Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Uranium distribution as a proxy for basin-scale fluid flow in distributive fluvial systems

This work was supported by the Fluvial Systems Research Group sponsors BG Group, BP, Chevron, ConocoPhilips, and Total. We thank reviews from Martin Stokes, an anonymous reviewer and Editor Stuart Jones.

Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature-induced grain coarsening. Despite this large-scale trend, a considerable grain size scatter occurs on the outcrop-scale indicating local influence of second-order effects such as thermal perturbations, fluid flow and second-phase particles. Second-phase particles, whose sizes range from nano- to the micron-scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron-scale second-phase particles, or to samples containing a large number of nano-scale particles. The obtained data set suggests that the second phases induce a temperature-controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D 1⁄4 C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second-phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian-type distributions of the pure samples. In contrast, fluid-enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable d18O and d13C isotope ratios in fluid-affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid-induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1.To summarize, considerable grain size variations of up to one order of magnitude can locally result from second-order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

EGCL: an extended G-Code Language with flow control, functions and mnemonic variables

In the context of computer numerical control (CNC) and computer aided manufacturing (CAM), the capabilities of programming languages such as symbolic and intuitive programming, program portability and geometrical portfolio have special importance -- They allow to save time and to avoid errors during part programming and permit code re-usage -- Our updated literature review indicates that the current state of art presents voids in parametric programming, program portability and programming flexibility -- In response to this situation, this article presents a compiler implementation for EGCL (Extended G-code Language), a new, enriched CNC programming language which allows the use of descriptive variable names, geometrical functions and flow-control statements (if-then-else, while) -- Our compiler produces low-level generic, elementary ISO-compliant Gcode, thus allowing for flexibility in the choice of the executing CNC machine and in portability -- Our results show that readable variable names and flow control statements allow a simplified and intuitive part programming and permit re-usage of the programs -- Future work includes allowing the programmer to define own functions in terms of EGCL, in contrast to the current status of having them as library built-in functions

On-demand generation and removal of alginate biocompatible microvalves for flow control in microfluidics

[EN] This paper describes, for the first time, the use of alginate hydrogels as miniaturised microvalves within microfluidic devices. These biocompatible and biodegradable microvalves are generated in situ and on demand, allowing for microfluidic flow control. The microfluidic devices were fabricated using an origami inspired technique of folding several layers of cyclic olefin polymer followed by thermocompression bonding. The hydrogels can be dehydrated at mild temperatures, 37◦C, to slightly open the microvalve and chemically erased using an ethylenediaminetetraacetic acid disodium salt (EDTA) solution, to completely open the channel, ensuring the reusability of the whole device and removal of damaged or defective valves for subsequent regeneration.

Methane and sulphate consumption in different habitats of the Håkon Mosby Mud Volcano (HMMV)

The Hakon Mosby Mud Volcano is a highly active methane seep hosting different chemosynthetic communities such as thiotrophic bacterial mats and siboglinid tubeworm assemblages. This study focuses on in situ measurements of methane fluxes to and from these different habitats, in comparison to benthic methane and oxygen consumption rates. By quantifying in situ oxygen, methane, and sulfide fluxes in different habitats, a spatial budget covering areas of 10-1000 -m diameter was established. The range of dissolved methane efflux (770-2 mmol m-2 d-1) from the center to the outer rim was associated with a decrease in temperature gradients from 46°C to < 1°C m-1, indicating that spatial variations in fluid flow control the distribution of benthic habitats and activities. Accordingly, total oxygen uptake (TOU) varied between the different habitats by one order of magnitude from 15 mmol m-2 d-1 to 161 mmol m-2 d-1. High fluid flow rates appeared to suppress benthic activities by limiting the availability of electron acceptors. Accordingly, the highest TOU was associated with the lowest fluid flow and methane efflux. This was most likely due to the aerobic oxidation of methane, which may be more relevant as a sink for methane as previously considered in submarine ecosystems.

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Shock control bump design optimization on Natural Laminar Aerofoil

The chapter investigates Shock Control Bumps (SCB) on a Natural Laminar Flow (NLF) aerofoil; RAE 5243 for Active Flow Control (AFC). A SCB approach is used to decelerate supersonic flow on the suction/pressure sides of transonic aerofoil that leads delaying shock occurrence or weakening of shock strength. Such an AFC technique reduces significantly the total drag at transonic speeds. This chapter considers the SCB shape design optimisation at two boundary layer transition positions (0 and 45%) using an Euler software coupled with viscous boundary layer effects and robust Evolutionary Algorithms (EAs). The optimisation method is based on a canonical Evolution Strategy (ES) algorithm and incorporates the concepts of hierarchical topology and parallel asynchronous evaluation of candidate solution. Two test cases are considered with numerical experiments; the first test deals with a transition point occurring at the leading edge and the transition point is fixed at 45% of wing chord in the second test. Numerical results are presented and it is demonstrated that an optimal SCB design can be found to significantly reduce transonic wave drag and improves lift on drag (L/D) value when compared to the baseline aerofoil design.

Thermal profile of a near-adiabatic compression process in a cylindrical tube and establishment of critical control elements for repeatable process control

The compressed gas industry and government agencies worldwide utilize "adiabatic compression" testing for qualifying high-pressure valves, regulators, and other related flow control equipment for gaseous oxygen service. This test methodology is known by various terms including adiabatic compression testing, gaseous fluid impact testing, pneumatic impact testing, and BAM testing as the most common terms. The test methodology will be described in greater detail throughout this document but in summary it consists of pressurizing a test article (valve, regulator, etc.) with gaseous oxygen within 15 to 20 milliseconds (ms). Because the driven gas1 and the driving gas2 are rapidly compressed to the final test pressure at the inlet of the test article, they are rapidly heated by the sudden increase in pressure to sufficient temperatures (thermal energies) to sometimes result in ignition of the nonmetallic materials (seals and seats) used within the test article. In general, the more rapid the compression process the more "adiabatic" the pressure surge is presumed to be and the more like an isentropic process the pressure surge has been argued to simulate. Generally speaking, adiabatic compression is widely considered the most efficient ignition mechanism for directly kindling a nonmetallic material in gaseous oxygen and has been implicated in many fire investigations. Because of the ease of ignition of many nonmetallic materials by this heating mechanism, many industry standards prescribe this testing. However, the results between various laboratories conducting the testing have not always been consistent. Research into the test method indicated that the thermal profile achieved (i.e., temperature/time history of the gas) during adiabatic compression testing as required by the prevailing industry standards has not been fully modeled or empirically verified, although attempts have been made. This research evaluated the following questions: 1) Can the rapid compression process required by the industry standards be thermodynamically and fluid dynamically modeled so that predictions of the thermal profiles be made, 2) Can the thermal profiles produced by the rapid compression process be measured in order to validate the thermodynamic and fluid dynamic models; and, estimate the severity of the test, and, 3) Can controlling parameters be recommended so that new guidelines may be established for the industry standards to resolve inconsistencies between various test laboratories conducting tests according to the present standards?

Similarity solutions for flow and heat transfer of non-Newtonian fluid over a stretching surface

Similarity solutions are carried out for flow of power law non-Newtonian fluid film on unsteady stretching surface subjected to constant heat flux. Free convection heat transfer induces thermal boundary layer within a semi-infinite layer of Boussinesq fluid. The nonlinear coupled partial differential equations (PDE) governing the flow and the boundary conditions are converted to a system of ordinary differential equations (ODE) using two-parameter groups. This technique reduces the number of independent variables by two, and finally the obtained ordinary differential equations are solved numerically for the temperature and velocity using the shooting method. The thermal and velocity boundary layers are studied by the means of Prandtl number and non-Newtonian power index plotted in curves.

Unsteady-Flow Of A Viscous-Fluid Between 2 Parallel Disks With A Time-Varying Gap Width And A Magnetic-Field

The unsteady incompressible viscous fluid flow between two parallel infinite disks which are located at a distance h(t*) at time t* has been studied. The upper disk moves towards the lower disk with velocity h'(t*). The lower disk is porous and rotates with angular velocity Omega(t*). A magnetic field B(t*) is applied perpendicular to the two disks. It has been found that the governing Navier-Stokes equations reduce to a set of ordinary differential equations if h(t*), a(t*) and B(t*) vary with time t* in a particular manner, i.e. h(t*) = H(1 - alpha t*)(1/2), Omega(t*) = Omega(0)(1 - alpha t*)(-1), B(t*) = B-0(1 - alpha t*)(-1/2). These ordinary differential equations have been solved numerically using a shooting method. For small Reynolds numbers, analytical solutions have been obtained using a regular perturbation technique. The effects of squeeze Reynolds numbers, Hartmann number and rotation of the disk on the flow pattern, normal force or load and torque have been studied in detail

Non-disruptive and null-deflection mass flow measurement by a pressure compensation technique

Accurate mass flow measurement is very important in various monitoring and control applications. This paper proposes a novel method of fluid flow measurement by compensating the pressure drop across the ends of measuring unit using a compensating pump. The pressure drop due to the flow is balanced by a feedback control loop. This is a null-deflection type of measurement. As the insertion of such a measuring unit does not affect the functioning of the systems, this is also a non-disruptive flow measurement method. The implementation and design of such a unit are discussed. The system is modeled and simulated using the bond graph technique and it is experimentally validated. (C) 2009 Elsevier Ltd. All rights reserved.