Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

Ambipolar diffusion in the solar atmosphere

The process of non-linear ambipolar diffusion in the region overlying the solar surface can be an effective mechanism for producing sharp magnetic structures and current sheets. These may be the sites responsible for the occurrence of connectivity of magnetic field lines, and the subsequent explosive input of energy for heating of some of the features in the atmosphere of the Sun..     

MONITORING THE DIFFUSION OF SUSPENDED SEDIMENTS AND STABILITY OF TIDAL RADIAL SAND RIDGES AREA USING MULTI-SATELLITE REMOTE SENSING DATA

Nine Landsat TM tapes and images and MSS images, 10 NOAA tapes and images and 1 SAR image from 1973 to 1997 were used to analyse the diffusion of suspended sediments and the change of tidal radial sand ridges in the northern part of the Changjiang River delta, the South Yellow Sea. The results showed that the diffusion of suspended sediments was controlled by the tide, net current, and submarine topography in this area. The distribution of suspended sediments had close relationship with thesubmarine topography. The old Huanghe River delta and the Changjiang River comprise the main sediment supply for the formation of radial sand ridges, whose evolution can be divided into three stages since the Huangbe River changed its course and flowed northward into the Bohai Sea.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Asymmetrically zoned reaction rims: assessment of grain boundary diffusivities and growth rates related to natural diffusion-controlled mineral reactions

This study explores garnet coronas around hedenbergite, which were formed by the reaction plagioclase + hedenbergite→garnet + quartz, to derive information about diffusion paths that allowed for material redistribution during reaction progress. Whereas quartz forms disconnected single grains along the garnet/hedenbergite boundaries, garnet forms ～20‐μm‐wide continuous polycrystalline rims along former plagioclase/hedenbergite phase boundaries. Individual garnet crystals are separated by low‐angle grain boundaries, which commonly form a direct link between the reaction interfaces of the plagioclase|garnet|hedenbergite succession. Compositional variations in garnet involve: (i) an overall asymmetric compositional zoning in Ca, Fe²⁺, Fe³⁺ and Al across the garnet layer; and (ii) micron‐scale compositional variations in the near‐grain boundary regions and along plagioclase/garnet phase boundaries. These compositional variations formed during garnet rim growth. Thereby, transfer of the chemical components occurred by a combination of fast‐path diffusion along grain boundaries within the garnet rim, slow diffusion through the interior of the garnet grains, and by fast diffusion along the garnet/plagioclase and the garnet/hedenbergite phase boundaries. Numerical simulation indicates that diffusion of Ca, Al and Fe²⁺ occurred about three to four, four and six to seven orders of magnitude faster along the grain boundaries than through the interior of the garnet grains. Fast‐path diffusion along grain boundaries contributed substantially to the bulk material transfer across the growing garnet rim. Despite the contribution of fast‐path diffusion, bulk diffusion through the garnet rim was too slow to allow for chemical equilibration of the phases involved in garnet rim formation even on a micrometre scale. Based on published garnet volume diffusion data the growth interval of a 20‐μm‐wide garnet rim is estimated at ～10³–10⁴ years at the inferred reaction conditions of 760 ± 50 °C at 7.6 kbar. Using the same parameterization of the growth law, 100‐μm‐ and 1‐mm‐thick garnet rims would grow within 10⁵–10⁶ and 10⁶–10⁷ years respectively.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Possible Density Change of the Origin of North Pacific Intermediate Water Due to Double Diffusion in the Mixed Water Region

In this study we test Talley's hypothesis that Oyashio winter mixed-layer water (26.5–26.6σ _θ) increases its density to produce the North Pacific Intermediate Water (NPIW) salinity minimum (26.7– 26.8σ_θ) in the Mixed Water Region, assuming a combination of cabbeling and double diffusion. The possible density change of Oyashio winter mixed-layer water is discussed using an instantaneous ratio of the change of temperature and salinity along any particular intrusion (R _l ). We estimate the range of R _l ^DD required to convert Oyashio winter mixed-layer water to the NPIW salinity minimum due to double diffusion, and then assume double-diffusive intrusions as this conversion mechanism. A double-diffusive intrusion model is used to estimate R _l ^DD in a situation where salt fingering dominates vertical mixing, as well as to determine whether Oyashio winter mixed-layer water can become the NPIW salinity minimum. Possible density changes are estimated from the model R _l ^DD by assuming the amount of density change due to cabbeling. From these results, we conclude that Oyashio winter mixed-layer water contributes to a freshening of the lighter layer of the NPIW salinity minimum (around 26.70σ_θ) in the MWR.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.