Splitting of a two-dimensional liquid plug at an airway bifurcation

2016 ◽  
Vol 793 ◽  
pp. 1-20 ◽  
Author(s):  
Benjamin L. Vaughan ◽  
James B. Grotberg

Certain medical treatments involve the introduction of exogenous liquids in the lungs. These liquids can form plugs within the airways. The plugs propagate throughout the branching network in the lungs being forced by airflow. They leave a deposited film on the airway walls and split at bifurcations. Understanding the resulting distribution of liquid throughout the lungs is important for the effective administration of the prescribed treatments. In this paper, we investigate numerically the splitting of a liquid plug by a two-dimensional pulmonary bifurcation under the influence of a transverse gravitational field. The splitting is characterized by the splitting ratio, which is the ratio of volume of the liquid plug in the daughter channels and depends on the capillary number and the orientation of the bifurcation plane with respect to a three-dimensional gravitational field. It is observed that gravity induces asymmetry in the splitting, causing the splitting ratio to be reduced. This effect is mitigated as the capillary number is increased. It is also observed that there exists a critical capillary number where the plug will not split and will instead propagate entirely into the gravitationally favoured daughter channel. We also compute the wall stresses on the bifurcation walls and observe the locations where stresses and their gradients are the highest in magnitude.

2005 ◽  
Vol 127 (5) ◽  
pp. 798-806 ◽  
Author(s):  
Y. Zheng ◽  
J. C. Anderson ◽  
V. Suresh ◽  
J. B. Grotberg

Many medical therapies require liquid plugs to be instilled into and delivered throughout the pulmonary airways. Improving these treatments requires a better understanding of how liquid distributes throughout these airways. In this study, gravitational and surface mechanisms determining the distribution of instilled liquids are examined experimentally using a bench-top model of a symmetrically bifurcating airway. A liquid plug was instilled into the parent tube and driven through the bifurcation by a syringe pump. The effect of gravity was adjusted by changing the roll angle (ϕ) and pitch angle (γ) of the bifurcation (ϕ=γ=0deg was isogravitational). ϕ determines the relative gravitational orientation of the two daughter tubes: when ϕ≠0deg, one daughter tube was lower (gravitationally favored) compared to the other. γ determines the component of gravity acting along the axial direction of the parent tube: when γ≠0deg, a nonzero component of gravity acts along the axial direction of the parent tube. A splitting ratio Rs, is defined as the ratio of the liquid volume in the upper daughter to the lower just after plug splitting. We measured the splitting ratio, Rs, as a function of: the parent-tube capillary number (Cap); the Bond number (Bo); ϕ; γ; and the presence of pre-existing plugs initially blocking either daughter tube. A critical capillary number (Cac) was found to exist below which no liquid entered the upper daughter (Rs=0), and above which Rs increased and leveled off with Cap. Cac increased while Rs decreased with increasing ϕ, γ, and Bo for blocked and unblocked cases at a given Cap>Cac. Compared to the nonblockage cases, Rs decreased (increased) at a given Cap while Cac increased (decreased) with an upper (lower) liquid blockage. More liquid entered the unblocked daughter with a blockage in one daughter tube, and this effect was larger with larger gravity effect. A simple theoretical model that predicts Rs and Cac is in qualitative agreement with the experiments over a wide range of parameters.


Author(s):  
Orest Shardt ◽  
J. J. Derksen ◽  
Sushanta K. Mitra

When droplets collide in a shear flow, they may coalesce or remain separate after the collision. At low Reynolds numbers, droplets coalesce when the capillary number does not exceed a critical value. We present three-dimensional simulations of droplet coalescence in a simple shear flow. We use a free-energy lattice Boltzmann method (LBM) and study the collision outcome as a function of the Reynolds and capillary numbers. We study the Reynolds number range from 0.2 to 1.4 and capillary numbers between 0.1 and 0.5. We determine the critical capillary number for the simulations (0.19) and find that it is does not depend on the Reynolds number. The simulations are compared with experiments on collisions between confined droplets in shear flow. The critical capillary number in the simulations is about a factor of 25 higher than the experimental value.


1999 ◽  
Vol 383 ◽  
pp. 29-54 ◽  
Author(s):  
ANTHONY D. SCHLEIZER ◽  
ROGER T. BONNECAZE

