scholarly journals Image-based modelling of lateral magma flow: the Basement Sill, Antarctica

2017 ◽  
Vol 4 (5) ◽  
pp. 161083 ◽  
Author(s):  
Nick Petford ◽  
Seyed Mirhadizadeh
Keyword(s):  
Numerical Models ◽  
Mass Density ◽  
Three Dimensions ◽  
Dry Valleys ◽  
Element Mesh ◽  
Magma Flow ◽  
Continental Scale ◽  
Low Viscosity ◽  
High Particle ◽  
High Reynolds

The McMurdo Dry Valleys magmatic system, Antarctica, provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle microstructure of a congested magma slurry. We simulated the flow regime in two and three dimensions using numerical models built on a finite-element mesh derived from field data. The model captures the flow behaviour of the Basement Sill magma over a viscosity range of 1–10 4  Pa s where the higher end (greater than or equal to 10 2  Pa s) corresponds to a magmatic slurry with crystal fractions varying between 30 and 70%. A novel feature of the model is the discovery of transient, low viscosity (less than or equal to 50 Pa s) high Reynolds number eddies formed along undulating contacts at the floor and roof of the intrusion. Numerical tracing of particle orbits implies crystals trapped in eddies segregate according to their mass density. Recovered shear strain rates (10 −3 –10 −5  s −1 ) at viscosities equating to high particle concentrations (around more than 40%) in the Sill interior point to shear-thinning as an explanation for some types of magmatic layering there. Model transport rates for the Sill magmas imply a maximum emplacement time of ca 10 5 years, consistent with geochemical evidence for long-range lateral flow. It is a theoretically possibility that fast-flowing magma on a continental scale will be susceptible to planetary-scale rotational forces.

Ice-Phase Precipitation

2017 ◽  
Vol 58 ◽  
pp. 6.1-6.36 ◽  
Author(s):  
I. Gultepe ◽  
A. J. Heymsfield ◽  
P. R. Field ◽  
D. Axisa
Keyword(s):  
Collection Efficiency ◽  
Numerical Models ◽  
Mass Density ◽  
Three Dimensional ◽  
Phase Precipitation ◽  
Microphysical Parameters ◽  
Mass Size ◽  
Ice Water ◽  
Ice Phase

AbstractIce-phase precipitation occurs at Earth’s surface and may include various types of pristine crystals, rimed crystals, freezing droplets, secondary crystals, aggregates, graupel, hail, or combinations of any of these. Formation of ice-phase precipitation is directly related to environmental and cloud meteorological parameters that include available moisture, temperature, and three-dimensional wind speed and turbulence, as well as processes related to nucleation, cooling rate, and microphysics. Cloud microphysical parameters in the numerical models are resolved based on various processes such as nucleation, mixing, collision and coalescence, accretion, riming, secondary ice particle generation, turbulence, and cooling processes. These processes are usually parameterized based on assumed particle size distributions and ice crystal microphysical parameters such as mass, size, and number and mass density. Microphysical algorithms in the numerical models are developed based on their need for applications. Observations of ice-phase precipitation are performed using in situ and remote sensing platforms, including radars and satellite-based systems. Because of the low density of snow particles with small ice water content, their measurements and predictions at the surface can include large uncertainties. Wind and turbulence affecting collection efficiency of the sensors, calibration issues, and sensitivity of ground-based in situ observations of snow are important challenges to assessing the snow precipitation. This chapter’s goals are to provide an overview for accurately measuring and predicting ice-phase precipitation. The processes within and below cloud that affect falling snow, as well as the known sources of error that affect understanding and prediction of these processes, are discussed.

2-D and 3-D resistivity image reconstruction using crosshole data

Geophysics ◽  
1992 ◽  
Vol 57 (10) ◽  
pp. 1270-1281 ◽  
Author(s):  
Hiromasa Shima
Keyword(s):  
Field Test ◽  
Numerical Models ◽  
Electrical Potential ◽  
Nonlinear Least Squares ◽  
Three Dimensions ◽  
Inversion Method ◽  
Good Resolution ◽  
Inversion Algorithm ◽  
Linear Inversion ◽  
Resistivity Model

