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 [2017/12/01 11:44] – yannsunfluidh:2d_heat_driven_cavity_incomp_flow [2020/02/03 17:43] (Version actuelle) – yann
Ligne 5: Ligne 5:
  
                 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 49: 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-03  , +                   Reference_Temperature= 1.0    ,  Reference_Dynamic_Viscosity= 0.71D-02  , 
-                   Prandtl = 0.71               ,  Thermal_Expansion_Coefficient= 1.0 /+                   Reference_Heat_Capacity=  1.0 , 
 +                   Prandtl = 0.71                ,  Thermal_Expansion_Coefficient= 1.0 /
                                              
 =========================================================================== ===========================================================================
Ligne 72: Ligne 82:
                  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= 256 ,Cells_Number_J_Direction=256 ,Cells_Number_K_Direction= 1,+                 Cells_Number_I_Direction= 64 ,Cells_Number_J_Direction=64 ,Cells_Number_K_Direction= 1,
                  Regular_Mesh= .true. /                  Regular_Mesh= .true. /
  
Ligne 101: Ligne 111:
 &Numerical_Methods  NS_NumericalMethod= "BDF2-SchemeO2"                    ,       !--- BDF2 + 2nd order centered scheme &Numerical_Methods  NS_NumericalMethod= "BDF2-SchemeO2"                    ,       !--- BDF2 + 2nd order centered scheme
                     MomentumConvection_Scheme="Centered-O2-Conservative"   ,       !--- conservative form for solving the velocity (momentum) equation                     MomentumConvection_Scheme="Centered-O2-Conservative"   ,       !--- conservative form for solving the velocity (momentum) equation
-                    Poisson_NumericalMethod="Home-Multigrid-ConstantMatrixCoef"  / !--- SOR  + multigrid method (homemade release) for solving the Poisson's equation with constant coefficient matrix                    +                    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  &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)                             Relaxation_Coefficient= 1.7 ,          !--- Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2)
-                            Number_max_Grid= 5,                    !--- Number of grid levels +                            Number_max_Grid= 6,                    !--- Number of grid levels 
-                            Number_max_Cycle= 10,                  !--- Number of multigrid cycles +                            Number_max_Cycle= 1,                  !--- 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= 3, !--- number of SOR iterations applied on any grid level during the restriction step (before the coarsest grid computation)
-                            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_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                             Number_Iteration_CoarsestGrid= 15 ,    !--- number of SOR iterations applied on the coarsest grid
-                            Convergence_Criterion= 1.D-10 /        !--- convergence tolerance on the residu of the Poisson's equation                   +                            Convergence_Criterion= 1.D-08 /        !--- convergence tolerance on the residu of the Poisson's equation                   
  
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Ligne 121: Ligne 130:
 &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 
 +    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 &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     " , Recording_Enabled = .true.  /      First velocity component +&Instantaneous_Fields_Listing  Name_of_Field = "U     "  /      First velocity component 
-&Instantaneous_Fields_Listing  Name_of_Field = "V     " , Recording_Enabled = .true.  /      Second velocity component +&Instantaneous_Fields_Listing  Name_of_Field = "V     "  /      Second velocity component 
-&Instantaneous_Fields_Listing  Name_of_Field = "T     " , Recording_Enabled = .true.  /      Temperature+&Instantaneous_Fields_Listing  Name_of_Field = "T     "  /      Temperature
  
 END OF FILE END OF FILE
sunfluidh/2d_heat_driven_cavity_incomp_flow.1512125070.txt.gz · Dernière modification : 2017/12/01 11:44 de yann
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki