Lattice boltzmann modeling of fluid flow to determine the permeability of a karst specimen

A limestone sample was scanned using computed tomography (CT) and the hydraulic conductivity of the 3D reconstructed sample was determined using Lattice- Boltzmann methods (LBM) at varying scales. Due to the shape and size of the original sample, it was challenging to obtain a consistent rectilinear test sample. Through visual inspection however, 91 mm and 76 mm samples were digitally cut from the original. The samples had porosities of 58% and 64% and produced hydraulic conductivity values of K= 13.5 m/s and K=34.5 m/s, respectively. Both of these samples were re-sampled to 1/8 and 1/64 of their original size to produce new virtual samples at lower resolutions of 0.542 mm/lu and 1.084 mm/lu, while still representing the same physical dimensions. The hydraulic conductivity tended to increase slightly as the resolution became coarser. In order to determine an REV, the 91 mm sample was also sub-sampled into blocks that were 1/8 and 1/64 the size of the original. The results were consistent with analytical expectations such as those produced by the Kozeny-Carman equation. A definitive REV size was not reached, however, indicating the need for a larger sample. The methods described here demonstrate the ability of LBM to test rock structures and sizes not normally attainable.

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.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

An irregular lattice model to simulate crack paths in bonded granular assemblies

Peer reviewed

An irregular lattice model to simulate crack paths in bonded granular assemblies

Peer reviewed

Family packing car for camping trip

Inscriptions: Verso: [stamped] Photograph by Freda Leinwand. [463 West Street, Studio 229G, New York, NY 10014].

Galois lattice dynamics

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Lattice approach to the dynamics of the DPSK-encoded soliton trains

We study the dynamical properties of the RZ-DPSK encoded sequences, focusing on the instabilities in the soliton train leading to the distortions of the information transmitted. The problem is reformulated within the framework of complex Toda chain model which allows one to carry out the simplified description of the optical soliton dynamics. We elucidate how the bit composition of the pattern affects the initial (linear) stage of the train dynamics and explain the general mechanisms of the appearance of unstable collective soliton modes. Then we discuss the nonlinear regime using asymptotic properties of the pulse stream at large propagation distances and analyze the dynamical behavior of the train classifying different scenarios for the pattern instabilities. Both approaches are based on the machinery of Hermitian and non-Hermitian lattice analysis. © 2010 IEEE.

In-situ examination of the selective etching of an alkanethiol monolayer covered Au{111} surface

An examination of the selective etching mechanism of a 1-alkanethiol self-assembled monolayer (SAM) covered Au{111} surface using in-situ atomic force microscopy (AFM) and molecular resolution scanning tunnelling microscopy (STM) is presented. The monolayer self-assembles on a smooth Au{111} surface and typically contains nanoscale non-uniformities such as pinholes, domain boundaries and monatomic depressions. During etching in a ferri/ferrocyanide water-based etchant, selective and preferential etching occurs at SAM covered Au(111) terrace and step edges where a lower SAM packing density is observed, resulting in triangular islands being relieved. The triangular islands are commensurate with the Au(111) lattice with their long edges parallel to its [11-0] direction. Thus, SAM etching is selective and preferential attack is localized to defects and step edges at sites of high molecular density contrast.

Lattice Engineering of III-Nitride Heterostructures and Their Applications

III-Nitride materials have recently become a promising candidate for superior applications over the current technologies. However, certain issues such as lack of native substrates, and high defect density have to be overcome for further development of III-Nitride technology. This work presents research on lattice engineering of III-Nitride materials, and the structural, optical, and electrical properties of its alloys, in order to approach the ideal material for various applications. We demonstrated the non-destructive and quantitative characterization of composition modulated nanostructure in InAlN thin films with X-ray diffraction. We found the development of the nanostructure depends on growth temperature, and the composition modulation has impacts on carrier recombination dynamics. We also showed that the controlled relaxation of a very thin AlN buffer (20 ~ 30 nm) or a graded composition InGaN buffer can significantly reduce the defect density of a subsequent epitaxial layer. Finally, we synthesized an InAlGaN thin films and a multi-quantum-well structure. Significant emission enhancement in the UVB range (280 – 320 nm) was observed compared to AlGaN thin films. The nature of the enhancement was investigated experimentally and numerically, suggesting carrier confinement in the In localization centers.

Galois lattice dynamics

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

The Plastic Coupled Map Lattice: a Novel Image‐processing Paradigm

Coupled map lattices (CML) can describe many relaxation and optimization algorithms currently used in image processing. We recently introduced the ‘‘plastic‐CML’’ as a paradigm to extract (segment) objects in an image. Here, the image is applied by a set of forces to a metal sheet which is allowed to undergo plastic deformation parallel to the applied forces. In this paper we present an analysis of our ‘‘plastic‐CML’’ in one and two dimensions, deriving the nature and stability of its stationary solutions. We also detail how to use the CML in image processing, how to set the system parameters and present examples of it at work. We conclude that the plastic‐CML is able to segment images with large amounts of noise and large dynamic range of pixel values, and is suitable for a very large scale integration(VLSI) implementation.

Aspects of Chiral Symmetry Breaking in Lattice QCD

Thesis (Ph.D.)--University of Washington, 2016-08

Loess as a Collapsible Soil: Some Basic Particle Packing

Loess is the most important collapsible soil; possibly the only engineering soil in which real collapse occurs. A real collapse involves a diminution in volume - it would be an open metastable packing being reduced to a more closely packed, more stable structure. Metastability is at the heart of the collapsible soils problem. To envisage and to model the collapse process in a metastable medium, knowledge is required about the nature and shape of the particles, the types of packings they assume (real and ideal), and the nature of the collapse process - a packing transition upon a change to the effective stress in a media of double porosity. Particle packing science has made little progress in geoscience discipline - since the initial packing paradigms set by Graton and Fraser (1935) - nevertheless is relatively well-established in the soft matter physics discipline. The collapse process can be represented by mathematical modelling of packing – including the Monte Carlo simulations - but relating representation to process remains difficult. This paper revisits the problem of sudden packing transition from a micro-physico-mechanical viewpoint (i.e. collapse imetan terms of structure-based effective stress). This cross-disciplinary approach helps in generalization on collapsible soils to be made that suggests loess is the only truly collapsible soil, because it is only loess which is so totally influenced by the packing essence of the formation process.