Documentation du code de simulation numérique SUNFLUIDH

Outils pour utilisateurs

Outils du site

Piste :
sunfluidh:2d_channel_flow_with_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_bar_incomp_flow [2017/10/02 16:35] – yannsunfluidh:2d_channel_flow_with_bar_incomp_flow [2019/11/13 16:33] (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>
 =========================================================================== ===========================================================================
 =========================================================================== ===========================================================================
-            MAIN INPUT DATA FILE : 2D CHANNEL FLOW WITH A CONSTRICTION (A SQUARE BAR)+            MAIN INPUT DATA FILE : 2D CHANNEL FLOW WITH A CONSTRICTION (A RECTANGULAR BAR)
                                    INCOMPRESSIBLE FLOW                                    INCOMPRESSIBLE FLOW
                                    ISOTHERM                                    ISOTHERM
Ligne 19: Ligne 19:
  
                  Initialisation = uniform velocity field                  Initialisation = uniform velocity field
-                 inlet flowrate = uniform profil+                 inlet flowrate = parabolic profil
  
 INCOMPRESSIBLE DOWNSTEP FLOW INCOMPRESSIBLE DOWNSTEP FLOW
Ligne 25: Ligne 25:
  
 ------------------------------------------ ------------------------------------------
-                   |  | +                 |       | 
-                   |__|+                 |_______|
 ---> inflow                     outflow ---> ---> inflow                     outflow --->
  
Ligne 54: Ligne 54:
                  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 unique) immersed body
 +
 +&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= "Set1"/
 +                             
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
               DEFINITION OF BOUNDARY CONDITIONS               DEFINITION OF BOUNDARY CONDITIONS
Ligne 65: Ligne 76:
 --------------------------------------------------- ---------------------------------------------------
 first set of wall boundary condition first set of wall boundary condition
-(This set corresponds to the default wall boundary conditions for the velocity. It is just shown for example and could be removed)+(This set corresponds to the default wall boundary conditions for the velocity. It is just shown as an example and could be removed)
  
 &Velocity_Wall_Boundary_Condition_Setup &Velocity_Wall_Boundary_Condition_Setup
Ligne 72: Ligne 83:
  West_Wall_Velocity_J= 0.0   ,   East_Wall_Velocity_J= 0.0   ,   Back_Wall_Velocity_J= 0.0    , Front_Wall_Velocity_J= 0.0    ,  South_Wall_Velocity_J= 0.0    , North_Wall_Velocity_J= 0.0 ,  West_Wall_Velocity_J= 0.0   ,   East_Wall_Velocity_J= 0.0   ,   Back_Wall_Velocity_J= 0.0    , Front_Wall_Velocity_J= 0.0    ,  South_Wall_Velocity_J= 0.0    , North_Wall_Velocity_J= 0.0 ,
  West_Wall_Velocity_K= 0.0   ,   East_Wall_Velocity_K= 0.0   ,   Back_Wall_Velocity_K= 0.0    , Front_Wall_Velocity_K= 0.0    ,  South_Wall_Velocity_K= 0.0    , North_Wall_Velocity_K= 0.0 /  West_Wall_Velocity_K= 0.0   ,   East_Wall_Velocity_K= 0.0   ,   Back_Wall_Velocity_K= 0.0    , Front_Wall_Velocity_K= 0.0    ,  South_Wall_Velocity_K= 0.0    , North_Wall_Velocity_K= 0.0 /
- 
- 
---------------------------------------------------- 
-     GEOMETRY OF THE IMMERSED BODIES 
---------------------------------------------------- 
- 
-First (and unic) immersed body 
- 
-&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= "Set1"/ 
                                                            
 ============================================================================= =============================================================================
Ligne 94: Ligne 94:
  
 &Inlet_Boundary_Conditions  Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 , &Inlet_Boundary_Conditions  Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Flow_Direction= 1 ,
 +                            Define_Velocity_profile = 1 ,
                             Plan_Location_Coordinate=  0.0    ,                             Plan_Location_Coordinate=  0.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 ,
Ligne 113: Ligne 114:
                    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_CoarsestGrid= 15 ,    !--- number of SOR iterations applied on the coarsest grid 
 +                            Number_Iteration_CoarseToFineGrid= 15, !--- number of SOR iterations applied on any grid level during the prolongation step (after the Coarsest grid computation) 
 +                            Convergence_Criterion= 1.D-08 /        !--- convergence tolerance on the residu of the Poisson's equation                   
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
                    SIMULATION MANAGEMENT                    SIMULATION MANAGEMENT
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-The numerical time step is imposed+The numerical time step is dynamic and is estimated by the constant CFL coefficient
  
 &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+01  , 
-                          TimeStep_Type = 0 , +                          TimeStep_Type = 1 , 
-                          Timestep_Max = 1.D-03 ,+                          CFL_min= 0.25   , CFL_max= 0.25, Iterations_For_Timestep_Linear_Progress= 1 , 
 +                          Timestep_Max = 1.D+03 ,
                           Simulation_Backup_Rate                 = 1000   , Simulation_Checking_Rate = 101 /                           Simulation_Backup_Rate                 = 1000   , Simulation_Checking_Rate = 101 /
 ============================================================================= =============================================================================
Ligne 143: Ligne 148:
  
 &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 153: Ligne 160:
 !--- 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.   /      Second velocity component+    InstantaneousFields_TimeRecordingRate= 5.0E+00 , 
 +    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
  
 !--- Statistics !--- Statistics
Ligne 161: Ligne 172:
 &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.  /+&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_bar_incomp_flow.1506954915.txt.gz · Dernière modification : 2017/10/02 16:35 de yann
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki