Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

Interações entre a carboximetilcelulose, carbonato de cálcio e bentonita: repercussões sobre as propriedades dos fluidos de perfuração aquosos

The role of carboxymethylcellulose (CMC) in association to calcium carbonate particles (CaCO3) in most water-based drilling fluids is to reduce the fluid loss to the surrounding formation. Another essential function is to provide rheological properties capable of maintaining in suspension the cuttings during drilling operation. Therefore, it is absolutely essential to correlate the polymer chemical structure (degree of substitution, molecular weight and distribution of substituent) with the physical-chemical properties of CaCO3, in order to obtain the better result at lower cost. Another important aspect refers to the clay hydration inhibitive properties of carboxymethylcellulose (CMC) in drilling fluids systems. The clay swelling promotes an undesirable damage that reduces the formation permeability and causes serious problems during the drilling operation. In this context, this thesis consists of two main parts. The first part refers to understanding of interactions CMC-CaCO3, as well as the corresponding effects on the fluid properties. The second part is related to understanding of mechanisms by which CMC adsorption occurs onto the clay surface, where, certainly, polymer chemical structure, ionic strength, molecular weight and its solvency in the medium are responsible to affect intrinsically the clay layers stabilization. Three samples of carboximetilcellulose with different molecular weight and degree of substitution (CMC A (9 x 104 gmol DS 0.7), CMC B (2.5 x 105 gmol DS 0.7) e CMC C (2.5 x 105 gmol DS 1.2)) and three samples of calcite with different average particle diameter and particle size distribution were used. The increase of CMC degree of substitution contributed to increase of polymer charge density and therefore, reduced its stability in brine, promoting the aggregation with the increase of filtrate volume. On the other hand, the increase of molecular weight promoted an increase of rheological properties with reduction of filtrate volume. Both effects are directly associated to hydrodynamic volume of polymer molecule in the medium. The granulometry of CaCO3 particles influenced not only the rheological properties, due to adsorption of polymers, but also the filtration properties. It was observed that the lower filtrate volume was obtained by using a CaCO3 sample of a low average size particle with wide dispersion in size. With regards to inhibition of clay swelling, the CMC performance was compared to other products often used (sodium chloride (NaCl), potassium chloride (KCl) and quaternary amine-based commercial inhibitor). The low molecular weight CMC (9 x 104 g/mol) showed slightly lower swelling degree compared to the high molecular weight (2.5 x 105 g/mol) along to 180 minutes. In parallel, it can be visualized by Scanning Electron Microscopy (SEM) that the high molecular weight CMC (2.5 x 105 g/mol e DS 0.7) promoted a reduction in pores formation and size of clay compared to low molecular weight CMC (9.0 x 104 g/mol e DS 0.7), after 1000 minutes in aqueous medium. This behavior was attributed to dynamic of interactions between clay and the hydrodynamic volume of CMC along the time, which is result of strong contribution of electrostatic interactions and hydrogen bounds between carboxylate groups and hydroxyls located along the polymer backbone and ionic and polar groups of clay surface. CMC adsorbs on clay surface promoting the skin formation , which is responsible to minimize the migration of water to porous medium. With the increase of degree of substitution, it was observed an increase of pores onto clay, suggesting that the higher charge density on polymer is responsible to decrease its flexibility and adsorption onto clay surface. The joint evaluation of these results indicate that high molecular weight is responsible to better results on control of rheological, filtration and clay swelling properties, however, the contrary effect is observed with the increase of degree of substitution. On its turn, the calcite presents better results of rheological and filtration properties with the decrease of average viii particle diameter and increase of particle size distribution. According to all properties evaluated, it has been obvious the interaction of CMC with the minerals (CaCO3 and clay) in the aqueous medium

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Fluid dynamics of air in a packed bed: velocity profiles and the continuum model assumption

Air flow through packed beds was analyzed experimentally under conditions ranging from those that reinforce the effect of the wall on the void fraction to those that minimize it. The packing was spherical particles, with a tube-to-particle diameter ratio (D/dp) between 3 and 60. Air flow rates were maintained between 1.3 and 4.44 m3/min, and gas velocity was measured with a Pitot tube positioned above the bed exit. Measurements were made at various radial and angular coordinate values, allowing the distribution of air flow across the bed to be described in detail. Comparison of the experimentally observed radial profiles with those derived from published equations revealed that at high D/dp ratios the measured and calculated velocity profiles behaved similarly. At low ratios, oscillations in the velocity profiles agreed with those in the voidage profiles, signifying that treating the porous medium as a continuum medium is questionable in these cases.

Simulação numérica de trocadores de calor multitubulares equipados com fitas helicoidais externamente aos tubos

