Monte Carlo simulation of liquid bridge rupture: Application to lung physiology

被引:16
作者
Alencar, Adriano M.
Wolfe, Elie
Buldyrev, Sergey V.
机构
[1] Harvard Univ, Sch Publ Hlth, Boston, MA 02115 USA
[2] Yeshiva Univ, Dept Phys, New York, NY 10033 USA
来源
PHYSICAL REVIEW E | 2006年 / 74卷 / 02期
关键词
D O I
10.1103/PhysRevE.74.026311
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
In the course of certain lung diseases, the surface properties and the amount of fluids coating the airways changes and liquid bridges may form in the small airways blocking the flow of air, impairing gas exchange. During inhalation, these liquid bridges may rupture due to mechanical instability and emit a discrete sound event called pulmonary crackle, which can be heard using a simple stethoscope. We hypothesize that this sound is a result of the acoustical release of energy that had been stored in the surface of liquid bridges prior to its rupture. We develop a lattice gas model capable of describing these phenomena. As a step toward modeling this process, we address a simpler but related problem, that of a liquid bridge between two planar surfaces. This problem has been analytically solved and we use this solution as a validation of the lattice gas model of the liquid bridge rupture. Specifically, we determine the surface free energy and critical stability conditions in a system containing a liquid bridge of volume Omega formed between two parallel planes, separated by a distance 2h, with a contact angle Theta using both Monte Carlo simulation of a lattice gas model and variational calculus based on minimization of the surface area with the volume and the contact angle constraints. In order to simulate systems with different contact angles, we vary the parameters between the constitutive elements of the lattice gas. We numerically and analytically determine the phase diagram of the system as a function of the dimensionless parameters h Omega(-1/3) and Theta. The regions of this phase diagram correspond to the mechanical stability and thermodynamical stability of the liquid bridge. We also determine the conditions for the symmetrical versus asymmetrical rupture of the bridge. We numerically and analytically compute the release of free energy during rupture. The simulation results are in agreement with the analytical solution. Furthermore, we discuss the results in connection to the rupture of similar bridges that exist in diseased lungs.
引用
收藏
页数:13
相关论文
共 36 条
[1]  
Abramowitz M., 1970, HDB MATH FUNCTIONS
[2]  
Akhiezer N. I., 1990, Elements of the Theory of Elliptic Functions
[3]   Scaling behavior in crackle sound during lung inflation [J].
Alencar, AM ;
Hantos, Z ;
Peták, F ;
Tolnai, J ;
Asztalos, T ;
Zapperi, S ;
Andrade, JS ;
Buldyrev, SV ;
Stanley, HE ;
Suki, B .
PHYSICAL REVIEW E, 1999, 60 (04) :4659-4663
[4]   Avalanche dynamics of crackle sound in the lung [J].
Alencar, AM ;
Buldyrev, SV ;
Majumdar, A ;
Stanley, HE ;
Suki, B .
PHYSICAL REVIEW LETTERS, 2001, 87 (08) :88101-1
[5]   Physiology - Dynamic instabilities in the inflating lung [J].
Alencar, AM ;
Arold, SP ;
Buldyrev, SV ;
Majumdar, A ;
Stamenovic, D ;
Stanley, HE ;
Suki, B .
NATURE, 2002, 417 (6891) :809-811
[6]  
ATHANASSENAS M, 1987, J REINE ANGEW MATH, V377, P97
[7]  
BHATT NS, 2002, PHYS REV E, V65, P6306
[8]  
BINDER K, 1997, MONTE CARLO SIMULATI
[9]   THE FORCES AND BEHAVIOR OF FLUIDS CONSTRAINED BY SOLIDS [J].
CARTER, WC .
ACTA METALLURGICA, 1988, 36 (08) :2283-2292
[10]  
FINN R, 1985, EQUILIBRIUM CAPILLAR