Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones when they are heated uniformly from below. The fault zone is assumed to be more permeable than its surrounding rocks. In particular, we have derived exact analytical solutions to the total critical Rayleigh numbers of the double diffusion-driven convective flow. Using the corresponding total critical Rayleigh numbers, the double diffusion-driven convective instability of a fluid-saturated three-dimensional geological fault zone system has been investigated. The related theoretical analysis demonstrates that: (1) The relative higher concentration of the chemical species at the top of the three-dimensional geological fault zone system can destabilize the convective flow of the system, while the relative lower concentration of the chemical species at the top of the three-dimensional geological fault zone system can stabilize the convective flow of the system. (2) The double diffusion-driven convective flow modes of the three-dimensional geological fault zone system are very close each other and therefore, the system may have the similar chance to pick up different double diffusion-driven convective flow modes, especially in the case of the fault thickness to height ratio approaching 0. (3) The significant influence of the chemical species diffusion on the convective instability of the three-dimensional geological fault zone system implies that the seawater intrusion into the surface of the Earth is a potential mechanism to trigger the convective flow in the shallow three-dimensional geological fault zone system.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Self-organization with traveling waves: A case for a convective torus

A traveling wave of BaSO4 in the chlorite-thiourea reaction has shown concentric precipitation patterns upon being triggered by the autocatalyst HOCl. The precipitation patterns show circular rings of alternate null and full precipitation regions. This self-organization appears to be the result of the formation of a convective torus. The formation of the convective torus can be described as a Benard-Marangoni instability with lateral heating.

Surface stickiness of drops of carbohydrate and organic acid solutions during convective drying: Experiments and modeling

Drying kinetics of low molecular weight sugars such as fructose, glucose, sucrose and organic acid such as citric acid and high molecular weight carbohydrate such as maltodextrin (DE 6) were determined experimentally using single drop drying experiments as well as predicted numerically by solving the mass and heat transfer equations. The predicted moisture and temperature histories agreed with the experimental ones within 6% average relative (absolute) error and average difference of +/- 1degreesC, respectively. The stickiness histories of these drops were determined experimentally and predicted numerically based on the glass transition temperature (T-g) of surface layer. The model predicted the experimental observations with good accuracy. A nonsticky regime for these materials during spray drying is proposed by simulating a drop, initially 120 mum in diameter, in a spray drying environment.

Selective entrainment of sediment graded by size and density under waves

The influence of near-bed sorting processes on heavy mineral content in suspension is discussed. Sediment concentrations above a rippled bed of mixed quartz and heavy mineral sand were measured under regular nonbreaking waves in the laboratory. Using the traditional gradient diffusion process, settling velocity would be expected to strongly affect sediment distribution. This was not observed during present trials. In fact, the vertical gradients of time-averaged suspension concentrations were found to be similar for the light and heavy minerals, despite their different settling velocities. This behavior implies a convective rather than diffusive distribution mechanism. Between the nonmoving bed and the lowest suspension sampling point, fight and heavy mineral concentration differs by two orders of magnitude. This discrimination against the heavy minerals in the pickup process is due largely to selective entrainment at the ripple face. Bed-form dynamics and the nature of quartz suspension profiles are found to be little affected by the trialed proportion of overall heavy minerals in the bed (3.8-22.1%).

Convective patterns under the Indo-Atlantic << box >>

Using fluid mechanics, we reinterpret the mantle images obtained from global and regional tomography together with geochemical, geological and paleomagnetic observations, and attempt to unravel the pattern of convection in the Indo-Atlantic "box" and its temporal evolution over the last 260 Myr. The << box >> presently contains a) a broad slow seismic anomaly at the CMB which has a shape similar to Pangea 250 Myr ago, and which divides into several branches higher in the lower mantle, b) a "superswell, centered on the western edge of South Africa, c) at least 6 "primary hotspots" with long tracks related to traps, and d) numerous smaller hotspots. In the last 260 Myr, this mantle box has undergone 10 trap events, 7 of them related to continental breakup. Several of these past events are spatially correlated with present-day seismic anomalies and/or upwellings. Laboratory experiments show that superswells, long-lived hotspot tracks and traps may represent three evolutionary stages of the same phenomenon, i.e. episodic destabilization of a hot, chemically heterogeneous thermal boundary layer, close to the bottom of the mantle. When scaled to the Earth's mantle, its recurrence time is on the order of 100-200 Myr. At any given time, the Indo-Atlantic box should contain 3 to 9 of these instabilities at different stages of their development, in agreement with observations. The return flow of the downwelling slabs, although confined to two main << boxes >> (Indo-Atlantic and Pacific) by subduction zone geometry, may therefore not be passive, but rather take the form of active thermochemical instabilities. (c) 2005 Elsevier B.V. All rights reserved.

