update input file for pipe flow

xcompact3d · Jan 30, 2024 · d11dd5c · d11dd5c
1 parent 0b5d84f
commit d11dd5c
Show file tree

Hide file tree

Showing 2 changed files with 282 additions and 0 deletions.
diff --git a/examples/Pipe-Flow/input_DNS_Re1000_LR.i3d b/examples/Pipe-Flow/input_DNS_Re1000_LR.i3d
@@ -0,0 +1,141 @@
+! -*- mode: f90 -*-
+
+!===================
+&BasicParam
+!===================
+
+! Flow type (1=Lock-exchange, 2=TGV, 3=Channel, 4=Periodic hill, 5=Cylinder, 6=dbg-schemes)
+itype = 14
+
+! Domain decomposition
+p_row=0               ! Row partition
+p_col=0               ! Column partition
+
+! Mesh
+nx=128                ! X-direction nodes
+ny=64                ! Y-direction nodes
+nz=64                ! Z-direction nodes
+istret = 0            ! y mesh refinement (0:no, 1:center, 2:both sides, 3:bottom)
+beta = 0.0            ! Refinement parameter (beta)
+
+! Domain
+xlx = 5.            ! Lx (Size of the box in x-direction)
+yly = 1.6             ! Ly (Size of the box in y-direction)
+zlz = 1.6             ! Lz (Size of the box in z-direction)
+
+! Flow parameters
+re  = 1000.           ! nu=1/re (Kinematic Viscosity)
+
+! Time stepping
+dt = 0.00025           ! Time step
+ifirst = 1            ! First iteration
+ilast = 5000        ! Last iteration
+
+! Enable modelling tools
+ilesmod=0             ! if 0 then DNS
+numscalar = 1         ! 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 = 3.0      ! Turbulence intensity (1=100%) !! Initial condition
+inflow_noise = 0.0    ! Turbulence intensity (1=100%) !! Inflow condition
+
+nclx1 = 0
+nclxn = 0
+ncly1 = 0
+nclyn = 0
+nclz1 = 0
+nclzn = 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 = 2             ! interpolation scheme (1: classic, 2: optimized, 3: optimized agressive)
+
+! Time scheme
+itimescheme = 3         ! Time integration scheme (1->Euler,2->AB2, 3->AB3, 4->AB4,5->RK3,6->RK4, 7-->CN2+AB3)
+
+! Dissipation control
+nu0nu = 3.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 = 10000   ! Frequency for writing backup file
+ioutput = 50         ! Frequency for visualization
+ilist = 25            ! Frequency for writing to screen        
+nvisu = 1             ! Size for visualisation collection
+
+/End
+
+!=================
+&Statistics
+!=================
+
+wrotation = 0.12      ! rotation speed to trigger turbulence
+spinup_time = 5000    ! number of time steps with a rotation to trigger turbulence
+nstat = 1             ! Size arrays for statistic collection
+initstat = 100001     ! Time steps after which statistics are collected
+
+/End
+
+!########################
+! OPTIONAL PARAMETERS
+!#######################
+
+!================
+&ibmstuff
+!================
+
+rai=0.50      ! when simulating a pipe, inner radius
+rao=0.70      ! when simulating a pipe, outer radius (defines wall thickness)
+nraf=1000     !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
+
+/End
+
+!================
+&ScalarParam
+!================
+
+Sc(1) = 0.71          ! Schmidt number
+
+nclxS1 = 0
+nclxSn = 0
+nclyS1 = 0
+nclySn = 0
+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
+
+&CASE
+/End
diff --git a/examples/Pipe-Flow/input_DNS_Re5300_HR.i3d b/examples/Pipe-Flow/input_DNS_Re5300_HR.i3d
@@ -0,0 +1,141 @@
+! -*- mode: f90 -*-
+
+!===================
+&BasicParam
+!===================
+
+! Flow type (1=Lock-exchange, 2=TGV, 3=Channel, 4=Periodic hill, 5=Cylinder, 6=dbg-schemes)
+itype = 14
+
+! Domain decomposition
+p_row=0               ! Row partition
+p_col=0               ! Column partition
+
+! Mesh
+nx=640                ! X-direction nodes
+ny=320                ! Y-direction nodes
+nz=320                ! Z-direction nodes
+istret = 0            ! y mesh refinement (0:no, 1:center, 2:both sides, 3:bottom)
+beta = 0.0            ! Refinement parameter (beta)
+
+! Domain
+xlx = 12.5            ! Lx (Size of the box in x-direction)
+yly = 1.6             ! Ly (Size of the box in y-direction)
+zlz = 1.6             ! Lz (Size of the box in z-direction)
+
+! Flow parameters
+re  = 5300.           ! nu=1/re (Kinematic Viscosity)
+
+! Time stepping
+dt = 0.0005           ! Time step
+ifirst = 1            ! First iteration
+ilast = 750000        ! Last iteration
+
+! Enable modelling tools
+ilesmod=0             ! if 0 then DNS
+numscalar = 1         ! 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 = 3.0      ! Turbulence intensity (1=100%) !! Initial condition
+inflow_noise = 0.0    ! Turbulence intensity (1=100%) !! Inflow condition
+
+nclx1 = 0
+nclxn = 0
+ncly1 = 0
+nclyn = 0
+nclz1 = 0
+nclzn = 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 = 2             ! interpolation scheme (1: classic, 2: optimized, 3: optimized agressive)
+
+! Time scheme
+itimescheme = 3         ! Time integration scheme (1->Euler,2->AB2, 3->AB3, 4->AB4,5->RK3,6->RK4, 7-->CN2+AB3)
+
+! Dissipation control
+nu0nu = 3.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 = 10000   ! Frequency for writing backup file
+ioutput = 500         ! Frequency for visualization
+ilist = 25            ! Frequency for writing to screen        
+nvisu = 1             ! Size for visualisation collection
+
+/End
+
+!=================
+&Statistics
+!=================
+
+wrotation = 0.12      ! rotation speed to trigger turbulence
+spinup_time = 5000    ! number of time steps with a rotation to trigger turbulence
+nstat = 1             ! Size arrays for statistic collection
+initstat = 100001     ! Time steps after which statistics are collected
+
+/End
+
+!########################
+! OPTIONAL PARAMETERS
+!#######################
+
+!================
+&ibmstuff
+!================
+
+rai=0.50      ! when simulating a pipe, inner radius
+rao=0.70      ! when simulating a pipe, outer radius (defines wall thickness)
+nraf=1000     !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
+
+/End
+
+!================
+&ScalarParam
+!================
+
+Sc(1) = 0.71          ! Schmidt number
+
+nclxS1 = 0
+nclxSn = 0
+nclyS1 = 0
+nclySn = 0
+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
+
+&CASE
+/End