Skip to main content

Simple shock tube calculator

Project description


This module calculates analytical solutions for a Sod shock tube problem.

A more detailed description can be found here.


from shocktubecalc import sod

# left_state and right_state set p, rho and u
# geometry sets left boundary on 0., right boundary on 1
# and initial position of the shock xi on 0.5
# t is the time evolution for which positions and states in tube should be calculated
# gamma denotes specific heat
# npts is number of points to be calculated for rarefaction wave
# Note that gamma and npts are default parameters (1.4 and 500) in solve function.

positions, regions, values = sod.solve(left_state=(1, 1, 0), right_state=(0.1, 0.125, 0.),
                                       geometry=(0., 1., 0.5), t=0.2, gamma=1.4, npts=500)

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

shocktubecalc-0.11.tar.gz (24.8 kB view hashes)

Uploaded source

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring Fastly Fastly CDN Google Google Object Storage and Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page