Evidences for direct magnetic patterning via diffusive transformations using femtosecond laser interferometry

The application of femtosecond laser interferometry to direct patterning of thin-film magnetic alloys is demonstrated. The formation of stripe gratings with submicron periodicities is achieved in Fe1-xVx (x=18-34wt. %) layers, with a difference in magnetic moments up to Delta mu/mu similar to 20 between adjacent stripes but without any significant development of the topographical relief (<1% of the film thickness). The produced gratings exhibit a robust effect of their anisotropy shape on magnetization curves in the film plane. The obtained data witness ultrafast diffusive transformations associated with the process of spinodal decomposition and demonstrate an opportunity for producing magnetic nanostructures with engineered properties upon this basis.

Microphysical evidence of the transition between predominant convective/stratiform rainfall associated with the intraseasonal oscillation in the Southwest Amazon

The distinction between convective and stratiform precipitation profiles around various precipitating systems existent in tropical regions is very important to the global atmospheric circulation, which is extremely sensitive to vertical latent heat distribution. In South America, the convective activity responds to the Intraseasonal Oscillation (IOS). This paper analyzes a disdrometer and a radar profiler data, installed in the Ji-Paraná airport, RO, Brazil, for the field experiment WETAMC/LBA & TRMM/LBA, during January and February of 1999. The microphysical analysis of wind regimes associated with IOS showed a large difference in type, size and microphysical processes of hydrometeor growth in each wind regime: easterly regimes had more turbulence and consequently convective precipitation formation, and westerly regimes had a more stratiform precipitation formation.

Convective dynamics of traveling autocatalytic fronts in a modulated gravity field

Aus: Physical chemistry, chemical physics, Vol. 16.2014, H. 47, S. 26279 - 26287

Contractive metrics for a Boltzmann equation for granular gases: diffusive equilibriaThe Patterson-Sullivan embedding and minimal volume entropy for outer space

We quantify the long-time behavior of a system of (partially) inelastic particles in a stochastic thermostat by means of the contractivity of a suitable metric in the set of probability measures. Existence, uniqueness, boundedness of moments and regularity of a steady state are derived from this basic property. The solutions of the kinetic model are proved to converge exponentially as t→ ∞ to this diffusive equilibrium in this distance metrizing the weak convergence of measures. Then, we prove a uniform bound in time on Sobolev norms of the solution, provided the initial data has a finite norm in the corresponding Sobolev space. These results are then combined, using interpolation inequalities, to obtain exponential convergence to the diffusive equilibrium in the strong L¹-norm, as well as various Sobolev norms.

Runaway Greenhouse Effect in a Semigray Radiative-Convective Model

The effects of the nongray absorption (i.e., atmospheric opacity varying with wavelength) on the possible upper bound of the outgoing longwave radiation (OLR) emitted by a planetary atmosphere have been examined. This analysis is based on the semigray approach, which appears to be a reasonable compromise between the complexity of nongray models and the simplicity of the gray assumption (i.e., atmospheric absorption independent of wavelength). Atmospheric gases in semigray atmospheres make use of constant absorption coefficients in finite-width spectral bands. Here, such a semigray absorption is introduced in a one-dimensional (1D) radiative– convective model with a stratosphere in radiative equilibrium and a troposphere fully saturated with water vapor, which is the semigray gas. A single atmospheric window in the infrared spectrum has been assumed. In contrast to the single absolute limit of OLR found in gray atmospheres, semigray ones may also show a relative limit. This means that both finite and infinite runaway effects may arise in some semigray cases. Of particular importance is the finding of an entirely new branch of stable steady states that does not appear in gray atmospheres. This new multiple equilibrium is a consequence of the nongray absorption only. It is suspected that this new set of stable solutions has not been previously revealed in analyses of radiative–convective models since it does not appear for an atmosphere with nongray parameters similar to those for the earth’s current state