More Like this

1. Pressure-driven gas ow in viscously deformable porous media: application to lava domes

   https://doi.org/10.1017/jfm2019-211  

   Hyman, D. ; Bursik, M. ; Pitman, E.B ( June 2019 , Journal of fluid mechanics)

   The behaviour of low-viscosity, pressure-driven compressible pore fluid flows in viscously deformable porous media is studied here with specific application to gas flow in lava domes. The combined flow of gas and lava is shown to be governed by a two-equation set of nonlinear mixed hyperbolic–parabolic type partial differential equations describing the evolution of gas pore pressure and lava porosity. Steady state solution of this system is achieved when the gas pore pressure is magmastatic and the porosity profile accommodates the magmastatic pressure condition by increased compaction of the medium with depth. A one-dimensional (vertical) numerical linear stability analysis (LSA) is presented here. As a consequence of the pore-fluid compressibility and the presence of gravitation compaction, the gradients present in the steady-state solution cause variable coefficients in the linearized equations which generate instability in the LSA despite the diffusion-like and dissipative terms in the original system. The onset of this instability is shown to be strongly controlled by the thickness of the flow and the maximum porosity, itself a function of the mass flow rate of gas. Numerical solutions of the fully nonlinear system are also presented and exhibit nonlinear wave propagation features such as shock formation. As applied to gas flow within lava domes, the details of this dynamics help explain observations of cyclic lava dome extrusion and explosion episodes. Because the instability is stronger in thicker flows, the continued extrusion and thickening of a lava dome constitutes an increasing likelihood of instability onset, pressure wave growth and ultimately explosion. 

   more » « less
2. The relativistic Euler equations: ESI notes on their geo-analytic structures and implications for shocks in 1D and multi-dimensions

   https://doi.org/10.1088/1361-6382/ad059a  

   Abbrescia, Leonardo ; Speck, Jared ( November 2023 , Classical and Quantum Gravity)

   In this article, we provide notes that complement the lectures on the relativistic Euler equations and shocks that were given by the second author at the programMathematical Perspectives of Gravitation Beyond the Vacuum Regime, which was hosted by the Erwin Schrödinger International Institute for Mathematics and Physics in Vienna in February 2022. We set the stage by introducing a standard first-order formulation of the relativistic Euler equations and providing a brief overview of local well-posedness in Sobolev spaces. Then, using Riemann invariants, we provide the first detailed construction of a localized subset of the maximal globally hyperbolic developments of an open set of initially smooth, shock-forming isentropic solutions in 1D, with a focus on describing the singular boundary and the Cauchy horizon that emerges from the singularity. Next, we provide an overview of the new second-order formulation of the 3Drelativistic Euler equations derived in Disconzi and Speck (2019Ann. Henri Poincare202173–270), its rich geometric and analytic structures, their implications for the mathematical theory of shock waves, and their connection to the setup we use in our 1Danalysis of shocks. We then highlight some key prior results on the study of shock formation and related problems. Furthermore, we provide an overview of how the formulation of the flow derived in Disconzi and Speck (2019Ann. Henri Poincare202173–270) can be used to study shock formation in multiple spatial dimensions. Finally, we discuss various open problems tied to shocks.

    

   more » « less
3. Spatially quasi-periodic water waves of finite depth

   https://doi.org/10.1098/rspa.2023.0019  

   Wilkening, Jon ; Zhao, Xinyu ( April 2023 , Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)

   We present a numerical study of spatially quasi-periodic gravity-capillary waves of finite depth in both the initial value problem and travelling wave settings. We adopt a quasi-periodic conformal mapping formulation of the Euler equations, where one-dimensional quasi-periodic functions are represented by periodic functions on a higher-dimensional torus. We compute the time evolution of free surface waves in the presence of a background flow and a quasi-periodic bottom boundary and observe the formation of quasi-periodic patterns on the free surface. Two types of quasi-periodic travelling waves are computed: small-amplitude waves bifurcating from the zero-amplitude solution and larger-amplitude waves bifurcating from finite-amplitude periodic travelling waves. We derive weakly nonlinear approximations of the first type and investigate the associated small-divisor problem. We find that waves of the second type exhibit striking nonlinear behaviour, e.g. the peaks and troughs are shifted non-periodically from the corresponding periodic waves due to the activation of quasi-periodic modes.

    

   more » « less
4. Uniqueness and Stability for the Shock Reflection-Diffraction Problem for Potential Flow

   Chen, Gui-Qiang G ; Feldman, Mikhail ; Xiang, Wei ( January 2019 , Proceedings of the XVII international conference on hyperbolic problems: Theory, numerics, applications)

   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
5. A Spectral/hp-Based Stabilized Solver with Emphasis on the Euler Equations

   https://doi.org/10.3390/fluids9010018  

   Ranjan, Rakesh ; Catabriga, Lucia ; Araya, Guillermo ( January 2024 , Fluids)

   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