Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

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.

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.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier 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.

Incentivos e condi????es para o desempenho dos servidores p??blicos: conclus??es a partir de um estudo sobre a Austr??lia e o Brasil

Para a maior parte das sociedades contempor??neas tornou-se praticamente imposs??vel dissociar a boa governan??a da pr??tica de princ??pios democr??ticos, de um lado, e da produ????o de pol??ticas e servi??os p??blicos de qualidade, de outro. Essa vis??o sugere que a quest??o de melhorar continuamente o desempenho das m??quinas governamentais est?? agora no topo das preocupa????es desses governos. As recentes ondas de reforma da administra????o p??blica refletem isso. Do desempenho depende, em grande medida, a legitima????o desses governos (Peters& Pierre, 2010, ). Exatamente para obter melhor desempenho e maior legitimidade, os governos t??m procurado aproximar cada vez mais a burocracia dos cidad??os, para tornar os servi??os p??blicos em maior conformidade com suas necessidades e prefer??ncias, mais r??pidos e mais precisos. Entre as in??meras quest??es que emergem, no campo do desempenho das organiza????es p??blicas, est?? a da gest??o dos seus recursos humanos: quais s??o os elementos que contribuem para a motiva????o dos servidores p??blicos e quais incentivos s??o mais eficazes s??o perguntas sempre presentes, pois do sucesso dessas pol??ticas internas depende, em larga medida, a pr??pria percep????o que os cidad??os fazem de seus respectivos governos. Afinal, como argumentam Peters e Pierre (2010), os cidad??os se confrontam, com muito maior frequ??ncia, com a burocracia do que com os agentes pol??ticos. Esse estudo tem por objetivo colocar algumas hip??teses sobre quais s??o os incentivos- bem-sucedidos- ao desempenho dos servidores p??blicos, comparando as formas de gest??o dos governos federais da Austr??lia e do Brasil.

Occurrence of storm-generated bedforms along the inner continental shelf - Southeastern Brazil

Levantamentos sonográficos desenvolvidos em duas áreas da plataforma interna adjacente ao estado do Espírito Santo revelam formas de fundo cuja origem é fortemente relacionada a eventos de tempestades. Estas são caracterizadas por alternâncias de bandas de areias finas sobrepostas, de forma alternada, a um fundo de areias grossas, gerando manchas de areias grossas intercaladas abruptamente com faixas de areias finas a profundidades de 25-30 m para GUA e de 05-08 m para BES. Em ambas as áreas as faixas de areias grossas apresentam marcas de ondulação geradas por ondas com orientação levemente paralela a linha de costa. A sedimentologia de GUA é composta por areias finas a muito finas lamosas carbonáticas com cascalhos siliciclásticos e a de BES por areias grossas e médias com cascalhos biodetríticos. Desta forma, é clara a ação de correntes/ondas de tempestade no fundo marinho de GUA e de BES; porém, a classificação quanto ao comportamento hidrodinâmico e ao transporte de sedimento é necessária para as estimativas sobre a formação, transporte e manutenção das feições encontradas.

Caracterização morfodinâmica do litoral Norte Fluminense, RJ, Brasil

