scholarly journals Interaction between Crustal-Scale Darcy and Hydrofracture Fluid Transport: A Numerical Study

Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Tamara de Riese ◽  
Paul D. Bons ◽  
Enrique Gomez-Rivas ◽  
Till Sachau
Keyword(s):  
Fluid Flow ◽  
Numerical Model ◽  
Transport Mechanism ◽  
Large Scale ◽  
Fluid Transport ◽  
Numerical Study ◽  
Ballistic Transport ◽  
Accretionary Prism ◽  
Small Scale ◽  
Porous Flow

Crustal-scale fluid flow can be regarded as a bimodal transport mechanism. At low hydraulic head gradients, fluid flow through rock porosity is slow and can be described as diffusional. Structures such as hydraulic breccias and hydrothermal veins both form when fluid velocities and pressures are high, which can be achieved by localized fluid transport in space and time, via hydrofractures. Hydrofracture propagation and simultaneous fluid flow can be regarded as a “ballistic” transport mechanism, which is activated when transport by diffusion alone is insufficient to release the local fluid overpressure. The activation of a ballistic system locally reduces the driving force, through allowing the escape of fluid. We use a numerical model to investigate the properties of the two transport modes in general and the transition between them in particular. We developed a numerical model in order to study patterns that result from bimodal transport. When hydrofractures are activated due to low permeability relative to fluid flux, many hydrofractures form that do not extend through the whole system. These abundant hydrofractures follow a power-law size distribution. A Hurst factor of ~0.9 indicates that the system self-organizes. The abundant small-scale hydrofractures organize the formation of large-scale hydrofractures that ascend through the whole system and drain fluids in large bursts. As the relative contribution of porous flow increases, escaping fluid bursts become less frequent, but more regular in time and larger in volume. We propose that metamorphic rocks with abundant veins, such as in the Kodiak accretionary prism (Alaska) and Otago schists (New Zealand), represent regions with abundant hydrofractures near the fluid source, while hydrothermal breccias are formed by the large fluid bursts that can ascend the crust to shallower levels.

scholarly journals Numerical Study of Violent Impact Flow Using a CIP-Based Model

2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Qiao-ling Ji ◽  
Xi-zeng Zhao ◽  
Sheng Dong
Keyword(s):  
Free Surface ◽  
Numerical Model ◽  
Large Scale ◽  
Stokes Equations ◽  
Numerical Study ◽  
Three Dimensional ◽  
Small Scale ◽  
Surface Boundary ◽  
Two Phase ◽  
Constrained Interpolation

A two-phase flow model is developed to study violent impact flow problem. The model governed by the Navier-Stokes equations with free surface boundary conditions is solved by a Constrained Interpolation Profile (CIP)-based high-order finite difference method on a fixed Cartesian grid system. The free surface is immersed in the computation domain and expressed by a one-fluid density function. An accurate Volume of Fluid (VOF)-type scheme, the Tangent of Hyperbola for Interface Capturing (THINC), is combined for the free surface treatment. Results of another two free surface capturing methods, the original VOF and CIP, are also presented for comparison. The validity and utility of the numerical model are demonstrated by applying it to two dam-break problems: a small-scale two-dimensional (2D) and three-dimensional (3D) full scale simulations and a large-scale 2D simulation. Main attention is paid to the water elevations and impact pressure, and the numerical results show relatively good agreement with available experimental measurements. It is shown that the present numerical model can give a satisfactory prediction for violent impact flow.

A numerical study of the transition of jet diffusion flames

1990 ◽  
Vol 431 (1882) ◽  
pp. 301-314 ◽  
Keyword(s):  
Reynolds Number ◽  
Diffusion Coefficient ◽  
Large Scale ◽  
Numerical Study ◽  
Small Scale ◽  
Surface Model ◽  
Flame Surface ◽  
Transition Mechanism ◽  
Air Stream ◽  
Scale Fluctuation

