Documentation du code de simulation numérique SUNFLUIDH

Outils pour utilisateurs

Outils du site

Piste :
sunfluidh:2d_eddy_convection_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_eddy_convection_incomp_flow [2016/12/05 12:04] – witkosunfluidh:2d_eddy_convection_incomp_flow [2019/12/10 18:13] (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 10: Ligne 10:
  
 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 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+                                     access by the data file --->  Initial_Field_Option_For_Velocity_I
  
 =========================================================================== ===========================================================================
Ligne 30: Ligne 30:
                                       in the module named module_user_define_init_fields.f90                                       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 /+&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= 5 /
 =========================================================================== ===========================================================================
                 GRAVITY                 GRAVITY
Ligne 65: Ligne 65:
 The walls at the ends of the domain are here replaced by periodical conditions The walls at the ends of the domain are here replaced by periodical conditions
  
-&Border_Domain_Boundary_Conditions West_Border= 4 , East_Border= 4 , Back_Border= 4 , Front_Border= 4 , North_Border = 0 , South_Border = 0 /+&Border_Domain_Boundary_Conditions West_BC_Name= "Periodic" , East_BC_Name= "Periodic" , Back_BC_Name= "Periodic" , Front_BC_Name= "Periodic" /
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
                    NUMERICAL METHODS                    NUMERICAL METHODS
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
- Default initialization:+&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                    
  
- Numerical scheme for Navier-Stokes equations : 2nd centrered order in space, 2nd order in time (BDF2 : 2nd oder Backward Differentiation Formulation) +&HomeData_PoissonSolver  SolverName="SOR"                ,         !--- Successive Over-Relaxation (SOR) method 
-                  for Poisson equation        : SOR iterative method coupled with a multigrid approach +                            Relaxation_Coefficient= 1.8  ,         !--- Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2) 
- +                            Number_max_Grid= 6,                    !--- Number of grid levels 
-Parameters for solving Poisson equation (multigrid parameters) : +                            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 
-&Numerical_Methods  Numerical_Scheme= 1 , +                            Number_Iteration_FineToCoarseGrid= 15, !--- number of SOR iterations applied on any grid level during the restriction step (before the coarsest grid computation) 
-                    Convective_Flux_Discretization_Type = 1 +                            Number_Iteration_CoarseToFineGrid= 15, !--- number of SOR iterations applied on any grid level during the prolongation step (after the Coarsest grid computation) 
-                    Numerical_Method_Poisson_Equation= 1, +                            Number_Iteration_CoarsestGrid= 15 ,    !--- number of SOR iterations applied on the coarsest grid 
-                    Number_max_Grid= 6 , Number_max_Cycle= 5 , Number_Iteration= 15 , +                            Convergence_Criterion= 1.D-10 /        !--- convergence tolerance on the residu of the Poisson's equation   
-                    Relaxation_Coefficient = 1.80 , Convergence_Criterion = 1.D-10 /+
  
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Ligne 89: Ligne 90:
  
 &Simulation_Management    Restart_Parameter= 0 , &Simulation_Management    Restart_Parameter= 0 ,
-                          Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-22,+                          Steady_Flow_Stopping_Criterion_Enabled = .false. , Steady_Flow_Stopping_Criterion = 1.D-22,
                           Temporal_Iterations_Number = 50000                  , Final_Time = 2.0D+01  ,                           Temporal_Iterations_Number = 50000                  , Final_Time = 2.0D+01  ,
                           TimeStep_Type = 1 ,                           TimeStep_Type = 1 ,
                           CFL_min= 0.25   , CFL_max= 0.25, Iterations_For_Timestep_Linear_Progress= 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 ,+                          Timestep_Min = 1.00                           , Timestep_Max = 1.00 ,
                           Simulation_Backup_Rate                 =  200   , Simulation_Checking_Rate = 51  /                           Simulation_Backup_Rate                 =  200   , Simulation_Checking_Rate = 51  /
  
 ============================================================================= =============================================================================
-                   PROBES MANAGEMENT+                   PROBES SETUP
 ============================================================================= =============================================================================
                                                       Probes order    U    , V     , W      , T      , P      , RHO                                                       Probes order    U    , V     , W      , T      , P      , RHO
Ligne 103: Ligne 104:
  
 &Probe_Location  Xi= 0.0 , Xj= 0.0 , Xk= 0.0  / &Probe_Location  Xi= 0.0 , Xj= 0.0 , Xk= 0.0  /
-&Probe_Location  Xi= 0.1 , Xj= 0.1 , Xk= 0.0 , End_of_Data_Block= .true. / +&Probe_Location  Xi= 0.1 , Xj= 0.1 , Xk= 0.0  /
-&Simulation_Management   Probe_Recording_Rate = 1 /+
  
 +&Simulation_Management  Probe_TimeIterationRecordingRate= 1   , 
 +                        Probe_StartTimeIterationRecording= 0  ,
 +                        Probe_RecordingReset=.false. /
  
 ============================================================================= =============================================================================
-           FIELDS RECORDING DECLARATION+              INSTANTANEOUS FIELD 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 /+&Simulation_Management 
 +    InstantaneousFields_RecordingReset=.false.     ,        
 +    InstantaneousFields_TimeRecordingRate= 1.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 
 +&Instantaneous_Fields_Listing  Name_of_Field = "P" /   Pressure related to the velocity field
  
-&Instantaneous_Fields_Listing  Name_of_Field = "U" , Recording_Enabled = .true.  /      First velocity component +END OF FILE
-&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> </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_eddy_convection_incomp_flow.1480935876.txt.gz · Dernière modification : 2016/12/05 12:04 de witko
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki