TECHNICAL PAPERS: Fluids/Heat/Transport

Steady Propagation of a Liquid Plug in a Two-Dimensional Channel

[+] Author and Article Information
Hideki Fujioka, James B. Grotberg

Department of Biomedical Engineering, University of Michigan, Ann Arbor, Michigan 48109 USA

J Biomech Eng 126(5), 567-577 (Nov 23, 2004) (11 pages) doi:10.1115/1.1798051 History: Received March 04, 2004; Revised May 20, 2004; Online November 23, 2004
Copyright © 2004 by ASME
Your Session has timed out. Please sign back in to continue.


Long,  W. , 1991, “Effects of two Rescue Doses of a Synthetic Surfactant on Mortality Rate and Survival Without Bronchopulmonary Dysplasia in 700- to 1350-Gram Infants With Respiratory Distress Syndrome. The American Exosurf Neonatal Study Group I,” J. Pediatr. (Rio J), 118(4), (Pt 1), pp. 595–605.
Hirschl,  R. B. , 1995, “Improvement of gas Exchange, Pulmonary Function, and Lung Injury With Partial Liquid Ventilation. A Study Model in a Setting of Severe Respiratory Failure,” Chest, 108(2), pp. 500–508.
Shaffer,  T. H., and Wolfson,  M. R., 1996, “Liquid Ventilation: An Alternative Ventilation Strategy for Management of Neonatal Respiratory Distress,” Eur. J. Pediatr., 155, Suppl 2, pp. S30–34.
Baden,  H. P. , 1997, “High-Frequency Oscillatory Ventilation With Partial Liquid Ventilation in a Model of Acute Respiratory Failure,” Crit. Care Med., 25(2), pp. 299–302.
Leach,  C. L. , 1996, “Partial Liquid Ventilation With Perflubron in Premature Infants With Severe Respiratory Distress Syndrome. The LiquiVent Study Group,” N. Engl. J. Med., 335(11), pp. 761–767.
Weiss,  D. J., Bonneau,  L., and Liggitt,  D., 2001, “Use of Perfluorochemical Liquid Allows Earlier Detection of Gene Expression and use of Less Vector in Normal Lung and Enhances Gene Expression in Acutely Injured Lung,” Molecular Therapy,3(5), pp. 734–745.
Nakazawa,  K. , 2001, “Pulmonary Administration of Prostacyclin (PGI(2)) During Partial Liquid Ventilation in an Oleic Acid-Induced Lung Injury: Inhalation of Aerosol or Intratracheal Instillation?” Intensive Care Medicine,27(1), pp. 243–250.
Dickson,  E. W. , 2002, “Liquid Ventilation With Perflubron in the Treatment of Rats With Pneumococcal Pneumonia,” Crit. Care Med., 30(2), pp. 393–395.
Yu,  J., and Chien,  Y. W., 1997, “Pulmonary Drug Delivery: Physiologic and Mechanistic Aspects,” Crit. Rev. Ther. Drug Carrier Syst., 14(4), pp. 395–453.
Raczka,  E. , 1998, “The Effect of Synthetic Surfactant Exosurf on Gene Transfer in Mouse Lung in Vivo,” Gene Ther., 5(10), pp. 1333–1339.
Jobe,  A. , 1996, “Surfactant Enhances Adenovirus-Mediated Gene Expression in Rabbit Lungs,” Gene Ther., 3(9), pp. 775–779.
Jensen,  O. E., Halpern,  D., and Grotberg,  J. B., 1994, “Transport of a Passive Solute by Surfactant-Driven Flows,” Chem. Eng. Sci., 49(8), pp. 1107–1117.
Zhang,  Y. L., Matar,  O. K., and Craster,  R. V., 2003, “A Theoretical Study of Chemical Delivery Within the Lung Using Exogenous Surfactant,” Med. Eng. Phys., 25(2), pp. 115–132.
Cassidy,  K. J. , 2001, “A rat Lung Model of Instilled Liquid Transport in the Pulmonary Airways,” J. Appl. Phys., 90, pp. 1955–1967.
Hughes,  J. M. B., Rosenzweig,  D. Y., and Kivitz,  P. B., 1970, “Site of Airway Closure in Excised Dog Lungs: Histologic Demonstration,” J. Appl. Phys., 29, pp. 340–344.
Macklem,  P. T., Proctor,  D. F., and Hogg,  J. C., 1970, “The Stability of Peripheral Airways,” Respir. Physiol., 8, pp. 191–203.
Kamm,  R. D., and Schroter,  R. C., 1989, “Is Airway Closure Caused by a Thin Liquid Instability?” Respir. Physiol., 75, pp. 141–156.
Halpern,  D., and Grotberg,  J. B., 1992, “Fluid-Elastic Instabilities of Liquid-Lined Flexible Tubes,” J. Fluid Mech., 244, pp. 615–632.
Bilek,  A. M., Dee,  K. C., and Gaver,  D. P., 2003, “Mechanisms of Surface-Tension-Induced Epithelial Cell Damage in a Model of Pulmonary Airway Reopening,” J. Appl. Phys., 94(2), pp. 770–783.
Olbricht,  W. L., 1996, “Pore-Scale Prototypes of Multiphase Flow in Porous Media,” Annu. Rev. Fluid Mech., 28, pp. 187–213.
Joseph, D. D., and Renardy, Y., 1991, Fundamentals of Two-Fluid Dynamics, Vol. I: Mathematical Theory and Application; Vol. II.-Lubricated Transport, Drops, and Miscible Liquid, Springer, New York.
Cassidy,  K. J., Gavriely,  N., and Grotberg,  J. B., 2001, “Liquid Plug Flow in Straight and Bifurcating Tubes,” J. Biomech. Eng., 123(6), pp. 580–589.
Howell,  P. D., Waters,  S. L., and Grotberg,  J. B., 2000, “The Propagation of a Liquid Bolus Along a Liquid-Lined Flexible Tube,” J. Fluid Mech., 406, pp. 309–335.
Waters,  S. L., and Grotberg,  J. B., 2002, “The Propagation of a Surfactant Laden Liquid Plug in a Capillary Tube,” Phys. Fluids, 14(2), pp. 471–480.
Bretherton,  F. P., 1961, “The Motion of Long Bubbles in Tubes,” J. Fluid Mech., 10(2), pp. 166–188.
Cox,  B. G., 1962, “On Driving a Viscous Fluid out of a Tube,” J. Fluid Mech., 14(1), pp. 81–96.
Giavedoni,  M. D., and Saita,  F. A., 1997, “The Axisymmetric and Plane Cases of a Gas Phase Steadily Displacing a Newtonian Liquid—A Simultaneous Solution of the Governing Equations,” Phys. Fluids, 9(8), pp. 2420–2428.
Heil,  M., 2001, “Finite Reynolds Number Effects in the Bretherton Problem,” Phys. Fluids, 13(9), pp. 2517–2521.
Wassmuth,  F., Laidlaw,  W. G., and Coombe,  D. A., 1993, “Calculation of Interfacial Flows and Surfactant Redistribution as a Gas-Liquid Interface Moves Between 2 Parallel Plates,” Phys. Fluids A, 5(7), pp. 1533–1548.
Giavedoni,  M. D., and Saita,  F. A., 1999, “The Rear Meniscus of a Long Bubble Steadily Displacing a Newtonian Liquid in a Capillary Tube,” Phys. Fluids, 11(4), pp. 786–794.
Horsfield,  K., Dart,  G., and Olsen,  D. E., 1971, “Models of the Human Bronchial Tree,” J. Appl. Phys., 11(2), pp. 207–217.
Patankar, S. V., 1980, Numerical Heat Transfer and Fluid Flow, Hemisphere, London.
Thompson, J. F., Soni, B. K., and Weatherill, N. P., 1999, Handbook of Grid Generation, CRC Press, Boca Raton, Fla.
Wesseling,  P. , 1998, “Computing Flows on General Two-Dimensional Nonsmooth Staggered Grids,” J. Eng. Math., 34(1–2), pp. 21–44.
Thompson,  J. F., and Warsi,  Z. U., 1982, “Boundary-Fitted Coordinate System for Numerical Solution of Partial Differential Equations—A Review,” J. Comput. Phys., 47, pp. 1–108.
Deen, W. M., 1998, in Analysis of Transport Phenomena, edited by K. E. Gubbins, Oxford University Press, New York.
Muzaferija,  S., and Peric,  M., 1997, “Computation of Free-Surface Flows Using the Finite-Volume Method and Moving Grids,” Numer. Heat Transfer, Part B, 32(4), pp. 369–384.
Press, W. H. et al., 1992, Numer. Recipes in C, 2nd ed., Cambridge University Press, Cambridge, England.
Demmel,  J. W. , 1999, “A Supernodal Approach to Sparse Partial Pivoting,” SIAM J. Matrix Anal. Appl., 20(3), pp. 720–755.
Halpern,  D., and Gaver,  D. P., 1994, “Boundary-Element Analysis of the Time-Dependent Motion of a Semiinfinite Bubble in a Channel,” J. Comput. Phys., 115(2), pp. 366–375.
Heil,  M., 2000, “Finite Reynolds Number Effects in the Propagation of an Air Finger Into a Liquid-Filled Flexible-Walled Channel,” J. Fluid Mech., 424, pp. 21–44.


