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 ### Plot Options ### plotFileNumber = 0 plotFileIntervalTime = 0.1e-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 = "pres" ### 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 (Copper) OPACITY OPTIONS ### op_targAbsorb = "op_tabpa" op_targEmiss = "op_tabpe" op_targTrans = "op_tabro" op_targFileType = "ionmix4" op_targFileName = "cu_op.cn4" ############################ # # # LASER PARAMETERS # # # ############################ useEnergyDeposition = .true. # Turn on laser energy deposition: Activates laser model ed_maxRayCount = 10000 # Max. number of rays which can be stored on a given process at onxe ed_gradOrder = 2 # Trun on laser ray refraction: Controls how electron number density and temperature are interpolated within a cell; 2 = lin. interpolation # Activate 3D-in-2D ray trace: ed_laser3Din2D = .true. # Activates the 3D-in-2D ray trace algorithm; No effect on single beam centered on z-plane; only for cylindrical simulation ed_laser3Din2DwedgeAngle = 0.1 # sets the wedge angle (should be between 0.1 and 1.0 degrees) ### LASER IO OPTIONS ### ed_useLaserIO = .true. # Options tell FLASH to write out the trajectory of ed_laserIOMaxNumberOfRays to the plot files ed_laserIOMaxNumberOfPositions = 10000 ed_laserIOMaxNumberOfRays = 128 ### SETUP LASER PULSES ### ed_numberOfPulses = 1 # Define Pulse 1: ed_numberOfSections_1 = 3 ed_time_1_1 = 0.0 # laser time points in seconds ed_time_1_2 = 1.0e-09 ed_time_1_3 = 2.5e-09 ed_power_1_1 = 0.0 # laser power points in watt ed_power_1_2 = 9.064720e+06 ed_power_1_3 = 9.064720e+06 # ed_adjustBeamsTargetIntensity = .true. # Adjusting the laser power automatically from 3D to 2D/1D ### SETUP LASER BEAM ### ed_numberOfBeams = 1 # Setup Gaussian Beam: ed_lensX_1 = 1000e-04 # Define the x-coordinate of the center of the lense ed_lensY_1 = 0e-04 # Define the y-coordinate of the center of the lense ed_lensZ_1 = 1007e-04 #why not defined? # Define the z-coordinate of the center of the lense ed_lensSemiAxisMajor_1 = 20.0e-04 # Set radius on the lense ed_targetX_1 = 0.0e-04 # Define the x-coordinate of the center of the laser focal spot; x = R direction ed_targetY_1 = 0.0e-04 # Define the y-coordinate of the center of the laser focal spot; y direction in the R-z simulation ed_targetZ_1 = 7.0e-04 #why not defined? # Define the z-coordinate of the center of the laser focal spot; z direction pointing out of the computer screen ed_targetSemiAxisMajor_1 = 10.0e-04 # Set radius of the laser spot ed_targetSemiAxisMinor_1 = 10.0e-04 #only using for an elliptical laser shape ed_pulseNumber_1 = 1 ed_wavelength_1 = 0.527 # wavelength unit: mum ed_crossSectionFunctionType_1 = "gaussian2D" # Intensity-Function ed_gaussianExponent_1 = 2.0 # exponent gamma ed_gaussianRadiusMajor_1 = 2.4022448175729e-04 # R_x ed_gaussianRadiusMinor_1 = 2.4022448175729e-04 # only using this for an elliptical laser shap ed_numberOfRays_1 = 4096 # Flash launch this amount of rays for the first beam every time the EnergyDeposition subroutine is called (here: called once) ed_gridType_1 = "radial2D" ed_gridnRadialTics_1 = 64 ed_semiAxisMajorTorsionAngle_1= 0.0 ed_semiAxisMajorTorsionAxis_1 = "x" ################################# # # # CONDUCTION PARAMETERS # # # ################################# useDiffuse = .true. useConductivity = .true. diff_useEleCond = .true. diff_eleFlMode = "fl_larsen" diff_eleFlCoef = 0.05 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 = 6800 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 shockLowerCFL = .true. # Shock Detect for numerical stability use_hybridOrder = .true. # Enforce Riemann density jump # Hydro boundary conditions: xl_boundary_type = "reflect" xr_boundary_type = "outflow" yl_boundary_type = "outflow" yr_boundary_type = "outflow" zl_boundary_type = "reflect" zr_boundary_type = "reflect" ############################## # # # INITIAL CONDITIONS # # # ############################## # define space occupied by the target material sim_targetRadius = 200.0e-04 # 0 < R < sim_targetRadius sim_targetHeight = 5.0e-04 sim_vacuumHeight = 7.0e-04 # sim_vacuumHeight < z < sim_vacuumHeight + sim_targetHeight # Target material: sim_rhoTarg = 8.96 # Set target material density in g/cc sim_teleTarg = 6800 # Set initial electron temperature sim_tionTarg = 6800 # Set initial ion temperature sim_tradTarg = 6800 # Set initial radiation temperature ms_targA = 63.546 # Set average atomic mass in units of dalton (unified atomic mass unit u) ms_targZ = 29.0 # Set average atomic number/Ordnungszahl ms_targZMin = 0.02 # Set minimum allowed ionization level for targ eos_targEosType = "eos_tab" # Set EOS model to tabulated EOS eos_targSubType = "ionmix4" # Set EOS table file to IONMIX 4 format eos_targTableFile = "copper_3336.cn4" # Set tabulated EOS file name # Chamber material defaults set for Helium at pressure 1.6 mbar: sim_rhoCham = 1.0e-06 # Initial Helium density is chosen to be fairly low (vacuum) sim_teleCham = 6800 # Set initial electron temperature sim_tionCham = 6800 # Set initial ion temperature sim_tradCham = 6800 # Set initial radiation temperature ms_chamA = 4.002602 # Set average atomic mass in units of dalton (unified atomic mass unit u) ms_chamZ = 2.0 # Set average atomic number/Ordnungszahl eos_chamEosType = "eos_tab" # Set EOS model to tabulated EOS eos_chamSubType = "ionmix4" # Set EOS table file to IONMIX 4 format eos_chamTableFile = "he-imx-005.cn4" # Set tabulated EOS file name ########################### # # # 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 = 1.0e-09 dtmin = 1.0e-16 dtinit = 1.0e-15 dtmax = 3.0e-09 nend = 10000000 ########################### # # # MESH PARAMETERS # # # ########################### geometry = "cylindrical" # Domain size: xmin = 0.0e-04 xmax = 65.0e-04 ymin = -5.0e-04 ymax = 70.0e-04 # Total number of blocks: nblockx = 2 nblocky = 2 lrefine_max = 4 lrefine_min = 1 refine_var_1 = "dens" refine_var_2 = "tele"