compressible_rk package¶
A method-of-lines compressible hydrodynamics solver. Piecewise constant reconstruction is done in space and a Runge-Kutta time integration is used to advance the solutiion.
Subpackages¶
Submodules¶
compressible_rk.fluxes module¶
This is a 2nd-order PLM method for a method-of-lines integration (i.e., no characteristic tracing).
We wish to solve
we want U_{i+1/2} – the interface values that are input to the Riemann problem through the faces for each zone.
Taylor expanding in space only yields:
dU
U = U + 0.5 dx --
i+1/2,j,L i,j dx
- compressible_rk.fluxes.fluxes(my_data, rp, ivars, solid, tc)[source]¶
unsplitFluxes returns the fluxes through the x and y interfaces by doing an unsplit reconstruction of the interface values and then solving the Riemann problem through all the interfaces at once
currently we assume a gamma-law EOS
- Parameters
- my_dataCellCenterData2d object
The data object containing the grid and advective scalar that we are advecting.
- rpRuntimeParameters object
The runtime parameters for the simulation
- varsVariables object
The Variables object that tells us which indices refer to which variables
- tcTimerCollection object
The timers we are using to profile
- Returns
- outndarray, ndarray
The fluxes on the x- and y-interfaces
compressible_rk.simulation module¶
- class compressible_rk.simulation.Simulation(solver_name, problem_name, rp, timers=None, data_class=<class 'mesh.patch.CellCenterData2d'>)[source]¶
Bases:
compressible.simulation.SimulationThe main simulation class for the method of lines compressible hydrodynamics solver