doc and new input files

docs/pages/user_guide/list_files.rst

@@ -74,7 +74,32 @@ This directory contains a series of input files for various ready-to-run simulat
 
 ``Turbulent-Boundary-Layer``: Folder with input files for a spatially evolving turbulent boundary layer (zero pressure gradient).
 
+**input_LR_lowRe_implicit.i3d**: input fiel for a 3D DNS for low Reynolds numbers.
 
+``2D-hill``: Folder with input files for the turbulent flow around a 2D periodic hill.
+
+**input_LR_Re1000.i3d**: input file for a 3D DNS with a Reynolds number of 1000.
+
+
+``Cavity``: Folder with input files for a cavity flow.
+
+**input_2D_LR_Re14084.i3d**: input file for a 2D DNS with a Reynolds number of 14084.
+
+``Lock-exchange``: Folder with input files for gravity currents in the lock exchange set-up (finite release of heavy fluid in light ambient).
+
+**input_gravitycurrent_3DRe2236_LR.i3d**: input file for a 3D DNS with a Reynolds number of 2236.
+
+``Mixing-layer``: Folder with input files for a mixing-layer flow.
+
+**input_2D_LR_periodic.i3d**: input file for a 2D DNS with a Reynolds number of 400.
+
+``Wind-Turbine``: Folder with input files for wind turbines simulations.
+
+**two_turbines/input.i3d**: input files for a two-turbine set-up with aligned turbines.
+
+``ABL``: Folder with input files for atmospheric boundary layer simulations.
+
+**input_neutral.i3d**: input file for a 3D LES of a neutrally stable atmospheric boundary layer at 10 m/s.
 
 2DECOMP directory
 ------------------

examples/2D-hill/input_LR_Re1000.i3d

@@ -0,0 +1,158 @@
+! -*- mode: f90 -*-
+
+!===================
+&BasicParam
+!===================
+
+! Flow type (1=Lock-exchange, 2=TGV, 3=Channel, 4=Periodic hill, 5=Cylinder, 6=dbg-schemes)
+itype = 4
+
+! Domain decomposition
+p_row=0               ! Row partition
+p_col=0               ! Column partition
+
+! Mesh
+nx=128                ! X-direction nodes
+ny=129                 ! Y-direction nodes
+nz=64                 ! Z-direction nodes
+istret = 2            ! y mesh refinement (0:no, 1:center, 2:both sides, 3:bottom)
+beta = 2.    ! Refinement parameter (beta)
+
+! Domain
+xlx = 9.             ! Lx (Size of the box in x-direction)
+yly = 3.036              ! Ly (Size of the box in y-direction)
+zlz = 4.5              ! Lz (Size of the box in z-direction)
+
+! Flow parameters
+re  = 1000.           ! nu=1/re (Kinematic Viscosity)
+cpg = F               ! if cpg=T, then re is friction Reynolds number
+
+! Time stepping
+dt = 0.005            ! Time step
+ifirst = 1            ! First iteration
+ilast = 5000        ! Last iteration
+
+! Enable modelling tools
+ilesmod=0             ! if 0 then DNS
+numscalar = 0         ! How many scalars? (Set to zero to disable scalars)
+iibm=2                ! Flag for immersed boundary method
+
+! Boundary and initial conditions
+iin = 1               ! Inflow conditions (1: classic, 2: turbinit, 4: SEM)
+u1  = 1.              ! u1 (max velocity) (for inflow condition)
+u2  = 1.              ! u2 (min velocity) (for inflow condition)
+init_noise = 0.125    ! Turbulence intensity (1=100%) !! Initial condition
+inflow_noise = 0.0    ! Turbulence intensity (1=100%) !! Inflow condition
+idir_stream  = 1      ! Index of the streamwise direction (1=X, 3=Z, 2 not available)
+
+nclx1 = 0
+nclxn = 0
+ncly1 = 2
+nclyn = 2
+nclz1 = 0
+nclzn = 0
+
+/End
+
+!====================
+&NumOptions
+!====================
+
+! Spatial derivatives
+ifirstder = 4         ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact)
+isecondder = 5        ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact, 5->hyperviscous 6th)
+ipinter = 3           ! interpolation scheme (1: classic, 2: optimized, 3: optimized agressive)
+
+! Time scheme
+iimplicit = 0
+itimescheme = 5       ! Time integration scheme (1->Euler,2->AB2, 3->AB3, 4->AB4,5->RK3,6->RK4, 7-->CN2+AB3)
+
+! Dissipation control
+nu0nu = 4.0             ! Ratio between hyperviscosity/viscosity at nu
+cnu = 0.44               ! Ratio between hypervisvosity at k_m=2/3pi and k_c= pi
+
+/End
+
+!=================
+&InOutParam
+!=================
+
+! Basic I/O
+irestart = 0          ! Read initial flow field ?
+icheckpoint = 1000    ! Frequency for writing backup file
+ioutput = 250       ! Frequency for visualization
+ilist = 25            ! Frequency for writing to screen        
+nvisu = 1             ! Size for visualisation collection
+nprobes = 0           ! Number of probes
+
+/End
+
+!=================
+&Statistics
+!=================
+
+wrotation = 0.12      ! rotation speed to trigger turbulence
+spinup_time = 0    ! number of time steps with a rotation to trigger turbulence
+nstat = 1             ! Size arrays for statistic collection
+initstat = 40000      ! Time steps after which statistics are collected
+
+/End
+
+
+!########################
+! OPTIONAL PARAMETERS
+!#######################
+
+!================
+&ScalarParam
+!================
+
+Sc(1) = 1.0           ! Schmidt number
+
+nclxS1 = 0
+nclxSn = 0
+nclyS1 = 2
+nclySn = 2
+nclzS1 = 0
+nclzSn = 0
+
+/End
+
+!================
+&LESModel
+!================
+
+iles = 0             ! LES Model (1: Phys Smag, 2: Phys WALE, 3: Phys dyn. Smag, 4: iSVV, 5: dyn SEV)
+smagcst = 0.14        ! Smagorinsky constant
+SmagWallDamp = 0      ! 1: Mason and Thomson Damping function, otherwise OFF
+walecst = 0.5         ! WALES Model Coefficient
+iconserv = 0          ! Formulation SGS divergence (0: non conservative, 1: conservative)
+
+/End
+
+!================
+&ibmstuff
+!================
+
+cex=5.        !when simulating a cylinder x coordinate of the center
+cey=6.        !when simulating a cylinder y coordinate of the center
+ra=0.5        !when simulating a cylinder, radius
+nraf=10       !level of refinement for iibm==2 to find the surface of the immersed object
+nobjmax=2     !number of immersed objects (DO NOT USE ZERO OBJECTS)
+iforces=0
+nvol=1        !Number of volumes for computing force balance
+
+/End
+
+&ForceCVs
+
+xld(1) = 4.0 !X left for volume control
+xrd(1) = 6.0 !X right for volume control
+yld(1) = 5.0 !Y bottom for volume control
+yud(1) = 7.0 !Y top for volume control
+
+/End
+
+
+&CASE
+/End

