Fractal Behaviour and Irregularity of Wetting Fronts in Heterogeneous Porous Media

The wetting front is the zone where water invades and advances into an initially dry porous material and it plays a crucial role in solute transport through the unsaturated zone. Water is an essential part of the physiological process of all plants. Through water, necessary minerals are moved from the roots to the parts of the plants that require them. Water moves chemicals from one part of the plant to another. It is also required for photosynthesis, for metabolism and for transpiration. The leaching of chemicals by wetting fronts is influenced by two major factors, namely: the irregularity of the fronts and heterogeneity in the distribution of chemicals, both of which have been described by using fractal techniques. Soil structure can significantly modify infiltration rates and flow pathways in soils. Relations between features of soil structure and features of infiltration could be elucidated from the velocities and the structure of wetting fronts. When rainwater falls onto soil, it doesn?t just pool on surfaces. Water ?or another fluid- acts differently on porous surfaces. If the surface is permeable (porous) it seeps down through layers of soil, filling that layer to capacity. Once that layer is filled, it moves down into the next layer. In sandy soil, water moves quickly, while it moves much slower through clay soil. The movement of water through soil layers is called the the wetting front. Our research concerns the motion of a liquid into an initially dry porous medium. Our work presents a theoretical framework for studying the physical interplay between a stationary wetting front of fractal dimension D with different porous materials. The aim was to model the mass geometry interplay by using the fractal dimension D of a stationary wetting front. The plane corresponding to the image is divided in several squares (the minimum correspond to the pixel size) of size length ". We acknowledge the help of Prof. M. García Velarde and the facilities offered by the Pluri-Disciplinary Institute of the Complutense University of Madrid. We also acknowledge the help of European Community under project Multi-scale complex fluid flows and interfacial phenomena (PITN-GA-2008-214919). Thanks are also due to ERCOFTAC (PELNoT, SIG 14)

Accelerating CFD via local POD plus Galerkin projections

Se desarrollan varias técnicas basadas en descomposición ortogonal propia (DOP) local y proyección de tipo Galerkin para acelerar la integración numérica de problemas de evolución, de tipo parabólico, no lineales. Las ideas y métodos que se presentan conllevan un nuevo enfoque para la modelización de tipo DOP, que combina intervalos temporales cortos en que se usa un esquema numérico estándard con otros intervalos temporales en que se utilizan los sistemas de tipo Galerkin que resultan de proyectar las ecuaciones de evolución sobre la variedad lineal generada por los modos DOP, obtenidos a partir de instantáneas calculadas en los intervalos donde actúa el código numérico. La variedad DOP se construye completamente en el primer intervalo, pero solamente se actualiza en los demás intervalos según las dinámicas de la solución, aumentando de este modo la eficiencia del modelo de orden reducido resultante. Además, se aprovechan algunas propiedades asociadas a la dependencia débil de los modos DOP tanto en la variable temporal como en los posibles parámetros de que pueda depender el problema. De esta forma, se aumentan la flexibilidad y la eficiencia computacional del proceso. La aplicación de los métodos resultantes es muy prometedora, tanto en la simulación de transitorios en flujos laminares como en la construcción de diagramas de bifurcación en sistemas dependientes de parámetros. Las ideas y los algoritmos desarrollados en la tesis se ilustran en dos problemas test, la ecuación unidimensional compleja de Ginzburg-Landau y el problema bidimensional no estacionario de la cavidad. Abstract Various ideas and methods involving local proper orthogonal decomposition (POD) and Galerkin projection are presented aiming at accelerating the numerical integration of nonlinear time dependent parabolic problems. The proposed methods come from a new approach to the POD-based model reduction procedures, which combines short runs with a given numerical solver and a reduced order model constructed by expanding the solution of the problem into appropriate POD modes, which span a POD manifold, and Galerkin projecting some evolution equations onto that linear manifold. The POD manifold is completely constructed from the outset, but only updated as time proceeds according to the dynamics, which yields an adaptive and flexible procedure. In addition, some properties concerning the weak dependence of the POD modes on time and possible parameters in the problem are exploited in order to increase the flexibility and efficiency of the low dimensional model computation. Application of the developed techniques to the approximation of transients in laminar fluid flows and the simulation of attractors in bifurcation problems shows very promising results. The test problems considered to illustrate the various ideas and check the performance of the algorithms are the onedimensional complex Ginzburg-Landau equation and the two-dimensional unsteady liddriven cavity problem.

Order 10 4 speedup in global linear instability analysis using matrix formation

