[论文解读] Deformation and orientation of a capsule with viscosity contrast in linear flows: a theoretical study

Paul Regazzi, Marc Leonetti|arXiv (Cornell University)|Feb 17, 2026

Blood properties and coagulation被引用 0

一句话总结

该论文提出一个扰动理论，用以表征在线性流场中初始为球形、黏度对比存在的囊泡（capsule）的形状与取向， Deformation 的二阶项包括表面张力与弯曲刚度，并将结果与边界积分模拟进行验证。

ABSTRACT

We develop a perturbation theory to study the shape and the orientation of an initially spherical capsule of radius R with a viscosity contrast, a surface tension σ and a bending rigidity $κ$ in linear flows. The elastic mechanical response of membrane to deformations is described by three elastic constitutive law which are either Hookean, Neohookean or Skalak type leading to the introduction of a surface shear elastic modulus $G_s$ and the Poisson ratio (or analog quantities). At the leading order, the deformation, i.e. the so-called Taylor parameter is proportional to the elastic capillary number Ca which evaluates the ratio between the external viscous stress and the elastic membrane response. In this linear regime, the results do not depend on the elastic constitutive law as expected. Without surface tension and bending rigidity, we recover the results of Barthes-Biesel & Rallison (1981) and notably the fact that the Taylor parameter does not depend on the viscosity contrast $λ$ contrary to the case of a viscous droplet. In our more general model, the deformation does no longer depend on $λ$ at the upper order. Now, the Taylor parameter also depends on two other dimensionless numbers: the surface elastocapillary ratio $σ/G_s$ and the dimensionless bending rigidity $B= κ/G_sR^2$. At the further order, the angle of inclination of the capsule with the direction of the shear flow, the analog of the Chaffey and Brenner equation for droplets is determined in each case. The results are in excellent agreement with the numerical ones performed with a code based on the boundary integral method providing an useful method to valid numerical developments.

研究动机与目标

• 理解线性流中囊泡力学的模型以类比细胞与微胶囊的行为。
• 建立一个扰动框架，在变形二阶内计算形状与取向。
• 通过弹性本构定律（胡克、Neo-Hookean、Skalak）将黏度对比、表面张力与弯曲刚性纳入分析。
• 推导囊泡的一般化取向关系（Chaffey–Brenner 型），并与数值结果比较）。

提出的方法

• 对球体周围的微小变形进行展开（r = 1 + F^(1) + F^(2)）。
• 通过 Cauchy 应力张量与表面能函数 w(I1,I2) 将膜弹性建模为胡克、Neo-Hookean 与 Skalak 本构定律。
• 通过 Helfrich 弯曲能量及表面张力项引入弯曲与表面张力，得到界面上的 f_el、f_kappa、f_sigma 力。
• 用无量纲化的 Ca = η dot{ε} R / G_s，包含黏度对比 λ、 elastocapillary Sigma = σ / G_s 与弯曲比 B = κ/(G_s R^2)。
• 用边界-积分类的谐和展开（实部谐函数）求解 Stokes 流，直到二阶，得到形变张量 F^(1)、K^(1) 与取向角。
• 针对每种本构定律给出在剪切流中的显式解析表达式，并与边界积分数值结果进行验证。

Figure 1 : Comparison between the numerical results and the theoretical equation ( 88 ) in the case of a capsule obeying to the Skalak constitutive law and with a viscosity contrast $\lambda$ : deviation from $\pi/4$ of the angle $\phi_{TT}$ of the capsule’s orientation with the direction of the she

实验结果

研究问题

• RQ1 黏度对比（λ）在线性流中超出领头项时如何影响囊泡变形与取向？
• RQ2 表面张力（σ）与弯曲刚度（κ）如何与膜的弹性本构相互作用，改变变形与取向角？
• RQ3 在不同本构定律下，剪切流中囊泡取向角的二阶修正（类似 Chaffey–Brenner 对囊泡的结果）是怎样的？
• RQ4 分析预测是否与高精度的边界积分模拟在相关参数范围（Ca、λ、Σ、B）一致？
• RQ5 该框架是否能在极限情形（如 Barthes-Biesel 与 Rallison 1981）回收已知结果并扩展到更一般的膜物理？

主要发现

• 在主导项上，变形（Taylor 参数）与 Ca 成正比，并且与内部黏度对比 λ 无关（得到 Barthes-Biesel 与 Rallison 1981 的结果）。
• 在下一个量级，变形仍与 λ 无关，而由 Σ = σ/G_s 与 B = κ/(G_s R^2) 控制。
• 剪切流中的取向角 φ_TT 给出二阶表达式（例如 Skalak 对应的 Eq. (88)），在合适极限下简化为已知结果；取向偏离 π/4，随 Ca 的变化而产生受 λ、C、B、Σ 影响的修正。
• 分析理论预测的变形分量与取向与数值边界积分模拟在所研究的参数空间内高度一致（对 λ Ca << 1 有效）。
• Neo-Hookean、Hooke 与 Skalak 本构定律给出一致的一阶变形与取向修正，并对每种情况给出明确公式（如剪切中基于 E_ij 的表达式）。
• 该框架给出平面剪切中半轴长度 L、S、W 以及截面变形 D_12、D_13、D_23 的显式表达，并确认对囊泡的一般 Chaffey–Brenner 型关系。

Figure 2 : Comparison between the numerical results and the theoretical equation ( 113 ) of the deformation of a Neo-Hookean capsule with the surface elastocapillary number $\Sigma$ under shear flow. The results are determined in the limit $\lambda\,Ca\,<<\,1$ . The range of $\Sigma$ is chosen to va