The Finite Element Numerical Study of Transonic Flow of Moist Air with Nonequilibrium Condensation

Abstract

The finite element numerical simulation of the transonic flow with condensation is studied in this paper. The transformed weak form is numerically discretized by using the classical Galerkin finite element spatial discretization scheme and the backward Euler difference time discretization scheme. Because the coupled Euler equations are convection equations, the results obtained by the classical Galerkin method often have nonphysical numerical oscillations, so the stabilization method is introduced. Another important feature of compressible Euler equation is the existence of shock wave. Because the thickness of shock wave is small and difficult to capture, the shock wave can be better captured by adding artificial viscosity term. Here, the streamline upwind/Petrov–Galerkin (SUPG) finite element method with good stability and the isotropic diffusion method with shock capture ability are used to solve the problem of transonic flow with condensation. In the simulation of steady problems, in order to improve the calculation accuracy of the existing area of the shock wave, the adaptive mesh refinement technology is added, which refines the mesh in the region of the shock wave to improve the resolution of the shock wave. Numerical experiments verify the feasibility and stabilization of the numerical method. Finally, the simulation of transonic flow around NACA-0012 and parallel-jet nozzle A1 has obtained the better numerical results.