Skip to content
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

Merged
merged 16 commits into from
Jun 18, 2024
Merged

Create problem to study He flame movement in a double det #2870

merged 16 commits into from
Jun 18, 2024

Conversation

brady-ryan
Copy link
Contributor

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

  • test suite needs to be run on this PR
  • this PR will change answers in the test suite to more than roundoff level
  • all newly-added functions have docstrings as per the coding conventions
  • the CHANGES file has been updated, if appropriate
  • if appropriate, this change is described in the docs


// Set up Castro data logs for this problem

if (castro::sum_interval > 0 && amrex::ParallelDescriptor::IOProcessor()) {
Copy link
Member

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
Copy link
Member

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

# PROBLEM PARAMETERS
problem.model_name = "toy_subch.hse.tanh.delta_50.00km.dx_73.24km"

problem.apply_vel_field = 1
Copy link
Member

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
Copy link
Member

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
Copy link
Member

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

amr.derive_plot_vars = ALL

# PROBLEM PARAMETERS
problem.model_name = "toy_subch.hse.tanh.delta_50.00km.dx_73.24km"
Copy link
Member

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
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

EOS -> network

@@ -0,0 +1 @@
amrex::Real maxVal (const std::string& name, amrex::Real time);
Copy link
Member

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

@zingale zingale merged commit 709aec4 into AMReX-Astro:development Jun 18, 2024
24 of 25 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants