-
Palabos-npFEM: Software for the Simulation of Cellular Blood Flow (Digital Blood)
Authors:
Christos Kotsalos,
Jonas Latt,
Bastien Chopard
Abstract:
Palabos-npFEM is a computational framework for the simulation of blood flow with fully resolved constituents. The software resolves the trajectories and deformed state of blood cells, such as red blood cells and platelets, and the complex interaction between them. The tool combines the lattice Boltzmann solver Palabos for the simulation of blood plasma (fluid phase), a finite element method (FEM)…
▽ More
Palabos-npFEM is a computational framework for the simulation of blood flow with fully resolved constituents. The software resolves the trajectories and deformed state of blood cells, such as red blood cells and platelets, and the complex interaction between them. The tool combines the lattice Boltzmann solver Palabos for the simulation of blood plasma (fluid phase), a finite element method (FEM) solver for the resolution of blood cells (solid phase), and an immersed boundary method (IBM) for the coupling of the two phases. Palabos-npFEM provides, on top of a CPU-only version, the option to simulate the deformable bodies on GPUs, thus the code is tailored for the fastest supercomputers. The software is integrated in the Palabos core library, and is available on the Git repository https://gitlab.com/unigespc/palabos. It offers the possibility to simulate various setups, e.g. several geometries and blood parameters, and due to its modular design, it allows external solvers to readily replace the provided ones.
△ Less
Submitted 9 November, 2020;
originally announced November 2020.
-
Personalized pathology test for Cardio-vascular disease: Approximate Bayesian computation with discriminative summary statistics learning
Authors:
Ritabrata Dutta,
Karim Zouaoui-Boudjeltia,
Christos Kotsalos,
Alexandre Rousseau,
Daniel Ribeiro de Sousa,
Jean-Marc Desmet,
Alain Van Meerhaeghe,
Antonietta Mira,
Bastien Chopard
Abstract:
Cardio/cerebrovascular diseases (CVD) have become one of the major health issue in our societies. But recent studies show that the present pathology tests to detect CVD are ineffectual as they do not consider different stages of platelet activation or the molecular dynamics involved in platelet interactions and are incapable to consider inter-individual variability. Here we propose a stochastic pl…
▽ More
Cardio/cerebrovascular diseases (CVD) have become one of the major health issue in our societies. But recent studies show that the present pathology tests to detect CVD are ineffectual as they do not consider different stages of platelet activation or the molecular dynamics involved in platelet interactions and are incapable to consider inter-individual variability. Here we propose a stochastic platelet deposition model and an inferential scheme to estimate the biologically meaningful model parameters using approximate Bayesian computation with a summary statistic that maximally discriminates between different types of patients. Inferred parameters from data collected on healthy volunteers and different patient types help us to identify specific biological parameters and hence biological reasoning behind the dysfunction for each type of patients. This work opens up an unprecedented opportunity of personalized pathology test for CVD detection and medical treatment.
△ Less
Submitted 9 February, 2022; v1 submitted 13 October, 2020;
originally announced October 2020.
-
Anomalous Platelet Transport & Fat-Tailed Distributions
Authors:
Christos Kotsalos,
Karim Zouaoui Boudjeltia,
Ritabrata Dutta,
Jonas Latt,
Bastien Chopard
Abstract:
The transport of platelets in blood is commonly assumed to obey an advection-diffusion equation. Here we propose a disruptive view, by showing that the random part of their velocity is governed by a fat-tailed probability distribution, usually referred to as a Lévy flight. Although for small spatio-temporal scales, it is hard to distinguish it from the generally accepted "red blood cell enhanced"…
▽ More
The transport of platelets in blood is commonly assumed to obey an advection-diffusion equation. Here we propose a disruptive view, by showing that the random part of their velocity is governed by a fat-tailed probability distribution, usually referred to as a Lévy flight. Although for small spatio-temporal scales, it is hard to distinguish it from the generally accepted "red blood cell enhanced" Brownian motion, for larger systems this effect is dramatic as the standard approach may underestimate the flux of platelets by several orders of magnitude, compromising in particular the validity of current platelet function tests.
△ Less
Submitted 16 September, 2021; v1 submitted 21 June, 2020;
originally announced June 2020.
-
Digital Blood in Massively Parallel CPU/GPU Systems for the Study of Platelet Transport
Authors:
Christos Kotsalos,
Jonas Latt,
Joel Beny,
Bastien Chopard
Abstract:
We propose a highly versatile computational framework for the simulation of cellular blood flow focusing on extreme performance without compromising accuracy or complexity. The tool couples the lattice Boltzmann solver Palabos for the simulation of the blood plasma, a novel finite element method (FEM) solver for the resolution of the deformable blood cells, and an immersed boundary method for the…
▽ More
We propose a highly versatile computational framework for the simulation of cellular blood flow focusing on extreme performance without compromising accuracy or complexity. The tool couples the lattice Boltzmann solver Palabos for the simulation of the blood plasma, a novel finite element method (FEM) solver for the resolution of the deformable blood cells, and an immersed boundary method for the coupling of the two phases. The design of the tool supports hybrid CPU-GPU executions (fluid, fluid-solid interaction on CPUs, the FEM solver on GPUs), and is non-intrusive, as each of the three components can be replaced in a modular way. The FEM-based kernel for solid dynamics outperforms other FEM solvers and its performance is comparable to the state-of-the-art mass-spring systems. We perform an exhaustive performance analysis on Piz Daint at the Swiss National Supercomputing Centre and provide case studies focused on platelet transport. The tests show that this versatile framework combines unprecedented accuracy with massive performance, rendering it suitable for the upcoming exascale architectures.
△ Less
Submitted 8 November, 2019;
originally announced November 2019.
-
Bridging the computational gap between mesoscopic and continuum modeling of red blood cells for fully resolved blood flow
Authors:
Christos Kotsalos,
Jonas Latt,
Bastien Chopard
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…
▽ More
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.
△ Less
Submitted 27 August, 2019; v1 submitted 15 March, 2019;
originally announced March 2019.