The dynamic behaviour and stability of a two-dimensional immiscible droplet subject to shear or pressure-driven flow between parallel plates is studied under conditions of negligible inertial and gravitational forces. The droplet is attached to the lower plate and forms two contact lines that are either fixed or mobile. The boundary-integral method is used to numerically determine the flow along and dynamics of the free surface. For surfactant-free interfaces with fixed contact lines, the deformation of the interface is determined for a range of capillary numbers, droplet to displacing fluid viscosity ratios, droplet sizes and flow type. It is shown that as the capillary number or viscosity ratio or size of the droplet increases, the deformation of the interface increases and above critical values of the capillary number no steady shape exists. For small droplets, and at low capillary numbers, shear and pressure-driven flows are shown to yield similar steady droplet shapes. The effect of surfactants is studied assuming a fixed amount of surfactant that is subject to convective–diffusive transport along the interface and no transport to or from the bulk fluids. Increasing the surface Péclet number, the ratio of convective to diffusive transport, leads to an accumulation of surfactant at the downstream end of the droplet and creates Marangoni stresses that immobilize the interface and reduce deformation. The no-slip boundary condition is then relaxed and an integral form of the Navier-slip model is used to examine the effects of allowing the droplet to slip along the solid surface in a pressure-driven flow. For contact angles less than or equal to 90°, a stable droplet spreads along the wall until a steady shape is reached, when the droplet translates across the wall at a constant velocity. The critical capillary number is larger for these droplets compared to those with pinned contact lines. For contact angles greater than 90°, the wetted area between a stable droplet and the wall decreases until a steady shape is reached. The critical capillary number for these droplets is less than that for pinned droplets. Above the critical capillary number the droplet completely detaches for a contact angle of 120°, or part of it is pinched off leaving behind a smaller attached droplet for contact angles less than or equal to 90°.


2011 ◽  
Vol 10 (2) ◽  
pp. 453-473 ◽  
Author(s):  
Jian-Jun Xu ◽  
Zhilin Li ◽  
John Lowengrub ◽  
Hongkai Zhao

AbstractIn this paper, we numerically investigate the effects of surfactant on drop-drop interactions in a 2D shear flow using a coupled level-set and immersed interface approach proposed in (Xu et al., J. Comput. Phys., 212 (2006), 590-616). We find that surfactant plays a critical and nontrivial role in drop-drop interactions. In particular, we find that the minimum distance between the drops is a non-monotone function of the surfactant coverage and Capillary number. This non-monotonic behavior, which does not occur for clean drops, is found to be due to the presence of Marangoni forces along the drop interfaces. This suggests that there are non-monotonic conditions for coalescence of surfactant-laden drops, as observed in recent experiments of Leal and co-workers. Although our study is two-dimensional, we believe that drop-drop interactions in three-dimensional flows should be qualitatively similar as the Maragoni forces in the near contact region in 3D should have a similar effect.


Author(s):  
Michael Higley ◽  
Michael Siegel ◽  
Michael R. Booty

Elastic capsules occur in nature in the form of cells and vesicles and are manufactured for biomedical applications. They are widely modelled, but there are few analytical results. In this paper, complex variable techniques are used to derive semi-analytical solutions for the steady-state response and time-dependent evolution of two-dimensional elastic capsules with an inviscid interior in Stokes flow. This provides a complete picture of the steady response of initially circular capsules in linear strain and shear flows as a function of the capillary number Ca . The analysis is complemented by spectrally accurate numerical computations of the time-dependent evolution. An imposed nonlinear strain that models the far-field velocity in Taylor's four-roller mill is found to lead to cusped steady shapes at a critical capillary number Ca c for Hookean capsules. Numerical simulation of the time-dependent evolution for Ca > Ca c shows the development of finite-time cusp singularities. The dynamics immediately prior to cusp formation are asymptotically self-similar, and the similarity exponents are predicted analytically and confirmed numerically. This is compelling evidence of finite-time singularity formation in fluid flow with elastic interfaces.


2001 ◽  
Vol 123 (6) ◽  
pp. 580-589 ◽  
Author(s):  
K. J. Cassidy ◽  
N. Gavriely ◽  
J. B. Grotberg