Grahic Jump Location
Comparison with other studies for the trailing film thickness h1 as a function of Ca at Re=0. The square points represent the results of the present study for LP=2. The broken line 27, the dashed line 40, and the solid line 25 represents the trailing film thickness of leading meniscus for a semi-infinite bubble, i.e., LP=∞.
Grahic Jump Location
The streamlines and pressure fields inside the liquid plug at Re=80 and λ=1000, for (a) LP=0.25, (b) LP=1, and (c) LP=2
Grahic Jump Location
The dimensionless macroscopic pressure gradient ΔP/LP versus Re at λ=1000 for different LP
Grahic Jump Location
The trailing film thickness, h1 versus Ca for different LP at Re=0
Grahic Jump Location
Streamlines and pressure fields inside the liquid plug at Ca=0.05 and Re=0, for (a) LP=0.25, (b) LP=1, and (c) LP=2. The directed lines represent the flow field, and the dashed lines denote lines of constant pressure. S1,S2,S3, and S4 are the locations of the stagnation points.
Grahic Jump Location
Streamlines and pressure fields inside the liquid plug at Ca=0.4 and Re=0 for (a) LP=0.25, (b) LP=1, and (c) LP=2      
Grahic Jump Location
The effect of inertia (Re) on the trailing film thickness h1 at λ=1000 for different LP. The point at Re=50 represents h1 for a semi-infinite bubble propagation 27.
Grahic Jump Location
Axial velocity profiles at the middle cross-section of the liquid plug (x=0) at Re=50, λ=1000 for different values of LP. The curve of “parabolic” is u=−1+32(1−h1)(1−y2) with h1=h1 (LP=2).
Grahic Jump Location
The streamlines and pressure fields inside the liquid plug at Re=50 and λ=1000, for (a) LP=0.25, (b) LP=1, and (c) LP=2
Grahic Jump Location
A sketch of the liquid plug model. A pressure difference between the front and back air finger, ΔP*=P1*−P2*, drives the liquid plug of the length LP* with constant speed U within a two-dimensional channel of half width H lined by a precursor film of thickness h2*. For steady state h2* is equal to the trailing film thickness at the rear end boundary h1*.
Grahic Jump Location
An example of a two-dimensional grid generated in the domain. A staggered grid is used where the pressure p is stored at the center of each cell and the velocity components are stored at control surfaces.
Grahic Jump Location
The wall shear stress, τ distribution, versus x for LP=2, λ=1000, and Re=40, 60, and 80
Grahic Jump Location
Velocity vectors distribution in the capillary wave for LP=2, Re=50, and λ=1000



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In