An official website of the United States government
We are concerned with free boundary problems arising from the analysis of multidimensional transonic shock waves for the Euler equations in compressible fluid dynamics. In this expository paper, we survey some recent developments in the analysis of multidimensional transonic shock waves and corresponding free boundary problems for the compressible Euler equations and related nonlinear partial differential equations (PDEs) of mixed type. The nonlinear PDEs under our analysis include the steady Euler equations for potential flow, the steady full Euler equations, the unsteady Euler equations for potential flow, and related nonlinear PDEs of mixed elliptic–hyperbolic type. The transonic shock problems include the problem of steady transonic flow past solid wedges, the von Neumann problem for shock reflection–diffraction, and the Prandtl–Meyer problem for unsteady supersonic flow onto solid wedges. We first show how these longstanding multidimensional transonic shock problems can be formulated as free boundary problems for the compressible Euler equations and related nonlinear PDEs of mixed type. Then we present an effective nonlinear method and related ideas and techniques to solve these free boundary problems. The method, ideas, and techniques should be useful to analyze other longstanding and newly emerging free boundary problems for nonlinear PDEs.
more »
« less
Loss of Regularity of Solutions of the Lighthill Problem for Shock Diffraction for Potential Flow
We are concerned with the suitability of the main models of compressible fluid dynamics for the Lighthill problem for shock diffraction by a convex corned wedge, by studying the regularity of solutions of the problem, which can be formulated as a free boundary problem. In this paper, we prove that there is no regular solution that is subsonic up to the wedge corner for potential flow. This indicates that, if the solution is subsonic at the wedge corner, at least a characteristic discontinuity (vortex sheet or entropy wave) is expected to be generated, which is consistent with the experimental and computational results. Therefore, the potential flow equation is not suitable for the Lighthill problem so that the compressible Euler system must be considered. In order to achieve the nonexistence result, a weak maximum principle for the solution is established, and several other mathematical techniques are developed. The methods and techniques developed here are also useful to the other problems with similar difficulties.
more »
« less
- Award ID(s):
- 1764278
- PAR ID:
- 10157327
- Date Published:
- Journal Name:
- SIAM journal on mathematical analysis
- Volume:
- 52
- Issue:
- 2
- ISSN:
- 0036-1410
- Page Range / eLocation ID:
- 1096–1114
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Using numerical models, we compute the evolution of the mantle flow field and the crystal preferred orientation (CPO) of mineral aggregates in the mantle wedge of generic subduction systems from their nascent to mature stage and investigate shear wave splitting (SWS) through the forearc mantle wedge corner and overriding crust. Upon subduction initiation, the maximum depth of slab‐mantle decoupling (MDD) is relatively shallow (∼20 km depth), resulting in mantle flow and CPO development in the wedge corner. As subduction continues, the MDD deepens, the wedge corner cools and stagnates, and the olivine CPO becomes frozen‐in. In the cool wedge corner, antigorite can form if water is available. In non‐deforming mantle, antigorite CPO develops relative to the host olivine CPO through topotactic growth. We calculate splitting parameters of synthetic local S waves based on the model‐predicted A‐ and B‐type olivine CPOs and topotactically grown antigorite CPO that replaces A‐type olivine CPO in the wedge corner. The fast direction is trench‐normal for A‐type olivine and antigorite CPOs and trench‐parallel for B‐type. When the delay times are long enough (>0.1 s), we find them positively correlated with the thickness of the mantle wedge corner. In NE Japan, where the results of detailed analyses on the spatial variation of the SWS parameters are available, such correlation is not observationally reported. However, the addition of an anisotropic overriding crust provides delay times (∼0.1 s) and trench‐normal fast directions that are consistent with the local SWS observations.more » « less
-
Direct Numerical Simulation (DNS) of spatially-developing turbulent boundary layers (SDTBL) is performed over isothermal/adiabatic flat plates for incompressible and compressible-subsonic (M∞ = 0.5 and 0.8) flow regimes. Similar low Reynolds numbers are considered in all cases with the purpose of assessing modest flow compressibility on low/high order flow statistics of Zero Pressure Gradient (ZPG) flows. The considered molecular Prandtl number is 0.72. Additionally, temperature is regarded as a passive scalar in the incompressible SDTBL with the purpose to examine differences in the thermal transport phenomena of subsonic flows, i.e., passive vs. active scalar. It was found that the Van Driest transform and Morkovin scaling are able to collapse incompressible and subsonic quantities very well.more » « less
-
When a plane shock hits a two-dimensional wedge head on, it experiences a reflection-diffraction process, and then a self-similar reflected shock moves outward as the original shock moves forward in time. The experimental, computational, and asymptotic analysis has indicated that various patterns occur, including regular reflection and Mach reflection. The von Neumann's conjectures on the transition from regular to Mach reflection involve the existence, uniqueness, and stability of regular shock reflection-diffraction configurations, generated by concave cornered wedges for compressible flow. In this paper, we discuss some recent developments in the study of the von Neumann's conjectures. More specifically, we discuss the uniqueness and stability of regular shock reflection-diffraction configurations governed by the potential flow equation in an appropriate class of solutions. We first show that the transonic shocks in the global solutions obtained in Chen-Feldman [19] are convex. Then we establish the uniqueness of global shock reflection-diffraction configurations with convex transonic shocks for any wedge angle larger than the detachment angle or the critical angle. Moreover, the solution under consideration is stable with respect to the wedge angle. Our approach also provides an alternative way of proving the existence of the admissible solutions established first in [19].more » « less
-
The solution of compressible flow equations is of interest with many aerospace engineering applications. Past literature has focused primarily on the solution of Computational Fluid Dynamics (CFD) problems with low-order finite element and finite volume methods. High-order methods are more the norm nowadays, in both a finite element and a finite volume setting. In this paper, inviscid compressible flow of an ideal gas is solved with high-order spectral/hp stabilized formulations using uniform high-order spectral element methods. The Euler equations are solved with high-order spectral element methods. Traditional definitions of stabilization parameters used in conjunction with traditional low-order bilinear Lagrange-based polynomials provide diffused results when applied to the high-order context. Thus, a revision of the definitions of the stabilization parameters was needed in a high-order spectral/hp framework. We introduce revised stabilization parameters, τsupg, with low-order finite element solutions. We also reexamine two standard definitions of the shock-capturing parameter, δ: the first is described with entropy variables, and the other is the YZβ parameter. We focus on applications with the above introduced stabilization parameters and analyze an array of problems in the high-speed flow regime. We demonstrate spectral convergence for the Kovasznay flow problem in both L1 and L2 norms. We numerically validate the revised definitions of the stabilization parameter with Sod’s shock and the oblique shock problems and compare the solutions with the exact solutions available in the literature. The high-order formulation is further extended to solve shock reflection and two-dimensional explosion problems. Following, we solve flow past a two-dimensional step at a Mach number of 3.0 and numerically validate the shock standoff distance with results obtained from NASA Overflow 2.2 code. Compressible flow computations with high-order spectral methods are found to perform satisfactorily for this supersonic inflow problem configuration. We extend the formulation to solve the implosion problem. Furthermore, we test the stabilization parameters on a complex flow configuration of AS-202 capsule analyzing the flight envelope. The proposed stabilization parameters have shown robustness, providing excellent results for both simple and complex geometries.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
