Différences

Ci-dessous, les différences entre deux révisions de la page.

--- sunfluidh:2d_channel_flow_with_heated_bar_incomp_flow [2016/11/30 11:31] – yann
+++ sunfluidh:2d_channel_flow_with_heated_bar_incomp_flow [2019/11/13 16:32] (Version actuelle) – yann
@@ Ligne 1: / Ligne 1: @@
 [[sunfluidh:sunfluidh_examples_data_set|Click here to come back to the previous page]]
-<file>
+<file txt input3d.dat>
 ===========================================================================
@@ Ligne 25: / Ligne 25: @@
                  Initialisation = uniform velocity field
                  inlet flowrate = uniform profil
-INCOMPRESSIBLE DOWNSTEP FLOW
  Tc
 ------------------------------------------
@@ Ligne 50: / Ligne 48: @@
                      (DIMENSIONLESS)
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-&Fluid_Properties    Reference_Dynamic_Viscosity = 1.00D-02, Reference_Density= 1.0  , Reference_Temperature= 1.0 /
+&Fluid_Properties    Reference_Dynamic_Viscosity = 1.00D-02, Reference_Density= 1.0  ,
+                     Heat_Transfer_Flow = .true., Reference_Temperature= 1.0, Prandtl = 0.71 /
 &Velocity_Initialization  I_Velocity_Reference_Value  = 1.0 , J_Velocity_Reference_Value = 0.0 , K_Velocity_Reference_Value = 0.0 /
@@ Ligne 62: / Ligne 61: @@
                  Cells_Number_I_Direction= 256 ,Cells_Number_J_Direction= 64 ,Cells_Number_K_Direction= 1,
                  Regular_Mesh= .true. /
+-------------------------------------------------------
+       GEOMETRY OF THE IMMERSED BODIES
+-------------------------------------------------------
+First (and unic) immersed body (temperature Th= 2.Tc)
+&Polyhedral_Immersed_Bodies  Xi_1= 4.5  , Xj_1= 0.5 ,Xk_1= 0.0  , Xi_2= 5.5  , Xj_2= 0.5 ,Xk_2= 0.0 ,
+                             Xi_3= 5.5  , Xj_3= 1.0  ,Xk_3= 0.0 , Xi_4= 4.5  , Xj_4= 1.0  ,Xk_4= 0.0 ,
+                             Wall_BC_DataSetName ="Set2" /
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
               DEFINITION OF BOUNDARY CONDITIONS
@@ Ligne 68: / Ligne 77: @@
 =============================================================================
                  WALL BOUNDARY CONDITION SETUP
+                     (DIMENSIONLESS)
 =============================================================================
- No wall
+---------------------------------------------------
+DATA SET FOR THE WALL BOUNDARY CONDITIONS
+---------------------------------------------------
+First set of wall boundary conditions (applied to the walls of the channel) :
+                                        imposed temperature : Tc
+--------
+&Heat_Wall_Boundary_Condition_Setup
+     Wall_BC_DataSetName ="Set1",
+      West_Heat_BC_Option = 0   ,  East_Heat_BC_Option = 0  ,  Back_Heat_BC_Option = 0  ,  Front_Heat_BC_Option = 0  ,
+      West_Wall_BC_Value= 1.0   ,  East_Wall_BC_Value= 1.0  ,  Back_Wall_BC_Value= 1.0  ,  Front_Wall_BC_Value= 1.0  /
+The usual wall boundary conditions for the velocity are used (no-slip and impermeability conditions).
+As they are the conditions by default, they are not explicitly written
+--------
+Second set of wall boundary conditions (applied to the walls of the square bar) :
+                                        imposed temperature : Th
+--------
+&Heat_Wall_Boundary_Condition_Setup
+      Wall_BC_DataSetName ="Set2",
+      West_Heat_BC_Option = 0   ,  East_Heat_BC_Option = 0  ,  Back_Heat_BC_Option = 0  ,  Front_Heat_BC_Option = 0  ,
+      West_Wall_BC_Value= 2.0   ,  East_Wall_BC_Value= 2.0  ,  Back_Wall_BC_Value= 2.0  ,  Front_Wall_BC_Value= 2.0  /
+The usual wall boundary conditions for the velocity are used (no-slip and impermeability conditions).
+As they are the conditions by default, they are not explicitly written
 =============================================================================
-               INLET AND OUTLET BOUNDARY CONDITIONS
+               INLET AND OUTLET  BOUNDARY CONDITIONS
+                     (DIMENSIONLESS)
 =============================================================================
-No inlet, No outlet
+Keep in mind that the domain is enclosed by default.
+Here the inlet and outlet conditions are located at the ends of the domain. They replace the walls by default over the insterested areas.
+Inlet : Uniform flowrate profil
+&Inlet_Boundary_Conditions  Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 ,
+                            Plan_Location_Coordinate=  0.0    ,
+                            Start_Coordinate_of_First_Span =  0.00  ,  End_Coordinate_of_First_Span = 1.00 ,
+                            Start_Coordinate_of_Second_Span= 0.0    ,  End_Coordinate_of_Second_Span= 0.0  ,
+                            Normal_Velocity_Reference_Value= 1.0  ,
+                            Temperature_Reference_Value= 1.0  /
+&Outlet_Boundary_Conditions  Type_of_BC= "OUTLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 ,
+                             Plan_Location_Coordinate=  10.0   ,
+                             Start_Coordinate_of_First_Span =  0.00  ,  End_Coordinate_of_First_Span = 1.00 ,
+                             Start_Coordinate_of_Second_Span= 0.0    ,  End_Coordinate_of_Second_Span= 0.0  /
 =============================================================================
                BORDER BOUNDARY CONDITIONS
 =============================================================================
