shocktube1dcalc 0.0.4

1D shocktube caculator to provide analytic solutions

1D shocktube caculator

This tool provdes 1D Shock Tube analytic solutions.

Getting Started

Prerequisites

• Python
• numpy
• scipy

Usage

from shocktube1dcalc import solver_analytic

# by default it will create a the shock tube based on Sod's classic condition.
shocktube = solver_analytic.ShockTube()

import numpy as np
mesh = np.linspace(-0.5, 0.5, 50)

analytic_solution = shocktube.get_analytic_solution(
    mesh, t=0.4
)

You may customize the physical status of the shocktube via:

shocktube = solver_analytic.ShockTube(rho_left=1.0, u_left=0.0, p_left=1.0, rho_right=0.125, u_right=0.0, p_right=0.1)

Contributing

See Contributing

Authors

Taihsiang Ho (tai271828) [e-mail]

Created from Lee-W/cookiecutter-python-template version 0.7.1