A unified solution framework is presented for one-, two- or three-dimensional complex non-symmetric eigenvalue problems, respectively governing linear modal instability of incompressible fluid flows in rectangular domains having two, one or no homogeneous spatial directions. The solution algorithm is based on subspace iteration in which the spatial discretization matrix is formed, stored and inverted serially. Results delivered by spectral collocation based on the Chebyshev-Gauss-Lobatto (CGL) points and a suite of high-order finite-difference methods comprising the previously employed for this type of work Dispersion-Relation-Preserving (DRP) and Padé finite-difference schemes, as well as the Summationby- parts (SBP) and the new high-order finite-difference scheme of order q (FD-q) have been compared from the point of view of accuracy and efficiency in standard validation cases of temporal local and BiGlobal linear instability. The FD-q method has been found to significantly outperform all other finite difference schemes in solving classic linear local, BiGlobal, and TriGlobal eigenvalue problems, as regards both memory and CPU time requirements. Results shown in the present study disprove the paradigm that spectral methods are superior to finite difference methods in terms of computational cost, at equal accuracy, FD-q spatial discretization delivering a speedup of ð (10 4). Consequently, accurate solutions of the three-dimensional (TriGlobal) eigenvalue problems may be solved on typical desktop computers with modest computational effort.

Fresnel-based modular solar fields for performance/cost optimization in solar thermal power plants: A comparison with parabolic trough collectors.

Linear Fresnel collectors are identified as a technology that should play a main role in order to reduce cost of Concentrating Solar Power. An optical and thermal analysis of the different blocks of the solar power plant is carried out, where Fresnel arrays are compared with the most extended linear technology: parabolic trough collectors. It is demonstrated that the optical performance of Fresnel array is very close to that of PTC, with similar values of maximum flux intensities. In addition, if the heat carrier fluid flows in series by the tubes of the receiver, relatively high thermal efficiencies are achieved. Thus, an annual solar to electricity efficiency of 19% is expected, which is similar to the state of the art in PTCs; this is done with a reduction of costs, thanks to lighter structures, that drives to an estimation of LCOE of around 6.5 c€/kWh.

Stable isotope ratios of interstitial fluids from ODP Leg 110 sites

Delta18O values of pore waters from the northern Barbados accretionary prism range from -0.3 to -3.6? and reflect pervasive reaction of volcanic ash to form smectite within the sedimentary sequence and continued low temperature alteration of basalt in the underlying ocean crust with the overprint of diffusive exchange between water in the sediment pores and the open ocean. Delta D values of pore waters in sediments sampled seaward of the deformation front drop from +5? at the sediment surface to -6? at the deepest levels sampled. These changes may also be related to alteration processes but remain largely enigmatic. Sediment deformation caused by impingement of the Caribbean plate on the Atlantic plate has instigated migration of chemically and isotopically distinct fluid along faults and coarse-grained sedimentary beds; delta18O values of pore waters are also locally affected by thrust stacking which increases diffusive pathlengths and possibly modifies diagenetic reaction rates in Pleistocene sediments. Migrating fluids are distinguished by anomalous delta18O values that are as much as 1? higher than those of surrounding fluids. Uncertainties in hydrogen isotope fractionation resulting from processes occurring under these conditions hinder identification of the hydrogen isotope composition of expelled fluid. Stable isotope analyses of pore waters help constrain the fluid migration history of the accretionary prism by limiting the source of fluids, the paths along which fluid flows, and the timing of faulting and subsequent fluid flow.

Nonlinear dynamics of cilia and flagella

Cilia and flagella are hairlike extensions of eukaryotic cells which generate oscillatory beat patterns that can propel micro-organisms and create fluid flows near cellular surfaces. The evolutionary highly conserved core of cilia and flagella consists of a cylindrical arrangement of nine microtubule doublets, called the axoneme. The axoneme is an actively bending structure whose motility results from the action of dynein motor proteins cross-linking microtubule doublets and generating stresses that induce bending deformations. The periodic beat patterns are the result of a mechanical feedback that leads to self-organized bending waves along the axoneme. Using a theoretical framework to describe planar beating motion, we derive a nonlinear wave equation that describes the fundamental Fourier mode of the axonemal beat. We study the role of nonlinearities and investigate how the amplitude of oscillations increases in the vicinity of an oscillatory instability. We furthermore present numerical solutions of the nonlinear wave equation for different boundary conditions. We find that the nonlinear waves are well approximated by the linearly unstable modes for amplitudes of beat patterns similar to those observed experimentally.

Concentração do óleo de girassol em compostos insaturados utilizando destilação molecular

