-
Notifications
You must be signed in to change notification settings - Fork 98
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Create problem to study He flame movement in a double det #2870
Conversation
|
||
// Set up Castro data logs for this problem | ||
|
||
if (castro::sum_interval > 0 && amrex::ParallelDescriptor::IOProcessor()) { |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this entire if-block can be removed
state(i,j,k,UMY) = 0.0_rt; | ||
state(i,j,k,UMZ) = 0.0_rt; | ||
|
||
// Now add the velocity perturbation |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this entire if-block can be removed -- we will not do a velocity perturbation
Exec/science/subch_planar/inputs_2d
Outdated
# PROBLEM PARAMETERS | ||
problem.model_name = "toy_subch.hse.tanh.delta_50.00km.dx_73.24km" | ||
|
||
problem.apply_vel_field = 1 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
all of these problem.* parameters are not needed here
@@ -0,0 +1,19 @@ | |||
H_min real 1.e-4_rt y | |||
|
|||
cutoff_density real 50.0_rt y |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this parameter is not needed
|
||
model_name string "" y | ||
|
||
apply_vel_field integer 0 y |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this parameter and all following are not needed
Exec/science/subch_planar/inputs_2d
Outdated
amr.derive_plot_vars = ALL | ||
|
||
# PROBLEM PARAMETERS | ||
problem.model_name = "toy_subch.hse.tanh.delta_50.00km.dx_73.24km" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this is the wrong initial model file
# This sets the EOS directory in $(MICROPHYSICS_HOME)/eos | ||
EOS_DIR := helmholtz | ||
|
||
# This sets the EOS directory in $(MICROPHYSICS_HOME)/networks |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
EOS -> network
Exec/science/subch_planar/Problem.H
Outdated
@@ -0,0 +1 @@ | |||
amrex::Real maxVal (const std::string& name, amrex::Real time); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this file is not needed
PR summary
Begin creating a problem to study the lateral movement of a He flame in the accreted layer of a white dwarf.
PR motivation
We want to determine how various shock detection thresholds and burning treatments affect the initial flame movement in a double detonation model. In particular, we seek to see if the results of disabling/enabling burning in detected shock zones converges as resolution becomes finer.
PR checklist
CHANGES
file has been updated, if appropriate