Castro

Castro copied to clipboard

In the job_info file, we report the CPU time, e.g.: ``` CPU time used since start of simulation (CPU-hours): 181.756 ``` This is just # of MPI * # of...

zingale

SDC/MOL doesn't work with multilevel

3

we need to come up with a scheme to do the ghost cell fill for the intermediate stages (which is not trivial with RK integration). we also need to address...

zingale

axisymmetric in 4th order

1

There are a number of changes we need to do to support 4th order in 2-d axisymmetric coords: * the transformation from cell-centers to cell-averages changes * the interpolant we...

zingale

For the diffusion to 4th order, we can just set the diffusive flux on the domain boundary to 0, due to symmetry.

zingale

derived variables are not high-order

2

The derived variables in the plotfiles are only computed to second-order accuracy. We may want to handle this via AMReX directly.

zingale

Implement optimal subcycling

1

Optimal subcycling is done in Nyx, and I want to implement this as an option in Castro. This involves two steps: 1) during timestep estimation, call the Amr function which...

maxpkatz

allow radial sponge to work in plane parallel geometry

1

For plane parallel geometries, the radial sponge should only depend on the vertical direction.

zingale

We should look into adding subsonic outflow boundaries. Maybe we can adapt the code from here: https://ccse.lbl.gov/people/motheau/Manuscripts_website/2017_AIAA_CFD_Motheau.pdf

zingale

Is FillPatch needed in diffusion?

3

Currently in the construction of the diffusion term, we do a FillPatch on the full state. But for the predictor (with CTU) or with any MOL-based hydro, we come into...

zingale

compute average RHS calls for reactions

3

It would be a useful diagnostic to compute and store / report the average number of RHS steps taken by the reaction network integrator each step.

zingale