Documentation du code de simulation numérique SUNFLUIDH

Outils pour utilisateurs

Outils du site

Piste :
sunfluidh:2d_heat_driven_cavity_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_heat_driven_cavity_incomp_flow [2016/11/30 10:51] – yannsunfluidh:2d_heat_driven_cavity_incomp_flow [2020/02/03 17:43] (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 HEAT-DRIVEN CAVITY FLOW MAIN INPUT DATA FILE : 2D HEAT-DRIVEN CAVITY FLOW
  
                 DIMENSIONLESS FORM  :                 DIMENSIONLESS FORM  :
-                                       Ra= 1.D+06+                                       Ra= 1.D+04
                                        Density scale   : rho_0 (fluid density)                                        Density scale   : rho_0 (fluid density)
                                        Length scale  H : height of cavity                                        Length scale  H : height of cavity
-                                       Velocity scale ---> U0= (k/H).Ra**0.25 (k thermal diffusivity)+                                       Velocity scale ---> U0= (k/H).Ra**0.5 (k thermal diffusivity)
                                        Temperature scale Th - Tc (Th= T hot ; Tc= T cold)                                        Temperature scale Th - Tc (Th= T hot ; Tc= T cold)
                                        Dimensionless Velocity    U*= U/U0                                        Dimensionless Velocity    U*= U/U0
Ligne 16: Ligne 16:
                                        dimensionless buoyancy term      = Pr.T*                                        dimensionless buoyancy term      = Pr.T*
                                        dimensionless domain Lx/H= 1 , Ly/H= 1                                        dimensionless domain Lx/H= 1 , Ly/H= 1
 +                                       
 +                                       Reference results (De Wahl Davis, IJNMF , Vol 3, 1983): 
 +                                          Averaged Nusselt number at the wall               : 2.243
 +                                          Averaged Nusselt number at the vertical mid-plane : 2.243
 +                                          Maximum Nusselt number at the wall                : 3.528
 +                                          Minimum Nusselt number at the wall                : 0.586
 +                                          Maximum value of the vertical velocity (horizontal mid-plane) : 0.19617
 +                                          Maximum value of the horizontal velocity (vertital mid-plane) : 0.16178
 +                                          
  
                  Q0= 0 (adiabatic)                  Q0= 0 (adiabatic)
Ligne 43: Ligne 52:
 GENERAL LAYOUT GENERAL LAYOUT
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
- +&Version File_Version="VERSION2.0"/
-   &Version File_Version="VERSION2.0"/+
        
 =========================================================================== ===========================================================================
Ligne 50: Ligne 58:
                (DIMENSIONLESS FORM)                (DIMENSIONLESS FORM)
 =========================================================================== ===========================================================================
- +&Fluid_Properties  Heat_Transfer_Flow = .true.   ,  Reference_Density= 1.0, 
-   &Fluid_Properties  Heat_Transfer_Flow = .true.  ,  Reference_Density= 1.0, +                   Reference_Temperature= 1.0    ,  Reference_Dynamic_Viscosity= 0.71D-02  , 
-                      Reference_Temperature= 1.0   ,  Reference_Dynamic_Viscosity= 0.71D-03  , +                   Reference_Heat_Capacity=  1.0 , 
-                      Prandtl = 0.71               ,  Thermal_Expansion_Coefficient= 1.0 /+                   Prandtl = 0.71                ,  Thermal_Expansion_Coefficient= 1.0 /
                                              
 =========================================================================== ===========================================================================
Ligne 59: Ligne 67:
                   (DIMENSIONLESS FORM)                   (DIMENSIONLESS FORM)
 =========================================================================== ===========================================================================
- +&Velocity_Initialization I_Velocity_Reference_Value          = 0.0  , J_Velocity_Reference_Value          = 0.0  , K_Velocity_Reference_Value          = 0.0 / 
-   &Velocity_Initialization I_Velocity_Reference_Value          = 0.0  , J_Velocity_Reference_Value          = 0.0  , K_Velocity_Reference_Value          = 0.0 / +&Temperature_Initialization Temperature_Reference_Value      = 0.5 /
-    +
-   &Temperature_Initialization Temperature_Reference_Value      = 0.5 /+
        
 =========================================================================== ===========================================================================
Ligne 68: Ligne 74:
                   (DIMENSIONLESS FORM)                   (DIMENSIONLESS FORM)
 =========================================================================== ===========================================================================
- +&Gravity  Gravity_Enabled= .true. , Gravity_Angle_IJ= 90.0  , Gravity_Angle_IK= 90.0 , Reference_Gravity_Constant= 0.71/
-   &Gravity  Gravity_Enabled= .true. , Gravity_Angle_IJ= 90.0  , Gravity_Angle_IK= 90.0 , Reference_Gravity_Constant= 0.71/+
        
 =========================================================================== ===========================================================================
                     DOMAIN FEATURES                     DOMAIN FEATURES
 =========================================================================== ===========================================================================
- +&Domain_Features Start_Coordinate_I_Direction= 0.00 , End_Coordinate_I_Direction= 1.00, 
-   &Domain_Features Start_Coordinate_I_Direction= 0.00 , End_Coordinate_I_Direction= 1.00, +                 Start_Coordinate_J_Direction= 0.00 , End_Coordinate_J_Direction= 1.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, 
-                    Start_Coordinate_K_Direction= 0.00 , End_Coordinate_K_Direction= 0.00, +                 Cells_Number_I_Direction= 64 ,Cells_Number_J_Direction=64 ,Cells_Number_K_Direction= 1, 
-                    Cells_Number_I_Direction= 256 ,Cells_Number_J_Direction=256 ,Cells_Number_K_Direction= 1, +                 Regular_Mesh= .true. /
-                    Regular_Mesh= .true. / +
-                    +
  
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Ligne 88: Ligne 91:
 WALL BOUNDARY CONDITION SETUP WALL BOUNDARY CONDITION SETUP
 ============================================================================= =============================================================================
 +Keep in mind that the domain is enclosed by default.
 +No new boundary conditions are defined at the ends of the domain : the walls by default are preserved
  
-   &Heat_Wall_Boundary_Condition_Setup +&Heat_Wall_Boundary_Condition_Setup 
-      West_Heat_BC_Option = 0   ,  East_Heat_BC_Option = 0  ,  Back_Heat_BC_Option = 1  ,  Front_Heat_BC_Option = 1  , South_Heat_BC_Option = 0 , North_Heat_BC_Option = 0, +      Wall_BC_DataSetName ="Set1", 
-      West_Wall_BC_Value= 1.0   ,  East_Wall_BC_Value= 0.0  ,  Back_Wall_BC_Value= 0.0  ,  Front_Wall_BC_Value= 0.0  , South_Wall_BC_Value= 0.0 , North_Wall_BC_Value= 0.0 , +      West_Heat_BC_Option = 0    ,  East_Heat_BC_Option = 0  ,  Back_Heat_BC_Option = 1  ,  Front_Heat_BC_Option = 1  , South_Heat_BC_Option = 0 , North_Heat_BC_Option = 0, 
-      End_of_Data_Block= .true.   /+      West_Wall_BC_Value= 1.0    ,  East_Wall_BC_Value= 0.0  ,  Back_Wall_BC_Value= 0.0  ,  Front_Wall_BC_Value= 0.0  , South_Wall_BC_Value= 0.0 , North_Wall_BC_Value= 0.0 /
              
 ============================================================================= =============================================================================
-BORDER BOUNDARY CONDITIONS+BORDER BOUNDARY CONDITIONS : The walls located by default at the ends of the domain remain unchanged
 ============================================================================= =============================================================================
 +!--- 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
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-Navier-Stokes equations ---> BDF2 + 2nd order centered scheme , semi-implicit  +&Numerical_Methods  NS_NumericalMethod= "BDF2-SchemeO2"                    ,       !--- BDF2 + 2nd order centered scheme 
-convection flux         ---> conservative form +                    MomentumConvection_Scheme="Centered-O2-Conservative"   ,       !--- conservative form for solving the velocity (momentum) equation 
-Poisson                 ---> Relaxed Gauss-Seidel method + Multigrid +                    Poisson_NumericalMethod="Home-SORMultigrid-ConstantMatrixCoef"  / !--- SOR  + multigrid method (homemade release) for solving the Poisson's equation with constant coefficient matrix                    
  
-   &Numerical_Methods  Numerical_Scheme= 1 , +&HomeData_PoissonSolver  SolverName="SOR"               ,          !--- Successive Over-Relaxation (SOR) method based on the red-black algorithm  
-                    Convective_Flux_Discretization_Type = 1    , Temperature_Advective_Flux_Discretization_Type = 1 , +                            Relaxation_Coefficient= 1.7 ,          !--- Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2) 
-                    Explicit_Solving_of_Density         = 0    , Velocity_Correction_Enabled        = .true., +                            Number_max_Grid= 6,                    !--- Number of grid levels 
-                    Numerical_Method_Poisson_Equation   = 1 ,  Number_max_Grid= 5, Number_max_Cycle= 10 , Number_Iteration= 15 , +                            Number_max_Cycle= 1,                  !--- Number of multigrid cycles 
-                    Relaxation_Coefficient= 1.7 , Convergence_Criterion= 1.D-10 /+                            Number_Iteration_FineToCoarseGrid= 3, !--- 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                   
  
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Ligne 119: Ligne 128:
 The numerical time step is estimated by means of the CFL coefficient The numerical time step is estimated by means of the 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 = 1000000              , Final_Time = 5.D+02  , +                          Temporal_Iterations_Number = 10000              , Final_Time = 5.D+02  , 