-&Border_Domain_Boundary_Conditions West_Border= 4 , East_Border= 4 , Back_Border= 4 , Front_Border= 4 , North_Border = 0 , South_Border = 0 /
+!--- No new boundary conditions are defined at the ends of the domain : walls by default are preserved, the inlet and outlet previously are defined above)
+!--- As "None" is the default setting for this namelist, it can be removed
+&Border_Domain_Boundary_Conditions West_BC_Name= "None" , East_BC_Name= "None" , Back_BC_Name= "None" , Front_BC_Name= "None" , North_BC_Name= "None" , South_BC_Name= "None" /
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
                    NUMERICAL METHODS
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
- Default initialization:
+&Numerical_Methods  NS_NumericalMethod= "BDF2-SchemeO2"                    ,       !--- BDF2 + 2nd order centered scheme
- Numerical scheme for Navier-Stokes equations : 2nd centrered order in space, 2nd order in time (BDF2 : 2nd oder Backward Differentiation Formulation)
+                    MomentumConvection_Scheme="Centered-O2-Conservative"   ,       !--- conservative form for solving the velocity (momentum) equation
-                  for Poisson equation        : SOR iterative method coupled with a multigrid approach
+                    Poisson_NumericalMethod="Home-SORMultigrid-ConstantMatrixCoef"  / !--- SOR  + multigrid method (homemade release) for solving the Poisson's equation with constant coefficient matrix
-Parameters for solving Poisson equation (multigrid parameters) :
-&Numerical_Methods  Numerical_Scheme= 1 ,
-                    Convective_Flux_Discretization_Type = 1
-                    Numerical_Method_Poisson_Equation= 1,
-                    Number_max_Grid= 6 , Number_max_Cycle= 5 , Number_Iteration= 15 ,
-                    Relaxation_Coefficient = 1.80 , Convergence_Criterion = 1.D-10 /
+&HomeData_PoissonSolver  SolverName="SOR"               ,          !--- Successive Over-Relaxation (SOR) method based on the red-black algorithm
+                            Relaxation_Coefficient= 1.7 ,          !--- Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2)
+                            Number_max_Grid= 4,                    !--- Number of grid levels
+                            Number_max_Cycle= 10,                  !--- 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-08 /        !--- convergence tolerance on the residu of the Poisson's equation
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
                    SIMULATION MANAGEMENT
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+The numerical time step is imposed
-Time step calculated from a CFL coefficient
 &Simulation_Management    Restart_Parameter= 0 ,
-                          Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-22,
+                          Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-16,
-                          Temporal_Iterations_Number = 50000                  , Final_Time = 2.0D+01  ,
+                          Temporal_Iterations_Number = 100000                  , Final_Time = 5.D+00  ,
-                          TimeStep_Type = 1 ,
+                          TimeStep_Type = 0 ,
-                          CFL_min= 0.25   , CFL_max= 0.25, Iterations_For_Timestep_Linear_Progress= 1 ,
+                          Timestep_Max = 1.D-03 ,
-                          Timestep_Min = 2.D-01                           , Timestep_Max = 2.D-01 ,
+                          Simulation_Backup_Rate                 = 1000   , Simulation_Checking_Rate = 101 /
-                          Simulation_Backup_Rate                 =  200   , Simulation_Checking_Rate = 51  /
 =============================================================================
                    PROBES MANAGEMENT
 =============================================================================
-                                                      Probes order    U    , V     , W      , T      , P      , RHO
+                                                     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= 3.0 , Xj= 0.5 , Xk= 0.0  /
-&Probe_Location  Xi= 0.1 , Xj= 0.1 , Xk= 0.0 , End_of_Data_Block= .true. /
+&Probe_Location  Xi= 6.0 , Xj= 0.5 , Xk= 0.0  /
-&Simulation_Management   Probe_Recording_Rate = 1 /
+&Simulation_Management  Probe_TimeIterationRecordingRate= 10  ,
+                        Probe_StartTimeIterationRecording= 0  ,
+                        Probe_RecordingReset=.false. /
 =============================================================================
-           FIELDS RECORDING DECLARATION
+           FIELDS RECORDING SETUP
 =============================================================================
-&Field_Recording_Setup  Precision_On_Instantaneous_Fields= 2 /
+&Field_Recording_Setup    Precision_On_Instantaneous_Fields= 2 /
-&Simulation_Management  Fields_Recording_Rate = 1.0D-00 /
+!--- Snapshots
+&Simulation_Management
+    InstantaneousFields_RecordingReset=.false.     ,
+    InstantaneousFields_TimeRecordingRate= 5.0E-01 ,
+    InstantaneousFields_RecordingStartTime= 0.D-00  /
+&Instantaneous_Fields_Listing  Name_of_Field = "U     " /      First velocity component
+&Instantaneous_Fields_Listing  Name_of_Field = "V     " /      Second velocity component
+&Instantaneous_Fields_Listing  Name_of_Field = "T     " /      Temperature
+!--- Statistics
+&Simulation_Management    Start_Time_For_Statistics= 1.D+03               , Time_Range_Statistic_Calculation = 5.D+00  /
+&Statistical_Fields_Listing  Name_of_Field = "<U>   "/
+&Statistical_Fields_Listing  Name_of_Field = "<V>   "/
+END OF FILE
-&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.  /      Pressure related to the velocity field
 </file>
 [[sunfluidh:sunfluidh_examples_data_set|Click here to come back to the previous page]]