An efficient, localised approach for the simulation of elastic blood vessels using the lattice Boltzmann method.
Ontology highlight
ABSTRACT:
SUBMITTER: McCullough JWS
PROVIDER: S-EPMC8688478 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description: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.
| S-EPMC4890317 | biostudies-literature
Project description:Porous structures of shales are reconstructed using the markov chain monte carlo (MCMC) method based on scanning electron microscopy (SEM) images of shale samples from Sichuan Basin, China. Characterization analysis of the reconstructed shales is performed, including porosity, pore size distribution, specific surface area and pore connectivity. The lattice Boltzmann method (LBM) is adopted to simulate fluid flow and Knudsen diffusion within the reconstructed shales. Simulation results reveal that the tortuosity of the shales is much higher than that commonly employed in the Bruggeman equation, and such high tortuosity leads to extremely low intrinsic permeability. Correction of the intrinsic permeability is performed based on the dusty gas model (DGM) by considering the contribution of Knudsen diffusion to the total flow flux, resulting in apparent permeability. The correction factor over a range of Knudsen number and pressure is estimated and compared with empirical correlations in the literature. For the wide pressure range investigated, the correction factor is always greater than 1, indicating Knudsen diffusion always plays a role on shale gas transport mechanisms in the reconstructed shales. Specifically, we found that most of the values of correction factor fall in the slip and transition regime, with no Darcy flow regime observed.
| S-EPMC4308705 | biostudies-literature
Project description:In this paper we develop a lattice Boltzmann algorithm to simulate red blood cell (RBC) behavior in shear flows. The immersed boundary method is employed to incorporate the fluid-membrane interaction between the flow field and deformable cells. The cell membrane is treated as a neo-Hookean viscoelastic material and a Morse potential is adopted to model the intercellular interaction. Utilizing the available mechanical properties of RBCs, multiple cells have been studied in shear flows using a two-dimensional approximation. These cells aggregate and form a rouleau under the action of intercellular interaction. The equilibrium configuration is related to the interaction strength. The end cells exhibit concave shapes under weak interaction and convex shapes under strong interaction. In shear flows, such a rouleau-like aggregate will rotate or be separated, depending on the relative strengths of the intercellular interaction and hydrodynamic viscous forces. These behaviors are qualitatively similar to experimental observations and show the potential of this numerical scheme for future studies of blood flow in microvessels.
| S-EPMC2254372 | biostudies-literature
Project description:Analyzing the fluid behavior in complex porous media like gas diffusion layers (GDLs) in polymer electrolyte fuel cells (PEFCs) can be accurately done using the lattice Boltzmann method (LBM). This article shows the data obtained from a study in which diffusion parameters such as porosity, gas phase tortuosity and diffusibility are computed considering simulated porous media [1]. The data were computed when a water drop obstacle is placed inside the GDL domain and the size of the water-drop is varied. Additionally, figures showing the evolution of the flow velocity field are presented alongside graphics that presents the change in local and bulk porosity for each obstacle size. Finally, there is a detailed method explanation concerning the implementation of the lattice Boltzmann method and a general description of computational codes for the domain and obstacle generation as well as the boundary conditions simulation. Data and processes in this article can be exploited in new attempts to solve real case problems in complex mesoscale media.
| S-EPMC6838375 | biostudies-literature
Project description:Nanodrug-carrier delivery in the blood stream is strongly influenced by nanoparticle (NP) dispersion. This paper presents a numerical study on NP transport and dispersion in red blood cell (RBC) suspensions under shear and channel flow conditions, utilizing an immersed boundary fluid-structure interaction model with a lattice Boltzmann fluid solver, an elastic cell membrane model and a particle motion model driven by both hydrodynamic loading and Brownian dynamics. The model can capture the multiphase features of the blood flow. Simulations were performed to obtain an empirical formula to predict NP dispersion rate for a range of shear rates and cell concentrations. NP dispersion rate predictions from the formula were then compared to observations from previous experimental and numerical studies. The proposed formula is shown to accurately predict the NP dispersion rate. The simulation results also confirm previous findings that the NP dispersion rate is strongly influenced by local disturbances in the flow due to RBC motion and deformation. The proposed formula provides an efficient method for estimating the NP dispersion rate in modeling NP transport in large-scale vascular networks without explicit RBC and NP models.
| S-EPMC5302481 | biostudies-literature
Project description:The three-dimensional multiple-relaxation-time LB (MRT-LB) and Lagrangian particle tracking methods were applied to simulate turbulent airflow and particle dispersion in a ventilated room with a partition. The turbulent airflow was simulated by large eddy simulation (LES) using the MRT-LB method with the Smagorinsky model. This method was verified by comparing it with experimental and other numerical results. Good agreement was observed between airflow simulation and experimental data. It is also demonstrated that the LES carried out by the MRT-LB method can produce airflow results very similar to the RNG LES and provide better prediction than the standard and RNG k-ε models. In order to further improve the efficiency of the MRT-LB method, the multi-block grid refinement (MBGR) technique was used. The accuracy and efficiency of the MBGR and the consistency of physical quantities across the interface were investigated. In simulation of particle dispersion in the model room, particles with diameters of 1 and 10 µm were considered. It is shown that this model can successfully capture dispersion characteristics of particles.
| S-EPMC7127572 | biostudies-literature
Project description:Immersed interface method (IIM) is a promising high-accuracy numerical scheme for the Poisson-Boltzmann model that has been widely used to study electrostatic interactions in biomolecules. However, the IIM suffers from instability and slow convergence for typical applications. In this study, we introduced both analytical interface and surface regulation into IIM to address these issues. The analytical interface setup leads to better accuracy and its convergence closely follows a quadratic manner as predicted by theory. The surface regulation further speeds up the convergence for nontrivial biomolecules. In addition, uncertainties of the numerical energies for tested systems are also reduced by about half. More interestingly, the analytical setup significantly improves the linear solver efficiency and stability by generating more precise and better-conditioned linear systems. Finally, we implemented the bottleneck linear system solver on GPUs to further improve the efficiency of the method, so it can be widely used for practical biomolecular applications. © 2019 Wiley Periodicals, Inc.
| S-EPMC6422926 | biostudies-literature
Project description:A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on a numerical equilibrium that has been extensively used by discrete velocity methods since the late 1990s. Here, it is extended to reproduce an arbitrary number of moments of the Maxwell-Boltzmann distribution. These extensions to the standard 5-constraint (mass, momentum and energy) approach lead to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high Reynolds-number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the supersonic flow past a NACA0012 aerofoil confirms the excellent behaviour of this model in a low-viscosity and supersonic regime. The flow past a sphere is further simulated to investigate the 3D behaviour of our model in the low-viscosity supersonic regime. The proposed model is shown to be substantially more efficient than the previous 5-moment D3Q343 DDF-LBM for both CPU and GPU architectures. It then opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.
| S-EPMC7333948 | biostudies-literature
Project description:Permeability is a key parameter for investigating the flow ability of sedimentary rocks. The conventional model for calculating permeability is derived from Darcy's law, which is valid only for continuum flow in porous rocks. We discussed the feasibility of simulating methane transport characteristics in the organic nano-pores of shale through the Lattice Boltzmann method (LBM). As a first attempt, the effects of high Knudsen number and the associated slip flow are considered, whereas the effect of adsorption in the capillary tube is left for future work. Simulation results show that at small Knudsen number, LBM results agree well with Poiseuille's law, and flow rate (flow capacity) is proportional to the square of the pore scale. At higher Knudsen numbers, the relaxation time needs to be corrected. In addition, velocity increases as the slip effect causes non negligible velocities on the pore wall, thereby enhancing the flow rate inside the pore, i.e., the permeability. Therefore, the LBM simulation of gas flow characteristics in organic nano-pores provides an effective way of evaluating the permeability of gas-bearing shale.
| S-EPMC4007072 | biostudies-other
Project description:The methane desorption and diffusion characteristics in coal are key factors affecting coalbed methane productivity. In this paper, we developed a lattice Boltzmann model for methane migration in the multiscale porous media of coal. In the simulation, the diffusion of methane in macropores/fractures is assumed to follow Fick's law, and that in the coal matrix is treated as Knudsen diffusion. In addition, the Langmuir adsorption kinetics equation is employed to describe the dynamic process of methane adsorption and desorption. The results indicated the following: (1) The specific surface area and fracture proportion of the coal will increase with the employment of hydraulic fracturing, which may prompt the gas desorption-diffusion efficiency. (2) The flow and diffusion of methane are closely related to each other. When the gas diffusivity is poor, the desorption-diffusion can be effectively accelerated by increasing the drainage intensity, but when the gas diffusivity is fine, the flow velocity has little influence on the methane desorption. In practice, if the estimated methane diffusion coefficient is below the order of 10-5 m2/s, more attention should be paid to its accuracy; otherwise, the obtained results may have a large deviation from the real value. (3) In the typical range of average pore sizes of coal, gas desorption rate growth with the increase of pore size makes the low-rank coal more advantageous in exploitation due to its larger average pore size. With the decline of reservoir pressure, the low- and high-rank coals more easily desorb methane than medium-rank coal. (4) In the kinetic study of the coalbed methane desorption-diffusion process, the accuracy of the obtained results may depend on the adsorption and desorption rate constants if the desorption rate constant is less than 106 1/s.
| S-EPMC8340423 | biostudies-literature