Classical shell-and-tube heat exchangers are usually equipped with segmental baffles. These baffles serve two basic functions: (a) they provide tube supports, thereby preventing or reducing mechanical problems, such as sagging or vibration; (b) they direct the fluid flow over the tubes so as to introduce a cross-flow component, thereby increasing the heat transfer. Segmented baffles have several sources of performance loss, some due to various leakage flows and others caused by stagnation zones. A new concept of longitudinal flow heat exchanger - based on placing twisted tapes along the tube bundle subchannels - was developed to mitigate drawbacks of other types of tubular heat exchangers. In this paper, a numerical model has been implemented in order to simulate the thermal-hydraulic feature of tubular heat exchangers equipped either with segmental baffles or with subchannel twisted tapes. The tube bundle has been described by means of an equivalent porous medium type model, allowing a macroscopic description of the shell-side flow. The basic equations - continuity, momentum and energy - have been solved by using the finite volume method. Typical numerical results have been compared with experimental data, reaching a very good agreement. A comparative analysis of different types of heat exchangers has been carried out, revealing the satisfactory thermal-hydraulic efficiency level of the twisted tapes heat exchangers.

Estimativa de alguns parâmetros de transporte de contaminantes para solos Brasileiros

This paper introduces an indirect estimate for the coefficients of distribution, hydrodynamic dispersion and retardation for contaminants commonly encountered in sanitary landfills and their liners, such as Cu2+ and K+; this estimate is based on the relationship between concentration and certain physical characteristics of typical Brazilian soils. The results of previous studies investigating the migration of contaminants were used to develop mathematical expressions from multiple non-linear regressions. Using minimal squares regression, this transport was linked to various combinations of contaminant concentration and both structural and textural characteristics of the porous medium. Various combinations of characteristics and concentrations were investigated, with a mathematical expression obtained for each. The relationship between percentage of clay and the contaminant content proved to be the most closely correlated with actual transport parameters, with coefficients close to one.

Aplicação da técnica de elementos finitos na estratificação de armazenadores térmicos com barreiras de meios porosos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Escoamentos em meios porosos: uma análise teórica e experimental dos efeitos dos esforços capilares provocados por percolação de água em elementos de alvenaria

Pós-graduação em Engenharia Mecânica - FEG

Medidores inteligentes: o primeiro passo em direção as redes inteligentes

Smart grids are the focus of major study today because of the necessity of modernization in electrical systems and reduction of greenhouse gas emissions that increases global warming. Reaching the best deployment method, you must first of all know the current electrical system and how to use them for the benefit of this new technology. Preparing the action plan we should be aware of the main points of smart grids in each step of the electricity system - generation, transmission and distribution. Analyzed these topics, this work will focus on the first step in the implementation of the smart grids: the smart meters, tool which is already being implemented in Brazil. The main characteristics and applications of these devices, as well as their communication structure with the core distributors will be showed during the paper. Finally, we present a case study which will be discussed and analyzed based in the results obtained with the implementation of smart meters in the city of Vancouver, Canada, where we have a considerable savings already in the first year, with fully paying the initial investment and still have a profit

Simulação numérica de fluxo de águas subterrâneas na fazenda Cara Preta - Santa Rita do Passa Quatro - SP

The term model refers to any representation of a real system. The use of models in Hydrogeology can be valuable predictive tools for management of groundwater resources. The numeric models of groundwater flow, object of this study, consist on a set of differential equations that describe the water flow in the porous medium. In this context, numeric simulations were made for a sub-basin located at Cara Preta farm – Santa Rita do Passa Quatro – SP. The aquifer at the local is composed by rocks of Pirambóia Formation, which is part of Guarani Aquifer System. It was developed a conceptual model from previous studies in the area, and from that, simulations were made through the software Visual Modflow®. The conceptual model established previously was considered consistent through the results of simulation.

Condutividade hidráulica saturada e direcional em xisto: ensaios em laboratório. -

Pós-graduação em Geociências e Meio Ambiente - IGCE

Comportamento das curvas de saturação da água e LNAPL em solo da formação Rio Claro

Pós-graduação em Geociências e Meio Ambiente - IGCE

Modeling hydrothermal system: deriving observables and hydrothermal instability in volcanic and non-volcanic setting

