Skip to content

allow for any C/O composition in the core #76

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.

Sign up for GitHub

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

Open
wants to merge 1 commit into
base: main
Choose a base branch
from
Open

allow for any C/O composition in the core #76

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions sub_chandra/_parameters
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,9 @@ X_C12 real 0.00
# mass fraction of O16 to put in the accreted envelope
X_O16 real 0.00

# mass fraction of core O16 composition
core_O16 real 0.5

delta real 1.e-6

xmin real 0.0
Expand Down
11 changes: 4 additions & 7 deletions sub_chandra/init_1d.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -15,7 +15,6 @@

using namespace amrex;


// generate an initial model for an arbitrary-mass, isothermal C WD
// with an isentropic He envelope on the surface.

Expand All @@ -37,7 +36,6 @@ AMREX_INLINE void init_1d() {

const int MAX_ITER = 1000;


// convert the envelope and WD mass into CGS

problem_rp::M_tot *= C::M_solar;
Expand All @@ -62,13 +60,12 @@ AMREX_INLINE void init_1d() {
Real xn_core[NumSpec] = {network_rp::small_x};
Real xn_he[NumSpec] = {network_rp::small_x};;

if (problem_rp::mixed_co_wd) {
xn_core[ic12] = 0.5_rt - 0.5_rt * (NumSpec - 1) * network_rp::small_x;
xn_core[io16] = 0.5_rt - 0.5_rt * (NumSpec - 1) * network_rp::small_x;
} else {
xn_core[ic12] = 1.0_rt - (NumSpec - 1) * network_rp::small_x;
if (problem_rp::core_O16 > 0.0_rt) {
xn_core[io16] = problem_rp::core_O16;
}

xn_core[ic12] = 1.0_rt - xn_core[io16] - (NumSpec - 2) * network_rp::small_x;

if (problem_rp::X_N14 > 0.0_rt) {
xn_he[in14] = problem_rp::X_N14;
}
Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-0.9.M_He-0.1.CO.CNO
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 5.e6
problem.temp_core = 3.e7
problem.temp_base = 2.3e8

problem.mixed_co_wd = 1

problem.X_C12 = 0.05
problem.X_O16 = 0.05
problem.X_N14 = 0.009
Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-0.981.M_He-0.05.CO.CO
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 5.e6
problem.temp_core = 5.e7
problem.temp_base = 2.e8

problem.mixed_co_wd = 1

problem.X_C12 = 0.05
problem.X_O16 = 0.05

Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-0.981.M_He-0.05.CO.CO.isothermal
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 5.e6
problem.temp_core = 5.e7
problem.temp_base = 2.e8

problem.mixed_co_wd = 1

problem.X_C12 = 0.05
problem.X_O16 = 0.05

Expand Down
2 changes: 1 addition & 1 deletion sub_chandra/inputs.M_WD-1.1.M_He-0.05
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ problem.delta = 5.e6
problem.temp_core = 1.e7
problem.temp_base = 1.75e8

problem.mixed_co_wd = 0
problem.core_O16 = 0.0

problem.X_N14 = 0.01

Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-1.1.M_He-0.05.CO
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 5.e6
problem.temp_core = 1.e7
problem.temp_base = 1.75e8

problem.mixed_co_wd = 1

problem.X_N14 = 0.0

problem.low_density_cutoff = 1.e-4
Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-1.1.M_He-0.05.CO.N14
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 5.e6
problem.temp_core = 1.e7
problem.temp_base = 1.75e8

problem.mixed_co_wd = 1

problem.X_N14 = 0.01

problem.low_density_cutoff = 1.e-4
Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-1.1.M_He-0.05.CO.N14.thin
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 1.e6
problem.temp_core = 1.e7
problem.temp_base = 1.75e8

problem.mixed_co_wd = 1

problem.X_N14 = 0.01

problem.low_density_cutoff = 1.e-4
Expand Down
2 changes: 0 additions & 2 deletions sub_chandra/inputs.M_WD-1.1.M_He-0.05.CO.N14.thinner
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,8 +11,6 @@ problem.delta = 2.e5
problem.temp_core = 1.e7
problem.temp_base = 1.75e8

problem.mixed_co_wd = 1

problem.X_N14 = 0.01

problem.low_density_cutoff = 1.e-4
Expand Down