Ceci est une ancienne révision du document !
Click here to come back to the previous page <file txt input3d.dat>
MAIN INPUT DATA FILE : 2D EDDY ON A PERIODIC DOMAIN (DIMENSIONAL RELEASE)
INCOMPRESSIBLE FLOW
ISOTHERM
DIMENSIONAL FORMULATION
Initialization of the velocity field is carried out from the user's routine Eddy_Velocity_Field in the module named module_user_define_init_fields.f90
access by the data file ---> Initial_Field_Option_For_Velocity_I= 4
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
GENERAL LAYOUT
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ &Version File_Version=“VERSION2.0”/
FLUID PROPERTIES
&Fluid_Properties Reference_Density= 860.0 , Reference_Dynamic_Viscosity = 43.D-03 /
INITIALIZATION OF THE VELOCITY COMPONENTS, THE TEMPERATURE AND SPECIES
Initial_Field_Option_For_Velocity_I : Special initialisation of a 2D eddy defined by the user's routine Eddy_Velocity_Field
in the module named module_user_define_init_fields.f90
&Velocity_Initialization I_Velocity_Reference_Value= 0.00 , J_Velocity_Reference_Value= 0.00, K_Velocity_Reference_Value= 0.00 , Initial_Field_Option_For_Velocity_I= 4 /
GRAVITY
No gravity
DOMAIN FEATURES
&Domain_Features Geometric_Layout= 0,
Start_Coordinate_I_Direction=-0.50 , End_Coordinate_I_Direction= 0.50,
Start_Coordinate_J_Direction=-0.50 , End_Coordinate_J_Direction= 0.50,
Start_Coordinate_K_Direction= 0.00 , End_Coordinate_K_Direction= 0.00,
Cells_Number_I_Direction= 128 ,Cells_Number_J_Direction= 128 ,Cells_Number_K_Direction= 1,
Regular_Mesh= .true. /
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
DEFINITION OF BOUNDARY CONDITIONS
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
WALL BOUNDARY CONDITION SETUP
No wall
INLET AND OUTLET BOUNDARY CONDITIONS
No inlet, No outlet
BORDER BOUNDARY CONDITIONS
Keep in mind that the domain is enclosed by default. The walls at the ends of the domain are here replaced by periodical conditions
&Border_Domain_Boundary_Conditions West_BC_Name= “Periodic” , East_BC_Name= “Periodic” , Back_BC_Name= “Periodic” , Front_BC_Name= “Periodic” / ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
NUMERICAL METHODS
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ &Numerical_Methods NS_NumericalMethod= “BDF2-SchemeO2” , !— BDF2 + 2nd order centered scheme
MomentumConvection_Scheme="Centered-O2-Conservative" , !--- conservative form for solving the velocity (momentum) equation
Poisson_NumericalMethod="Home-SORMultigrid-ConstantMatrixCoef" / !--- SOR + multigrid method (homemade release) for solving the Poisson's equation with constant coefficient matrix
&HomeData_PoissonSolver SolverName=“SOR” , !— Successive Over-Relaxation (SOR) method based on the red-black algorithm
Relaxation_Coefficient= 1.8 , !--- Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2)
Number_max_Grid= 6, !--- Number of grid levels
Number_max_Cycle= 5, !--- Number of multigrid cycles
Number_Iteration= 0, !--- Maximum number of SOR iterations method applied for any grid level, if 0 (or removed) the 3 next data are considered
Number_Iteration_FineToCoarseGrid= 15, !--- number of SOR iterations applied on any grid level during the restriction step (before the coarsest grid computation)
Number_Iteration_CoarseToFineGrid= 15, !--- number of SOR iterations applied on any grid level during the prolongation step (after the Coarsest grid computation)
Number_Iteration_CoarsestGrid= 15 , !--- number of SOR iterations applied on the coarsest grid
Convergence_Criterion= 1.D-10 / !--- convergence tolerance on the residu of the Poisson's equation
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
SIMULATION MANAGEMENT
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Time step calculated from a CFL coefficient
&Simulation_Management Restart_Parameter= 0 ,
Steady_Flow_Stopping_Criterion_Enabled = .false. , Steady_Flow_Stopping_Criterion = 1.D-22,
Temporal_Iterations_Number = 50000 , Final_Time = 2.0D+01 ,
TimeStep_Type = 1 ,
CFL_min= 0.25 , CFL_max= 0.25, Iterations_For_Timestep_Linear_Progress= 1 ,
Timestep_Min = 2.D-01 , Timestep_Max = 2.D-01 ,
Simulation_Backup_Rate = 200 , Simulation_Checking_Rate = 51 /
PROBES MANAGEMENT
Probes order U , V , W , T , P , RHO
&Probe_Quantities_Enabled Temporal_Series_For_Quantity_Enabled(:) = .true., .true., .false., .false., .true. , .false. /
&Probe_Location Xi= 0.0 , Xj= 0.0 , Xk= 0.0 / &Probe_Location Xi= 0.1 , Xj= 0.1 , Xk= 0.0 / &Simulation_Management Probe_Recording_Rate = 1 /
FIELDS RECORDING DECLARATION
&Field_Recording_Setup Precision_On_Instantaneous_Fields= 2 / &Simulation_Management Fields_Recording_Rate = 1.0D-00 /
&Instantaneous_Fields_Listing Name_of_Field = “U” , Recording_Enabled = .true. / First velocity component &Instantaneous_Fields_Listing Name_of_Field = “V” , Recording_Enabled = .true. / Second velocity component &Instantaneous_Fields_Listing Name_of_Field = “W” , Recording_Enabled = .false. / Third velocity component &Instantaneous_Fields_Listing Name_of_Field = “P” , Recording_Enabled = .true. , End_of_Data_Block= .true./ Pressure related to the velocity field
</file txt> Click here to come back to the previous page
- fr