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.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Observations of anti-winds in a deep alpine valley, Lake Tekapo, New Zealand

This paper presents observations of summertime anti-winds monitored under ideal conditions in the Lake Tekapo hydro-catchment situated in the central Southern Alps, New Zealand. Onset and cessation of anti-winds was observed to coincide with the change in phase of the surface limbs of thermally generated valley and mountain winds under settled anti-cyclonic conditions. Anti-winds were best developed in the early morning before surface heating and associated convective mixing of the valley atmosphere began to mask the boundaries between the surface based limb of the mountain-valley wind and the corresponding anti-wind. By mid-day, the anti-valley wind exceeded the height of the surrounding ridgeline and became embedded in the topographically channeled gradient wind. Observations presented here have both theoretical and applied implications with regard to the development of thermally generated wind systems in deep alpine valleys, and their role in the dispersion of air pollution.

Asymptotic behaviour of two-dimensional time-delayed Navier-Stokes equations

We consider the two-dimensional Navier-Stokes equations with a time-delayed convective term and a forcing term which contains some hereditary features. Some results on existence and uniqueness of solutions are established. We discuss the asymptotic behaviour of solutions and we also show the exponential stability of stationary solutions.

Particulate urban air pollution affects the functional morphology of mouse placenta

in humans, adverse pregnancy outcomes (low birth weight, prematurity, and intrauterine growth retardation) are associated with exposure to urban air pollution. Experimental data have also shown that such exposure elicits adverse reproductive outcomes. We hypothesized that the effects of urban air pollution on pregnancy outcomes could be related to changes in functional morphology of the placenta. To test this, future dams were exposed during pregestational and gestational periods to filtered or nonfiltered air in exposure chambers. Placentas were collected from near-term pregnancies and prepared for microscopical examination. Fields of view on vertical uniform random tissue slices were analyzed using stereological methods. Volumes of placental compartments were estimated, and the labyrinth was analyzed further in terms of its maternal vascular spaces, fetal capillaries, trophoblast, and exchange surface areas. From these primary data, secondary quantities were derived: vessel calibers (expressed as diameters), trophoblast thickness (arithmetic mean), and total and mass-specific morphometric diffusive conductances for oxygen of the intervascular barrier. Two-way analysis of variance showed that both periods of exposure led to significantly smaller fetal weights. Pregestational exposure to nonfiltered air led to significant increases in fetal capillary surface area and in total and mass-specific conductances. However, the calibers of maternal blood spaces were reduced. Gestational exposure to nonfiltered air was associated with reduced volumes, calibers, and surface areas of maternal blood spaces and with greater fetal capillary surfaces and diffusive conductances. The findings indicate that urban air pollution affects placental functional morphology. Fetal weights are compromised despite attempts to improve diffusive transport across the placenta.

The application of electron microscopy and microanalysis in conjunction with organic petrology to further the understanding of organic-mineral association: examples from Mount Isa and McArthur Basins, Australia

Organic petrology supported by electron microscopical and micro-analytical techniques was applied to organic matter in Proterozoic sediments to better understand hydrothermal processes responsible for ore-grade mineralisation. It was shown that organic maturation was not only closely linked to the geological history of the sediments, but also highlighted heat transfer by convection as differentiated from conduction solely through sediment burial and step-wise subsidence. Water-rock ratios effect organic maturation in hydrothermal systems, and erratic reflectance profiles are indicators of convective heat transfer. Identification and characterisation of organic materials in terms of source rock and migrated hydrocarbons was shown to be a powerful tool in reconstructing the thermal history of sediments, identifying hydrothermal episodes, fluid pathways and heat source in the northern Australian Proterozoic basins. Higher reflectance of organic matter towards the central parts of the Mount Isa Basin and some of the most northerly parts point to proximity to higher heat flow at times, in contrast to relatively low temperatures prevailing in the western parts of the basin, next to the Murphy Inlier. A close correlation shown between peak organic reflectance values and super-sequence boundaries farther highlighted the valuable information to be gained from organic petrology, by allowing the separation of processes responsible for metal dissolution and transportation from those inducing precipitation. (C) 2001 Elsevier Science B.V All rights reserved.