Vegetable oils are characterized as important raw materials in the supplying of natural substances of interest pharmaceutical, food and cosmetic industry. Sunflower oil stands out for its important composition present in unsaturated fatty acids such as oleic acid (C18:1) and linoleic (C18:2), responsible for many health benefits. The main objective of this study is obtain enriched fractions in unsaturated compounds from refined sunflower oil. The oil used in this study was characterized by the determination of some properties, like iodine number, acid number and viscosity. A transesterification was done to transform the triglycerides into their corresponding methyl esters of fatty acids. These was submitted the molecular distillation process, for present as an efficient alternative to separation and purification of these substances, using high vacuum and low temperatures. Of the esters fractions that was obtained, were analyzed by gas chromatography. The experimental design technique was used to evaluate the influence of the temperature variation of evaporation and condensation system on the percentage obtained residue. The evaporator temperature proved to be the most influential variable on the studied response. The optimized conditions for the answer was studied at 100 °C for evaporator temperature and 10 °C for the condenser temperature. The graph of "split ratio" showed that for the lowest flow feed (1 mL/min) and higher evaporator temperature (110 °C) was obtained in the largest fraction of distillate. It also used the study of the influence of evaporator temperature on the concentration of unsaturated compounds. The best operating conditions for temperature was 90 °C reached 82.21 % of unsaturated compounds. Elimination curves of the unsaturated compounds present in the distillate stream were obtained. The simulation results of the molecular distillation process of sunflower oil showed the concentration profiles for three different feed flow rates. The speed, temperature and thickness profiles of the liquid film were obtained. The speed of the film increases as the fluid flows through the walls of the evaporator, reaching a maximum on length of 0.075 m. The film thickness decreases on the route, since many compounds are volatilized. The result of the temperature profile had to be consistent with the literature reproduced, being constant after reaching the maximum operating temperature in the length of 0.15 m. This study allowed characterizing and focusing, through experimental analysis, unsaturated compounds and observing the sunflower oil´s behavior through process simulation.

Concentração do óleo de girassol em compostos insaturados utilizando destilação molecular

Vegetable oils are characterized as important raw materials in the supplying of natural substances of interest pharmaceutical, food and cosmetic industry. Sunflower oil stands out for its important composition present in unsaturated fatty acids such as oleic acid (C18:1) and linoleic (C18:2), responsible for many health benefits. The main objective of this study is obtain enriched fractions in unsaturated compounds from refined sunflower oil. The oil used in this study was characterized by the determination of some properties, like iodine number, acid number and viscosity. A transesterification was done to transform the triglycerides into their corresponding methyl esters of fatty acids. These was submitted the molecular distillation process, for present as an efficient alternative to separation and purification of these substances, using high vacuum and low temperatures. Of the esters fractions that was obtained, were analyzed by gas chromatography. The experimental design technique was used to evaluate the influence of the temperature variation of evaporation and condensation system on the percentage obtained residue. The evaporator temperature proved to be the most influential variable on the studied response. The optimized conditions for the answer was studied at 100 °C for evaporator temperature and 10 °C for the condenser temperature. The graph of "split ratio" showed that for the lowest flow feed (1 mL/min) and higher evaporator temperature (110 °C) was obtained in the largest fraction of distillate. It also used the study of the influence of evaporator temperature on the concentration of unsaturated compounds. The best operating conditions for temperature was 90 °C reached 82.21 % of unsaturated compounds. Elimination curves of the unsaturated compounds present in the distillate stream were obtained. The simulation results of the molecular distillation process of sunflower oil showed the concentration profiles for three different feed flow rates. The speed, temperature and thickness profiles of the liquid film were obtained. The speed of the film increases as the fluid flows through the walls of the evaporator, reaching a maximum on length of 0.075 m. The film thickness decreases on the route, since many compounds are volatilized. The result of the temperature profile had to be consistent with the literature reproduced, being constant after reaching the maximum operating temperature in the length of 0.15 m. This study allowed characterizing and focusing, through experimental analysis, unsaturated compounds and observing the sunflower oil´s behavior through process simulation.

Automatic dynamic mask extraction for PIV images containing an unsteady interface, bubbles, and a moving structure

When performing Particle Image Velocimetry (PIV) measurements in complex fluid flows with moving interfaces and a two-phase flow, it is necessary to develop a mask to remove non-physical measurements. This is the case when studying, for example, the complex bubble sweep-down phenomenon observed in oceanographic research vessels. Indeed, in such a configuration, the presence of an unsteady free surface, of a solid–liquid interface and of bubbles in the PIV frame, leads to generate numerous laser reflections and therefore spurious velocity vectors. In this note, an image masking process is developed to successively identify the boundaries of the ship and the free surface interface. As the presence of the solid hull surface induces laser reflections, the hull edge contours are simply detected in the first PIV frame and dynamically estimated for consecutive ones. As for the unsteady surface determination, a specific process is implemented like the following: i) the edge detection of the gradient magnitude in the PIV frame, ii) the extraction of the particles by filtering high-intensity large areas related to the bubbles and/or hull reflections, iii) the extraction of the rough region containing these particles and their reflections, iv) the removal of these reflections. The unsteady surface is finally obtained with a fifth-order polynomial interpolation. The resulted free surface is successfully validated from the Fourier analysis and by visualizing selected PIV images containing numerous spurious high intensity areas. This paper demonstrates how this data analysis process leads to PIV images database without reflections and an automatic detection of both the free surface and the rigid body. An application of this new mask is finally detailed, allowing a preliminary analysis of the hydrodynamic flow.

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems