Theoretical and numerical analyses of convective instability in porous media with upward throughflow

Exact analytical solutions have been obtained for a hydrothermal system consisting of a horizontal porous layer with upward throughflow. The boundary conditions considered are constant temperature, constant pressure at the top, and constant vertical temperature gradient, constant Darcy velocity at the bottom of the layer. After deriving the exact analytical solutions, we examine the stability of the solutions using linear stability theory and the Galerkin method. It has been found that the exact solutions for such a hydrothermal system become unstable when the Rayleigh number of the system is equal to or greater than the corresponding critical Rayleigh number. For small and moderate Peclet numbers (Pe less than or equal to 6), an increase in upward throughflow destabilizes the convective flow in the horizontal layer. To confirm these findings, the finite element method with the progressive asymptotic approach procedure is used to compute the convective cells in such a hydrothermal system. Copyright (C) 1999 John Wiley & Sons, Ltd.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

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.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Kinetics of photoacclimation in corals

Traditional models describing the relationship between photosynthesis (P) and irradiance (I) do not account for photoacclimation to short-term variation in irradiance. Here we develop and test a model that predicts the rate of photosynthesis under fluctuating irradiances at the scale of days to weeks. Using oxygen respirometry, we measured the rates of change in the P-I model parameters P-max (maximum rate of gross photosynthesis) and I-k (sub-saturation irradiance) of the photo-symbiotic coral Turbinaria mesenterina (Lamarck) following large and small increases and decreases in growth irradiance. We analyse the behaviour of the dynamic P-I model in turbid-water conditions using a dataset of 3-month continuous irradiance as the input variable. In response to upward or downward changes in experimental growth irradiance, I-k values decreased or increased exponentially, reaching new and stable levels within 5-10 days. I-k responded 4 times stronger than P-max to changes in growth irradiance. The kinetics of I-k did not show hysteresis, and changed in similar ways when irradiance was increased or decreased in small or large amounts. This suggests that mechanisms associated with photo-protection during increases in irradiance, and the maximisation of photosynthetic efficiency during decreases in irradiance, are equally potent. On the scale of months, the dynamic P-I model did not predict higher rates of photosynthesis than the static P-I model, but buffered the variation in photosynthesis during periods of reduced irradiance. Fourier analysis indicated that the kinetics of I-k closely matches the main periodicities in daily irradiance (1-2 weeks). The recorded kinetics of photoacclimation in the Turbinaria-zooxanthella symbiosis is comparable to that of free-living phytoplankton and faster than that of higher plants.

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.

Age changes in personality traits and their heritabilities during the adult years: evidence from Australian Twin Registry samples

Short versions of four Eysenck personality scales had been included in questionnaires given to several adult samples from the Australian Twin Registry, comprising altogether some 5400 pairs. Means and regressions with age are compared for three samples at average ages of 23, 37, and 61 years, and for two samples of retested individuals, one tested twice at average ages of 29 and 37 years, and one tested three times at average ages of 45, 56, and 62 years, For both males and females the trends for Psychoticism (P), Extraversion (E), and Neuroticism (N) were generally downward with age, and for Lie (L), upward. However, in the longitudinal sample between ages 56 and 62 the trends for P, E, and I stopped or reversed, although N continued downward. Heritabilities were reasonably stable across age for P, E, and N, and the effects of shared environments negligible, but L showed some influence of shared environment as well as genes in all but the oldest age group. (C) 2001 Elsevier Science Ltd. All rights reserved.

Facies architecture of the CRP-3 drillhole, Victoria Land Basin, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diametites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ("Motif A"). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as "Motif B. Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ("o", "p" and "q") correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

The entrapment of particles by windbreaks

A theory is developed for calculating the entrapment of particles by a windbreak, with four results. (1) The fraction of particles in the oncoming flow which pass through the windbreak, or transmittance of the windbreak for particles (sigma), is related to the optical porosity (tau). The very simple approximation sigma=tau works well for most applications involving the interception of spray droplets by windbreaks. Results from a field experiment agree with the theoretical predictions. (2) A new equation for the bulk drag coefficient of a windbreak is tested against numerical, wind tunnel and field experiments. This enables the bleed velocity for the flow through the windbreak to be predicted in terms of the screen pressure coefficient (k) of the barrier. (3) The relationship between k and tau is different for a vegetative barrier than for a screen across a confined duct, implying a lower Fc for given tau. (4) The total deposition of particles to a windbreak is determined by a trade-off between particle absorption and throughflow, implying an optimum value of tau for maximum total deposition. For particles larger than 30 mum and vegetation elements smaller than 30 mm, this occurs near tau = 0.2. (C) 2001 Elsevier Science Ltd. All rights reserved.

Cytological basis of photoresponsive behavior in a sponge larva

Ontogenetic changes in the photoresponse of larvae from the demosponge Reneira sp. were studied by analyzing the swimming paths of individual larvae exposed to diffuse white light. Larvae swam upward upon release from the adult, but were negatively phototactic until at least 12 hours after release. The larval photoreceptors are presumed to be a posterior ring of columnar monociliated epithelial cells that possess 120-mum-long cilia and pigment-filled protrusions. A sudden increase in light intensity caused these cilia to become rigidly straight. If the light intensity remained high, the cilia gradually bent over the pigmented vesicles in the adjacent cytoplasm, and thus covered one entire pole of the larva. The response was reversed upon a sudden decrease in light intensity. The ciliated cells were sensitive to changes in light intensity in larvae of all ages. This response is similar to the shadow response in tunicate larvae or the shading of the photoreceptor in Euglena and is postulated to allow the larvae to steer away from brighter light to darker areas, such as under coral rubble-the preferred site of the adult sponge on the reef flat. In the absence of a coordinating system in cellular sponges, the spatial organization and autonomous behavior of the pigmented posterior cells control the rapid responses to light shown by these larvae.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.