A numerical study on the transition from laminar to turbulent of two-dimensional fuel jet flames developed in a co-flowing air stream was made by adopting the flame surface model of infinite chemical reaction rate and unit Lewis number. The time dependent compressible Navier–Stokes equation was solved numerically with the equation for coupling function by using a finite difference method. The temperature-dependence of viscosity and diffusion coefficient were taken into account so as to study effects of increases of these coefficients on the transition. The numerical calculation was done for the case when methane is injected into a co-flowing air stream with variable injection Reynolds number up to 2500. When the Reynolds number was smaller than 1000 the flame, as well as the flow, remained laminar in the calculated domain. As the Reynolds number was increased above this value, a transition point appeared along the flame, downstream of which the flame and flow began to fluctuate. Two kinds of fluctuations were observed, a small scale fluctuation near the jet axis and a large scale fluctuation outside the flame surface, both of the same origin, due to the Kelvin–Helmholtz instability. The radial distributions of density and transport coefficients were found to play dominant roles in this instability, and hence in the transition mechanism. The decreased density in the flame accelerated the instability, while the increase in viscosity had a stabilizing effect. However, the most important effect was the increase in diffusion coefficient. The increase shifted the flame surface, where the large density decrease occurs, outside the shear layer of the jet and produced a thick viscous layer surrounding the jet which effectively suppressed the instability.

scholarly journals Stochastic modelling and diffusion modes for proper orthogonal decomposition models and small-scale flow analysis

2017 ◽  
Vol 826 ◽  
pp. 888-917 ◽  
Author(s):  
Valentin Resseguier ◽  
Etienne Mémin ◽  
Dominique Heitz ◽  
Bertrand Chapron
Keyword(s):  
Fluid Flow ◽  
Velocity Component ◽  
Large Scale ◽  
Stochastic Modelling ◽  
Wake Flow ◽  
Small Scale ◽  
Cylinder Wake ◽  
Diffusion Term ◽  
Advection Term ◽  
Dimensional Models

We present here a new stochastic modelling approach in the constitution of fluid flow reduced-order models. This framework introduces a spatially inhomogeneous random field to represent the unresolved small-scale velocity component. Such a decomposition of the velocity in terms of a smooth large-scale velocity component and a rough, highly oscillating component gives rise, without any supplementary assumption, to a large-scale flow dynamics that includes a modified advection term together with an inhomogeneous diffusion term. Both of those terms, related respectively to turbophoresis and mixing effects, depend on the variance of the unresolved small-scale velocity component. They bring an explicit subgrid term to the reduced system which enables us to take into account the action of the truncated modes. Besides, a decomposition of the variance tensor in terms of diffusion modes provides a meaningful statistical representation of the stationary or non-stationary structuration of the small-scale velocity and of its action on the resolved modes. This supplies a useful tool for turbulent fluid flow data analysis. We apply this methodology to circular cylinder wake flow at Reynolds numbers $Re=100$ and $Re=3900$. The finite-dimensional models of the wake flows reveal the energy and the anisotropy distributions of the small-scale diffusion modes. These distributions identify critical regions where corrective advection effects, as well as structured energy dissipation effects, take place. In providing rigorously derived subgrid terms, the proposed approach yields accurate and robust temporal reconstruction of the low-dimensional models.

scholarly journals The new design of the wood stove based on the numerical analysis and experimental research

2020 ◽  
Vol 154 ◽  
pp. 04001
Author(s):  
Przemysław Motyl ◽  
Marcin Wikło ◽  
Julita Bukalska ◽  
Bartosz Piechnik ◽  
Rafał Kalbarczyk
Keyword(s):  
Fluid Flow ◽  
Heat Exchange ◽  
Numerical Analysis ◽  
Heat Release ◽  
Experimental Research ◽  
Numerical Study ◽  
Combustion Process ◽  
Small Scale ◽  
Wood Stove ◽  
Wood Stoves

In Europe, especially in Poland, wood-fired stoves remain one of the most popular renewable household heating. The use of wood logs in small-scale units stoves are expected to increase substantially. The work proposes a comprehensive approach to modify the design of wood stoves with heating power up to 20 kW, including design works, simulations, and experimental research. The article also presents the numerical study of a combustion process including fluid flow, chemical combustion reaction, and heat exchange in the wood stove. The retrofit enhanced a more stable heat release from the wood stove, which increased efficiency and reduction of the harmful components of combustion.

Multiscale 3D stress field modelling for the URL 'Reiche Zeche' using a discontinuum model approach

2020 ◽  
Author(s):  
Sebastian Rehde ◽  
Prof. Dr.-Ing. habil. Heinz Konietzky
Keyword(s):  
Numerical Model ◽  
Stress Field ◽  
Numerical Modelling ◽  
Large Scale ◽  
Small Scale ◽  
Hydraulic Stimulation ◽  
Field Modelling ◽  
Slip Tendency ◽  
Project Site ◽  
Model Approach

