Temperature effect in potassium and nitrate ions in soil transport

When doing researches on solute dynamics in porous medium, the knowledge of medium characteristics and percolating liquids, as well as of external factors is very important. An important external factor is temperature and, in this sense, our purpose was determining potassium and nitrate transport parameters for different values of temperature, in miscible displacement experiments. Evaluated parameters were retardation factor (R), diffusion/dispersion coefficient (D) and dispersivity, at ambient temperature (25 up to 28 ºC), 40 ºC and 50 ºC. Salts used were potassium nitrate and potassium chlorate, prepared in a solution made up of 5 ppm nitrate and 2.000 ppm potassium, with Red-Yellow Latosol porous medium. Temperature exhibited a positive influence upon porous medium solution and upon dispersion coefficient.

Liquid transportation and distribution during (re)wetting : impact on coated papers in printing

In many industrial applications, such as the printing and coatings industry, wetting of porous materials by liquids includes not only imbibition and permeation into the bulk but also surface spreading and evaporation. By understanding these phenomena, valuable information can be obtained for improved process control, runnability and printability, in which liquid penetration and subsequent drying play important quality and economic roles. Knowledge of the position of the wetting front and the distribution/degree of pore filling within the structure is crucial in describing the transport phenomena involved. Although exemplifying paper as a porous medium in this work, the generalisation to dynamic liquid transfer onto a surface, including permeation and imbibition into porous media, is of importance to many industrial and naturally occurring environmental processes. This thesis explains the phenomena in the field of heatset web offset printing but the content and the analyses are applicable in many other printing methods and also other technologies where water/moisture monitoring is crucial in order to have a stable process and achieve high quality end products. The use of near-infrared technology to study the water and moisture response of porous pigmented structures is presented. The use of sensitive surface chemical and structural analysis, as well as the internal structure investigation of a porous structure, to inspect liquid wetting and distribution, complements the information obtained by spectroscopic techniques. Strong emphasis has been put on the scale of measurement, to filter irrelevant information and to understand the relationship between interactions involved. The near-infrared spectroscopic technique, presented here, samples directly the changes in signal absorbance and its variation in the process at multiple locations in a print production line. The in-line non-contact measurements are facilitated by using several diffuse reflectance probes, giving the absolute water/moisture content from a defined position in the dynamic process in real-time. The nearinfrared measurement data illustrate the changes in moisture content as the paper is passing through the printing nips and dryer, respectively, and the analysis of the mechanisms involved highlight the roles of the contacting surfaces and the relative liquid carrier properties of both non-image and printed image areas. The thesis includes laboratory studies on wetting of porous media in the form of coated paper and compressed pigment tablets by mono-, dual-, and multi-component liquids, and paper water/moisture content analysis in both offline and online conditions, thus also enabling direct sampling of temporal water/moisture profiles from multiple locations. One main focus in this thesis was to establish a measurement system which is able to monitor rapid changes in moisture content of paper. The study suggests that near-infrared diffuse reflectance spectroscopy can be used as a moisture sensitive system and to provide accurate online qualitative indicators, but, also, when accurately calibrated, can provide quantification of water/moisture levels, its distribution and dynamic liquid transfer. Due to the high sensitivity, samples can be measured with excellent reproducibility and good signal to noise ratio. Another focus of this thesis was on the evolution of the moisture content, i.e. changes in moisture content referred to (re)wetting, and liquid distribution during printing of coated paper. The study confirmed different wetting phases together with the factors affecting each phase both for a single droplet and a liquid film applied on a porous substrate. For a single droplet, initial capillary driven imbibition is followed by equilibrium pore filling and liquid retreat by evaporation. In the case of a liquid film applied on paper, the controlling factors defining the transportation were concluded to be the applied liquid volume in relation to surface roughness, capillarity and permeability of the coating giving the liquid uptake capacity. The printing trials confirmed moisture gradients in the printed sheet depending on process parameters such as speed, fountain solution dosage and drying conditions as well as the printed layout itself. Uneven moisture distribution in the printed sheet was identified to be one of the sources for waving appearance and the magnitude of waving was influenced by the drying conditions.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