examples/Cavity/input_2D_LR_Re14084.i3d

@@ -0,0 +1,126 @@
+! -*- mode: f90 -*-
+
+!===================
+&BasicParam
+!===================
+
+! Flow type (1=Lock-exchange, 2=TGV, 3=Channel, 4=Periodic hill, 5=Cylinder, 6=dbg-schemes)
+itype = 13
+
+! Domain decomposition
+p_row=0               ! Row partition
+p_col=0               ! Column partition
+
+! Mesh
+nx=257                ! X-direction nodes
+ny=129                ! Y-direction nodes
+nz=1                  ! Z-direction nodes
+istret = 0            ! y mesh refinement (0:no, 1:center, 2:both sides, 3:bottom)
+
+! Domain
+xlx = 1.              ! Lx (Size of the box in x-direction)
+yly = 1.              ! Ly (Size of the box in y-direction)
+zlz = 1.              ! Lz (Size of the box in z-direction)
+
+! Flow parameters
+re  = 14084.5         ! nu=1/re (Kinematic Viscosity)
+
+! Time stepping
+dt = 0.002            ! Time step
+ifirst = 1            ! First iteration
+ilast = 5000          ! Last iteration
+
+! Enable modelling tools
+numscalar = 1         ! How many scalars? (Set to zero to disable scalars)
+
+! Velocity boundary conditions
+nclx1 = 2
+nclxn = 2
+ncly1 = 2
+nclyn = 2
+nclz1 = 0
+nclzn = 0
+
+! Gravity vector
+gravx = 0.
+gravy = 1.
+gravz = 0.
+
+/End
+
+!====================
+&NumOptions
+!====================
+
+! Spatial derivatives
+ifirstder = 4         ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact)
+isecondder = 4        ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact, 5->hyperviscous 6th)
+ipinter = 1           ! interpolation scheme (1: classic, 2: optimized, 3: optimized agressive)
+
+! Time scheme
+itimescheme = 2       ! Time integration scheme (1->Euler,2->AB2, 3->AB3, 4->AB4,5->RK3,6->RK4, 7-->CN2+AB3)
+
+! Dissipation control
+nu0nu = 4.0             ! Ratio between hyperviscosity/viscosity at nu
+cnu = 0.44              ! Ratio between hypervisvosity at k_m=2/3pi and k_c= pi
+
+/End
+
+!=================
+&InOutParam
+!=================
+
+! Basic I/O
+irestart = 0          ! Read initial flow field ?
+icheckpoint = 400000  ! Frequency for writing backup file
+ioutput = 500         ! Frequency for visualization
+ilist = 25            ! Frequency for writing to screen        
+nvisu = 1             ! Size for visualisation collection
+nprobes = 1           ! Number of probes
+
+/End
+
+!=================
+&Statistics
+!=================
+
+nstat = 1             ! Size arrays for statistic collection
+initstat = 300000     ! Time steps after which statistics are collected
+
+/End
+
+&CASE
+/End
+
+!=================
+&ProbesParam
+!=================
+
+flag_all_digits = F     ! By default, only 6 digits, 16 when True
+flag_extra_probes = F   ! By default, gradients are not monitored
+xyzprobes(1,1) = 0.1    ! X position of probe 1
+xyzprobes(2,1) = 0.1    ! Y position of probe 1
+xyzprobes(3,1) = 0.1    ! Z position of probe 1
+
+/End
+
+!########################
+! OPTIONAL PARAMETERS
+!#######################
+
+!================
+&ScalarParam
+!================
+
+Sc(1) = 0.71          ! Schmidt number
+sc_skew(1) = F        ! Skew-symetric convection
+ri(1) = 0.71          ! Richardson number
+
+nclxS1 = 2
+nclxSn = 2
+nclyS1 = 2
+nclySn = 2
+nclzS1 = 0
+nclzSn = 0
+
+/End