=========================================================================== =========================================================================== MAIN INPUT DATA FILE : SUCCESSIVE COMPRESSION CYCLES OF A 2D QUASI-ADIABATIC CAVITY LMN APPORACH Wall Qh=0 ------------------------------------------------------------| | | | | <-- | --> | Wall Qh= 0 T= T0 | | | | ------------------------------------------------------------| Symmetric plan DIMENSIONAL SETUP fluid : air at P= 101325 Pa and T= 294.146 K (initial condition at t=0) cavity size= 0.10m * 0.015m volume variation : A.sin(2pi.f.t- w0) w0= pi/2 A= 0.01 m f= 50 Hz =========================================================================== =========================================================================== &Version File_Version="VERSION2.0"/ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ GENERAL LAYOUT (DIMENSIONLESS) ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ &Fluid_Properties Variable_Density = .true. , Constant_Mass_Flow = .false. , Variable_Fluid_Volume= .true. , Molecular_Mass= 0.02896421357024561 , Heat_Transfer_Flow = .true. , Reference_Dynamic_Viscosity = 1.795D-05, Reference_Density = 1.2 , Reference_Temperature = 294.146 , Prandtl = 0.726 , Reference_Heat_Capacity = 1004.7093960142244 , Heat_Capacity_Ratio = 1.4 / &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= 294.146 , Initial_Field_Option_For_Temperature= 0 / =========================================================================== DOMAIN FEATURES (DIMENSIONLESS) =========================================================================== &Domain_Features Start_Coordinate_I_Direction=-0.020 , End_Coordinate_I_Direction= 0.10, Start_Coordinate_J_Direction= 0.00 , End_Coordinate_J_Direction= 0.015, Start_Coordinate_K_Direction= 0.00 , End_Coordinate_K_Direction= 0.00, Cells_Number_I_Direction= 512 ,Cells_Number_J_Direction= 64 ,Cells_Number_K_Direction= 1, Regular_Mesh= .false. / ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ DEFINITION OF BOUNDARY CONDITIONS ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ============================================================================= WALL BOUNDARY CONDITION SETUP (DIMENSIONLESS) ============================================================================= ============================================================================= INLET AND OUTLET BOUNDARY CONDITIONS (DIMENSIONLESS) ============================================================================= &Inlet_Boundary_Conditions Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Plan_Location_Coordinate= -0.0200 , Start_Coordinate_of_First_Span = 0.00 , End_Coordinate_of_First_Span = 0.015 , Start_Coordinate_of_Second_Span= 0.0 , End_Coordinate_of_Second_Span= 0.0 , Flow_Direction= 1 , Define_Velocity_Profile= 0, Normal_Velocity_Reference_Value= 0.0 , Temperature_Reference_Value= 294.146 , Density_Reference_Value= 1.2 , Temperature_BC_Type="Neumann", Density_BC_Type="Neumann", Variable_Flowrate= 0/ !Time_Fct_Name= "Sinus_zero_average" , Time_Fct_Threshold= 0.0 , !Time_Fct_Time_Scale= 2.0E-02 , Time_Fct_Magnitude= 4.712 / ============================================================================= IMMERSED BOUNDARY METHODS : PENALIZATION METHOD ============================================================================= &Immersed_Boundary_Methods Immersed_Boundary_Method_Type= 1 , Name_of_Solid="RECTANGLE", Linked_IBM_Inlet_Number = 1 , Gravity_IBM_Enabled = .false. , StrongPenaltyForVelocity_Enabled=.true., Coordinate_Gravity_Center_At_Rest_I=-0.04 , Coordinate_Gravity_Center_At_Rest_J= 0.0075 , Coordinate_Gravity_Center_At_Rest_K= 0.0 , Coordinate_Gravity_Center_Init_I= -0.05 , Coordinate_Gravity_Center_Init_J= 0.0075 , Coordinate_Gravity_Center_Init_K= 0.0 , Coordinate_Restricted_Motion_I=-0.04 , Coordinate_Restricted_Motion_J= 0.0 , Coordinate_Restricted_Motion_K= 0.0 , Size_Object_I= 0.08 , Size_Object_J= 0.015 , Size_Object_K= 0.0 , Reference_Velocity_I= 0.0 , Reference_Velocity_J= 0.0 , Reference_Velocity_K= 0.0 , Motion_Magnitude_I= 0.010 , Motion_Magnitude_J= 0.00 , Motion_Magnitude_K= 0.0 , Mass= 1.2 , Spring_Stiffness_Constant_I= 0.0 , Spring_Stiffness_Constant_J= 0.0 , Spring_Stiffness_Constant_K= 0.0 , Transitional_Time= 0.00, Forced_Frequency= 5.0D+01 , OffSet= 0.5, Heat_Transfer_Type= 0, Reference_Temperature= 294.146 , Reference_Heat_Flux= 0.0 , Material_Thermal_Conductivity= 2.4841D-02 , Material_Mass_Heat_Capacity= 1004.7094, Material_Density= 1.2 / ============================================================================= ENDS BOUNDARY CONDITIONS ============================================================================= &Border_Domain_Boundary_Conditions West_BC_Name = "None" , East_BC_Name = "None" , Back_BC_Name = "Symmetric" , Front_BC_Name = "None" , North_BC_Name = "None" , South_BC_Name = "None" / ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ NUMERICAL METHODS ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ &Numerical_Methods !!!NS_NumericalMethod= "BDF2-SchemeO2", NS_NumericalMethod= "CN-SchemeO2", !!!NS_NumericalMethod= "CN-SchemeO2-SpecialLowMachFlow", !!MomentumConvection_Scheme="Centered-O2-Conservative" , MomentumConvection_Scheme="Centered-O2-Convective_2" , !!!MomentumConvection_Scheme="Upwind-O1-Convective" , TemperatureAdvection_Scheme="Centered-O2-Convective_2" , !!!TemperatureAdvection_Scheme="Upwind-O1-Convective" , !!!TemperatureAdvection_Scheme="Centered-O2-Convective" , !Poisson_NumericalMethod="Home-PartialDiagonalization" , !Poisson_NumericalMethod="Home-SORMultigrid-ConstantMatrixCoef" , Poisson_NumericalMethod="Home-SORMultigrid-VariableMatrixCoef" , Off_Set_Poisson_Source_term= .true./ !HomeData_PoissonSolver Direction_1= 2, Direction_2= 1, Direction_3= 3/ &HomeData_PoissonSolver SolverName="SOR" ,!Successive Over-Relaxation (SOR) method based on the red-black algorithm MultiGrid_Type="V_Cycle", Relaxation_Coefficient= 1.45 ,!Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2) Number_max_Grid= 7 ,!Number of grid levels Number_max_Cycle= 9 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= 30 ,!number of SOR iterations applied on the coarsest grid ResidualNormalisation_Enabled= .true., ConvergenceStrengthening_Enabled= .true. , Convergence_Criterion= 1.D-04 /!convergence tolerance on the residu of the Poisson's equation ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ SIMULATION MANAGEMENT ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ The numerical time step is dynamic and is estimated by the constant CFL coefficient &Simulation_Management Restart_Parameter= 0 , Steady_Flow_Stopping_Criterion_Enabled = .true. , Steady_Flow_Stopping_Criterion = 1.D-16, Temporal_Iterations_Number = 4000000 , Final_Time = 4.00D-00 , TimeStep_Type = 0 , TimeStep_Max = 1.0D-05, CFL_Max = 0.3 , Simulation_Backup_Rate = 1000 , Simulation_Checking_Rate = 101 / ============================================================================= PROBES MANAGEMENT ============================================================================== &Simulation_Management Probe_Recording_Rate = 100 / &Probe_Quantities_Enabled Temporal_Series_For_Quantity_Enabled(:)= .true. , .true., .false., .false., .false., .false. / &Probe_Location Xi=-0.011 , Xj= 0.007 , Xk= 0.0 / ============================================================================= FIELDS RECORDING DECLARATION ============================================================================= !------------------------------ !--- Instantaneous fields !------------------------------ &Field_Recording_Setup Precision_On_Instantaneous_Fields= 1 / &Simulation_Management InstantaneousFields_TimeRecordingRate = 1.0D-03 InstantaneousFields_RecordingStartTime=0.0 / &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 " / &Instantaneous_Fields_Listing Name_of_Field = "T " / &Instantaneous_Fields_Listing Name_of_Field = "TRACE " / &Instantaneous_Fields_Listing Name_of_Field = "RHO " / !Instantaneous_Fields_Listing Name_of_Field = "divU " / !Instantaneous_Fields_Listing Name_of_Field = "divRU " / !Instantaneous_Fields_Listing Name_of_Field = "PHI " / !Instantaneous_Fields_Listing Name_of_Field = "drho " / !Instantaneous_Fields_Listing Name_of_Field = "SRC_P " / !------------------------------ !--- Statistic fields !------------------------------ &Simulation_Management Start_Time_For_Statistics= 2.D-01 , Time_Range_Statistic_Calculation = 3.8D-01 / &Statistical_Fields_Listing Name_of_Field = " " / &Statistical_Fields_Listing Name_of_Field = " " / &Statistical_Fields_Listing Name_of_Field = " " / &Statistical_Fields_Listing Name_of_Field = "

" / &Statistical_Fields_Listing Name_of_Field = " " / !------------------------------ !--- Time series !------------------------------ &Time_Series_Data Title= "MassFlowBalance" , Filename= "check_mass.dat" , iter_rec= 10 , Time_Start_Mean_Calculation= 40. ,Range_Mean_Calculation= 80. Restart_Enabled= .false. / &Time_Series_Data Title= "VolumeFlowBalance" , Filename= "check_vol.dat" , iter_rec= 10 , Time_Start_Mean_Calculation= 40. ,Range_Mean_Calculation= 80. Restart_Enabled= .false. / &Time_Series_Data Title= "LMN_Data" , Filename= "check_lmndata.dat" , iter_rec= 10 , Time_Start_Mean_Calculation= 40. ,Range_Mean_Calculation= 80. Restart_Enabled= .false. / &Time_Series_Data Title= "KE_Average" , Filename= "check_kin_energy.dat" , iter_rec= 10 , Time_Start_Mean_Calculation= 40. ,Range_Mean_Calculation= 80. Restart_Enabled= .false. /