Source author record

Benjamin Mauroy

Benjamin Mauroy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Biological Physics physics.comp-ph physics.flu-dyn Tissues and Organs math.AP physics.med-ph

Catalog footprint

What is connected

6works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Optimal efficiency of high frequency chest wall oscillations and links with resistance and compliance in a model of the lung

Chest physiotherapy is a set of techniques used to help the draining of the mucus from the lung in pathological situations. The choice of the techniques, and their adjustment to the patients or to the pathologies, remains as of today largely empirical. High Frequency Chest Wall Oscillation (HFCWO) is one of these techniques, performed with a device that applies oscillating pressures on the chest. However, there is no clear understanding of how HFCWO devices interact with the lung biomechanics. Hence, we study idealised HFCWO manipulations applied to a mathematical and numerical model of the biomechanics of the lung. The lung is represented by a fluid--structure interaction model based on an airway tree that is coupled to an homogeneous elastic medium. We show that our model is driven by two dimensionless numbers that drive the effect of the idealised HFCWO manipulation on the model of the lung. Our model allow to analyze the stress applied to an idealised mucus by the air--mucus interaction and by the airway walls deformation. This stress behaves as a buffer and has the effect of reducing the stress needed to overcome the idealised mucus yield stress. Moreover, our model predicts the existence of an optimal range of the working frequencies of HFCWO. This range is in agreement with the frequencies actually used by practitioners during HFCWO maneuvers. Finally, our model suggests that analyzing the mouth airflow during HFCWO maneuvers could allow to estimate the compliance and the hydrodynamic resistance of the lung of a patient.

preprint2015arXiv

Murray's law revisited: Quémada's fluid model and fractal trees

In 1926, Murray proposed the first law for the optimal design of blood vessels. He minimized the power dissipation arising from the trade-off between fluid circulation and blood maintenance. The law, based on a constant fluid viscosity, states that in the optimal configuration the fluid flow rate inside the vessel is proportional to the cube of the vessel radius, implying that wall shear stress is not dependent on the vessel radius. Murray's law has been found to be true in blood macrocirculation, but not in microcirculation. In 2005, Alarcón et al took into account the non monotonous dependence of viscosity on vessel radius - Fåhræus-Lindqvist effect - due to phase separation effect of blood. They were able to predict correctly the behavior of wall shear stresses in microcirculation. One last crucial step remains however: to account for the dependence of blood viscosity on shear rates. In this work, we investigate how viscosity dependence on shear rate affects Murray's law. We extended Murray's optimal design to the whole range of Quémada's fluids, that models pseudo-plastic fluids such as blood. Our study shows that Murray's original law is not restricted to Newtonian fluids, it is actually universal for all Quémada's fluid as long as there is no phase separation effect. When phase separation effect occurs, then we derive an extended version of Murray's law. Our analyses are very general and apply to most of fluids with shear dependent rheology. Finally, we study how these extended laws affect the optimal geometries of fractal trees to mimic an idealized arterial network.

preprint2014arXiv

Functional optimization of the arterial network

We build an evolutionary scenario that explains how some crucial physiological constraints in the arterial network of mammals - i.e. hematocrit, vessels diameters and arterial pressure drops - could have been selected by evolution. We propose that the arterial network evolved while being constrained by its function as an organ. To support this hypothesis, we focus our study on one of the main function of blood network: oxygen supply to the organs. We consider an idealized organ with a given oxygen need and we optimize blood network geometry and hematocrit with the constraint that it must fulfill the organ oxygen need. Our model accounts for the non-Newtonian behavior of blood, its maintenance cost and Fåhræus effects (decrease in average concentration of red blood cells as the vessel diameters decrease). We show that the mean shear rates (relative velocities of fluid layers) in the tree vessels follow a scaling law related to the multi-scale property of the tree network, and we show that this scaling law drives the behavior of the optimal hematocrit in the tree. We apply our scenario to physiological data and reach results fully compatible with the physiology: we found an optimal hematocrit of 0.43 and an optimal ratio for diameter decrease of about 0.79. Moreover our results show that pressure drops in the arterial network should be regulated in order for oxygen supply to remain optimal, suggesting that the amplitude of the arterial pressure drop may have co-evolved with oxygen needs.