Theoretical changes in the distribution of electrical potential near subsurface resistivity anomalies have been studied using two resistivity models. The results suggest that the greatest response from such anomalies can be observed with buried electrodes, and that the resistivity model of a volume between boreholes can be accurately reconstructed by using crosshole data. The distributive properties of crosshole electrical potential data obtained by the pole‐pole array method have also been examined using the calculated partial derivative of the observed apparent resistivity with respect to a small cell within a given volume. The results show that for optimum two‐dimensional (2-D) and three‐dimensional (3-D) target imaging, in‐line data and crossline data should be combined, and an area outside the zone of exploration should be included in the analysis. In this paper, the 2-D and 3-D resistivity images presented are reconstructed from crosshole data by the combination of two inversion algorithms. The first algorithm uses the alpha center method for forward modeling and reconstructs a resistivity model by a nonlinear least‐squares inversion. Alpha centers express a continuously varying resistivity model, and the distribution of the electrical potential from the model can be calculated quickly. An initial general model is determined by the resistivity backprojection technique (RBPT) prior to the first inversion step. The second process uses finite elements and a linear inversion algorithm to improve the resolution of the resistivity model created by the first step. Simple 2-D and 3-D numerical models are discussed to illustrate the inversion method used in processing. Data from several field studies are also presented to demonstrate the capabilities of using crosshole resistivity exploration techniques. The numerical experiments show that by using the combined reconstruction algorithm, thin conductive layers can be imaged with good resolution for 2-D and 3-D cases. The integration of finite‐element computations is shown to improve the image obtained by the alpha center inversion process for 3-D applications. The first field test uses horizontal galleries to evaluate complex 2-D features of a zinc mine. The second field test illustrates the use of three boreholes at a dam site to investigate base rock features and define the distribution of an altered zone in three dimensions.

On the Calculation of the Boundary Layer

1940 ◽  
Vol 44 (359) ◽  
pp. 795-801 ◽  
Author(s):  
L. Prandtl Gottingen
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
High Reynolds Number ◽  
Main Flow ◽  
Friction Zone ◽  
Narrow Zone ◽  
Low Viscosity ◽  
Effective Components ◽  
High Reynolds

The field of viscosity forces in fluid currents of low viscosity—i.e., high Reynolds number—is frequently confined to a narrow zone in the vicinity of the surface; consequently, calculation can be simplified by the following assumptions :—1. In the friction zone, customarily called the “ boundary layer,” only the most effective components of the viscosity force μΔω, are considered—e.g., for the x component of this force, μΔu (calling the direction of main flow parallel to the surface, the x direction, and that perpendicular thereto the y direction) the expressions μ (∂2u/∂x2) and μ (∂2u/∂z2) can be neglected for μ (∂2u/∂y2).

scholarly journals Data-Driven, Multi-Model Workflow Suggests Strong Influence from Hurricanes on the Generation of Turbidity Currents in the Gulf of Mexico

2020 ◽  
Vol 8 (8) ◽  
pp. 586
Author(s):  
Courtney Harris ◽  
Jaia Syvitski ◽  
H.G. Arango ◽  
E.H. Meiburg ◽  
Sagy Cohen ◽  
...  
Keyword(s):  
Gulf Of Mexico ◽  
Continental Shelf ◽  
Numerical Models ◽  
Ocean Dynamics ◽  
Turbidity Currents ◽  
Regional Ocean Modeling System ◽  
Transport Pathway ◽  
Seabed Sediment ◽  
Continental Scale ◽  
Sediment Routing