A morfologia das praias é reflexo da ação hidrodinfunica com o tipo de sedimento disponível. Esta interação se dá a partir da base de ação das ondas, tornando-se mais efetiva na zona de arrebentação onde a energia é dissipada, originando na zona de surfe que se estende em direção à linha de costa até Q limite do espraiamento na face de praia. O objetivo deste trabalho é apresentar as características morfodinâmicas do litoral NE do Estado do Rio de Janeiro, baseando-se na análise da variação espaço Hemporal de perfis de praia, na morfologia da plataforma continental e em fotos aéreas. Os perfis de praia foram levantados durante cinco anos pela PETROBRÁS S/A, em oito estações ao longo do litoral. O comportamento morfodinâmico distinto de setores do litoral estudado é função do padrão de ondas e direção da linha de costa, morfologia da plataforma interna controlando a extensão da zona de surfe e o número de zonas de arrebentação, características dos sedimentos e gradiente da face de praia. Foram reconhecidos quatro compartimentos morfodinâmicos distintos, com base na morfologia dos perfis praiais e plataforma interna, no índice de mobilidade da linha de praia, granulometria x gradiente da face de praia e o parâmetro Q. Os quatro compartimentos são: Atafona/Foz do Rio Paraíba do Sul - estágio morfodinâmico intermediário a dissipativo, com o maior índice de mobilidade e elevada taxa de erosão; Sul de Atafona-Cabo de São Tomé - estágio intermediário a refletivo com baixo índice de mobilidade da praia; Cabo de São Tomé - estágio refletivo a intermediário com alto índice de mobilidade da praia; Cabo de São Tomé-Cabiúnas - estágio refletivo com baixo índice de mobilidade.

Adapted Gaussian basis sets for atoms from Li through Xe generated with the generator coordinate Hartree-Fock method

Utiliza-se o método coordenada geradora Hartree-Fock para gerar bases Gaussianas adaptadas para os átomos de Li (Z=3) até Xe (Z=54). Neste método, integram-se as equações de Griffin-Hill-Wheeler-Hartree-Fock através da técnica de discretização integral. Comparam-se as funções de ondas geradas neste trabalho com as funções de ondas Roothaan-Hartree-Fock de Clementi e Roetti (1974) e com outros conjuntos de bases relatados na literatura. Para os átomos estudados aqui, os erros em nossas energias totais relativos aos limites numéricos Hartree-Fock são sempre menores que 7,426 milihartree.

Gravitational waves : a 100-year tool applied to the dark energy problem

Recent observations from type Ia Supernovae and from cosmic microwave background (CMB) anisotropies have revealed that most of the matter of the Universe interacts in a repulsive manner, composing the so-called dark energy constituent of the Universe. Determining the properties of dark energy is one of the most important tasks of modern cosmology and this is the main motivation for this work. The analysis of cosmic gravitational waves (GW) represents, besides the CMB temperature and polarization anisotropies, an additional approach in the determination of parameters that may constrain the dark energy models and their consistence. In recent work, a generalized Chaplygin gas model was considered in a flat universe and the corresponding spectrum of gravitational waves was obtained. In the present work we have added a massless gas component to that model and the new spectrum has been compared to the previous one. The Chaplygin gas is also used to simulate a L-CDM model by means of a particular combination of parameters so that the Chaplygin gas and the L-CDM models can be easily distinguished in the theoretical scenarios here established. We find that the models are strongly degenerated in the range of frequencies studied. This degeneracy is in part expected since the models must converge to each other when some particular combinations of parameters are considered.

Metodologia para obtenção do hidrograma de escoamento superficial em encostas e canais. Parte I: desenvolvimento e avaliação

Desenvolveu-se uma metodologia que permite obter o hidrograma de escoamento superficial e a vazão máxima para qualquer posição ao longo de uma encosta e para seções transversais de canais utilizando o modelo de ondas cinemáticas. A área da encosta é dividida num sistema matricial composto por 100 linhas e 100 colunas. Na encosta, considera-se que o escoamento ocorre na direção da declividade e que a vazão de cada pixel é a soma da vazão produzida nesse com a vazão advinda dos pixels que contribuem com o escoamento superficial para o pixel em análise. No canal, a vazão é calculada pela soma dos hidrogramas advindos das colunas do sistema reticulado. A comparação entre os valores de lâmina e vazão máxima de escoamento superficial obtidas experimentalmente e calculadas em duas condições (encosta e bacia) permitiu evidenciar que a metodologia, comparada aos métodos Racional e do Número da Curva, ofereceu boas estimativas tanto da lâmina quanto da vazão máxima de escoamento superficial.