[FLASH-USERS] FLASH: low laser power spherical target interaction
Majid Masnavi
majidmasnavi at creol.ucf.edu
Tue Jul 21 17:11:09 EDT 2015
Hello FLASH-users,
The FLASH code gives an unexpected low electron temperature (~ 10 eV) plasma from a laser-spherical (solid or liquid) Al (and perhaps other elements) target interaction. The maximum laser power is ~ (3-4)x10^6 W. The laser rise time is ~ 17-20 ns. The experimental results demonstrated plasma with temperature > 20 eV. Is there any mistake in the following file? I would be very happy to recieve your comments.
Best
Majid
# Run the spherical target
# setup call:
# ./setup ... -auto -2d +hdf5typeio species=cham,targ +mtmmmt +laser +uhd3t +mgd mgd_meshgroups=6 -parfile=flash.par -objdir=...
#
run_comment = "test simulation of ..."
log_file = "....log"
basenm = "..._"
# ...
#########################
# OUTPUT PARAMETERS #
#########################
### Checkpoint Options ###
checkpointFileIntervalTime = 1.0
checkpointFileIntervalStep = 4
### Plot Options ###
plotFileNumber = 0
plotFileIntervalStep = 10
plot_var_1 = "dens"
plot_var_2 = "depo"
plot_var_3 = "tele"
plot_var_4 = "tion"
plot_var_5 = "ye "
plot_var_6 = "sumy"
### 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 = "vacuum"
rt_mgdXrBoundaryType = "vacuum"
rt_mgdYlBoundaryType = "vacuum"
rt_mgdYrBoundaryType = "vacuum"
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 = "al-imx-003.cn4"
############################
# #
# LASER PARAMETERS #
# #
############################
ed_enforcePositiveTele = .true.
ed_enforcePositiveNele = .true.
useEnergyDeposition = .true. #turn on the laser
ed_maxRayCount = 2000
ed_gradOrder = 2
### LASER IO OPTIONS ###
ed_useLaserIO = .true.
ed_laserIOMaxNumberOfPositions = 10000
ed_laserIOMaxNumberOfRays = 128
### SETUP LASER PULSES ###
ed_numberOfPulses = 1
# Define Pulse 1:
# in nano second scale (e-09)
ed_numberOfSections_1 =3
ed_time_1_1 = 0.0
ed_time_1_2 = 20.0e-09
ed_time_1_3 = 2.10e-09
# in scale (e+ watt)
ed_power_1_1 = 4.0e+3 (4.0e+2, ...)
ed_power_1_2 = 4.0e+6
ed_power_1_3 = 0
### SETUP LASER BEAM ###
ed_numberOfBeams = 1
ed_printBeams = .true.
ed_printPulses = .true.
ed_printRays = .true.
ed_printMain = .true.
# Setup Gaussian Beam 1:
ed_lensX_1 = 500.0e-04 #x-position of L
ed_lensY_1 = 0 #y-position of L
ed_targetX_1 = 20.0e-04 #x-position of T
ed_targetY_1 = 0 #y-position of T
ed_lensSemiAxisMajor_1 = 10.0e-04 (...) #
ed_targetSemiAxisMajor_1 = 10.0e-04 (...)
ed_targetSemiAxisMinor_1 = 10.0e-04 (...)
ed_pulseNumber_1 = 1
ed_wavelength_1 = 1.053
ed_crossSectionFunctionType_1 = "gaussian1D"
ed_gaussianExponent_1 = 1.0 (...) #
ed_gaussianRadiusMajor_1 = 5.50e-04 (...)
ed_gaussianRadiusMinor_1 = 5.50e-04 (...)
ed_numberOfRays_1 = 512
ed_gridType_1 = "regular1D"
#ed_semiAxisMajorTorsionAngle_1 = 0.0
#ed_semiAxisMajorTorsionAxis_1 = "x"
#################################
# #
# CONDUCTION PARAMETERS #
# #
#################################
useDiffuse = .true.
useConductivity = .true.
diff_useEleCond = .true.
diff_eleFlMode = "fl_minmax"
diff_eleFlCoef = 0.06
diff_thetaImplct = 1.0
diff_useIonCond = .true.
diff_ionFlMode = "fl_harmonic"
diff_ionFlCoef = 0.8
diff_ionthetaImplct = 1.0
diff_eleXlBoundaryType = "outflow"
diff_eleXrBoundaryType = "outflow"
diff_eleYlBoundaryType = "outflow"
diff_eleYrBoundaryType = "outflow"
diff_eleZlBoundaryType = "outflow"
diff_eleZrBoundaryType = "outflow"
diff_ionXlBoundaryType = "outflow"
diff_ionXrBoundaryType = "outflow"
diff_ionYlBoundaryType = "outflow"
diff_ionYrBoundaryType = "outflow"
diff_ionZlBoundaryType = "outflow"
diff_ionZrBoundaryType = "outflow"
####################################
# #
# HEAT EXCHANGE PARAMETERS #
# #
####################################
useHeatexchange = .true.
##########################
# #
# EOS PARAMETERS #
# #
##########################
eosModeInit = "dens_temp_gather"
eosMode = "dens_ie_recal_gather"
smallt = 1.0
smallx = 1.0e-99
eos_useLogTables = .false.
############################
# #
# HYDRO PARAMETERS #
# #
############################
# Use second order hybrid solver with minmod slope limiter. This
# essentially eliminates any Carbuncle instability.
order = 3 # Interpolation order (first/second/third/fifth order)
slopeLimiter = "mc" # 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 = .true. # Entropy fix for the Roe solver
shockDetect = .true. # Shock Detect for numerical stability
addThermalFlux = .false.
use_hybridOrder = .true.
# Hydro boundary conditions:
xl_boundary_type = "outflow"
xr_boundary_type = "outflow"
yl_boundary_type = "outflow"
yr_boundary_type = "outflow"
zl_boundary_type = "outflow"
zr_boundary_type = "outflow"
##########################
# INITIAL CONDITIONS #
##########################
sim_initGeom = "sphere"
sim_targetRadius = 0.0
sim_shellThickness = 20.0e-04
# Target material defaults set for Aluminum at room temperature:
sim_rhoTarg = 2.7
sim_teleTarg = 290.11375
sim_tionTarg = 290.11375
sim_tradTarg = 290.11375
ms_targA = 26.9
ms_targZ = 13
ms_targZMin = 0.02
eos_targEosType = "eos_tab"
eos_targSubType = "ionmix4"
eos_targTableFile = "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
eos_chamEosType = "eos_tab"
eos_chamSubType = "ionmix4"
eos_chamTableFile = "he-imx-005.cn4"
#######################
# TIME PARAMETERS #
#######################
tstep_change_factor = 1.10
cfl = 0.2
dt_diff_factor = 1.0e+100 # Disable diffusion dt
rt_dtFactor = 0.02
hx_dtFactor = 1.0e+100
tmax = 21.0e-09
dtmin = 1.0e-16
dtinit = 1.0e-15
dtmax = 3.0e-09
nend = 100000000
#######################
# MESH PARAMETERS #
#######################
geometry = "cartesian"
# Domain size:
xmin = -100.0e-04
xmax = 100.0e-04
ymin = -100.0e-04
ymax = 100.0e-04
# Total number of blocks:
#nblockx = 2
#nblocky = 5
nblockx = 1
nblocky = 1
lrefine_max = 4
lrefine_min = 2
refine_var_1 = "dens"
refine_var_2 = "tele"
refine_var_3 = "tion"
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20150721/2818e8ca/attachment.htm>
More information about the flash-users
mailing list