Turbidity currents deliver sediment rapidly from the continental shelf to the slope and beyond; and can be triggered by processes such as shelf resuspension during oceanic storms; mass failure of slope deposits due to sediment- and wave-pressure loadings; and localized events that grow into sustained currents via self-amplifying ignition. Because these operate over multiple spatial and temporal scales, ranging from the eddy-scale to continental-scale; coupled numerical models that represent the full transport pathway have proved elusive though individual models have been developed to describe each of these processes. Toward a more holistic tool, a numerical workflow was developed to address pathways for sediment routing from terrestrial and coastal sources, across the continental shelf and ultimately down continental slope canyons of the northern Gulf of Mexico, where offshore infrastructure is susceptible to damage by turbidity currents. Workflow components included: (1) a calibrated simulator for fluvial discharge (Water Balance Model - Sediment; WBMsed); (2) domain grids for seabed sediment textures (dbSEABED); bathymetry, and channelization; (3) a simulator for ocean dynamics and resuspension (the Regional Ocean Modeling System; ROMS); (4) A simulator (HurriSlip) of seafloor failure and flow ignition; and (5) A Reynolds-averaged Navier–Stokes (RANS) turbidity current model (TURBINS). Model simulations explored physical oceanic conditions that might generate turbidity currents, and allowed the workflow to be tested for a year that included two hurricanes. Results showed that extreme storms were especially effective at delivering sediment from coastal source areas to the deep sea, at timescales that ranged from individual wave events (~hours), to the settling lag of fine sediment (~days).

scholarly journals How do basin margins record long-term tectonic and climatic changes?

Geology ◽  
2020 ◽  
Vol 48 (9) ◽  
pp. 893-897
Author(s):  
Jinyu Zhang ◽  
Zoltán Sylvester ◽  
Jacob Covault
Keyword(s):  
Sea Level ◽  
Numerical Models ◽  
Source Area ◽  
Climatic Changes ◽  
Steady States ◽  
Sediment Supply ◽  
Relative Sea Level ◽  
Continental Scale ◽  
Sediment Routing ◽  
Basin Margin

Abstract A long-standing goal of sedimentary geoscience is to understand how tectonic and climatic changes are reflected in basin fill. Here, we use 14 numerical models of continental-scale sediment-routing systems spanning millions of years to investigate the responses of sediment supply and basin sedimentation to changes in uplift and precipitation in the source area. We also investigate the extent to which these signals can be altered by relative sea level (the sum of subsidence and eustasy). In cases of constant relative sea level, sediment supply and margin progradation have similar responses because nearly all of the sediment is transported beyond the coastal plain and continental shelf to the basin margin. Thus, margin progradation can be used as a proxy for sediment supply. However, changes in uplift and precipitation result in different erosional patterns in the source area and different basin-margin depositional patterns. Changes in uplift result in gradual (over several million years) adjustment to new steady states of source-area erosion and margin progradation, whereas changes in precipitation result in abrupt changes in erosion and progradation followed by a return to the initial steady states. In cases of changing relative sea level, sediment storage on the shelf attenuates signals of uplift, but signals of precipitation change can be interpreted in the basin-margin record because climate-induced sediment supply changes are large enough to influence margin progradation. Understanding the relationship between sediment supply and basin-margin progradation, and their linked responses to forcings, improves our ability to interpret signals of environmental change in the stratigraphic record.

scholarly journals The vortex-entrainment sheet in an inviscid fluid: theory and separation at a sharp edge

2019 ◽  
Vol 866 ◽  
pp. 660-688 ◽  
Author(s):  
A. C. DeVoria ◽  
K. Mohseni
Keyword(s):  
Mass Density ◽  
Normal Component ◽  
Vortex Sheet ◽  
Three Dimensions ◽  
Sharp Edge ◽  
Pressure Jump ◽  
Inviscid Fluid ◽  
Fluid Theory ◽  
Viscous Boundary ◽  
Singularity Removal

