run_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 = 1.0 checkpointFileIntervalStep = 1000 ### Plot Options ### plotFileNumber = 0 plotFileIntervalStep = 6000 plotFileIntervalTime = 0.2e-09 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 = "volx" plot_var_11 = "voly" plot_var_12 = "volz" plot_var_13 = "velx" plot_var_14 = "vely" plot_var_15 = "velz" ### Restart Options ### restart = .false. checkpointFileNumber = 0 ######################################## # # # RADIATION/OPACITY PARAMETERS # # # ######################################## rt_useMGD = .true. rt_mgdNumGroups = 6 rt_mgdBounds_1 = 1.0e-01 rt_mgdBounds_2 = 1.0e+00 rt_mgdBounds_3 = 1.0e+01 rt_mgdBounds_4 = 1.0e+02 rt_mgdBounds_5 = 1.0e+03 rt_mgdBounds_6 = 1.0e+04 rt_mgdBounds_7 = 1.0e+05 rt_mgdFlMode = "fl_harmonic" rt_mgdFlCoef = 1.0 rt_mgdXlBoundaryType = "reflecting" rt_mgdXrBoundaryType = "vacuum" rt_mgdYlBoundaryType = "vacuum" rt_mgdYrBoundaryType = "reflecting" rt_mgdZlBoundaryType = "reflecting" rt_mgdZrBoundaryType = "reflecting" useOpacity = .true. ### SET CHAMBER (HELIUM) OPACITY OPTIONS ### op_chamAbsorb = "op_tabpa" op_chamEmiss = "op_tabpe" op_chamTrans = "op_tabro" op_chamFileType = "ionmix4" op_chamFileName = "he-imx-005.cn4" ### SET TARGET (ALUMINUM) OPACITY OPTIONS ### op_targAbsorb = "op_tabpa" op_targEmiss = "op_tabpe" op_targTrans = "op_tabro" op_targFileType = "ionmix4" op_targFileName = "Sn.cn4" ############################ # # # LASER PARAMETERS # # # ############################ useEnergyDeposition = .true. ed_maxRayCount = 100000 ed_gradOrder = 2 # Activate 3D-in-2D ray trace: ed_laser3Din2D = .false. ed_laser3Din2DwedgeAngle = 0.1 # Activate Beams Target Intensity adjust Yaoxing Wu ed_adjustBeamsTargetIntensity = .ture. ### LASER IO OPTIONS ### ed_useLaserIO = .true. ed_laserIOMaxNumberOfPositions = 10000 ed_laserIOMaxNumberOfRays = 128 ### SETUP LASER PULSES ### ed_numberOfPulses = 1 # Define Pulse 1: ed_numberOfSections_1 = 33 ed_time_1_1 = 0 ed_time_1_2 = 5.000000e-10 ed_time_1_3 = 1.000000e-09 ed_time_1_4 = 1.500000e-09 ed_time_1_5 = 2.000000e-09 ed_time_1_6 = 2.500000e-09 ed_time_1_7 = 3.000000e-09 ed_time_1_8 = 3.500000e-09 ed_time_1_9 = 4.000000e-09 ed_time_1_10 = 4.500000e-09 ed_time_1_11 = 5.000000e-09 ed_time_1_12 = 5.500000e-09 ed_time_1_13 = 6.000000e-09 ed_time_1_14 = 6.500000e-09 ed_time_1_15 = 7.000000e-09 ed_time_1_16 = 7.500000e-09 ed_time_1_17 = 8.000000e-09 ed_time_1_18 = 8.500000e-09 ed_time_1_19 = 9.000000e-09 ed_time_1_20 = 9.500000e-09 ed_time_1_21 = 1.000000e-08 ed_time_1_22 = 1.050000e-08 ed_time_1_23 = 1.100000e-08 ed_time_1_24 = 1.150000e-08 ed_time_1_25 = 1.200000e-08 ed_time_1_26 = 1.250000e-08 ed_time_1_27 = 1.300000e-08 ed_time_1_28 = 1.350000e-08 ed_time_1_29 = 1.400000e-08 ed_time_1_30 = 1.450000e-08 ed_time_1_31 = 1.500000e-08 ed_time_1_32 = 1.550000e-08 ed_time_1_33 = 1.600000e-08 ed_power_1_1 =1.912586e+05 ed_power_1_2 =3.555360e+05 ed_power_1_3 =6.350010e+05 ed_power_1_4 =1.089666e+06 ed_power_1_5 =1.796556e+06 ed_power_1_6 =2.845877e+06 ed_power_1_7 =4.331316e+06 ed_power_1_8 =6.333616e+06 ed_power_1_9 =8.898398e+06 ed_power_1_10 =1.201158e+07 ed_power_1_11 =1.557818e+07 ed_power_1_12 =1.941161e+07 ed_power_1_13 =2.323992e+07 ed_power_1_14 =2.673227e+07 ed_power_1_15 =2.954372e+07 ed_power_1_16 =3.137060e+07 ed_power_1_17 =3.200433e+07 ed_power_1_18 =3.137060e+07 ed_power_1_19 =2.954372e+07 ed_power_1_20 =2.673227e+07 ed_power_1_21 =2.323992e+07 ed_power_1_22 =1.941161e+07 ed_power_1_23 =1.557818e+07 ed_power_1_24 =1.201158e+07 ed_power_1_25 =8.898398e+06 ed_power_1_26 =6.333616e+06 ed_power_1_27 =4.331316e+06 ed_power_1_28 =2.845877e+06 ed_power_1_29 =1.796556e+06 ed_power_1_30 =1.089666e+06 ed_power_1_31 =6.350010e+05 ed_power_1_32 =3.555360e+05 ed_power_1_33 =1.912586e+05 ### SETUP LASER BEAM ### ed_numberOfBeams = 1 # Setup Gaussian Beam1: ed_lensX_1 = 0.0e-04 ed_lensY_1 = 0.0e-04 ed_lensZ_1 = -6000.0e-04 ed_lensSemiAxisMajor_1 = 62.5e-04 ed_targetX_1 = 0.0e-04 ed_targetY_1 = 0.0e-04 ed_targetZ_1 = 0.0e-04 ed_targetSemiAxisMajor_1 = 62.5e-04 ed_targetSemiAxisMinor_1 = 62.5e-04 ed_pulseNumber_1 = 1 ed_wavelength_1 = 1.064 ed_crossSectionFunctionType_1 = "gaussian1D" ed_gaussianExponent_1 = 2.0 ed_gaussianRadiusMajor_1 = 62.5e-04 ed_gaussianRadiusMinor_1 = 62.5e-04 ed_numberOfRays_1 = 4096 ed_gridType_1 = "regular1D" ed_gridnRadialTics_1 = 64 ed_semiAxisMajorTorsionAngle_1= 0.0 ed_semiAxisMajorTorsionAxis_1 = "x" # Setup Gaussian Beam2: # ed_lensX_2 = 0.0e-04 # ed_lensY_2 = 0.0e-04 # ed_lensZ_2 = -6000.0e-04 # ed_lensSemiAxisMajor_2 = 125e-04 # ed_targetX_2 = 0.0e-04 # ed_targetY_2 = 0.0e-04 # ed_targetZ_2 = 0.0e-04 # ed_targetSemiAxisMajor_2 = 125e-04 # ed_targetSemiAxisMinor_2 = 125e-04 # ed_pulseNumber_2 = 1 # ed_wavelength_2 = 1.064 # ed_crossSectionFunctionType_2 = "gaussian2D" # ed_gaussianExponent_2 = 2.0 # ed_gaussianRadiusMajor_2 = 125e-04 # ed_gaussianRadiusMinor_2 = 125e-04 # ed_numberOfRays_2 = 4096 # ed_gridType_2 = "radial2D" # ed_gridnRadialTics_2 = 64 # ed_semiAxisMajorTorsionAngle_2= 0.0 # ed_semiAxisMajorTorsionAxis_2 = "x" ################################# # # # CONDUCTION PARAMETERS # # # ################################# useDiffuse = .true. useConductivity = .true. diff_useEleCond = .true. diff_eleFlMode = "fl_larsen" diff_eleFlCoef = 0.08 diff_thetaImplct = 1.0 diff_eleXlBoundaryType = "neumann" diff_eleXrBoundaryType = "neumann" diff_eleYlBoundaryType = "neumann" diff_eleYrBoundaryType = "neumann" diff_eleZlBoundaryType = "neumann" diff_eleZrBoundaryType = "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" #xl_boundary_type = "outflow" xr_boundary_type = "outflow" yl_boundary_type = "outflow" yr_boundary_type = "outflow" zl_boundary_type = "reflect" zr_boundary_type = "reflect" ############################## # # # INITIAL CONDITIONS # # # ############################## sim_targetRadius = 100.0e-04 #target radius sim_targetHeight = 100.0e-04 #target thickness sim_vacuumHeight = 0.0e-04 #target position # Target material defaults set for Aluminum at room temperature: sim_rhoTarg = 7.28 #target density sim_teleTarg = 290.11375 #intial tem(K) sim_tionTarg = 290.11375 sim_tradTarg = 290.11375 ms_targA = 118.71 #target atom weight ms_targZ = 50.0 #target Z ms_targZMin = 0.02 eos_targEosType = "eos_tab" eos_targSubType = "ionmix4" eos_targTableFile = "Sn.cn4" # Chamber material defaults set for Helium at pressure 1.6 mbar:1.6 mbar = 1.0e-6, 3e-6 torr = 1.7e-10 sim_rhoCham = 1.0e-8 sim_teleCham = 290.11375 sim_tionCham = 290.11375 sim_tradCham = 290.11375 ms_chamA = 4.002602 ms_chamZ = 2.0 eos_chamEosType = "eos_tab" eos_chamSubType = "ionmix4" eos_chamTableFile = "he-imx-005.cn4" ########################### # # # TIME PARAMETERS # # # ########################### tstep_change_factor = 1.10 cfl = 0.4 dt_diff_factor = 1.0e+100 # Disable diffusion dt rt_dtFactor = 1.0e+100 hx_dtFactor = 1.0e+100 tmax = 16.0e-09 dtmin = 1.0e-16 dtinit = 1.0e-15 dtmax = 3.0e-09 nend = 10000000 ########################### # # # MESH PARAMETERS # # # ########################### geometry = "cartesian" # Domain size: xmin = 0.0 xmax = 2000.0e-04 #unit 2000 um ymin = -3000.0e-04 #unit ymax = 1000.0e-04 # Total number of blocks: nblockx = 1 #intial mesh nblocky = 2 lrefine_max = 4 # max split mesh lrefine_min = 1 refine_var_1 = "dens" refine_var_2 = "tele"