Engineering Computational Mechanics

Lattice Boltzmann Simulation Studies

Description

This cluster of papers focuses on the application and development of the Lattice Boltzmann Method (LBM) for simulating complex flows, including multiphase flow, fluid-structure interaction, and heat transfer in porous media. It also explores the integration of LBM with immersed boundary methods and its use in modeling binary fluid systems and kinetic theory representation of hydrodynamics.

Keywords

Lattice Boltzmann Method; Immersed Boundary Method; Multiphase Flow; Complex Geometries; Fluid-Structure Interaction; Porous Media; Kinetic Theory; Thermal Properties; Binary Fluid Systems; Immersion Techniques

Since the publication of the first edition of Multiphase Flow with Droplets and Particles, there have been significant advances in science and engineering applications of multiphase fluid flow. Maintaining the … Since the publication of the first edition of Multiphase Flow with Droplets and Particles, there have been significant advances in science and engineering applications of multiphase fluid flow. Maintaining the pedagogical approach that made the first edition so popular, this second edition provides a background in this important area of fluid mecha
A new technique is described for the numerical investigation of the time-dependent flow of an incompressible fluid, the boundary of which is partially confined and partially free. The full Navier-Stokes … A new technique is described for the numerical investigation of the time-dependent flow of an incompressible fluid, the boundary of which is partially confined and partially free. The full Navier-Stokes equations are written in finite-difference form, and the solution is accomplished by finite-time-step advancement. The primary dependent variables are the pressure and the velocity components. Also used is a set of marker particles which move with the fluid. The technique is called the marker and cell method. Some examples of the application of this method are presented. All non-linear effects are completely included, and the transient aspects can be computed for as much elapsed time as desired.
We propose the lattice BGK models, as an alternative to lattice gases or the lattice Boltzmann equation, to obtain an efficient numerical scheme for the simulation of fluid dynamics. With … We propose the lattice BGK models, as an alternative to lattice gases or the lattice Boltzmann equation, to obtain an efficient numerical scheme for the simulation of fluid dynamics. With a properly chosen equilibrium distribution, the Navier-Stokes equation is obtained from the kinetic BGK equation at the second-order of approximation. Compared to lattice gases, the present model is noise-free, has Galileian invariance and a velocity-independent pressure. It involves a relaxation parameter that influences the stability of the new scheme. Numerical simulations are shown to confirm the speed of sound and the shear viscosity.
Abstract : A continuum mechanics approach to two-phase flow is reviewed. An averaging procedure is discussed and applied to the exact equations of motion. Constitutive equations are supplied and discussed … Abstract : A continuum mechanics approach to two-phase flow is reviewed. An averaging procedure is discussed and applied to the exact equations of motion. Constitutive equations are supplied and discussed for the stresses, pressure differences and the interfacial force. The nature of the resulting equations is studied. (Author)
A lattice Boltzmann scheme able to model the hydrodynamics of phase separation and two-phase flow is described. Thermodynamic consistency is ensured by introducing a nonideal pressure tensor directly into the … A lattice Boltzmann scheme able to model the hydrodynamics of phase separation and two-phase flow is described. Thermodynamic consistency is ensured by introducing a nonideal pressure tensor directly into the collision operator. We also show how an external chemical potential can be used to supplement standard boundary conditions in order to investigate the effect of wetting on phase separation and fluid flow in confined geometries. The approach has the additional advantage of reducing many of the unphysical discetrization problems common to previous lattice Boltzmann methods.
It is known that the Frisch-Hasslacher-Pomeau lattice-gas automaton model and related models possess some rather unphysical effects. These are (1) a non-Galilean invariance caused by a density-dependent coefficient in the … It is known that the Frisch-Hasslacher-Pomeau lattice-gas automaton model and related models possess some rather unphysical effects. These are (1) a non-Galilean invariance caused by a density-dependent coefficient in the convection term, and (2) a velocity-dependent equation of state. In this paper, we show that both of these effects can be eliminated exactly in a lattice Boltzmann-equation model.
We discuss an alternative technique to the lattice-gas automata for the study of hydrodynamic properties, namely, we propose to model the lattice gas with a Boltzmann equation. This approach completely … We discuss an alternative technique to the lattice-gas automata for the study of hydrodynamic properties, namely, we propose to model the lattice gas with a Boltzmann equation. This approach completely eliminates the statistical noise that plagues the usual lattice-gas simulations and therefore permits simulations that demand much less computer time. It is estimated to be more efficient than the lattice-gas automata for intermediate to low Reynolds number $R\ensuremath{\lesssim}100$.
We present the details of a lattice Boltzmann approach to phase separation in nonideal one- and two-component fluids. The collision rules are chosen such that the equilibrium state corresponds to … We present the details of a lattice Boltzmann approach to phase separation in nonideal one- and two-component fluids. The collision rules are chosen such that the equilibrium state corresponds to an input free energy and the bulk flow is governed by the continuity, Navier-Stokes, and, for the binary fluid, a convection-diffusion equation. Numerical results are compared to simple analytic predictions to confirm that the equilibrium state is indeed thermodynamically consistent and that the kinetics of the approach to equilibrium lie within the expected universality classes. The approach is compared to other lattice Boltzmann simulations of nonideal systems. \textcopyright{} 1996 The American Physical Society.
A lattice Boltzmann model is developed which has the ability to simulate flows containing multiple phases and components. Each of the components can be immiscible with the others and can … A lattice Boltzmann model is developed which has the ability to simulate flows containing multiple phases and components. Each of the components can be immiscible with the others and can have different mass values. The equilibrium state of each component can have a nonideal gas equation of state at a prescribed temperature exhibiting thermodynamic phase transitions. The scheme incorporated in this model is the introduction of an interparticle potential. The dynamical rules in this model are local so it is highly efficient to compute on massively parallel computers. This model has many applications in large-scale numerical simulations of various types of fluid flows.
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 … 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.
We introduce a lattice Boltzmann model for simulating immiscible binary fluids in two dimensions. The model, based on the Boltzmann equation of lattice-gas hydrodynamics, incorporates features of a previously introduced … We introduce a lattice Boltzmann model for simulating immiscible binary fluids in two dimensions. The model, based on the Boltzmann equation of lattice-gas hydrodynamics, incorporates features of a previously introduced discrete immiscible lattice-gas model. A theoretical value of the surface-tension coefficient is derived and found to be in excellent agreement with values obtained from simulations. The model serves as a numerical method for the simulation of immiscible two-phase flow; a preliminary application illustrates a simulation of flow in a two-dimensional microscopic model of a porous medium. Extension of the model to three dimensions appears straightforward.
A finite-difference method for solving the time-dependent NavierStokes equations for an incompressible fluid is introduced. This method uses the primitive variables, i.e. the velocities and the pressure, and is equally … A finite-difference method for solving the time-dependent NavierStokes equations for an incompressible fluid is introduced. This method uses the primitive variables, i.e. the velocities and the pressure, and is equally applicable to problems in two and three space dimensions. Test problems are solved, and an application to a three-dimensional convection problem is presented.
We show that discrete lattice effects must be considered in the introduction of a force into the lattice Boltzmann equation. A representation of the forcing term is then proposed. With … We show that discrete lattice effects must be considered in the introduction of a force into the lattice Boltzmann equation. A representation of the forcing term is then proposed. With the representation, the Navier-Stokes equation is derived from the lattice Boltzmann equation through the Chapman-Enskog expansion. Several other existing force treatments are also examined.
The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in … The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyperviscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters that optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some two-dimensional LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long-wavelength limit (wave vector k=0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.
This paper is concerned with the mathematical structure of the immersed boundary (IB) method, which is intended for the computer simulation of fluid–structure interaction, especially in biological fluid dynamics. The … This paper is concerned with the mathematical structure of the immersed boundary (IB) method, which is intended for the computer simulation of fluid–structure interaction, especially in biological fluid dynamics. The IB formulation of such problems, derived here from the principle of least action, involves both Eulerian and Lagrangian variables, linked by the Dirac delta function. Spatial discretization of the IB equations is based on a fixed Cartesian mesh for the Eulerian variables, and a moving curvilinear mesh for the Lagrangian variables. The two types of variables are linked by interaction equations that involve a smoothed approximation to the Dirac delta function. Eulerian/Lagrangian identities govern the transfer of data from one mesh to the other. Temporal discretization is by a second-order Runge–Kutta method. Current and future research directions are pointed out, and applications of the IB method are briefly discussed.
We describe in detail a recently proposed lattice-Boltzmann model [X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993)] for simulating flows with multiple phases and components. In particular, … We describe in detail a recently proposed lattice-Boltzmann model [X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993)] for simulating flows with multiple phases and components. In particular, the focus is on the modeling of one-component fluid systems which obey nonideal gas equations of state and can undergo a liquid-gas-type phase transition. The model is shown to be momentum conserving. From the microscopic mechanical stability condition, the densities in bulk liquid and gas phases are obtained as functions of a temperaturelike parameter. Comparisons with the thermodynamic theory of phase transitions show that the lattice-Boltzmann-equation model can be made to correspond exactly to an isothermal process. The density profile in the liquid-gas interface is also obtained as a function of the temperaturelike parameter and is shown to be isotropic. The surface tension, which can be changed independently, is calculated. The analytical conclusions are verified by numerical simulations.
This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of 15-velocity and 19-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in … This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of 15-velocity and 19-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.
We study the velocity boundary condition for curved boundaries in the lattice Boltzmann equation (LBE). We propose a LBE boundary condition for moving boundaries by combination of the “bounce-back” scheme … We study the velocity boundary condition for curved boundaries in the lattice Boltzmann equation (LBE). We propose a LBE boundary condition for moving boundaries by combination of the “bounce-back” scheme and spatial interpolations of first or second order. The proposed boundary condition is a simple, robust, efficient, and accurate scheme. Second-order accuracy of the boundary condition is demonstrated for two cases: (1) time-dependent two-dimensional circular Couette flow and (2) two-dimensional steady flow past a periodic array of circular cylinders (flow through the porous media of cylinders). For the former case, the lattice Boltzmann solution is compared with the analytic solution of the Navier–Stokes equation. For the latter case, the lattice Boltzmann solution is compared with a finite-element solution of the Navier–Stokes equation. The lattice Boltzmann solutions for both flows agree very well with the solutions of the Navier–Stokes equations. We also analyze the torque due to the momentum transfer between the fluid and the boundary for two initial conditions: (a) impulsively started cylinder and the fluid at rest, and (b) uniformly rotating fluid and the cylinder at rest.
▪ Abstract We present an overview of the lattice Boltzmann method (LBM), a parallel and efficient algorithm for simulating single-phase and multiphase fluid flows and for incorporating additional physical complexities. … ▪ Abstract We present an overview of the lattice Boltzmann method (LBM), a parallel and efficient algorithm for simulating single-phase and multiphase fluid flows and for incorporating additional physical complexities. The LBM is especially useful for modeling complicated boundary conditions and multiphase interfaces. Recent extensions of this method are described, including simulations of fluid turbulence, suspension flows, and reaction diffusion systems.
A new and very general technique for simulating solid–fluid suspensions is described; its most important feature is that the computational cost scales linearly with the number of particles. The method … A new and very general technique for simulating solid–fluid suspensions is described; its most important feature is that the computational cost scales linearly with the number of particles. The method combines Newtonian dynamics of the solid particles with a discretized Boltzmann equation for the fluid phase; the many-body hydrodynamic interactions are fully accounted for, both in the creeping-flow regime and at higher Reynolds numbers. Brownian motion of the solid particles arises spontaneously from stochastic fluctuations in the fluid stress tensor, rather than from random forces or displacements applied directly to the particles. In this paper, the theoretical foundations of the technique are laid out, illustrated by simple analytical and numerical examples; in a companion paper (Part 2), extensive numerical tests of the method, for stationary flows, time-dependent flows, and finite-Reynolds-number flows, are reported.
With its roots in kinetic theory and the cellular automaton concept, the lattice-Boltzmann (LB) equation can be used to obtain continuum flow quantities from simple and local update rules based … With its roots in kinetic theory and the cellular automaton concept, the lattice-Boltzmann (LB) equation can be used to obtain continuum flow quantities from simple and local update rules based on particle interactions. The simplicity of formulation and its versatility explain the rapid expansion of the LB method to applications in complex and multiscale flows. We review many significant developments over the past decade with specific examples. Some of the most active developments include the entropic LB method and the application of the LB method to turbulent flow, multiphase flow, and deformable particle and fiber suspensions. Hybrid methods based on the combination of the Eulerian lattice with a Lagrangian grid system for the simulation of moving deformable boundaries show promise for more efficient applications to a broader class of problems. We also discuss higher-order boundary conditions and the simulation of microchannel flow with finite Knudsen number. Additionally, the remarkable scalability of the LB method for parallel processing is shown with examples. Teraflop simulations with the LB method are routine, and there is no doubt that this method will be one of the first candidates for petaflop computational fluid dynamics in the near future.
The authors develop finite difference methods for elliptic equations of the form \[ \nabla \cdot (\beta (x)\nabla u(x)) + \kappa (x)u(x) = f(x)\] in a region $\Omega $ in one … The authors develop finite difference methods for elliptic equations of the form \[ \nabla \cdot (\beta (x)\nabla u(x)) + \kappa (x)u(x) = f(x)\] in a region $\Omega $ in one or two space dimensions. It is assumed that $\Omega $ is a simple region (e.g., a rectangle) and that a uniform rectangular grid is used. The situation is studied in which there is an irregular surface $\Gamma $ of codimension 1 contained in $\Omega $ across which $\beta ,\kappa $, and f may be discontinuous, and along which the source f may have a delta function singularity. As a result, derivatives of the solution u may be discontinuous across $\Gamma $. The specification of a jump discontinuity in u itself across $\Gamma $ is allowed. It is shown that it is possible to modify the standard centered difference approximation to maintain second order accuracy on the uniform grid even when $\Gamma $ is not aligned with the grid. This approach is also compared with a discrete delta function approach to handling singular sources, as used in Peskin’s immersed boundary method.
The field of fluid mechanics is rapidly advancing, driven by unprecedented volumes of data from experiments, field measurements, and large-scale simulations at multiple spatiotemporal scales. Machine learning (ML) offers a … The field of fluid mechanics is rapidly advancing, driven by unprecedented volumes of data from experiments, field measurements, and large-scale simulations at multiple spatiotemporal scales. Machine learning (ML) offers a wealth of techniques to extract ...Read More
A new and very general technique for simulating solid–fluid suspensions has been described in a previous paper (Part 1); the most important feature of the new method is that the … A new and very general technique for simulating solid–fluid suspensions has been described in a previous paper (Part 1); the most important feature of the new method is that the computational cost scales linearly with the number of particles. In this paper (Part 2), extensive numerical tests of the method are described; results are presented for creeping flows, both with and without Brownian motion, and at finite Reynolds numbers. Hydrodynamic interactions, transport coefficients, and the short-time dynamics of random dispersions of up to 1024 colloidal particles have been simulated.
Pressure (density) and velocity boundary conditions are studied for 2-D and 3-D lattice Boltzmann BGK models (LBGK) and a new method to specify these conditions is proposed. These conditions are … Pressure (density) and velocity boundary conditions are studied for 2-D and 3-D lattice Boltzmann BGK models (LBGK) and a new method to specify these conditions is proposed. These conditions are constructed in consistency with the wall boundary condition, based on the idea of bounceback of the non-equilibrium distribution. When these conditions are used together with the incompressible LBGK model [J. Stat. Phys. 81, 35 (1995)] the simulation results recover the analytical solution of the plane Poiseuille flow driven by a pressure (density) difference. The half-way wall bounceback boundary condition is also used with the pressure (density) inlet/outlet conditions proposed in this paper and in Phys. Fluids 8, 2527 (1996) to study 2-D Poiseuille flow and 3-D square duct flow. The numerical results are approximately second-order accurate. The magnitude of the error of the half-way wall bounceback boundary condition is comparable with that of other published boundary conditions and it has better stability behavior.
Abstract In recent years, certain forms of the Boltzmann equation--now going by the name of "Lattice Boltzmann equation" (LBE)--have emerged which relinquish most mathematical complexities of the true Boltzmann equation … Abstract In recent years, certain forms of the Boltzmann equation--now going by the name of "Lattice Boltzmann equation" (LBE)--have emerged which relinquish most mathematical complexities of the true Boltzmann equation without sacrificing physical fidelity in the description of complex fluid motion. This book provides the first detailed survey of LBE theory and its major applications to date. Accessible to a broad audience of scientists dealing with complex system dynamics, the book also portrays future developments in allied areas of science where fluid motion plays a distinguished role.
Abstract The reduction in permeability of sediments due to blockages caused by suspended fine particles is a common concern for the extraction processes of oil, natural gas, or methane gas … Abstract The reduction in permeability of sediments due to blockages caused by suspended fine particles is a common concern for the extraction processes of oil, natural gas, or methane gas from methane hydrate. In this study, the permeability reduction caused by suspended fine particles was newly modelled. Solid–water two-phase flow in frame sand sediment was numerically simulated by a three-dimensional Lattice Boltzmann method. For frame sand, shapes of real sand grains were extracted by series expansion of spherical harmonics from CT-scan images and packed in a microscopic computational domain. For each fine particle, a motion equation is solved using the pressure integrated on its surface with considering its collision to the frame sand surfaces. The calculated relative permeability could not be modelled solely by the volume saturation of the fine particles, but also their specific surface area was required.
Though ubiquitous in many engineering applications, including drug delivery, the compound droplet hydrodynamics in confined geometries have been barely surveyed. For the first time, this study thoroughly investigates the hydrodynamics … Though ubiquitous in many engineering applications, including drug delivery, the compound droplet hydrodynamics in confined geometries have been barely surveyed. For the first time, this study thoroughly investigates the hydrodynamics of a ferrofluid compound droplet (FCD) during its migration in a microchannel under the presence of a pressure-driven flow and a uniform external magnetic field (UEMF) to manipulate its morphology and retard its breakup. Finite difference and phase-field multiple-relaxation time lattice Boltzmann approaches are coupled to determine the magnetic field and ternary flow system, respectively. First, the influence of the magnetic Bond number ( ${Bo}_m$ ) on the FCD morphology is explored depending on whether the core or shell is ferrofluid when the UEMF is applied along $\alpha =0^\circ$ and $\alpha =90^\circ$ relative to the fluid flow. It is ascertained that imposing the UEMF at $\alpha =0^\circ$ when the shell is ferrofluid can postpone the breakup. Intriguingly, when the core is ferrofluid, strengthening the UEMF enlarges the shell deformation. Afterwards, the effects of the capillary number ( ${Ca}$ ), density ratio, viscosity ratio, radius ratio and surface tension coefficients are scrutinised on the FCD deformation and breakup. The results indicate that augmenting the core-to-shell viscosity and density ratios accelerates the breakup process. Additionally, surface tension between the core and shell suppresses the core deformation. Moreover, increasing the ${Ca}$ intensifies the viscous drag force exerted on the shell, flattening its rear side, which causes a triangular-like configuration. Ultimately, by varying ${Bo}_m$ and ${Ca}$ , five distinct regimes are observed, whose regime map is established.
In this work, a systematic study is carried out concerning the dynamic behaviour of finite-size spheroidal particles in non-isothermal shear flows between parallel plates. The simulations rely on a hybrid … In this work, a systematic study is carried out concerning the dynamic behaviour of finite-size spheroidal particles in non-isothermal shear flows between parallel plates. The simulations rely on a hybrid method combining the lattice Boltzmann method with a finite-difference solver. Fluid–particle and heat–particle interactions are accounted for by using the immersed boundary method. The effect of particle Reynolds number ( $\textit{Re}_p=1{-}90$ ), Grashof number ( ${Gr}=0{-}200$ ), initial position and initial orientation of the particle are thoroughly examined. For the isothermal prolate particle, we observed that above a certain Reynolds number, the particle undergoes a pitchfork bifurcation; at an even higher Reynolds number, it returns to the centre position. In contrast, the hot particle behaves differently, with no pitchfork bifurcation. Instead, the Reynolds and Grashof numbers can induce oscillatory tumbling or log-rolling motions in either the lower or upper half of the channel. Heat transfer also plays an important role: at low Grashof numbers, the particle settles near the lower wall, while increasing the Grashof number shifts it towards the upper side. Moreover, the presence of thermal convection increases the rotational speed of the particle. Surprisingly, beyond the first critical Reynolds number, the equilibrium position of the thermal particle shifts closer to the centreline compared with that of a neutrally buoyant isothermal particle. Moreover, higher Grashof numbers can cause the particle to transition from tumbling to log-rolling or even a no-rotation mode. The initial orientation has a stronger influence at low Grashof numbers, while the initial position shows no strong effect in non-isothermal cases.
Microfluidic devices offer well-defined physical environments that are suitable for effective cell seeding and in vitro three-dimensional (3D) cell culture experiments. These platforms have been employed to model in vivo … Microfluidic devices offer well-defined physical environments that are suitable for effective cell seeding and in vitro three-dimensional (3D) cell culture experiments. These platforms have been employed to model in vivo conditions for studying mechanical forces, cell–extracellular matrix (ECM) interactions, and to elucidate transport mechanisms in 3D tissue-like structures, such as tumor and lymph node organoids. Studies have shown that fluid flow behavior in microfluidic slides (µ-slides) directly influences shear stress, which has emerged as a key factor affecting cell proliferation and differentiation. This study investigates fluid flow in the porous channel of a µ-slide using computational fluid dynamics (CFD) techniques to analyze the impact of perfusion flow rate and porous properties on resulting shear stresses. The model of the µ-slide filled with a permeable biomaterial is considered. Porous media fluid flow in the channel is characterized by adding a momentum loss term to the standard Navier–Stokes equations, with a physiological range of permeability values. Numerical simulations are conducted to obtain data and contour plots of the filtration velocity and flow-induced shear stress distributions within the device channel. The filtration flow is subsequently measured by performing protein perfusions into the slide embedded with native human-derived ECM, while the flow rate is controlled using a syringe pump. The relationships between inlet flow rate and shear stress, as well as filtration flow and ECM permeability, are analyzed. The findings provide insights into the impact of shear stress, informing the optimization of perfusion conditions for studying tissues and cells under fluid flow.
The dramatic increase in domestic and industrial waste over recent centuries has significantly polluted water bodies, threatening aquatic life and human activities such as drinking, recreation, and commerce. Understanding pollutant … The dramatic increase in domestic and industrial waste over recent centuries has significantly polluted water bodies, threatening aquatic life and human activities such as drinking, recreation, and commerce. Understanding pollutant dispersion is essential for designing effective waste management systems, employing both experimental and computational techniques. Among Computational Fluid Dynamics (CFD) techniques, the Lattice Boltzmann Method (LBM) has emerged as a novel approach based on a discretized Boltzmann equation. The versatility and parallelization capability of this method makes it particularly attractive for fluid dynamics simulations using high-performance computing. Motivated by its successful application across various scientific disciplines, this study explores the potential of LBM to model pollutant mixing and dilution from outfalls into surface water bodies, focusing specifically on vertical dense jets in crossflow (JICF), a key scenario for the diffusion of brine from desalination plants. A full-LBM scheme is employed to model both the hydrodynamics and the transport of the saline concentration field, and Large Eddy Simulations (LES) are employed in the framework of LBM to reduce computational costs typically associated with turbulence modeling, together with a recursive regularization procedure for the collision operator to achieve greater stability. Several key aspects of vertical dense JICF are considered. The simulations successfully capture general flow characteristics corresponding to jets with varying crossflow parameter urF and most of the typical vortical structures associated with JICF. Relevant quantities such as the terminal rise height, the impact distance, the dilution at the terminal rise height, and the dilution at the impact point are compared with experimental results and semi-empirical relations. The results show a systematic underestimation of these quantities, but the key trends are successfully captured, highlighting LBM’s promise as a tool for simulating wastewater dispersion in aquatic environments.
<title>Abstract</title> We introduce neural network architectures that combine physics-informed neural networks (PINNs) with the immersed boundary method (IBM) to solve fluid-structure interaction (FSI) problems. Our approach features two distinct architectures: … <title>Abstract</title> We introduce neural network architectures that combine physics-informed neural networks (PINNs) with the immersed boundary method (IBM) to solve fluid-structure interaction (FSI) problems. Our approach features two distinct architectures: a Single-FSI network with a unified parameter space, and an innovative Eulerian-Lagrangian network that maintains separate parameter spaces for fluid and structure domains. We study each architecture using standard Tanh and adaptive B-spline activation functions. Empirical studies on a 2D cavity flow problem involving a moving solid structure show that the Eulerian-Lagrangian architecture performs significantly better. The adaptive B-spline activation further enhances accuracy by providing locality-aware representation near boundaries. While our methodology shows promising results in predicting the velocity field, pressure recovery remains challenging due to the absence of explicit force-coupling constraints in the current formulation. Our findings underscore the importance of domain-specific architectural design and adaptive activation functions for modeling FSI problems within the PINN framework.
Abhay | INTERANTIONAL JOURNAL OF SCIENTIFIC RESEARCH IN ENGINEERING AND MANAGEMENT
Abstract—Machine learning (ML) models deployed in dy- namic, real-world environments are susceptible to performance degradation over time due to concept drift—the phenomenon where the underlying data distribution changes. This poses … Abstract—Machine learning (ML) models deployed in dy- namic, real-world environments are susceptible to performance degradation over time due to concept drift—the phenomenon where the underlying data distribution changes. This poses significant challenges to maintaining model reliability and pre- dictive accuracy in production systems. In this project, we propose a fully automated pipeline for drift detection and model retraining, designed to ensure sustained model performance with minimal human intervention. The pipeline leverages statistical drift monitoring techniques through Evidently AI to detect distributional changes in incoming data, generate actionable drift reports, and trigger retraining only when drift is significant and impacts performance. A Boolean logic-based trigger mechanism is used to initiate model retraining using recent data, followed by rigorous evaluation and comparison with the incumbent model. If the retrained model demonstrates improved performance, it is deployed into production using a controlled update strategy. The entire system is modular, scalable, and integrates seamlessly with MLOps workflows. This automated approach not only reduces operational overhead but also enhances model resilience, making it well-suited for applications in real-time analytics, IoT, and adaptive decision systems. Index Terms—Concept drift, model retraining, automated machine learning, drift detection, MLOps, real-time monitoring, adaptive learning systems, Evidently AI, data stream mining, performance-aware retraining, model lifecycle management, edge computing, unsupervised drift detection.
We explore theoretical aspects of boundary conditions for lattice Boltzmann methods, focusing on a toy two-velocities scheme. By mapping lattice Boltzmann schemes to Finite Difference schemes, we facilitate rigorous consistency … We explore theoretical aspects of boundary conditions for lattice Boltzmann methods, focusing on a toy two-velocities scheme. By mapping lattice Boltzmann schemes to Finite Difference schemes, we facilitate rigorous consistency and stability analyses. We develop kinetic boundary conditions for inflows and outflows, highlighting the trade-off between accuracy and stability, which we successfully overcome. Stability is assessed using GKS (Gustafsson, Kreiss, and Sundstr{\"o}m) analysis and -- when this approach fails on coarse meshes -- spectral and pseudo-spectral analyses of the scheme's matrix that explain effects germane to low resolutions.