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_data : CellCenterData2d object
The data object containing the grid and advective scalar that we are advecting.
- rp : RuntimeParameters object
The runtime parameters for the simulation
- vars : Variables object
The Variables object that tells us which indices refer to which variables
- tc : TimerCollection object
The timers we are using to profile
Returns: - out : ndarray, 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