<p>Underneath the small town of Freiberg, Saxony, stretches the ore mine complex 'Reiche Zeche'. The underground laboratory (URL) inside the mine was inaugurated in 1919 and is an internationally acknowledged institution for experimental work of variable scales and subjects. Our work is part of the Stimtec project, which aims on improving planning and conducting hydraulic stimulation in anisotropic, crystalline rocks. The project comprises numerical modelling and field work inside the URL. Prior to the numerical analysis, we implemented a tool to perform a slip tendency analysis of faults that were mapped along the tunnel walls at the project site. It allows to assess the slip tendency of arbitrarily oriented faults and stress fields. The tool is used for preselection of stimulation intervals, enabling identification of faults which are likely to be reactivated by hydraulic stimulation. <br>We perform the stress field modelling using a multiscale numerical model approach. Therefore, we set up three different sized models deriving from a large scale 3D geomodel. The geomodel contains the topography, drifts and 47 fault structures taken from mine maps. The project site and measurement points are positioned in the center of the model. From the large scale geomodel, we developed a simplified numerical model geometry with 12 major faults, disregarding the galleries. We use the distinct element code 3DEC for discontinuous numerical modelling of the stress field. This allows to take into account discrete displacements along the faults. Far field stress is taken from previous investigations and literature as boundary and initial conditions. The resulting stress  field provides the stress tensors for calculating the corresponding forces for each gridpoint at the model boundaries of the small scale model. The small scale numerical model is smaller by a factor of 10, including two major fault segments, the galleries and mapped local faults. Hydraulic fracturing stress measurements taken during the field tests indicate that the stress field is strongly distorted in the vicinity of the tunnels and excavations along the ore veins. Hence, we developed a third model approach, a 2.5D slice model, to investigate the influence of the assumed excavation damage zones.<br>With this work, we provide an approach to predict the stress field inside the complex, anisotropic rock volume. Within the framework of the Stimtec project, we developed a workflow for planning hydraulic stimulation tests and 3D geological models for a diverse set of further appliations in the URL 'Reiche Zeche'.</p>

scholarly journals Impact Assessment of Orography on the Extreme Precipitation Event of July 2010 over Pakistan: A Numerical Study

2015 ◽  
Vol 2015 ◽  
pp. 1-20 ◽  
Author(s):  
Khan Muhammad Tahir ◽  
Yan Yin ◽  
Yong Wang ◽  
Zaheer A. Babar ◽  
Dong Yan
Keyword(s):  
Large Scale ◽  
Extreme Event ◽  
Numerical Study ◽  
Deep Convection ◽  
Temporal Distribution ◽  
Precipitation Event ◽  
Windward Side ◽  
Small Scale ◽  
Extreme Rainfall Event ◽  
Synoptic Conditions

The topography influences monsoon precipitation and gives rise to significant rainfall events in South Asia. The physical mechanism involved in such events includes mechanical uplifting, thermodynamics, small scale cloud processes, and large scale atmospheric circulations. The investigation into orographic precipitation is pursued by synoptic and model analysis. Deep convection occurs as warm moist airflow is channeling over steep mountains. WRF model coupled with Morrison double moment scheme is used to assess the relative impact of topography on extreme rainfall event of 26–30 July 2010 in Pakistan. Two sensitivity tests with full topography (CTL) and reduced topography by 50% (LOW) are carried out. Two distinct precipitation zones over Hindukush and Himalaya mountains are identified. The topographic changes significantly affect moisture divergence and spatial and temporal distribution of precipitation. A low level jet is created on windward side of big mountains, yielding enhanced moisture flux and instability. Eddy kinetic energy significantly changes with orographic height. Energy flux created further unstabilized atmosphere and deep convection, producing wide spread heavy rainfall in the area in Himalaya foothills. Under the set synoptic conditions, orographic orientation enhanced the moisture accumulation and deep convection, resulting in occurrence of this extreme event.

Liquid Fuel Spray Characterization for the Design of the Dual-Fuel PGT10B Combustor