Industrial radioactive barite scale: suppression of radium uptake by introduction of competing ions

Incorporation of radioactive isotopes during the formation of barite mineral scale is a widespread phenomenon occurring within the oil, mining and process industries. In a series of experiments radioactive barite/celestite solid solutions (SSBarite-Celcstite) have been synthesized under controlled conditions by the counter diffusion of Ra-226, Ba2+, Sr24+ and SO42- ions through a porous medium (silica gel), to investigate inhibiting effects in Ra uptake associated with the introduction of a competing ion (Sr2+). From characterization studies, the particle size and the morphology of the crystals appear to be related to the initial [Sr]/[Ba] molar ratio of the starting solution. Typically, systems richer in Sr produce smaller sized crystals and clusters characterized by a lower degree of order. The activity introduced to the system is mainly incorporated in the crystals generated from the barite/celestite solid solution as suggested by the activity profiles of the hydrogel columns analysed by gamma-spectrometry. There is a relationship between the initial [Sr]/[Ba] molar ratio of the starting solution and the activity exhibited by the synthesized crystals. An effective inhibition of the Ra-226 uptake during formation of the crystals (SSBarite-Celestite) was obtained through the introduction of a competing ion (Sr2+): the higher the initial [Sr]/[Ba] molar ratio of the starting solution, the lower the intensity of the activity peak in the crystals. (C) 2003 Published by Elsevier Ltd.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Velocity-based moving mesh methods for nonlinear partial differential equations

This article describes a number of velocity-based moving mesh numerical methods formultidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

Direct Computation of Stochastic Flow in Reservoirs with Uncertain Parameters

A direct method is presented for determining the uncertainty in reservoir pressure, flow, and net present value (NPV) using the time-dependent, one phase, two- or three-dimensional equations of flow through a porous medium. The uncertainty in the solution is modelled as a probability distribution function and is computed from given statistical data for input parameters such as permeability. The method generates an expansion for the mean of the pressure about a deterministic solution to the system equations using a perturbation to the mean of the input parameters. Hierarchical equations that define approximations to the mean solution at each point and to the field covariance of the pressure are developed and solved numerically. The procedure is then used to find the statistics of the flow and the risked value of the field, defined by the NPV, for a given development scenario. This method involves only one (albeit complicated) solution of the equations and contrasts with the more usual Monte-Carlo approach where many such solutions are required. The procedure is applied easily to other physical systems modelled by linear or nonlinear partial differential equations with uncertain data.

Sea ice is a mushy layer

[1] Sea ice is a two-phase, two-component, reactive porous medium: an example of what is known in other contexts as a mushy layer. The fundamental conservation laws underlying the mathematical description of mushy layers provide a robust foundation for the prediction of sea-ice evolution. Here we show that the general equations describing mushy layers reduce to the model of Maykut and Untersteiner (1971) under the same approximations employed therein.

Modelagem e previsão da perda de injetividade em poços canhoneados

Waterflooding is a technique largely applied in the oil industry. The injected water displaces oil to the producer wells and avoid reservoir pressure decline. However, suspended particles in the injected water may cause plugging of pore throats causing formation damage (permeability reduction) and injectivity decline during waterflooding. When injectivity decline occurs it is necessary to increase the injection pressure in order to maintain water flow injection. Therefore, a reliable prediction of injectivity decline is essential in waterflooding projects. In this dissertation, a simulator based on the traditional porous medium filtration model (including deep bed filtration and external filter cake formation) was developed and applied to predict injectivity decline in perforated wells (this prediction was made from history data). Experimental modeling and injectivity decline in open-hole wells is also discussed. The injectivity of modeling showed good agreement with field data, which can be used to support plan stimulation injection wells

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.