Bridging the computational gap between mesoscopic and continuum modeling of red blood cells for fully resolved blood flow

DOI10.1016/J.JCP.2019.108905MaRDI QIDQ2222595zbMATH OpenOpenAlexFDO

Authors Christos Kotsalos, J. Latt, B. Chopard

Publication date 27 January 2021

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1903.06479

zbMATH Keywords

extended projective dynamicsfully resolved blood flowimmersed boundarylattice Boltzmann methodnodal projective FEMred blood cells

Mathematics Subject Classification ID

Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.) (74F10) Finite element methods applied to problems in fluid mechanics (76M10) Particle methods and lattice-gas methods (76M28) Suspensions (76T20) Physiological flows (76Z05) Physiological flow (92C35)

Abstract: We present a computational framework for the simulation of blood flow with fully resolved red blood cells (RBCs) using a modular approach that consists of a lattice Boltzmann solver for the blood plasma, a novel finite element based solver for the deformable bodies and an immersed boundary method for the fluid-solid interaction. For the RBCs, we propose a nodal projective FEM (npFEM) solver which has theoretical advantages over the more commonly used mass-spring systems (mesoscopic modeling), such as an unconditional stability, versatile material expressivity, and one set of parameters to fully describe the behavior of the body at any mesh resolution. At the same time, the method is substantially faster than other FEM solvers proposed in this field, and has an efficiency that is comparable to the one of mesoscopic models. At its core, the solver uses specially defined potential energies, and builds upon them a fast iterative procedure based on quasi-Newton techniques. For a known material, our solver has only one free parameter that demands tuning, related to the body viscoelasticity. In contrast, state-of-the-art solvers for deformable bodies have more free parameters, and the calibration of the models demands special assumptions regarding the mesh topology, which restrict their generality and mesh independence. We propose as well a modification to the potential energy proposed by Skalak et al. 1973 for the red blood cell membrane, which enhances the strain hardening behavior at higher deformations. Our viscoelastic model for the red blood cell, while simple enough and applicable to any kind of solver as a post-convergence step, can capture accurately the characteristic recovery time and tank-treading frequencies. The framework is validated using experimental data, and it proves to be scalable for multiple deformable bodies.

Recommendations

• A new two-component approach in modeling red blood cells
• The finite element immersed boundary method for the numerical simulation of the motion of red blood cells in microfluidic flows
• Coupling the lattice-Boltzmann and spectrin-link methods for the direct numerical simulation of cellular blood flow
• A lattice Boltzmann fictitious domain method for modeling red blood cell deformation and multiple-cell hydrodynamic interactions in flow
• Multiscale modelling of erythrocytes in Stokes flow

Cites work

• A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems
• An efficient multilevel solution scheme for large scale non‐linear systems
• Coupling the lattice-Boltzmann and spectrin-link methods for the direct numerical simulation of cellular blood flow
• Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method
• Flow patterns around heart valves: A numerical method
• scientific article; zbMATH DE number 3144204 (Why is no real title available?)
• scientific article; zbMATH DE number 5060482 (Why is no real title available?)
• IMMERSED BOUNDARY METHODS
• Lift generation by a two-dimensional symmetric flapping wing: immersed boundary-lattice Boltzmann simulations
• Nonlinear Continuum Mechanics for Finite Element Analysis
• Numerical models for differential problems. Translated by Silvia Quarteroni.
• On the bending algorithms for soft objects in flows
• Projective dynamics: fusing constraint projections for fast simulation
• Simulating deformable particle suspensions using a coupled lattice-Boltzmann and finite-element method
• Some practical procedures for the solution of nonlinear finite element equations
• The immersed boundary method
• The lattice Boltzmann method. Principles and practice

Cited in

(10)

• A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows
• A curvilinear lattice Boltzmann scheme for thermal flows
• Palabos: parallel lattice Boltzmann solver
• Parameterizing the Morse potential for coarse-grained modeling of blood plasma
• An Improved Coarse-Grained Model to Accurately Predict Red Blood Cell Morphology and Deformability
• Efficient supersonic flow simulations using lattice Boltzmann methods based on numerical equilibria
• npFEM
• Heat transfer and fluid flow around starfish-inspired cylinder by immersed boundary-simplified thermal lattice Boltzmann method
• Eulerian formulation of the tensor-based morphology equations for strain-based blood damage modeling
• Determination of platelets properties in blood: combining discrete numerical models with experiments

Describes a project that uses

Uses Software

• Hemocell

This page was built for publication: Bridging the computational gap between mesoscopic and continuum modeling of red blood cells for fully resolved blood flow

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222595)