A finite-length liquid plug may be present in an airway due to disease, airway closure, or by direct instillation for medical therapy. Air forced by ventilation propagates the plug through the airways, where it deposits fluid onto the airway walls. The plug may encounter single or bifurcating airways, an airway surface liquid, and other liquid plugs in nearby airways. In order to understand how these flow situations influence plug transport, benchtop experiments are performed for liquid plug flow in: Case (i) straight dry tubes, Case (ii) straight pre-wetted tubes, Case (iii) bifurcating dry tubes, and Case (iv) bifurcating tubes with a liquid blockage in one daughter. Data are obtained for the trailing film thickness and plug splitting ratio as a function of capillary number and plug volumes. For Case (i), the finite length plug in a dry tube has similar behavior to a semi-infinite plug. For Case (ii), the trailing film thickness is dependent upon the plug capillary number (Ca) and not the precursor film thickness, although the shortening or lengthening of the liquid plug is influenced by the precursor film. For Case (iii), the plug splits evenly between the two daughters and the deposited film thickness depends on the local plug Ca, except for a small discrepancy that may be due to an entrance effect or from curvature of the tubes. For Case (iv), a plug passing from the parent to daughters will deliver more liquid to the unblocked daughter (nearly double, consistently) and then the plug will then travel at greater Ca in the unblocked daughter as the blocked. The flow asymmetry is enhanced for a larger blockage volume and diminished for a larger parent plug volume and parent-Ca.


1966 ◽  
Vol 25 ◽  
pp. 46-48 ◽  
Author(s):  
M. Lecar

“Dynamical mixing”, i.e. relaxation of a stellar phase space distribution through interaction with the mean gravitational field, is numerically investigated for a one-dimensional self-gravitating stellar gas. Qualitative results are presented in the form of a motion picture of the flow of phase points (representing homogeneous slabs of stars) in two-dimensional phase space.


Author(s):  
H.A. Cohen ◽  
T.W. Jeng ◽  
W. Chiu

This tutorial will discuss the methodology of low dose electron diffraction and imaging of crystalline biological objects, the problems of data interpretation for two-dimensional projected density maps of glucose embedded protein crystals, the factors to be considered in combining tilt data from three-dimensional crystals, and finally, the prospects of achieving a high resolution three-dimensional density map of a biological crystal. This methodology will be illustrated using two proteins under investigation in our laboratory, the T4 DNA helix destabilizing protein gp32*I and the crotoxin complex crystal.


Author(s):  
B. Ralph ◽  
A.R. Jones

In all fields of microscopy there is an increasing interest in the quantification of microstructure. This interest may stem from a desire to establish quality control parameters or may have a more fundamental requirement involving the derivation of parameters which partially or completely define the three dimensional nature of the microstructure. This latter categorey of study may arise from an interest in the evolution of microstructure or from a desire to generate detailed property/microstructure relationships. In the more fundamental studies some convolution of two-dimensional data into the third dimension (stereological analysis) will be necessary.In some cases the two-dimensional data may be acquired relatively easily without recourse to automatic data collection and further, it may prove possible to perform the data reduction and analysis relatively easily. In such cases the only recourse to machines may well be in establishing the statistical confidence of the resultant data. Such relatively straightforward studies tend to result from acquiring data on the whole assemblage of features making up the microstructure. In this field data mode, when parameters such as phase volume fraction, mean size etc. are sought, the main case for resorting to automation is in order to perform repetitive analyses since each analysis is relatively easily performed.


Author(s):  
Yu Liu

The image obtained in a transmission electron microscope is the two-dimensional projection of a three-dimensional (3D) object. The 3D reconstruction of the object can be calculated from a series of projections by back-projection, but this algorithm assumes that the image is linearly related to a line integral of the object function. However, there are two kinds of contrast in electron microscopy, scattering and phase contrast, of which only the latter is linear with the optical density (OD) in the micrograph. Therefore the OD can be used as a measure of the projection only for thin specimens where phase contrast dominates the image. For thick specimens, where scattering contrast predominates, an exponential absorption law holds, and a logarithm of OD must be used. However, for large thicknesses, the simple exponential law might break down due to multiple and inelastic scattering.


Sign in / Sign up

Export Citation Format

Share Document