In this paper a model for viscous boundary and shear layers in three dimensions is introduced and termed a vortex-entrainment sheet. The vorticity in the layer is accounted for by a conventional vortex sheet. The mass and momentum in the layer are represented by a two-dimensional surface having its own internal tangential flow. Namely, the sheet has a mass density per-unit-area making it dynamically distinct from the surrounding outer fluid and allowing the sheet to support a pressure jump. The mechanism of entrainment is represented by a discontinuity in the normal component of the velocity across the sheet. The velocity field induced by the vortex-entrainment sheet is given by a generalized Birkhoff–Rott equation with a complex sheet strength. The model was applied to the case of separation at a sharp edge. No supplementary Kutta condition in the form of a singularity removal is required as the flow remains bounded through an appropriate balance of normal momentum with the pressure jump across the sheet. A pressure jump at the edge results in the generation of new vorticity. The shedding angle is dictated by the normal impulse of the intrinsic flow inside the bound sheets as they merge to form the free sheet. When there is zero entrainment everywhere the model reduces to the conventional vortex sheet with no mass. Consequently, the pressure jump must be zero and the shedding angle must be tangential so that the sheet simply convects off the wedge face. Lastly, the vortex-entrainment sheet model is demonstrated on several example problems.

NEW FICTITIOUS DOMAIN METHODS: FORMULATION AND ANALYSIS

2005 ◽  
Vol 15 (10) ◽  
pp. 1575-1594 ◽  
Author(s):  
IVO BABUŠKA ◽  
EUGENE G. PODNOS ◽  
GREGORY J. RODIN
Keyword(s):  
Model Problem ◽  
Nonlinear Problems ◽  
Theory Of Elasticity ◽  
Two Dimensions ◽  
Three Dimensions ◽  
Fictitious Domain ◽  
Element Mesh ◽  
Multiple Interfaces ◽  
Fictitious Domain Methods ◽  
The Cost

Two closely-related fictitious domain methods for solving problems involving multiple interfaces are introduced. Like other fictitious domain methods, the proposed methods simplify the task of finite element mesh generation and provide access to solvers that can take advantage of uniform structured grids. The proposed methods do not involve the Lagrange multipliers, which makes them quite different from existing fictitious domain methods. This difference leads to an advantageous form of the inf–sup condition, and allows one to avoid time-consuming integration over curvilinear surfaces. In principle, the proposed methods have the same rate of convergence as existing fictitious domain methods. Nevertheless it is shown that, at the cost of introducing additional unknowns, one can improve the quality of the solution near the interfaces. The methods are presented using a two-dimensional model problem formulated in the context of linearized theory of elasticity. The model problem is sufficient for presenting method details and mathematical foundations. Although the model problem is formulated in two dimensions and involves only one interface, there are no apparent conceptual difficulties to extending the methods to three dimensions and multiple interfaces. Further, it is possible to extend the methods to nonlinear problems involving multiple interfaces.

scholarly journals Coupling ice flow models of varying orders of complexity with the Tiling method

2012 ◽  
Vol 58 (210) ◽  
pp. 776-786 ◽  
Author(s):  
Helene Seroussi ◽  
Hachmi Ben Dhia ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Eric Rignot ◽  
...  
Keyword(s):  
Numerical Models ◽  
Computational Cost ◽  
Higher Order ◽  
Ice Flow ◽  
Flow Models ◽  
Continental Scale ◽  
Critical Areas ◽  
Warming Climate ◽  
Computationally Intensive ◽  
A New Technique

AbstractIce flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models, instead of the traditional shallow-ice approximation, whose assumptions are not valid in certain critical areas. These higher-order models are, however, computationally intensive and difficult to use at the continental scale. Here we present a new technique, the Tiling method, to couple ice flow models of varying orders of complexity. The goal of the method is to limit the spatial extent of where higherorder models are applied to reduce the computational cost, while maintaining the model precision. We apply this method on synthetic geometries to demonstrate its practical use. We first use a geometry for which all models yield the same results to check the consistency of the method. Then we apply our method to a geometry for which a full-Stokes model is required in the vicinity of the ice front. Our results show that the hybrid models present significant improvements over mono-model approaches and reduce computational times.

Sign in / Sign up

Export Citation Format

Share Document