[FLASH-USERS] Strong differences between 4.7.1 and 4.8
Dr. Arthur Hirsch-Passicos
a.hirsch at hzdr.de
Tue Sep 3 03:24:40 EDT 2024
Dear Flash users,
I noticed some strong differences between my FLASH cases between 4.7.1
and 4.8. As I had changed quite a few parameters from the base of
LaserSlab to my case, I came back to LaserSlab and applies variations of
my parameters progressively.
While the base LaserSlab case presented some small variations (of the
order of a few % in ion temperature in the joined images), the
differences became massive when I switched from cylindrical to cartesian
grid, with a factor 4 increase of the ion temperature in the joined
images, increasing to a factor 10 at the end of the simulation (Both set
of images joined).
I have seen the same kind of results on short pulse simulations and on
lower timescale, making me think it is not due to the change on
HeatExchange and PlasmaState decribed in the changelog. As I see it,
they would improve the propagation of absorbed energy through the target
and not increase the level of absorbed energy. But I have not find a
change in the energy deposition that would explain any of the observed
results.
Does any of you have an explanation of what would be the reason for such
changes?
I joined images of the ion temperature in both versions for the
cartesian and cylindrical simulations as well as the par file for the
cartesian case.
Thank you in advance for your help
Best regards,
--
Dr. Arthur Hirsch-Passicos
Post-doctoral Researcher
Laser Particle Acceleration Division
Institute of Radiation Physics
Helmholtz-Zentrum Dresden-Rossendorf e.V.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cartesian_0P18ns_Flash48.png
Type: image/png
Size: 29414 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0008.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cartesian_0P18ns_Flash471.png
Type: image/png
Size: 30077 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0009.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cartesian_2ns_Flash48.png
Type: image/png
Size: 38383 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0010.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cartesian_2ns_Flash471.png
Type: image/png
Size: 38689 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0011.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cylindrical_0P18ns_Flash48.png
Type: image/png
Size: 55969 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0012.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cylindrical_0P18ns_Flash471.png
Type: image/png
Size: 53614 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0013.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cylindrical_2ns_Flash48.png
Type: image/png
Size: 54293 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0014.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cylindrical_2ns_Flash471.png
Type: image/png
Size: 54587 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240903/cc704c25/attachment-0015.png>
-------------- next part --------------
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 = 100
plotFileIntervalTime = 0.01e-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"
### 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 = "al-imx-003.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 = 0.1e-09
ed_time_1_3 = 1.0e-09
ed_time_1_4 = 1.1e-09
ed_power_1_1 = 0.0
ed_power_1_2 = 1.0e+09
ed_power_1_3 = 1.0e+09
ed_power_1_4 = 0.0
### SETUP LASER BEAM ###
ed_numberOfBeams = 1
# Setup Gaussian Beam:
ed_lensX_1 = 1000.0e-04
ed_lensY_1 = 0.0e-04
ed_lensZ_1 = -1000.0e-04
ed_lensSemiAxisMajor_1 = 10.0e-04
ed_targetX_1 = 0.0e-04
ed_targetY_1 = 0.0e-04
ed_targetZ_1 = 60.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 = 4.0
ed_gaussianRadiusMajor_1 = 7.5e-04
ed_gaussianRadiusMinor_1 = 7.5e-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 = 20.0e-04
sim_vacuumHeight = 60.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.9815386
ms_targZ = 13.0
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.4
dt_diff_factor = 1.0e+100 # Disable diffusion dt
rt_dtFactor = 1.0e+100
hx_dtFactor = 1.0e+100
tmax = 2.0e-09
dtmin = 1.0e-16
dtinit = 1.0e-15
dtmax = 3.0e-09
nend = 10000000
###########################
# #
# MESH PARAMETERS #
# #
###########################
geometry = "cylindrical"
# Domain size:
xmin = 0.0
xmax = 40.0e-04
ymin = 0.0e-04
ymax = 80.0e-04
# Total number of blocks:
nblockx = 1
nblocky = 2
lrefine_max = 4
lrefine_min = 1
refine_var_1 = "dens"
refine_var_2 = "tele"
More information about the flash-users
mailing list