2001 ◽  
Author(s):  
F. A. Tap ◽  
R. Modi ◽  
J. P. Van Buijtenen
Keyword(s):  
Liquid Film ◽  
Gas Turbine ◽  
Large Scale ◽  
Drop Size ◽  
Numerical Study ◽  
Fuel Spray ◽  
Small Scale ◽  
Dual Fuel ◽  
Test Results ◽  
Liquid Concentration

The Dry Low NOx (DLN) silo combustor of the Nuovo Pignone PGT10B gas turbine is being redesigned to meet Dual-Fuel capability. A prototype with specially designed fuel injectors, placed on airfoil-shaped elements, was tested at cold conditions (using water instead of Diesel fuel) to map the spray mass distribution at the premixer exit. The resulting profile showed high concentrations of liquid near the premixer centerline and on the premixer wall. Parallel to this test, a small-scale experimental and numerical study was made of a single atomizer of the fuel system, placed in cross flow position. This small-scale study was launched in order to gain insight in the behavior of the spray, as well as to assess the relative importance of spray modeling parameters. The PDPA experiments and 2D CFD simulations of these experiments showed fair agreement on the average drop size distribution and drop size-velocity correlation. The flow visualization also revealed liquid film formation on the surface of the airfoil, behind the injector, due to the low atomization pressure differential at cold conditions. Using this modeling experience, the spray patternation test with the prototype combustor has been modeled using an existing 3D CFD model of the premixer. The model also showed high liquid concentration on the wall, but not near the centerline. From the results of the small-scale study it is concluded that the measured high concentration near the premixer centerline is not a result of the flow field. It is assumed that in the complete assembly the liquid film from the injector vanes accumulates on the center body, resulting in a high liquid concentration downstream on the premixer centerline. Overall, the application of CFD analyses on the tests performed proved to be a very useful tool to evaluate the test results. The modeling experience identified the important factors in modeling the fuel spray in a gas turbine environment, but further evolution of computer resources is required before large-scale test results will be reproducible with CFD models.

2D and 3D Multiscale Computational Modeling of Dynamic Microorgan Devices as Drug Screening Platforms

2015 ◽  
Author(s):  
Filippos Tourlomousis ◽  
Robert C. Chang
Keyword(s):  
Large Scale ◽  
Numerical Study ◽  
Numerical Results ◽  
Scale Model ◽  
Small Scale ◽  
Optimum Process ◽  
Shear Stresses ◽  
Enzymatic Reactions ◽  
Modeling Approach

The ability to incorporate three-dimensional (3D) hepatocyte-laden hydrogel constructs using layered fabrication approaches into devices that can be perfused with drugs enables the creation of dynamic microorgan devices (DMDs) that offer an optimal analog of the in vivo liver metabolism scenario. The dynamic nature of such in vitro metabolism models demands reliable numerical tools to determine the optimum process, material, and geometric parameters for the most effective metabolic conversion of the perfused drug into the liver microenvironment. However, there is a current lack of literature that integrates computational approaches to guide the optimum design of such devices. The groundwork of the present numerical study has been laid by our previous study [1], where the authors modeled in 2D an in vitro DMD of arbitrary dimensions and identified the modeling challenges towards meaningful results. These constructs are hosted in the chamber of the microfluidic device serving as walls of the microfluidic array of channels through which a fluorescent drug substrate is perfused into the microfluidic printed channel walls at a specified volumetric flow rate assuring Stokes flow conditions (Re<<1). Due to the porous nature of the hydrogel walls, a metabolized drug product is collected at the outlet port. A rigorous FEM based modeling approach is presented for a single channel parallel model geometry (1 free flow channel with 2 porous walls), where the hydrodynamics, mass transfer and pharmacokinetics equations are solved numerically in order to yield the drug metabolite concentration profile at the DMD outlet. The fluid induces shear stresses are assessed both in 3D, with only 27 cells modeled as single compartment voids, where all of the enzymatic reactions are assumed to take place. In this way, the mechanotransduction effect that alters the hepatocyte metabolic activity is assessed for a small scale model. This approach overcomes the numerical limitations imposed by the cell density (∼1012 cells/m3) of the large scale DMD device. In addition, a compartmentalization technique is proposed in order to assess the metabolism process at the subcellular level. The numerical results are validated with experiments to reveal the robustness of the proposed modeling approach and the necessity of scaling the numerical results by preserving dynamic and biochemical similarity between the small and large scale model.

Sign in / Sign up

Export Citation Format

Share Document