Continuous quantum measurement of two coupled quantum dots using a point contact: A quantum trajectory approach

We obtain the finite-temperature unconditional master equation of the density matrix for two coupled quantum dots (CQD's) when one dot is subjected to a measurement of its electron occupation number using a point contact (PC). To determine how the CQD system state depends on the actual current through the PC device, we use the so-called quantum trajectory method to derive the zero-temperature conditional master equation. We first treat the electron tunneling through the PC barrier as a classical stochastic point process (a quantum-jump model). Then we show explicitly that our results can be extended to the quantum-diffusive limit when the average electron tunneling rate is very large compared to the extra change of the tunneling rate due to the presence of the electron in the dot closer to the PC. We find that in both quantum-jump and quantum-diffusive cases, the conditional dynamics of the CQD system can be described by the stochastic Schrodinger equations for its conditioned state vector if and only if the information carried away from the CQD system by the PC reservoirs can be recovered by the perfect detection of the measurements.

Dynamics of a mesoscopic charge quantum bit under continuous quantum measurement

We present the conditional quantum dynamics of an electron tunneling between two quantum dots subject to a measurement using a low transparency point contact or tunnel junction. The double dot system forms a single qubit and the measurement corresponds to a continuous in time readout of the occupancy of the quantum dot. We illustrate the difference between conditional and unconditional dynamics of the qubit. The conditional dynamics is discussed in two regimes depending on the rate of tunneling through the point contact: quantum jumps, in which individual electron tunneling current events can be distinguished, and a diffusive dynamics in which individual events are ignored, and the time-averaged current is considered as a continuous diffusive variable. We include the effect of inefficient measurement and the influence of the relative phase between the two tunneling amplitudes of the double dot/point contact system.

A comparison of thermal condition between pilot- and full-scale furnaces for studying slagging and fouling propensity in PF boilers

The Australian Coal Industry Research Laboratory (ACIRL) furnace is scaled to simulate slagging and fouling in operating boilers. This requires that the gas and target temperatures, the heat flux, and the flow pattern be the same as those in real boilers. The gas and target temperatures are maintained by insulating the wall and cooling the target respectively. The flow pattern of a small burner cannot be the same as a large furnace. However, this flow pattern is partially compensated for by placing the slagging panels in three vertical locations. The paper develops the models of radiant heat transfer from the flame to the deposits both in pilot-scale and full-scale furnaces. They are used to compare the effective radiant heat transfer of the pilot- and full-scale furnaces. The experimental data both from the pilot- and full-scale furnaces are used to verify the incident heat flux and temperature profiles in the pilot- and full-scale furnaces. The results showed that the thermal condition in the pilot-scale furnace meets the requirements for studying the slagging regarding the gas temperature and the incident heat flux, particularly for the panel #1. The gas temperature in the convective section also meets the requirement for studying the fouling.

Rebreathing potential of infant mattresses and bedcovers

Objective : To establish the CO2 dispersion and retention properties of some mattresses and bed coverings commercially available in Australia. Methods : Five mattresses were studied in (i) an in vivo model in which an infant's head was covered by a headbox, rebreathing was allowed to occur, and the final steady state CO2 concentration was measured; and (ii) an in vitro model in which 5% CO2 in a headbox was allowed to disperse, and the time taken for the concentration to reach 1% was measured. Five types of bedcover were studied in (i) an in vivo model in which an infant's head was covered by a bedcover and the final steady state CO2 concentration was measured; and (ii) an in vitro model in which 5% CO2 under a bedcover was allowed to disperse, and the time taken for the concentration to reach 1% was measured. Results : The steady state CO2 concentrations ranged from 0.6% to 3.0% for the mattresses (P < 0.05). The time for CO2 to disperse ranged from 5.5 min to 30.4 min (P < 0.05). Steady state CO2 concentrations ranged from 2.5% to 3.6% for the bedcoverings (P > 0.05). The time for CO2 to disperse ranged from 5.4 min to 7.7 min (P > 0.05). Conclusions : Some commercial cot mattresses and bedcoverings allow high concentrations of CO2 to accumulate in rebreathing environments. Some mattress types studied were more diffusive to CO2 , whereas there was no difference between the bedcovers studied. This may have implications for vulnerable infants at risk of sudden infant death syndrome.

Techniques for detailed heat transfer measurements in cold supersonic blowdown tunnels using thermochromic liquid crystals

The paper presents methods for measurement of convective heat transfer distributions in a cold flow, supersonic blowdown wind tunnel. The techniques involve use of the difference between model surface temperature and adiabatic wall temperature as the driving temperature difference for heat transfer and no active heating or cooling of the test gas or model is required. Thermochromic liquid crystals are used for surface temperature indication and results presented from experiments in a Mach 3 flow indicate that measurements of the surface heat transfer distribution under swept shock wave boundary layer interactions can be made. (C) 2002 Elsevier Science Ltd. All rights reserved.

Instabilities and drop formation in cylindrical liquid jets in reduced gravity

The effects of convective and absolute instabilities on the formation of drops formed from cylindrical liquid jets of glycerol/water issuing into still air were investigated. Medium-duration reduced gravity tests were conducted aboard NASA's KC-135 and compared to similar tests performed under normal gravity conditions to aid in understanding the drop formation process. In reduced gravity, the Rayleigh-Chandrasekhar Equation was found to accurately predict the transition between a region of absolute and convective instability as defined by a critical Weber number. Observations of the physics of the jet, its breakup, and subsequent drop dynamics under both gravity conditions and the effects of the two instabilities on these processes are presented. All the normal gravity liquid jets investigated, in regions of convective or absolute instability, were subject to significant stretching effects, which affected the subsequent drop and associated geometry and dynamics. These effects were not displayed in reduced gravity and, therefore, the liquid jets would form drops which took longer to form (reduction in drop frequency), larger in size, and more spherical (surface tension effects). Most observed changes, in regions of either absolute or convective instabilities, were due to a reduction in the buoyancy force and an increased importance of the surface tension force acting on the liquid contained in the jet or formed drop. Reduced gravity environments allow better investigations to be performed into the physics of liquid jets, subsequently formed drops, and the effects of instabilities on these systems. In reduced gravity, drops form up to three times more slowly and as a consequence are up to three times larger in volume in the theoretical absolute instability region than in the theoretical convective instability region. This difference was not seen in the corresponding normal gravity tests due to the masking effects of gravity. A drop is shown to be able to form and detach in a region of absolute instability, and spanning the critical Weber number (from a region of convective to absolute instability) resulted in a marked change in dynamics and geometry of the liquid jet and detaching drops. (C) 2002 American Institute of Physics.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

The effect of heat transfer mode on heart rate responses and hysteresis during heating and cooling in the estuarine crocodile Crocodylus porosus

The effect of heating and cooling on heart rate in the estuarine crocodile Crocodylus porosus was studied in response to different heat transfer mechanisms and heat loads. Three heating treatments were investigated. C. porosus were: (1) exposed to a radiant heat source under dry conditions; (2) heated via radiant energy while half-submerged in flowing water at 23degreesC and (3) heated via convective transfer by increasing water temperature from 23degreesC to 35degreesC. Cooling was achieved in all treatments by removing the heat source and with C. porosus half-submerged in flowing water at 23degreesC. In all treatments, the heart rate of C. porosus increased markedly in response to heating and decreased rapidly with the removal of the heat source. Heart rate during heating was significantly faster than during cooling at any given body temperature, i.e. there was a significant heart rate hysteresis. There were two identifiable responses to heating and cooling. During the initial stages of applying or removing the heat source, there was a dramatic increase or decrease in heart rate ('rapid response'), respectively, indicating a possible cardiac reflex. This rapid change in heart rate with only a small change or no change in body temperature (

Estudo da camada limite atmosférica em regiões metropolitanas costeiras com simulações de brisa marítima

O principal objetivo deste trabalho foi identificar e caracterizar a evolução diária da Camada Limite Atmosférica (CLA) na Região da Grande Vitória (RGV), Estado do Espírito Santo, Brasil e na Região de Dunkerque (RD), Departamento Nord Pas-de-Calais, França, avaliando a acurácia de parametrizações usadas no modelo meteorológico Weather Research and Forecasting (WRF) em detectar a formação e atributos da Camada Limite Interna (CLI) que é formada pelas brisas marítimas. A RGV tem relevo complexo, em uma região costeira de topografia acidentada e uma cadeia de montanhas paralela à costa. A RD tem relevo simples, em uma região costeira com pequenas ondulações que não chegam a ultrapassar 150 metros, ao longo do domínio de estudos. Para avaliar os resultados dos prognósticos feitos pelo modelo, foram utilizados os resultados de duas campanhas: uma realizada na cidade de Dunkerque, no norte da França, em Julho de 2009, utilizando um sistema light detection and ranging (LIDAR), um sonic detection and ranging (SODAR) e dados de uma estação meteorológica de superfície (EMS); outra realizada na cidade de Vitória – Espírito Santo, no mês de julho de 2012, também usando um LIDAR, um SODAR e dados de uma EMS. Foram realizadas simulações usando três esquemas de parametrizações para a CLA, dois de fechamento não local, Yonsei University (YSU) e Asymmetric Convective Model 2 (ACM2) e um de fechamento local, Mellor Yamada Janjic (MYJ) e dois esquemas de camada superficial do solo (CLS), Rapid Update Cycle (RUC) e Noah. Tanto para a RGV quanto para a RD, foram feitas simulações com as seis possíveis combinações das três parametrizações de CLA e as duas de CLS, para os períodos em que foram feitas as campanhas, usando quatro domínios aninhados, sendo os três maiores quadrados com dimensões laterais de 1863 km, 891 km e 297 km, grades de 27 km, 9 km e 3 km, respectivamente, e o domínio de estudo, com dimensões de 81 km na direção Norte-Sul e 63 km na Leste-Oeste, grade de 1 km, com 55 níveis verticais, até um máximo de, aproximadamente, 13.400 m, mais concentrados próximos ao solo. Os resultados deste trabalho mostraram que: a) dependendo da configuração adotada, o esforço computacional pode aumentar demasiadamente, sem que ocorra um grande aumento na acurácia dos resultados; b) para a RD, a simulação usando o conjunto de parametrizações MYJ para a CLA com a parametrização Noah produziu a melhor estimativa captando os fenômenos da CLI. As simulações usando as parametrizações ACM2 e YSU inferiram a entrada da brisa com atraso de até três horas; c) para a RGV, a simulação que usou as parametrizações YSU para a CLA em conjunto com a parametrização Noah para CLS foi a que conseguiu fazer melhores inferências sobre a CLI. Esses resultados sugerem a necessidade de avaliações prévias do esforço computacional necessário para determinadas configurações, e sobre a acurácia de conjuntos de parametrizações específicos para cada região pesquisada. As diferenças estão associadas com a capacidade das diferentes parametrizações em captar as informações superficiais provenientes das informações globais, essenciais para determinar a intensidade de mistura turbulenta vertical e temperatura superficial do solo, sugerindo que uma melhor representação do uso de solo é fundamental para melhorar as estimativas sobre a CLI e demais parâmetros usados por modelos de dispersão de poluentes atmosféricos.