preprint2012arXiv

The camera method, or how to track numerically a deformable particle moving in a fluid network

The goal of this work is to follow the displacement and possible deformation of a free particle in a fluid flow in 2D axi-symmetry, 2D or 3D using the classical finite elements method without the usual drawbacks finite elements bring for fluid-structure interaction, i.e. huge numerical problems and strong mesh distortions. Working with finite elements is a choice motivated by the fact that finite elements are well known by a large majority of researchers and are easy to manipulate. The method we describe in this paper, called the camera method, is well adapted to the study of a single particle in a network and most particularly when the study focuses on the particle behaviour. The camera method is based on two principles: 1/ the fluid structure interaction problem is restricted to a neighbourhood of the particle, thus reducing drastically the number of degrees of freedom of the problem; 2/ the neighbourhood mesh moves and rotates with the particle, thus avoiding most of the mesh distortions that occur in a standard ALE method. In this article, we present the camera method and the conditions under which it can be used. Then we apply it to several examples from the literature in 2D axi-symmetry, 2D and 3D.

preprint2010arXiv

Shape minimization of the dissipated energy in dyadic trees

In this paper, we study the role of boundary conditions on the optimal shape of a dyadic tree in which flows a Newtonian fluid. Our optimization problem consists in finding the shape of the tree that minimizes the viscous energy dissipated by the fluid with a constrained volume, under the assumption that the total flow of the fluid is conserved throughout the structure. These hypotheses model situations where a fluid is transported from a source towards a 3D domain into which the transport network also spans. Such situations could be encountered in organs like for instance the lungs and the vascular networks. Two fluid regimes are studied: (i) low flow regime (Poiseuille) in trees with an arbitrary number of generations using a matricial approach and (ii) non linear flow regime (Navier-Stokes, moderate regime with a Reynolds number 100) in trees of two generations using shape derivatives in an augmented Lagrangian algorithm coupled with a 2D/3D finite elements code to solve Navier-Stokes equations. It relies on the study of a finite dimensional optimization problem in the case (i) and on a standard shape optimization problem in the case (ii). We show that the behaviours of both regimes are very similar and that the optimal shape is highly dependent on the boundary conditions of the fluid applied at the leaves of the tree.

preprint2007arXiv

Following red blood cells in a pulmonary capillary

The red blood cells or erythrocytes are biconcave shaped cells and consist mostly in a membrane delimiting a cytosol with a high concentration in hemoglobin. This membrane is highly deformable and allows the cells to go through narrow passages like the capillaries which diameters can be much smaller than red blood cells one. They carry oxygen thanks to hemoglobin, a complex molecule that have very high affinity for oxygen. The capacity of erythrocytes to load and unload oxygen is thus a determinant factor in their efficacy. In this paper, we will focus on the pulmonary capillary where red blood cells capture oxygen. We propose a camera method in order to numerically study the behavior of the red blood cell along a whole capillary. Our goal is to understand how erythrocytes geometrical changes along the capillary can affect its capacity to capture oxygen. The first part of this document presents the model chosen for the red blood cells along with the numerical method used to determine and follow their shapes along the capillary. The membrane of the red blood cell is complex and has been modelled by an hyper-elastic approach coming from Mills et al (2004). This camera method is then validated and confronted with a standard ALE method. Some geometrical properties of the red blood cells observed in our simulations are then studied and discussed. The second part of this paper deals with the modeling of oxygen and hemoglobin chemistry in the geometries obtained in the first part. We have implemented a full complex hemoglobin behavior with allosteric states inspired from Czerlinski et al (1999).