Chromodynamic lattice Boltzmann method for the simulation of drops, erythrocytes, and other vesicles

SPENDLOVE, James, XU, Xu, SCHENKEL, Torsten, GUNN, Julian and HALLIDAY, Ian (2023). Chromodynamic lattice Boltzmann method for the simulation of drops, erythrocytes, and other vesicles. Communications in computational physics, 33 (1), 283-309.

Chromodynamic_Lattice_Boltzmann_Method_for_the_Simulation_of_Drops__Erythr.pdf - Accepted Version
All rights reserved.

Download (2MB) | Preview
Official URL:
Link to published version::


Recently, we have validated a three-dimensional, single framework multicomponent lattice Boltzmann method, modified to generate vesicles (rather than drops) [“Three-dimensional single framework multicomponent lattice Boltzmann equation method for vesicle hydrodynamics,” Phys. Fluids 33, 077110 (2021)]. This approach implements an immersed boundary force distribution, characterised by bending rigidity, surface tension, preferred curvature and conserved membrane area, in which work we successfully validated isolated vesicle flows against other methodologies and experiment. Like most immersed boundary algorithms, our method relies on numerical computation of high-order spatial derivatives and an intricate body force density. The next step is to verify that it has sufficient numerical stability to address the anticipated application of high volume fraction flows of highly deformable objects in intimate interaction. It is this in silico verification – of both the class of fluid object attainable and the stability of the later in strong, straining and shearing flows which is at issue, here. We extend our method to simulate multiple variously deflated vesicles and multiple liquid droplets still within a single framework, from which our fluid objects emerge as particular parameterisations. We present data from simulations containing up to four vesicles (five immiscible fluid species), which threshold verifies that simulations containing unlimited fluid objects are possible [“Modeling the flow of dense suspensions of deformable particles in three dimensions,” Phys. Rev. E 75, 066707 (2007)]. These data also assure the ability of our immersed boundary forcing to preserve the character and integrity of fluid objects in interactions characterised by large local velocity gradients (intimate squeezing, shearing and elongational straining). Throughout, we take interfacial or membrane area, A, as a proxy for stability and physical veracity.

Item Type: Article
Additional Information: Published online date taken from pdf properties information
Uncontrolled Keywords: Applied Mathematics
Identification Number:
Page Range: 283-309
SWORD Depositor: Symplectic Elements
Depositing User: Symplectic Elements
Date Deposited: 17 Oct 2022 11:04
Last Modified: 17 Feb 2024 01:18

Actions (login required)

View Item View Item


Downloads per month over past year

View more statistics