-                             TimeStep_Type = 0 , +                          TimeStep_Type = 1 , 
-                             Timestep_Min = 1.D-03                           , Timestep_Max = 1.D-03 ,+                             Timestep_Min = 5.D-02                           , Timestep_Max = 5.D-02 ,
                              CFL_Min      = 0.5                              , CFL_Max      = 0.5 ,                              CFL_Min      = 0.5                              , CFL_Max      = 0.5 ,
                             Iterations_For_Timestep_Linear_Progress= 1,                             Iterations_For_Timestep_Linear_Progress= 1,
Ligne 135: Ligne 144:
            
 ============================================================================= =============================================================================
-FIELDS RECORDING DECLARATION+FIELDS RECORDING SETUP
 ============================================================================= =============================================================================
-   &Simulation_Management  Fields_Recording_Rate = 5.D+01 / +&Simulation_Management 
-   &Field_Recording_Setup     Check_Special_Features= "Heat_Driven_Cavity_Flow",  Precision_On_Instantaneous_Fields= 2 /   Here, a special variable devoted to results of heat driven cavity flows is active+    InstantaneousFields_RecordingReset=.false.     ,        
 +    InstantaneousFields_TimeRecordingRate= 5.0E+00 , 
 +    InstantaneousFields_RecordingStartTime= 0.D-00  / 
 +&Field_Recording_Setup     Check_Special_Features= "Heat_Driven_Cavity_Flow",  Precision_On_Instantaneous_Fields= 2 /   Here, a special variable devoted to results of heat driven cavity flows is active 
 + 
 +&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 
 + 
 +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 = "T     " , Recording_Enabled = .true. , End_of_Data_Block= .true. /      Temperature 
 </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_heat_driven_cavity_incomp_flow.1480499489.txt.gz · Dernière modification : 2016/11/30 10:51 de yann
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki