pyro.compressible_rk.problems package

pyro.compressible_rk.problems package#

Submodules#

pyro.compressible_rk.problems.acoustic_pulse module#

The acoustic pulse problem described in McCorquodale & Colella 2011. This uses a uniform background and a small pressure perturbation that drives a low Mach number soundwave. This problem is useful for testing convergence of a compressible solver.

pyro.compressible_rk.problems.acoustic_pulse.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.acoustic_pulse.init_data(myd, rp)[source]#: initialize the acoustic_pulse problem. This comes from McCorquodale & Coella 2011

pyro.compressible_rk.problems.advect module#

A simple advection test. A density perturbation is set with a constant pressure in the domain and a velocity field is set to advect the profile across the domain. This is useful for testing convergence.

pyro.compressible_rk.problems.advect.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.advect.init_data(my_data, rp)[source]#: initialize a smooth advection problem for testing convergence

pyro.compressible_rk.problems.bubble module#

A buoyant perturbation (bubble) is placed in an isothermal hydrostatic atmosphere (plane-parallel). It will rise and deform (due to shear)

pyro.compressible_rk.problems.bubble.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.bubble.init_data(my_data, rp)[source]#: initialize the bubble problem

pyro.compressible_rk.problems.gresho module#

The Gresho vortex problem sets up a toroidal velocity field that is balanced by a radial pressure gradient. This is in equilibrium and the state should remain unchanged in time. This version of the problem is based on Miczek, Roepke, and Edelmann 2014.

pyro.compressible_rk.problems.gresho.finalize()[source]#: print out any information to the userad at the end of the run

pyro.compressible_rk.problems.gresho.init_data(my_data, rp)[source]#: initialize the Gresho vortex problem

pyro.compressible_rk.problems.hse module#

Initialize an isothermal hydrostatic atmosphere. It should remain static. This is a test of our treatment of the gravitational source term.

pyro.compressible_rk.problems.hse.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.hse.init_data(my_data, rp)[source]#: initialize the HSE problem

pyro.compressible_rk.problems.kh module#

A Kelvin-Helmholtz shear problem. There are 2 shear layers, with the and an optional vertical bulk velocity. This can be used to test the numerical dissipation in the solver. This setup is based on McNally et al. 2012.

pyro.compressible_rk.problems.kh.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.kh.init_data(my_data, rp)[source]#: initialize the Kelvin-Helmholtz problem

pyro.compressible_rk.problems.logo module#

Generate the pyro logo! The word “pyro” is written in the center of the domain and perturbations are placed in the 4 corners to drive converging shocks inward to scramble the logo.

pyro.compressible_rk.problems.logo.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.logo.init_data(my_data, rp)[source]#: initialize the logo problem

pyro.compressible_rk.problems.quad module#

The quadrant problem from Shulz-Rinne et al. 1993; Lax and Lui 1998. Four different states are initialized in the quadrants of the domain, driving shocks and other hydrodynamic waves at the interfaces. This can be used to test the symmetry of the solver.

pyro.compressible_rk.problems.quad.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.quad.init_data(my_data, rp)[source]#: initialize the quadrant problem

pyro.compressible_rk.problems.ramp module#

A shock hitting a ramp at an oblique angle. This is based on Woodward & Colella 1984.

pyro.compressible_rk.problems.ramp.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.ramp.init_data(my_data, rp)[source]#: initialize the double Mach reflection problem

pyro.compressible_rk.problems.rt module#

A single-mode Rayleigh-Taylor instability.

pyro.compressible_rk.problems.rt.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.rt.init_data(my_data, rp)[source]#: initialize the rt problem

pyro.compressible_rk.problems.rt2 module#

A RT problem with two distinct modes: short wavelength on the left and long wavelength on the right. This allows one to see how the growth rate depends on wavenumber.

pyro.compressible_rk.problems.rt2.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.rt2.init_data(my_data, rp)[source]#: initialize the rt problem

pyro.compressible_rk.problems.sedov module#

The classic Sedov problem.

pyro.compressible_rk.problems.sedov.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.sedov.init_data(my_data, rp)[source]#: initialize the sedov problem

pyro.compressible_rk.problems.sod module#

A general shock tube problem for comparing the solver to an exact Riemann solution.

pyro.compressible_rk.problems.sod.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.sod.init_data(my_data, rp)[source]#: initialize the sod problem

pyro.compressible_rk.problems.test module#

A setup intended for unit testing.

pyro.compressible_rk.problems.test.finalize()[source]#: print out any information to the user at the end of the run

pyro.compressible_rk.problems.test.init_data(my_data, rp)[source]#: an init routine for unit testing