SUMMARY The occurrence of plate tectonics on Earth is rooted in the physics of lithospheric ductile weakening and shear-localization. The pervasiveness of mylonites at lithospheric shear zones is a key piece of evidence that localization correlates with reduction in mineral grain size. Most lithospheric mylonites are polymineralic and the interaction between mineral phases, such as olivine and pyroxene, especially through Zener pinning, impedes normal grain growth while possibly enhancing grain damage, both of which facilitate grain size reduction and weakening, as evident in lab experiments and field observations. The efficacy of pinning, however, relies on the mineral phases being mixed and dispersed at the grain scale, where well-mixed states lead to greater mylonitization. To model grain mixing between different phases at the continuum scale, we previously developed a theory treating grain-scale processes as diffusion between phases, but driven by imposed compressive stresses acting on the boundary between phases. Here we present a new model for shearing rock that combines our theory for diffusive grain mixing, 2-D non-Newtonian flow and two-phase grain damage. The model geometry is designed specifically for comparison to torsional shear-deformation experiments. Deformation is either forced by constant velocity or constant stress boundary conditions. As the layer is deformed, mixing zones between different mineralogical units undergo enhanced grain size reduction and weakening, especially at high strains. For constant velocity boundary experiments, stress drops towards an initial piezometric plateau by a strain of around 4; this is also typical of monophase experiments for which this initial plateau is the final steady state stress. However, polyphase experiments can undergo a second large stress drop at strains of 10–20, and which is associated with enhanced phase mixing and resultant grain size reduction and weakening. Model calculations for polyphase media with grain mixing and damage capture the experimental behaviour when damage to the interface between phases is moderately slower or less efficient than damage to the grain boundaries. Other factors such as distribution and bulk fraction of the secondary phase, as well as grain-mixing diffusivity also influence the timing of the second stress drop. For constant stress boundary conditions, the strain rate increases during weakening and localization. For a monophase medium, there is theoretically one increase in strain rate to a piezometric steady state. But for the polyphase model, the strain rate undergoes a second abrupt increase, the timing for which is again controlled by interface damage and grain mixing. The evolution of heterogeneity through mixing and deformation, and that of grain size distributions also compare well to experimental observations. In total, the comparison of theory to deformation experiments provides a framework for guiding future experiments, scaling microstructural physics to geodynamic applications and demonstrates the importance of grain mixing and damage for the formation of plate tectonic boundaries.
more »
« less
Hydraulic control of flow in a multi-passage system connecting two basins
When a fluid stream in a conduit splits in order to pass around an obstruction, it is possible that one branch will be critically controlled while the other remains not so. This is apparently the situation in Pacific Ocean abyssal circulation, where most of the northward flow of Antarctic bottom water passes through the Samoan Passage, where it is hydraulically controlled, while the remainder is diverted around the Manihiki Plateau and is not controlled. These observations raise a number of questions concerning the dynamics necessary to support such a regime in the steady state, the nature of upstream influence and the usefulness of rotating hydraulic theory to predict the partitioning of volume transport between the two paths, which assumes the controlled branch is inviscid. Through the use of a theory for constant potential vorticity flow and accompanying numerical model, we show that a steady-state regime similar to what is observed is dynamically possible provided that sufficient bottom friction is present in the uncontrolled branch. In this case, the upstream influence that typically exists for rotating channel flow is transformed into influence into how the flow is partitioned. As a result, the partitioning of volume flux can still be reasonably well predicted with an inviscid theory that exploits the lack of upstream influence.
more »
« less
- Award ID(s):
- 1658027
- NSF-PAR ID:
- 10346122
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 940
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We vary the inflow properties in a finite-volume solver to investigate their effects on the computed cyclonic motion in a right-cylindrical vortex chamber. The latter comprises eight tangential injectors through which steady-state air is introduced under incompressible and inviscid conditions. To minimize cell skewness around injectors, a fine tetrahedral mesh is implemented first and then converted into polyhedral elements, namely, to improve convergence characteristics and precision. Once convergence is achieved, our principal variables are evaluated and compared using a range of inflow parameters. These include the tangential injector speed, count, diameter, and elevation. The resulting computations show that well-resolved numerical simulations can properly predict the forced vortex behavior that dominates in the core region as well as the free vortex tail that prevails radially outwardly, beyond the point of peak tangential speed. It is also shown that augmenting the mass influx by increasing the number of injectors, injector size, or average injection speed further amplifies the vortex strength and all peak velocities while shifting the mantle radially inwardly. Overall, the axial velocity is found to be the most sensitive to vertical displacements of the injection plane. By raising the injection plane to the top half portion of the chamber, the flow character is markedly altered, and an axially unidirectional vortex is engendered, particularly, with no upward motion or mantle formation. Conversely, the tangential and radial velocities are found to be axially independent and together with the pressure distribution prove to be the least sensitive to injection plane relocations.more » « less
-
We vary the inflow properties in a finite-volume solver to investigate their effects on the computed cyclonic motion in a right-cylindrical vortex chamber. The latter comprises eight tangential injectors through which steady-state air is introduced under incompressible and inviscid conditions. To minimize cell skewness around injectors, a fine tetrahedral mesh is implemented first and then converted into polyhedral elements, namely, to improve convergence characteristics and precision. Once convergence is achieved, our principal variables are evaluated and compared using a range of inflow parameters. These include the tangential injector speed, count, diameter, and elevation. The resulting computations show that well-resolved numerical simulations can properly predict the forced vortex behavior that dominates in the core region as well as the free vortex tail that prevails radially outwardly, beyond the point of peak tangential speed. It is also shown that augmenting the mass influx by increasing the number of injectors, injector size, or average injection speed, further amplifies the vortex strength and all peak velocities while shifting the mantle radially inwardly. Overall, the axial velocity is found to be the most sensitive to vertical displacements of the injection plane. By raising the injection plane to the top half portion of the chamber, the flow character is markedly altered, and an axially unidirectional vortex is engendered, particularly, with no upward motion or mantle formation. Conversely, the tangential and radial velocities are found to be axially independent and together with the pressure distribution prove to be the least sensitive to injection plane relocations.more » « less
-
more » « less
Abstract The aquatic bladderwort
Utricularia gibba captures zooplankton in mechanically triggered underwater traps. With characteristic dimensions <1 mm, the trapping structures are among the smallest known that work by suction—a mechanism that would not be effective in the creeping‐flow regime. To understand the adaptations that make suction feeding possible on this small scale, we have measured internal flow speeds during artificially triggered feeding strikes in the absence of prey. These data are compared with complementary analytical models of the suction event: an inviscid model of the jet development in time and a steady‐state model incorporating friction. The initial dynamics are well described by a time‐dependent Bernoulli equation in which the action of the trap door is represented by a step increase in driving pressure. According to this model, the observed maximum flow speed (5.2 m/s) depends only on the pressure difference, whereas the initial acceleration (3 × 104 m/s2) is determined by pressure difference and channel length. Because the terminal speed is achieved quickly (~0.2 ms) and the channel is short, the remainder of the suction event (~2.0 ms) is effectively an undeveloped viscous steady state. The steady‐state model predicts that only 17% of power is lost to friction. The energy efficiency and steady‐state fluid speed decrease rapidly with decreasing channel diameter, setting a lower limit on practical bladderwort size. -
null (Ed.)This work focuses on the use of a finite-volume solver to describe the wall-bounded cyclonic flowfield that evolves in a swirl-driven thrust chamber. More specifically, a non-reactive, cold-flow simulation is carried out using an idealized chamber configuration of a square-shaped, right-cylindrical enclosure with eight tangential injectors and a variable nozzle size. For simplicity, we opt for air as the working fluid and perform our simulations under steady, incompressible, and inviscid flow conditions. First, a meticulously developed mesh that consists of tetrahedral elements is generated in a manner to minimize the overall skewness, especially near injectors; second, this mesh is converted into a polyhedral grid to improve convergence characteristics and accuracy. After achieving convergence in all variables, our three velocity components are examined and compared to an existing analytical solution obtained by Vyas and Majdalani (Vyas, A. B., and Majdalani, J., “Exact Solution of the Bidirectional Vortex,” AIAA Journal, Vol. 44, No. 10, 2006, pp. 2208-2216). We find that the numerical model is capable of predicting the expected forced vortex behavior in the inner core region as well as the free vortex tail in the inviscid region. Moreover, the results appear to be in fair agreement with the Vyas–Majdalani solution derived under similarly inviscid conditions, and thus resulting in a quasi complex-lamellar profile. In this work, we are able to ascertain the axial independence of the swirl velocity no matter the value of the outlet radius, which confirms the key assumption made in most analytical models of wall-bounded cyclonic motions. Moreover, the pressure distribution predicted numerically is found to be in fair agreement with both theoretical formulations and experimental measurements of cyclone separators. The bidirectional character of the flowfield is also corroborated by the axial and radial velocity distributions, which are found to be concurrent with theory. Then using parametric trade studies, the sensitivity of the numerical simulations to the outlet diameter of the chamber is explored to determine the influence of outlet nozzle variations on the mantle location and the number of mantles. Since none of the cases considered here promote the onset of multiple mantles, we are led to believe that more factors are involved in producing more mantles than one. Besides the exit diameter, the formation of a multiple mantle structure may be influenced by the physical boundary conditions, nozzle radius, inlet curvature, and length. In closing, we show that the latter plays a significant role in controlling the development of backflow regions inside the chamber.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2022.212
