Documentation du code de simulation numérique SUNFLUIDH

Outils pour utilisateurs

Outils du site

Piste : • 2d_channel_flow_with_heated_bar_incomp_flow
sunfluidh:2d_channel_flow_with_heated_bar_incomp_flow

Différences

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

Lien vers cette vue comparative

Les deux révisions précédentesRévision précédente
Prochaine révision		Révision précédente
sunfluidh:2d_channel_flow_with_heated_bar_incomp_flow [2016/11/30 12:25] – yannsunfluidh: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]] [[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                  Initialisation = uniform velocity field
                  inlet flowrate = uniform profil                  inlet flowrate = uniform profil
- +                 
-INCOMPRESSIBLE DOWNSTEP FLOW +
  Tc  Tc
 ------------------------------------------ ------------------------------------------
Ligne 50: Ligne 48:
                      (DIMENSIONLESS)                      (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 / &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,                  Cells_Number_I_Direction= 256 ,Cells_Number_J_Direction= 64 ,Cells_Number_K_Direction= 1,
                  Regular_Mesh= .true. /                  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               DEFINITION OF BOUNDARY CONDITIONS
Ligne 77: Ligne 87:
  
 &Heat_Wall_Boundary_Condition_Setup &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_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  /       West_Wall_BC_Value= 1.0   ,  East_Wall_BC_Value= 1.0  ,  Back_Wall_BC_Value= 1.0  ,  Front_Wall_BC_Value= 1.0  /
Ligne 89: Ligne 100:
  
 &Heat_Wall_Boundary_Condition_Setup &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_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  , +      West_Wall_BC_Value= 2.0   ,  East_Wall_BC_Value= 2.0  ,  Back_Wall_BC_Value= 2.0  ,  Front_Wall_BC_Value= 2.0  /
-      End_of_Data_Block= .true.   /+
  
 The usual wall boundary conditions for the velocity are used (no-slip and impermeability conditions). 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 As they are the conditions by default, they are not explicitly written
- 
-------------------------------------------------------- 
-       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 , 
-                             Identifier_Value= 2, 
-                             End_of_Data_Block= .true. / 
  
 ============================================================================= =============================================================================
Ligne 121: Ligne 121:
                             Start_Coordinate_of_Second_Span= 0.0    ,  End_Coordinate_of_Second_Span= 0.0  ,                             Start_Coordinate_of_Second_Span= 0.0    ,  End_Coordinate_of_Second_Span= 0.0  ,
                             Normal_Velocity_Reference_Value= 1.0  ,                             Normal_Velocity_Reference_Value= 1.0  ,
-                            Temperature_Reference_Value= 1.0  , +                            Temperature_Reference_Value= 1.0  /
-                            End_of_Data_Block= .true. /+
  
 &Outlet_Boundary_Conditions  Type_of_BC= "OUTLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 , &Outlet_Boundary_Conditions  Type_of_BC= "OUTLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 ,
                              Plan_Location_Coordinate=  10.0   ,                              Plan_Location_Coordinate=  10.0   ,
                              Start_Coordinate_of_First_Span =  0.00  ,  End_Coordinate_of_First_Span = 1.00 ,                              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  , +                             Start_Coordinate_of_Second_Span= 0.0    ,  End_Coordinate_of_Second_Span= 0.0  /
-                             End_of_Data_Block= .true. /+
 ============================================================================= =============================================================================
                BORDER BOUNDARY CONDITIONS                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) !--- 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_Border= 0 , East_Border= 0 , Back_Border= 0 , Front_Border= 0 , North_Border = 0 , South_Border = 0 /+&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
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-          For Navier-Stokes Equations  : Numerical scheme --> BDF2 - Spatial discretization --> 2nd order centered scheme +                     
-                                         convection flux ---> conservative form +&Numerical_Methods  NS_NumericalMethod= "BDF2-SchemeO2"                    ,       !--- BDF2 + 2nd order centered scheme 
-          For Poisson's equation       : Iterative method --> SOR coupled with a multigrid method.+                    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                    
  
-Parameters of the multigrid method : +&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) 
-&Numerical_Methods  Numerical_Scheme= 1 , +                            Number_max_Grid= 4,                    !--- Number of grid levels 
-                    Convective_Flux_Discretization_Type = 1 +                            Number_max_Cycle= 10,                  !--- Number of multigrid cycles 
-                    Numerical_Method_Poisson_Equation= 1, +                            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_max_Grid= 4 , Number_max_Cycle= 10 , Number_Iteration= 15 , +                            Number_Iteration_FineToCoarseGrid= 15, !--- number of SOR iterations applied on any grid level during the restriction step (before the coarsest grid computation) 
-                    Relaxation_Coefficient = 1.80 , Convergence_Criterion = 1.D-08 /+                            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                    SIMULATION MANAGEMENT
Ligne 156: Ligne 158:
 &Simulation_Management    Restart_Parameter= 0 , &Simulation_Management    Restart_Parameter= 0 ,
                           Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-16,                           Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-16,
-                          Temporal_Iterations_Number = 100000                  , Final_Time = 5.D+02  ,+                          Temporal_Iterations_Number = 100000                  , Final_Time = 5.D+00  ,
                           TimeStep_Type = 0 ,                           TimeStep_Type = 0 ,
                           Timestep_Max = 1.D-03 ,                           Timestep_Max = 1.D-03 ,
Ligne 168: Ligne 170:
  
 &Probe_Location  Xi= 3.0 , Xj= 0.5 , Xk= 0.0  / &Probe_Location  Xi= 3.0 , Xj= 0.5 , Xk= 0.0  /
-&Probe_Location  Xi= 6.0 , Xj= 0.5 , Xk= 0.0 , End_of_Data_Block= .true. / +&Probe_Location  Xi= 6.0 , Xj= 0.5 , Xk= 0.0  / 
-&Simulation_Management   Probe_Recording_Rate = 10 /+&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 /
Ligne 178: Ligne 182:
 !--- Snapshots !--- Snapshots
  
-&Simulation_Management    Fields_Recording_Rate = 5.D+01  / +&Simulation_Management 
-&Instantaneous_Fields_Listing  Name_of_Field = "U     " , Recording_Enabled = .true.  /      First velocity component +    InstantaneousFields_RecordingReset=.false.     ,        
-&Instantaneous_Fields_Listing  Name_of_Field = "V     " , Recording_Enabled = .true. , End_of_Data_Block= .true  /      Second velocity component+    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 !--- Statistics
Ligne 186: Ligne 195:
 &Simulation_Management    Start_Time_For_Statistics= 1.D+03               , Time_Range_Statistic_Calculation = 5.D+00  / &Simulation_Management    Start_Time_For_Statistics= 1.D+03               , Time_Range_Statistic_Calculation = 5.D+00  /
  
-&Statistical_Fields_Listing  Name_of_Field = "<U>   " , Recording_Enabled = .true.  / +&Statistical_Fields_Listing  Name_of_Field = "<U>   "/ 
-&Statistical_Fields_Listing  Name_of_Field = "<V>   " , Recording_Enabled = .true. , End_of_Data_Block= .true /+&Statistical_Fields_Listing  Name_of_Field = "<V>   "/ 
 + 
 + 
 +END OF FILE
  
 </file> </file>
 [[sunfluidh:sunfluidh_examples_data_set|Click here to come back to the previous page]] [[sunfluidh:sunfluidh_examples_data_set|Click here to come back to the previous page]]
sunfluidh/2d_channel_flow_with_heated_bar_incomp_flow.1480505122.txt.gz · Dernière modification : 2016/11/30 12:25 de yann
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki