Skip to content

[WIP] clean-up the transverse flux correction #1920

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
zingale wants to merge 10 commits into
base: development
Choose a base branch
from
Open

[WIP] clean-up the transverse flux correction #1920

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
3a5d173
work in progress
zingale Jun 25, 2021
c3ee56c
first cut
zingale Jun 25, 2021
b1de06d
fix undefined T
zingale Jun 25, 2021
1a32a8e
fix internal energyt
zingale Jun 25, 2021
91667ff
fix 2-d update
zingale Jun 25, 2021
1847e49
Merge branch 'development' into transverse
zingale Jun 25, 2021
0da12e6
Merge branch 'development' into transverse
zingale Jul 27, 2021
e85d45c
Merge branch 'development' into transverse
zingale Sep 1, 2021
07cc554
Merge branch 'development' into transverse
zingale Sep 11, 2021
e341a51
Merge branch 'development' into transverse
zingale May 10, 2022
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
44 changes: 22 additions & 22 deletions Source/hydro/Castro_ctu_hydro.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -569,9 +569,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
vol_arr,
hdt, hdtdy);

reset_edge_state_thermo(xbx, ql.array());
//reset_edge_state_thermo(xbx, ql.array());

reset_edge_state_thermo(xbx, qr.array());
//reset_edge_state_thermo(xbx, qr.array());

// solve the final Riemann problem axross the x-interfaces

Expand Down Expand Up @@ -613,9 +613,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
vol_arr,
hdt, hdtdx);

reset_edge_state_thermo(ybx, ql.array());
//reset_edge_state_thermo(ybx, ql.array());

reset_edge_state_thermo(ybx, qr.array());
//reset_edge_state_thermo(ybx, qr.array());


// solve the final Riemann problem axross the y-interfaces
Expand Down Expand Up @@ -697,9 +697,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdx);

reset_edge_state_thermo(tyxbx, qmyx.array());
//reset_edge_state_thermo(tyxbx, qmyx.array());

reset_edge_state_thermo(tyxbx, qpyx.array());
//reset_edge_state_thermo(tyxbx, qpyx.array());

// [lo(1), lo(2)-1, lo(3)], [hi(1), hi(2)+1, hi(3)+1]
const Box& tzxbx = amrex::grow(zbx, IntVect(AMREX_D_DECL(0,1,0)));
Expand All @@ -725,9 +725,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdx);

reset_edge_state_thermo(tzxbx, qmzx.array());
//reset_edge_state_thermo(tzxbx, qmzx.array());

reset_edge_state_thermo(tzxbx, qpzx.array());
//reset_edge_state_thermo(tzxbx, qpzx.array());

// compute F^y
// [lo(1)-1, lo(2), lo(3)-1], [hi(1)+1, hi(2)+1, hi(3)+1]
Expand Down Expand Up @@ -773,9 +773,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdy);

reset_edge_state_thermo(txybx, qmxy.array());
//reset_edge_state_thermo(txybx, qmxy.array());

reset_edge_state_thermo(txybx, qpxy.array());
//reset_edge_state_thermo(txybx, qpxy.array());

// [lo(1)-1, lo(2), lo(3)], [hi(1)+1, hi(2), lo(3)+1]
const Box& tzybx = amrex::grow(zbx, IntVect(AMREX_D_DECL(1,0,0)));
Expand Down Expand Up @@ -804,9 +804,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdy);

reset_edge_state_thermo(tzybx, qmzy.array());
//reset_edge_state_thermo(tzybx, qmzy.array());

reset_edge_state_thermo(tzybx, qpzy.array());
//reset_edge_state_thermo(tzybx, qpzy.array());

// compute F^z
// [lo(1)-1, lo(2)-1, lo(3)], [hi(1)+1, hi(2)+1, hi(3)+1]
Expand Down Expand Up @@ -852,9 +852,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdz);

reset_edge_state_thermo(txzbx, qmxz.array());
//reset_edge_state_thermo(txzbx, qmxz.array());

reset_edge_state_thermo(txzbx, qpxz.array());
//reset_edge_state_thermo(txzbx, qpxz.array());

// [lo(1)-1, lo(2), lo(3)], [hi(1)+1, hi(2)+1, lo(3)]
const Box& tyzbx = amrex::grow(ybx, IntVect(AMREX_D_DECL(1,0,0)));
Expand Down Expand Up @@ -883,9 +883,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdt, cdtdz);

reset_edge_state_thermo(tyzbx, qmyz.array());
//reset_edge_state_thermo(tyzbx, qmyz.array());

reset_edge_state_thermo(tyzbx, qpyz.array());
//reset_edge_state_thermo(tyzbx, qpyz.array());

// we now have q?zx, q?yx, q?zy, q?xy, q?yz, q?xz

Expand Down Expand Up @@ -946,9 +946,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp2_arr,
hdtdy, hdtdz);

reset_edge_state_thermo(xbx, ql.array());
//reset_edge_state_thermo(xbx, ql.array());

reset_edge_state_thermo(xbx, qr.array());
//reset_edge_state_thermo(xbx, qr.array());

#ifdef SIMPLIFIED_SDC
#ifdef REACTIONS
Expand Down Expand Up @@ -1025,9 +1025,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp1_arr,
hdtdx, hdtdz);

reset_edge_state_thermo(ybx, ql.array());
//reset_edge_state_thermo(ybx, ql.array());

reset_edge_state_thermo(ybx, qr.array());
//reset_edge_state_thermo(ybx, qr.array());

#ifdef SIMPLIFIED_SDC
#ifdef REACTIONS
Expand Down Expand Up @@ -1106,9 +1106,9 @@ Castro::construct_ctu_hydro_source(Real time, Real dt)
qgdnvtmp2_arr,
hdtdx, hdtdy);

reset_edge_state_thermo(zbx, ql.array());
//reset_edge_state_thermo(zbx, ql.array());

reset_edge_state_thermo(zbx, qr.array());
//reset_edge_state_thermo(zbx, qr.array());

#ifdef SIMPLIFIED_SDC
#ifdef REACTIONS
Expand Down
Loading