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.

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

Jump to bottom

update mom_flux_has_p to address with 2D spherical #2958

Merged
merged 11 commits into from
Sep 16, 2024
Merged

update mom_flux_has_p to address with 2D spherical #2958

Changes from 1 commit
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
3395e2e
update mom_flux_has_p in Castro_util
zhichen3 Sep 11, 2024
ba68c91
Merge branch 'development' into 2d_spherical_mom_flux_has_p
zhichen3 Sep 11, 2024
79562a8
change places associated with mom_flux_has_p except castro_ctu/mol_hy…
zhichen3 Sep 11, 2024
ca04f9e
fix incorrect calculation
zhichen3 Sep 11, 2024
442b5d8
remove redundant else and fix CoordSys::spherical -> CoordSys:SPHERICAL
zhichen3 Sep 11, 2024
8606de8
missing colon
zhichen3 Sep 11, 2024
3631380
fix typo
zhichen3 Sep 11, 2024
f9df3a7
add preprocessor to ensure we're in 2d or more
zhichen3 Sep 11, 2024
727f49f
fix compilation
zhichen3 Sep 11, 2024
0bd9707
fix typo
zhichen3 Sep 12, 2024
abe9edd
Merge branch 'development' into 2d_spherical_mom_flux_has_p
zingale Sep 16, 2024
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
6 changes: 3 additions & 3 deletions Source/hydro/Castro_ctu.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,14 +69,14 @@ Castro::consup_hydro(const Box& bx,
// Add gradp term to momentum equation -- only for axisymmetric
// coords (and only for the radial flux).

U_new(i,j,k,UMX) += - dt * (qx(i+1,j,k,GDPRES) - qx(i,j,k,GDPRES)) / geomdata.CellSize()[0];
U_new(i,j,k,UMX) += - dt * (qx(i+1,j,k,GDPRES) - qx(i,j,k,GDPRES)) / geomdata.CellSize(1);
zingale marked this conversation as resolved.
Show resolved Hide resolved

#if AMREX_SPACEDIM >= 2
} else if (n == UMY && !mom_flux_has_p(1, 1, geomdata.Coord())) {
// Add gradp term to polar(theta) momentum equation for Spherical 2D geometry

Real r = geom.ProbLoArray()[0] + (static_cast<Real>(i) + 0.5_rt) * geomdata.CellSize()[0];
U_new(i,j,k,UMY) += - dt * (qy(i,j+1,k,GDPRES) - qy(i,j,k,GDPRES)) / (r * geomdata.CellSize()[1]);
Real r = geomdata.ProbLo(0) + (static_cast<Real>(i) + 0.5_rt) * geomdata.CellSize(0);
U_new(i,j,k,UMY) += - dt * (qy(i,j+1,k,GDPRES) - qy(i,j,k,GDPRES)) / (r * geomdata.CellSize(1));
#endif
}
});
Expand Down