=========================================================================== =========================================================================== MAIN INPUT DATA FILE : LMN APPORACH FOR THERMOACOUSTIC PROBLEMS MODELING THE ATCHLEY'S EXPERIMENT (Atchley et al., JASA 1990) DIMENSIONAL FORM symmetrical plan ------------------------------------------------------- | plate | --------------------- <---> <---> Vl.sin(wt) Vr.sin(wt) --------------------- | plate | ------------------------------------------------------- symmetrical plan D= 0.5 kx=0.35 =========================================================================== =========================================================================== &Version File_Version="VERSION2.0"/ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ GENERAL LAYOUT (DIMENSIONLESS) ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ &Fluid_Properties Variable_Density = .true. , Constant_Mass_Flow = .false. , Variable_Fluid_Volume = .false. , Molecular_Mass= 4.D-03 , Heat_Transfer_Flow = .true. , Reference_Dynamic_Viscosity = 1.985000D-05, Reference_Density = 1.833106D-01 , Reference_Temperature = 298.4 , Prandtl = 6.781648D-01 , Reference_Heat_Capacity = 5193.0 , Heat_Capacity_Ratio = 1.67 / &Velocity_Initialization I_Velocity_Reference_Value = 1.D-04 , J_Velocity_Reference_Value = 0.0 , K_Velocity_Reference_Value = 0.0 / &Temperature_Initialization Temperature_Reference_Value= 298.4 , Initial_Field_Option_For_Temperature= 0 / =========================================================================== DOMAIN FEATURES (DIMENSIONLESS) =========================================================================== &Domain_Features Start_Coordinate_I_Direction=-0.0153945 , End_Coordinate_I_Direction= 0.0153945, Start_Coordinate_J_Direction= 0.00 , End_Coordinate_J_Direction= 1.710500D-03, Start_Coordinate_K_Direction= 0.00 , End_Coordinate_K_Direction= 0.00, Cells_Number_I_Direction= 256 ,Cells_Number_J_Direction= 64 ,Cells_Number_K_Direction= 1, Regular_Mesh= .false. / ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ DEFINITION OF BOUNDARY CONDITIONS ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ============================================================================= WALL BOUNDARY CONDITION SETUP ============================================================================= &Heat_Wall_Boundary_Condition_Setup Wall_BC_DataSetName ="Set1", West_Heat_BC_Option = 0 , East_Heat_BC_Option = 0 , Back_Heat_BC_Option = 0 , Front_Heat_BC_Option = 0 , South_Heat_BC_Option = 0 , North_Heat_BC_Option = 0 / &Heat_Wall_Boundary_Condition_Setup Wall_BC_DataSetName ="Set2", West_Heat_BC_Option = 2 , East_Heat_BC_Option = 2 , Back_Heat_BC_Option = 2 , Front_Heat_BC_Option = 2 , South_Heat_BC_Option = 2 , North_Heat_BC_Option = 2, Material_Thermal_Conductivity= 7.296 , Material_Mass_Heat_Capacity = 1.67540D+03 , Material_Density = 1.0D+02 / &Polyhedral_Immersed_Bodies Xi_1=-3.425D-03, Xj_1= 0.0000 , Xk_1= 0.0 , Xi_2= 3.425D-03, Xj_2= 0.0000 , Xk_2= 0.0 , Xi_3= 3.425D-03, Xj_3= 9.525D-05 , Xk_3= 0.0 , Xi_4=-3.425D-03, Xj_4= 9.525D-05 , Xk_4= 0.0 , Wall_BC_DataSetName="Set2" / &Polyhedral_Immersed_Bodies Xi_1=-3.425D-03, Xj_1= 1.61525D-03 ,Xk_1= 0.0 , Xi_2= 3.425D-03, Xj_2= 1.61525D-03 ,Xk_2= 0.0 , Xi_3= 3.425D-03, Xj_3= 1.71050D-03 ,Xk_3= 0.0, Xi_4=-3.425D-03, Xj_4= 1.71050D-03 ,Xk_4= 0.0 , Wall_BC_DataSetName="Set2" / ============================================================================= INLET AND OUTLET BOUNDARY CONDITIONS ============================================================================= &Inlet_Boundary_Conditions Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Plan_Location_Coordinate= -0.0153945 , Start_Coordinate_of_First_Span = 0.00 , End_Coordinate_of_First_Span = 1.710500D-03 , Flow_Direction= 1 , Define_Velocity_Profile= 0, Normal_Velocity_Reference_Value= 2.282681 , Temperature_Reference_Value= 298.4 , Density_Reference_Value= 1.833106D-01 , Temperature_BC_Type="Neumann", Density_BC_Type="Neumann", Variable_Flowrate= 1 , Time_Fct_Name= "Sinus_zero_average" , Time_Fct_Threshold= 0.0 , Time_Fct_Time_Scale= 1.76173611342339D-03 , Time_Fct_Magnitude= 2.195749 / &Inlet_Boundary_Conditions Type_of_BC= "INLET", Direction_Normal_Plan= 1 , Plan_Location_Coordinate= 0.0153945 , Start_Coordinate_of_First_Span = 0.00 , End_Coordinate_of_First_Span = 1.710500D-03 , Flow_Direction=-1 , Define_Velocity_Profile= 0, Normal_Velocity_Reference_Value= 2.051303 , Temperature_Reference_Value= 298.4 , Density_Reference_Value= 1.833106D-01 , Temperature_BC_Type="Neumann", Density_BC_Type="Neumann", Variable_Flowrate= 1 , Time_Fct_Name= "Sinus_zero_average" , Time_Fct_Threshold= 0.0 , Time_Fct_Time_Scale= 1.76173611342339D-03 , Time_Fct_Magnitude= 2.195749 / ============================================================================= BORDER BOUNDARY CONDITIONS ============================================================================= &Border_Domain_Boundary_Conditions West_BC_Name = "None" , East_BC_Name = "None" , Back_BC_Name = "Symmetric" , Front_BC_Name = "Symmetric" , 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 SolverName="SOR" ,!Successive Over-Relaxation (SOR) method based on the red-black algorithm MultiGrid_Type="V_Cycle", Relaxation_Coefficient= 1.75 ,!Relaxation coefficient of the SOR method ( 1 <= Relaxation_Coefficient < 2) Number_max_Grid= 7 ,!Number of grid levels Number_max_Cycle= 5 Number_Iteration_FineToCoarseGrid=1 ,!number of SOR iterations applied on any grid level during the restriction step (before the coarsest grid computation) Number_Iteration_CoarseToFineGrid= 10 ,!number of SOR iterations applied on any grid level during the prolongation step (after the Coarsest grid computation) Number_Iteration_CoarsestGrid= 20 ,!number of SOR iterations applied on the coarsest grid ResidualNormalisation_Enabled= .true., !ConvergenceStrengthening_Enabled= .true. , Prolongation_Type="Linear" , Convergence_Criterion= 1.D-06 /!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 = 200000 , Final_Time = 7.5D-01 , TimeStep_Type = 0 , !TimeStep_Max = 1.76173611342339D-04, TimeStep_Max = 5.D-06, CFL_Max = 0.3 , Iterations_For_Timestep_Linear_Progress= 1 , Simulation_Backup_Rate = 1000 , Simulation_Checking_Rate = 101 / ============================================================================= PROBES MANAGEMENT ============================================================================== &Simulation_Management Probe_Recording_Rate = 21 / &Probe_Quantities_Enabled Temporal_Series_For_Quantity_Enabled(:)= .true. , .true., .false., .True., .false., .false. / &Probe_Location Xi= -3.420D-03, Xj= 9.5000D-05 , Xk= 0.0 / &Probe_Location Xi= 3.420D-03, Xj= 9.5000D-05 , Xk= 0.0 / &Probe_Location Xi= -3.420D-03, Xj= 1.6155D-03 , Xk= 0.0 / &Probe_Location Xi= 3.420D-03, Xj= 1.6155D-03 , Xk= 0.0 / &Probe_Location Xi= -3.420D-03, Xj= 8.5525D-04 , Xk= 0.0 / &Probe_Location Xi= 3.420D-03, Xj= 8.5525D-04 , Xk= 0.0 / ============================================================================= FIELDS RECORDING DECLARATION ============================================================================= !------------------------------ !--- Instantaneous fields !------------------------------ &Field_Recording_Setup Precision_On_Instantaneous_Fields= 1 / &Simulation_Management InstantaneousFields_TimeRecordingRate = 1.76D-04 InstantaneousFields_RecordingStartTime=0.71 / &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 = "RHO " / !Instantaneous_Fields_Listing Name_of_Field = "TRACE " / !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+05 , 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= "VolumeFlowBalance" , Filename= "check_volflux.dat" , iter_rec= 20 , Time_Start_Mean_Calculation= 2.D+05 ,Range_Mean_Calculation= 1.76D-01 Restart_Enabled= .false. / &Time_Series_Data Title= "LMN_Data" , Filename= "check_lmndata.dat" , iter_rec= 20 , Time_Start_Mean_Calculation= 2.D+05 ,Range_Mean_Calculation= 1.76D-01 Restart_Enabled= .false. / &Time_Series_Data Title= "KE_Average" , Filename= "check_kin_energy.dat" , iter_rec= 20 , Time_Start_Mean_Calculation= 2.D+05 ,Range_Mean_Calculation= 1.76D-01 Restart_Enabled= .false. /