Hydrothermal fluids are a fundamental resource for understanding and monitoring volcanic and non-volcanic systems. This thesis is focused on the study of hydrothermal system through numerical modeling with the geothermal simulator TOUGH2. Several simulations are presented, and geophysical and geochemical observables, arising from fluids circulation, are analyzed in detail throughout the thesis. In a volcanic setting, fluids feeding fumaroles and hot spring may play a key role in the hazard evaluation. The evolution of the fluids circulation is caused by a strong interaction between magmatic and hydrothermal systems. A simultaneous analysis of different geophysical and geochemical observables is a sound approach for interpreting monitored data and to infer a consistent conceptual model. Analyzed observables are ground displacement, gravity changes, electrical conductivity, amount, composition and temperature of the emitted gases at surface, and extent of degassing area. Results highlight the different temporal response of the considered observables, as well as the different radial pattern of variation. However, magnitude, temporal response and radial pattern of these signals depend not only on the evolution of fluid circulation, but a main role is played by the considered rock properties. Numerical simulations highlight differences that arise from the assumption of different permeabilities, for both homogeneous and heterogeneous systems. Rock properties affect hydrothermal fluid circulation, controlling both the range of variation and the temporal evolution of the observable signals. Low temperature fumaroles and low discharge rate may be affected by atmospheric conditions. Detailed parametric simulations were performed, aimed to understand the effects of system properties, such as permeability and gas reservoir overpressure, on diffuse degassing when air temperature and barometric pressure changes are applied to the ground surface. Hydrothermal circulation, however, is not only a characteristic of volcanic system. Hot fluids may be involved in several mankind problems, such as studies on geothermal engineering, nuclear waste propagation in porous medium, and Geological Carbon Sequestration (GCS). The current concept for large-scale GCS is the direct injection of supercritical carbon dioxide into deep geological formations which typically contain brine. Upward displacement of such brine from deep reservoirs driven by pressure increases resulting from carbon dioxide injection may occur through abandoned wells, permeable faults or permeable channels. Brine intrusion into aquifers may degrade groundwater resources. Numerical results show that pressure rise drives dense water up to the conduits, and does not necessarily result in continuous flow. Rather, overpressure leads to new hydrostatic equilibrium if fluids are initially density stratified. If warm and salty fluid does not cool passing through the conduit, an oscillatory solution is then possible. Parameter studies delineate steady-state (static) and oscillatory solutions.

Critical review of coupled flux formulations for clay membranes based on nonequilibrium thermodynamics

Extensive research conducted over the past several decades has indicated that semipermeable membrane behavior (i.e., the ability of a porous medium to restrict the passage of solutes) may have a significant influence on solute migration through a wide variety of clay-rich soils, including both natural clay formations (aquitards, aquicludes) and engineered clay barriers (e.g., landfill liners and vertical cutoff walls). Restricted solute migration through clay membranes generally has been described using coupled flux formulations based on nonequilibrium (irreversible) thermodynamics. However, these formulations have differed depending on the assumptions inherent in the theoretical development, resulting in some confusion regarding the applicability of the formulations. Accordingly, a critical review of coupled flux formulations for liquid, current, and solutes through a semipermeable clay membrane under isothermal conditions is undertaken with the goals of explicitly resolving differences among the formulations and illustrating the significance of the differences from theoretical and practical perspectives. Formulations based on single-solute systems (i.e., uncharged solute), single-salt systems, and general systems containing multiple cations or anions are presented. Also, expressions relating the phenomenological coefficients in the coupled flux equations to relevant soil properties (e.g., hydraulic conductivity and effective diffusion coefficient) are summarized for each system. A major difference in the formulations is shown to exist depending on whether counter diffusion or salt diffusion is assumed. This difference between counter and salt diffusion is shown to affect the interpretation of values for the effective diffusion coefficient in a clay membrane based on previously published experimental data. Solute transport theories based on both counter and salt diffusion then are used to re-evaluate previously published column test data for the same clay membrane. The results indicate that, despite the theoretical inconsistency between the counter-diffusion assumption and the salt-diffusion conditions of the experiments, the predictive ability of solute transport theory based on the assumption of counter diffusion is not significantly different from that based on the assumption of salt diffusion, provided that the input parameters used in each theory are derived under the same assumption inherent in the theory. Nonetheless, salt-diffusion theory is fundamentally correct and, therefore, is more appropriate for problems involving salt diffusion in clay membranes. Finally, the fact that solute diffusion cannot occur in an ideal or perfect membrane is not explicitly captured in any of the theoretical expressions for total solute flux in clay membranes, but rather is generally accounted for via inclusion of an effective porosity, n(e), or a restrictive tortuosity factor, tau(r), in the formulation of Fick's first law for diffusion. Both n(e) and tau(r) have been correlated as a linear function of membrane efficiency. This linear correlation is supported theoretically by pore-scale modeling of solid-liquid interactions, but experimental support is limited. Additional data are needed to bolster the validity of the linear correlation for clay membranes. Copyright 2012 Elsevier B.V. All rights reserved.