Lattice Boltzmann Method for 3-D Flows with Curved Boundary
NASA Technical Reports Server (NTRS)
Mei, Renwei; Shyy, Wei; Yu, Dazhi; Luo, Li-Shi
2002-01-01
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
Multiple-Relaxation-Time Lattice Boltzmann Models in 3D
NASA Technical Reports Server (NTRS)
dHumieres, Dominique; Ginzburg, Irina; Krafczyk, Manfred; Lallemand, Pierre; Luo, Li-Shi; Bushnell, Dennis M. (Technical Monitor)
2002-01-01
This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of fifteen-velocity and nineteen-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in three dimensions at Re=500 and 2000 is performed. The results clearly demonstrate the superior numerical stability of the multiple-relaxation-time lattice Boltzmann equation over the popular lattice Bhatnagar-Gross-Krook equation.
A lattice-Boltzmann scheme of the Navier-Stokes equations on a 3D cuboid lattice
NASA Astrophysics Data System (ADS)
Min, Haoda; Peng, Cheng; Wang, Lian-Ping
2015-11-01
The standard lattice-Boltzmann method (LBM) for fluid flow simulation is based on a square (in 2D) or cubic (in 3D) lattice grids. Recently, two new lattice Boltzmann schemes have been developed on a 2D rectangular grid using the MRT (multiple-relaxation-time) collision model, by adding a free parameter in the definition of moments or by extending the equilibrium moments. Here we developed a lattice Boltzmann model on 3D cuboid lattice, namely, a lattice grid with different grid lengths in different spatial directions. We designed our MRT-LBM model by matching the moment equations from the Chapman-Enskog expansion with the Navier-Stokes equations. The model guarantees correct hydrodynamics. A second-order term is added to the equilibrium moments in order to restore the isotropy of viscosity on a cuboid lattice. The form and the coefficients of the extended equilibrium moments are determined through an inverse design process. An additional benefit of the model is that the viscosity can be adjusted independent of the stress-moment relaxation parameter, thus improving the numerical stability of the model. The resulting cuboid MRT-LBM model is then validated through benchmark simulations using laminar channel flow, turbulent channel flow, and the 3D Taylor-Green vortex flow.
Development and validation of a 3D Lattice Boltzmann model for volcano aeroacoustics
NASA Astrophysics Data System (ADS)
Brogi, Federico; Bonadonna, Costanza; Ripepe, Maurizio; Chopard, Bastien; Malaspinas, Orestis; Latt, Jonas; Falcone, Jean-Luc
2015-04-01
Infrasound measurements have a great potential for the real time characterization of volcanic plume source parameters [Ripepe et al., 2013]. Nonetheless many shortcomings have been highlighted in the understanding of the infrasound monitoring. In particular, the application of the classical acoustic source models to volcanic explosive eruptions has shown to be challenging and a better knowledge of the link between the acoustic radiation and actual volcanic fluid dynamics processes is required. New insights into this subject could be given by the study of realistic aeroacoustic numerical simulations of a volcanic jet. Our work mainly focuses on developing and validating such numerical model to determine when and if classical model source theory can be applied to explain volcanic infrasound data. Lattice Boltzmann strategies (LB) provide the opportunity to develop an accurate, computationally fast, 3D physical model for a volcanic jet and wave propagation. In the field of aeroacoustic applications, dedicated LB schemes has been proven to have the low dispersion and dissipative properties needed for capturing the weak acoustic pressure fluctuations. However, when dealing with simulations of realistic flows, artificial boundaries are defined around the flow region. The reflected waves from these boundaries can have significant influence on the flow field and overwhelm the acoustic field of interest. A special absorbing boundary layer has been implemented in our model to suppress the reflected waves [Xu et al., 2013]. In addition, for highly multi-scale turbulent flows, such as volcanic plumes, the number of grid points needed to represent the smallest scales might become intractable and the most complicated physics happen only in small portions of the computational domain. The implementation of the grid refinement, in our model allow us to insert local finer grids only where is actually needed [Lagrava et al., 2012] and to increase the size of the computational domain
Lattice Boltzmann simulation of dynamics of plunge and pitch of 3D flexible wing
NASA Astrophysics Data System (ADS)
Qi, Dewei; Shyy, Wei
2008-11-01
The method of lattice Boltzmann (LB) simulation has been used to simulate fluid structures and motion of a flexible insect wing in a 3D space. In the method, a beam has been discretized into a chain of rigid segments. Each segment is connected through ball and socket joints at its ends. One segment may be bent and twisted with its neighboring segment. A constraint force is applied to each joint to ensure the solid structure moving as a whole flexible elastic body.We have demonstrated that the LB method is suitable for modeling of aerodynamics of insects flight at low Reynolds numbers. First, a simulation of plunging and pitching of a rigid wing is performed at Re=75 in a 2D space and the results of lift forces and flow structures are in excellent agreement with the previous results. Second, plunging and pitching of a flexible wing in span-wise direction is simulated at Re=136 in a 3D space. We found that when twisting elasticity is large enough the twisting angle could be controlled at a level of smaller than 0.2 degree. It is shown that as bending and twisting elasticity is large enough, the motion of flexible wing approaches that of a rigid membrane wing. The simulation results show that the optimization of flexibility in span-wise direction will benefit thrust and an intermediate level is favorable. The results are consistent with experimental finding.
Lattice Boltzmann Model of 3D Multiphase Flow in Artery Bifurcation Aneurysm Problem
Abas, Aizat; Mokhtar, N. Hafizah; Ishak, M. H. H.; Abdullah, M. Z.; Ho Tian, Ang
2016-01-01
This paper simulates and predicts the laminar flow inside the 3D aneurysm geometry, since the hemodynamic situation in the blood vessels is difficult to determine and visualize using standard imaging techniques, for example, magnetic resonance imaging (MRI). Three different types of Lattice Boltzmann (LB) models are computed, namely, single relaxation time (SRT), multiple relaxation time (MRT), and regularized BGK models. The results obtained using these different versions of the LB-based code will then be validated with ANSYS FLUENT, a commercially available finite volume- (FV-) based CFD solver. The simulated flow profiles that include velocity, pressure, and wall shear stress (WSS) are then compared between the two solvers. The predicted outcomes show that all the LB models are comparable and in good agreement with the FVM solver for complex blood flow simulation. The findings also show minor differences in their WSS profiles. The performance of the parallel implementation for each solver is also included and discussed in this paper. In terms of parallelization, it was shown that LBM-based code performed better in terms of the computation time required. PMID:27239221
Lattice Boltzmann Model of 3D Multiphase Flow in Artery Bifurcation Aneurysm Problem.
Abas, Aizat; Mokhtar, N Hafizah; Ishak, M H H; Abdullah, M Z; Ho Tian, Ang
2016-01-01
This paper simulates and predicts the laminar flow inside the 3D aneurysm geometry, since the hemodynamic situation in the blood vessels is difficult to determine and visualize using standard imaging techniques, for example, magnetic resonance imaging (MRI). Three different types of Lattice Boltzmann (LB) models are computed, namely, single relaxation time (SRT), multiple relaxation time (MRT), and regularized BGK models. The results obtained using these different versions of the LB-based code will then be validated with ANSYS FLUENT, a commercially available finite volume- (FV-) based CFD solver. The simulated flow profiles that include velocity, pressure, and wall shear stress (WSS) are then compared between the two solvers. The predicted outcomes show that all the LB models are comparable and in good agreement with the FVM solver for complex blood flow simulation. The findings also show minor differences in their WSS profiles. The performance of the parallel implementation for each solver is also included and discussed in this paper. In terms of parallelization, it was shown that LBM-based code performed better in terms of the computation time required. PMID:27239221
Dual FIB-SEM 3D imaging and lattice boltzmann modeling of porosimetry and multiphase flow in chalk.
Rinehart, Alex; Petrusak, Robin; Heath, Jason E.; Dewers, Thomas A.; Yoon, Hongkyu
2010-12-01
Mercury intrusion porosimetry (MIP) is an often-applied technique for determining pore throat distributions and seal analysis of fine-grained rocks. Due to closure effects, potential pore collapse, and complex pore network topologies, MIP data interpretation can be ambiguous, and often biased toward smaller pores in the distribution. We apply 3D imaging techniques and lattice-Boltzmann modeling in interpreting MIP data for samples of the Cretaceous Selma Group Chalk. In the Mississippi Interior Salt Basin, the Selma Chalk is the apparent seal for oil and gas fields in the underlying Eutaw Fm., and, where unfractured, the Selma Chalk is one of the regional-scale seals identified by the Southeast Regional Carbon Sequestration Partnership for CO2 injection sites. Dual focused ion - scanning electron beam and laser scanning confocal microscopy methods are used for 3D imaging of nanometer-to-micron scale microcrack and pore distributions in the Selma Chalk. A combination of image analysis software is used to obtain geometric pore body and throat distributions and other topological properties, which are compared to MIP results. 3D data sets of pore-microfracture networks are used in Lattice Boltzmann simulations of drainage (wetting fluid displaced by non-wetting fluid via the Shan-Chen algorithm), which in turn are used to model MIP procedures. Results are used in interpreting MIP results, understanding microfracture-matrix interaction during multiphase flow, and seal analysis for underground CO2 storage.
LBM-EP: Lattice-Boltzmann method for fast cardiac electrophysiology simulation from 3D images.
Rapaka, S; Mansi, T; Georgescu, B; Pop, M; Wright, G A; Kamen, A; Comaniciu, Dorin
2012-01-01
Current treatments of heart rhythm troubles require careful planning and guidance for optimal outcomes. Computational models of cardiac electrophysiology are being proposed for therapy planning but current approaches are either too simplified or too computationally intensive for patient-specific simulations in clinical practice. This paper presents a novel approach, LBM-EP, to solve any type of mono-domain cardiac electrophysiology models at near real-time that is especially tailored for patient-specific simulations. The domain is discretized on a Cartesian grid with a level-set representation of patient's heart geometry, previously estimated from images automatically. The cell model is calculated node-wise, while the transmembrane potential is diffused using Lattice-Boltzmann method within the domain defined by the level-set. Experiments on synthetic cases, on a data set from CESC'10 and on one patient with myocardium scar showed that LBM-EP provides results comparable to an FEM implementation, while being 10 - 45 times faster. Fast, accurate, scalable and requiring no specific meshing, LBM-EP paves the way to efficient and detailed models of cardiac electrophysiology for therapy planning. PMID:23286029
Dual FIB-SEM 3D Imaging and Lattice Boltzmann Modeling of Porosimetry and Multiphase Flow in Chalk
NASA Astrophysics Data System (ADS)
Rinehart, A. J.; Yoon, H.; Dewers, T. A.; Heath, J. E.; Petrusak, R.
2010-12-01
Mercury intrusion porosimetry (MIP) is an often-applied technique for determining pore throat distributions and seal analysis of fine-grained rocks. Due to closure effects, potential pore collapse, and complex pore network topologies, MIP data interpretation can be ambiguous, and often biased toward smaller pores in the distribution. We apply 3D imaging techniques and lattice-Boltzmann modeling in interpreting MIP data for samples of the Cretaceous Selma Group Chalk. In the Mississippi Interior Salt Basin, the Selma Chalk is the apparent seal for oil and gas fields in the underlying Eutaw Fm., and, where unfractured, the Selma Chalk is one of the regional-scale seals identified by the Southeast Regional Carbon Sequestration Partnership for CO2 injection sites. Dual focused ion - scanning electron beam and laser scanning confocal microscopy methods are used for 3D imaging of nanometer-to-micron scale microcrack and pore distributions in the Selma Chalk. A combination of image analysis software is used to obtain geometric pore body and throat distributions and other topological properties, which are compared to MIP results. 3D data sets of pore-microfracture networks are used in Lattice Boltzmann simulations of drainage (wetting fluid displaced by non-wetting fluid via the Shan-Chen algorithm), which in turn are used to model MIP procedures. Results are used in interpreting MIP results, understanding microfracture-matrix interaction during multiphase flow, and seal analysis for underground CO2 storage. This work was supported by the US Department of Energy, Office of Basic Energy Sciences as part of an Energy Frontier Research Center. Sandia National Laboratories is a multi-program laboratory operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Company, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000.
Lattice Boltzmann Stokesian dynamics.
Ding, E J
2015-11-01
Lattice Boltzmann Stokesian dynamics (LBSD) is presented for simulation of particle suspension in Stokes flows. This method is developed from Stokesian dynamics (SD) with resistance and mobility matrices calculated using the time-independent lattice Boltzmann algorithm (TILBA). TILBA is distinguished from the traditional lattice Boltzmann method (LBM) in that a background matrix is generated prior to the calculation. The background matrix, once generated, can be reused for calculations for different scenarios, thus the computational cost for each such subsequent calculation is significantly reduced. The LBSD inherits the merits of the SD where both near- and far-field interactions are considered. It also inherits the merits of the LBM that the computational cost is almost independent of the particle shape. PMID:26651812
Lattice Boltzmann Stokesian dynamics
NASA Astrophysics Data System (ADS)
Ding, E. J.
2015-11-01
Lattice Boltzmann Stokesian dynamics (LBSD) is presented for simulation of particle suspension in Stokes flows. This method is developed from Stokesian dynamics (SD) with resistance and mobility matrices calculated using the time-independent lattice Boltzmann algorithm (TILBA). TILBA is distinguished from the traditional lattice Boltzmann method (LBM) in that a background matrix is generated prior to the calculation. The background matrix, once generated, can be reused for calculations for different scenarios, thus the computational cost for each such subsequent calculation is significantly reduced. The LBSD inherits the merits of the SD where both near- and far-field interactions are considered. It also inherits the merits of the LBM that the computational cost is almost independent of the particle shape.
NASA Astrophysics Data System (ADS)
Kajzer, A.; Pozorski, J.; Szewc, K.
2014-08-01
In the paper we present Large-eddy simulation (LES) results of 3D Taylor- Green vortex obtained by the three different computational approaches: Smoothed Particle Hydrodynamics (SPH), Lattice Boltzmann Method (LBM) and Finite Volume Method (FVM). The Smagorinsky model was chosen as a subgrid-scale closure in LES for all considered methods and a selection of spatial resolutions have been investigated. The SPH and LBM computations have been carried out with the use of the in-house codes executed on GPU and compared, for validation purposes, with the FVM results obtained using the open-source CFD software OpenFOAM. A comparative study in terms of one-point statistics and turbulent energy spectra shows a good agreement of LES results for all methods. An analysis of the GPU code efficiency and implementation difficulties has been made. It is shown that both SPH and LBM may offer a significant advantage over mesh-based CFD methods.
Lattice Boltzmann morphodynamic model
NASA Astrophysics Data System (ADS)
Zhou, Jian Guo
2014-08-01
Morphological change due to sediment transport is a common natural phenomenon in real flows. It involves complex processes of erosion and deposition such as those along beaches and in river beds, imposing a strong strain on human beings. Studying and understanding morphodynamic evolution are essential to protect living environment. Although there are conventional numerical methods like finite difference method and finite volume method for forecast of morphological change by solving flow and morphodynamic equations, the methods are too complex/inefficient to be applied to a real large scale problem. To overcome this, a lattice Boltzmann method is developed to simulate morphological evolution under flows. It provides an alternative way of studying morphodynamics at the full advantages of the lattice Boltzmann methodology. The model is verified by applications to the evolution of one and two dimensional sand dunes under shallow water flows.
NASA Astrophysics Data System (ADS)
Paradis, Hedvig; Andersson, Martin; Sundén, Bengt
2015-09-01
A 3D model at microscale by the lattice Boltzmann method (LBM) is proposed for part of an anode of a solid oxide fuel cell (SOFC) to analyze the interaction between the transport and reaction processes and structural parameters. The equations of charge, momentum, heat and mass transport are simulated in the model. The modeling geometry is created with randomly placed spheres to resemble the part of the anode structure close to the electrolyte. The electrochemical reaction processes are captured at specific sites where spheres representing Ni and YSZ materials are present with void space. This work focuses on analyzing the effect of structural parameters such as porosity, and percentage of active reaction sites on the ionic current density and concentration of H2 using LBM. It is shown that LBM can be used to simulate an SOFC anode at microscale and evaluate the effect of structural parameters on the transport processes to improve the performance of the SOFC anode. It was found that increasing the porosity from 30 to 50 % decreased the ionic current density due to a reduction in the number of reaction sites. Also the consumption of H2 decreased with increasing porosity. When the percentage of active reaction sites was increased while the porosity was kept constant, the ionic current density increased. However, the H2 concentration was slightly reduced when the percentage of active reaction sites was increased. The gas flow tortuosity decreased with increasing porosity.
NASA Astrophysics Data System (ADS)
Paradis, Hedvig; Andersson, Martin; Sundén, Bengt
2016-08-01
A 3D model at microscale by the lattice Boltzmann method (LBM) is proposed for part of an anode of a solid oxide fuel cell (SOFC) to analyze the interaction between the transport and reaction processes and structural parameters. The equations of charge, momentum, heat and mass transport are simulated in the model. The modeling geometry is created with randomly placed spheres to resemble the part of the anode structure close to the electrolyte. The electrochemical reaction processes are captured at specific sites where spheres representing Ni and YSZ materials are present with void space. This work focuses on analyzing the effect of structural parameters such as porosity, and percentage of active reaction sites on the ionic current density and concentration of H2 using LBM. It is shown that LBM can be used to simulate an SOFC anode at microscale and evaluate the effect of structural parameters on the transport processes to improve the performance of the SOFC anode. It was found that increasing the porosity from 30 to 50 % decreased the ionic current density due to a reduction in the number of reaction sites. Also the consumption of H2 decreased with increasing porosity. When the percentage of active reaction sites was increased while the porosity was kept constant, the ionic current density increased. However, the H2 concentration was slightly reduced when the percentage of active reaction sites was increased. The gas flow tortuosity decreased with increasing porosity.
Lattice gas and lattice Boltzmann computational physics
Chen, S.
1993-05-01
Recent developments of the lattice gas automata method and its extension to the lattice Boltzmann method have provided new computational schemes for solving a variety of partial differential equations and modeling different physics systems. The lattice gas method, regarded as the simplest microscopic and kinetic approach which generates meaningful macroscopic dynamics, is fully parallel and can be easily programmed on parallel machines. In this talk, the author will review basic principles of the lattice gas and lattice Boltzmann method, its mathematical foundation and its numerical implementation. A detailed comparison of the lattice Boltzmann method with the lattice gas technique and other traditional numerical schemes, including the finite-difference scheme and the pseudo-spectral method, for solving the Navier-Stokes hydrodynamic fluid flows, will be discussed. Recent achievements of the lattice gas and the the lattice Boltzmann method and their applications in surface phenomena, spinodal decomposition and pattern formation in chemical reaction-diffusion systems will be presented.
He, Xiaoyi; Lou, Li-Shi Lou, Li-Shi
1997-12-01
In this paper, the lattice Boltzmann equation is directly derived from the Boltzmann equation. It is shown that the lattice Boltzmann equation is a special discretized form of the Boltzmann equation. Various approximations for the discretization of the Boltzmann equation in both time and phase space are discussed in detail. A general procedure to derive the lattice Boltzmann model from the continuous Boltzmann equation is demonstrated explicitly. The lattice Boltzmann models derived include the two-dimensional 6-bit, 7-bit, and 9-bit, and three-dimensional 27-bit models. {copyright} {ital 1997} {ital The American Physical Society}
Crystallographic Lattice Boltzmann Method.
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-01-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows. PMID:27251098
Crystallographic Lattice Boltzmann Method
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-01-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows. PMID:27251098
Crystallographic Lattice Boltzmann Method
NASA Astrophysics Data System (ADS)
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-06-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows.
Fully relativistic lattice Boltzmann algorithm
Romatschke, P.; Mendoza, M.; Succi, S.
2011-09-15
Starting from the Maxwell-Juettner equilibrium distribution, we develop a relativistic lattice Boltzmann (LB) algorithm capable of handling ultrarelativistic systems with flat, but expanding, spacetimes. The algorithm is validated through simulations of a quark-gluon plasma, yielding excellent agreement with hydrodynamic simulations. The present scheme opens the possibility of transferring the recognized computational advantages of lattice kinetic theory to the context of both weakly and ultrarelativistic systems.
Lattice Boltzmann methods in Geosciences
NASA Astrophysics Data System (ADS)
Huber, Christian; Parmigiani, Andrea; Su, Yanqing
2014-05-01
Numerical models often offer the only possible approach to study the complex non-linear dynamics of geodynamical processes that are difficult or impossible to scale for laboratory experiments. The development of improved computer resources has allowed the emergence of large-scale parallel computations in Earth Sciences. These resources have lead to an increasing complexity in models where a greater number of adjustable parameters arise. Although the increasing number of free parameters offers a greater flexibility to fit satisfyingly the set of available constraints (e.g. geochemical, structural) it also provides new challenges in terms of the size of the parameter space and non-uniqueness of model solutions. Another significant challenge associated with state-of-the-art models is that their complexity is in general associated with the addition of parameterizations of the unresolved (small) scale processes. This trend calls for the development of complementary high-performance models to constrain the physics at small-scales where mass, momentum and energy exchanges at interfaces between different phases control the dynamics in heterogeneous media. We argue that more attention should be devoted to the development of multiphase numerical modeling at the granular (pore) scale to investigate the dynamical behavior of heterogeneous media and the emergence of feedbacks that influence the response of these media at much greater scales. The lattice Boltzmann method is a paradigm that emerged almost three decades ago. It is based on kinetic theory and follows a bottom-up approach that contrast the top-down strategy of standard methods such as Finite Volumes, FEM and Finite Differences. Lattice Boltzmann is ideally suited to handle the complex dynamics of multiphase systems at small spatial scales and is very efficient for parallel programing. In this presentation, we discuss the development of different lattice Boltzmann models developed in our group over the last years
Lattice Boltzmann equation for relativistic quantum mechanics.
Succi, Sauro
2002-03-15
Relativistic versions of the quantum lattice Boltzmann equation are discussed. It is shown that the inclusion of nonlinear interactions requires the standard collision operator to be replaced by a pair of dynamic fields coupling to the relativistic wave function in a way which can be described by a multicomponent complex lattice Boltzmann equation. PMID:16210189
Recent advances in lattice Boltzmann methods
Chen, S.; Doolen, G.D.; He, X.; Nie, X.; Zhang, R.
1998-12-31
In this paper, the authors briefly present the basic principles of lattice Boltzmann method and summarize recent advances of the method, including the application of the lattice Boltzmann method for fluid flows in MEMS and simulation of the multiphase mixing and turbulence.
Lattice Boltzmann modeling of phonon transport
NASA Astrophysics Data System (ADS)
Guo, Yangyu; Wang, Moran
2016-06-01
A novel lattice Boltzmann scheme is proposed for phonon transport based on the phonon Boltzmann equation. Through the Chapman-Enskog expansion, the phonon lattice Boltzmann equation under the gray relaxation time approximation recovers the classical Fourier's law in the diffusive limit. The numerical parameters in the lattice Boltzmann model are therefore rigorously correlated to the bulk material properties. The new scheme does not only eliminate the fictitious phonon speed in the diagonal direction of a square lattice system in the previous lattice Boltzmann models, but also displays very robust performances in predicting both temperature and heat flux distributions consistent with analytical solutions for diverse numerical cases, including steady-state and transient, macroscale and microscale, one-dimensional and multi-dimensional phonon heat transport. This method may provide a powerful numerical tool for deep studies of nonlinear and nonlocal heat transports in nanosystems.
A Lattice Boltzmann Method for Turbomachinery Simulations
NASA Technical Reports Server (NTRS)
Hsu, A. T.; Lopez, I.
2003-01-01
Lattice Boltzmann (LB) Method is a relatively new method for flow simulations. The start point of LB method is statistic mechanics and Boltzmann equation. The LB method tries to set up its model at molecular scale and simulate the flow at macroscopic scale. LBM has been applied to mostly incompressible flows and simple geometry.
Lattice Boltzmann approach for complex nonequilibrium flows.
Montessori, A; Prestininzi, P; La Rocca, M; Succi, S
2015-10-01
We present a lattice Boltzmann realization of Grad's extended hydrodynamic approach to nonequilibrium flows. This is achieved by using higher-order isotropic lattices coupled with a higher-order regularization procedure. The method is assessed for flow across parallel plates and three-dimensional flows in porous media, showing excellent agreement of the mass flow with analytical and numerical solutions of the Boltzmann equation across the full range of Knudsen numbers, from the hydrodynamic regime to ballistic motion. PMID:26565365
Higher Order Thermal Lattice Boltzmann Model
NASA Astrophysics Data System (ADS)
Sorathiya, Shahajhan; Ansumali, Santosh
2013-03-01
Lattice Boltzmann method (LBM) modelling of thermal flows, compressible and micro flows requires an accurate velocity space discretization. The sub optimality of Gauss-Hermite quadrature in this regard is well known. Most of the thermal LBM in the past have suffered from instability due to lack of proper H-theorem and accuracy. Motivated from these issues, the present work develops along the two works and and imposes an eighth higher order moment to get correct thermal physics. We show that this can be done by adding just 6 more velocities to D3Q27 model and obtain a ``multi-speed on lattice thermal LBM'' with 33 velocities in 3D and calO (u4) and calO (T4) accurate fieq with a consistent H-theorem and inherent numerical stability. Simulations for Rayleigh-Bernard as well as velocity and temperature slip in micro flows matches with analytical results. Lid driven cavity set up for grid convergence is studied. Finally, a novel data structure is developed for HPC. The authors express their gratitude for computational resources and financial support provide by Jawaharlal Nehru Centre for Advanced Scientific Research (JNCASR), Bangalore, India.
Entropic lattice Boltzmann model for compressible flows.
Frapolli, N; Chikatamarla, S S; Karlin, I V
2015-12-01
We present a lattice Boltzmann model (LBM) that covers the entire range of fluid flows, from low Mach weakly compressible to transonic and supersonic flows. One of the most restrictive limitations of the lattice Boltzmann method, the low Mach number limit, is overcome here by three fundamental changes to the LBM scheme: use of an appropriately chosen multispeed lattice, accurate evaluation of the equilibrium, and the entropic relaxation for the collision. The range of applications is demonstrated through the simulation of a bow shock in front of an airfoil and the simulation of decaying compressible turbulence with shocklets. PMID:26764625
Entropic lattice Boltzmann model for compressible flows
NASA Astrophysics Data System (ADS)
Frapolli, N.; Chikatamarla, S. S.; Karlin, I. V.
2015-12-01
We present a lattice Boltzmann model (LBM) that covers the entire range of fluid flows, from low Mach weakly compressible to transonic and supersonic flows. One of the most restrictive limitations of the lattice Boltzmann method, the low Mach number limit, is overcome here by three fundamental changes to the LBM scheme: use of an appropriately chosen multispeed lattice, accurate evaluation of the equilibrium, and the entropic relaxation for the collision. The range of applications is demonstrated through the simulation of a bow shock in front of an airfoil and the simulation of decaying compressible turbulence with shocklets.
Visualization of 3D optical lattices
NASA Astrophysics Data System (ADS)
Lee, Hoseong; Clemens, James
2016-05-01
We describe the visualization of 3D optical lattices based on Sisyphus cooling implemented with open source software. We plot the adiabatic light shift potentials found by diagonalizing the effective Hamiltonian for the light shift operator. Our program incorporates a variety of atomic ground state configurations with total angular momentum ranging from j = 1 / 2 to j = 4 and a variety of laser beam configurations including the two-beam lin ⊥ lin configuration, the four-beam umbrella configuration, and four beams propagating in two orthogonal planes. In addition to visualizing the lattice the program also evaluates lattice parameters such as the oscillation frequency for atoms trapped deep in the wells. The program is intended to help guide experimental implementations of optical lattices.
Lattice Boltzmann model for compressible fluids
NASA Technical Reports Server (NTRS)
Alexander, F. J.; Chen, H.; Chen, S.; Doolen, G. D.
1992-01-01
A lattice Boltzmann model is derived which simulates compressible fluids. By choosing the parameters of the equilibrium distribution appropriately, the sound speed (which may be set arbitrarily low), bulk viscosity, and kinematic viscosity can be selected. This model simulates compressible flows and can include shocks. With a proper rescaling and zero-sound speed, this model simulates Burgers's equation. The viscosity determined by a Chapman-Enskog expansion compares well with that measured form simulations. The exact solutions of Burgers's equation on the unit circle are compared to solutions of lattice Boltzmann model finding reasonable agreement.
Fast Lattice Boltzmann Solver for Relativistic Hydrodynamics
Mendoza, M.; Herrmann, H. J.; Boghosian, B. M.; Succi, S.
2010-07-02
A lattice Boltzmann formulation for relativistic fluids is presented and numerically validated through quantitative comparison with recent hydrodynamic simulations of relativistic fluids. In order to illustrate its capability to handle complex geometries, the scheme is also applied to the case of a three-dimensional relativistic shock wave, generated by a supernova explosion, impacting on a massive interstellar cloud. This formulation opens up the possibility of exporting the proven advantages of lattice Boltzmann methods, namely, computational efficiency and easy handling of complex geometries, to the context of (mildly) relativistic fluid dynamics at large, from quark-gluon plasmas up to supernovae with relativistic outflows.
Connection Between the Lattice Boltzmann Equation and the Beam Scheme
NASA Technical Reports Server (NTRS)
Xu, Kun; Luo, Li-Shi
1999-01-01
In this paper we analyze and compare the lattice Boltzmann equation with the beam scheme in details. We notice the similarity and differences between the lattice Boltzmann equation and the beam scheme. We show that the accuracy of the lattice Boltzmann equation is indeed second order in space. We discuss the advantages and limitations of lattice Boltzmann equation and the beam scheme. Based on our analysis, we propose an improved multi-dimensional beam scheme.
Boundary conditions for free interfaces with the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Bogner, Simon; Ammer, Regina; Rüde, Ulrich
2015-09-01
In this paper we analyze the boundary treatment of the lattice Boltzmann method (LBM) for simulating 3D flows with free surfaces. The widely used free surface boundary condition of Körner et al. [27] is shown to be first order accurate. The article presents a new free surface boundary scheme that is suitable for second order accurate simulations based on the LBM. The new method takes into account the free surface position and its orientation with respect to the computational lattice. Numerical experiments confirm the theoretical findings and illustrate the different behavior of the original method and the new method.
Physical symmetry and lattice symmetry in the lattice Boltzmann method
Cao, N.; Chen, S.; Jin, S.; Martinez, D.
1997-01-01
The lattice Boltzmann method (LBM) is regarded as a specific finite difference discretization for the kinetic equation of the discrete velocity distribution function. We argue that for finite sets of discrete velocity models, such as LBM, the physical symmetry is necessary for obtaining the correct macroscopic Navier-Stokes equations. In contrast, the lattice symmetry and the Lagrangian nature of the scheme, which is often used in the lattice gas automaton method and the existing lattice Boltzmann methods and directly associated with the property of particle dynamics, is not necessary for recovering the correct macroscopic dynamics. By relaxing the lattice symmetry constraint and introducing other numerical discretization, one can also obtain correct hydrodynamics. In addition, numerical simulations for applications, such as nonuniform meshes and thermohydrodynamics can be easily carried out and numerical stability can be ensured by the Courant-Friedricks-Lewey condition and using the semi-implicit collision scheme. {copyright} {ital 1997} {ital The American Physical Society}
Stability analysis of lattice Boltzmann methods
Sterling, J.D.; Chen, Shiyi
1996-01-01
The lattice Boltzmann equation describes the evolution of the velocity distribution function on a lattice in a manner that macroscopic fluid dynamical behavior is recovered. Although the equation is a derivative of lattice gas automata, it may be interpreted as a Lagrangian finite-difference method for the numerical simulation of the discrete-velocity Boltzmann equation that makes use of a BGK collision operator. As a result, it is not surprising that numericaI instability of lattice Boltzmann methods have been frequently encountered by researchers. We present an analysis of the stability of perturbations of the particle populations linearized about equilibrium values corresponding to a constant-density uniform mean flow. The linear stability depends on the following parameters: the distribution of the mass at a site between the different discrete speeds, the BGK relaxation time, the mean velocity, and the wave-number of the perturbations. This parameter space is too large to compute the complete stability characteristics. We report some stability results for a subset of the parameter space for a 7-velocity hexagonal lattice, a 9-velocity square lattice, and a 15-velocity cubic lattice. Results common to all three lattices are (1) the BGK relaxation time {tau} must be greater than 1/2 corresponding to positive shear viscosity, (2) there exists a maximum stable mean velocity for fixed values of theother parameters, and (3) as {tau} is increased from 1/2 the maximum stable velocity increases monotonically until some fixed velocity is reached which does not change for larger {tau}.
Structural stability of Lattice Boltzmann schemes
NASA Astrophysics Data System (ADS)
David, Claire; Sagaut, Pierre
2016-02-01
The goal of this work is to determine classes of traveling solitary wave solutions for Lattice Boltzmann schemes by means of a hyperbolic ansatz. It is shown that spurious solitary waves can occur in finite-difference solutions of nonlinear wave equation. The occurrence of such a spurious solitary wave, which exhibits a very long life time, results in a non-vanishing numerical error for arbitrary time in unbounded numerical domain. Such a behavior is referred here to have a structural instability of the scheme, since the space of solutions spanned by the numerical scheme encompasses types of solutions (solitary waves in the present case) that are not solutions of the original continuous equations. This paper extends our previous work about classical schemes to Lattice Boltzmann schemes (David and Sagaut 2011; 2009a,b; David et al. 2007).
Lattice Boltzmann model for wave propagation.
Zhang, Jianying; Yan, Guangwu; Shi, Xiubo
2009-08-01
A lattice Boltzmann model for two-dimensional wave equation is proposed by using the higher-order moment method. The higher-order moment method is based on the solution of a series of partial differential equations obtained by using multiscale technique and Chapman-Enskog expansion. In order to obtain the lattice Boltzmann model for the wave equation with higher-order accuracy of truncation errors, we removed the second-order dissipation term and the third-order dispersion term by employing the moments up to fourth order. The reversibility in time appears owing to the absence of the second-order dissipation term and the third-order dispersion term. As numerical examples, some classical examples, such as interference, diffraction, and wave passing through a convex lens, are simulated. The numerical results show that this model can be used to simulate wave propagation. PMID:19792280
Lattice Boltzmann model for simulation of magnetohydrodynamics
NASA Technical Reports Server (NTRS)
Chen, Shiyi; Chen, Hudong; Martinez, Daniel; Matthaeus, William
1991-01-01
A numerical method, based on a discrete Boltzmann equation, is presented for solving the equations of magnetohydrodynamics (MHD). The algorithm provides advantages similar to the cellular automaton method in that it is local and easily adapted to parallel computing environments. Because of much lower noise levels and less stringent requirements on lattice size, the method appears to be more competitive with traditional solution methods. Examples show that the model accurately reproduces both linear and nonlinear MHD phenomena.
Lattice Boltzmann method and channel flow
NASA Astrophysics Data System (ADS)
Stensholt, Sigvat; Mongstad Hope, Sigmund
2016-07-01
Lattice Boltzmann methods are presented at an introductory level with a focus on fairly simple simulations that can be used to test and illustrate the model’s capabilities. Two scenarios are presented. The first is a simple laminar flow in a straight channel driven by a pressure gradient (Poiseuille flow). The second is a more complex, including a wedge where Moffatt vortices may be induced if the wedge is deep enough. Simulations of the Poiseuille flow scenario accurately capture the theoretical velocity profile. The experiment shows the location of the fluid-wall boundary and the effects viscosity has on the velocity and convergence time. The numerical capabilities of the lattice Boltzmann model are tested further by simulating the more complex Moffatt vortex scenario. The method reproduces with high accuracy the theoretical predction that Moffat vortices will not form in a wedge if the vertex angle exceeds 146°. Practical issues limitations of the lattice Boltzmann method are discussed. In particular the accuracy of the bounce-back boundary condition is first order dependent on the grid resolution.
Consistent lattice Boltzmann equations for phase transitions.
Siebert, D N; Philippi, P C; Mattila, K K
2014-11-01
Unlike conventional computational fluid dynamics methods, the lattice Boltzmann method (LBM) describes the dynamic behavior of fluids in a mesoscopic scale based on discrete forms of kinetic equations. In this scale, complex macroscopic phenomena like the formation and collapse of interfaces can be naturally described as related to source terms incorporated into the kinetic equations. In this context, a novel athermal lattice Boltzmann scheme for the simulation of phase transition is proposed. The continuous kinetic model obtained from the Liouville equation using the mean-field interaction force approach is shown to be consistent with diffuse interface model using the Helmholtz free energy. Density profiles, interface thickness, and surface tension are analytically derived for a plane liquid-vapor interface. A discrete form of the kinetic equation is then obtained by applying the quadrature method based on prescribed abscissas together with a third-order scheme for the discretization of the streaming or advection term in the Boltzmann equation. Spatial derivatives in the source terms are approximated with high-order schemes. The numerical validation of the method is performed by measuring the speed of sound as well as by retrieving the coexistence curve and the interface density profiles. The appearance of spurious currents near the interface is investigated. The simulations are performed with the equations of state of Van der Waals, Redlich-Kwong, Redlich-Kwong-Soave, Peng-Robinson, and Carnahan-Starling. PMID:25493907
Lattice Boltzmann model for numerical relativity.
Ilseven, E; Mendoza, M
2016-02-01
In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems. PMID:26986435
Lattice Boltzmann model for numerical relativity
NASA Astrophysics Data System (ADS)
Ilseven, E.; Mendoza, M.
2016-02-01
In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
Multireflection boundary conditions for lattice Boltzmann models.
Ginzburg, Irina; d'Humières, Dominique
2003-12-01
We present a general framework for several previously introduced boundary conditions for lattice Boltzmann models, such as the bounce-back rule and the linear and quadratic interpolations. The objectives are twofold: first to give theoretical tools to study the existing link-type boundary conditions and their corresponding accuracy; second to design boundary conditions for general flows which are third-order kinetic accurate. Using these new boundary conditions, Couette and Poiseuille flows are exact solutions of the lattice Boltzmann models for a Reynolds number Re=0 (Stokes limit) for arbitrary inclination with the lattice directions. Numerical comparisons are given for Stokes flows in periodic arrays of spheres and cylinders, linear periodic array of cylinders between moving plates, and for Navier-Stokes flows in periodic arrays of cylinders for Re<200. These results show a significant improvement of the overall accuracy when using the linear interpolations instead of the bounce-back reflection (up to an order of magnitude on the hydrodynamics fields). Further improvement is achieved with the new multireflection boundary conditions, reaching a level of accuracy close to the quasianalytical reference solutions, even for rather modest grid resolutions and few points in the narrowest channels. More important, the pressure and velocity fields in the vicinity of the obstacles are much smoother with multireflection than with the other boundary conditions. Finally the good stability of these schemes is highlighted by some simulations of moving obstacles: a cylinder between flat walls and a sphere in a cylinder. PMID:14754343
Lattice Boltzmann approach to thermal transpiration
Sofonea, Victor
2006-11-15
Diffuse reflection boundary conditions are introduced in a thermal lattice Boltzmann model to allow for variable fluid density and temperature along the walls. The capability of this model to capture the main characteristics of the thermal transpiration phenomenon in a box at nonvanishing Knudsen numbers is demonstrated. The thermal creep velocity is found to be proportional to the temperature gradient imposed at the wall, whereas the accuracy of the simulation results are found to be of first or second order, depending on the numerical scheme.
Lattice-Boltzmann-based Simulations of Diffusiophoresis
NASA Astrophysics Data System (ADS)
Castigliego, Joshua; Kreft Pearce, Jennifer
We present results from a lattice-Boltzmann-base Brownian Dynamics simulation on diffusiophoresis and the separation of particles within the system. A gradient in viscosity that simulates a concentration gradient in a dissolved polymer allows us to separate various types of particles by their deformability. As seen in previous experiments, simulated particles that have a higher deformability react differently to the polymer matrix than those with a lower deformability. Therefore, the particles can be separated from each other. This simulation, in particular, was intended to model an oceanic system where the particles of interest were zooplankton, phytoplankton and microplastics. The separation of plankton from the microplastics was achieved.
Thermal lattice Boltzmann method for complex microflows
NASA Astrophysics Data System (ADS)
Yasuoka, Haruka; Kaneda, Masayuki; Suga, Kazuhiko
2016-07-01
A methodology to simulate thermal fields in complex microflow geometries is proposed. For the flow fields, the regularized multiple-relaxation-time lattice Boltzmann method (LBM) is applied coupled with the diffusive-bounce-back boundary condition for wall boundaries. For the thermal fields, the regularized lattice Bhatnagar-Gross-Krook model is applied. For the thermal wall boundary condition, a newly developed boundary condition, which is a mixture of the diffuse scattering and constant temperature conditions, is applied. The proposed set of schemes is validated by reference data in the Fourier flows and square cylinder flows confined in a microchannel. The obtained results confirm that it is essential to apply the regularization to the thermal LBM for avoiding kinked temperature profiles in complex thermal flows. The proposed wall boundary condition is successful to obtain thermal jumps at the walls with good accuracy.
Multi-Species Thermal Lattice Boltzmann Models
NASA Astrophysics Data System (ADS)
Wah, Darren; Vahala, George; Vahala, Linda; Pavlo, Pavol; Carter, Jonathan
1998-11-01
Thermal Lattice Boltzmann models (TLBM) are ideal for simulating nonlinear macroscopic conservation systems because of their inherent parallelizeability (nearly all operations are purely local). The TLBM solves a linear BGK-like kinetic equation so that the standard nonlinear convective terms in the standard fluid codes are now replaced by a simple shift operator (linear advection) at the kinetic level. Here we extend our previous TLBM to handle a two-species system, utilizing the models of Morse (1964),Greene (1973) and Kotelnikov & Montgomery (1997). Each kinetic equation now has 2 BGK-like relaxation terms : the first is due to self-collisions and the other is due to different- species collisions. The relaxation rates used are appropriate for electron-ion collisions. Certain constraints can be imposed on the relaxed distribution functions so that the cross-species momentum and energy evolutions relax at the rate determined from the full nonlinear Boltzmann integral collision operator. Ionization and recombination processes will also be examined. Both hexagonal and octagonal lattices are studied.
Immersed boundary method implemented in lattice Boltzmann GPU code
NASA Astrophysics Data System (ADS)
Devincentis, Brian; Smith, Kevin; Thomas, John
2015-11-01
Lattice Boltzmann is well suited to efficiently utilize the rapidly increasing compute power of GPUs to simulate viscous incompressible flows. Fluid-structure interaction with solids of arbitrarily complex geometry can be modeled in this framework with the immersed boundary method (IBM). In IBM a solid is modeled by its surface which applies a force at the neighboring lattice points. The majority of published IBMs require solving a linear system in order to satisfy the no-slip condition. However, the method presented by Wang et al. (2014) is unique in that it produces equally accurate results without solving a linear system. Furthermore, the algorithm can be applied in a parallel manner over the immersed boundary and is, therefore, well suited for GPUs. Here, a 2D and 3D version of their algorithm is implemented in Sailfish CFD, a GPU-based open source lattice Boltzmann code. One issue unaddressed by most published work is how to correct force and torque calculated from IBM for translating and rotating solids. These corrections are necessary because the fluid inside the solid affects its inertia in a non-trivial manner. Therefore, this implementation uses the Lagrangian points approximation correction shown by Suzuki and Inamuro (2011) to be accurate.
Thermal equation of state for lattice Boltzmann gases
NASA Astrophysics Data System (ADS)
Ran, Zheng
2009-06-01
The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.
Lattice-Boltzmann Simulation of Tablet Disintegration
NASA Astrophysics Data System (ADS)
Jiang, Jiaolong; Sun, Ning; Gersappe, Dilip
Using the lattice-Boltzmann method, we developed a 2D model to study the tablet disintegration involving the swelling and wicking mechanisms. The surface area and disintegration profile of each component were obtained by tracking the tablet structure in the simulation. Compared to pure wicking, the total surface area is larger for swelling and wicking, which indicates that the swelling force breaks the neighboring bonds. The disintegration profiles show that the tablet disintegrates faster than pure wicking, and there are more wetted active pharmaceutical ingredient particles distributed on smaller clusters. Our results indicate how the porosity would affect the disintegration process by changing the wetting area of the tablet as well as by changing the swelling force propagation.
Lattice Boltzmann simulations of lymphatic pumping
NASA Astrophysics Data System (ADS)
Kunert, Christian; Padera, Tim P.; Munn, Lance L.
2012-02-01
Lymphatic flow plays an important role in the progress of many diseases, including lymphedema and metastasis. However lymphatic pumping and flow is poorly understood. Here, we present a computer model that is based on biological observations of lymphatic pumping. Fluid flow is simulated by a D2Q9 lattice Boltzmann model. The boundary of the vessels moves according to shear-induced nitric oxide production, and wall motion transfers momentum to the fluid to induce flow. Because the model only includes local properties, it can be highly parallelized. In our case we utilize graphic processors (GPU) to achieve high performance computation. We show that the model provides stable pumping over a wide range of parameter values, with optimum flow achieved in the biological ranges. Furthermore, we investigate the efficiency by analyzing the flow rate and pumping frequency in order to compare the behavior of the model with existing in vivo data.
Lattice Boltzmann formulation for Braginskii magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Dellar, Paul
2012-03-01
We present a lattice Boltzmann formulation of the Braginskii magnetohydrodynamic equations that describe large-scale motions in strongly magnetised plasmas. Fluid quantities, density, velocity and stress, are represented by a finite set of distribution functions associated with particles moving on a square or cubic lattice. Equilibrium distributions are constructed from Hermite moment expansions, so slowly varying solutions of the discrete kinetic equation exactly satisfy the Navier--Stokes or MHD momentum equations. Electromagnetic quantities are represented by a second kinetic equation for a set of vector-valued distribution functions. Maxwell's equations and the resistive MHD induction equation may be recovered from slowly varying solutions using different scalings. The resulting algorithm, comprising only local operations at grid points and data copying between adjacent points, readily lends itself to large-scale parallel computations. We modify the collision operator to apply different relaxation times to components of the stress parallel and perpendicular to the local magnetic field, simulating a form of the Braginskii MHD equations encountered in astrophysics. Large shears develop in simulations where the fluid velocity perpendicular to the field lines reverses.
Modeling adsorption with lattice Boltzmann equation
Guo, Long; Xiao, Lizhi; Shan, Xiaowen; Zhang, Xiaoling
2016-01-01
The research of adsorption theory has recently gained renewed attention due to its critical relevance to a number of trending industrial applications, hydrogen storage and shale gas exploration for instance. The existing theoretical foundation, laid mostly in the early twentieth century, was largely based on simple heuristic molecular interaction models and static interaction potential which, although being insightful in illuminating the fundamental mechanisms, are insufficient for computations with realistic adsorbent structure and adsorbate hydrodynamics, both critical for real-life applications. Here we present and validate a novel lattice Boltzmann model incorporating both adsorbate-adsorbate and adsorbate-adsorbent interactions with hydrodynamics which, for the first time, allows adsorption to be computed with real-life details. Connection with the classic Ono-Kondo lattice theory is established and various adsorption isotherms, both within and beyond the IUPAC classification are observed as a pseudo-potential is varied. This new approach not only enables an important physical to be simulated for real-life applications, but also provides an enabling theoretical framework within which the fundamentals of adsorption can be studied. PMID:27256325
Modeling adsorption with lattice Boltzmann equation.
Guo, Long; Xiao, Lizhi; Shan, Xiaowen; Zhang, Xiaoling
2016-01-01
The research of adsorption theory has recently gained renewed attention due to its critical relevance to a number of trending industrial applications, hydrogen storage and shale gas exploration for instance. The existing theoretical foundation, laid mostly in the early twentieth century, was largely based on simple heuristic molecular interaction models and static interaction potential which, although being insightful in illuminating the fundamental mechanisms, are insufficient for computations with realistic adsorbent structure and adsorbate hydrodynamics, both critical for real-life applications. Here we present and validate a novel lattice Boltzmann model incorporating both adsorbate-adsorbate and adsorbate-adsorbent interactions with hydrodynamics which, for the first time, allows adsorption to be computed with real-life details. Connection with the classic Ono-Kondo lattice theory is established and various adsorption isotherms, both within and beyond the IUPAC classification are observed as a pseudo-potential is varied. This new approach not only enables an important physical to be simulated for real-life applications, but also provides an enabling theoretical framework within which the fundamentals of adsorption can be studied. PMID:27256325
Convection in Multiphase Fluid Flows Using Lattice Boltzmann Methods
NASA Astrophysics Data System (ADS)
Biferale, L.; Perlekar, P.; Sbragaglia, M.; Toschi, F.
2012-03-01
We present high-resolution numerical simulations of convection in multiphase flows (boiling) using a novel algorithm based on a lattice Boltzmann method. We first study the thermodynamical and kinematic properties of the algorithm. Then, we perform a series of 3D numerical simulations changing the mean properties in the phase diagram and compare convection with and without phase coexistence at Rayleigh number Ra˜107. We show that in the presence of nucleating bubbles non-Oberbeck-Boussinesq effects develop, the mean temperature profile becomes asymmetric, and heat-transfer and heat-transfer fluctuations are enhanced, at all Ra studied. We also show that small-scale properties of velocity and temperature fields are strongly affected by the presence of the buoyant bubble leading to high non-Gaussian profiles in the bulk.
Lattice Boltzmann model for resistive relativistic magnetohydrodynamics.
Mohseni, F; Mendoza, M; Succi, S; Herrmann, H J
2015-08-01
In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1/2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere. PMID:26382548
Lattice Boltzmann model for resistive relativistic magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Mohseni, F.; Mendoza, M.; Succi, S.; Herrmann, H. J.
2015-08-01
In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1 / 2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.
Lattice Boltzmann algorithm for continuum multicomponent flow.
Halliday, I; Hollis, A P; Care, C M
2007-08-01
We present a multicomponent lattice Boltzmann simulation for continuum fluid mechanics, paying particular attention to the component segregation part of the underlying algorithm. In the principal result of this paper, the dynamics of a component index, or phase field, is obtained for a segregation method after U. D'Ortona [Phys. Rev. E 51, 3718 (1995)], due to Latva-Kokko and Rothman [Phys. Rev. E 71 056702 (2005)]. The said dynamics accord with a simulation designed to address multicomponent flow in the continuum approximation and underwrite improved simulation performance in two main ways: (i) by reducing the interfacial microcurrent activity considerably and (ii) by facilitating simulational access to regimes of flow with a low capillary number and drop Reynolds number [I. Halliday, R. Law, C. M. Care, and A. Hollis, Phys. Rev. E 73, 056708 (2006)]. The component segregation method studied, used in conjunction with Lishchuk's method [S. V. Lishchuk, C. M. Care, and I. Halliday, Phys. Rev. E 67, 036701 (2003)], produces an interface, which is distributed in terms of its component index; however, the hydrodynamic boundary conditions which emerge are shown to support the notion of a sharp, unstructured, continuum interface. PMID:17930175
Lattice Boltzmann Simulations of Peristaltic Particle Transport
NASA Astrophysics Data System (ADS)
Connington, Kevin; Kang, Qinjun; Viswanathan, Hari; Chen, Shiyi; Abdel-Fattah, Amr
2008-11-01
A peristaltic flow occurs when a tube or channel with flexible walls transports the contained fluid by progressing a series of contraction or expansion waves along the length of those walls. It is a mechanism used to transport fluid and immersed solid particles when it is ineffective or impossible to impose a favorable pressure gradient or desirous to avoid contact between the transported mixture and mechanical moving parts. Peristaltic transport occurs in many physiological situations and has myriad industrial applications. We focus our study on the peristaltic transport of a macroscopic particle in a two dimensional channel using the Lattice Boltzmann Method (LBM). We systematically investigate the effect of variation of the relevant non-dimensional parameters of the system on the particle transport. We examine the particle behavior when the system exhibits the peculiar phenomenon of fluid ``trapping.'' Finally, we analyze how the particle presence affects stress, pressure, and dissipation in the fluid in hopes of determining preferred working conditions for peristaltic transport of shear-sensitive particles.
Three-dimensional lattice Boltzmann model for magnetic reconnection
Mendoza, M.; Munoz, J. D.
2008-02-15
We develop a three-dimensional (3D) lattice Boltzmann model that recovers in the continuous limit the two-fluids theory for plasmas, and consequently includes the generalized Ohm's law. The model reproduces the magnetic reconnection process just by giving the right initial equilibrium conditions in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, the plasma is handled similar to two fluids with an interaction term, each one with distribution functions associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered as a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model can simulate either viscous fluids in the incompressible limit or nonviscous compressible fluids, and successfully reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. The reconnection rate in the magnetotail obtained with this model lies between R=0.062 and R=0.073, in good agreement with the observations.
Three-dimensional lattice Boltzmann model for magnetic reconnection.
Mendoza, M; Muñoz, J D
2008-02-01
We develop a three-dimensional (3D) lattice Boltzmann model that recovers in the continuous limit the two-fluids theory for plasmas, and consequently includes the generalized Ohm's law. The model reproduces the magnetic reconnection process just by giving the right initial equilibrium conditions in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, the plasma is handled similar to two fluids with an interaction term, each one with distribution functions associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered as a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model can simulate either viscous fluids in the incompressible limit or nonviscous compressible fluids, and successfully reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. The reconnection rate in the magnetotail obtained with this model lies between R=0.062 and R=0.073, in good agreement with the observations. PMID:18352154
LUDWIG: A parallel Lattice-Boltzmann code for complex fluids
NASA Astrophysics Data System (ADS)
Desplat, Jean-Christophe; Pagonabarraga, Ignacio; Bladon, Peter
2001-03-01
This paper describes Ludwig, a versatile code for the simulation of Lattice-Boltzmann (LB) models in 3D on cubic lattices. In fact, Ludwig is not a single code, but a set of codes that share certain common routines, such as I/O and communications. If Ludwig is used as intended, a variety of complex fluid models with different equilibrium free energies are simple to code, so that the user may concentrate on the physics of the problem, rather than on parallel computing issues. Thus far, Ludwig's main application has been to symmetric binary fluid mixtures. We first explain the philosophy and structure of Ludwig which is argued to be a very effective way of developing large codes for academic consortia. Next we elaborate on some parallel implementation issues such as parallel I/O, and the use of MPI to achieve full portability and good efficiency on both MPP and SMP systems. Finally, we describe how to implement generic solid boundaries, and look in detail at the particular case of a symmetric binary fluid mixture near a solid wall. We present a novel scheme for the thermodynamically consistent simulation of wetting phenomena, in the presence of static and moving solid boundaries, and check its performance.
RNA folding on the 3D triangular lattice
2009-01-01
Background Difficult problems in structural bioinformatics are often studied in simple exact models to gain insights and to derive general principles. Protein folding, for example, has long been studied in the lattice model. Recently, researchers have also begun to apply the lattice model to the study of RNA folding. Results We present a novel method for predicting RNA secondary structures with pseudoknots: first simulate the folding dynamics of the RNA sequence on the 3D triangular lattice, next extract and select a set of disjoint base pairs from the best lattice conformation found by the folding simulation. Experiments on sequences from PseudoBase show that our prediction method outperforms the HotKnot algorithm of Ren, Rastegari, Condon and Hoos, a leading method for RNA pseudoknot prediction. Our method for RNA secondary structure prediction can be adapted into an efficient reconstruction method that, given an RNA sequence and an associated secondary structure, finds a conformation of the sequence on the 3D triangular lattice that realizes the base pairs in the secondary structure. We implemented a suite of computer programs for the simulation and visualization of RNA folding on the 3D triangular lattice. These programs come with detailed documentation and are accessible from the companion website of this paper at http://www.cs.usu.edu/~mjiang/rna/DeltaIS/. Conclusion Folding simulation on the 3D triangular lattice is effective method for RNA secondary structure prediction and lattice conformation reconstruction. The visualization software for the lattice conformations of RNA structures is a valuable tool for the study of RNA folding and is a great pedagogic device. PMID:19891777
An axisymmetric multiple-relaxation-time lattice Boltzmann scheme
NASA Astrophysics Data System (ADS)
Xie, Wenjun
2015-01-01
A multiple-relaxation-time (MRT) lattice Boltzmann (LB) scheme developed for axisymmetric flows recovers the complete continuity and Navier-Stokes equations. This scheme follows the strategy of the standard D2Q9 model by using a single particle distribution function and a simple "collision-streaming" updating rule. The extra terms related to axisymmetry in the macroscopic equations are recovered by adding source terms into the LB equation, which are simple and involve no gradients. The compressible effect retained in the Navier-Stokes equations is recovered by introducing a term related to the reversed transformation matrix for MRT collision operator, so as to produce a correct bulk viscosity, making it suitable for compressible flows with high frequency and low Mach number. The validity of the scheme is demonstrated by testing the Hagen-Poiseuille flow and 3D Womersley flow, as well as the standing acoustic waves in a closed cylindrical chamber. The numerical experiments show desirable stability at low viscosities, enabling to simulate a standing ultrasound field in centimeters space.
NASA Technical Reports Server (NTRS)
Mei, Ren-Wei; Shyy, Wei; Yu, Da-Zhi; Luo, Li-Shi; Rudy, David (Technical Monitor)
2001-01-01
The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is
The Lattice Boltzmann Method applied to neutron transport
Erasmus, B.; Van Heerden, F. A.
2013-07-01
In this paper the applicability of the Lattice Boltzmann Method to neutron transport is investigated. One of the main features of the Lattice Boltzmann method is the simultaneous discretization of the phase space of the problem, whereby particles are restricted to move on a lattice. An iterative solution of the operator form of the neutron transport equation is presented here, with the first collision source as the starting point of the iteration scheme. A full description of the discretization scheme is given, along with the quadrature set used for the angular discretization. An angular refinement scheme is introduced to increase the angular coverage of the problem phase space and to mitigate lattice ray effects. The method is applied to a model problem to investigate its applicability to neutron transport and the results are compared to a reference solution calculated, using MCNP. (authors)
The lattice Boltzmann method and the problem of turbulence
Djenidi, L.
2015-03-10
This paper reports a brief review of numerical simulations of homogeneous isotopic turbulence (HIT) using the lattice Boltzmann method (LBM). The LBM results shows that the details of HIT are well captured and in agreement with existing data. This clearly indicates that the LBM is as good as current Navier-Stokes solvers and is very much adequate for investigating the problem of turbulence.
Lattice Boltzmann method for weakly ionized isothermal plasmas
Li Huayu; Ki, Hyungson
2007-12-15
In this paper, a lattice Boltzmann method (LBM) for weakly ionized isothermal plasmas is presented by introducing a rescaling scheme for the Boltzmann transport equation. Without using this rescaling, we found that the nondimensional relaxation time used in the LBM is too large and the LBM does not produce physically realistic results. The developed model was applied to the electrostatic wave problem and the diffusion process of singly ionized helium plasmas with a 1-3% degree of ionization under an electric field. The obtained results agree well with theoretical values.
Zalzale, M.; McDonald, P.J.
2012-12-15
The lattice Boltzmann method is used to investigate the permeability of microstructures of cement pastes generated using the numerical models CEMHYD3D (Bentz, 1997) and {mu}IC (Bishnoi and Scrivener, 2009). Results are reported as a function of paste water-to-cement ratio and degree of hydration. The permeability decreases with increasing hydration and decreasing water-to-cement ratio in agreement with experiment. However the permeability is larger than the experimental data recorded using beam bending methods (Vichit-Vadakan and Scherer, 2002). Notwithstanding, the lattice Boltzmann results compare favourably with alternate numerical methods of permeability calculation for cement model microstructures. In addition, we show early results for the liquid/vapour capillary adsorption and desorption isotherms in the same model {mu}IC structures. The broad features of the experimental capillary porosity isotherm are reproduced, although further work is required to adequately parameterise the model.
Entropic multirelaxation lattice Boltzmann models for turbulent flows.
Bösch, Fabian; Chikatamarla, Shyam S; Karlin, Ilya V
2015-10-01
We present three-dimensional realizations of a class of lattice Boltzmann models introduced recently by the authors [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014)] and review the role of the entropic stabilizer. Both coarse- and fine-grid simulations are addressed for the Kida vortex flow benchmark. We show that the outstanding numerical stability and performance is independent of a particular choice of the moment representation for high-Reynolds-number flows. We report accurate results for low-order moments for homogeneous isotropic decaying turbulence and second-order grid convergence for most assessed statistical quantities. It is demonstrated that all the three-dimensional lattice Boltzmann realizations considered herein converge to the familiar lattice Bhatnagar-Gross-Krook model when the resolution is increased. Moreover, thanks to the dynamic nature of the entropic stabilizer, the present model features less compressibility effects and maintains correct energy and enstrophy dissipation. The explicit and efficient nature of the present lattice Boltzmann method renders it a promising candidate for both engineering and scientific purposes for highly turbulent flows. PMID:26565366
Entropic multirelaxation lattice Boltzmann models for turbulent flows
NASA Astrophysics Data System (ADS)
Bösch, Fabian; Chikatamarla, Shyam S.; Karlin, Ilya V.
2015-10-01
We present three-dimensional realizations of a class of lattice Boltzmann models introduced recently by the authors [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014), 10.1103/PhysRevE.90.031302] and review the role of the entropic stabilizer. Both coarse- and fine-grid simulations are addressed for the Kida vortex flow benchmark. We show that the outstanding numerical stability and performance is independent of a particular choice of the moment representation for high-Reynolds-number flows. We report accurate results for low-order moments for homogeneous isotropic decaying turbulence and second-order grid convergence for most assessed statistical quantities. It is demonstrated that all the three-dimensional lattice Boltzmann realizations considered herein converge to the familiar lattice Bhatnagar-Gross-Krook model when the resolution is increased. Moreover, thanks to the dynamic nature of the entropic stabilizer, the present model features less compressibility effects and maintains correct energy and enstrophy dissipation. The explicit and efficient nature of the present lattice Boltzmann method renders it a promising candidate for both engineering and scientific purposes for highly turbulent flows.
Lattice Boltzmann Model for Electronic Structure Simulations
NASA Astrophysics Data System (ADS)
Mendoza, M.; Herrmann, H. J.; Succi, S.
2015-09-01
Recently, a new connection between density functional theory and kinetic theory has been proposed. In particular, it was shown that the Kohn-Sham (KS) equations can be reformulated as a macroscopic limit of the steady-state solution of a suitable single-particle kinetic equation. By using a discrete version of this new formalism, the exchange and correlation energies of simple atoms and the geometrical configuration of the methane molecule were calculated accurately. Here, we discuss the main ideas behind the lattice kinetic approach to electronic structure computations, offer some considerations for prospective extensions, and also show additional numerical results, namely the geometrical configuration of the water molecule.
Lattice Boltzmann method for the Saint-Venant equations
NASA Astrophysics Data System (ADS)
Liu, Haifei; Wang, Hongda; Liu, Shu; Hu, Changwei; Ding, Yu; Zhang, Jie
2015-05-01
The Saint-Venant equations represent the hydrodynamic principles of unsteady flows in open channel network through a set of non-linear partial differential equations. In this paper, a new lattice Boltzmann approach to solving the one-dimensional Saint-Venant equations (LABSVE) is developed, demonstrating the variation of discharge and sectional area with external forces, such as bed slope and bed friction. Our research recovers the Saint-Venant equations through deducing the Chapman-Enskog expansion on the lattice Boltzmann equation, which is a mesoscopic technique, bridging the molecular movement and macroscopic physical variables. It is also a fully explicit process, providing simplicity for programming. The model is verified by three benchmark tests: (i) a one-dimensional subcritical gradient flow; (ii) a dam-break wave flow; (iii) a flood event on the Yongding River. The results showed the accuracy of the proposed method and its good applicability in solving Saint-Venant problems.
Conjugate heat transfer with the entropic lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Pareschi, G.; Frapolli, N.; Chikatamarla, S. S.; Karlin, I. V.
2016-07-01
A conjugate heat-transfer model is presented based on the two-population entropic lattice Boltzmann method. The present approach relies on the extension of Grad's boundary conditions to the two-population model for thermal flows, as well as on the appropriate exact conjugate heat-transfer condition imposed at the fluid-solid interface. The simplicity and efficiency of the lattice Boltzmann method (LBM), and in particular of the entropic multirelaxation LBM, are retained in the present approach, thus enabling simulations of turbulent high Reynolds number flows and complex wall boundaries. The model is validated by means of two-dimensional parametric studies of various setups, including pure solid conduction, conjugate heat transfer with a backward-facing step flow, and conjugate heat transfer with the flow past a circular heated cylinder. Further validations are performed in three dimensions for the case of a turbulent flow around a heated mounted cube.
Conjugate heat transfer with the entropic lattice Boltzmann method.
Pareschi, G; Frapolli, N; Chikatamarla, S S; Karlin, I V
2016-07-01
A conjugate heat-transfer model is presented based on the two-population entropic lattice Boltzmann method. The present approach relies on the extension of Grad's boundary conditions to the two-population model for thermal flows, as well as on the appropriate exact conjugate heat-transfer condition imposed at the fluid-solid interface. The simplicity and efficiency of the lattice Boltzmann method (LBM), and in particular of the entropic multirelaxation LBM, are retained in the present approach, thus enabling simulations of turbulent high Reynolds number flows and complex wall boundaries. The model is validated by means of two-dimensional parametric studies of various setups, including pure solid conduction, conjugate heat transfer with a backward-facing step flow, and conjugate heat transfer with the flow past a circular heated cylinder. Further validations are performed in three dimensions for the case of a turbulent flow around a heated mounted cube. PMID:27575234
Polar Coordinate Lattice Boltzmann Kinetic Modeling of Detonation Phenomena
NASA Astrophysics Data System (ADS)
Lin, Chuan-Dong; Xu, Ai-Guo; Zhang, Guang-Cai; Li, Ying-Jun
2014-11-01
A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier—Stokes equations, including contribution of chemical reaction, via the Chapman—Enskog expansion. For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.
Modified lattice Boltzmann method for compressible fluid simulations.
Hinton, F L; Rosenbluth, M N; Wong, S K; Lin-Liu, Y R; Miller, R L
2001-06-01
A modified lattice Boltzmann algorithm is shown to have much better stability to growing temperature perturbations, when compared with the standard lattice Boltzmann algorithm. The damping rates of long-wavelength waves, which determine stability, are derived using a collisional equilibrium distribution function which has the property that the Euler equations are obtained exactly in the limit of zero time step. Using this equilibrium distribution function, we show that our algorithm has inherent positive hyperviscosity and hyperdiffusivity, for very small values of viscosity and thermal diffusivity, which are lacking in the standard algorithm. Short-wavelength modes are shown to be stable for temperatures greater than a lower limit. Results from a computer code are used to compare these algorithms, and to confirm the damping rate predictions made analytically. Finite amplitude sound waves in the simulated fluid steepen, as expected from gas dynamic theory. PMID:11415085
Lattice Boltzmann method on unstructured grids: further developments.
Ubertini, S; Bella, G; Succi, S
2003-07-01
We discuss further developments of the finite-volume lattice Boltzmann formulation on unstructured grids. It is shown that the method tolerates significant grid distortions without showing any appreciable numerical viscosity effects at second order in the mesh size. A theoretical argument of plausibility for such a property is presented. In addition, a set of boundary conditions which permit to handle flows with open boundaries is also introduced and numerically demonstrated for the case of channel flows and driven cavity flows. PMID:12935281
Deviations from Boltzmann-Gibbs Statistics in Confined Optical Lattices.
Dechant, Andreas; Kessler, David A; Barkai, Eli
2015-10-23
We investigate the semiclassical phase-space probability distribution P(x,p) of cold atoms in a Sisyphus cooling lattice with an additional harmonic confinement. We pose the question of whether this nonequilibrium steady state satisfies the equivalence of energy and probability. This equivalence is the foundation of Boltzmann-Gibbs and generalized thermostatic statistics, and a prerequisite for the description in terms of a temperature. At large energies, P(x,p) depends only on the Hamiltonian H(x,p) and the answer to the question is yes. In distinction to the Boltzmann-Gibbs state, the large-energy tails are power laws P(x,p)∝H(x,p)(-1/D), where D is related to the depth of the optical lattice. At intermediate energies, however, P(x,p) cannot be expressed as a function of the Hamiltonian and the equivalence between energy and probability breaks down. As a consequence the average potential and kinetic energy differ and no well-defined temperature can be assigned. The Boltzmann-Gibbs state is regained only in the limit of deep optical lattices. For strong confinement relative to the damping, we derive an explicit expression for the stationary phase-space distribution. PMID:26551114
Zhang, Yimeng; Li, Xiong; Samonds, Jason M; Lee, Tai Sing
2016-03-01
Bayesian theory has provided a compelling conceptualization for perceptual inference in the brain. Central to Bayesian inference is the notion of statistical priors. To understand the neural mechanisms of Bayesian inference, we need to understand the neural representation of statistical regularities in the natural environment. In this paper, we investigated empirically how statistical regularities in natural 3D scenes are represented in the functional connectivity of disparity-tuned neurons in the primary visual cortex of primates. We applied a Boltzmann machine model to learn from 3D natural scenes, and found that the units in the model exhibited cooperative and competitive interactions, forming a "disparity association field", analogous to the contour association field. The cooperative and competitive interactions in the disparity association field are consistent with constraints of computational models for stereo matching. In addition, we simulated neurophysiological experiments on the model, and found the results to be consistent with neurophysiological data in terms of the functional connectivity measurements between disparity-tuned neurons in the macaque primary visual cortex. These findings demonstrate that there is a relationship between the functional connectivity observed in the visual cortex and the statistics of natural scenes. They also suggest that the Boltzmann machine can be a viable model for conceptualizing computations in the visual cortex and, as such, can be used to predict neural circuits in the visual cortex from natural scene statistics. PMID:26712581
Three-dimensional lattice Boltzmann model for compressible flows.
Sun, Chenghai; Hsu, Andrew T
2003-07-01
A three-dimensional compressible lattice Boltzmann model is formulated on a cubic lattice. A very large particle-velocity set is incorporated in order to enable a greater variation in the mean velocity. Meanwhile, the support set of the equilibrium distribution has only six directions. Therefore, this model can efficiently handle flows over a wide range of Mach numbers and capture shock waves. Due to the simple form of the equilibrium distribution, the fourth-order velocity tensors are not involved in the formulation. Unlike the standard lattice Boltzmann model, no special treatment is required for the homogeneity of fourth-order velocity tensors on square lattices. The Navier-Stokes equations were recovered, using the Chapman-Enskog method from the Bhatnagar-Gross-Krook (BGK) lattice Boltzmann equation. The second-order discretization error of the fluctuation velocity in the macroscopic conservation equation was eliminated by means of a modified collision invariant. The model is suitable for both viscous and inviscid compressible flows with or without shocks. Since the present scheme deals only with the equilibrium distribution that depends only on fluid density, velocity, and internal energy, boundary conditions on curved wall are easily implemented by an extrapolation of macroscopic variables. To verify the scheme for inviscid flows, we have successfully simulated a three-dimensional shock-wave propagation in a box and a normal shock of Mach number 10 over a wedge. As an application to viscous flows, we have simulated a flat plate boundary layer flow, flow over a cylinder, and a transonic flow over a NACA0012 airfoil cascade. PMID:12935242
NASA Astrophysics Data System (ADS)
Pingen, Georg
The objective of this work is the development of a formal design approach for fluidic systems, providing conceptually novel design layouts with the provision of only boundary conditions and some basic parameters. The lattice Boltzmann method (LBM) is chosen as a flow model due to its simplicity, inherent use of immersed boundary methods, parallelizability, and general flexibility. Immersed Boundary Methods in the form of a Brinkmann penalization are used to continuously vary the flow from fluid to solid, leading to a material distribution based boundary representation. An analytical adjoint sensitivity analysis is derived for the lattice Boltzmann method, enabling the combination of the lattice Boltzmann method with optimization techniques. This results in the first application of design optimization with the lattice Boltzmann method. In particular, the first LBM topology optimization framework for 2D and 3D problems is developed and validated with numerical design optimization problems for drag and pressure drop minimization. To improve the parallel scalability of the LBM sensitivity analysis and permit the solution of large 2D and 3D problems, iterative solvers are studied and a parallel GMRES Schur Complement method is applied to the solution of the linear adjoint problem in the LBM sensitivity analysis. This leads to improved parallel scalability through reduced memory use and algorithmic speedup. The potential of the developed design approach for fluidic systems is illustrated with the optimization of a 3D dual-objective fixed-geometry valve. The use of a parametric level-set method coupled with the LBM material distribution based topology optimization framework is shown to provide further versatility for design applications. Finally, the use of a penalty formulation of the fluid volume constraint permits the topology optimization of flows at moderate Reynolds numbers for a steady-state pipe bend application. Concluding, this work has led to the development of
McHardy, Christopher; Horneber, Tobias; Rauh, Cornelia
2016-07-25
Based on the kinetic theory of photons, a new lattice Boltzmann method for the simulation of 3D radiation transport is presented. The method was successfully validated with Monte Carlo simulations of radiation transport in optical thick absorbing and non-absorbing turbid media containing either isotropic or anisotropic scatterers. Moreover, for the approximation of Mie-scattering, a new iterative algebraic approach for the discretization of the scattering phase function was developed, ensuring full conservation of energy and asymmetry after discretization. It was found that the main error sources of the method are caused by linearization and ray effects and suggestions for further improvement of the method are made. PMID:27464152
Lattice Boltzmann Equation On a 2D Rectangular Grid
NASA Technical Reports Server (NTRS)
Bouzidi, MHamed; DHumieres, Dominique; Lallemand, Pierre; Luo, Li-Shi; Bushnell, Dennis M. (Technical Monitor)
2002-01-01
We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitrasy inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.
Lattice Boltzmann method for the fractional advection-diffusion equation
NASA Astrophysics Data System (ADS)
Zhou, J. G.; Haygarth, P. M.; Withers, P. J. A.; Macleod, C. J. A.; Falloon, P. D.; Beven, K. J.; Ockenden, M. C.; Forber, K. J.; Hollaway, M. J.; Evans, R.; Collins, A. L.; Hiscock, K. M.; Wearing, C.; Kahana, R.; Villamizar Velez, M. L.
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β , the fractional order α , and the single relaxation time τ , the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering.
Lattice Boltzmann simulation of chemical dissolution in porous media.
Kang, Qinjun; Zhang, Dongxiao; Chen, Shiyi; He, Xiaoyi
2002-03-01
In this paper, we develop a lattice Boltzmann model for simulating the transport and reaction of fluids in porous media. To simulate such a system, we account for the interaction of forced convection, molecular diffusion, and surface reaction. The problem is complicated by the evolution of the porous media geometry due to chemical reactions, which may significantly and continuously modify the hydrologic properties of the media. The particular application that motivates the present study is acid stimulation, a common technique used to increase production from petroleum reservoirs. This technique involves the injection of acid (e.g., hydrochloric acid, HCl, acetic acid, HAc) into the formation to dissolve minerals comprising the rock. As acid is injected, highly conductive channels or "wormholes" may be formed. The dissolution of carbonate rocks in 0.5M HCl and 0.5M HAc is simulated with the lattice Boltzmann model developed in this study. The dependence of dissolution process and the geometry of the final wormhole pattern on the acid type and the injection rate is studied. The results agree qualitatively with the experimental and theoretical analyses of others and substantiate the previous finding that there exists an optimal injection rate at which the wormhole is formed as well as the number of pore volumes of the injected fluid to break through is minimized. This study also confirms the experimentally observed phenomenon that the optimal injection rate decreases and the corresponding minimized number of pore volumes to break through increases as the acid is changed from HCl to HAc. Simulations suggest that the proposed lattice Boltzmann model may serve as an alternative reliable quantitative approach to study chemical dissolution in porous media. PMID:11909255
Lattice Boltzmann method for the fractional advection-diffusion equation.
Zhou, J G; Haygarth, P M; Withers, P J A; Macleod, C J A; Falloon, P D; Beven, K J; Ockenden, M C; Forber, K J; Hollaway, M J; Evans, R; Collins, A L; Hiscock, K M; Wearing, C; Kahana, R; Villamizar Velez, M L
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β, the fractional order α, and the single relaxation time τ, the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering. PMID:27176431
Accuracy of non-Newtonian Lattice Boltzmann simulations
NASA Astrophysics Data System (ADS)
Conrad, Daniel; Schneider, Andreas; Böhle, Martin
2015-11-01
This work deals with the accuracy of non-Newtonian Lattice Boltzmann simulations. Previous work for Newtonian fluids indicate that, depending on the numerical value of the dimensionless collision frequency Ω, additional artificial viscosity is introduced, which negatively influences the accuracy. Since the non-Newtonian fluid behavior is incorporated through appropriate modeling of the dimensionless collision frequency, a Ω dependent error EΩ is introduced and its influence on the overall error is investigated. Here, simulations with the SRT and the MRT model are carried out for power-law fluids in order to numerically investigate the accuracy of non-Newtonian Lattice Boltzmann simulations. A goal of this accuracy analysis is to derive a recommendation for an optimal choice of the time step size and the simulation Mach number, respectively. For the non-Newtonian case, an error estimate for EΩ in the form of a functional is derived on the basis of a series expansion of the Lattice Boltzmann equation. This functional can be solved analytically for the case of the Hagen-Poiseuille channel flow of non-Newtonian fluids. With the help of the error functional, the prediction of the global error minimum of the velocity field is excellent in regions where the EΩ error is the dominant source of error. With an optimal simulation Mach number, the simulation is about one order of magnitude more accurate. Additionally, for both collision models a detailed study of the convergence behavior of the method in the non-Newtonian case is conducted. The results show that the simulation Mach number has a major impact on the convergence rate and second order accuracy is not preserved for every choice of the simulation Mach number.
A Unified Theory of Non-Ideal Gas Lattice Boltzmann Models
NASA Technical Reports Server (NTRS)
Luo, Li-Shi
1998-01-01
A non-ideal gas lattice Boltzmann model is directly derived, in an a priori fashion, from the Enskog equation for dense gases. The model is rigorously obtained by a systematic procedure to discretize the Enskog equation (in the presence of an external force) in both phase space and time. The lattice Boltzmann model derived here is thermodynamically consistent and is free of the defects which exist in previous lattice Boltzmann models for non-ideal gases. The existing lattice Boltzmann models for non-ideal gases are analyzed and compared with the model derived here.
Lattice Boltzmann Modeling of Thrombosis in Giant Aneurysms
NASA Astrophysics Data System (ADS)
Chopard, B.; Ouared, R.; Ruefenacht, D. A.; Yilmaz, H.
We propose a numerical model of blood flow and blood clotting whose purpose is to describe thrombus formation in cerebral aneurysms. We identify possible mechanisms that can cause occurence of spontaneous thrombosis in unruptured giant intracranial aneurysms. Our main claim is that, under normal conditions, there is a low shear rate threshold below which thrombosis starts and growths. This assumption is supported by several evidences from literature. The proposed mechanisms are incorporated into a Lattice Boltzmann (LB) model for blood flow and platelets adhesion and aggregation. Numerical simulations show that the low shear rate threshold assumption together with aneurysm geometry account well for the observations.
Thrombosis modeling in intracranial aneurysms: a lattice Boltzmann numerical algorithm
NASA Astrophysics Data System (ADS)
Ouared, R.; Chopard, B.; Stahl, B.; Rüfenacht, D. A.; Yilmaz, H.; Courbebaisse, G.
2008-07-01
The lattice Boltzmann numerical method is applied to model blood flow (plasma and platelets) and clotting in intracranial aneurysms at a mesoscopic level. The dynamics of blood clotting (thrombosis) is governed by mechanical variations of shear stress near wall that influence platelets-wall interactions. Thrombosis starts and grows below a shear rate threshold, and stops above it. Within this assumption, it is possible to account qualitatively well for partial, full or no occlusion of the aneurysm, and to explain why spontaneous thrombosis is more likely to occur in giant aneurysms than in small or medium sized aneurysms.
Derivation of the lattice Boltzmann model for relativistic hydrodynamics
Mendoza, M.; Herrmann, H. J.; Boghosian, B. M.; Succi, S.
2010-11-15
A detailed derivation of the lattice Boltzmann scheme for relativistic fluids recently proposed in M. Mendoza, B. Boghosian, H. Herrmann, and S. Succi, Phys. Rev. Lett. 105, 014502 (2010) is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid-dynamic problems, namely, shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast wave on massive interstellar clouds. Close to second-order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time steps, are reported.
Thermodynamic-consistent lattice Boltzmann model for nonideal fluids
NASA Astrophysics Data System (ADS)
Wen, Binghai; Qin, Zhangrong; Zhang, Chaoying; Fang, Haiping
2015-11-01
A lattice Boltzmann model to simulate phase separation and two-phase flow is proposed. The nonideal force in multiphase flow is directly computed from the free energy. Thermodynamic consistency and Galilean invariance are theoretically analyzed and numerically verified. Remarkably, the theoretical simplicity endues the model with the advantages of high efficiency and easy implementation. We also find that it can work well together with various equations of state in order to simulate different kinds of multiphase flows. Importantly, it has a tunable parameter κ, which can be used to reduce the effect of spurious current and adjust the surface tension to meet the requirements of researches.
An implicit Lagrangian lattice Boltzmann method for the compressible flows
NASA Astrophysics Data System (ADS)
Yan, Guangwu; Dong, Yinfeng; Liu, Yanhong
2006-08-01
In this paper, we propose a new Lagrangian lattice Boltzmann method (LBM) for simulating the compressible flows. The new scheme simulates fluid flows based on the displacement distribution functions. The compressible flows, such as shock waves and contact discontinuities are modelled by using Lagrangian LBM. In this model, we select the element in the Lagrangian coordinate to satisfy the basic fluid laws. This model is a simpler version than the corresponding Eulerian coordinates, because the convection term of the Euler equations disappears. The numerical simulations conform to classical results.
Static contact angle in lattice Boltzmann models of immiscible fluids.
Latva-Kokko, M; Rothman, Daniel H
2005-10-01
We study numerically the capillary rise between two horizontal plates and in a rectangular tube, using a lattice Boltzmann (LB) method. We derive an equation for the static fluid-solid contact angle as a function of the wetting tendency of the walls and test its validity. We show that the generalized Laplace law with two independent radii of curvature is followed in capillary rise in rectangular tubes. Our method removes the history dependence of the fluid-solid contact angle that had been present in earlier LB schemes. PMID:16383561
Permeability of Partially Molten Rocks from Lattice-Boltzmann Modeling
NASA Astrophysics Data System (ADS)
Garapic, G.; Faul, U.
2013-12-01
Timescales of melt transport at mid-ocean ridges from mantle source to the surface depend on permeability of the partially molten mantle. The permeability is usually predicted indirectly from experimental observations based on porosities that are much higher than the porosities inferred for the partially molten mantle. Low porosities are for example predicted by geochemical models from the onset of melt migration. Since melting starts at the grain scale, permeability of the partially molten mantle will depend on the grain-scale melt distribution. We reconstructed a 3-D view of melt geometry of two experimentally produced samples of partially molten olivine which demonstrates that melt exists in thin layers on two-grain boundaries (Garapić et al.,G3, 2013). The wetted two-grain boundaries have a width about 100 times smaller than the average grain size. Additionally, the pore space consists of a network of triple-junction tubules at all porosities, and large 'melt pools'. Due to the relative size of the wetted two-grain boundaries as well as the size of the triple-junction network compared to the grain size imagining and numerical analyses of partially molten samples require high resolution. Since no direct experimental permeability measurements are possible on partially molten aggregates, we investigate numerically the permeability as a function of porosity for this system. We simulate porous flow through an artificial pore volume using the lattice-Boltzmann method (LBM) and Palabos LB code. Flow simulations were done on a computer cluster on three or four 125 GB nodes with 16 processors per node. With the available memory and allowed run time the maximum size of our pore structure was 1100 voxels per edge. In its simplest form the pore structure consists of a network of cylinders within a matrix of cubic grains. To approximate the observed 3-D melt geometry we added randomly distributed sheets on cube faces ('wetted two-grain boundaries') as well as randomly
Lattice Boltzmann model for a steady radiative transfer equation.
Yi, Hong-Liang; Yao, Feng-Ju; Tan, He-Ping
2016-08-01
A complete lattice Boltzmann model (LBM) is proposed for the steady radiative transfer equation (RTE). The RTE can be regarded as a pure convection equation with a source term. To derive the expressions for the equilibrium distribution function and the relaxation time, an artificial isotropic diffusion term is introduced to form a convection-diffusion equation. When the dimensionless relaxation time has a value of 0.5, the lattice Boltzmann equation (LBE) is exactly applicable to the original steady RTE. We also perform a multiscale analysis based on the Chapman-Enskog expansion to recover the macroscopic RTE from the mesoscopic LBE. The D2Q9 model is used to solve the LBE, and the numerical results obtained by the LBM are comparable to the results obtained by other methods or analytical solutions, which demonstrates that the proposed model is highly accurate and stable in simulating multidimensional radiative transfer. In addition, we find that the convergence rate of the LBM depends on the transport properties of RTE: for diffusion-dominated RTE with a large optical thickness, the LBM shows a second-order convergence rate in space, while for convection-dominated RTE with a small optical thickness, a lower convergence rate is observed. PMID:27627417
Lattice Boltzmann formulation for conjugate heat transfer in heterogeneous media
NASA Astrophysics Data System (ADS)
Karani, Hamid; Huber, Christian
2015-02-01
In this paper, we propose an approach for studying conjugate heat transfer using the lattice Boltzmann method (LBM). The approach is based on reformulating the lattice Boltzmann equation for solving the conservative form of the energy equation. This leads to the appearance of a source term, which introduces the jump conditions at the interface between two phases or components with different thermal properties. The proposed source term formulation conserves conductive and advective heat flux simultaneously, which makes it suitable for modeling conjugate heat transfer in general multiphase or multicomponent systems. The simple implementation of the source term approach avoids any correction of distribution functions neighboring the interface and provides an algorithm that is independent from the topology of the interface. Moreover, our approach is independent of the choice of lattice discretization and can be easily applied to different advection-diffusion LBM solvers. The model is tested against several benchmark problems including steady-state convection-diffusion within two fluid layers with parallel and normal interfaces with respect to the flow direction, unsteady conduction in a three-layer stratified domain, and steady conduction in a two-layer annulus. The LBM results are in excellent agreement with analytical solution. Error analysis shows that our model is first-order accurate in space, but an extension to a second-order scheme is straightforward. We apply our LBM model to heat transfer in a two-component heterogeneous medium with a random microstructure. This example highlights that the method we propose is independent of the topology of interfaces between the different phases and, as such, is ideally suited for complex natural heterogeneous media. We further validate the present LBM formulation with a study of natural convection in a porous enclosure. The results confirm the reliability of the model in simulating complex coupled fluid and thermal dynamics
Lattice Boltzmann Methods to Address Fundamental Boiling and Two-Phase Problems
Uddin, Rizwan
2012-01-01
This report presents the progress made during the fourth (no cost extension) year of this three-year grant aimed at the development of a consistent Lattice Boltzmann formulation for boiling and two-phase flows. During the first year, a consistent LBM formulation for the simulation of a two-phase water-steam system was developed. Results of initial model validation in a range of thermo-dynamic conditions typical for Boiling Water Reactors (BWRs) were shown. Progress was made on several fronts during the second year. Most important of these included the simulation of the coalescence of two bubbles including the surface tension effects. Work during the third year focused on the development of a new lattice Boltzmann model, called the artificial interface lattice Boltzmann model (AILB model) for the 3 simulation of two-phase dynamics. The model is based on the principle of free energy minimization and invokes the Gibbs-Duhem equation in the formulation of non-ideal forcing function. This was reported in detail in the last progress report. Part of the efforts during the last (no-cost extension) year were focused on developing a parallel capability for the 2D as well as for the 3D codes developed in this project. This will be reported in the final report. Here we report the work carried out on testing the AILB model for conditions including the thermal effects. A simplified thermal LB model, based on the thermal energy distribution approach, was developed. The simplifications are made after neglecting the viscous heat dissipation and the work done by pressure in the original thermal energy distribution model. Details of the model are presented here, followed by a discussion of the boundary conditions, and then results for some two-phase thermal problems.
NASA Astrophysics Data System (ADS)
Yu, Huidan; Chen, Xi; Wang, Zhiqiang; Deep, Debanjan; Lima, Everton; Zhao, Ye; Teague, Shawn D.
2014-06-01
In this paper, we develop a mass-conserved volumetric lattice Boltzmann method (MCVLBM) for numerically solving fluid dynamics with willfully moving arbitrary boundaries. In MCVLBM, fluid particles are uniformly distributed in lattice cells and the lattice Boltzmann equations deal with the time evolution of the particle distribution function. By introducing a volumetric parameter P (x,y,z,t) defined as the occupation of solid volume in the cell, we distinguish three types of lattice cells in the simulation domain: solid cell (pure solid occupation, P =1), fluid cell (pure fluid occupation, P =0), and boundary cell (partial solid and partial fluid, 0
lattice Boltzmann equations are self-regularized through P and consist of three parts: (1) collision taking into account the momentum exchange between the willfully moving boundary and the flow; (2) streaming accompanying a volumetric bounce-back procedure in boundary cells; and (3) boundary-induced volumetric fluid migration moving the residual fluid particles into the flow domain when the boundary swipes over a boundary cell toward a solid cell. The MCVLBM strictly satisfies mass conservation and can handle irregular boundary orientation and motion with respect to the mesh. Validation studies are carried out in four cases. The first is to simulate fluid dynamics in syringes focusing on how MCVLBM captures the underlying physics of flow driven by a willfully moving piston. The second and third cases are two-dimensional (2D) peristaltic flow and three-dimensional (3D) pipe flow, respectively. In each case, we compare the MCVLBM simulation result with the analytical solution and achieve quantitatively good agreements. The fourth case is to simulate blood flow in human aortic arteries with a very complicated irregular boundary. We study steady flow in two dimensions and unsteady flow via the pulsation of the cardiac cycle in three dimensions. In the 2D case, both vector (velocity) and
Large-scale lattice-Boltzmann simulations over lambda networks
NASA Astrophysics Data System (ADS)
Saksena, R.; Coveney, P. V.; Pinning, R.; Booth, S.
Amphiphilic molecules are of immense industrial importance, mainly due to their tendency to align at interfaces in a solution of immiscible species, e.g., oil and water, thereby reducing surface tension. Depending on the concentration of amphiphiles in the solution, they may assemble into a variety of morphologies, such as lamellae, micelles, sponge and cubic bicontinuous structures exhibiting non-trivial rheological properties. The main objective of this work is to study the rheological properties of very large, defect-containing gyroidal systems (of up to 10243 lattice sites) using the lattice-Boltzmann method. Memory requirements for the simulation of such large lattices exceed that available to us on most supercomputers and so we use MPICH-G2/MPIg to investigate geographically distributed domain decomposition simulations across HPCx in the UK and TeraGrid in the US. Use of MPICH-G2/MPIg requires the port-forwarder to work with the grid middleware on HPCx. Data from the simulations is streamed to a high performance visualisation resource at UCL (London) for rendering and visualisation. Lighting the Blue Touchpaper for UK e-Science - Closing Conference of ESLEA Project March 26-28 2007 The George Hotel, Edinburgh, UK
Chebyshev collocation spectral lattice Boltzmann method for simulation of low-speed flows.
Hejranfar, Kazem; Hajihassanpour, Mahya
2015-01-01
In this study, the Chebyshev collocation spectral lattice Boltzmann method (CCSLBM) is developed and assessed for the computation of low-speed flows. Both steady and unsteady flows are considered here. The discrete Boltzmann equation with the Bhatnagar-Gross-Krook approximation based on the pressure distribution function is considered and the space discretization is performed by the Chebyshev collocation spectral method to achieve a highly accurate flow solver. To provide accurate unsteady solutions, the time integration of the temporal term in the lattice Boltzmann equation is made by the fourth-order Runge-Kutta scheme. To achieve numerical stability and accuracy, physical boundary conditions based on the spectral solution of the governing equations implemented on the boundaries are used. An iterative procedure is applied to provide consistent initial conditions for the distribution function and the pressure field for the simulation of unsteady flows. The main advantage of using the CCSLBM over other high-order accurate lattice Boltzmann method (LBM)-based flow solvers is the decay of the error at exponential rather than at polynomial rates. Note also that the CCSLBM applied does not need any numerical dissipation or filtering for the solution to be stable, leading to highly accurate solutions. Three two-dimensional (2D) test cases are simulated herein that are a regularized cavity, the Taylor vortex problem, and doubly periodic shear layers. The results obtained for these test cases are thoroughly compared with the analytical and available numerical results and show excellent agreement. The computational efficiency of the proposed solution methodology based on the CCSLBM is also examined by comparison with those of the standard streaming-collision (classical) LBM and two finite-difference LBM solvers. The study indicates that the CCSLBM provides more accurate and efficient solutions than these LBM solvers in terms of CPU and memory usage and an exponential
Lattice Boltzmann simulations of a viscoelastic shear-thinning fluid
NASA Astrophysics Data System (ADS)
Papenkort, S.; Voigtmann, Th.
2015-07-01
We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell model that combines the effects of (linear) viscoelasticity and (nonlinear) shear thinning. We discuss the transient dynamics of velocities, shear stresses, and normal stress differences in planar pressure-driven channel flow, after switching on (startup) and off (cessation) of the driving pressure. This transient dynamics depends nontrivially on the channel width due to an interplay between hydrodynamic momentum diffusion and slow structural relaxation.
Lattice Boltzmann simulations of a viscoelastic shear-thinning fluid.
Papenkort, S; Voigtmann, Th
2015-07-28
We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell model that combines the effects of (linear) viscoelasticity and (nonlinear) shear thinning. We discuss the transient dynamics of velocities, shear stresses, and normal stress differences in planar pressure-driven channel flow, after switching on (startup) and off (cessation) of the driving pressure. This transient dynamics depends nontrivially on the channel width due to an interplay between hydrodynamic momentum diffusion and slow structural relaxation. PMID:26233150
Gauss Quadratures - the Keystone of Lattice Boltzmann Models
NASA Astrophysics Data System (ADS)
Piaud, Benjamin; Blanco, Stéphane; Fournier, Richard; Ambruş, Victor Eugen; Sofonea, Victor
2014-01-01
In this paper, we compare two families of Lattice Boltzmann (LB) models derived by means of Gauss quadratures in the momentum space. The first one is the HLB(N;Qx,Qy,Qz) family, derived by using the Cartesian coordinate system and the Gauss-Hermite quadrature. The second one is the SLB(N;K,L,M) family, derived by using the spherical coordinate system and the Gauss-Laguerre, as well as the Gauss-Legendre quadratures. These models order themselves according to the maximum order N of the moments of the equilibrium distribution function that are exactly recovered. Microfluidics effects (slip velocity, temperature jump, as well as the longitudinal heat flux that is not driven by a temperature gradient) are accurately captured during the simulation of Couette flow for Knudsen number (kn) up to 0.25.
Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation
NASA Technical Reports Server (NTRS)
Orr, Jeb S.; Powers, Joseph F.; Yang, Hong Q.
2015-01-01
A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments relies upon the global momentum conservation of the fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. Numerical examples illustrate the method's application to predicting bulk fluid motion including lateral propellant slosh in low-g conditions.
Lattice-Boltzmann hydrodynamics of anisotropic active matter.
de Graaf, Joost; Menke, Henri; Mathijssen, Arnold J T M; Fabritius, Marc; Holm, Christian; Shendruk, Tyler N
2016-04-01
A plethora of active matter models exist that describe the behavior of self-propelled particles (or swimmers), both with and without hydrodynamics. However, there are few studies that consider shape-anisotropic swimmers and include hydrodynamic interactions. Here, we introduce a simple method to simulate self-propelled colloids interacting hydrodynamically in a viscous medium using the lattice-Boltzmann technique. Our model is based on raspberry-type viscous coupling and a force/counter-force formalism, which ensures that the system is force free. We consider several anisotropic shapes and characterize their hydrodynamic multipolar flow field. We demonstrate that shape-anisotropy can lead to the presence of a strong quadrupole and octupole moments, in addition to the principle dipole moment. The ability to simulate and characterize these higher-order moments will prove crucial for understanding the behavior of model swimmers in confining geometries. PMID:27059561
Lattice-Boltzmann hydrodynamics of anisotropic active matter
NASA Astrophysics Data System (ADS)
de Graaf, Joost; Menke, Henri; Mathijssen, Arnold J. T. M.; Fabritius, Marc; Holm, Christian; Shendruk, Tyler N.
2016-04-01
A plethora of active matter models exist that describe the behavior of self-propelled particles (or swimmers), both with and without hydrodynamics. However, there are few studies that consider shape-anisotropic swimmers and include hydrodynamic interactions. Here, we introduce a simple method to simulate self-propelled colloids interacting hydrodynamically in a viscous medium using the lattice-Boltzmann technique. Our model is based on raspberry-type viscous coupling and a force/counter-force formalism, which ensures that the system is force free. We consider several anisotropic shapes and characterize their hydrodynamic multipolar flow field. We demonstrate that shape-anisotropy can lead to the presence of a strong quadrupole and octupole moments, in addition to the principle dipole moment. The ability to simulate and characterize these higher-order moments will prove crucial for understanding the behavior of model swimmers in confining geometries.
Finite-difference lattice-Boltzmann methods for binary fluids.
Xu, Aiguo
2005-06-01
We investigate two-fluid Bhatnagar-Gross-Krook (BGK) kinetic methods for binary fluids. The developed theory works for asymmetric as well as symmetric systems. For symmetric systems it recovers Sirovich's theory and is summarized in models A and B. For asymmetric systems it contributes models C, D, and E which are especially useful when the total masses and/or local temperatures of the two components are greatly different. The kinetic models are discretized based on an octagonal discrete velocity model. The discrete-velocity kinetic models and the continuous ones are required to describe the same hydrodynamic equations. The combination of a discrete-velocity kinetic model and an appropriate finite-difference scheme composes a finite-difference lattice Boltzmann method. The validity of the formulated methods is verified by investigating (i) uniform relaxation processes, (ii) isothermal Couette flow, and (iii) diffusion behavior. PMID:16089910
Full Eulerian lattice Boltzmann model for conjugate heat transfer
NASA Astrophysics Data System (ADS)
Hu, Yang; Li, Decai; Shu, Shi; Niu, Xiaodong
2015-12-01
In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently. The present conjugate interface treatment is validated by several steady and unsteady numerical tests, including pure heat conduction, forced convection, and natural convection problems. Both flat and curved interfaces are also involved. The obtained results show good agreement with the analytical and/or finite volume results.
Two relaxation time lattice Boltzmann model for rarefied gas flows
NASA Astrophysics Data System (ADS)
Esfahani, Javad Abolfazli; Norouzi, Ali
2014-01-01
In this paper, the lattice Boltzmann equation (LBE) with two relaxation times (TRT) is implemented in order to study gaseous flow through a long micro/nano-channel. A new relation is introduced for the reflection factor in the bounce-back/specular reflection (BSR) boundary condition based on the analytical solution of the Navier-Stokes equations. The focus of the present study is on comparing TRT with the other LBE models called multiple relaxation times (MRT) and single relaxation time (SRT) in simulation of rarefied gas flows. After a stability analysis for the TRT and SRT models, the numerical results are presented and validated by the analytical solution of the Navier-Stokes equations with slip boundary condition, direct simulation of Monte Carlo (DSMC) and information preservation (IP) method. The effect of various gases on flow behavior is also investigated by using the variable hard sphere (VHS) model through the symmetrical relaxation time.
Lattice Boltzmann method for mixtures at variable Schmidt number
NASA Astrophysics Data System (ADS)
Monteferrante, Michele; Melchionna, Simone; Marconi, Umberto Marini Bettolo
2014-07-01
When simulating multicomponent mixtures via the Lattice Boltzmann Method, it is desirable to control the mutual diffusivity between species while maintaining the viscosity of the solution fixed. This goal is herein achieved by a modification of the multicomponent Bhatnagar-Gross-Krook evolution equations by introducing two different timescales for mass and momentum diffusion. Diffusivity is thus controlled by an effective drag force acting between species. Numerical simulations confirm the accuracy of the method for neutral binary and charged ternary mixtures in bulk conditions. The simulation of a charged mixture in a charged slit channel show that the conductivity and electro-osmotic mobility exhibit a departure from the Helmholtz-Smoluchowski prediction at high diffusivity.
Lattice Boltzmann method for mixtures at variable Schmidt number.
Monteferrante, Michele; Melchionna, Simone; Marconi, Umberto Marini Bettolo
2014-07-01
When simulating multicomponent mixtures via the Lattice Boltzmann Method, it is desirable to control the mutual diffusivity between species while maintaining the viscosity of the solution fixed. This goal is herein achieved by a modification of the multicomponent Bhatnagar-Gross-Krook evolution equations by introducing two different timescales for mass and momentum diffusion. Diffusivity is thus controlled by an effective drag force acting between species. Numerical simulations confirm the accuracy of the method for neutral binary and charged ternary mixtures in bulk conditions. The simulation of a charged mixture in a charged slit channel show that the conductivity and electro-osmotic mobility exhibit a departure from the Helmholtz-Smoluchowski prediction at high diffusivity. PMID:25005272
Modeling of urban traffic networks with lattice Boltzmann model
NASA Astrophysics Data System (ADS)
Meng, Jian-ping; Qian, Yue-hong; Dai, Shi-qiang
2008-02-01
It is of great importance to uncover the characteristics of traffic networks. However, there have been few researches concerning kinetics models for urban traffic networks. In this work, a lattice Boltzmann model (LBM) for urban traffic networks is proposed by incorporating the ideas of the Biham-Middleton-Levine (BML) model into the LBM for road traffic. In the present model, situations at intersections with the red and green traffic signals are treated as a kind of boundary conditions varying with time. Thus, the urban traffic network could be described in the mesoscopic level. By performing numerical simulations under the periodic boundary conditions, the behavior of average velocity is investigated in detail. The numerical results agree quite well with those given by the Chowdhury-Schadschneider (ChSch) model (Chowdhury D. and Schadschneider A., Phys. Rev. E, 59 (1999) R1311). Furthermore, the statistical noise is reduced in this discrete kinetics model, thus, the present model has considerably high computational efficiency.
A lattice Boltzmann model for the Burgers-Fisher equation.
Zhang, Jianying; Yan, Guangwu
2010-06-01
A lattice Boltzmann model is developed for the one- and two-dimensional Burgers-Fisher equation based on the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different time scales. In order to obtain the two-dimensional Burgers-Fisher equation, vector sigma(j) has been used. And in order to overcome the drawbacks of "error rebound," a new assumption of additional distribution is presented, where two additional terms, in first order and second order separately, are used. Comparisons with the results obtained by other methods reveal that the numerical solutions obtained by the proposed method converge to exact solutions. The model under new assumption gives better results than that with second order assumption. PMID:20590325
Lattice Boltzmann model for generalized nonlinear wave equations
NASA Astrophysics Data System (ADS)
Lai, Huilin; Ma, Changfeng
2011-10-01
In this paper, a lattice Boltzmann model is developed to solve a class of the nonlinear wave equations. Through selecting equilibrium distribution function and an amending function properly, the governing evolution equation can be recovered correctly according to our proposed scheme, in which the Chapman-Enskog expansion is employed. We validate the algorithm on some problems where analytic solutions are available, including the second-order telegraph equation, the nonlinear Klein-Gordon equation, and the damped, driven sine-Gordon equation. It is found that the numerical results agree well with the analytic solutions, which indicates that the present algorithm is very effective and can be used to solve more general nonlinear problems.
Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation
NASA Technical Reports Server (NTRS)
Orr, Jeb S.; Powers, Joseph F.; Yang, Hong Q
2015-01-01
A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments is based upon a novel approach that relies on the global momentum conservation of the closed fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. A numerical example illustrates the method's application to prediction of bulk fluid behavior during a spacecraft ullage settling maneuver.
Lattice-Boltzmann Simulation of Coalescence-Driven Island Coarsening
Hakan Basagaoglu; Christopher T. Green; Paul Meakin; Benjamin J. McCoy
2004-10-01
A two-dimensional lattice-Boltzmann model (LBM) with fluid-fluid interactions was used to simulate first-order phase separation in a thin fluid film. The intermediate asymptotic time dependence of the mean island size, island number concentration, and polydispersity were determined and compared with the predictions of the distribution-kinetics model. The comparison revealed that the combined effects of growth, coalescence, and Ostwald ripening control the phase transition process in the LBM simulations. However, the overall process is dominated by coalescence, which is independent of island mass. As the phase transition advances, the mean island size increases, the number of islands decrease, and the polydispersity approaches unity, which conforms to the predictions of the distribution-kinetics model. The effects of the domain size on the intermediate asymptotic island size distribution, scaling form of the island size distribution, and the crossover to the long-term asymptotic behavior were elucidated. (C) 2004 American Institute of Physics.
Full Eulerian lattice Boltzmann model for conjugate heat transfer.
Hu, Yang; Li, Decai; Shu, Shi; Niu, Xiaodong
2015-12-01
In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently. The present conjugate interface treatment is validated by several steady and unsteady numerical tests, including pure heat conduction, forced convection, and natural convection problems. Both flat and curved interfaces are also involved. The obtained results show good agreement with the analytical and/or finite volume results. PMID:26764851
A Lattice Boltzmann Model for Oscillating Reaction-Diffusion
NASA Astrophysics Data System (ADS)
Rodríguez-Romo, Suemi; Ibañez-Orozco, Oscar; Sosa-Herrera, Antonio
2016-07-01
A computational algorithm based on the lattice Boltzmann method (LBM) is proposed to model reaction-diffusion systems. In this paper, we focus on how nonlinear chemical oscillators like Belousov-Zhabotinsky (BZ) and the chlorite-iodide-malonic acid (CIMA) reactions can be modeled by LBM and provide with new insight into the nature and applications of oscillating reactions. We use Gaussian pulse initial concentrations of sulfuric acid in different places of a bidimensional reactor and nondiffusive boundary walls. We clearly show how these systems evolve to a chaotic attractor and produce specific pattern images that are portrayed in the reactions trajectory to the corresponding chaotic attractor and can be used in robotic control.
Influence of asperities on fluid and thermal flow in a fracture: A coupled lattice Boltzmann study
NASA Astrophysics Data System (ADS)
Neuville, A.; Flekkøy, E. G.; Toussaint, R.
2013-07-01
The characteristics of the hydro-thermal flow which occurs when a cold fluid is injected into a hot fractured bedrock depend on the morphology of the fracture. We consider a sharp triangular asperity, invariant in one direction, perturbing an otherwise flat fracture. We investigate its influence on the macroscopic hydraulic transmissivity and heat transfer efficiency, at fixed low Reynolds number. In this study, numerical simulations are done with a coupled lattice Boltzmann method that solves both the complete Navier-Stokes and advection-diffusion equations in three dimensions. The results are compared with those obtained under lubrication approximations which rely on many hypotheses and neglect the three-dimensional (3-D) effects. The lubrication results are obtained by analytically solving the Stokes equation and a two-dimensional (integrated over the thickness) advection-diffusion equation. We use a lattice Boltzmann method with a double distribution (for mass and energy transport) on hypercubic and cubic lattices. Beyond some critical slope for the boundaries, the velocity profile is observed to be far from a quadratic profile in the vicinity of the sharp asperity: the fluid within the triangular asperity is quasi-static. We find that taking account of both the 3-D effects and the cooling of the rock, are important for the thermal exchange. Neglecting these effects with lubrication approximations results in overestimating the heat exchange efficiency. The evolution of the temperature over time, toward steady state, also shows complex behavior: some sites alternately reheat and cool down several times, making it difficult to forecast the extracted heat.
A Three-Dimensional Multi-Mesh Lattice Boltzmann Model for Multiphysics Simulations
NASA Astrophysics Data System (ADS)
Hashemi, Amirreza; Eshraghi, Mohsen; Felicelli, Sergio
2015-11-01
The lattice Boltzmann method (LBM) is known as an attractive computational method for modeling fluid flow and, more recently, transport phenomena. As any numerical method, the computational cost of LBM simulations depends on the density of the computational grids. The cost of simulations can become enormous when multiple equations are solved in three dimensions. In this work, the development of a multi-block multi-grid LBM model is discussed for three-dimensional (3D) multiphysics simulations. In a system of multiple coupled equations with different length scales, a multi-block mesh with different grids for each model would enhance the computational efficiency and stability of the model. Embedded-type grids facilitate the transfer of information between lattices while allowing larger time steps. In addition, a non-uniform mesh is considered within each mode that allows mesh refinement within each physical model when required. The multi-mesh method was developed to solve for transport phenomena including fluid flow, mass and heat transfer. The huge memory demands of LBM simulations in 3D was significantly reduced using this scheme. Moreover, by reducing the number of lattice points, cost communication in parallel processing was largely decreased.
Towards Full Aircraft Airframe Noise Prediction: Lattice Boltzmann Simulations
NASA Technical Reports Server (NTRS)
Khorrami, Mehdi R.; Fares, Ehab; Casalino, Damiano
2014-01-01
Computational results for an 18%-scale, semi-span Gulfstream aircraft model are presented. Exa Corporation's lattice Boltzmann PowerFLOW(trademark) solver was used to perform time-dependent simulations of the flow field associated with this high-fidelity aircraft model. The simulations were obtained for free-air at a Mach number of 0.2 with the flap deflected at 39 deg (landing configuration). We focused on accurately predicting the prominent noise sources at the flap tips and main landing gear for the two baseline configurations, namely, landing flap setting without and with gear deployed. Capitalizing on the inherently transient nature of the lattice Boltzmann formulation, the complex time-dependent flow features associated with the flap were resolved very accurately and efficiently. To properly simulate the noise sources over a broad frequency range, the tailored grid was very dense near the flap inboard and outboard tips. Extensive comparison of the computed time-averaged and unsteady surface pressures with wind tunnel measurements showed excellent agreement for the global aerodynamic characteristics and the local flow field at the flap inboard and outboard tips and the main landing gear. In particular, the computed fluctuating surface pressure field for the flap agreed well with the measurements in both amplitude and frequency content, indicating that the prominent airframe noise sources at the tips were captured successfully. Gear-flap interaction effects were remarkably well predicted and were shown to affect only the inboard flap tip, altering the steady and unsteady pressure fields in that region. The simulated farfield noise spectra for both baseline configurations, obtained using a Ffowcs-Williams and Hawkings acoustic analogy approach, were shown to be in close agreement with measured values.
Study of three-dimensional electro-osmotic flow with curved boundary via lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Chen, Q.; Zhang, X. B.; Li, Q.; Jiang, X. S.; Zhou, H. P.
2016-01-01
A three-dimensional (3D) lattice Boltzmann model and boundary method is developed to simulate electro-osmotic flow (EOF) with a charged spherical particle immersed in an electrolyte solution. The general governing equations for electro-osmotic transport are Navier-Stokes equations for fluid flow and the Poisson-Boltzmann equation for electric potential distribution around the particle. Two sets of D3Q19 lattice structure with curved boundary conditions are implemented. The simulation results are compared with analytical predictions and are found to be in excellent agreement. The potential distribution appears circularly symmetric and the flow velocity decreases with the cross-sectional area for flow passage increasing due to the mass conservation. The effects of the ionic concentration, the sphere radius, electric potential and external electric field on the velocity profiles are investigated. The flow velocity increases with both the electric potential and the external electric field. However, the variation in flow velocity with the ionic concentration and the sphere radius is complex due to the change in electrical double layer (EDL) thickness.
One-dimensional transient radiative transfer by lattice Boltzmann method.
Zhang, Yong; Yi, Hongliang; Tan, Heping
2013-10-21
The lattice Boltzmann method (LBM) is extended to solve transient radiative transfer in one-dimensional slab containing scattering media subjected to a collimated short laser irradiation. By using a fully implicit backward differencing scheme to discretize the transient term in the radiative transfer equation, a new type of lattice structure is devised. The accuracy and computational efficiency of this algorithm are examined firstly. Afterwards, effects of the medium properties such as the extinction coefficient, the scattering albedo and the anisotropy factor, and the shapes of laser pulse on time-resolved signals of transmittance and reflectance are investigated. Results of the present method are found to compare very well with the data from the literature. For an oblique incidence, the LBM results in this paper are compared with those by Monte Carlo method generated by ourselves. In addition, transient radiative transfer in a two-Layer inhomogeneous media subjected to a short square pulse irradiation is investigated. At last, the LBM is further extended to study the transient radiative transfer in homogeneous medium with a refractive index discontinuity irradiated by the short pulse laser. Several trends on the time-resolved signals different from those for refractive index of 1 (i.e. refractive-index-matched boundary) are observed and analysed. PMID:24150298
Volumetric lattice Boltzmann simulation for blood flow in aorta arteries
NASA Astrophysics Data System (ADS)
Deep, Debanjan; Yu, Huidan (Whitney); Teague, Shawn
2012-11-01
Complicated moving boundaries pose a major challenge in computational fluid dynamics for complex flows, especially in the biomechanics of both blood flow in the cardiovascular system and air flow in the respiratory system where the compliant nature of the vessels can have significant effects on the flow rate and wall shear stress. We develop a computation approach to treat arbitrarily moving boundaries using a volumetric representation of lattice Boltzmann method, which distributes fluid particles inside lattice cells. A volumetric bounce-back procedure is applied in the streaming step while momentum exchange between the fluid and moving solid boundary are accounted for in the collision sub-step. Additional boundary-induced migration is introduced to conserve fluid mass as the boundary moves across fluid cells. The volumetric LBM (VLBM) is used to simulate blood flow in both normal and dilated aorta arteries. We first compare flow structure and pressure distribution in steady state with results from Navier-Stokes based solver and good agreements are achieved. Then we focus on wall stress within the aorta for different heart pumping condition and present quantitative measurement of wall shear and normal stress.
Modeling groundwater flow by lattice Boltzmann method in curvilinear coordinates
NASA Astrophysics Data System (ADS)
Budinski, Ljubomir; Fabian, Julius; Stipic, Matija
2015-07-01
In order to promote the use of the lattice Boltzmann method (LBM) for the simulation of isotropic groundwater flow in a confined aquifer with arbitrary geometry, Poisson's equation was transformed into a curvilinear coordinate system. With the metric function between the physical and the computational domain established, Poisson's equation written in Cartesian coordinates was transformed in curvilinear coordinates. Following, the appropriate equilibrium function for the D2Q9 square lattice has been defined. The resulting curvilinear formulation of the LBM for groundwater flow is capable of modeling flow in domains of complex geometry with the opportunity of local refining/coarsening of the computational mesh corresponding to the complexity of the flow pattern and the required accuracy. Since the proposed form of the LBM uses the transformed equation of flow implemented in the equilibrium function, finding a solution does not require supplementary procedures along the curvilinear boundaries, nor in the zones requiring mesh density adjustments. Thus, the basic concept of the LBM is completely maintained. The improvement of the proposed LBM over the previously published classical methods is completely verified by three examples with analytical solutions. The results demonstrate the advantages of the proposed curvilinear LBM in modeling groundwater flow in complex flow domains.
Evaluation of the Lattice-Boltzmann Equation Solver PowerFLOW for Aerodynamic Applications
NASA Technical Reports Server (NTRS)
Lockard, David P.; Luo, Li-Shi; Singer, Bart A.; Bushnell, Dennis M. (Technical Monitor)
2000-01-01
A careful comparison of the performance of a commercially available Lattice-Boltzmann Equation solver (Power-FLOW) was made with a conventional, block-structured computational fluid-dynamics code (CFL3D) for the flow over a two-dimensional NACA-0012 airfoil. The results suggest that the version of PowerFLOW used in the investigation produced solutions with large errors in the computed flow field; these errors are attributed to inadequate resolution of the boundary layer for reasons related to grid resolution and primitive turbulence modeling. The requirement of square grid cells in the PowerFLOW calculations limited the number of points that could be used to span the boundary layer on the wing and still keep the computation size small enough to fit on the available computers. Although not discussed in detail, disappointing results were also obtained with PowerFLOW for a cavity flow and for the flow around a generic helicopter configuration.
Lattice Uehling-Uhlenbeck Boltzmann-Bhatnagar-Gross-Krook hydrodynamics of quantum gases
NASA Astrophysics Data System (ADS)
Yang, Jaw-Yen; Hung, Li-Hsin
2009-05-01
We present a semiclassical lattice Boltzmann method based on quantum kinetic theory. The method is directly derived by projecting the Uehling-Uhlenbeck Boltzmann-Bhatnagar-Gross-Krook equations onto the tensor Hermite polynomials following Grad’s moment expansion method. The intrinsic discrete nodes of the Gauss-Hermite quadrature provide the natural lattice velocities for the semiclassical lattice Boltzmann method. Gases of particles of arbitrary statistics can be considered. Simulation of one-dimensional compressible gas flow and two-dimensional hydrodynamic flows are shown. The results indicate distinct characteristics of the effects of quantum statistics.
Study of hydrodynamic instabilities with a multiphase lattice Boltzmann model
NASA Astrophysics Data System (ADS)
Velasco, Ali Mauricio; Muñoz, José Daniel
2015-10-01
Rayleigh-Taylor and Kelvin-Helmholtz hydrodynamic instabilities are frequent in many natural and industrial processes, but their numerical simulation is not an easy challenge. This work simulates both instabilities by using a lattice Boltzmann model on multiphase fluids at a liquid-vapour interface, instead of multicomponent systems like the oil-water one. The model, proposed by He, Chen and Zhang (1999) [1] was modified to increase the precision by computing the pressure gradients with a higher order, as proposed by McCracken and Abraham (2005) [2]. The resulting model correctly simulates both instabilities by using almost the same parameter set. It also reproduces the relation γ ∝√{ A} between the growing rate γ of the Rayleigh-Taylor instability and the relative density difference between the fluids (known as the Atwood number A), but including also deviations observed in experiments at low density differences. The results show that the implemented model is a useful tool for the study of hydrodynamic instabilities, drawing a sharp interface and exhibiting numerical stability for moderately high Reynolds numbers.
Force Evaluation in the Lattice Boltzmann Method Involving Curved Geometry
NASA Technical Reports Server (NTRS)
Mei, Renwei; Yu, Dazhi; Shyy, Wei; Luo, Li-Shi; Bushnell, Dennis M. (Technical Monitor)
2002-01-01
The present work investigates two approaches for force evaluation in the lattice Boltzmann equation: the momentum- exchange method and the stress-integration method on the surface of a body. The boundary condition for the particle distribution functions on curved geometries is handled with second order accuracy based on our recent works. The stress-integration method is computationally laborious for two-dimensional flows and in general difficult to implement for three-dimensional flows, while the momentum-exchange method is reliable, accurate, and easy to implement for both two-dimensional and three-dimensional flows. Several test cases are selected to evaluate the present methods, including: (i) two-dimensional pressure-driven channel flow; (ii) two-dimensional uniform flow past a column of cylinders; (iii) two-dimensional flow past a cylinder asymmetrically placed in a channel (with vortex shedding); (iv) three-dimensional pressure-driven flow in a circular pipe; and (v) three-dimensional flow past a sphere. The drag evaluated by using the momentum-exchange method agrees well with the exact or other published results.
Lattice Boltzmann modeling of directional wetting: comparing simulations to experiments.
Jansen, H Patrick; Sotthewes, Kai; van Swigchem, Jeroen; Zandvliet, Harold J W; Kooij, E Stefan
2013-07-01
Lattice Boltzmann Modeling (LBM) simulations were performed on the dynamic behavior of liquid droplets on chemically striped patterned surfaces, ultimately with the aim to develop a predictive tool enabling reliable design of future experiments. The simulations accurately mimic experimental results, which have shown that water droplets on such surfaces adopt an elongated shape due to anisotropic preferential spreading. Details of the contact line motion such as advancing of the contact line in the direction perpendicular to the stripes exhibit pronounced similarities in experiments and simulations. The opposite of spreading, i.e., evaporation of water droplets, leads to a characteristic receding motion first in the direction parallel to the stripes, while the contact line remains pinned perpendicular to the stripes. Only when the aspect ratio is close to unity, the contact line also starts to recede in the perpendicular direction. Very similar behavior was observed in the LBM simulations. Finally, droplet movement can be induced by a gradient in surface wettability. LBM simulations show good semiquantitative agreement with experimental results of decanol droplets on a well-defined striped gradient, which move from high- to low-contact angle surfaces. Similarities and differences for all systems are described and discussed in terms of the predictive capabilities of LBM simulations to model direction wetting. PMID:23944550
Lattice Boltzmann modeling of three-phase incompressible flows.
Liang, H; Shi, B C; Chai, Z H
2016-01-01
In this paper, based on multicomponent phase-field theory we intend to develop an efficient lattice Boltzmann (LB) model for simulating three-phase incompressible flows. In this model, two LB equations are used to capture the interfaces among three different fluids, and another LB equation is adopted to solve the flow field, where a new distribution function for the forcing term is delicately designed. Different from previous multiphase LB models, the interfacial force is not used in the computation of fluid velocity, which is more reasonable from the perspective of the multiscale analysis. As a result, the computation of fluid velocity can be much simpler. Through the Chapman-Enskog analysis, it is shown that the present model can recover exactly the physical formulations for the three-phase system. Numerical simulations of extensive examples including two circular interfaces, ternary spinodal decomposition, spreading of a liquid lens, and Kelvin-Helmholtz instability are conducted to test the model. It is found that the present model can capture accurate interfaces among three different fluids, which is attributed to its algebraical and dynamical consistency properties with the two-component model. Furthermore, the numerical results of three-phase flows agree well with the theoretical results or some available data, which demonstrates that the present LB model is a reliable and efficient method for simulating three-phase flow problems. PMID:26871191
Consistent lattice Boltzmann methods for incompressible axisymmetric flows.
Zhang, Liangqi; Yang, Shiliang; Zeng, Zhong; Yin, Linmao; Zhao, Ya; Chew, Jia Wei
2016-08-01
In this work, consistent lattice Boltzmann (LB) methods for incompressible axisymmetric flows are developed based on two efficient axisymmetric LB models available in the literature. In accord with their respective original models, the proposed axisymmetric models evolve within the framework of the standard LB method and the source terms contain no gradient calculations. Moreover, the incompressibility conditions are realized with the Hermite expansion, thus the compressibility errors arising in the existing models are expected to be reduced by the proposed incompressible models. In addition, an extra relaxation parameter is added to the Bhatnagar-Gross-Krook collision operator to suppress the effect of the ghost variable and thus the numerical stability of the present models is significantly improved. Theoretical analyses, based on the Chapman-Enskog expansion and the equivalent moment system, are performed to derive the macroscopic equations from the LB models and the resulting truncation terms (i.e., the compressibility errors) are investigated. In addition, numerical validations are carried out based on four well-acknowledged benchmark tests and the accuracy and applicability of the proposed incompressible axisymmetric LB models are verified. PMID:27627407
Thermal multicomponent lattice Boltzmann model for catalytic reactive flows.
Kang, Jinfen; Prasianakis, Nikolaos I; Mantzaras, John
2014-06-01
Catalytic reactions are of great interest in many applications related to power generation, fuel reforming and pollutant abatement, as well as in various biochemical processes. A recently proposed lattice Boltzmann model for thermal binary-mixture gas flows [J. Kang, N. I. Prasianakis, and J. Mantzaras, Phys. Rev. E. 87, 053304 (2013)] is revisited and extended for the simulation of multispecies flows with catalytic reactions. The resulting model can handle flows with large temperature and concentration gradients. The developed model is presented in detail and validated against a finite volume Navier-Stokes solver in the case of channel-flow methane catalytic combustion. The surface chemistry is treated with a one-step global reaction for the catalytic total oxidation of methane on platinum. In order to take into account thermal effects, the catalytic boundary condition of S. Arcidiacono, J. Mantzaras, and I. V. Karlin [Phys. Rev. E 78, 046711 (2008)] is adapted to account for temperature variations. Speed of sound simulations further demonstrate the physical integrity and unique features of the model. PMID:25019915
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Pazdniakou, A.; Adler, P. M.
2013-12-01
The lattice Boltzmann method (LBM) is applied to calculate the dynamic permeability K(ω) of porous media; an oscillating macroscopic pressure gradient is imposed in order to generate oscillating flows. The LBM simulation yields the time dependent seepage velocity of amplitude A and phase shift B which are used to calculate K(ω). The procedure is validated for plane Poiseuille flows where excellent agreement with the analytical solution is obtained. The limitations of the method are discussed. When the ratio between the kinematic viscosity and the characteristic size of the pores is high, the corresponding Knudsen number Kn is high and the numerical values of K(ω) are incorrect with a positive imaginary part; it is only when Kn is small enough that correct values are obtained. The influence of the time discretization of the oscillating body force is studied; simulation results are influenced by an insufficient discretization, i.e., it is necessary to avoid using too high frequencies. The influence of absolute errors in the seepage velocity amplitude δA and the phase shift δB on K(ω) shows that for high ω even small errors in B can cause drastic errors in ReK(ω). The dynamic permeability of reconstructed and real (sandstone) porous media is calculated for a large range of frequencies and the universal scaling behavior is verified. Very good correspondences with the theoretical predictions are observed.
Lattice Boltzmann modeling of three-phase incompressible flows
NASA Astrophysics Data System (ADS)
Liang, H.; Shi, B. C.; Chai, Z. H.
2016-01-01
In this paper, based on multicomponent phase-field theory we intend to develop an efficient lattice Boltzmann (LB) model for simulating three-phase incompressible flows. In this model, two LB equations are used to capture the interfaces among three different fluids, and another LB equation is adopted to solve the flow field, where a new distribution function for the forcing term is delicately designed. Different from previous multiphase LB models, the interfacial force is not used in the computation of fluid velocity, which is more reasonable from the perspective of the multiscale analysis. As a result, the computation of fluid velocity can be much simpler. Through the Chapman-Enskog analysis, it is shown that the present model can recover exactly the physical formulations for the three-phase system. Numerical simulations of extensive examples including two circular interfaces, ternary spinodal decomposition, spreading of a liquid lens, and Kelvin-Helmholtz instability are conducted to test the model. It is found that the present model can capture accurate interfaces among three different fluids, which is attributed to its algebraical and dynamical consistency properties with the two-component model. Furthermore, the numerical results of three-phase flows agree well with the theoretical results or some available data, which demonstrates that the present LB model is a reliable and efficient method for simulating three-phase flow problems.
Multiblock approach for the passive scalar thermal lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Huang, Rongzong; Wu, Huiying
2014-04-01
A multiblock approach for the passive scalar thermal lattice Boltzmann method (TLBM) with multiple-relaxation-time collision scheme is proposed based on the Chapman-Enskog analysis. The interaction between blocks is executed in the moment space directly and an external force term is considered. Theoretical analysis shows that all the nonequilibrium parts of the nonconserved moments should be rescaled, while the nonequilibrium parts of the conserved moments can be calculated directly. Moreover, a local scheme based on the pseudoparticles for computing heat flux is proposed with no need to calculate temperature gradient based on the finite-difference scheme. In order to validate the multiblock approach and local scheme for computing heat flux, thermal Couette flow with wall injection is simulated and good results are obtained, which show that the adoption of the multiblock approach does not deteriorate the convergence rate of TLBM and the local scheme for computing heat flux has second-order convergence rate. Further application of the present approach is the simulation of natural convection in a square cavity with the Rayleigh number up to 109.
Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows.
Li, Q; Luo, K H; Li, X J
2012-07-01
The pseudopotential lattice Boltzmann (LB) model is a widely used multiphase model in the LB community. In this model, an interaction force, which is usually implemented via a forcing scheme, is employed to mimic the molecular interactions that cause phase segregation. The forcing scheme is therefore expected to play an important role in the pseudoepotential LB model. In this paper, we aim to address some key issues about forcing schemes in the pseudopotential LB model. First, theoretical and numerical analyses will be made for Shan-Chen's forcing scheme [Shan and Chen, Phys. Rev. E 47, 1815 (1993)] and the exact-difference-method forcing scheme [Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)]. The nature of these two schemes and their recovered macroscopic equations will be shown. Second, through a theoretical analysis, we will reveal the physics behind the phenomenon that different forcing schemes exhibit different performances in the pseudopotential LB model. Moreover, based on the analysis, we will present an improved forcing scheme and numerically demonstrate that the improved scheme can be treated as an alternative approach to achieving thermodynamic consistency in the pseudopotential LB model. PMID:23005565
Force evaluation in the lattice Boltzmann method involving curved geometry
NASA Astrophysics Data System (ADS)
Mei, Renwei; Yu, Dazhi; Shyy, Wei; Luo, Li-Shi
2002-04-01
The present work investigates two approaches for force evaluation in the lattice Boltzmann equation: the momentum-exchange method and the stress-integration method on the surface of a body. The boundary condition for the particle distribution functions on curved geometries is handled with second-order accuracy based on our recent works [Mei et al., J. Comput. Phys. 155, 307 (1999); ibid. 161, 680 (2000)]. The stress-integration method is computationally laborious for two-dimensional flows and in general difficult to implement for three-dimensional flows, while the momentum-exchange method is reliable, accurate, and easy to implement for both two-dimensional and three-dimensional flows. Several test cases are selected to evaluate the present methods, including: (i) two-dimensional pressure-driven channel flow; (ii) two-dimensional uniform flow past a column of cylinders; (iii) two-dimensional flow past a cylinder asymmetrically placed in a channel (with vortex shedding); (iv) three-dimensional pressure-driven flow in a circular pipe; and (v) three-dimensional flow past a sphere. The drag evaluated by using the momentum-exchange method agrees well with the exact or other published results.
Lattice Boltzmann simulations of settling behaviors of irregularly shaped particles
NASA Astrophysics Data System (ADS)
Zhang, Pei; Galindo-Torres, S. A.; Tang, Hongwu; Jin, Guangqiu; Scheuermann, A.; Li, Ling
2016-06-01
We investigated the settling dynamics of irregularly shaped particles in a still fluid under a wide range of conditions with Reynolds numbers Re varying between 1 and 2000, sphericity ϕ and circularity c both greater than 0.5, and Corey shape factor (CSF) less than 1. To simulate the particle settling process, a modified lattice Boltzmann model combined with a turbulence module was adopted. This model was first validated using experimental data for particles of spherical and cubic shapes. For irregularly shaped particles, two different types of settling behaviors were observed prior to particles reaching a steady state: accelerating and accelerating-decelerating, which could be distinguished by a critical CSF value of approximately 0.7. The settling dynamics were analyzed with a focus on the projected areas and angular velocities of particles. It was found that a minor change in the starting projected area, an indicator of the initial particle orientation, would not strongly affect the settling velocity for low Re. Periodic oscillations developed for all simulated particles when Re>100 . The amplitude of these oscillations increased with Re. However, the periods were not sensitive to Re. The critical Re that defined the transition between the steady and periodically oscillating behaviors depended on the inertia tensor. In particular, the maximum eigenvalue of the inertia tensor played a major role in signaling this transition in comparison to the intermediate and minimum eigenvalues.
Lattice Boltzmann Modeling of Micro-fluidic Devices
Clague, D S
2002-01-28
The results to date do indeed show that the lattice Boltzmann method accurately solves relevant, non-trivial flow problems. The parallelization of both the fluid and the mobile species in flow has enhanced this capability such that it is useful for solving relevant problems in a timely fashion. The initial studies of stationary or capture species revealed evidence of hydrodynamic screening between upstream and downstream particles. Numerical studies reveal that the critical length for which the test particle is hydrodynamically decoupled from upstream and downstream particles is on the order of 30 sphere radii. For mobile species, the LB capability was shown to be naturally suited for predicting the hydrodynamic lift phenomenon (inertial lift). A conversion factor was developed based on scaling arguments to include relevant forces generated by external fields. Using this conversion, an analytic solution for the Dielectrophoretic force was included into the LB capability which enabled the study of Dielectrophoretic particle capture. The Non-Newtonian enhancements have expanded the applicability of the LB capability to more physical systems. Specifically, with the bead-n-spring representation of macromolecules researchers will be able to study chain dynamics in micro-, physiological and Bio-MEMS environments. Furthermore, the ability to capture the shear thinning behavior, without any increase in computational time, positions this capability to be applied to a whole host of new problems involving biofluids.
Thermal multicomponent lattice Boltzmann model for catalytic reactive flows
NASA Astrophysics Data System (ADS)
Kang, Jinfen; Prasianakis, Nikolaos I.; Mantzaras, John
2014-06-01
Catalytic reactions are of great interest in many applications related to power generation, fuel reforming and pollutant abatement, as well as in various biochemical processes. A recently proposed lattice Boltzmann model for thermal binary-mixture gas flows [J. Kang, N. I. Prasianakis, and J. Mantzaras, Phys. Rev. E. 87, 053304 (2013), 10.1103/PhysRevE.87.053304] is revisited and extended for the simulation of multispecies flows with catalytic reactions. The resulting model can handle flows with large temperature and concentration gradients. The developed model is presented in detail and validated against a finite volume Navier-Stokes solver in the case of channel-flow methane catalytic combustion. The surface chemistry is treated with a one-step global reaction for the catalytic total oxidation of methane on platinum. In order to take into account thermal effects, the catalytic boundary condition of S. Arcidiacono, J. Mantzaras, and I. V. Karlin [Phys. Rev. E 78, 046711 (2008), 10.1103/PhysRevE.78.046711] is adapted to account for temperature variations. Speed of sound simulations further demonstrate the physical integrity and unique features of the model.
Lattice Boltzmann simulations of drops colliding with solid surfaces
NASA Astrophysics Data System (ADS)
Jia, X.; McLaughlin, J. B.; Kontomaris, K.
2009-04-01
Video images of drops colliding with solid surfaces shown by Rioboo et al. (2002) reveal that, for large drop velocities, the drops flatten and form a ring structure before receding and, in some cases, rebounding from the surface. They described the sequence of events in terms of four distinct regimes. During the initial kinematic phase, the dimensionless wetting radius of the drop follows a universal form if the drop Weber and Reynolds numbers are sufficiently large. In the second phase, the drop becomes highly flattened and the values of the Weber and Reynolds numbers influence the time evolution of the dimensionless wetting radius and its maximum value. This is followed by a third phase in which the wetting radius begins to decrease with time and the wettability of the surface influences the dynamics. This paper presents simulation results for the early stages of drop impact and spreading on a partially wetting solid surface. The simulations were performed with a modified version of the lattice Boltzmann method (LBM) developed by Inamuro et al. (2004) for a liquid-gas density ratio of 1000. The Inamuro et al. version of the LBM was modified by incorporating rigid, no-slip boundary conditions and incorporating a boundary condition on the normal derivative of the order parameter to impose the desired equilibrium contact angle.
Massively parallel simulations of multiphase flows using Lattice Boltzmann methods
NASA Astrophysics Data System (ADS)
Ahrenholz, Benjamin
2010-03-01
In the last two decades the lattice Boltzmann method (LBM) has matured as an alternative and efficient numerical scheme for the simulation of fluid flows and transport problems. Unlike conventional numerical schemes based on discretizations of macroscopic continuum equations, the LBM is based on microscopic models and mesoscopic kinetic equations. The fundamental idea of the LBM is to construct simplified kinetic models that incorporate the essential physics of microscopic or mesoscopic processes so that the macroscopic averaged properties obey the desired macroscopic equations. Especially applications involving interfacial dynamics, complex and/or changing boundaries and complicated constitutive relationships which can be derived from a microscopic picture are suitable for the LBM. In this talk a modified and optimized version of a Gunstensen color model is presented to describe the dynamics of the fluid/fluid interface where the flow field is based on a multi-relaxation-time model. Based on that modeling approach validation studies of contact line motion are shown. Due to the fact that the LB method generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallelization. Hence, it is possible to perform efficient simulations in complex geometries at a large scale by massively parallel computations. Here, the results of drainage and imbibition (Degree of Freedom > 2E11) in natural porous media gained from microtomography methods are presented. Those fully resolved pore scale simulations are essential for a better understanding of the physical processes in porous media and therefore important for the determination of constitutive relationships.
Lattice Boltzmann method for one-dimensional vector radiative transfer.
Zhang, Yong; Yi, Hongliang; Tan, Heping
2016-02-01
A one-dimensional vector radiative transfer (VRT) model based on lattice Boltzmann method (LBM) that considers polarization using four Stokes parameters is developed. The angular space is discretized by the discrete-ordinates approach, and the spatial discretization is conducted by LBM. LBM has such attractive properties as simple calculation procedure, straightforward and efficient handing of boundary conditions, and capability of stable and accurate simulation. To validate the performance of LBM for vector radiative transfer, four various test problems are examined. The first case investigates the non-scattering thermal-emitting atmosphere with no external collimated solar. For the other three cases, the external collimated solar and three different scattering types are considered. Particularly, the LBM is extended to solve VRT in the atmospheric aerosol system where the scattering function contains singularities and the hemisphere space distributions for the Stokes vector are presented and discussed. The accuracy and computational efficiency of this algorithm are discussed. Numerical results show that the LBM is accurate, flexible and effective to solve one-dimensional polarized radiative transfer problems. PMID:26906779
Large Eddy Simulations using Lattice Boltzmann algorithms. Final report
Serling, J.D.
1993-09-28
This report contains the results of a study performed to implement eddy-viscosity models for Large-Eddy-Simulations (LES) into Lattice Boltzmann (LB) algorithms for simulating fluid flows. This implementation requires modification of the LB method of simulating the incompressible Navier-Stokes equations to allow simulation of the filtered Navier-Stokes equations with some subgrid model for the Reynolds stress term. We demonstrate that the LB method can indeed be used for LES by simply locally adjusting the value of the BGK relaxation time to obtain the desired eddy-viscosity. Thus, many forms of eddy-viscosity models including the standard Smagorinsky model or the Dynamic model may be implemented using LB algorithms. Since underresolved LB simulations often lead to instability, the LES model actually serves to stabilize the method. An alternative method of ensuring stability is presented which requires that entropy increase during the collision step of the LB method. Thus, an alternative collision operator is locally applied if the entropy becomes too low. This stable LB method then acts as an LES scheme that effectively introduces its own eddy viscosity to damp short wavelength oscillations.
Multiple anisotropic collisions for advection-diffusion Lattice Boltzmann schemes
NASA Astrophysics Data System (ADS)
Ginzburg, Irina
2013-01-01
This paper develops a symmetrized framework for the analysis of the anisotropic advection-diffusion Lattice Boltzmann schemes. Two main approaches build the anisotropic diffusion coefficients either from the anisotropic anti-symmetric collision matrix or from the anisotropic symmetric equilibrium distribution. We combine and extend existing approaches for all commonly used velocity sets, prescribe most general equilibrium and build the diffusion and numerical-diffusion forms, then derive and compare solvability conditions, examine available anisotropy and stable velocity magnitudes in the presence of advection. Besides the deterioration of accuracy, the numerical diffusion dictates the stable velocity range. Three techniques are proposed for its elimination: (i) velocity-dependent relaxation entries; (ii) their combination with the coordinate-link equilibrium correction; and (iii) equilibrium correction for all links. Two first techniques are also available for the minimal (coordinate) velocity sets. Even then, the two-relaxation-times model with the isotropic rates often gains in effective stability and accuracy. The key point is that the symmetric collision mode does not modify the modeled diffusion tensor but it controls the effective accuracy and stability, via eigenvalue combinations of the opposite parity eigenmodes. We propose to reduce the eigenvalue spectrum by properly combining different anisotropic collision elements. The stability role of the symmetric, multiple-relaxation-times component, is further investigated with the exact von Neumann stability analysis developed in diffusion-dominant limit.
Consistent lattice Boltzmann methods for incompressible axisymmetric flows
NASA Astrophysics Data System (ADS)
Zhang, Liangqi; Yang, Shiliang; Zeng, Zhong; Yin, Linmao; Zhao, Ya; Chew, Jia Wei
2016-08-01
In this work, consistent lattice Boltzmann (LB) methods for incompressible axisymmetric flows are developed based on two efficient axisymmetric LB models available in the literature. In accord with their respective original models, the proposed axisymmetric models evolve within the framework of the standard LB method and the source terms contain no gradient calculations. Moreover, the incompressibility conditions are realized with the Hermite expansion, thus the compressibility errors arising in the existing models are expected to be reduced by the proposed incompressible models. In addition, an extra relaxation parameter is added to the Bhatnagar-Gross-Krook collision operator to suppress the effect of the ghost variable and thus the numerical stability of the present models is significantly improved. Theoretical analyses, based on the Chapman-Enskog expansion and the equivalent moment system, are performed to derive the macroscopic equations from the LB models and the resulting truncation terms (i.e., the compressibility errors) are investigated. In addition, numerical validations are carried out based on four well-acknowledged benchmark tests and the accuracy and applicability of the proposed incompressible axisymmetric LB models are verified.
Lattice Boltzmann simulation for forced two-dimensional turbulence.
Xia, YuXian; Qian, YueHong
2014-08-01
The direct numerical simulations of forced two-dimensional turbulent flow are presented by using the lattice Boltzmann method. The development of an energy-enstrophy double cascade is investigated in the two cases of external force of two-dimensional turbulence, Gaussian force and Kolmogorov force. It is found that the friction force is a necessary condition of the occurrence of a double cascade. The energy spectrum k(-3) in the enstrophy inertial range is in accord with the classical Kraichnan theory for both external forces. The energy spectrum of the Gaussian force case in an inverse cascade is k(-2); however, the Kolmogorov force drives the k(-5/3) energy in a backscatter cascade. The result agrees with Scott's standpoint, which describes nonrobustness of the two-dimensional turbulent inverse cascade. Also, intermittency is found for the enstrophy cascade in two cases of the external force form. Intermittency refers to the nonuniform distribution of saddle points in the two-dimensional turbulent flow. PMID:25215817
Lattice Boltzmann simulation for forced two-dimensional turbulence
NASA Astrophysics Data System (ADS)
Xia, YuXian; Qian, YueHong
2014-08-01
The direct numerical simulations of forced two-dimensional turbulent flow are presented by using the lattice Boltzmann method. The development of an energy-enstrophy double cascade is investigated in the two cases of external force of two-dimensional turbulence, Gaussian force and Kolmogorov force. It is found that the friction force is a necessary condition of the occurrence of a double cascade. The energy spectrum k-3 in the enstrophy inertial range is in accord with the classical Kraichnan theory for both external forces. The energy spectrum of the Gaussian force case in an inverse cascade is k-2; however, the Kolmogorov force drives the k-5/3 energy in a backscatter cascade. The result agrees with Scott's standpoint, which describes nonrobustness of the two-dimensional turbulent inverse cascade. Also, intermittency is found for the enstrophy cascade in two cases of the external force form. Intermittency refers to the nonuniform distribution of saddle points in the two-dimensional turbulent flow.
Lattice Boltzmann Simulation Optimization on Leading Multicore Platforms
Williams, Samuel; Carter, Jonathan; Oliker, Leonid; Shalf, John; Yelick, Katherine
2008-02-01
We present an auto-tuning approach to optimize application performance on emerging multicore architectures. The methodology extends the idea of search-based performance optimizations, popular in linear algebra and FFT libraries, to application-specific computational kernels. Our work applies this strategy to a lattice Boltzmann application (LBMHD) that historically has made poor use of scalar microprocessors due to its complex data structures and memory access patterns. We explore one of the broadest sets of multicore architectures in the HPC literature, including the Intel Clovertown, AMD Opteron X2, Sun Niagara2, STI Cell, as well as the single core Intel Itanium2. Rather than hand-tuning LBMHD for each system, we develop a code generator that allows us identify a highly optimized version for each platform, while amortizing the human programming effort. Results show that our auto-tuned LBMHD application achieves up to a 14x improvement compared with the original code. Additionally, we present detailed analysis of each optimization, which reveal surprising hardware bottlenecks and software challenges for future multicore systems and applications.
Treatment of moving boundaries in lattice-Boltzmann simulations.
NASA Astrophysics Data System (ADS)
Indireshkumar, K.; Pal, A.; Brasseur, J. G.
2000-11-01
We consider the treatment of moving boundaries with the lattice-Boltzmann (LB) technique, where the treatment of the boundary often does not precisely conserve mass and spurious fluctuations in density/pressure result from boundary motion through fixed grids. First, we applied the extrapolation method proposed by Chen et. al.(S. Y. Chen, D. Martinez, and R Mei, Phys. Fluids) 8, 2527 (1996) to incompressible flow induced by the movement of a piston in a 2D ``cylinder'' with mass flow out of or into the cylinder. In these simulations, the velocity of the boundary nodes is set equal to the (known) velocity of the boundary (piston) in the equilibrium distribution function (Method I). In a second set of simulations, the boundary node velocities are obtained by interpolating between interior nodes and the boundary, thus including the effect of boundary position more precisely (Method II). Comparison of LB predictions with simulations using FIDAP show pressure agreement to witnin 2 %. The total mass is conserved to within 0.1% with Method I and improves to within 0.02 % using method II. Spurious fluctuations in density/pressure due to boundary movement is about 0.9% with Method I, which improves significantly to about 0.3% with Method II. The application of these simple techniques to more complex geometries and wall (and fluid) motions in a stomach during gastric emptying will be presented.
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Adler, P.; Pazdniakou, A.
2012-04-01
The main objective of our work is to determine the dynamic permeability of three dimensional porous media by means of the Lattice Boltzmann method (LBM). The Navier-Stokes equation can be numerically solved by LBM which is widely used to address various fluid dynamics problems. Space is discretized by a three-dimensional cubic lattice and time is discretized as well. The generally accepted notation for lattice Boltzmann models is DdQq where D stands for space dimension and Q for the number of discrete velocities. The present model is denoted by D3Q19. Moreover, the Two Relaxation Times variant of the Multi Relaxation Times model is implemented. Bounce back boundary conditions are used on the solid-fluid interfaces. The porous medium is spatially periodic. Reconstructed media were used; they are obtained by imposing a porosity and a correlation function characterized by a correlation length. Real samples can be obtained by MicroCT. In contrast with other previous contributions, the dynamic permeability K(omega) which is a complex number, is derived by imposing an oscillating body force of pulsation omega on the unit cell and by deriving the amplitude and the phase shift of the resulting time dependent seepage velocity. The influence of two limiting parameters, namely the Knudsen number Kn and the discretization for high frequencies, on K(omega) is carefully studied for the first time. Kn is proportional to nu/(cs H) where nu is the kinematic viscosity, cs the speed of sound in the fluid and H a characteristic length scale of the porous medium. Several porous media such as the classical plane Poiseuille flow and the reconstructed media are used to show that it is only for small enough values of Kn that reliable results are obtained. Otherwise, the data depend on Kn and may even be totally unphysical. However, it should be noticed that the limiting value of Kn could not be derived in general since it depends very much on the structure of the medium. Problems occur at
Dynamic permeability of porous media by the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Pazdniakou, A.; Adler, P. M.
2011-12-01
The main objective of our work is to determine the dynamic permeability of three dimensional porous media by means of the Lattice Boltzmann method (LBM). The Navier-Stokes equation can be numerically solved by LBM which is widely used to address various fluid dynamics problems. Space is discretized by a three-dimensional cubic lattice and time is discretized as well. The generally accepted notation for lattice Boltzmann models is DdQq where D stands for space dimension and Q for the number of discrete velocities. The present model is denoted by D3Q19. Moreover, the Two Relaxation Times variant of the Multi Relaxation Times model is implemented. Bounce back boundary conditions are used on the solid-fluid interfaces. The porous medium is spatially periodic. Reconstructed media were used; they are obtained by imposing a porosity and a correlation function characterized by a correlation length. Real samples can be obtained by MicroCT. In contrast with other previous contributions, the dynamic permeability K(omega) which is a complex number, is derived by imposing an oscillating body force of pulsation omega on the unit cell and by deriving the amplitude and the phase shift of the resulting time dependent seepage velocity. The influence of two limiting parameters, namely the Knudsen number Kn and the discretization for high frequencies, on K(omega) is carefully studied for the first time. Kn is proportional to nu/(c_s H) where nu is the kinematic viscosity, c_s the speed of sound in the fluid and H a characteristic length scale of the porous medium. Several porous media such as the classical plane Poiseuille flow and the reconstructed media are used to show that it is only for small enough values of Kn that reliable results are obtained. Otherwise, the data depend on Kn and may even be totally unphysical. However, it should be noticed that the limiting value of Kn could not be derived in general since it depends very much on the structure of the medium. Problems occur
NASA Astrophysics Data System (ADS)
Zhang, Jianying; Yan, Guangwu
2016-04-01
A lattice Boltzmann model for solving the (2+1) dimensional cubic-quintic complex Ginzburg-Landau equation (CQCGLE) is proposed. Different from the classic lattice Boltzmann models, this lattice Boltzmann model is based on uniformly distributed lattice points in a two-dimensional space, and the evolution of the model is about a spatial axis rather than time. The algorithm provides advantages similar to the lattice Boltzmann method in that it is easily adapted to complex Ginzburg-Landau equations. Numerical results reproduce the phenomena of the fusion of necklace-ring pattern and the effect of non-linearity on the soliton in the CQCGLE.
NASA Astrophysics Data System (ADS)
Zhang, Jianying; Yan, Guangwu
2015-12-01
A spatiotemporal lattice Boltzmann model for solving the three-dimensional cubic-quintic complex Ginzburg-Landau equation (CQCGLE) is proposed. Different from the classic lattice Boltzmann models, this lattice Boltzmann model is based on uniformly distributed lattice points in a three-dimensional spatiotemporal space, and the evolution of the model is about a spatial axis rather than time. The algorithm possesses advantages similar to the lattice Boltzmann method in that it is easily adapted to complex Ginzburg-Landau equations. Examples show that the model reproduces the phenomena in the CQCGLE accurately.
Hybrid lattice-Boltzmann and finite-difference simulation of electroosmotic flow in a microchannel
NASA Astrophysics Data System (ADS)
Masilamani, Kannan; Ganguly, Suvankar; Feichtinger, Christian; Rüde, Ulrich
2011-04-01
A three-dimensional (3D) transient mathematical model is developed to simulate electroosmotic flows (EOFs) in a homogeneous, square cross-section microchannel, with and without considering the effects of axial pressure gradients. The general governing equations for electroosmotic transport are incompressible Navier-Stokes equations for fluid flow and the nonlinear Poisson-Boltzmann (PB) equation for electric potential distribution within the channel. In the present numerical approach, the hydrodynamic equations are solved using a lattice-Boltzmann (LB) algorithm and the PB equation is solved using a finite-difference (FD) method. The hybrid LB-FD numerical scheme is implemented on an iterative framework solving the system of coupled time-dependent partial differential equations subjected to the pertinent boundary conditions. Transient behavior of the EOF and effects due to the variations of different physicochemical parameters on the electroosmotic velocity profile are investigated. Transport characteristics for the case of combined electroosmotic- and pressure-driven microflows are also examined with the present model. For the sake of comparison, the cases of both favorable and adverse pressure gradients are considered. EOF behaviors of the non-Newtonian fluid are studied through implementation of the power-law model in the 3D LB algorithm devised for the fluid flow analysis. Numerical simulations reveal that the rheological characteristic of the fluid changes the EOF pattern to a considerable extent and can have significant consequences in the design of electroosmotically actuated bio-microfluidic systems. To improve the performance of the numerical solver, the proposed algorithm is implemented for parallel computing architectures and the overall parallel performance is found to improve with the number of processors.
Multigrid lattice Boltzmann method for accelerated solution of elliptic equations
NASA Astrophysics Data System (ADS)
Patil, Dhiraj V.; Premnath, Kannan N.; Banerjee, Sanjoy
2014-05-01
A new solver for second-order elliptic partial differential equations (PDEs) based on the lattice Boltzmann method (LBM) and the multigrid (MG) technique is presented. Several benchmark elliptic equations are solved numerically with the inclusion of multiple grid-levels in two-dimensional domains at an optimal computational cost within the LB framework. The results are compared with the corresponding analytical solutions and numerical solutions obtained using the Stone's strongly implicit procedure. The classical PDEs considered in this article include the Laplace and Poisson equations with Dirichlet boundary conditions, with the latter involving both constant and variable coefficients. A detailed analysis of solution accuracy, convergence and computational efficiency of the proposed solver is given. It is observed that the use of a high-order stencil (for smoothing) improves convergence and accuracy for an equivalent number of smoothing sweeps. The effect of the type of scheduling cycle (V- or W-cycle) on the performance of the MG-LBM is analyzed. Next, a parallel algorithm for the MG-LBM solver is presented and then its parallel performance on a multi-core cluster is analyzed. Lastly, a practical example is provided wherein the proposed elliptic PDE solver is used to compute the electro-static potential encountered in an electro-chemical cell, which demonstrates the effectiveness of this new solver in complex coupled systems. Several orders of magnitude gains in convergence and parallel scaling for the canonical problems, and a factor of 5 reduction for the multiphysics problem are achieved using the MG-LBM.
Lattice Boltzmann description of magnetization in porous media
NASA Astrophysics Data System (ADS)
Guyer, R. A.; McCall, K. R.
2000-08-01
The magnetic moments of fluid particles filling the pore space of a porous material (1) reside in a complex space, (2) are carried by the particles in diffusive exploration of the pore space, and (3) relax when the particles approach relaxation sites on the walls of the pore space. Further, when the magnetic moments are the object of a nuclear magnetic resonance experiment, they are manipulated by rf magnetic fields, internal magnetic field gradients, and applied magnetic field gradient pulses. In this paper, a lattice Boltzmann computational procedure is described that accounts for all of the vagaries in the experience of a magnetic moment in a porous material. The time evolution of the longitudinal and transverse magnetization, in a variety of experimental situations, can be simulated with this computational procedure. The z component of the magnetization, the longitudinal magnetization, is described by a set of coarse grained distribution functions for a scalar fluid. The time evolution of these distribution functions involves a scattering process (to account for diffusion) and a probability of transmission out of the pore space (to account for surface relaxation). A numerical example, involving a pore adjacent to a microporous region, is examined in detail. The transverse magnetization is a vector. Its x and y magnetization components are carried by separate scalar fluids. There is a set of coarse grained distribution functions for each fluid. Radio frequency magnetic fields, internal magnetic field gradients, applied magnetic field gradient pulses, etc., represent conversion processes in which the two fluids transform into one another. Two examples, one involving a periodic field gradient and a Hahn echo, and the other involving an isolated pore and a PFG sequence, are examined in detail.
Force method in a pseudo-potential lattice Boltzmann model
NASA Astrophysics Data System (ADS)
Hu, Anjie; Li, Longjian; Uddin, Rizwan
2015-08-01
Single component pseudo-potential lattice Boltzmann models have been widely studied due to their simplicity and stability in multiphase simulations. While numerous models have been proposed, comparative analysis and advantages and disadvantages of different force schemes are often lacking. A pseudo-potential model to simulate large density ratios proposed by Kupershtokh et al. [1] is analyzed in detail in this work. Several common used force schemes are utilized and results compared. Based on the numerical results, the relatively most accurate force scheme proposed by Guo et al. [2] is selected and applied to improve the accuracy of Kupershtokh et al.'s model. Results obtained using the modified Kupershtokh et al.'s model [1] for different value of τ are compared with those obtained using Li et al.'s model [3]. Effect of relaxation time τ on the accuracy of the results is reported. Moreover, it is noted that the error in the density ratio predicted by the model is directly correlated with the magnitude of the spurious velocities on (curved) interfaces. Simulation results show that, the accuracy of Kupershtokh et al.'s model can be improved with Guo et al.'s force scheme [2]. However, the errors and τ's effects are still noticeable when density ratios are large. To improve the accuracy of the pseudo-potential model and to reduce the effects of τ, two possible methods were discussed in the present work. Both, a rescaling of the equation of state and multi-relaxation time, are applied and are shown to improve the prediction of the density ratios.
Polar-coordinate lattice Boltzmann modeling of compressible flows
NASA Astrophysics Data System (ADS)
Lin, Chuandong; Xu, Aiguo; Zhang, Guangcai; Li, Yingjun; Succi, Sauro
2014-01-01
We present a polar coordinate lattice Boltzmann kinetic model for compressible flows. A method to recover the continuum distribution function from the discrete distribution function is indicated. Within the model, a hybrid scheme being similar to, but different from, the operator splitting is proposed. The temporal evolution is calculated analytically, and the convection term is solved via a modified Warming-Beam (MWB) scheme. Within the MWB scheme a suitable switch function is introduced. The current model works not only for subsonic flows but also for supersonic flows. It is validated and verified via the following well-known benchmark tests: (i) the rotational flow, (ii) the stable shock tube problem, (iii) the Richtmyer-Meshkov (RM) instability, and (iv) the Kelvin-Helmholtz instability. As an original application, we studied the nonequilibrium characteristics of the system around three kinds of interfaces, the shock wave, the rarefaction wave, and the material interface, for two specific cases. In one of the two cases, the material interface is initially perturbed, and consequently the RM instability occurs. It is found that the macroscopic effects due to deviating from thermodynamic equilibrium around the material interface differ significantly from those around the mechanical interfaces. The initial perturbation at the material interface enhances the coupling of molecular motions in different degrees of freedom. The amplitude of deviation from thermodynamic equilibrium around the shock wave is much higher than those around the rarefaction wave and material interface. By comparing each component of the high-order moments and its value in equilibrium, we can draw qualitatively the main behavior of the actual distribution function. These results deepen our understanding of the mechanical and material interfaces from a more fundamental level, which is indicative for constructing macroscopic models and other kinds of kinetic models.
Peristaltic particle transport using the Lattice Boltzmann method
Connington, Kevin William; Kang, Qinjun; Viswanathan, Hari S; Abdel-fattah, Amr; Chen, Shiyi
2009-01-01
Peristaltic transport refers to a class of internal fluid flows where the periodic deformation of flexible containing walls elicits a non-negligible fluid motion. It is a mechanism used to transport fluid and immersed solid particles in a tube or channel when it is ineffective or impossible to impose a favorable pressure gradient or desirous to avoid contact between the transported mixture and mechanical moving parts. Peristaltic transport occurs in many physiological situations and has myriad industrial applications. We focus our study on the peristaltic transport of a macroscopic particle in a two-dimensional channel using the lattice Boltzmann method. We systematically investigate the effect of variation of the relevant dimensionless parameters of the system on the particle transport. We find, among other results, a case where an increase in Reynolds number can actually lead to a slight increase in particle transport, and a case where, as the wall deformation increases, the motion of the particle becomes non-negative only. We examine the particle behavior when the system exhibits the peculiar phenomenon of fluid trapping. Under these circumstances, the particle may itself become trapped where it is subsequently transported at the wave speed, which is the maximum possible transport in the absence of a favorable pressure gradient. Finally, we analyze how the particle presence affects stress, pressure, and dissipation in the fluid in hopes of determining preferred working conditions for peristaltic transport of shear-sensitive particles. We find that the levels of shear stress are most hazardous near the throat of the channel. We advise that shear-sensitive particles should be transported under conditions where trapping occurs as the particle is typically situated in a region of innocuous shear stress levels.
An immersed boundary-lattice Boltzmann method for single- and multi-component fluid flows
NASA Astrophysics Data System (ADS)
Li, Zhe; Favier, Julien; D'Ortona, Umberto; Poncet, Sébastien
2016-01-01
The paper presents a numerical method to simulate single- and multi-component fluid flows around moving/deformable solid boundaries, based on the coupling of Immersed Boundary (IB) and Lattice Boltzmann (LB) methods. The fluid domain is simulated with LB method using the single relaxation time BGK model, in which an interparticle potential model is applied for multi-component fluid flows. The IB-related force is directly calculated with the interpolated definition of the fluid macroscopic velocity on the Lagrangian points that define the immersed solid boundary. The present IB-LB method can better ensure the no-slip solid boundary condition, thanks to an improved spreading operator. The proposed method is validated through several 2D/3D single- and multi-component fluid test cases with a particular emphasis on wetting conditions on solid wall. Finally, a 3D two-fluid application case is given to show the feasibility of modeling the fluid transport via a cluster of beating cilia.
Application of the lattice Boltzmann/lattice gas technique to multi-fluid flow in porous media
Soll, W.E.; Chen, S.Y.; Eggert, K.G.; Grunau, D.W.; Janecky, D.R.
1994-03-01
The lattice Boltzmann approach to modeling fluid flow provides an efficient and reliable method for solving the Navier-Stokes equations and studying multi-fluid flow problems. In this paper, we report state of the art capabilities of our lattice Boltzmann simulator for single- and two-fluid flows in two- and three-dimensional problems. We review the development of the code and present some of the latest results. Some of the flexibility available in the model includes arbitrary pore space descriptions, wettability effects, surface tension relations, and chemical reactivity. Simulations of two-fluid flow through a digitized micromodel geometry, and through a high resolution, digitized sample of Berea sandstone are presented. Relative permeability as a function of wettability and capillary number is discussed. Integration of the lattice Boltzmann approach into larger scale models to build a more powerful tool for analyzing constitutive behavior is considered.
NASA Astrophysics Data System (ADS)
Abdelaziz, Ramadan; Sussumu Komori, Fabio
2015-04-01
Recently, Lattice Boltzmann Modelling (LBM) techniques attract many scientists in various fields of research. This work shows the capability for LBM to simulate the fluid flow and solute transport in porous and fracture media, additionally, how to study behavior of nanofluids submitted to a temperature gradient, which it is an important process in natural aquatic environments, water treatment, and other water related technologies. LBSim is used in this work as Lattice Boltzmann Model simulator software. In this article, a series of cases using the lattice Boltzmann method are presented, showing the capability of the method in simulating phenomena with fluid flow and heat transfer in porous media. Results show that the lattice Boltzmann method delivers reliable and helpful simulations for the analyses of processes in water related technologies. Thus, LBSim is a recommended tool for simulating fluid flow at laminar and turbulent condition, and heat and mass transport under complex geometry and boundary condition. parameter values. Keywords: Lattice Boltzmann Model, LBSim, Fractures Media, Porous Media, nanofluids
A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows
NASA Astrophysics Data System (ADS)
Li, Kai; Zhong, Cheng-Wen
2015-05-01
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiple-relaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time (SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version. Project supported by the Innovation Fund for Aerospace Science and Technology of China (Grant No. 2009200066) and the Aeronautical Science Fund of China (Grant No. 20111453012).
NASA Technical Reports Server (NTRS)
Chau, Jessica Furrer; Or, Dani; Sukop, Michael C.; Steinberg, S. L. (Principal Investigator)
2005-01-01
Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity.
NASA Astrophysics Data System (ADS)
Biciuşcă, Tonino; Horga, Adrian; Sofonea, Victor
2015-10-01
We use a two-dimensional Lattice Boltzmann model to investigate the liquid-vapour phase separation in an isothermal van der Waals fluid. The model is based on the expansion of the distribution function up to the third order in terms of Hermite polynomials. In two dimensions, this model is an off-lattice one and has 16 velocities. The Corner Transport Upwind Scheme is used to evolve the corresponding distribution functions on a square lattice. The resulting code allows one to follow the liquid-vapour phase separation on lattices up to 4096 × 4096 nodes using a Tesla M2090 Graphics Processing Unit.
Fully-Lagrangian and Lattice-Boltzmann Methods for Solving Systems of Conservation Equations
NASA Astrophysics Data System (ADS)
Ancona, M. G.
1994-11-01
A class of "fully-Lagrangian" methods for solving systems of conservation equations is defined. The key step in formulating these methods is the definition of a new set of field variables for which Lagrangian discretization is trivial. Recently popular lattice-Boltzmann simulation schemes for solving such systems are shown to be a useful sub-class of these fully-Lagrangian methods in which (a) the conservation laws are satisfied at each grid point, (b) the Lagrangian variables are expanded perturbatively, and (c) discretization error is used to represent physics. Such schemes are typically derived using methods of kinetic theory. Our numerical analysis approach shows that the conventional physical derivation, while certainly valid and fruitful, is not essential, that it often confuses physics and numerics and that it can be unnecessarily constraining. For example, we show that lattice-Boltzmann-like methods can be non-perturbative and can be made higher-order, implicit and/or with non-uniform grids. Furthermore, our approach provides new perspective on the relationship between lattice-Boltzmann methods and finite-difference techniques. Among other things, we show that the lattice-Boltzmann schemes are only conditionally consistent and in some cases are identical to the well-known Dufort-Frankel method. Through this connection, the lattice-Boltzmann method provides a rational basis for understanding Dufort-Frankel and gives a pathway for its generalization. At the same time, that Dufort Frankel is no longer much used suggests that the lattice-Boltzmann approach might also share this fate.
Chiavazzo, Eliodoro; Karlin, Iliya V.; Gorban, Alexander N.; Boulouchos, Konstantinos
2010-10-15
A new framework of simulation of reactive flows is proposed based on a coupling between accurate reduced reaction mechanism and the lattice Boltzmann representation of the flow phenomena. The model reduction is developed in the setting of slow invariant manifold construction, and the simplest lattice Boltzmann equation is used in order to work out the procedure of coupling of the reduced model with the flow solver. Practical details of constructing slow invariant manifolds of a reaction system under various thermodynamic conditions are reported. The proposed method is validated with the two-dimensional simulation of a premixed counterflow flame in the hydrogen-air mixture. (author)
Lattice Boltzmann model for Coulomb-driven flows in dielectric liquids
NASA Astrophysics Data System (ADS)
Luo, Kang; Wu, Jian; Yi, Hong-Liang; Tan, He-Ping
2016-02-01
In this paper, we developed a unified lattice Boltzmann model (LBM) to simulate electroconvection in a dielectric liquid induced by unipolar charge injection. Instead of solving the complex set of coupled Navier-Stokes equations, the charge conservation equation, and the Poisson equation of electric potential, three consistent lattice Boltzmann equations are formulated. Numerical results are presented for both strong and weak injection regimes, and different scenarios for the onset and evolution of instability, bifurcation, and chaos are tracked. All LBM results are found to be highly consistent with the analytical solutions and other numerical work.
Lattice Boltzmann model for Coulomb-driven flows in dielectric liquids.
Luo, Kang; Wu, Jian; Yi, Hong-Liang; Tan, He-Ping
2016-02-01
In this paper, we developed a unified lattice Boltzmann model (LBM) to simulate electroconvection in a dielectric liquid induced by unipolar charge injection. Instead of solving the complex set of coupled Navier-Stokes equations, the charge conservation equation, and the Poisson equation of electric potential, three consistent lattice Boltzmann equations are formulated. Numerical results are presented for both strong and weak injection regimes, and different scenarios for the onset and evolution of instability, bifurcation, and chaos are tracked. All LBM results are found to be highly consistent with the analytical solutions and other numerical work. PMID:26986441
Lattice Boltzmann simulation of a fluid flow around a triangular unit of three isothermal cylinders
NASA Astrophysics Data System (ADS)
Alinejad, J.
2016-01-01
The lattice Boltzmann method is employed to simulate heat transfer in the flow past three arrangements of elliptical and circular cylinders under an isothermal boundary condition. The lattice Boltzmann equations and the Bhatnagar-Gross-Krook model are used to simulate two-dimensional forced convection at 30 ≤ Re ≤ 100 and Pr = 0.71. Pressure distributions, isotherms, and streamlines are obtained. Vortex shedding maps are observed in detail for several cases. The present results are in good agreement with available experimental and numerical data.
Lattice Boltzmann Hydrodynamic and Transport Modeling of Everglades Mangrove Estuaries
NASA Astrophysics Data System (ADS)
Sukop, M. C.; Engel, V.
2010-12-01
Lattice Boltzmann methods are being developed and applied to simulate groundwater and surface water flows, and heat, solute, and particle transport. Their ability to solve Navier-Stokes, St. Venant, or Darcy equations with closely coupled solute transport and density-dependent flow effects in geometrically complex domains is attractive for inverse modeling of tracer release data and forward modeling of carbon transport in mangrove estuaries under various future conditions. Key physical processes to be simulated include tidal cycles, storm surge, sea level change, variable upstream stage, subsurface groundwater inputs, and precipitation/recharge and their effects on estuary salinity and carbon transport in the estuaries and groundwater beneath the mangroves. Carbon sources and storage in the aquifer and exchanges at the mangrove-estuary interface and carbon transformations in the water column also need to be simulated. Everglades tidal mangrove estuaries are characterized by relatively high velocity (approaching 1 m s-1) tidal flows. The channels are generally less than 2 m in depth. Tidal fluctuations approach 2 m leading to significant areas of periodic inundation and emergence of oyster beds, shell beaches, mangrove root masses, and sandy beaches. Initial models are two-dimensional, although a three-dimensional model explicitly incorporating bathymetry, density-dependent flow, and wind-driven circulation could be developed. Preliminary work highlights some of the abilities of early models. A satellite image of a 64-km2 area surrounding a CO2 flux tower is used to provide the model geometry. Model resolution is 15 m per grid node. A sinusoidal tidal stage variation and constant, high salinity are applied to the Gulf side of the model while a constant stage (corresponding to mean tide), zero salinity boundary is applied on the inland side. The Navier-Stokes equations coupled with the advection-diffusion equation are solved in the open channels. The mangrove areas
Twisted 3D N=4 supersymmetric YM on deformed A{sub 3}{sup *} lattice
Saidi, El Hassan
2014-01-15
We study a class of twisted 3D N=4 supersymmetric Yang-Mills (SYM) theory on particular 3-dimensional lattice L{sub 3D} formally denoted as L{sub 3D}{sup su{sub 3}×u{sub 1}} and given by non-trivial fibration L{sub 1D}{sup u{sub 1}}×L{sub 2D}{sup su{sub 3}} with base L{sub 2D}{sup su{sub 3}}=A{sub 2}{sup *}, the weight lattice of SU(3). We first, develop the twisted 3D N=4 SYM in continuum by using superspace method where the scalar supercharge Q is manifestly exhibited. Then, we show how to engineer the 3D lattice L{sub 3D}{sup su{sub 3}×u{sub 1}} that host this theory. After that we build the lattice action S{sub latt} invariant under the following three points: (i) U(N) gauge invariance, (ii) BRST symmetry, (iii) the S{sub 3} point group symmetry of L{sub 3D}{sup su{sub 3}×u{sub 1}}. Other features such as reduction to twisted 2D supersymmetry with 8 supercharges living on L{sub 2D}≡L{sub 2D}{sup su{sub 2}×u{sub 1}}, the extension to twisted maximal 5D SYM with 16 supercharges on lattice L{sub 5D}≡L{sub 5D}{sup su{sub 4}×u{sub 1}} as well as the relation with known results are also given.
The Quantum Dynamics of a Dilute Gas in a 3D BCC Optical Lattice
NASA Astrophysics Data System (ADS)
Reichl, Linda; Boretz, Yingyue
2015-03-01
The classical and quantum dynamics of a dilute gas of rubidium atoms, in a 3D body-centered cubic optical lattice, is studied for a range of polarizations of the laser beams forming the lattice. The relative polarization of the lasers determines the the structure of the potential energy seen by the rubidium atoms. If three pairs of in-phase mutually perpendicular laser beams, with the same wavelength, form the lattice, only a limited range of possible couplings can be realized in the lab. We have determined the band structure of the BCC optical lattice for all theoretically possible couplings, and find that the band structure for lattices realizable in the lab, differs significantly from that expected for a BCC crystal. As coupling is increased, the lattice becomes increasingly chaotic and it becomes possible to produce band structure that has qualitative similarity to a BCC. Welch Foundation
Chen, Yu; Navarro, Laurent; Wang, Yan; Courbebaisse, Guy
2014-01-01
Computed Tomography Angiography (CTA) plays an essential role in the diagnosis, treatment evaluation, and monitoring of cerebral aneurysms. Segmentation of CTA medical images of giant intracranial aneurysms (GIA) provides quantitative measurements of thrombus and aneurysms geometrical characteristics allowing 3D reconstruction. In fact, GIA demonstrated neuroradiological features and propensity of partial or total spontaneous intra-aneurysmal thrombosis generating a thrombus. Despite intensive researches on medical image segmentation, aneurysm (Lumen, Thrombus, and Parent Blood Vessels) segmentation remains as a difficult problem that has not been yet resolved. In this paper, we proposed a Lattice Boltzmann Geodesic Active Contour Method (LBGM) for aneurysm segmentation in CTA images in order to estimate both the volumes of the thrombus and the aneurysm. Although the noise in the CTA images is very strong and the edges of the thrombus are not so different than the surrounding tissues, the aneurysms are segmented effectively. Based on these results, a method using a dome-neck aspect ratio (AR) parameter for the evaluation of the Spontaneous Thrombosis (ST) phenomena demonstrates the promising potentiality of this LBGM for clinical applications. PMID:24077409
Study of Gas Flow Characteristics in Tight Porous Media with a Microscale Lattice Boltzmann Model
Zhao, Jianlin; Yao, Jun; Zhang, Min; Zhang, Lei; Yang, Yongfei; Sun, Hai; An, Senyou; Li, Aifen
2016-01-01
To investigate the gas flow characteristics in tight porous media, a microscale lattice Boltzmann (LB) model with the regularization procedure is firstly adopted to simulate gas flow in three-dimensional (3D) digital rocks. A shale digital rock and a sandstone digital rock are reconstructed to study the effects of pressure, temperature and pore size on microscale gas flow. The simulation results show that because of the microscale effect in tight porous media, the apparent permeability is always higher than the intrinsic permeability, and with the decrease of pressure or pore size, or with the increase of temperature, the difference between apparent permeability and intrinsic permeability increases. In addition, the Knudsen numbers under different conditions are calculated and the results show that gas flow characteristics in the digital rocks under different Knudsen numbers are quite different. With the increase of Knudsen number, gas flow in the digital rocks becomes more uniform and the effect of heterogeneity of the porous media on gas flow decreases. Finally, two commonly used apparent permeability calculation models are evaluated by the simulation results and the Klinkenberg model shows better accuracy. In addition, a better proportionality factor in Klinkenberg model is proposed according to the simulation results. PMID:27587293
Study of Gas Flow Characteristics in Tight Porous Media with a Microscale Lattice Boltzmann Model.
Zhao, Jianlin; Yao, Jun; Zhang, Min; Zhang, Lei; Yang, Yongfei; Sun, Hai; An, Senyou; Li, Aifen
2016-01-01
To investigate the gas flow characteristics in tight porous media, a microscale lattice Boltzmann (LB) model with the regularization procedure is firstly adopted to simulate gas flow in three-dimensional (3D) digital rocks. A shale digital rock and a sandstone digital rock are reconstructed to study the effects of pressure, temperature and pore size on microscale gas flow. The simulation results show that because of the microscale effect in tight porous media, the apparent permeability is always higher than the intrinsic permeability, and with the decrease of pressure or pore size, or with the increase of temperature, the difference between apparent permeability and intrinsic permeability increases. In addition, the Knudsen numbers under different conditions are calculated and the results show that gas flow characteristics in the digital rocks under different Knudsen numbers are quite different. With the increase of Knudsen number, gas flow in the digital rocks becomes more uniform and the effect of heterogeneity of the porous media on gas flow decreases. Finally, two commonly used apparent permeability calculation models are evaluated by the simulation results and the Klinkenberg model shows better accuracy. In addition, a better proportionality factor in Klinkenberg model is proposed according to the simulation results. PMID:27587293
Lattice Boltzmann Method of Different BGA Orientations on I-Type Dispensing Method
Gan, Z. L.; Ishak, M. H. H.; Abdullah, M. Z.; Khor, Soon Fuat
2016-01-01
This paper studies the three dimensional (3D) simulation of fluid flows through the ball grid array (BGA) to replicate the real underfill encapsulation process. The effect of different solder bump arrangements of BGA on the flow front, pressure and velocity of the fluid is investigated. The flow front, pressure and velocity for different time intervals are determined and analyzed for potential problems relating to solder bump damage. The simulation results from Lattice Boltzmann Method (LBM) code will be validated with experimental findings as well as the conventional Finite Volume Method (FVM) code to ensure highly accurate simulation setup. Based on the findings, good agreement can be seen between LBM and FVM simulations as well as the experimental observations. It was shown that only LBM is capable of capturing the micro-voids formation. This study also shows an increasing trend in fluid filling time for BGA with perimeter, middle empty and full orientations. The perimeter orientation has a higher pressure fluid at the middle region of BGA surface compared to middle empty and full orientation. This research would shed new light for a highly accurate simulation of encapsulation process using LBM and help to further increase the reliability of the package produced. PMID:27454872
Lattice Boltzmann Method of Different BGA Orientations on I-Type Dispensing Method.
Abas, Aizat; Gan, Z L; Ishak, M H H; Abdullah, M Z; Khor, Soon Fuat
2016-01-01
This paper studies the three dimensional (3D) simulation of fluid flows through the ball grid array (BGA) to replicate the real underfill encapsulation process. The effect of different solder bump arrangements of BGA on the flow front, pressure and velocity of the fluid is investigated. The flow front, pressure and velocity for different time intervals are determined and analyzed for potential problems relating to solder bump damage. The simulation results from Lattice Boltzmann Method (LBM) code will be validated with experimental findings as well as the conventional Finite Volume Method (FVM) code to ensure highly accurate simulation setup. Based on the findings, good agreement can be seen between LBM and FVM simulations as well as the experimental observations. It was shown that only LBM is capable of capturing the micro-voids formation. This study also shows an increasing trend in fluid filling time for BGA with perimeter, middle empty and full orientations. The perimeter orientation has a higher pressure fluid at the middle region of BGA surface compared to middle empty and full orientation. This research would shed new light for a highly accurate simulation of encapsulation process using LBM and help to further increase the reliability of the package produced. PMID:27454872
Lattice Boltzmann simulation of asymmetric flow in nematic liquid crystals with finite anchoring
NASA Astrophysics Data System (ADS)
Zhang, Rui; Roberts, Tyler; Aranson, Igor S.; de Pablo, Juan J.
2016-02-01
Liquid crystals (LCs) display many of the flow characteristics of liquids but exhibit long range orientational order. In the nematic phase, the coupling of structure and flow leads to complex hydrodynamic effects that remain to be fully elucidated. Here, we consider the hydrodynamics of a nematic LC in a hybrid cell, where opposite walls have conflicting anchoring boundary conditions, and we employ a 3D lattice Boltzmann method to simulate the time-dependent flow patterns that can arise. Due to the symmetry breaking of the director field within the hybrid cell, we observe that at low to moderate shear rates, the volumetric flow rate under Couette and Poiseuille flows is different for opposite flow directions. At high shear rates, the director field may undergo a topological transition which leads to symmetric flows. By applying an oscillatory pressure gradient to the channel, a net volumetric flow rate is found to depend on the magnitude and frequency of the oscillation, as well as the anchoring strength. Taken together, our findings suggest several intriguing new applications for LCs in microfluidic devices.
Lattice Boltzmann simulation of asymmetric flow in nematic liquid crystals with finite anchoring.
Zhang, Rui; Roberts, Tyler; Aranson, Igor S; de Pablo, Juan J
2016-02-28
Liquid crystals (LCs) display many of the flow characteristics of liquids but exhibit long range orientational order. In the nematic phase, the coupling of structure and flow leads to complex hydrodynamic effects that remain to be fully elucidated. Here, we consider the hydrodynamics of a nematic LC in a hybrid cell, where opposite walls have conflicting anchoring boundary conditions, and we employ a 3D lattice Boltzmann method to simulate the time-dependent flow patterns that can arise. Due to the symmetry breaking of the director field within the hybrid cell, we observe that at low to moderate shear rates, the volumetric flow rate under Couette and Poiseuille flows is different for opposite flow directions. At high shear rates, the director field may undergo a topological transition which leads to symmetric flows. By applying an oscillatory pressure gradient to the channel, a net volumetric flow rate is found to depend on the magnitude and frequency of the oscillation, as well as the anchoring strength. Taken together, our findings suggest several intriguing new applications for LCs in microfluidic devices. PMID:26931724
Lattice-Boltzmann simulations of three-dimensional fluid flow on a desktop computer.
Brewster, Jeffrey D
2007-04-01
The lattice-Boltzmann (LB) method is a cellular automaton approach to simulating fluid flow with many advantages over conventional methods based on the Navier-Stokes equations. It is conceptually simple, amenable to a wide array of boundary conditions, and can be adapted to handle thermal, density, miscibility, and other effects. The LB approach has been used to model a number of fluid systems of interest to analytical chemists, including chromatography columns, micromixers, and electroosmotic pumps. However, widespread use of this tool has been limited, in part because virtually all large-scale 3D simulations in the literature have been executed on supercomputers. This work demonstrates that such simulations can be executed in reasonable periods of time (hours) on a desktop computer using a cross-platform software package that is easy to learn and use. This package incorporates several improvements that enhance the utility of the LB approach, including an algorithm for speeding common calculations by 2 orders of magnitude and a scheme for handling convection-diffusion equations of interest in electrochemical and surface reaction studies. PMID:17319648
Extending a CAD-Based Cartesian Mesh Generator for the Lattice Boltzmann Method
Cantrell, J Nathan; Inclan, Eric J; Joshi, Abhijit S; Popov, Emilian L; Jain, Prashant K
2012-01-01
This paper describes the development of a custom preprocessor for the PaRAllel Thermal Hydraulics simulations using Advanced Mesoscopic methods (PRATHAM) code based on an open-source mesh generator, CartGen [1]. PRATHAM is a three-dimensional (3D) lattice Boltzmann method (LBM) based parallel flow simulation software currently under development at the Oak Ridge National Laboratory. The LBM algorithm in PRATHAM requires a uniform, coordinate system-aligned, non-body-fitted structured mesh for its computational domain. CartGen [1], which is a GNU-licensed open source code, already comes with some of the above needed functionalities. However, it needs to be further extended to fully support the LBM specific preprocessing requirements. Therefore, CartGen is being modified to (i) be compiler independent while converting a neutral-format STL (Stereolithography) CAD geometry to a uniform structured Cartesian mesh, (ii) provide a mechanism for PRATHAM to import the mesh and identify the fluid/solid domains, and (iii) provide a mechanism to visually identify and tag the domain boundaries on which to apply different boundary conditions.
NASA Astrophysics Data System (ADS)
Geneva, Nicholas; Wang, Lian-Ping
2015-11-01
In the past 25 years, the mesoscopic lattice Boltzmann method (LBM) has become an increasingly popular approach to simulate incompressible flows including turbulent flows. While LBM solves more solution variables compared to the conventional CFD approach based on the macroscopic Navier-Stokes equation, it also offers opportunities for more efficient parallelization. In this talk we will describe several different algorithms that have been developed over the past 10 plus years, which can be used to represent the two core steps of LBM, collision and streaming, more effectively than standard approaches. The application of these algorithms spans LBM simulations ranging from basic channel to particle laden flows. We will cover the essential detail on the implementation of each algorithm for simple 2D flows, to the challenges one faces when using a given algorithm for more complex simulations. The key is to explore the best use of data structure and cache memory. Two basic data structures will be discussed and the importance of effective data storage to maximize a CPU's cache will be addressed. The performance of a 3D turbulent channel flow simulation using these different algorithms and data structures will be compared along with important hardware related issues.
Increasing the filling of ultracold KRb molecules in a 3D optical lattice
NASA Astrophysics Data System (ADS)
Moses, Steven; Covey, Jacob; Gadway, Bryce; Yan, Bo; Miecnikowski, Matthew; Ye, Jun; Jin, Deborah
2015-05-01
Ultracold polar molecules, with their long-range electric dipolar interactions, offer new opportunities for studying quantum magnetism and many-body physics. Recently, our group observed spin exchange interactions between KRb molecules in a 3D optical lattice, which is one of the first steps towards studying lattice spin models with polar molecules. The lattice fillings were about 10% or less in these experiments. Future experiments will benefit greatly from lower entropies and higher lattice fillings. Here, we have investigated the molecular creation process in a 3D optical lattice with the goal of maximizing the filling fraction. We start by loading a BEC of Rb and a degenerate Fermi gas of K into a 3D optical lattice. In the absence of K, Rb is a Mott insulator. We study how the Mott insulator and the filling of Rb are affected by the presence of K and develop a strategy to maintain high Rb filling throughout the molecule production process. We also find that we can convert a large fraction of these Rb to molecules when we operate with low Rb numbers. We acknowledge funding from DARPA, DOE, NIST, NSF, AFOSR, and the NDSEG Graduate Fellowship.
Boltzmann Fluctuations in Numerical Simulations of Nonequilibrium Lattice Threshold Systems
Rundle, J.B.; Klein, W.; Gross, S.; Turcotte, D.L.
1995-08-21
Nonequilibrium threshold systems such as slider blocks are now used to model a variety of dynamical systems, including earthquake faults, driven neural networks, and sliding charge density waves. We show that for general mean field models driven at low rates fluctuations in the internal energy field are characterized by Boltzmann statistics. Numerical simulations confirm this prediction. Our results indicate that mean field models can be effectively treated as equilibrium systems.
McGhee, J.M.; Roberts, R.M.; Morel, J.E.
1997-06-01
A spherical harmonics research code (DANTE) has been developed which is compatible with parallel computer architectures. DANTE provides 3-D, multi-material, deterministic, transport capabilities using an arbitrary finite element mesh. The linearized Boltzmann transport equation is solved in a second order self-adjoint form utilizing a Galerkin finite element spatial differencing scheme. The core solver utilizes a preconditioned conjugate gradient algorithm. Other distinguishing features of the code include options for discrete-ordinates and simplified spherical harmonics angular differencing, an exact Marshak boundary treatment for arbitrarily oriented boundary faces, in-line matrix construction techniques to minimize memory consumption, and an effective diffusion based preconditioner for scattering dominated problems. Algorithm efficiency is demonstrated for a massively parallel SIMD architecture (CM-5), and compatibility with MPP multiprocessor platforms or workstation clusters is anticipated.
Pore-scale discretisation limits of multiphase lattice-Boltzmann methods
NASA Astrophysics Data System (ADS)
Li, Z.; Middleton, J.; Varslot, T.; Sheppard, A.
2015-12-01
Lattice-Boltzmann (LB) modeling is a popular method for the numerical solution of the Navier-Stokes equations and several multi-component LB models are widely used to simulate immiscible two-phase fluid flow in porous media. However, there has been relatively little study of the models' ability to make optimal use of 3D imagery by considering the minimum number of grid points that are needed to represent geometric features such as pore throats. This is of critical importance since 3D images of geological samples are a compromise between resolution and field of view. In this work we explore the discretisation limits of LB models, their behavior near these limits, and the consequences of this behavior for simulations of drainage and imbibition. We quantify the performance of two commonly used multiphase LB models: Shan-Chen (SC) and Rothman-Keller (RK) models in a set of tests, including simulations of bubbles in bulk fluid, on flat surfaces, confined in flat/tilted tubes, and fluid invasion into single tubes. Simple geometries like these allow better quantification of model behavior and better understanding of breakdown mechanisms. In bulk fluid, bubble radii less than 2.5 grid units (image voxels) cause numerical instability in SC model; the RK model is stable to a radius of 2.5 units and below, but with poor agreement with the Laplace's law. When confined to a flat duct, the SC model can simulate similar radii to RK model, but with higher interface spurious currents than the RK model and some risk of instability. In tilted ducts with 'staircase' voxel-level roughness, the SC model seems to average the roughness, whereas for RK model only the 'peaks' of the surface are relevant. Overall, our results suggest that LB models can simulate fluid capillary pressure corresponding to interfacial radii of just 1.5 grid units, with the RK model exhibiting significantly better stability.
Lattice Boltzmann simulation of three-dimensional Rayleigh-Taylor instability
NASA Astrophysics Data System (ADS)
Liang, H.; Li, Q. X.; Shi, B. C.; Chai, Z. H.
2016-03-01
In this paper, the three-dimensional (3D) Rayleigh-Taylor instability (RTI) with low Atwood number (At=0.15 ) in a long square duct (12 W ×W ×W ) is studied by using a multiple-relaxation-time lattice Boltzmann (LB) multiphase model. The effect of the Reynolds number on the interfacial dynamics and bubble and spike amplitudes at late time is investigated in detail. The numerical results show that at sufficiently large Reynolds numbers, a sequence of stages in the 3D immiscible RTI can be observed, which includes the linear growth, terminal velocity growth, reacceleration, and chaotic development stages. At late stage, the RTI induces a very complicated topology structure of the interface, and an abundance of dissociative drops are also observed in the system. The bubble and spike velocities at late stage are unstable and their values have exceeded the predictions of the potential flow theory [V. N. Goncharov, Phys. Rev. Lett. 88, 134502 (2002), 10.1103/PhysRevLett.88.134502]. The acceleration of the bubble front is also measured and it is found that the normalized acceleration at late time fluctuates around a constant value of 0.16. When the Reynolds number is reduced to small values, some later stages cannot be reached sequentially. The interface becomes relatively smoothed and the bubble velocity at late time is approximate to a constant value, which coincides with the results of the extended Layzer model [S.-I. Sohn, Phys. Rev. E 80, 055302(R) (2009), 10.1103/PhysRevE.80.055302] and the modified potential theory [R. Banerjee, L. Mandal, S. Roy, M. Khan, and M. R. Guptae, Phys. Plasmas 18, 022109 (2011), 10.1063/1.3555523]. In our simulations, the Graphics Processing Unit (GPU) parallel computing is also used to relieve the massive computational cost.
Lattice Boltzmann simulation of three-dimensional Rayleigh-Taylor instability.
Liang, H; Li, Q X; Shi, B C; Chai, Z H
2016-03-01
In this paper, the three-dimensional (3D) Rayleigh-Taylor instability (RTI) with low Atwood number (A(t)=0.15) in a long square duct (12W × W × W) is studied by using a multiple-relaxation-time lattice Boltzmann (LB) multiphase model. The effect of the Reynolds number on the interfacial dynamics and bubble and spike amplitudes at late time is investigated in detail. The numerical results show that at sufficiently large Reynolds numbers, a sequence of stages in the 3D immiscible RTI can be observed, which includes the linear growth, terminal velocity growth, reacceleration, and chaotic development stages. At late stage, the RTI induces a very complicated topology structure of the interface, and an abundance of dissociative drops are also observed in the system. The bubble and spike velocities at late stage are unstable and their values have exceeded the predictions of the potential flow theory [V. N. Goncharov, Phys. Rev. Lett. 88, 134502 (2002)]. The acceleration of the bubble front is also measured and it is found that the normalized acceleration at late time fluctuates around a constant value of 0.16. When the Reynolds number is reduced to small values, some later stages cannot be reached sequentially. The interface becomes relatively smoothed and the bubble velocity at late time is approximate to a constant value, which coincides with the results of the extended Layzer model [S.-I. Sohn, Phys. Rev. E 80, 055302(R) (2009)] and the modified potential theory [R. Banerjee, L. Mandal, S. Roy, M. Khan, and M. R. Guptae, Phys. Plasmas 18, 022109 (2011)]. In our simulations, the Graphics Processing Unit (GPU) parallel computing is also used to relieve the massive computational cost. PMID:27078453
NASA Astrophysics Data System (ADS)
Brogi, F.; Malaspinas, O.; Bonadonna, C.; Chopard, B.; Ripepe, M.
2015-12-01
Low frequency (< 20Hz) acoustic measurements have a great potential for the real time characterization of volcanic plume source parameters. Using the classical source theory, acoustic data can be related to the exit velocity of the volcanic jet and to mass eruption rate, based on the geometric constrain of the vent and the mixture density. However, the application of the classical acoustic source models to volcanic explosive eruptions has shown to be challenging and a better knowledge of the link between the acoustic radiation and actual volcanic fluid dynamics processes is required. New insights into this subject could be given by the study of realistic aeroacoustic numerical simulations of a volcanic jet. Lattice Boltzmann strategies (LBS) provide the opportunity to develop an accurate, computationally fast, 3D physical model for a volcanic jet. In the field of aeroacoustic applications, dedicated LBS has been proven to have the low dissipative properties needed for capturing the weak acoustic pressure fluctuations. However, due to the big disparity in magnitude between the flow and the acoustic disturbances, even weak spurious noise sources in simulations can ruin the accuracy of the acoustic predictions. Reflected waves from artificial boundaries defined around the flow region can have significant influence on the flow field and overwhelm the acoustic field of interest. In addition, for highly multiscale turbulent flows, such as volcanic plumes, the number of grid points needed to represent the smallest scales might become intractable and the most complicated physics happen only in small portions of the computational domain. The implementation of the grid refinement, in our model allow us to insert local finer grids only where is actually needed and to increase the size of the computational domain for running more realistic simulations. 3D LBS model simulations for turbulent jet aeroacoustics have been accurately validated. Both mean flow and acoustic results
The lattice Boltzmann model for the second-order Benjamin-Ono equations
NASA Astrophysics Data System (ADS)
Lai, Huilin; Ma, Changfeng
2010-04-01
In this paper, in order to extend the lattice Boltzmann method to deal with more complicated nonlinear equations, we propose a 1D lattice Boltzmann scheme with an amending function for the second-order (1 + 1)-dimensional Benjamin-Ono equation. With the Taylor expansion and the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The equilibrium distribution function and the amending function are obtained. Numerical simulations are carried out for the 'good' Boussinesq equation and the 'bad' one to validate the proposed model. It is found that the numerical results agree well with the analytical solutions. The present model can be used to solve more kinds of nonlinear partial differential equations.
Topology optimization in thermal-fluid flow using the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Yaji, Kentaro; Yamada, Takayuki; Yoshino, Masato; Matsumoto, Toshiro; Izui, Kazuhiro; Nishiwaki, Shinji
2016-02-01
This paper proposes a topology optimization method for thermal-fluid flow problems using the lattice Boltzmann method (LBM). The design sensitivities are derived based on the adjoint lattice Boltzmann method (ALBM), whose basic idea is that the adjoint problem is first formulated using a continuous adjoint approach, and the adjoint problem is then solved using the LBM. In this paper, the discrete velocity Boltzmann equation, in which only the particle velocities are discretized, is introduced to the ALBM to deal with the various boundary conditions in the LBM. The novel sensitivity analysis is applied in two flow channel topology optimization problems: 1) a pressure drop minimization problem, and 2) a heat exchange maximization problem. Several numerical examples are provided to confirm the utility of the proposed method.
Developing extensible lattice-Boltzmann simulationsfor general-purpose graphics-programming units
Walsh, S C; Saar, M O
2011-10-27
Lattice-Boltzmann methods are versatile numerical modeling techniques capable of reproducing a wide variety of fluid-mechanical behavior. These methods are well suited to parallel implementation, particularly on the single-instruction multiple data (SIMD) parallel processing environments found in computer graphics processing units (GPUs). Although more recent programming tools dramatically improve the ease with which GPU programs can be written, the programming environment still lacks the flexibility available to more traditional CPU programs. In particular, it may be difficult to develop modular and extensible programs that require variable on-device functionality with current GPU architectures. This paper describes a process of automatic code generation that overcomes these difficulties for lattice-Boltzmann simulations. It details the development of GPU-based modules for an extensible lattice-Boltzmann simulation package - LBHydra. The performance of the automatically generated code is compared to equivalent purpose written codes for both single-phase, multiple-phase, and multiple-component flows. The flexibility of the new method is demonstrated by simulating a rising, dissolving droplet in a porous medium with user generated lattice-Boltzmann models and subroutines.
Developing extensible lattice-Boltzmann simulators for general-purpose graphics-processing units
Walsh, S C; Saar, M O
2011-12-21
Lattice-Boltzmann methods are versatile numerical modeling techniques capable of reproducing a wide variety of fluid-mechanical behavior. These methods are well suited to parallel implementation, particularly on the single-instruction multiple data (SIMD) parallel processing environments found in computer graphics processing units (GPUs). Although more recent programming tools dramatically improve the ease with which GPU programs can be written, the programming environment still lacks the flexibility available to more traditional CPU programs. In particular, it may be difficult to develop modular and extensible programs that require variable on-device functionality with current GPU architectures. This paper describes a process of automatic code generation that overcomes these difficulties for lattice-Boltzmann simulations. It details the development of GPU-based modules for an extensible lattice-Boltzmann simulation package - LBHydra. The performance of the automatically generated code is compared to equivalent purpose written codes for both single-phase, multiple-phase, and multiple-component flows. The flexibility of the new method is demonstrated by simulating a rising, dissolving droplet in a porous medium with user generated lattice-Boltzmann models and subroutines.
Reassessing the single relaxation time Lattice Boltzmann method for the simulation of Darcy’s flows
NASA Astrophysics Data System (ADS)
Prestininzi, Pietro; Montessori, Andrea; La Rocca, Michele; Succi, Sauro
2016-09-01
It is shown that the single relaxation time (SRT) version of the Lattice Boltzmann (LB) equation permits to compute the permeability of Darcy’s flows in porous media within a few percent accuracy. This stands in contrast with previous claims of inaccuracy, which we relate to the lack of recognition of the physical dependence of the permeability on the Knudsen number.
Revised lattice Boltzmann model for traffic flow with equilibrium traffic pressure
NASA Astrophysics Data System (ADS)
Shi, Wei; Lu, Wei-Zhen; Xue, Yu; He, Hong-Di
2016-02-01
A revised lattice Boltzmann model concerning the equilibrium traffic pressure is proposed in this study to tackle the phase transition phenomena of traffic flow system. The traditional lattice Boltzmann model has limitation to investigate the complex traffic phase transitions due to its difficulty for modeling the equilibrium velocity distribution. Concerning this drawback, the equilibrium traffic pressure is taken into account to derive the equilibrium velocity distribution in the revised lattice Boltzmann model. In the proposed model, a three-dimensional velocity-space is assumed to determine the equilibrium velocity distribution functions and an alternative, new derivative approach is introduced to deduct the macroscopic equations with the first-order accuracy level from the lattice Boltzmann model. Based on the linear stability theory, the stability conditions of the corresponding macroscopic equations can be obtained. The outputs indicate that the stability curve is divided into three regions, i.e., the stable region, the neutral stability region, and the unstable region. In the stable region, small disturbance appears in the initial uniform flow and will vanish after long term evolution, while in the unstable region, the disturbance will be enlarged and finally leads to the traffic system entering the congested state. In the neutral stability region, small disturbance does not vanish with time and maintains its amplitude in the traffic system. Conclusively, the stability of traffic system is found to be enhanced as the equilibrium traffic pressure increases. Finally, the numerical outputs of the proposed model are found to be consistent with the recognized, theoretical results.
NASA Technical Reports Server (NTRS)
Menon, Suresh
2003-01-01
This report summarizes the progress made in the first 8 to 9 months of this research. The Lattice Boltzmann Equation (LBE) methodology for Large-eddy Simulations (LES) of microblowing has been validated using a jet-in-crossflow test configuration. In this study, the flow intake is also simulated to allow the interaction to occur naturally. The Lattice Boltzmann Equation Large-eddy Simulations (LBELES) approach is capable of capturing not only the flow features associated with the flow, such as hairpin vortices and recirculation behind the jet, but also is able to show better agreement with experiments when compared to previous RANS predictions. The LBELES is shown to be computationally very efficient and therefore, a viable method for simulating the injection process. Two strategies have been developed to simulate multi-hole injection process as in the experiment. In order to allow natural interaction between the injected fluid and the primary stream, the flow intakes for all the holes have to be simulated. The LBE method is computationally efficient but is still 3D in nature and therefore, there may be some computational penalty. In order to study a large number or holes, a new 1D subgrid model has been developed that will simulate a reduced form of the Navier-Stokes equation in these holes.
NASA Astrophysics Data System (ADS)
Wong, T.; Sun, W.
2012-12-01
Microcomputed tomography can be used to characterize the geometry of the pore space of a sedimentary rock, with resolution that is sufficiently refined for the realistic simulation of physical properties based on the 3D image. Significant advances have been made on the characterization of pore size distribution and connectivity, development of techniques such as lattice Boltzmann method to simulate permeability, and its upscaling. Sun, Andrade and Rudnicki (2011) recently introduced a multiscale method that dynamically links these three aspects, which were often treated separately in previous computational schemes. In this study, we improve the efficiency of this multiscale method by introducing a flood-fill algorithm to determine connectivity of the pores, followed by a multiscale lattice Boltzmann/finite element calculation to obtain homogenized effective anisotropic permeability. The improved multiscale method also includes new capacity to consistently determine electrical conductivity and formation factor from CT images. Furthermore, we also introduce a level set based method that transforms pore geometry to finite element mesh and thus enables direct simulation of pore-scale flow with finite element method. When applied to the microCT data acquired by Lindquist et al. (2000) for four Fontainebleau sandstone samples with porosities ranging from 7.5% to 22%, this multiscale method has proved to be computationally efficient and our simulations has provided new insights into the relation among permeability, pore geometry and connectivity.
Loading mode dependent effective properties of octet-truss lattice structures using 3D-printing
NASA Astrophysics Data System (ADS)
Challapalli, Adithya
Cellular materials, often called lattice materials, are increasingly receiving attention for their ultralight structures with high specific strength, excellent impact absorption, acoustic insulation, heat dissipation media and compact heat exchangers. In alignment with emerging additive manufacturing (AM) technology, realization of the structural applications of the lattice materials appears to be becoming faster. Considering the direction dependent material properties of the products with AM, by directionally dependent printing resolution, effective moduli of lattice structures appear to be directionally dependent. In this paper, a constitutive model of a lattice structure, which is an octet-truss with a base material having an orthotropic material property considering AM is developed. In a case study, polyjet based 3D printing material having an orthotropic property with a 9% difference in the principal direction provides difference in the axial and shear moduli in the octet-truss by 2.3 and 4.6%. Experimental validation for the effective properties of a 3D printed octet-truss is done for uniaxial tension and compression test. The theoretical value based on the micro-buckling of truss member are used to estimate the failure strength. Modulus value appears a little overestimate compared with the experiment. Finite element (FE) simulations for uniaxial compression and tension of octettruss lattice materials are conducted. New effective properties for the octet-truss lattice structure are developed considering the observed behavior of the octet-truss structure under macroscopic compression and tension trough simulations.
Increasing the filling fraction of ultracold polar KRb molecules in a 3D optical lattice
NASA Astrophysics Data System (ADS)
Moses, Steven; Gadway, Bryce; Yan, Bo; Covey, Jacob; Jin, Deborah; Ye, Jun
2013-05-01
Gases of ultracold polar molecules with sufficiently low entropy are an ideal experimental scenario to look for signatures of long-range many-body interactions. Having a high filling fraction in a 3D lattice is one way to achieve a low entropy state. In earlier work, we showed that preformed pairs of K and Rb in a 3D lattice (sites that have exactly one K and one Rb) are converted to KRb Feshbach molecules with nearly 100% efficiency. Since the STIRAP transfer from Feshbach molecules to ground-state molecules is 90-100% efficient, loading a 3D lattice with a large fraction of preformed pairs will lead to a large filling fraction of ground-state molecules. Our scheme is to load a Mott insulator of Rb and band insulator of K. After we have loaded a lattice with a high filling fraction, we will investigate effects of dipolar interactions in a many-body system. We acknowledge funding from NIST, NSF, AFOSR-MURI, and the NDSEG Graduate Fellowship.
Lattice-Boltzmann simulation of laser interaction with weakly ionized helium plasmas
Li Huayu; Ki, Hyungson
2010-07-15
This paper presents a lattice Boltzmann method for laser interaction with weakly ionized plasmas considering electron impact ionization and three-body recombination. To simulate with physical properties of plasmas, the authors' previous work on the rescaling of variables is employed and the electromagnetic fields are calculated from the Maxwell equations by using the finite-difference time-domain method. To calculate temperature fields, energy equations are derived separately from the Boltzmann equations. In this way, we attempt to solve the full governing equations for plasma dynamics. With the developed model, the continuous-wave CO{sub 2} laser interaction with helium is simulated successfully.
NASA Astrophysics Data System (ADS)
Chen, Lin-Jie; Ma, Chang-Feng
2010-01-01
This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut + αuux + βunux + γuxx + δuxxx + ζuxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.
Lin, X.
1991-01-01
This paper reports the development of an object-oriented programming methodology for particle simulations. It is established on the (m reductionist) view that many physical phenomena cana be reduced to many-body problems. By doing the reduction, many seemly unrelated physical phenomena can be simulated in a systematic way and a high-level programming system can be constructed to facilitate the programming and the solution of the simulations. In the object-oriented particle simulation methodology, a hierarchy of abstract particles is defined to represent a variety of characteristics in physical system simulations. A simulation program is constructed from particles derived from the abstract particles. The object- oriented particle simulation methodology provides a unifying modeling and simulation framework for a variety of simulation applications with the use of particle methods. It allows easy composition of simulation programs from predefined software modules and facilitates software reusability. It greatly increase the productivity of simulation program constructions. Boltzmann (after Ludwig Boltzmann, 1844-1906) is a prototype programming system in the object-oriented particle simulation methodology. Boltzmann is implemented in C++ and the X Window System. It contains a library of data types and functions that support simulations in particle methods. Moreover, it provides a visualization window to support friendly user-computer interaction. Examples of the application of the Boltzmann programming system are presented. The effectiveness of the object-oriented particle simulation methodology is demonstrated. A user's manual is included in the appendix.
Finite-difference lattice Boltzmann simulation on acoustics-induced particle deposition
NASA Astrophysics Data System (ADS)
Fu, Sau-Chung; Yuen, Wai-Tung; Wu, Chili; Chao, Christopher Yu-Hang
2015-10-01
Particle manipulation by acoustics has been investigated for many years. By a proper design, particle deposition can be induced by the same principle. The use of acoustics can potentially be developed into an energy-efficient technique for particle removal or filtration system as the pressure drop due to acoustic effects is low and the flow velocity is not necessary to be high. Two nonlinear acoustic effects, acoustic streaming and acoustic radiation pressure, are important. Acoustic streaming introduces vortices and stagnation points on the surface of an air duct and removes the particles by deposition. Acoustic radiation pressure causes particles to form agglomerates and enhances inertial impaction and/or gravitational sedimentation. The objective of this paper is to develop a numerical model to investigate the particle deposition induced by acoustic effects. A three-step approach is adopted and lattice Boltzamnn technique is employed as the numerical method. This is because the lattice Boltzmann equation is hyperbolic and can be solved locally, explicitly, and efficiently on parallel computers. In the first step, the acoustic field and its mean square fluctuation values are calculated. Due to the advantage of the lattice Boltzmann technique, a simple, stable and fast lattice Boltzmann method is proposed and verified. The result of the first step is input into the second step to solve for acoustic streaming. Another finite difference lattice Boltzmann method, which has been validated by a number of flows and benchmark cases in the literature, is used. The third step consists in tracking the particle's motion by a Lagrangian approach where the acoustic radiation pressure is considered. The influence of the acoustics effects on particle deposition is explained. The numerical result matches with an experiment. The model is a useful tool for optimizing the design and helps to further develop the technique.
Lattice-Boltzmann simulation of coalescence-driven island coarsening
Basagaoglu, H.; Green, C.T.; Meakin, P.; McCoy, B.J.
2004-01-01
The first-order phase separation in a thin fluid film was simulated using a two-dimensional lattice-Boltzman model (LBM) with fluid-fluid interactions. The effects of the domain size on the intermediate asymptotic island size distribution were also discussed. It was observed that the overall process is dominated by coalescence which is independent of island mass. The results show that the combined effects of growth, coalescence, and Ostwald ripening control the phase transition process in the LBM simulations.
Enhanced hybrid search algorithm for protein structure prediction using the 3D-HP lattice model.
Zhou, Changjun; Hou, Caixia; Zhang, Qiang; Wei, Xiaopeng
2013-09-01
The problem of protein structure prediction in the hydrophobic-polar (HP) lattice model is the prediction of protein tertiary structure. This problem is usually referred to as the protein folding problem. This paper presents a method for the application of an enhanced hybrid search algorithm to the problem of protein folding prediction, using the three dimensional (3D) HP lattice model. The enhanced hybrid search algorithm is a combination of the particle swarm optimizer (PSO) and tabu search (TS) algorithms. Since the PSO algorithm entraps local minimum in later evolution extremely easily, we combined PSO with the TS algorithm, which has properties of global optimization. Since the technologies of crossover and mutation are applied many times to PSO and TS algorithms, so enhanced hybrid search algorithm is called the MCMPSO-TS (multiple crossover and mutation PSO-TS) algorithm. Experimental results show that the MCMPSO-TS algorithm can find the best solutions so far for the listed benchmarks, which will help comparison with any future paper approach. Moreover, real protein sequences and Fibonacci sequences are verified in the 3D HP lattice model for the first time. Compared with the previous evolutionary algorithms, the new hybrid search algorithm is novel, and can be used effectively to predict 3D protein folding structure. With continuous development and changes in amino acids sequences, the new algorithm will also make a contribution to the study of new protein sequences. PMID:23824509
Flow force and torque on submerged bodies in lattice-Boltzmann methods via momentum exchange.
Giovacchini, Juan P; Ortiz, Omar E
2015-12-01
We review the momentum exchange method to compute the flow force and torque on a submerged body in lattice-Boltzmann methods by presenting an alternative derivation. Our derivation does not depend on a particular implementation of the boundary conditions at the body surface, and it relies on general principles. After the introduction of the momentum exchange method in lattice-Boltzmann methods, some formulations were introduced to compute the fluid force on static and moving bodies. These formulations were introduced in a rather intuitive, ad hoc way. In our derivation, we recover the proposals most frequently used, in some cases with minor corrections, gaining some insight into the two most used formulations. At the end, we present some numerical tests to compare different approaches on a well-known benchmark test that support the correctness of the formulas derived. PMID:26764848
Conservative phase-field lattice Boltzmann model for interface tracking equation.
Geier, Martin; Fakhari, Abbas; Lee, Taehun
2015-06-01
Based on the phase-field theory, we propose a conservative lattice Boltzmann method to track the interface between two different fluids. The presented model recovers the conservative phase-field equation and conserves mass locally and globally. Two entirely different approaches are used to calculate the gradient of the phase field, which is needed in computation of the normal to the interface. One approach uses finite-difference stencils similar to many existing lattice Boltzmann models for tracking the two-phase interface, while the other one invokes central moments to calculate the gradient of the phase field without any finite differences involved. The former approach suffers from the nonlocality of the collision operator while the latter is entirely local making it highly suitable for massive parallel implementation. Several benchmark problems are carried out to assess the accuracy and stability of the proposed model. PMID:26172824
Theory of the Lattice Boltzmann Equation: Symmetry properties of Discrete Velocity Sets
NASA Technical Reports Server (NTRS)
Rubinstein, Robert; Luo, Li-Shi
2007-01-01
In the lattice Boltzmann equation, continuous particle velocity space is replaced by a finite dimensional discrete set. The number of linearly independent velocity moments in a lattice Boltzmann model cannot exceed the number of discrete velocities. Thus, finite dimensionality introduces linear dependencies among the moments that do not exist in the exact continuous theory. Given a discrete velocity set, it is important to know to exactly what order moments are free of these dependencies. Elementary group theory is applied to the solution of this problem. It is found that by decomposing the velocity set into subsets that transform among themselves under an appropriate symmetry group, it becomes relatively straightforward to assess the behavior of moments in the theory. The construction of some standard two- and three-dimensional models is reviewed from this viewpoint, and procedures for constructing some new higher dimensional models are suggested.
NASA Astrophysics Data System (ADS)
Held, M.; Kendl, A.
2015-10-01
A lattice Boltzmann method (LBM) approach to the Charney-Hasegawa-Mima (CHM) model for adiabatic drift wave turbulence in magnetised plasmas is implemented. The CHM-LBM model contains a barotropic equation of state for the potential, a force term including a cross-product analogous to the Coriolis force in quasigeostrophic models, and a density gradient source term. Expansion of the resulting lattice Boltzmann model equations leads to cold-ion fluid continuity and momentum equations, which resemble CHM dynamics under drift ordering. The resulting numerical solutions of standard test cases (monopole propagation, stable drift modes and decaying turbulence) are compared to results obtained by a conventional finite difference scheme that directly discretizes the CHM equation. The LB scheme resembles characteristic CHM dynamics apart from an additional shear in the density gradient direction. The occurring shear reduces with the drift ratio and is ascribed to the compressible limit of the underlying LBM.
NASA Astrophysics Data System (ADS)
Ding, E. J.
2015-06-01
The time-independent lattice Boltzmann algorithm (TILBA) is developed to calculate the hydrodynamic interactions between two particles in a Stokes flow. The TILBA is distinguished from the traditional lattice Boltzmann method in that a background matrix (BGM) is generated prior to the calculation. The BGM, once prepared, can be reused for calculations for different scenarios, and the computational cost for each such calculation will be significantly reduced. The advantage of the TILBA is that it is easy to code and can be applied to any particle shape without complicated implementation, and the computational cost is independent of the shape of the particle. The TILBA is validated and shown to be accurate by comparing calculation results obtained from the TILBA to analytical or numerical solutions for certain problems.
Momentum-exchange method in lattice Boltzmann simulations of particle-fluid interactions
NASA Astrophysics Data System (ADS)
Chen, Yu; Cai, Qingdong; Xia, Zhenhua; Wang, Moran; Chen, Shiyi
2013-07-01
The momentum exchange method has been widely used in lattice Boltzmann simulations for particle-fluid interactions. Although proved accurate for still walls, it will result in inaccurate particle dynamics without corrections. In this work, we reveal the physical cause of this problem and find that the initial momentum of the net mass transfer through boundaries in the moving-boundary treatment is not counted in the conventional momentum exchange method. A corrected momentum exchange method is then proposed by taking into account the initial momentum of the net mass transfer at each time step. The method is easy to implement with negligible extra computation cost. Direct numerical simulations of a single elliptical particle sedimentation are carried out to evaluate the accuracy for our method as well as other lattice Boltzmann-based methods by comparisons with the results of the finite element method. A shear flow test shows that our method is Galilean invariant.
NASA Astrophysics Data System (ADS)
Obliger, Amaël; Duvail, Magali; Jardat, Marie; Coelho, Daniel; Békri, Samir; Rotenberg, Benjamin
2013-07-01
We report the calculation of all the transfer coefficients which couple the solvent and ionic fluxes through a charged pore under the effect of pressure, electrostatic potential, and concentration gradients. We use a combination of analytical calculations at the Poisson-Nernst-Planck and Navier-Stokes levels of description and mesoscopic lattice simulations based on kinetic theory. In the absence of added salt, i.e., when the only ions present in the fluid are the counterions compensating the charge of the surface, exact analytical expressions for the fluxes in cylindrical pores allow us to validate a new lattice-Boltzmann electrokinetics (LBE) scheme which accounts for the osmotic contribution to the transport of all species. The influence of simulation parameters on the numerical accuracy is thoroughly investigated. In the presence of an added salt, we assess the range of validity of approximate expressions of the fluxes computed from the linearized Poisson-Boltzmann equation by a systematic comparison with LBE simulations.
NASA Astrophysics Data System (ADS)
Kang, XiuYing; Su, YanPing
2012-10-01
Cross-flows around two, three and four circular cylinders in tandem, side-by-side, isosceles triangle and square arrangements are simulated using the incompressible lattice Boltzmann method with a second-order accurate curved boundary condition at Reynolds number 200 and the cylinder center-to-center transverse or/and longitudinal spacing 1.5 D, where D is the identical circular cylinder diameter. The wake patterns, pressure and force distributions on the cylinders and mechanism of flow dynamics are investigated and compared among the four cases. The results also show that flows around the three or four cylinders significantly differ from those of the two cylinders in the tandem and side-by-side arrangements although there are some common features among the four cases due to their similarity of structures, which are interesting, complex and useful for practical applications. This study provides a useful database to validate the simplicity, accuracy and robustness of the Lattice Boltzmann method.
Flow force and torque on submerged bodies in lattice-Boltzmann methods via momentum exchange
NASA Astrophysics Data System (ADS)
Giovacchini, Juan P.; Ortiz, Omar E.
2015-12-01
We review the momentum exchange method to compute the flow force and torque on a submerged body in lattice-Boltzmann methods by presenting an alternative derivation. Our derivation does not depend on a particular implementation of the boundary conditions at the body surface, and it relies on general principles. After the introduction of the momentum exchange method in lattice-Boltzmann methods, some formulations were introduced to compute the fluid force on static and moving bodies. These formulations were introduced in a rather intuitive, ad hoc way. In our derivation, we recover the proposals most frequently used, in some cases with minor corrections, gaining some insight into the two most used formulations. At the end, we present some numerical tests to compare different approaches on a well-known benchmark test that support the correctness of the formulas derived.
Regularized lattice Boltzmann model for a class of convection-diffusion equations.
Wang, Lei; Shi, Baochang; Chai, Zhenhua
2015-10-01
In this paper, a regularized lattice Boltzmann model for a class of nonlinear convection-diffusion equations with variable coefficients is proposed. The main idea of the present model is to introduce a set of precollision distribution functions that are defined only in terms of macroscopic moments. The Chapman-Enskog analysis shows that the nonlinear convection-diffusion equations can be recovered correctly. Numerical tests, including Fokker-Planck equations, Buckley-Leverett equation with discontinuous initial function, nonlinear convection-diffusion equation with anisotropic diffusion, are carried out to validate the present model, and the results show that the present model is more accurate than some available lattice Boltzmann models. It is also demonstrated that the present model is more stable than the traditional single-relaxation-time model for the nonlinear convection-diffusion equations. PMID:26565368
Prediction of sound absorption in rigid porous media with the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
da Silva, Andrey Ricardo; Mareze, Paulo; Brandão, Eric
2016-02-01
In this work, sound absorption phenomena associated with the viscous shear stress within rigid porous media is investigated with a simple isothermal lattice Boltzmann BGK model. Simulations are conducted for different macroscopic material properties such as sample thickness and porosity and the results are compared with the exact analytical solution for materials with slit-like structure in terms of acoustic impedance and sound absorption coefficient. The numerical results agree very well with the exact solution, particularly for the sound absorption coefficient. The small deviations found in the low frequency limit for the real part of the acoustic impedance are attributed to the ratio between the thicknesses of the slit and the viscous boundary layer. The results suggest that the lattice Boltzmann method can be a very compelling numerical tool for simulating viscous sound absorption phenomena in the time domain, particularly due to its computational simplicity when compared to traditional continuum based techniques.
Rossby vortex simulation on a paraboloidal coordinate system using the lattice Boltzmann method.
Yu, H; Zhao, K
2001-11-01
In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc., can be expected. PMID:11736137
A study of multiphase flow in fractured porous media using a microscale lattice Boltzmann approach
Soll, W.E.; Eggert, K.E.; Grunau, D.W.; Schafer-Perini, A.L.
1994-02-01
The lattice Boltzmann technique has been shown to be an efficient and reliable approach to modeling single- and multi-fluid flow in porous media systems. The flexibility of this approach in discretizing the pore/solid space means it is particularly well suited to capturing fluid behavior, fluid-fluid interactions, and fluid-solid interactions at the scale of the individual pores. Such flexibility readily lends itself to studying processes occurring at physical interfaces, such as between a fracture and the surrounding porous matrix. Here we present pore-level simulations of fluid flow through a fracture embedded in an unsaturated matrix. Simulations are run on the massively parallel Connection Machine 5 (CM-5) using the two-fluid, two-dimensional lattice Boltzmann flow simulator developed at Los Alamos National Laboratory. We look at the effect of pressure gradients and initial matrix saturation on infiltration into the matrix and fluid flow along the fracture.
A Lattice-Boltzmann model for suspensions of self-propelling colloidal particles.
Ramachandran, S; Sunil Kumar, P B; Pagonabarraga, I
2006-06-01
We present a Lattice-Boltzmann method for simulating self-propelling (active) colloidal particles in two dimensions. Active particles with symmetric and asymmetric force distribution on their surface are considered. The velocity field generated by a single active particle, changing its orientation randomly, and the different time scales involved are characterized in detail. The steady-state speed distribution in the fluid, resulting from the activity, is shown to deviate considerably from the equilibrium distribution. PMID:16779527
A Lattice-Boltzmann model for suspensions of self-propelling colloidal particles
NASA Astrophysics Data System (ADS)
Ramachandran, S.; Kumar, P. B. Sunil; Pagonabarraga, I.
2006-06-01
We present a Lattice-Boltzmann method for simulating self-propelling (active) colloidal particles in two dimensions. Active particles with symmetric and asymmetric force distribution on their surface are considered. The velocity field generated by a single active particle, changing its orientation randomly, and the different time scales involved are characterized in detail. The steady-state speed distribution in the fluid, resulting from the activity, is shown to deviate considerably from the equilibrium distribution.
Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method
NASA Technical Reports Server (NTRS)
Chen, Hudong; Chen, Shiyi; Matthaeus, William H.
1992-01-01
A lattice Boltzmann model is presented which gives the complete Navier-Stokes equation and may provide an efficient parallel numerical method for solving various fluid problems. The model uses the single-time relaxation approximation and a particular Maxwell-type distribution. The model eliminates exactly (1) the non-Galilean invariance caused by a density-dependent coefficient in the convection term and (2) a velocity-dependent equation of state.
On the lattice Boltzmann method for multiphase flows with large density ratios
NASA Astrophysics Data System (ADS)
Kim, Seung Hyun; Pitsch, Heinz
2015-12-01
An analysis of the lattice Boltzmann (LB) method for multiphase flows with large density ratios is presented. It is shown that for incompressible, multiphase LB methods, the divergence-free condition is not satisfied within the formal accuracy of the LB method, when the density ratio between the two phases is large enough. The discrete differentiation-by-parts rule is responsible for this error. A new multiphase LB method to resolve this issue is proposed.
Entropic Lattice Boltzmann Methods for Fluid Mechanics: Thermal, Multi-phase and Turbulence
NASA Astrophysics Data System (ADS)
Chikatamarla, Shyam; Boesch, F.; Frapolli, N.; Mazloomi, A.; Karlin, I.
2014-11-01
With its roots in statistical mechanics and kinetic theory, the lattice Boltzmann method (LBM) is a paradigm-changing innovation, offering for the first time an intrinsically parallel CFD algorithm. Over the past two decades, LBM has achieved numerous results in the field of CFD and is now in a position to challenge state-of-the art CFD techniques. Major restyling of LBM resulted in an unconditionally stable entropic LBM which restored Second Law (Boltzmann H theorem) in the LBM kinetics and thus enabled affordable direct simulations of fluid turbulence. In this talk, we shall review recent advances in ELBM as a practical, modeling-free tool for simulation of complex flow phenomenon. We shall present recent simulations of fluid turbulence including turbulent channel flow, flow past a circular cylinder, creation and dynamics of vortex tubes, and flow past a surface mounted cube. Apart from its achievements in turbulent flow simulations, ELBM has also presented us the opportunity to extend lattice Boltzmann method to higher order lattices which shall be employed for turbulent, multi-phase and thermal flow simulations. A new class of entropy functions are proposed to handle non-ideal equation of state and surface tension terms in multi-phase flows. It is shown the entropy principle brings unconditional stability and thermodynamic consistency to all the three flow regimes considered here. Acknowledgements: ERC Advanced Grant ``ELBM'' and CSCS grant s437 are deeply acknowledged. References:
Modeling flue pipes: Subsonic flow, lattice Boltzmann, and parallel distributed computers
NASA Astrophysics Data System (ADS)
Skordos, Panayotis A.
1995-01-01
The problem of simulating the hydrodynamics and the acoustic waves inside wind musical instruments such as the recorder the organ, and the flute is considered. The problem is attacked by developing suitable local-interaction algorithms and a parallel simulation system on a cluster of non-dedicated workstations. Physical measurements of the acoustic signal of various flue pipes show good agreement with the simulations. Previous attempts at this problem have been frustrated because the modeling of acoustic waves requires small integration time steps which make the simulation very compute-intensive. In addition, the simulation of subsonic viscous compressible flow at high Reynolds numbers is susceptible to slow-growing numerical instabilities which are triggered by high-frequency acoustic modes. The numerical instabilities are mitigated by employing suitable explicit algorithms: lattice Boltzmann method, compressible finite differences, and fourth-order artificial-viscosity filter. Further, a technique for accurate initial and boundary conditions for the lattice Boltzmann method is developed, and the second-order accuracy of the lattice Boltzmann method is demonstrated. The compute-intensive requirements are handled by developing a parallel simulation system on a cluster of non-dedicated workstations. The system achieves 80 percent parallel efficiency (speedup/processors) using 20 HP-Apollo workstations. The system is built on UNIX and TCP/IP communication routines, and includes automatic process migration from busy hosts to free hosts.
Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications
NASA Astrophysics Data System (ADS)
Wu, J.; Shu, C.
2009-04-01
A version of immersed boundary-lattice Boltzmann method (IB-LBM) is proposed in this work. It is based on the lattice Boltzmann equation with external forcing term proposed by Guo et al. [Z. Guo, C. Zheng, B. Shi, Discrete lattice effects on the forcing term in the lattice Boltzmann method, Phys. Rev. E 65 (2002) 046308], which can well consider the effect of external force to the momentum and momentum flux as well as the discrete lattice effect. In this model, the velocity is contributed by two parts. One is from the density distribution function and can be termed as intermediate velocity, and the other is from the external force and can be considered as velocity correction. In the conventional IB-LBM, the force density (external force) is explicitly computed in advance. As a result, we cannot manipulate the velocity correction to enforce the non-slip boundary condition at the boundary point. In the present work, the velocity corrections (force density) at all boundary points are considered as unknowns which are computed in such a way that the non-slip boundary condition at the boundary points is enforced. The solution procedure of present IB-LBM is exactly the same as the conventional IB-LBM except that the non-slip boundary condition can be satisfied in the present model while it is only approximately satisfied in the conventional model. Numerical experiments for the flows around a circular cylinder and an airfoil show that there is no any penetration of streamlines to the solid body in the present results. This is not the case for the results obtained by the conventional IB-LBM. Another advantage of the present method is its simple calculation of force on the boundary. The force can be directly calculated from the relationship between the velocity correction and the force density.
Design of Chern and Mott insulators in buckled 3 d oxide honeycomb lattices
NASA Astrophysics Data System (ADS)
Doennig, David; Baidya, Santu; Pickett, Warren E.; Pentcheva, Rossitza
2016-04-01
Perovskite (La X O3 )2/(LaAlO3)4(111) superlattices with X spanning the entire 3 d transition-metal series combine the strongly correlated, multiorbital nature of electrons in transition-metal oxides with a honeycomb lattice as a key feature. Based on density functional theory calculations including strong interaction effects, we establish trends in the evolution of electronic states as a function of several control parameters: band filling, interaction strength, spin-orbit coupling (SOC), and lattice instabilities. Competition between local pseudocubic and global trigonal symmetry as well as the additional flexibility provided by the magnetic and spin degrees of freedom of 3 d ions lead to a broad array of distinctive broken-symmetry ground states not accessible for the (001)-growth direction, offering a platform to design two-dimensional electronic functionalities. Constraining the symmetry between the two triangular sublattices causes X =Mn , Co, and Ti to emerge as Chern insulators driven by SOC. For X =Mn we illustrate how interaction strength and lattice distortions can tune these systems between a Dirac semimetal, a Chern and a trivial Mott insulator.
NASA Astrophysics Data System (ADS)
Xia, Yong; Lu, De-Tang; Liu, Yang; Xu, You-Sheng
2009-03-01
The multiple-relaxation-time lattice Boltzmann method (MRT-LBM) is implemented to numerically simulate the cross flow over a longitudinal vibrating circular cylinder. This research is carried out on a three-dimensional (3D) finite cantilevered cylinder to investigate the effect of forced vibration on the wake characteristics and the 3D effect of a cantilevered cylinder. To meet the accuracy of this method, the present calculation is carried out at a low Reynolds number Re = 100, as well as to make the vibration obvious, we make the vibration strong enough. The calculation results indicate that the vibration has significant influence on the wake characteristics. When the vibrating is big enough, our early works show that the 2D vortex shedding would be locked up by vibration. Contrarily, this phenomenon would not appear in the present 3D case because of the end effect of the cantilevered cylinder.
NASA Astrophysics Data System (ADS)
Zhang, Jianying; Yan, Guangwu; Wang, Moran
2016-02-01
A lattice Boltzmann model for solving the three-dimensional cubic-quintic complex Ginzburg-Landau equation (CQCGLE) is proposed. Differently from the classic lattice Boltzmann models, this lattice Boltzmann model is based on uniformly distributed lattice points in a three-dimensional space, and the evolution of the model is about a spatial axis rather than time. The algorithm provides advantages similar to the lattice Boltzmann method in that it is easily adapted to complex Ginzburg-Landau equations. Examples show that the model accurately reproduces the vortex tori pattern in the CQCGLE.
Linearized lattice Boltzmann method for micro- and nanoscale flow and heat transfer
NASA Astrophysics Data System (ADS)
Shi, Yong; Yap, Ying Wan; Sader, John E.
2015-07-01
Ability to characterize the heat transfer in flowing gases is important for a wide range of applications involving micro- and nanoscale devices. Gas flows away from the continuum limit can be captured using the Boltzmann equation, whose analytical solution poses a formidable challenge. An efficient and accurate numerical simulation of the Boltzmann equation is thus highly desirable. In this article, the linearized Boltzmann Bhatnagar-Gross-Krook equation is used to develop a hierarchy of thermal lattice Boltzmann (LB) models based on half-space Gaussian-Hermite (GH) quadrature ranging from low to high algebraic precision, using double distribution functions. Simplified versions of the LB models in the continuum limit are also derived, and are shown to be consistent with existing thermal LB models for noncontinuum heat transfer reported in the literature. Accuracy of the proposed LB hierarchy is assessed by simulating thermal Couette flows for a wide range of Knudsen numbers. Effects of the underlying quadrature schemes (half-space GH vs full-space GH) and continuum-limit simplifications on computational accuracy are also elaborated. The numerical findings in this article provide direct evidence of improved computational capability of the proposed LB models for modeling noncontinuum flows and heat transfer at small length scales.
Linearized lattice Boltzmann method for micro- and nanoscale flow and heat transfer.
Shi, Yong; Yap, Ying Wan; Sader, John E
2015-07-01
Ability to characterize the heat transfer in flowing gases is important for a wide range of applications involving micro- and nanoscale devices. Gas flows away from the continuum limit can be captured using the Boltzmann equation, whose analytical solution poses a formidable challenge. An efficient and accurate numerical simulation of the Boltzmann equation is thus highly desirable. In this article, the linearized Boltzmann Bhatnagar-Gross-Krook equation is used to develop a hierarchy of thermal lattice Boltzmann (LB) models based on half-space Gaussian-Hermite (GH) quadrature ranging from low to high algebraic precision, using double distribution functions. Simplified versions of the LB models in the continuum limit are also derived, and are shown to be consistent with existing thermal LB models for noncontinuum heat transfer reported in the literature. Accuracy of the proposed LB hierarchy is assessed by simulating thermal Couette flows for a wide range of Knudsen numbers. Effects of the underlying quadrature schemes (half-space GH vs full-space GH) and continuum-limit simplifications on computational accuracy are also elaborated. The numerical findings in this article provide direct evidence of improved computational capability of the proposed LB models for modeling noncontinuum flows and heat transfer at small length scales. PMID:26274307
Zu, Y Q; He, S
2013-04-01
A lattice Boltzmann model (LBM) is proposed based on the phase-field theory to simulate incompressible binary fluids with density and viscosity contrasts. Unlike many existing diffuse interface models which are limited to density matched binary fluids, the proposed model is capable of dealing with binary fluids with moderate density ratios. A new strategy for projecting the phase field to the viscosity field is proposed on the basis of the continuity of viscosity flux. The new LBM utilizes two lattice Boltzmann equations (LBEs): one for the interface tracking and the other for solving the hydrodynamic properties. The LBE for interface tracking can recover the Chan-Hilliard equation without any additional terms; while the LBE for hydrodynamic properties can recover the exact form of the divergence-free incompressible Navier-Stokes equations avoiding spurious interfacial forces. A series of 2D and 3D benchmark tests have been conducted for validation, which include a rigid-body rotation, stationary and moving droplets, a spinodal decomposition, a buoyancy-driven bubbly flow, a layered Poiseuille flow, and the Rayleigh-Taylor instability. It is shown that the proposed method can track the interface with high accuracy and stability and can significantly and systematically reduce the parasitic current across the interface. Comparisons with momentum-based models indicate that the newly proposed velocity-based model can better satisfy the incompressible condition in the flow fields, and eliminate or reduce the velocity fluctuations in the higher-pressure-gradient region and, therefore, achieve a better numerical stability. In addition, the test of a layered Poiseuille flow demonstrates that the proposed scheme for mixture viscosity performs significantly better than the traditional mixture viscosity methods. PMID:23679542
NASA Astrophysics Data System (ADS)
Zu, Y. Q.; He, S.
2013-04-01
A lattice Boltzmann model (LBM) is proposed based on the phase-field theory to simulate incompressible binary fluids with density and viscosity contrasts. Unlike many existing diffuse interface models which are limited to density matched binary fluids, the proposed model is capable of dealing with binary fluids with moderate density ratios. A new strategy for projecting the phase field to the viscosity field is proposed on the basis of the continuity of viscosity flux. The new LBM utilizes two lattice Boltzmann equations (LBEs): one for the interface tracking and the other for solving the hydrodynamic properties. The LBE for interface tracking can recover the Chan-Hilliard equation without any additional terms; while the LBE for hydrodynamic properties can recover the exact form of the divergence-free incompressible Navier-Stokes equations avoiding spurious interfacial forces. A series of 2D and 3D benchmark tests have been conducted for validation, which include a rigid-body rotation, stationary and moving droplets, a spinodal decomposition, a buoyancy-driven bubbly flow, a layered Poiseuille flow, and the Rayleigh-Taylor instability. It is shown that the proposed method can track the interface with high accuracy and stability and can significantly and systematically reduce the parasitic current across the interface. Comparisons with momentum-based models indicate that the newly proposed velocity-based model can better satisfy the incompressible condition in the flow fields, and eliminate or reduce the velocity fluctuations in the higher-pressure-gradient region and, therefore, achieve a better numerical stability. In addition, the test of a layered Poiseuille flow demonstrates that the proposed scheme for mixture viscosity performs significantly better than the traditional mixture viscosity methods.
The sign-factor of the 3D Ising model on dual BCC lattice
NASA Astrophysics Data System (ADS)
Khachatryan, Sh.; Sedrakyan, A.
2002-01-01
We modify the two-dimensional model for the sign-factor of the regular 3D Ising model (3DIM) presented by Kavalov and Sedrakyan (Phys. Lett. 173B (1986) 449 and Nucl. Phys. 285B (1987) 264) for the case of dual to body centered cubic (DBCC) three-dimensional lattice. The advantage of this lattice is in an absence of self-intersections of the two-dimensional surfaces embedded there. We investigate simpler case of the model with scalar fermions (instead of SU(2) needed for 3DIM) and have found it's spectrum, which appeared to be massless. We reformulate the model by use of R-matrix formalism and a new interesting structure appears in a necessity to introduce three-particle R(3)ijk-matrices. We formulate the integrability property of the model for more general case.
Deconfinement Phase Transition in a 3D Nonlocal U(1) Lattice Gauge Theory
Arakawa, Gaku; Ichinose, Ikuo; Matsui, Tetsuo; Sakakibara, Kazuhiko
2005-06-03
We introduce a 3D compact U(1) lattice gauge theory having nonlocal interactions in the temporal direction, and study its phase structure. The model is relevant for the compact QED{sub 3} and strongly correlated electron systems like the t-J model of cuprates. For a power-law decaying long-range interaction, which simulates the effect of gapless matter fields, a second-order phase transition takes place separating the confinement and deconfinement phases. For an exponentially decaying interaction simulating matter fields with gaps, the system exhibits no signals of a second-order transition.
Lattice Boltzmann simulation of electrostatic double layer interaction force for nanoparticles
NASA Astrophysics Data System (ADS)
Shi, Grace X.; Jin, Yan; Lazouskaya, Volha; Wang, Chao; Wang, Lian-Ping
2011-11-01
Modeling the transport and retention of nanoparticles (NPs) through soil porous media requires an accurate description of the electrostatic interaction force between a nanoparticle and soil grain. In this study, we apply the lattice Boltzmann method to directly solve the nonlinear Poisson Boltzmann (PB) equation for several geometric configurations including plate-plate, NP-plate, and NP-NP interactions, for any surface potentials and interaction distances and for different boundary conditions. Interaction energy and force are then derived from the simulations. For the case of plate-plate interaction, the simulation results are compared to the exact solution of the nonlinear PB equation. It is shown that the linear PB solution is valid when the nondimensional surface potential is less than one, and that the linear PB solution over-predicts the interaction force for intermediate gap distances but under-predicts the force for small gap distances. For NP-plate and NP-NP interactions, an axisymmetric lattice Boltzmann formulation is developed to solve the governing equations. The results will be compared to the classic approximate expressions of interaction force to evaluate their validity and to study the effect of nanoparticle size. Work supported by NSF and USDA.
Theory of the lattice Boltzmann method: acoustic and thermal properties in two and three dimensions.
Lallemand, Pierre; Luo, Li-Shi
2003-09-01
The focus of the present work is to provide an analysis for the acoustic and thermal properties of the energy-conserving lattice Boltzmann models, and a solution to the numerical defects and instability associated with these models in two and three dimensions. We discover that a spurious algebraic coupling between the shear and energy modes of the linearized evolution operator is a defect universal to the energy-conserving Boltzmann models in two and three dimensions. This spurious mode coupling is highly anisotropic and may occur at small values of wave number k along certain directions, and it is a direct consequence of the following key features of the lattice Boltzmann equation: (1) its simple spatial-temporal dynamics, (2) the linearity of the relaxation modeling for collision operator, and (3) the energy-conservation constraint. To eliminate the spurious mode coupling, we propose a hybrid thermal lattice Boltzmann equation (HTLBE) in which the mass and momentum conservation equations are solved by using the multiple-relaxation-time model due to d'Humières, whereas the diffusion-advection equation for the temperature is solved separately by using finite-difference technique (or other means). Through the Chapman-Enskog analysis we show that the hydrodynamic equations derived from the proposed HTLBE model include the equivalent effect of gamma=C(P)/C(V) in both the speed and attenuation of sound. Appropriate coupling between the energy and velocity field is introduced to attain correct acoustics in the model. The numerical stability of the HTLBE scheme is analyzed by solving the dispersion equation of the linearized collision operator. We find that the numerical stability of the lattice Boltzmann scheme improves drastically once the spurious mode coupling is removed. It is shown that the HTLBE scheme is far superior to the existing thermal LBE schemes in terms of numerical stability, flexibility, and possible generalization for complex fluids. We also present
Discovery of a 3d-transition-metal-based ferromagnetic Kondo lattice system
NASA Astrophysics Data System (ADS)
Us Saleheen, Ahmad; Samanta, Tapas; Lepkowski, Daniel; Shankar, Alok; Prestigiacomo, Joseph; Dubenko, Igor; Quetz, Abdiel; McDougald, Roy, Jr.; McCandless, Gregory; Chan, Julia; Adams, Philip; Young, David; Ali, Naushad; Stadler, Shane
2015-03-01
The formation of a Kondo lattice results in a wide variety of exotic phenomena associated with the competition between the Kondo effect and the RKKY interaction, such as heavy fermions, non-Fermi liquid behavior, unconventional superconductivity, and so on. A quantum critical point (QCP) has been frequently observed at the boundaries of competing phases for antiferromagnetic materials. However, the existence of a ferromagnetic (FM) QCP is unclear. Moreover, FM Kondo lattices are rare. Here we report the discovery of a FM Kondo lattice system Mn1-xFexCoGe, which is the first example of a 3d-metal-based system (i.e., not rare-earth-based). Resistivity, magnetic susceptibility, heat capacity and thermopower studies on a single crystal sample indicate that the anisotropic FM kondo lattice has formed along c-axis. The signature of a spin density wave transition was also observed above the Kondo minimum, below which the resistivity follows a log(T) behavior. This work was supported by the U.S. Department of Energy (Grant Nos. DE-FG02-13ER46946 and DE-FG02-06ER46291).
NASA Astrophysics Data System (ADS)
Fiorentino, Eve-Agnès; Toussaint, Renaud; Jouniaux, Laurence
2014-05-01
We study the coupling between hydraulic and electric flows in a porous medium at small scale using the Lattice Boltzmann method. This method is a computational fluid dynamics technique that is used for advection and diffusion modeling. We implement a coupled Lattice Boltzmann algorithm that solves both the mass transport and the electric field arising from charges displacements. The streaming potential and electroosmosis phenomena occur in a variety of situations and derive from this coupling. We focus on the streaming potential which is described using the ratio between the created potential difference and the applied pressure gradient. The streaming potential is assumed to be a linear function of the fluid conductivity, but experimental results highlight anomalous behaviors at low and high salinity. We try to account for them by setting extreme conditions that are likely to generate non-linearities. Several pore radii are tested so as to determine what is the effect of a radius that is comparable to the Debye length, the screening length of the electric potential, due to the ions in the electrolyte. The volumetric integral of the electrical current is calculated for comparison with the 2D simulations. High values of zeta potential are tested to verify if the discrepancy regarding the theoretical result is concentration-dependent. We try to include a surface conductivity term in the coefficient formulation. Some tests including a rugosity on the channel walls are performed. All of these attempts show a normal behaviour of the streaming potential at high salinity. We observe a decrease of the ratio at low conductivity, showing that this ratio is modified when the pore radius becomes negligible compared with the Debye length, which is physically meaningful in little pores at low concentrations. References : S. Pride. Governing equations for the coupled electromagnetics and acoustics of porous media. Physical Review B, 50 : 15678-15696, 1994. D. A. Wolf
Boundary conditions of the lattice Boltzmann method for convection-diffusion equations
NASA Astrophysics Data System (ADS)
Huang, Juntao; Yong, Wen-An
2015-11-01
In this paper, we employ an asymptotic analysis technique and construct two boundary schemes accompanying the lattice Boltzmann method for convection-diffusion equations with general Robin boundary conditions. One scheme is for straight boundaries, with the boundary points locating at any distance from the lattice nodes, and has second-order accuracy. The other is for curved boundaries, has only first-order accuracy and is much simpler than the existing schemes. Unlike those in the literature, our schemes involve only the current lattice node. Such a "single-node" boundary schemes are highly desirable for problems with complex geometries. The two schemes are validated numerically with a number of examples. The numerical results show the utility of the constructed schemes and very well support our theoretical predications.
Lattice Boltzmann technique for heat transport phenomena coupled with melting process
NASA Astrophysics Data System (ADS)
Ibrahem, A. M.; El-Amin, M. F.; Mohammadein, A. A.; Gorla, Rama Subba Reddy
2016-04-01
In this work, the heat transport phenomena coupled with melting process are studied by using the enthalpy-based lattice Boltzmann method (LBM). The proposed model is a modified version of thermal LB model, where could avoid iteration steps and ensures high accuracy. The Bhatnagar-Gross-Krook (BGK) approximation with a D1Q2 lattice was used to determine the temperature field for one-dimensional melting by conduction and multi-distribution functions (MDF) with D2Q9 lattice was used to determine the density, velocity and temperature fields for two-dimensional melting by natural convection. Different boundary conditions including Dirichlet, adiabatic and bounce-back boundary conditions were used. The influence of increasing Rayleigh number (from 103 to 105) on temperature distribution and melting process is studied. The obtained results show that a good agreement with the analytical solution for melting by conduction case and with the benchmark solution for melting by convection.
Lattice-Boltzmann method for the simulation of transport phenomena in charged colloids.
Horbach, J; Frenkel, D
2001-12-01
We present a simulation scheme based on the lattice-Boltzmann method to simulate the dynamics of charged colloids in an electrolyte. In our model we describe the electrostatics on the level of a Poisson-Boltzmann equation and the hydrodynamics of the fluid by the linearized Navier-Stokes equations. We verify our simulation scheme by means of a Chapman-Enskog expansion. Our method is applied to the calculation of the reduced sedimentation velocity U/U(0) for a cubic array of charged spheres in an electrolyte. We show that we recover the analytical solution first derived by Booth [F. Booth, J. Chem. Phys. 22, 1956 (1954)] for a weakly charged, isolated sphere in an unbounded electrolyte. The present method makes it possible to go beyond the Booth theory, and we discuss the dependence of the sedimentation velocity on the charge of the spheres. Finally we compare our results to experimental data. PMID:11736191
Radiative or neutron transport modeling using a lattice Boltzmann equation framework
NASA Astrophysics Data System (ADS)
Bindra, H.; Patil, D. V.
2012-07-01
In this paper, the lattice Boltzmann equation (LBE)-based framework is used to obtain the solution for the linear radiative or neutron transport equation. The LBE framework is devised for the integrodifferential forms of these equations which arise due to the inclusion of the scattering terms. The interparticle collisions are neglected, hence omitting the nonlinear collision term. Furthermore, typical representative examples for one-dimensional or two-dimensional geometries and inclusion or exclusion of the scattering term (isotropic and anisotropic) in the Boltzmann transport equation are illustrated to prove the validity of the method. It has been shown that the solution from the LBE methodology is equivalent to the well-known Pn and Sn methods. This suggests that the LBE can potentially provide a more convenient and easy approach to solve the physical problems of neutron and radiation transport.