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.

Caractérisation des occupations du sol en milieu urbain par imagerie radar

Cette étude vise à tester la pertinence des images RSO - de moyenne et de haute résolution - à la caractérisation des types d’occupation du sol en milieu urbain. Elle s’est basée sur des approches texturales à partir des statistiques de deuxième ordre. Plus spécifiquement, on recherche les paramètres de texture les plus pertinents pour discriminer les objets urbains. Il a été utilisé à cet égard des images Radarsat-1 en mode fin en polarisation HH et Radarsat-2 en mode fin en double et quadruple polarisation et en mode ultrafin en polarisation HH. Les occupations du sol recherchées étaient le bâti dense, le bâti de densité moyenne, le bâti de densité faible, le bâti industriel et institutionnel, la végétation de faible densité, la végétation dense et l’eau. Les neuf paramètres de textures analysés ont été regroupés, en familles selon leur définition mathématique. Les paramètres de ressemblance/dissemblance regroupent l’Homogénéité, le Contraste, la Similarité et la Dissimilarité. Les paramètres de désordre sont l’Entropie et le Deuxième Moment Angulaire. L’Écart-Type et la Corrélation sont des paramètres de dispersion et la Moyenne est une famille à part. Il ressort des expériences que certaines combinaisons de paramètres de texture provenant de familles différentes utilisés dans les classifications donnent de très bons résultants alors que d’autres associations de paramètres de texture de définition mathématiques proches génèrent de moins bons résultats. Par ailleurs on constate que si l’utilisation de plusieurs paramètres de texture améliore les classifications, la performance de celle-ci plafonne à partir de trois paramètres. Malgré la bonne performance de cette approche basée sur la complémentarité des paramètres de texture, des erreurs systématiques dues aux effets cardinaux subsistent sur les classifications. Pour pallier à ce problème, il a été développé un modèle de compensation radiométrique basé sur la section efficace radar (SER). Une simulation radar à partir du modèle numérique de surface du milieu a permis d'extraire les zones de rétrodiffusion des bâtis et d'analyser les rétrodiffusions correspondantes. Une règle de compensation des effets cardinaux fondée uniquement sur les réponses des objets en fonction de leur orientation par rapport au plan d'illumination par le faisceau du radar a été mise au point. Des applications de cet algorithme sur des images RADARSAT-1 et RADARSAT-2 en polarisations HH, HV, VH, et VV ont permis de réaliser de considérables gains et d’éliminer l’essentiel des erreurs de classification dues aux effets cardinaux.

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.

A review of biochars’ potential role in the remediation, revegetation and restoration of contaminated soils

Biochars are biological residues combusted under low oxygen conditions, resulting in a porous, low density carbon rich material. Their large surface areas and cation exchange capacities, determined to a large extent by source materials and pyrolysis temperatures, enables enhanced sorption of both organic and inorganic contaminants to their surfaces, reducing pollutant mobility when amending contaminated soils. Liming effects or release of carbon into soil solution may increase arsenic mobility, whilst low capital but enhanced retention of plant nutrients can restrict revegetation on degraded soils amended only with biochars; the combination of composts, manures and other amendments with biochars could be their most effective deployment to soils requiring stabilisation by revegetation. Specific mechanisms of contaminant-biochar retention and release over time and the environmental impact of biochar amendments on soil organisms remain somewhat unclear but must be investigated to ensure that the management of environmental pollution coincides with ecological sustainability.

P-T-Fluid evolution and graphite deposition during retrograde metamorphism in Ribeira Fold Belt, SE Brazil: Oxygen fugacity, fluid inclusions and C-O-H isotopic evidence

