An official website of the United States government
Abstract Two‐phase flow, a system where Stokes flow and Darcy flow are coupled, is of great importance in the Earth's interior, such as in subduction zones, mid‐ocean ridges, and hotspots. However, it remains challenging to solve the two‐phase equations accurately in the zero‐porosity limit, for example, when melt is fully frozen below solidus temperature. Here we propose a new three‐field formulation of the two‐phase system, with solid velocity (vs), total pressure (Pt), and fluid pressure (Pf) as unknowns, and present a robust finite‐element implementation, which can be used to solve problems in which domains of both zero porosity and non‐zero porosity are present. The reformulated equations include regularization to avoid singularities and exactly recover to the standard single‐phase incompressible Stokes problem at zero porosity. We verify the correctness of our implementation using the method of manufactured solutions and analytic solutions and demonstrate that we can obtain the expected convergence rates in both space and time. Example experiments, such as self‐compaction, falling block, and mid‐ocean ridge spreading show that this formulation can robustly resolve zero‐ and non‐zero‐porosity domains simultaneously, and can be used for a large range of applications in various geodynamic settings.
more »
« less
Pressure-driven gas ow in viscously deformable porous media: application to lava domes
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
- PAR ID:
- 10111295
- Date Published:
- Journal Name:
- Journal of fluid mechanics
- Volume:
- 869
- ISSN:
- 1469-7645
- Page Range / eLocation ID:
- 85-109
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The interfacial instability in a two-phase mixing layers between parallel gas and liquid streams is important to two-phase atomization. Depending on the inflow conditions and fluid properties, interfacial instability can be convective or absolute. The goal of the present study is to investigate the impact of gas viscosity on the interfacial instability. Both interface-resolved simulations and linear stability analysis (LSA) have been conducted. In LSA, the Orr–Sommerfeld equation is solved to analyze the spatio-temporal viscous modes. When the gas viscosity decreases, the Reynold number (Re) increases accordingly. The LSA demonstrates that when Re is higher than a critical threshold, the instability transitions from the absolute to the convective (A/C) regimes. Such a Re-induced A/C transition is also observed in the numerical simulations, though the critical Re observed in simulations is significantly lower than that predicted by LSA. The LSA results indicate that the temporal growth rate decreases with Re. When the growth rate reaches zero, the A/C transition will occur. The Re-induced A/C transition is observed in both confined and unconfined mixing layers and also in cases with low and high gas-to-liquid density ratios. In the transition from typical absolute and convective regimes, a weak absolute regime is identified in the simulations, for which the spectrograms show both the absolute and convective modes. The dominant frequency in the weak absolute regime can be influenced by the perturbation introduced at the inlet. The simulation results also show that the wave propagation speed can vary in space. In the absolute instability regime, the wave propagation speed agrees well with the absolute mode celerity near the inlet and increases to the Dimotakis speed further downstream.more » « less
-
Abstract. Fluid and melt transport in the solid mantle can be modeled as a two-phase flow in which the liquid flow is resisted by the compaction of the viscously deforming solid mantle. Given the wide impact of liquid transport on the geodynamical and geochemical evolution of the Earth, the so-called “compaction equations” are increasingly being incorporated into geodynamical modeling studies. When implementing these equations, it is common to use a regularization technique to handle the porosity singularity in the dry mantle. Moreover, it is also common to enforce a positive porosity (liquid fraction) to avoid unphysical negative values of porosity. However, the effects of this “capped” porosity on the liquid flow and mass conservation have not been quantitatively evaluated. Here, we investigate these effects using a series of 1- and 2-dimensional numerical models implemented using the commercial finite-element package COMSOL Multiphysics®. The results of benchmarking experiments against a semi-analytical solution for 1- and 2-D solitary waves illustrate the successful implementation of the compaction equations. We show that the solutions are accurate when the element size is smaller than half of the compaction length. Furthermore, in time-evolving experiments where the solid is stationary (immobile), we show that the mass balance errors are similarly low for both the capped and uncapped (i.e., allowing negative porosity) experiments. When Couette flow, convective flow, or subduction corner flow of the solid mantle is assumed, the capped porosity leads to overestimations of the mass of liquid in the model domain and the mass flux of liquid across the model boundaries, resulting in intrinsic errors in mass conservation even if a high mesh resolution is used. Despite the errors in mass balance, however, the distributions of the positive porosity and peaks (largest positive liquid fractions) in both the uncapped and capped experiments are similar. Hence, the capping of porosity in the compaction equations can be reasonably used to assess the main pathways and first-order distribution of fluids and melts in the mantle.more » « less
-
Purpose The purpose of this study is to develop nonlinear and linearized models of DW printing dynamics that capture the complexity of DW while remaining integrable into control schemes. Control of material metering in extrusion-based additive manufacturing modalities, such as positive displacement direct-write (DW), is critical for manufacturing accuracy. However, in DW, transient flows are poorly controlled due to capacitive pressure dynamics – pressure is stored and slowly released over time from the build material and other compliant system elements, adversely impacting flow rate start-ups and stops. Thus far, modeling of these dynamics has ranged from simplistic, potentially omitting key contributors to the observed phenomena, to highly complex, making usage in control schemes difficult. Design/methodology/approach The authors present nonlinear and linearized models that seek to both capture the capacitive and nonlinear resistive fluid elements of DW systems and to pose them as ordinary differential equations for integration into control schemes. The authors validate the theoretical study with experimental flow rate and material measurements across a range of extrusion nozzle sizes and materials. The authors explore the contribution of the system and build material bulk modulus to these dynamics. Findings The authors show that all tested models accurately describe the measured dynamics, facilitating ease of integration into future control systems. Additionally, the authors show that system bulk modulus may be substantially reduced through appropriate system design. However, the remaining build material bulk modulus is sufficient to require feedback control for accurate material delivery. Originality/value This study presents new nonlinear and linear models for DW printing dynamics. The authors show that linear models are sufficient to describe the dynamics, with small errors between nonlinear and linear models. The authors demonstrate control is necessary for accurate material delivery in DW.more » « less
-
Abstract Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm2019-211
