Rusanov Solver

examples/advection_rusanov/GNUmakefile

@@ -0,0 +1,50 @@
+# NGA location if not yet defined
+NGA_HOME ?= ../..
+
+# Compilation parameters
+PRECISION  = DOUBLE
+USE_MPI    = TRUE
+USE_FFTW   = FALSE
+USE_HYPRE  = FALSE
+USE_LAPACK = TRUE
+PROFILE    = FALSE
+DEBUG      = TRUE
+COMP       = gnu
+EXEBASE    = nga
+
+# Directories that contain user-defined code
+Udirs   := src
+
+# Include user-defined sources
+Upack   += $(foreach dir, $(Udirs), $(wildcard $(dir)/Make.package))
+Ulocs   += $(foreach dir, $(Udirs), $(wildcard $(dir)))
+include $(Upack)
+INCLUDE_LOCATIONS += $(Ulocs)
+VPATH_LOCATIONS   += $(Ulocs)
+
+# Define external libraries - this can also be in .profile
+HDF5_DIR =  /opt/hdf5
+HYPRE_DIR = /opt/hypre
+LAPACK_DIR = /opt/lapack-3.10.1
+FFTW_DIR = /opt/fftw-3.3.10
+
+# NGA compilation definitions
+include $(NGA_HOME)/tools/GNUMake/Make.defs
+
+# Include NGA base code
+Bdirs   := particles core hyperbolic data solver config grid libraries
+Bpack   += $(foreach dir, $(Bdirs), $(NGA_HOME)/src/$(dir)/Make.package)
+include $(Bpack)
+
+# Inform user of Make.packages used
+ifdef Ulocs
+   $(info Taking user code from: $(Ulocs))
+endif
+$(info Taking base code from: $(Bdirs))
+
+# Target definition
+all: $(executable)
+	@echo COMPILATION SUCCESSFUL
+
+# NGA compilation rules
+include $(NGA_HOME)/tools/GNUMake/Make.rules

examples/advection_rusanov/README

@@ -0,0 +1,2 @@
+standard advection test case advecting different shapes over a periodic domain
+to evaluate accuracy and diffusion

examples/advection_rusanov/input

@@ -0,0 +1,24 @@
+# Parallelization
+Partition :	2 1 1
+
+# Mesh definition
+Lx :    1.0
+nx :    128
+Ly :    1.0
+ny :    128
+Lz :    1.0
+nz :    1
+Periodic : true
+
+# Initialization
+Shapes : c
+
+# Direction
+Advection direction : 1 1 0
+
+# Time integration
+Max cfl number :    0.5
+
+# Ensight output
+Ensight output period : 0.01
+

examples/advection_rusanov/src/Make.package

@@ -0,0 +1,2 @@
+# List here the extra files here
+f90EXE_sources += geometry.f90 simulation.f90

examples/advection_rusanov/src/geometry.f90

@@ -0,0 +1,80 @@
+!> Various definitions and tools for initializing NGA2 config
+module geometry
+   use config_class, only: config
+   use precision,    only: WP
+   implicit none
+   private
+
+   !> Single config
+   type(config), public :: cfg
+
+   public :: geometry_init
+
+contains
+
+   !> Initialization of problem geometry
+   subroutine geometry_init
+      use sgrid_class, only: sgrid
+      use param,       only: param_read
+      implicit none
+      type(sgrid) :: grid
+
+      ! Create a grid from input params
+      create_grid: block
+         use sgrid_class, only: cartesian
+         logical :: periodic
+         integer :: i, j, k, nx, ny, nz
+         real(WP) :: Lx, Ly, Lz
+         real(WP), dimension(:), allocatable :: x, y, z
+
+         ! read in grid definition
+         call param_read('Lx', Lx)
+         call param_read('nx', nx)
+         call param_read('Ly', Ly)
+         call param_read('ny', ny)
+         call param_read('Lz', Lz)
+         call param_read('nz', nz)
+         call param_read('Periodic', periodic)
+
+         ! allocate
+         allocate(x(nx+1), y(ny+1), z(nz+1))
+
+         ! create simple rectilinear grid
+         do i = 1, nx+1
+            x(i) = real(i-1,WP) / real(nx, WP) * Lx
+         end do
+         do j = 1, ny+1
+            y(j) = real(j-1,WP) / real(ny, WP) * Ly
+         end do
+         do k = 1, nz+1
+            z(k) = real(k-1,WP) / real(nz, WP) * Lz
+         end do
+
+         ! generate serial grid object
+         grid=sgrid(coord=cartesian, no=1, x=x, y=x, z=x, xper=periodic,      &
+           & yper=periodic, zper=periodic)
+
+      end block create_grid
+
+      ! create a config from that grid on our entire group
+      create_cfg: block
+         use parallel, only: group
+         integer, dimension(3) :: partition
+
+         ! read in partition
+         call param_read('Partition', partition, short='p')
+
+         ! create partitioned grid
+         cfg = config(grp=group, decomp=partition, grid=grid)
+
+      end block create_cfg
+
+      ! Create masks for this config
+      create_walls: block
+         cfg%VF=1.0_WP
+      end block create_walls
+
+   end subroutine geometry_init
+
+end module geometry
+