[FLASH-USERS] multispecies ionmix not found
Aodhan McIlvenny
amcilvenny01 at qub.ac.uk
Wed Feb 12 01:30:27 EST 2020
Hi all,
I have 2 questions which I hope are easy to answer; I am a bit stuck
when trying to change the material in the laserslab example.
* I am trying to use the polystyrene ionmix4 file which is in the same
format as the aluminium one however it returns an error that it
can't be found. Any ideas why this is so?
* Secondly, can I track the carbon and proton species separately from
this one material? Ideally I would like to return the position and
temperature of carbon and protons separately. Is there a partciular
syntax for looking at each ion species from 1 material?
I have attached my .par file. I would much appreciate some help! I am
using FLASH.4.6.2 and setup command was
./setup -auto LaserSlab -2d +cylindrical -nxb=16 -nyb=16 +hdf5typeio
species=cham,targ +mtmmmt +laser +uhd3t +mgd mgd_meshgroups=6
-parfile=VulcanTest2.par
Thanks in advance,
Aodhan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20200212/aab5ee3f/attachment.htm>
-------------- next part --------------
run_comment = "VulcanTestrun"
log_file = "VulcanTest.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 = 100
plotFileIntervalTime = 50e-15
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"
### 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 = "vacuum"
rt_mgdZrBoundaryType = "vacuum"
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 = "polystyrene-imx-002.cn4"
############################
# #
# LASER PARAMETERS #
# #
############################
useEnergyDeposition = .true.
ed_maxRayCount = 10000
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 = 4
ed_time_1_1 = 0.0
ed_time_1_2 = 1e-12
ed_time_1_3 = 15e-12
ed_time_1_4 = 20e-12
ed_power_1_1 = 0.0
ed_power_1_2 = 5.35e+12
ed_power_1_3 = 5.35e+13
ed_power_1_4 = 5.35e+14
### SETUP LASER BEAM ###
ed_numberOfBeams = 1
# Setup Gaussian Beam:
ed_lensX_1 = 15.0e-04
ed_lensY_1 = 0.0e-04
ed_lensZ_1 = -1000.0e-04
ed_lensSemiAxisMajor_1 = 10.0e-04
ed_targetX_1 = 15.0e-04
ed_targetY_1 = 0.0e-04
ed_targetZ_1 = 17.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 = "gaussian2D"
ed_gaussianExponent_1 = 2.0
ed_gaussianRadiusMajor_1 = 4.27e-04
ed_gaussianRadiusMinor_1 = 4.27e-04
ed_numberOfRays_1 = 4096
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.06
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"
xr_boundary_type = "outflow"
yl_boundary_type = "outflow"
yr_boundary_type = "outflow"
zl_boundary_type = "reflect"
zr_boundary_type = "reflect"
##############################
# #
# INITIAL CONDITIONS #
# #
##############################
sim_targetRadius = 200.0e-04
sim_targetHeight = 1.0e-04
sim_vacuumHeight = 17.0e-04
# Target material defaults set for Aluminum at room temperature:
sim_rhoTarg = 1.2
sim_teleTarg = 290.11375
sim_tionTarg = 290.11375
sim_tradTarg = 290.11375
ms_targA = 12
ms_targZ = 6
ms_targZMin = 0.02
eos_targEosType = "eos_tab"
eos_targSubType = "ionmix4"
eos_targTableFile = "polystyrene-imx-002.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.4
dt_diff_factor = 1.0e+100 # Disable diffusion dt
rt_dtFactor = 1.0e+100
hx_dtFactor = 1.0e+100
dtmin = 1.0e-16
dtinit = 1.0e-15
dtmax = 3.0e-09
tmax = 20.0e-12
nend = 10000000
###########################
# #
# MESH PARAMETERS #
# #
###########################
geometry = "cylindrical"
# Domain size:
xmin = 0.0e-4
xmax = 30.0e-04
ymin = 0.0e-04
ymax = 20.0e-04
# Total number of blocks:
nblockx = 1
nblocky = 2
lrefine_max = 4
lrefine_min = 2
refine_var_1 = "dens"
refine_var_2 = "tele"
More information about the flash-users
mailing list