w2run_comment = "Laser Slab Example Simulation" log_file = "lasslab.log" basenm = "lasslab_" # This particular parfile is used as an example that is described in # detail in the users guide. ########################## # # # I/O PARAMETERS # # # ########################## ### Checkpoint Options ### checkpointFileIntervalTime = 10.e-12 #2.5e-12 #checkpointFileIntervalStep = 100 ### Plot Options ### plotFileNumber = 0 #plotFileIntervalStep = 100 plotFileIntervalTime = 10.e-12 #2.5e-12 #10.e-12 plot_var_1 = "dens" plot_var_2 = "depo" plot_var_3 = "tele" plot_var_4 = "tion" plot_var_5 = "trad" plot_var_6 = "ye " plot_var_7 = "sumy" plot_var_8 = "cham" plot_var_9 = "targ" plot_var_10 = "fpos" plot_var_11 = "fvel" plot_var_12 = "phas" ### Restart Options ### restart = .false. #checkpointFileNumber = 27 #plotFileNumber = 28 #./setup -auto LaserSlab_FOAM -2d +cylindrical -nxb=100 -nyb=250 +hdf5typeio species=cham,targ +mtmmmt +laser +uhd3t -unit=physics/sourceTerms/Cool/CoolMain/Foam_hybrid3 ######################################## # # # FOAM SUBSCALE PARAMETERS # # # ######################################## useFoam = .false. density_cylinder = 0. length_pore = 0.e-4 wavelength = 0. ######################################## # # # RADIATION/OPACITY PARAMETERS # # # ######################################## rt_useMGD = .false. rt_mgdNumGroups = 20 rt_mgdBounds_1 = 1.0e-01 rt_mgdBounds_2 = 1.0e+00 rt_mgdBounds_3 = 5.0e+00 rt_mgdBounds_4 = 1.0e+01 rt_mgdBounds_5 = 2.0e+01 rt_mgdBounds_6 = 4.0e+01 rt_mgdBounds_7 = 1.0e+02 rt_mgdBounds_8 = 2.0e+02 rt_mgdBounds_9 = 3.2e+02 rt_mgdBounds_10 = 4.0e+02 rt_mgdBounds_11 = 6.8e+02 rt_mgdBounds_12 = 9.0e+02 rt_mgdBounds_13 = 1.1e+03 rt_mgdBounds_14 = 1.55e+03 rt_mgdBounds_15 = 2.5e+03 rt_mgdBounds_16 = 4.0e+03 rt_mgdBounds_17 = 6.5e+03 rt_mgdBounds_18 = 1.065e+04 rt_mgdBounds_19 = 1.7e+04 rt_mgdBounds_20 = 2.5e+04 rt_mgdBounds_21 = 4.0e+04 rt_mgdFlMode = "fl_harmonic" rt_mgdFlCoef = 1.0 rt_mgdXlBoundaryType = "vacuum" #"reflecting" rt_mgdXrBoundaryType = "vacuum" rt_mgdYlBoundaryType = "vacuum" rt_mgdYrBoundaryType = "reflecting" rt_mgdZlBoundaryType = "reflecting" rt_mgdZrBoundaryType = "reflecting" useOpacity = .false. ### SET CHAMBER (HELIUM) OPACITY OPTIONS ### op_chamAbsorb = "op_tabpa" op_chamEmiss = "op_tabpe" op_chamTrans = "op_tabro" op_chamFileType = "ionmix4" op_chamFileName = "he-imx-20g.cn4" #"he-imx-005.cn4"#"he-imx-20g.cn4" #"ar-imx-20g.cn4" ### SET TARGET (FOAM) OPACITY OPTIONS ### op_targAbsorb = "op_tabpa" op_targEmiss = "op_tabpe" op_targTrans = "op_tabro" op_targFileType = "ionmix4" op_targFileName = "he-imx-20g.cn4" #"al-imx-003.cn4"# "cu-imx-20g.cn4" ############################ # # # LASER PARAMETERS # # # ############################ useEnergyDeposition = .true. ed_maxRayCount = 200000 ed_gradOrder = 2 # Activate 3D-in-2D ray trace: ed_laser3Din2D = .true. ed_laser3Din2DwedgeAngle = 0.1 ### LASER IO OPTIONS ### ed_useLaserIO = .true. ed_laserIOMaxNumberOfPositions = 10000 ed_laserIOMaxNumberOfRays = 128 ### SETUP LASER PULSES ### ed_numberOfPulses = 1 # Define Pulse 1: ed_numberOfSections_1 = 10 ed_time_1_1 = 0.00000000 ed_time_1_2 = 1.11111106E-10 ed_time_1_3 = 2.22222213E-10 ed_time_1_4 = 3.33333305E-10 ed_time_1_5 = 4.44444426E-10 ed_time_1_6 = 5.55555546E-10 ed_time_1_7 = 6.66666611E-10 ed_time_1_8 = 7.77777731E-10 ed_time_1_9 = 8.88888851E-10 ed_time_1_10 = 9.99999972E-10 ed_power_1_1 = 283173280. ed_power_1_2 = 5.93503744E+09 ed_power_1_3 = 5.81368668E+10 ed_power_1_4 = 2.66156540E+11 ed_power_1_5 = 5.69481953E+11 ed_power_1_6 = 5.69481953E+11 ed_power_1_7 = 2.66156671E+11 ed_power_1_8 = 5.81369078E+10 ed_power_1_9 = 5.93503744E+09 ed_power_1_10 = 283173280. #ed_power_1_5 = 1.113e12 #1.25e12 #2.3e14 #ed_power_1_6 = 1.113e12 #1.25e12 #2.3e14 #ed_adjustBeamsTargetIntensity = .false. ### SETUP LASER BEAM ### ed_numberOfBeams = 1 # Setup Gaussian Beam: ed_lensX_1 = 0.e-04 #25.e-04 ed_lensY_1 = 0.e-04 ed_lensZ_1 = 6400.e-04 ed_lensSemiAxisMajor_1 = 500.e-04 #25.0e-04 ed_targetX_1 = 0.e-04 #25.0e-04 ed_targetY_1 = 0.0e-04 ed_targetZ_1 = 12.5e-04 #400.0e-04 ed_targetSemiAxisMajor_1 = 500.0e-04 #25.0e-04 ed_targetSemiAxisMinor_1 = 500.0e-04 #25.0e-04 ed_pulseNumber_1 = 1 ed_wavelength_1 = 0.438 #0.381 #0.35 #in micron ed_crossSectionFunctionType_1 = "gaussian2D" #"gaussian2D" ed_gaussianExponent_1 = 1. ed_gaussianRadiusMajor_1 = 120.e-04 #25.e-04 ed_gaussianRadiusMinor_1 = 120.e-04 #25.e-04 ed_numberOfRays_1 = 100000#8400#4200 #4096 ed_gridType_1 = "radial2D" #"statistical2D" #"radial2D" ed_gridnRadialTics_1 = 600 #64#600 ed_semiAxisMajorTorsionAngle_1= 0.0 ed_semiAxisMajorTorsionAxis_1 = "x" ed gaussianCenterMajor_1 = 0. ed_gaussiancenterMinor_1 = 0. ################################# # # # CONDUCTION PARAMETERS # # # ################################# useDiffuse = .true. useConductivity = .true. diff_useEleCond = .true. diff_eleFlMode = "fl_minmax" #"fl_larsen" diff_eleFlCoef = 0.05 #0.03#0.05 #0.06 #0.1 diff_thetaImplct = 1.0 diff_eleXlBoundaryType = "neumann" diff_eleXrBoundaryType = "neumann" diff_eleYlBoundaryType = "neumann" diff_eleYrBoundaryType = "neumann" diff_eleZlBoundaryType = "neumann" diff_eleZrBoundaryType = "neumann" diff_useIonCond = .true. diff_ionFlMode = "fl_minmax" #"fl_larsen" diff_ionFlCoef = 0.05 #0.03#0.05 #0.06 #0.1 diff_thetaImplct = 1.0 diff_ionXlBoundaryType = "neumann" diff_ionXrBoundaryType = "neumann" diff_ionYlBoundaryType = "neumann" diff_ionYrBoundaryType = "neumann" diff_ionZlBoundaryType = "neumann" diff_ionZrBoundaryType = "neumann" #################################### # # # HEAT EXCHANGE PARAMETERS # # # #################################### useHeatexchange = .true. ########################## # # # EOS PARAMETERS # # # ########################## eosModeInit = "dens_temp_gather" smallt = 1.0 smallx = 1.0e-99 eos_useLogTables = .false. ############################ # # # HYDRO PARAMETERS # # # ############################ useHydro = .true. order = 3 # Interpolation order (first/second/third/fifth order) slopeLimiter = "minmod" # Slope limiters (minmod, mc, vanLeer, hybrid, limited) LimitedSlopeBeta = 1. # Slope parameter for the "limited" slope by Toro charLimiting = .true. # Characteristic limiting vs. Primitive limiting use_avisc = .true. # use artificial viscosity (originally for PPM) cvisc = 0.1 # coefficient for artificial viscosity use_flattening = .false. # use flattening (dissipative) (originally for PPM) use_steepening = .false. # use contact steepening (originally for PPM) use_upwindTVD = .false. # use upwind biased TVD slope for PPM (need nguard=6) RiemannSolver = "hllc" # Roe, HLL, HLLC, LLF, Marquina, hybrid entropy = .false. # Entropy fix for the Roe solver shockDetect = .false. # Shock Detect for numerical stability use_hybridOrder = .true. # Enforce Riemann density jump # Hydro boundary conditions: xl_boundary_type = "reflect" #"reflect" xr_boundary_type = "outflow" yl_boundary_type = "outflow" yr_boundary_type = "outflow" zl_boundary_type = "reflect" zr_boundary_type = "reflect" ############################## # # # INITIAL CONDITIONS # # # ############################## sim_targetRadius = 500.0e-04 #50.00e-04 100.00e-04 sim_targetHeight = 12.5e-04#400.0e-04 sim_vacuumHeight = 0.0e-04 # Target material defaults set for saran (C2H2Cl2) at room temperature: sim_rhoTarg = 1.701#1.e-2 #1.e-2#8.97 sim_teleTarg = 300.1#290.11375 sim_tionTarg = 300.1#290.11375 sim_tradTarg = 300.1#290.11375 ms_targA = 16.2 #7.2 #8.73#6.510 ms_targZ = 8. #3.85 #4.54#3.5 ms_targZMin = 0.02 #0.02#3.85#4.54#1.0 eos_targEosType = "eos_tab"#"eos_gam" #"eos_tab" #eos_singlespeciesA = 16.2 #7.2 #8.73 #eos_singlespeciesZ = 8 #3.85 #4.54 eos_targSubType = "ionmix4" eos_targTableFile = "saran-imx-30g.cn4" #"polystyrene-imx-008.cn4" # "he-imx-20g.cn4" #"cu-imx-20g.cn4" #"al-imx-003.cn4" # Chamber material defaults set for Helium at pressure 1.6 mbar: sim_rhoCham = 1.0e-06 sim_teleCham = 290.11375 sim_tionCham = 290.11375 sim_tradCham = 290.11375 ms_chamA = 4.002602 ms_chamZ = 2.0 ms_chamZMin = 0.02 #2.#0.02 eos_chamEosType = "eos_tab"#"eos_gam" #"eos_tab" eos_chamSubType = "ionmix4" eos_chamTableFile = "he-imx-20g.cn4" ########################### # # # TIME PARAMETERS # # # ########################### tstep_change_factor = 1.10 cfl = 0.2#0.4 dt_diff_factor = 1.0e+100 # Disable diffusion dt rt_dtFactor = 1.0e+100 hx_dtFactor = 1.0e+100 tmax = 1.e-9 #0.6e-9#0.5e-9 #1000.0e-12 dtmin = 1.0e-16 dtinit = 1.0e-15 dtmax = 1.e-13 #0.1e-12# 1.e-14 0.05e-12 nend = 10000000 ########################### # # # MESH PARAMETERS # # # ########################### geometry = "cylindrical" #"cylindrical" # Domain size: xmin = 0. xmax = 500.e-04 #50.0e-04 100.00e-04 ymin = 0.0e-04 ymax = 500.0e-04 # Total number of blocks: nblockx = 1 nblocky = 1 lrefine_max = 1 lrefine_min = 1 refine_var_1 = "dens" refine_var_2 = "tele"