from pygasflow.interactive.diagrams import DeLavalDiagram
from pygasflow.nozzles import CD_TOP_Nozzle
from pygasflow.solvers import De_Laval_Solver
nozzle = CD_TOP_Nozzle(Ri=0.2, Rt=0.1, Re=0.6, K=0.8)
solver = De_Laval_Solver(
    gamma=1.4, R=287.05, T0=500, P0=101325,
    Pb_P0_ratio=0.2, nozzle=nozzle
)
DeLavalDiagram(solver=solver)