Numerical Modeling of Phase Transition in Shock Tube Problem
In this work, we developed a numerical model for phase transition in shock-tube problem. This work focused on simulating shock-tube problem and the numerical framework is based on a two-phase flow model. Our simulation tool is an in-house FORTRAN multiphase compressible flow solver, RocSDT . Specifically, this involves the Euler equations along with the Stiffened Gas Equation of State (EOS), an interface of liquid-gas type, and the effects of heat and mass transfer due to phase transition. This thesis described the implementation and incorporation of a phase transition model into a system of five partial differential equation model and the associated thermodynamic theory. The thermal and chemical source terms for phase transition modeling are handled via a sequence of relaxation processes that drive the mixture to the desired equilibrium conditions. In particular, a simple algebraic system of equations is used to solve for the equilbrium conditions. This algebraic relaxation technique along with our five-equation model are shown to be simple and effective.