Combined fluid inclusion (FI) microthermometry, Raman spectroscopy, X-ray diffraction, C-O-H isotopes and oxygen fugacities of granulites from central Ribeira Fold Belt, SE Brazil, provided the following results: i) Magnetite-Hematite fO(2) estimates range from 10(-11.5) bar (QFM + 1) to 10(-18.3) bar (QFM - 1) for the temperature range of 896 degrees C-656 degrees C, implying fO(2) decrease from metamorphic peak temperatures to retrograde conditions; ii) 5 main types of fluid inclusions were observed: a) CO(2) and CO(2)-N(2) (0-11 mol%) high to medium density (1.01-0.59 g/cm(3)) FI; b) CO(2) and CO(2)-N(2) (0-36 mol%) low density (0.19-0.29 g/cm(3)) FI; c) CO(2) (94-95 mol%)-N(2) (3 mol%)-CH(4) (2-3 mol%)-H(2)O (water phi(v) (25 degrees C) = 0.1) FI; d) low-salinity H(2)O-CO(2) FI; and e) late low-salinity H(2)O FI; iii) Raman analyses evidence two graphite types in khondalites: an early highly ordered graphite (T similar to 450 degrees C) overgrown by a disordered kind (T similar to 330 degrees C); iv) delta(18)O quartz results of 10.3-10.7 parts per thousand, imply high-temperature CO(2) delta(18)O values of 14.4-14.8 parts per thousand, suggesting the involvement of a metamorphic fluid, whereas lower temperature biotite delta(18)O and delta D results of 7.5-8.5 parts per thousand and -54 to -67 parts per thousand respectively imply H(2)O delta(18)O values of 10-11 parts per thousand and delta D(H2O) of -23 to -36 parts per thousand suggesting delta(18)O depletion and increasing fluid/rock ratio from metamorphic peak to retrograde conditions. Isotopic results are compatible with low-temperature H(2)O influx and fO(2) decrease that promoted graphite deposition in retrograde granulites, simultaneous with low density CO(2), CO(2)-N(2) and CO(2)-N(2)-CH(4)-H(2)O fluid inclusions at T = 450-330 degrees C. Graphite delta(13)C results of -10.9 to -11.4 parts per thousand imply CO(2) delta(13)C values of -0.8 to -1.3 parts per thousand suggesting decarbonation of Cambrian marine carbonates with small admixture of lighter biogenic or mantle derived fluids. Based on these results, it is suggested that metamorphic fluids from the central segment of Ribeira Fold Belt evolved to CO(2)-N(2) fluids during granulitic metamorphism at high fO(2), followed by rapid pressure drop at T similar to 400-450 degrees C during late exhumation that caused fO(2) reduction induced by temperature decrease and water influx, turning carbonic fluids into CO(2)-H(2)O (depleting biotite delta(18)O and delta D values), and progressively into H(2)O. When fO(2) decreased substantially by mixture of carbonic and aqueous fluids, graphite deposited forming khondalites. (C) 2010 Elsevier Ltd. All rights reserved.

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.

População de Metarhizium anisopliae em diferentes tipos e graus de compactação do solo

Este trabalho objetivou investigar a influência do tipo e compactação do solo na sobrevivência do fungo Metarhizium anisopliae. A sobrevivência do fungo foi determinada em quatro tipos de solos: Latossolo Vermelho textura argilosa, Latossolo Vermelho textura média, Argissolo Vermelho Amarelo textura arenosa média e Argissolo Vermelho Amarelo textura areno-argilosa, com maior teor de matéria orgânica. Para determinar o efeito da compactação na sobrevivência do fungo usaram-se os três primeiros tipos de solos nas densidades de 1,12, 1,32, 1,50g cm-3; 1,22, 1,44, 1,65g cm-3; 1,30, 1,50, 1,70g cm-3, respectivamente. Por meio da contagem de unidades formadoras de colônias (UFC) em placas de Petri, fizeram-se avaliações da sobrevivência do fungo, após zero, 20, 40, 60, 80, 100 e 120 dias de incubação a 27 ± 1ºC. Houve influência significativa do tipo de solo e do grau de compactação na sobrevivência do fungo, obtendo-se maior quantidade de UFC no solo textura areno-argilosa. Entre os demais solos, a maior sobrevivência ocorreu no solo textura arenosa e a menor no solo textura argilosa. O efeito da compactação foi significativo para o tipo de solo, exceto no solo textura arenosa. Independentemente do tipo de solo, a maior sobrevivência foi observada nos valores médios de densidade. A compactação teve maior impacto no solo textura média, onde ocorreu queda mais acentuada na quantidade de UFC em todas as densidades.

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.