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.

Residence time distributions of gas flowing through rotating drum bioreactors

Residence time distribution studies of gas through a rotating drum bioreactor for solid-state fermentation were performed using carbon monoxide as a tracer gas. The exit concentration as a function of time differed considerably from profiles expected for plug flow, plug flow with axial dispersion, and continuous stirred tank reactor (CSTR) models. The data were then fitted by least-squares analysis to mathematical models describing a central plug flow region surrounded by either one dead region (a three-parameter model) or two dead regions (a five-parameter model). Model parameters were the dispersion coefficient in the central plug flow region, the volumes of the dead regions, and the exchange rates between the different regions. The superficial velocity of the gas through the reactor has a large effect on parameter values. Increased superficial velocity tends to decrease dead region volumes, interregion transfer rates, and axial dispersion. The significant deviation from CSTR, plug flow, and plug flow with axial dispersion of the residence time distribution of gas within small-scale reactors can lead to underestimation of the calculation of mass and heat transfer coefficients and hence has implications for reactor design and scaleup. (C) 2001 John Wiley & Sons, Inc.

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.

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.

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.

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.

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 (

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%).

Uma solução da equação difusão-advecção com o termo contragradiente

Neste trabalho, apresenta­se uma solução para equação de difusão­advecção considerando o  termo  contragradiente  que  é um termo  adicional. Esse termo  adicional  contém informações  sobre a assimetria, escala de tempo Lagrangeana e velocidade turbulenta vertical. A solução  da equação foi obtida pela utilização da técnica de Transformada de Laplace, considerando a  Camada  Limite  Planetária  (CLP)  como  um  sistema  de  multicamadas.  Os  parâmetros  turbulentos foram derivados da teoria de difusão estatística de Taylor, combinada com a teoria  da similaridade. Assim, são apresentadas simulações para diferentes valores de assimetria, o  que  propiciou  a  obtenção  de  uma  concentração  de  contaminantes  em  diferentes  alturas,  em  uma  camada  limite  convectiva.  A  avaliação  do  desempenho  do  modelo,  considerando  a  assimetria  no  processo  de  dispersão  de  poluentes  atmosféricos, foi realizada  através  de  um  experimento  de  tanque  convectivo  tradicional.  Nesse  experimento,  o  termo  contragradiente  influenciou a concentração de poluentes para uma camada limite convectiva. Entretanto, com  as  parametrizações  utilizadas,  o  modelo  não  conseguiu  captar  de  forma  eficiente  o  comportamento da concentração em pontos mais distantes da fonte.

Comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde: estabilidade de tensão

Neste trabalho é efectuado, não só o diagnóstico em regime permanente, mas também o estudo, simulação e análise do comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde. Os estudos de estabilidade transitória desempenham um importante papel, tanto no planeamento como na operação dos sistemas de potência. Tais estudos são realizados, em grande parte, através de simulação digital no domínio do tempo, utilizando integração numérica para resolver as equações não-lineares que modelam a dinâmica do sistema e dependem da existência de registos reais de perturbação (ex: osciloperturbografia). O objectivo do trabalho será também verificar a aplicabilidade dos requisitos técnicos que as unidades geradoras devem ter, no que concerne ao controlo de tensão, estabelecidos na futura regulamentação europeia desenvolvida pela ENTSO-E (European Network Transmission System Operator for Electricity). De entre os requisitos analisou-se a capacidade das máquinas existentes suportarem cavas de tensão decorrentes de curto-circuitos trifásicos simétricos, Fault Ride Through, no ponto de ligação à rede. Identificaram-se para o efeito os factores que influenciam a estabilidade desta rede, em regime perturbado nomeadamente: (i) duração do defeito, (ii) caracterização da carga, com e sem a presença do sistema de controlo de tensão (AVR) em unidades de geração síncronas. Na ausência de registos reais sobre o comportamento do sistema, conclui-se que este é sensível à elasticidade das cargas em particular do tipo potência constante, existindo risco de perda de estabilidade, neste caso, para defeitos superiores a 5ms sem AVR. A existência de AVR nesta rede afigura-se como indispensável para garantir estabilidade de tensão sendo contudo necessário proceder a uma correcta parametrização.

Distribuição de vídeo sobre uma rede IEEE 802.11n

A recente norma IEEE 802.11n oferece um elevado débito em redes locais sem fios sendo por isso esperado uma adopção massiva desta tecnologia substituindo progressivamente as redes 802.11b/g. Devido à sua elevada capacidade esta recente geração de redes sem fios 802.11n permite um crescimento acentuado de serviços audiovisuais. Neste contexto esta dissertação procura estudar a rede 802.11n, caracterizando o desempenho e a qualidade associada a um serviço de transmissão de vídeo, recorrendo para o efeito a uma arquitectura de simulação da rede 802.11n. Desta forma é caracterizado o impacto das novas funcionalidades da camada MAC introduzidas na norma 801.11n, como é o caso da agregação A-MSDU e A-MPDU, bem como o impacto das novas funcionalidades da camada física como é o caso do MIMO; em ambos os casos uma optimização da parametrização é realizada. Também se verifica que as principais técnicas de codificação de vídeo H.264/AVC para optimizar o processo de distribuição de vídeo, permitem optimizar o desempenho global do sistema de transmissão. Aliando a optimização e parametrização da camada MAC, da camada física, e do processo de codificação, é possível propor um conjunto de configurações que permitem obter o melhor desempenho na qualidade de serviço da transmissão de conteúdos de vídeo numa rede 802.11n. A arquitectura de simulação construída nesta dissertação é especificamente adaptada para suportar as técnicas de agregação da camada MAC, bem como para suportar o encapsulamento em protocolos de rede que permitem a transmissão dos pacotes de vídeo RTP, codificados em H.264/AVC.