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:01] – créée 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>
+===========================================================================
+===========================================================================
+            MAIN INPUT DATA FILE : 2D CHANNEL FLOW WITH A SQUARE BAR (OBSTACLE)
+                                   INCOMPRESSIBLE FLOW
+                                   HEATED SQUARE BAR  AT TEMPERATURE Th = 2.Tc
+                                   OTHER WALLS AT Tc
+                                   TEMPERATURE OF THE INFLOW : Tc
+                                   GRAVITY & BUOYANCY ARE NEGLECTED
+            DIMENSIONLESS LAYOUT :
+                 Length scale     : h (the channel height)
+                 Reynolds number    Re_h= rho_0.U_0.h/nu= 100
+                 dimensionless quantities :
+                   velocity U*      = U/U_0
+                   temperature T*= T/T0     with T0= Tc ---> Tc*= 1 and Th*= 2
+                   kinetic viscosity= 1/Re_h
+                 dimensionless domain : Lx/h= 10
+                 Initialisation = uniform velocity field
+                 inlet flowrate = uniform profil
+ Tc
+------------------------------------------
+                   |Th|
+                   |__|
+---> inflow                     outflow --->
+      (Tc)
+------------------------------------------
+ Tc
+J
+^
+|
+|
+---->I
+===========================================================================
+===========================================================================
+&Version File_Version="VERSION2.0"/
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+                     GENERAL LAYOUT
+                     (DIMENSIONLESS)
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+&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 /
+===========================================================================
+                    DOMAIN FEATURES
+                     (DIMENSIONLESS)
+===========================================================================
+&Domain_Features Start_Coordinate_I_Direction= 0.00 , End_Coordinate_I_Direction= 10.00,
+                 Start_Coordinate_J_Direction= 0.00 , End_Coordinate_J_Direction= 1.00,
+                 Start_Coordinate_K_Direction= 0.00 , End_Coordinate_K_Direction= 0.00,
+                 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
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+=============================================================================
+                 WALL BOUNDARY CONDITION SETUP
+                     (DIMENSIONLESS)
+=============================================================================
+---------------------------------------------------
+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
+                     (DIMENSIONLESS)
+=============================================================================
+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
+=============================================================================
+!--- 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
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+&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.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
+&Simulation_Management    Restart_Parameter= 0 ,
+                          Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-16,
+                          Temporal_Iterations_Number = 100000                  , Final_Time = 5.D+00  ,
+                          TimeStep_Type = 0 ,
+                          Timestep_Max = 1.D-03 ,
+                          Simulation_Backup_Rate                 = 1000   , Simulation_Checking_Rate = 101 /
+=============================================================================
+                   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= 3.0 , Xj= 0.5 , Xk= 0.0  /
+&Probe_Location  Xi= 6.0 , Xj= 0.5 , Xk= 0.0  /
+&Simulation_Management  Probe_TimeIterationRecordingRate= 10  ,
+                        Probe_StartTimeIterationRecording= 0  ,
+                        Probe_RecordingReset=.false. /
+=============================================================================
+           FIELDS RECORDING SETUP
+=============================================================================
+&Field_Recording_Setup    Precision_On_Instantaneous_Fields= 2 /
+!--- 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
 </file>
 [[sunfluidh:sunfluidh_examples_data_set|Click here to come back to the previous page]]