Micro-objets déformables biomimétiques sous écoulement
Dynamique des fluides et transfert en microgravité
Axe de recherche :
Micro-objets déformables sous forçage hydrodynamique (voir les détails sur la page dédiée)
Publications scientifiques au M2P2
Sudip Das, Marc Jaeger, Marc Leonetti, Rochish M. Thaokar, Paul G. Chen. Effect of pulse width on the dynamics of a deflated vesicle in unipolar and bipolar pulsed electric fields. Physics of Fluids, American Institute of Physics, 2021, 33 (8), pp.081905. ⟨10.1063/5.0057168⟩. ⟨hal-03317441⟩ Plus de détails...
Giant unilamellar vesicles subjected to pulsed direct-current (pulsed-DC) fields are promising biomimetic systems to investigate the electroporation of cells. In strong electric fields, vesicles undergo significant deformation, which strongly alters the transmembrane potential, consequently the electroporation. Previous theoretical studies investigated the electrodeformation of vesicles in DC fields (which are not pulsed). In this work, we computationally investigate the deformation of a deflated vesicle under unipolar, bipolar, and two-step unipolar pulses and show sensitive dependence of intermediate shapes on type of pulse and the pulse width. Starting with the stress-free initial shape of a deflated vesicle, which is similar to a prolate spheroid, the analysis is presented for the cases with higher and lower conductivities of the inner fluid medium relative to the outer fluid medium. For the ratio of inner to outer fluid conductivity, σ r = 10, the shape always remains prolate, including when the field is turned off. For σ r = 0.1, several complex dynamics are observed, such as the prolate-to-oblate (PO), prolate-to-oblate-to-prolate (POP) shape transitions in time depending upon the strength of the field and the pulse properties. In this case, on turning off the field, a metastable oblate equilibrium shape is seen, that seems to be a characteristics of a deflated vesicle leading to POPO transitions. When a two-step unipolar pulse (a combination of a strong and a weak subpulse) is applied, a vesicle can reach an oblate or a prolate final shape depending upon the relative durations of the two subpulses. This study suggests that the transmembrane potential can be regulated using a bipolar pulsed-DC field. It also shows that the shapes admitted in the dynamics of a vesicle depends upon whether the pulse is unipolar or bipolar. Parameters are suggested wherein, the simulation results can be demonstrated in experiments.
Sudip Das, Marc Jaeger, Marc Leonetti, Rochish M. Thaokar, Paul G. Chen. Effect of pulse width on the dynamics of a deflated vesicle in unipolar and bipolar pulsed electric fields. Physics of Fluids, American Institute of Physics, 2021, 33 (8), pp.081905. ⟨10.1063/5.0057168⟩. ⟨hal-03317441⟩
Jinming Lyu, Paul G. Chen, G. Boedec, M. Leonetti, Marc Jaeger. An isogeometric boundary element method for soft particles flowing in microfluidic channels. Computers and Fluids, Elsevier, 2021, 214, pp.104786. ⟨10.1016/j.compfluid.2020.104786⟩. ⟨hal-02476945v2⟩ Plus de détails...
Understanding the flow of deformable particles such as liquid drops, synthetic capsules and vesicles, and biological cells confined in a small channel is essential to a wide range of potential chemical and biomedical engineering applications. Computer simulations of this kind of fluid-structure (mem-brane) interaction in low-Reynolds-number flows raise significant challenges faced by an intricate interplay between flow stresses, complex particles' in-terfacial mechanical properties, and fluidic confinement. Here, we present an isogeometric computational framework by combining the finite-element method (FEM) and boundary-element method (BEM) for an accurate prediction of the deformation and motion of a single soft particle transported in microfluidic channels. The proposed numerical framework is constructed consistently with the isogeometric analysis paradigm; Loop's subdivision elements are used not only for the representation of geometry but also for the membrane mechanics solver (FEM) and the interfacial fluid dynamics solver (BEM). We validate our approach by comparison of the simulation results with highly accurate benchmark solutions to two well-known examples available in the literature, namely a liquid drop with constant surface tension in a circular tube and a capsule with a very thin hyperelastic membrane in a square channel. We show that the numerical method exhibits second-order convergence in both time and space. To further demonstrate the accuracy and long-time numerically stable simulations of the algorithm, we perform hydrodynamic computations of a lipid vesicle with bending stiffness and a red blood cell with a composite membrane in capillaries. The present work offers some possibilities to study the deformation behavior of confining soft particles, especially the particles' shape transition and dynamics and their rheological signature in channel flows.
Jinming Lyu, Paul G. Chen, G. Boedec, M. Leonetti, Marc Jaeger. An isogeometric boundary element method for soft particles flowing in microfluidic channels. Computers and Fluids, Elsevier, 2021, 214, pp.104786. ⟨10.1016/j.compfluid.2020.104786⟩. ⟨hal-02476945v2⟩
Paul G. Chen, J Lyu, M Jaeger, M. Leonetti. Shape transition and hydrodynamics of vesicles in tube flow. Physical Review Fluids, American Physical Society, 2020, 5 (4), pp.043602. ⟨10.1103/PhysRevFluids.5.043602⟩. ⟨hal-02415320v2⟩ Plus de détails...
The steady motion and deformation of a lipid-bilayer vesicle translating through a circular tube in low Reynolds number pressure-driven flow are investigated numerically using an axisymmetric boundary element method. This fluid-structure interaction problem is determined by three dimen-sionless parameters: reduced volume (a measure of the vesicle asphericity), geometric confinement (the ratio of the vesicle effective radius to the tube radius), and capillary number (the ratio of viscous to bending forces). The physical constraints of a vesicle--fixed surface area and enclosed volume when it is confined in a tube--determine critical confinement beyond which it cannot pass through without rupturing its membrane. The simulated results are presented in a wide range of reduced volumes [0.6, 0.98] for different degrees of confinement; the reduced volume of 0.6 mimics red blood cells. We draw a phase diagram of vesicle shapes and propose a shape transition line separating the parachutelike shape region from the bulletlike one in the reduced volume versus confinement phase space. We show that the shape transition marks a change in the behavior of vesicle mobility, especially for highly deflated vesicles. Most importantly, high-resolution simulations make it possible for us to examine the hydrodynamic interaction between the wall boundary and the vesicle surface at conditions of very high confinement, thus providing the limiting behavior of several quantities of interest, such as the thickness of lubrication film, vesicle mobility and its length, and the extra pressure drop due to the presence of the vesicle. This extra pressure drop holds implications for the rheology of dilute vesicle suspensions. Furthermore, we present various correlations and discuss a number of practical applications. The results of this work may serve as a benchmark for future studies and help devise tube-flow experiments.
Paul G. Chen, J Lyu, M Jaeger, M. Leonetti. Shape transition and hydrodynamics of vesicles in tube flow. Physical Review Fluids, American Physical Society, 2020, 5 (4), pp.043602. ⟨10.1103/PhysRevFluids.5.043602⟩. ⟨hal-02415320v2⟩
Wenjun Liu, Paul G. Chen, Jalil Ouazzani, Qiusheng Liu. Thermocapillary flow transition in an evaporating liquid layer in a heated cylindrical cell. International Journal of Heat and Mass Transfer, Elsevier, 2020, 153, pp.119587. ⟨10.1016/j.ijheatmasstransfer.2020.119587⟩. ⟨hal-02488876⟩ Plus de détails...
Motived by recent ground-based and microgravity experiments investigating the interfacial dynamics of a volatile liquid (FC-72, P r = 12.34) contained in a heated cylindrical cell, we numerically study the thermocapillary-driven flow in such an evaporating liquid layer. Particular attention is given to the prediction of the transition of the axisymmetric flow to fully three-dimensional patterns when the applied temperature increases. The numerical simulations rely on an improved one-sided model of evaporation by including heat and mass transfer through the gas phase via the heat transfer Biot number and the evaporative Biot number. We present the axisymmetric flow characteristics, show the variation of the transition points with these Biot numbers, and more importantly elucidate the twofold role of the latent heat of evaporation in the stability; evaporation not only destabilizes the flow but also stabilizes it, depending upon the place where the evaporation-induced thermal gradients come into play. We also show that buoyancy in the liquid layer has a stabilizing effect, though its effect is insignificant. At high Marangoni numbers, the numerical simulations revealed smaller-scale thermal patterns formed on the surface of a thinner evaporating layer, in qualitative agreement with experimental observations. The present work helps to gain a better understanding of the role of a phase change in the thermocapillary instability of an evaporating liquid layer.
Wenjun Liu, Paul G. Chen, Jalil Ouazzani, Qiusheng Liu. Thermocapillary flow transition in an evaporating liquid layer in a heated cylindrical cell. International Journal of Heat and Mass Transfer, Elsevier, 2020, 153, pp.119587. ⟨10.1016/j.ijheatmasstransfer.2020.119587⟩. ⟨hal-02488876⟩
Journal: International Journal of Heat and Mass Transfer
Xue Chen, Xun Wang, Paul G. Chen, Qiusheng Liu. Determination of Diffusion Coefficient in Droplet Evaporation Experiment Using Response Surface Method. Microgravity Science and Technology, Springer, 2018, 30, pp.675-682. ⟨10.1007/s12217-018-9645-2⟩. ⟨hal-02112826⟩ Plus de détails...
Evaporation of a liquid droplet resting on a heated substrate is a complex free-surface advection-diffusion problem, in which the main driving force of the evaporation is the vapor concentration gradient across the droplet surface. Given the uncertainty associated with the diffusion coefficient of the vapor in the atmosphere during space evaporation experiments due to the environmental conditions, a simple and accurate determination of its value is of paramount importance for a better understanding of the evaporation process. Here we present a novel approach combining numerical simulations and experimental results to address this issue. Specifically, we construct a continuous function of output using a Kriging-based response surface method, which allows to use the numerical results as a black-box with a limited number of inputs and outputs. Relevant values of the diffusion coefficient can then be determined by solving an inverse problem which is based on accessible experimental data and the proposed response surface. In addition, on the basis of our numerical simulation results, we revisit a widely used formula for the prediction of the evaporation rate in the literature and propose a refined expression for the droplets evaporating on a heated substrate.
Xue Chen, Xun Wang, Paul G. Chen, Qiusheng Liu. Determination of Diffusion Coefficient in Droplet Evaporation Experiment Using Response Surface Method. Microgravity Science and Technology, Springer, 2018, 30, pp.675-682. ⟨10.1007/s12217-018-9645-2⟩. ⟨hal-02112826⟩
Jinming Lyu, Paul G. Chen, Gwenn Boedec, Marc Leonetti, Marc Jaeger. Hybrid continuum–coarse-grained modeling of erythrocytes. Comptes Rendus Mécanique, Elsevier Masson, 2018, 346, pp.439-448. ⟨10.1016/j.crme.2018.04.015⟩. ⟨hal-01785429⟩ Plus de détails...
The red blood cell (RBC) membrane is a composite structure, consisting of a phospholipid bilayer and an underlying membrane-associated cytoskeleton. Both continuum and particle-based coarse-grained RBC models make use of a set of vertices connected by edges to represent the RBC membrane, which can be seen as a triangular surface mesh for the former and a spring network for the latter. Here, we present a modeling approach combining an existing continuum vesicle model with a coarse-grained model for the cytoskeleton. Compared to other two-component approaches, our method relies on only one mesh, representing the cytoskeleton, whose velocity in the tangential direction of the membrane may be different from that of the lipid bilayer. The finitely extensible nonlinear elastic (FENE) spring force law in combination with a repulsive force defined as a power function (POW), called FENE-POW, is used to describe the elastic properties of the RBC membrane. The mechanical interaction between the lipid bilayer and the cytoskeleton is explicitly computed and incorporated into the vesicle model. Our model includes the fundamental mechanical properties of the RBC membrane, namely fluidity and bending rigidity of the lipid bilayer, and shear elasticity of the cytoskeleton while maintaining surface-area and volume conservation constraint. We present three simulation examples to demonstrate the effectiveness of this hybrid continuum--coarse-grained model for the study of RBCs in fluid flows.
Jinming Lyu, Paul G. Chen, Gwenn Boedec, Marc Leonetti, Marc Jaeger. Hybrid continuum–coarse-grained modeling of erythrocytes. Comptes Rendus Mécanique, Elsevier Masson, 2018, 346, pp.439-448. ⟨10.1016/j.crme.2018.04.015⟩. ⟨hal-01785429⟩
Xue Chen, Xun Wang, Paul G. Chen, Qiusheng Liu. Thermal effects of substrate on Marangoni flow in droplet evaporation: Response surface and sensitivity analysis. International Journal of Heat and Mass Transfer, Elsevier, 2017, 113, pp.354 - 365. ⟨10.1016/j.ijheatmasstransfer.2017.05.076⟩. ⟨hal-01532757⟩ Plus de détails...
In this paper, the evaporation of sessile droplets resting on a substrate with different thermal properties is numerically investigated. Computations are based on a transient axisymmetric numerical model. Special attention is paid to evaluate thermal effects of substrate on the structure of bulk fluid flow in the course of evaporation. Numerical results reveal that Marangoni convection induced by non-uniform distribution of temperature along the interface exhibits three distinctly different behaviours: inward flow, multicellular flow and outward flow, consequently resulting in different particle depositions. It is highlighted that three factors (i.e. relative thermal conductivity, relative substrate thickness and relative substrate temperature) strongly affect the flow pattern. In order to further investigate the coupling effects of different influential factors, a Kriging-based response surface method is introduced. We model the flow behaviour as a function of continuous influential factors using a limited number of computations corresponding to discrete values of the inputs. The sensitivities of the Marangoni flow are also analysed using Sobol’ index to study the coupling mechanisms of influential factors. The proposed method can be used to forecast the flow patterns for any input parameter without additional sophisticated computer simulation, and allows to confidently estimate an unknown environmental parameter.
Xue Chen, Xun Wang, Paul G. Chen, Qiusheng Liu. Thermal effects of substrate on Marangoni flow in droplet evaporation: Response surface and sensitivity analysis. International Journal of Heat and Mass Transfer, Elsevier, 2017, 113, pp.354 - 365. ⟨10.1016/j.ijheatmasstransfer.2017.05.076⟩. ⟨hal-01532757⟩
Journal: International Journal of Heat and Mass Transfer
Xue Chen, Paul G. Chen, Jalil Ouazzani, Qiusheng Liu. Numerical simulations of sessile droplet evaporating on heated substrate. The European Physical Journal. Special Topics, EDP Sciences, 2017, 226 (6), pp.1325-1335. ⟨10.1140/epjst/e2016-60203-y⟩. ⟨hal-01509843⟩ Plus de détails...
Motivated by the space project EFILE, a 2D axisymmetric numerical model in the framework of ALE method is developed to investigate the coupled physical mechanism during the evaporation of a pinned drop that partially wets on a heated substrate. The model accounts for mass transport in surrounding air, Marangoni convection inside the drop and heat conduction in the substrate as well as moving interface. Numerical results predict simple scaling laws for the evaporation rate which scales linearly with drop radius but follows a power-law with substrate temperature. It is highlighted that thermal effect of the substrate has a great impact on the temperature profile at the drop surface, which leads to a multicellular thermocapillary flow pattern. In particular, the structure of the multicellular flow behavior induced within a heated drop is mainly controlled by a geometric parameter (aspect ratio). A relationship between the number of thermal cells and the aspect ratio is proposed。
Xue Chen, Paul G. Chen, Jalil Ouazzani, Qiusheng Liu. Numerical simulations of sessile droplet evaporating on heated substrate. The European Physical Journal. Special Topics, EDP Sciences, 2017, 226 (6), pp.1325-1335. ⟨10.1140/epjst/e2016-60203-y⟩. ⟨hal-01509843⟩
Journal: The European Physical Journal. Special Topics
E Alekseenko, B Roux, D Fougere, Paul G. Chen. The effect of wind induced bottom shear stress and salinity on Zostera noltii replanting in a Mediterranean coastal lagoon. Estuarine, Coastal and Shelf Science, Elsevier, 2017, 187, pp.293-305. ⟨10.1016/j.ecss.2017.01.010⟩. ⟨hal-01453377⟩ Plus de détails...
The paper concerns the wind influence on bottom shear stress and salinity levels in a Mediterranean semi-enclosed coastal lagoon (Etang de Berre), with respect to a replanting program of Zostera noltii . The MARS3D numerical model is used to analyze the 3D current, salinity and temperature distribution induced by three meteorological, oceanic and anthropogenic forcings in this lagoon. The numerical model has been carefully validated by comparison with daily observations of the vertical salinity and temperature profiles at three mooring stations, for one year. Then, two modelling scenarios are considered. The first scenario (scen.## 1), starting with an homogeneous salinity of S = 20 PSU and without wind forcing, studies a stratification process under the influence of a periodic seawater inflow and a strong freshwater inflow from an hydropower plant (250 m3/s). Then, in the second scenario (scen.## 2), we study how a strong wind of 80 km/h can mix the haline stratification obtained at the end of scen.## 1. The most interesting results concern four nearshore replanting areas; two are situated on the eastern side of EB and two on the western side. The results of scen.## 2 show that all these areas are subject to a downwind coastal jet. Concerning bottom salinity, the destratification process is very beneficial; it always remains greater than 12 PSU for a N-NW wind of 80 km/h and an hydropower runoff of 250 m3/s. Special attention is devoted to the bottom shear stress (BSS) for different values of the bottom roughness parameter (for gravels, sands and silts), and to the bottom salinity. Concerning BSS, it presents a maximum near the shoreline and decreases along transects perpendicular to the shoreline. There exists a zone, parallel to the shoreline, where BSS presents a minimum (close to zero). When comparing the BSS value at the four replanting areas with the critical value, BSScr, at which the sediment mobility would occur, we see that for the smaller roughness values (ranging from z0 = 3.5 × 10-4 mm, to 3.5 × 10-2 mm) BSS largely surpasses this critical value. For a N-NW wind speed of 40 km/h (which is blowing for around 100 days per year), BSS still largely surpasses BSScr - at least for the silt sediments (ranging from z0 = 3.5 × 10-4 mm, to 3.5 × 10-3 mm). This confirms the possibility that the coastal jet could generate sediment mobility which could have a negative impact for SAV replanting.
E Alekseenko, B Roux, D Fougere, Paul G. Chen. The effect of wind induced bottom shear stress and salinity on Zostera noltii replanting in a Mediterranean coastal lagoon. Estuarine, Coastal and Shelf Science, Elsevier, 2017, 187, pp.293-305. ⟨10.1016/j.ecss.2017.01.010⟩. ⟨hal-01453377⟩
Achim Guckenberger, Marcel P. Schrame, Paul G. Chen, Marc Leonetti, Stephan Gekle. On the bending algorithms for soft objects in flows. Computer Physics Communications, Elsevier, 2016, 207, pp.1-23. ⟨10.1016/j.cpc.2016.04.018⟩. ⟨hal-01314722⟩ Plus de détails...
One of the most challenging aspects in the accurate simulation of three-dimensional soft objects such as vesicles or biological cells is the computation of membrane bending forces. The origin of this difficulty stems from the need to numerically evaluate a fourth order derivative on the discretized surface geometry. Here we investigate six different algorithms to compute membrane bending forces, including regularly used methods as well as novel ones. All are based on the same physical model (due to Canham and Helfrich) and start from a surface discretization with flat triangles. At the same time, they differ substantially in their numerical approach. We start by comparing the numerically obtained mean curvature, the Laplace-Beltrami operator of the mean curvature and finally the surface force density to analytical results for the discocyte resting shape of a red blood cell. We find that none of the considered algorithms converges to zero error at all nodes and that for some algorithms the error even diverges. There is furthermore a pronounced influence of the mesh structure: Discretizations with more irregular triangles and node connectivity present serious difficulties for most investigated methods. To assess the behavior of the algorithms in a realistic physical application, we investigate the deformation of an initially spherical capsule in a linear shear flow at small Reynolds numbers. To exclude any influence of the flow solver, two conceptually very different solvers are employed: the Lattice-Boltzmann and the Boundary Integral Method. Despite the largely different quality of the bending algorithms when applied to the static red blood cell, we find that in the actual flow situation most algorithms give consistent results for both hydrodynamic solvers. Even so, a short review of earlier works reveals a wide scattering of reported results for, e.g., the Taylor deformation parameter. Besides the presented application to biofluidic systems, the investigated algorithms are also of high relevance to the computer graphics and numerical mathematics communities.
Achim Guckenberger, Marcel P. Schrame, Paul G. Chen, Marc Leonetti, Stephan Gekle. On the bending algorithms for soft objects in flows. Computer Physics Communications, Elsevier, 2016, 207, pp.1-23. ⟨10.1016/j.cpc.2016.04.018⟩. ⟨hal-01314722⟩
Bo Xun, Kai Li, Paul G. Chen, Wen-Rui Hu. Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge. International Journal of Heat and Mass Transfer, Elsevier, 2009, 52 (19-20), pp.4211-4220. ⟨10.1016/j.ijheatmasstransfer.2009.04.008⟩. ⟨hal-01307193⟩ Plus de détails...
In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and ''thermal " energy change of the critical disturbances were investigated
Bo Xun, Kai Li, Paul G. Chen, Wen-Rui Hu. Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge. International Journal of Heat and Mass Transfer, Elsevier, 2009, 52 (19-20), pp.4211-4220. ⟨10.1016/j.ijheatmasstransfer.2009.04.008⟩. ⟨hal-01307193⟩
Journal: International Journal of Heat and Mass Transfer
Jian-Kang Zhang, Bo Xun, Paul G. Chen. A Continuation Method Applied to the Study of Thermocapillary Instabilities in Liquid Bridges. Microgravity Science & Technology, 2009, 21 (Suppl 1), pp.111-117. ⟨10.1007/s12217-009-9111-2⟩. ⟨hal-01307220⟩ Plus de détails...
A continuation method is applied to investigate the linear stability of the steady, axisymmetric thermocapillary flows in liquid bridges. The method is based upon an appropriate extended system of perturbation equations depending on the nature of transition of the basic flow. The dependence of the critical Reynolds number and corresponding azimuthal wavenumber on serval parameters is presented for both cylindrical and non-cylindrical liquid bridges.
Jian-Kang Zhang, Bo Xun, Paul G. Chen. A Continuation Method Applied to the Study of Thermocapillary Instabilities in Liquid Bridges. Microgravity Science & Technology, 2009, 21 (Suppl 1), pp.111-117. ⟨10.1007/s12217-009-9111-2⟩. ⟨hal-01307220⟩
B Xun, Paul G. Chen, K Li, Z Yin, W.R. Hu. A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges. Advances in Space Research, Elsevier, 2008, 41 (12), pp.2094-2100. ⟨10.1016/j.asr.2007.07.016⟩. ⟨hal-01307161⟩ Plus de détails...
A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4 ⩽ Pr ⩽ 50). We focus on the relationships between the critical Reynolds number Re c , the azimuthal wave number m, the aspect ratio C and the Prandtl number Pr. A detailed Re c–Pr stability diagram is given for liquid bridges with various C. In the region of Pr > 1, which has been less studied previously and where Re c has been usually believed to decrease with the increase of Pr, we found Re c exhibits an early increase for liquid bridges with C around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocap-illary convection at large Prandtl number, and that the stability property of the basic flow is determined by the ''effective'' part of liquid bridge.
B Xun, Paul G. Chen, K Li, Z Yin, W.R. Hu. A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges. Advances in Space Research, Elsevier, 2008, 41 (12), pp.2094-2100. ⟨10.1016/j.asr.2007.07.016⟩. ⟨hal-01307161⟩