Enzyme Catalysis Causes Fluid Flow, Motility, and Directional Transport on Supported Lipid Bilayers

The dynamic interplay between the composition of lipid membranes and the behavior of membrane-bound enzymes is critical to the understanding of cellular function and viability, and the design of membrane-based biosensing platforms. While there is a significant body of knowledge about how lipid composition and dynamics affect membrane-bound enzymes, little is known about how enzyme catalysis influences the motility and lateral transport on lipid membranes. Using enzyme-attached lipids in supported bilayers (SLBs), we provide direct evidence of catalysis-induced fluid flow that underlies the observed motility on SLBs. Additionally, by using active enzyme patches, we demonstrate the directional transport of tracer particles on SLBs. As expected, enhancing the membrane viscosity by incorporating cholesterol into the bilayer suppresses the overall movement. These are the first steps in understanding diffusion and transport on lipid membranes due to active, out-of-equilibrium processes that are the hallmark of living systems. In general, our study demonstrates how active enzymes can be used to control diffusion and transport in confined 2-D environments.

Critical conditions for development of a second pair of Dean vortices in curved microfluidic channels

The centrifugal force in flow through a curved channel initiates a hydrodynamic instability that results in the development of Dean vortices, a pair of counter-rotating roll cells across the channel that deflect the high velocity fluid in the center toward the outer (concave) wall. When this secondary flow toward the concave (outer) wall is too strong to be dissipated by viscous effects, an additional pair of vortices emerges near the outer wall. Combining numerical simulation and dimensional analysis, we find that the critical condition for the onset of the second vortex pair depends on γ^{1/2}Dn (γ: channel aspect ratio; Dn: Dean number). We also investigate the development length for the additional vortex pair in channels with different aspect ratios and curvatures. The larger centrifugal force at higher Dean numbers produces the additional vortices further upstream, with the required development length being inversely proportional to the Reynolds number and increasing linearly with the radius of curvature of the channel.

Simultaneous Screening of Nitrofuran Metabolites in Honey Using Biochip Array Technology: Validation Study According to the Decision 2002/657/EC of the European Union

Background

Although no authorization is available for antibiotics to treat bee diseases, some veterinary compounds are used by beekeepers, and each country sets its own thresholds. Inappropriate and excessive use of these drugs can cause allergic reactions and antibiotic resistance in humans who consume the remaining antibiotic residues in honey and its products. It is, therefore, relevant to monitor the presence of antibiotic residues in this matrix.

Objectives

A rapid method for the simultaneous screening of nitrofuran metabolite residues in honey was validated according to Commission Decision 2002/657/EC (C.D 657) and the European guideline for the validation of screening methods for veterinary medicines.

Methods

This multi-analytical screening method enables the simultaneous determination of four nitrofuran metabolites [3-amino-2-oxazolidone (AOZ), 3-amino-5-morpholinomethyl-2-oxazolidinone (AMOZ), 1-Aminohydantoin HCl (ADH), and semicarbazide (SEC)] from a single honey sample. Thirty-five honey samples were collected randomly as real samples for screening from Tehran, IR Iran, Germany, and the Netherlands in 2018.

Results

For all four antibiotic residues, the positivity threshold T was higher than the cut-off value Fm, and no false-positive results were obtained for three antibiotics (AOZ, AMOZ, and SEC). Detection capabilities (CCβ) of all compounds were under the minimum required performance limit (MRPL) authorized by the European Commission (currently 1 μg/kg). The screening results of 15 domestic and 20 imported honey samples showed that the levels of AOZ in 6.66% and 10% of the samples, the level of AMOZ in 13.33% and 0% of the samples, and the level of SEC in 33.33% and 40% of the samples were less than the cut-off ([in relative light units (RLUs)], respectively.

Conclusions

This study found that this technique is valid for detecting and quantifying three antibiotic residues in honey samples at the measured validation levels. This method was simple, rapid, and capable of simultaneously screening three nitrofuran metabolites from a single honey sample.

Photothermal Atomic Force Microscopy Coupled with Infrared Spectroscopy (AFM-IR) Analysis of High Extinction Coefficient Materials: A Case Study with Silica and Silicate Glasses

Photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) brings significant value as a spatially resolved surface analysis technique for disordered oxide materials such as glasses, but additional development and fundamental understanding of governing principles is needed to interpret AFM-IR spectra, since the existing theory described for organic materials does not work for materials with high extinction coefficients for infrared (IR) absorption. This paper describes theoretical calculation of a transient temperature profile inside the IR-absorbing material considering IR refraction at the interface as well as IR adsorption and heat transfer inside the sample. This calculation explains the differences in peak positions and amplitudes of AFM-IR spectra from those of specular reflectance and extinction coefficient spectra. It also addresses the information depth of the AFM-IR characterization of bulk materials. AFM-IR applied to silica and silicate glass surfaces has demonstrated novel capability of characterizing subsurface structural changes and surface heterogeneity due to mechanical stresses from physical contacts, as well as chemical alterations manifested in surface layers through aqueous corrosion.

Scale-up Issues for Commercial Depth Filters in Bioprocessing

Significant increases in cell density and product titer have led to renewed interest in the application of depth filtration for initial clarification of cell culture fluid in antibody production. The performance of these depth filters will depend on the local pressure and velocity distribution within the filter capsule, but these are very difficult to probe experimentally, leading to challenges in both process design and scale-up. We have used a combination of carefully designed experimental studies and computational fluid dynamics (CFD) to examine these issues in both lab-scale (Supracap^TM 50) and pilot-scale (Stax^TM ) depth filter modules, both employing dual-layer lenticular PDH4 media containing diatomaceous earth. The Supracap^TM 50 showed a more rapid increase in transmembrane pressure and a more rapid DNA breakthrough during filtration of a Chinese Hamster Ovary cell culture fluid. These results were explained using CFD calculations which showed very different flow distributions within the modules. CFD predictions were further validated using measurements of the residence time distribution and dye binding in both the lab-scale and pilot-plant modules. These results provide important insights into the factors controlling the performance and scale-up of these commercially important depth filters as well as a framework that can be broadly applied to develop more effective depth filters and depth filtration processes. This article is protected by copyright. All rights reserved.

A Volume-Corrected Wenzel Model

The Wenzel model, commonly used for predicting the equilibrium contact angle (CA) of drops which penetrate the asperities of a rough surface, does not account for the liquid volume stored in the asperities. Interestingly, many previous experimental and molecular dynamics studies have noted discrepancies between observed CAs and those predicted by the Wenzel model because of this neglected liquid volume. Here, we apply a thermodynamic model to wetting of periodically patterned surfaces to derive a volume-corrected Wenzel equation in the limit of small pattern wavelength (compared to drop size). We show that the corrected equilibrium CA is smaller than that predicted by the Wenzel equation and that the reduction in CA can be significant when the liquid volume within the asperities becomes non-negligible compared to the total droplet volume. In such cases, the corrected CAs agree reasonably well with experimental observations and results of molecular dynamics simulations reported in previous studies.

Conducting polymer microcontainers for biomedical applications

Advancement in the development of metallic-based implantable micro-scale bioelectronics has been limited by low signal to noise ratios and low charge injection at electrode-tissue interfaces. Further, implantable electrodes lose their long-term functionality because of unfavorable reactive tissue responses. Thus, substantial incentive exists to produce bioelectronics capable of delivering therapeutic compounds while improving electrical performance. Here, we have produced hollow poly(pyrrole) microcontainers (MCs) using poly(lactic-co-glycolic) acid (PLGA) as degradable templates. We demonstrate that the effective surface area of the electrode increases significantly as deposition charge density is increased, resulting in a 91% decrease in impedance and an 85% increase in charge storage capacity versus uncoated gold electrodes. We also developed an equivalent circuit model to quantify the effect of conducting polymer film growth on impedance. These MC-modified electrodes offer the potential to improve the electrical properties of implantable bioelectronics, as well as provide potential controlled release avenues for drug delivery applications.

Conducting Polymer Microcups for Organic Bioelectronics and Drug Delivery Applications

An ideal neural device enables long-term, sensitive, and selective communication with the nervous system. To accomplish this task, the material interface should mimic the biophysical and the biochemical properties of neural tissue. By contrast, microfabricated neural probes utilize hard metallic conductors, which hinder their long-term performance because these materials are not intrinsically similar to soft neural tissue. This study reports a method for the fabrication of monodisperse conducting polymer microcups. It is demonstrated that the physical surface properties of conducting polymer microcups can be precisely modulated to control electrical properties and drug-loading/release characteristics.

Effect of Gravity on the Configuration of Droplets on Two-Dimensional Physically Patterned Surfaces

Wetting of solid surfaces is important for many potential applications, including the design of low-drag and antifouling/self-cleaning surfaces, and it is usually quantified by the contact angle and by contact angle hysteresis. Both the chemistry and the physical patterning of the surface are known to affect the contact angle. In studying the wetting of such surfaces, most models focus on the small Bond number (Bo) limit in which the effect of gravity is negligible, which simplifies free energy calculations. In this work, we employ a thermodynamic model for surfaces patterned with two-dimensional asperities, which remains applicable for nonzero Bo. We employ two versions of the model: one in which we require the liquid-vapor interface to remain a circular cap, and another in which we allow the liquid-vapor interface to deform. We find that the effects of gravity are twofold. First, drops with larger Bo tend to flatten and spread across the surface relative to the same size drops with Bo = 0. Second, gravity makes it more favorable for drops to penetrate surface asperities compared to the case of Bo = 0, which also tends to lower the contact angles. The main effect of droplet deformation is to produce larger contact angles for the same wetting configuration. Finally, we compare our model predictions with relevant experimental observations. We find very close agreement with the experiments, thereby validating our theoretical model.

Effects of Hierarchical Surface Roughness on Droplet Contact Angle

Superhydrophobic surfaces often incorporate roughness on both micron and nanometer length scales, although a satisfactory understanding of the role of this hierarchical roughness in causing superhydrophobicity remains elusive. We present a two-dimensional thermodynamic model to describe wetting on hierarchically grooved surfaces by droplets for which the influence of gravity is negligible. By creating wetting phase diagrams for droplets on surfaces with both single-scale and hierarchical roughness, we find that hierarchical roughness leads to greatly expanded superhydrophobic domains in phase space over those for a single scale of roughness. Our results indicate that an important role of the nanoscale roughness is to increase the effective Young's angle of the microscale features, leading to smaller required aspect ratios (height to width) for the surface structures. We then show how this idea may be used to design a hierarchically rough surface with optimally high contact angles.

Kinetics of droplet wetting mode transitions on grooved surfaces: forward flux sampling

The wetting configuration of a liquid droplet on a rough or physically patterned surface is typically characterized by either the Cassie wetting mode, in which the droplet resides on top of the roughness, or the Wenzel mode, in which the droplet penetrates into the roughness. For a fixed surface topology and droplet size, one of these modes corresponds to the global free-energy minimum. However, the other state is often metastable and long-lived due to a free-energy barrier that hinders the transition between the two wetting states. Metastable wetting states have been observed experimentally, and we also observe them in molecular dynamics (MD) simulations of a droplet on a grooved surface. Using forward flux sampling, we study the kinetics of the Cassie to Wenzel and Wenzel to Cassie transitions for two-dimensional droplets on periodically grooved substrates. The global-minimum wetting states that emerge from our nanoscale MD approach are consistent with those predicted by a macroscopic model based on free energy minimization. We find that the free-energy barriers for these transitions depend on the droplet size and surface topology. A committor analysis indicates that the transition-state ensemble consists of droplets that are on the verge of initiating/breaking contact with the substrate at the bottom of the grooves.

Locomotion of microorganisms near a no-slip boundary in a viscoelastic fluid

Locomotion of microorganisms plays a vital role in most of their biological processes. In many of these processes, microorganisms are exposed to complex fluids while swimming in confined domains, such as spermatozoa in mucus of mammalian reproduction tracts or bacteria in extracellular polymeric matrices during biofilm formation. Thus, it is important to understand the kinematics of propulsion in a viscoelastic fluid near a no-slip boundary. We use a squirmer model with a time-reversible body motion to analytically investigate the swimming kinematics in an Oldroyd-B fluid near a wall. Analysis of the time-averaged motion of the swimmer shows that both pullers and pushers in a viscoelastic fluid swim towards the no-slip boundary if they are initially located within a small domain of "attraction" in the vicinity of the wall. In contrast, neutral swimmers always move towards the wall regardless of their initial distance from the wall. Outside the domain of attraction, pullers and pushers are both repelled from the no-slip boundary. Time-averaged locomotion is most pronounced at a Deborah number of unity. We examine the swimming trajectories of different types of swimmers as a function of their initial orientation and distance from the no-slip boundary.

Kinematic matrix theory and universalities in self-propellers and active swimmers

We describe an efficient and parsimonious matrix-based theory for studying the ensemble behavior of self-propellers and active swimmers, such as nanomotors or motile bacteria, that are typically studied by differential-equation-based Langevin or Fokker-Planck formalisms. The kinematic effects for elementary processes of motion are incorporated into a matrix, called the "kinematrix," from which we immediately obtain correlators and the mean and variance of angular and position variables (and thus effective diffusivity) by simple matrix algebra. The kinematrix formalism enables us recast the behaviors of a diverse range of self-propellers into a unified form, revealing universalities in their ensemble behavior in terms of new emergent time scales. Active fluctuations and hydrodynamic interactions can be expressed as an additive composition of separate self-propellers.

Nanomotor mechanisms and motive force distributions from nanorotor trajectories

Nanomotors convert chemical energy into mechanical motion. For a given motor type, the underlying chemical reaction that enables motility is typically well known, but the detailed, quantitative mechanism by which this reaction breaks symmetry and converts chemical energy to mechanical motion is often less clear, since it is difficult experimentally to measure important parameters such as the spatial distribution of chemical species around the nanorotor during operation. Without this information on how motor geometry affects motor function, it is difficult to control and optimize nanomotor behavior. Here we demonstrate how one easily observable characteristic of nanomotor operation-the visible trajectory of a nanorotor-can provide quantitative information about the role of asymmetry in nanomotor operation, as well as insights into the spatial distribution of motive force along the surface of the nanomotor, the motive torques, and the effective diffusional motion.

Wetting on physically patterned solid surfaces: the relevance of molecular dynamics simulations to macroscopic systems

We used molecular dynamics (MD) simulations to study the wetting of Lennard-Jones cylindrical droplets on surfaces patterned with grooves. By scaling the surface topography parameters with the droplet size, we find that the preferred wetting modes and contact angles become independent of the droplet size. This result is in agreement with a mathematical model for the droplet free energy at small Bond numbers for which the effects of gravity are negligible. The MD contact angles for various wetting modes are in good agreement with those predicted by the mathematical model. We construct phase diagrams of the dependence of the wetting modes observed in the MD simulations on the topography of the surface. Depending on the topographical parameters characterizing the surface, multiple wetting modes can be observed, as is also seen experimentally. Thus, our studies indicate that MD simulations can yield insight into the large-length-scale behavior of droplets on patterned surfaces.

Microencapsulation of chemotherapeutics into monodisperse and tunable biodegradable polymers via electrified liquid jets: control of size, shape, and drug release

This paper describes microencapsulation of antitumor agent 1,3-bis(2-chloroethyl)-1-nitrosourea (BCNU, Carmustine) into biodegradable polymer poly(lactic-co-glycolic) acid (PLGA) using an electrojetting technique. The resulting BCNU-loaded PLGA microcapsules have significantly higher drug encapsulation efficiency, more tunable drug loading capacity, and (3) narrower size distribution than those generated using other encapsulation methods.

Chiral diffusion of rotary nanomotors

Neither a purely deterministic rotary nanomotor nor a purely orientational diffuser exhibits long-term translational motion, but coupling rotation to orientational diffusion yields translational diffusion. We demonstrate that this effective translational diffusion can easily dominate the ordinary thermal translational diffusion for experimentally relevant nanomotors, and that this effective diffusion is chiral. Unpowered chiral particles do not exhibit chiral diffusion, but a nanorotor has both handedness and an instantaneous direction of powered motion, thus-unlike an unpowered particle-its diffusional motion can distinguish left from right.

A theory for the morphological dependence of wetting on a physically patterned solid surface

We present a theoretical model for predicting equilibrium wetting configurations of two-dimensional droplets on periodically grooved hydrophobic surfaces. The main advantage of our model is that it accounts for pinning/depinning of the contact line at step edges, a feature that is not captured by the Cassie and Wenzel models. We also account for the effects of gravity (via the Bond number) on various wetting configurations that can occur. Using free-energy minimization, we construct phase diagrams depicting the dependence of the wetting modes (including the number of surface grooves involved in the wetting configuration) and their corresponding contact angles on the geometrical parameters characterizing the patterned surface. In the limit of vanishing Bond number, the predicted wetting modes and contact angles become independent of drop size if the geometrical parameters are scaled with drop radius. Contact angles predicted by our continuum-level theoretical model are in good agreement with corresponding results from nanometer-scale molecular dynamics simulations. Our theoretical predictions are also in good agreement with experimentally measured contact angles of small drops, for which gravitational effects on interface deformation are negligible. We show that contact-line pinning is important for superhydrophobicity and that the contact angle is maximized when the droplet size is comparable to the length scale of the surface pattern.

Ozone Uptake During Inspiratory Flow in a Model of the Larynx, Trachea and Primary Bronchial Bifurcation

Three-dimensional simulations of the transport and uptake of a reactive gas such as O(3) were compared between an idealized model of the larynx, trachea, and first bifurcation and a second "control" model in which the larynx was replaced by an equivalent, cylindrical, tube segment. The Navier-Stokes equations, Spalart-Allmaras turbulence equation, and convection-diffusion equation were implemented at conditions reflecting inhalation into an adult human lung. Simulation results were used to analyze axial velocity, turbulent viscosity, local fractional uptake, and regional uptake. Axial velocity data revealed a strong laryngeal jet with a reattachment point in the proximal trachea. Turbulent viscosity data indicated that jet turbulence occurred only at high Reynolds numbers and was attenuated by the first bifurcation. Local fractional uptake data affirmed hotspots previously reported at the first carina, and suggested additional hotspots at the glottal constriction and jet reattachment point in the proximal trachea. These laryngeal effects strongly depended on inlet Reynolds number, with maximal effects (approaching 15%) occurring at maximal inlet flow rates. While the increase in the regional uptake caused by the larynx subsided by the end of the model, the effect of the larynx on cumulative uptake persisted further downstream. These results suggest that with prolonged exposure to a reactive gas, entire regions of the larynx and proximal trachea could show signs of tissue injury.

Gas Uptake in a Three-Generation Model Geometry During Steady Expiration: Comparison of Axisymmetric and Three-Dimensional Models

Mass transfer coefficients were predicted and compared for uptake of a formaldehyde-air gas system using an axisymmetric single path model (ASPM) and a three-dimensional computational fluid dynamics model (CFDM) in three-generation model geometry at steady expiratory flow. The flow and concentration fields in the ASPM were solved using Galerkin's finite-element method and in the CFDM using a commercial finite-element software, FIDAP. Numerical results were compared for two different inlet flow rates, wall mass transfer coefficients, and bifurcation angles. The mass transfer coefficients variation with bifurcation unit from the ASPM and CFDM compared qualitatively and quantitatively closely at all flows and lower wall mass transfer coefficients for both 40 degrees and 70 degrees bifurcation angles. However, at higher wall mass transfer coefficients, quantitatively they were within 40% for both the bifurcation angles. Also, at higher flow and wall mass transfer coefficients, they were off qualitatively for a 70 degrees bifurcation angle although the uptake compared qualitatively. This is due to the normalization of uptake within a bifurcation unit with the average of inlet and outlet average concentrations. Both CFDM and ASPM predict the same trends of increase in mass transfer coefficients with inlet flow and wall mass transfer coefficients. Also, the local values of the mass transfer coefficients compared closely at all conditions. These results validate the simplified ASPM and the complex CFDM. Mass transfer coefficients increase with bifurcation angles and with a flat inlet velocity profile compared to a parabolic velocity profile since the flow is non-fully developed and hence, the uptake increases.

Gas Uptake in a Three-Generation Model Geometry with a Flat Inlet Velocity During Steady Inspiration: Comparison of Axisymmetric and Three-Dimensional Models

Mass transfer coefficients were predicted and compared for uptake of reactive gas system using an axisymmetric single-path model (ASPM) with experimentally predicted values in a two-generation geometry and with a three-dimensional computational fluid dynamics model (CFDM) in a three-generation model geometry at steady inspiratory flow with a flat inlet velocity profile. The flow and concentration fields in the ASPM were solved using Galerkin's finite element method and in the CFDM using a commercial finite element software FIDAP. ASPM predicted average gas phase mass transfer coefficients within 25% of the experimental values. Numerical results in terms of overall mass transfer coefficients from the two models within each bifurcation unit were compared for two different inlet flow rates, wall mass transfer coefficients, and bifurcation angles. The overall mass transfer coefficients variation with bifurcation unit from the ASPM and CFDM compared qualitatively and quantitatively closely at lower wall mass transfer coefficients for both 40 degree and 70 degree bifurcation angles. But at higher wall mass transfer coefficients, quantitatively they were off in the range of 2-10% for 40 degree bifurcation angle and in the range of 4-15% for 70 degree bifurcation angle. Both CFDM and ASPM predict the same trends of increase in mass transfer coefficients with inlet flow, wall mass transfer coefficients, and during inspiration compared to expiration. Higher mass transfer coefficients were obtained with a flat velocity profile compared to a parabolic velocity profile using ASPM. These results validate the simplified ASPM and the complex CFDM.

Comparison of Axisymmetric and Three-Dimensional Models for Gas Uptake in a Single Bifurcation During Steady Expiration

Reactive gas uptake is predicted and compared in a single bifurcation at steady expiratory flow in terms of Sherwood number using an axisymmetric single-path model (ASPM) and a three-dimensional computational fluid dynamics model (CFDM). ASPM is validated in a two-generation geometry by comparing the average gas-phase mass transfer coefficients with the experimental values. ASPM predicted mass transfer coefficients within 20% of the experimental values. The flow and concentration variables in the ASPM were solved using Galerkin finite element method and in the CFDM using commercial finite element software FIDAP. The simulations were performed for reactive gas flowing at Reynolds numbers ranging from 60 to 350 in both symmetric bifurcation for three bifurcation angles, 30deg, 70deg, and 90deg, and in an asymmetric bifurcation. The numerical models compared with each other qualitatively but quantitatively they were within 0.4–8% due to nonfully developed flow in the parent branch predicted by the CFDM. The radially averaged concentration variation along the axial location matched qualitatively between the CFDM and ASPM but quantitatively they were within 32% due to differences in the flow field. ASPM predictions compared well with the CFDM predictions for an asymmetric bifurcation. These results validate the simplified ASPM and the complex CFDM. ASPM predicts higher Sherwood number with a flat velocity inlet profile compared to a parabolic inlet velocity profile. Sherwood number increases with the inlet average velocity, wall mass transfer coefficient, and bifurcation angle since the boundary layer grows slower in the parent and daughter branches.

Thermocapillary convection in double-layer fluid structures within a two-dimensional open cavity

Thermocapillary convection within a differentially-heated open rectangular cavity containing two immiscible liquid layers is considered in the absence of gravitational effects. The temperature and flow fields in the two layers are computed using domain mapping in conjunction with a finite-difference scheme on a staggered grid. The melt-encapsulant and air-encapsulant interfaces are allowed to deform, with the contact lines pinned on the solid boundaries. The presence of a free surface at the top leads to increased convection in the encapsulant phase while retarding thermocapillary flow in the melt. The intensity of thermocapillary convection in the encapsulated layer is reduced as the viscosity of the encapsulant is increased or the thickness of the encapsulant layer is decreased. Choosing an encapsulant with a greater sensitivity of interfacial tension to temperature (as compared to that of the melt phase) can almost completely suppress thermocapillary convection in the melt. Deformations of the melt-encapsulant interface in an open cavity are found to be larger than those in a closed cavity with a rigid top surface, due to higher pressure gradients realized in the encapsulant phase. In contrast to interface deformation behavior reported earlier for a double-layer system in a closed cavity, the shape of the melt-encapsulant interface is qualitatively similar for all values of the viscosity ratio, with the interface dipping into the melt near the cold wall, and into the encapsulant near the hot wall. For the double-layers considered in this study, a free surface at the top of the encapsulant layer was found to be more effective than a rigid top in reducing the intensity of thermocapillary convection in the melt.

Three-Dimensional Simulations of Reactive Gas Uptake in Single Airway Bifurcations

The pattern of lung injury induced by the inhalation of ozone (O(3)) depends on the dose delivered to different tissues in the airways. This study examined the distribution of O(3) uptake in a single, symmetrically branched airway bifurcation. Reaction in the epithelial lining fluid was assumed to be so rapid that O(3) concentration was negligible along the entire surface of the bifurcation wall. Three-dimensional numerical solutions of the continuity, Navier-Stokes and convection-diffusion equations were obtained for steady inspiratory and expiratory flows at Reynolds numbers ranging from 100 to 500. The total rate of O(3) uptake was found to increase with increasing flow rate during both inspiration and expiration. Hot spots of O(3) flux appeared at the carina of the bifurcation for virtually all inspiratory and expiratory Reynolds numbers considered in the simulations. At the lowest expiratory Reynolds number, however, the location of the maximum flux was shifted to the outer wall of the daughter branch. For expiratory flow, additional hot spots of flux were found on the parent branch wall just downstream of the branching region. In all cases, O(3) uptake in the single bifurcation was larger than that in a straight tube of equal inlet radius and wall surface area. This study provides insight into the effect of flow conditions on O(3) uptake and dose distribution in individual bifurcations.

Numerical simulation of the spontaneous penetration of liquids into cylindrical capillaries

Dynamics of spontaneous capillary penetration of a liquid into a cylindrical pore is studied numerically over the entire duration of an experiment, including the initial stages of penetration during which inertial effects are dominant. Partial slip in the vicinity of the moving contact-line is allowed by using an empirical constitutive relation between the dynamic contact angle and contact line speed in order to avoid the stress singularity arising from the presence of the moving contact line on the solid wall. A finite-difference scheme on a staggered body-fitted grid is used to solve for the time-dependent flow field and to determine the time evolution of the shape of the advancing meniscus. The results of dynamic simulations of capillary rise under both normal and microgravity conditions are compared with the reported experimental observations. The simulation results are found to be in good agreement with the experimental measurements in both cases. Numerical simulations capture the different flow regimes identified in previous studies of spontaneous capillary penetration. The velocity dependence of the dynamic contact line is found to have a significant effect on kinetics of wetting in the intermediate-time flow regime, wherein the capillary force is balanced by convective losses.

Double-layer thermocapillary convection in a differentially heated cavity

Many materials-processing applications such as crystal growth from the melt involve thermocapillary flows that can affect the quality of the final product, particularly under microgravity conditions where the influence of buoyancy-driven convection is minimized. When the melt contains volatile components, as in the production of III-V semiconductor crystals, it is often encapsulated in a low-melting point amorphous molten glass phase such as boron oxide or pyrolytic boron nitride in order to prevent evaporation of the volatile components. The addition of the encapsulant layer and the melt-encapsulant interface in such cases can alter the thermocapillary flow in the melt. In this study, thermocapillary convection within a differentially heated rectangular cavity containing two immiscible liquid layers is considered in the absence of gravity. Domain mapping is used in conjunction with a finite difference scheme on a staggered grid to solve for the temperature and flow fields. The melt-encapsulant and the air-encapsulant interfaces are allowed to deform, with the contact lines pinned on the solid boundaries. The computed flow fields are compared to the corresponding results for a cavity with a rigid top surface. The presence of a free surface at the top leads to increased convection in the encapsulant phase while suppressing the thermocapillary flow in the melt phase. The flow pattern in the encapsulated layer is strongly dependent on the viscosity of the encapsulant layer. The intensity of the thermocapillary flow within the melt is significantly reduced as the viscosity of the encapsulant layer is increased. However, for a higher encapsulant viscosity, the retarding effect of the free top surface on thermocapillary convection in the melt is weakened.

Spontaneous spreading of surfactant-bearing drops in the sorption-controlled limit

Axisymmetric spreading of a liquid drop containing a soluble surfactant on a smooth solid substrate is numerically investigated for the case in which surfactant mass transfer between the interface and the bulk liquid is sorption/kinetic controlled. The fastest spreading rate is achieved by drops with O(1) values of Biot number for which the rate of surface convection is comparable to the sorption rate, and the surfactant molecules transferred to the interface are effectively convected to the contact line region.

Coalescence of Bubbles Translating through a Tube

The results of an experimental study of the interaction and coalescence of two air bubbles translating in a cylindrical tube are presented. Both pressure- and buoyancy-driven motion of the two bubbles in a Newtonian suspending fluid within the tube are considered. The close approach of the two bubbles is examined using image analysis, and measurements of the coalescence time are reported for various bubble size ratios and capillary numbers. For pressure-driven motion of bubbles, coalescence is found to occur in an axisymmetric configuration for all bubble size ratios considered in the experiments. For buoyancy-driven motion, on the other hand, the disturbance flow behind the leading bubble causes the trailing bubble to move radially out toward the tube wall when the trailing bubble size becomes very small compared to the size of the leading bubble. In that case, coalescence occurs in a nonaxisymmetric configuration, with a time scale for coalescence that is substantially larger than that for coalescence in the axisymmetric configuration. When the imposed flow is in the direction of the buoyancy force, coalescence time is independent of bubble size ratio, and decreases as the capillary number increases. Experimental measurements of the radius of the thin liquid film separating the two bubbles are used in conjunction with a simple film drainage model to predict the dependence of the coalescence time on the bubble size ratio.

Pressure-driven motion of surfactant-laden drops through cylindrical capillaries: effect of surfactant solubility

The effect of bulk-soluble surfactants on the dynamics of a drop translating through a cylindrical tube under low-Reynolds-number conditions is investigated. Interfacial surfactant adsorption/desorption is modeled according to the Frumkin adsorption framework, and the bulk-insoluble surfactant limit is recovered as the rate of surfactant sorption becomes large compared to that of bulk diffusion. As the equilibrium surface coverage is increased, the mechanism by which drop mobility is reduced changes from uniform retardation at low surface coverage to the formation of a stagnant cap at high surface coverage. For large capillary numbers, the drop does not achieve a steady shape, and eventually it breaks up either through the formation of a penetrating viscous jet of suspending fluid, or by continuous elongation and pinch-off. Surfactants have a destabilizing effect on transient drop shapes by accelerating the formation and development of the penetrating viscous jet that leads to drop breakup. The critical conditions for drop breakup, as well as the mode of breakup, depend on the manner in which the strength of the flow (i.e., the capillary number) is increased.

Thermocapillary flow in double-layer fluid structures: An effective single-layer model

Thermocapillary flows are of considerable technological importance in materials processing applications such as crystal growth from the melt, particularly under microgravity conditions where the influence of buoyancy-driven convection is minimized. In this study, thermally driven convection within a differentially heated rectangular cavity containing two immiscible liquid layers is considered in the absence of gravity. The introduction of a more viscous encapsulant layer leads to a significant reduction in the intensity of the thermocapillary flow within the encapsulated layer. Interface deformations are small when the contact line of the interface is pinned on the solid boundaries. The higher viscosity of the encapsulant layer gives rise to a larger pressure gradient in that layer, thereby resulting in interface deformations that are qualitatively different from those observed at the free surface in the absence of the encapsulant layer. The flow pattern in the encapsulated layer and the resulting interface deformations are strongly dependent on both the thickness and the viscosity of the encapsulant layer. It is shown that the flow within the encapsulated layer may be closely approximated by simply considering the single-layer problem with a modified stress condition at the interface. The modified tangential stress balance for the effective single-layer model is derived based on asymptotic results for small-aspect-ratio double-layer systems and the insight gained from double-layer computations for finite-aspect-ratio systems. It is shown that the single-layer model accurately predicts the flow in the double-layer system even for large aspect-ratios.

Surfactant-assisted spreading of a liquid drop on a smooth solid surface

Axisymmetric spreading of a liquid drop covered with an insoluble surfactant monolayer on a smooth solid substrate is numerically investigated. As the drop spreads, the adsorbed surfactant molecules are constantly redistributed along the air-liquid interface by convection and diffusion, leading to nonuniformities in surface tension along the interface. The resulting Marangoni stresses affect the spreading rate by altering the surface flow and the drop shape. In addition, surfactant accumulation in the vicinity of the moving contact line affects the spreading rate by altering the balance of line forces. Two different models for the constitutive relation at the moving contact line are used, in conjunction with a surface equation of state based on the Frumkin adsorption framework, to probe the surfactant influence. The coupled evolution equations for the drop shape and monolayer concentration profile are integrated using a pseudospectral method to determine the rate of surfactant-assisted spreading over a wide range of the dimensionless parameters governing the spreading process. The insoluble monolayer enhances spreading through two mechanisms; a reduction in the equilibrium contact angle, and an increase in the magnitude of the radial pressure gradient within the drop due to the formation of positive surface curvature near the moving contact line. Both mechanisms are driven by the accumulation of surfactant at the contact line due to surface convection. Although the Marangoni stresses induced at the air-liquid interface reduce the rate of spreading during the initial stages of spreading, their retarding effect is overwhelmed by the favorable effects of the aforementioned mechanisms to lead to an overall enhancement in the rate of spreading in most cases. The spreading rate is found to be higher for bulkier surfactants with stronger repulsive interactions. With the exception of monolayers with strong cohesive interactions which tend to retard the spreading process, the overall effect of an insoluble monolayer is to increase the rate of drop spreading. Simulation results for small Bond numbers indicate the existence of a power-law region for the time-dependence of the basal radius of the drop, consistent with experimental measurements.

Surfactant Effect on the Buoyancy-Driven Motion of Bubbles and Drops in a Tube

The effect of surfactants on the buoyancy-driven motion of bubbles and drops in a vertical tube is experimentally examined. The terminal velocities of fluid particles are measured and their steady shapes are quantitatively characterized in systems with various bulk-phase concentrations of surfactant. In the case of air bubbles, the presence of surfactant retards the motion of small bubbles due to the development of adverse Marangoni stresses, whereas it enhances the motion of large bubbles by allowing them to deform away from the tube wall more easily. For viscous drops, the surfactant-enhanced regime of particle motion becomes more pronounced in the sense that the terminal velocity becomes more sensitive to surfactant concentration, whereas the surfactant effect in the surfactant-retarded regime becomes weaker.

Grade progression and regression in recurrent urothelial cancer

PURPOSE

Recurrent urothelial cancers are reported to have characteristics similar to those of the primary tumor, with 10% to 25% of low grade tumors recurring as high grade disease. We determined how often grade progression and regression occur and whether abnormalities in p53 protein expression in original tumors are preserved in recurrences.

MATERIALS AND METHODS

Two groups of patients treated for recurrent stages Ta/T1 urothelial bladder cancers with at least 1 tumor-free examination between the index and recurrent tumors were reviewed. Group 1 included 115 patients in whom the first available tumor was compared with the last recurrence and group 2 included 42 in whom the initial tumor was compared with the first recurrence. Immunohistochemical analysis of p53 expression was performed on a subset of 34 tumor pairs.

RESULTS

In group 1, 33 grade 3 tumors (45%) recurred as grade 1 or 2 tumors, while 9 of 82 grades 1 and 2 tumors (11%) recurred as grade 3 tumors. Five of 7 group 2 grade 3 tumors (71%) recurred as grade 1 or 2 disease, while 1 of 35 grades 1 and 2 tumors (3%) recurred as grade 3 disease. In the 34 pairs studied immunohistochemically 6 of 14 grade 3 tumors recurred at lower grades. Nuclear p53 over expression occurred in 21 index tumors (12 of 14 grade 3, 8 of 17 grade 2 and 1 of 3 grade 1) and in 9 recurrences (6 of 10 grade 3, 2 of 17 grade 2 and 1 of 7 grade 1). Only 7 of 21 p53 positive and 2 of 12 p53 negative index tumors were p53 positive on recurrence.

CONCLUSIONS

While progression from low to high grade occurred in less than 15% of patients, grade regression was observed in almost 50%. The loss of p53 positivity in regressing tumors indicates that these recurrences are molecularly distinct from the corresponding initial tumor.

A Numerical Study of the Effect of Insoluble Surfactants on the Stability of a Viscous Drop Translating in a Hele–Shaw Cell

A circular drop is a linearly stable solution for the buoyancy-driven motion of drops in a Hele-Shaw cell [Gupta et al. J. Colloid Interface Sci.218(1), 338 (1999)]. In the absence of surface-active agents, an initially prolate drop always goes to a steady circular shape while initially oblate drops exhibit complex dynamics [Gupta et al. J. Colloid Interface Sci.222, 107 (2000)]. In this study, the effect of insoluble surfactant impurities on the critical conditions for drop breakup is explored by using the Langmuir adsorption framework in conjunction with a physically based expression for the depth-averaged tangential stress exerted on a two-phase interface in a Hele-Shaw cell. It is shown that the presence of surfactants can have both a stabilizing and a destabilizing effect on the shape of the drop, depending on the Bond number, the magnitude of the initial perturbation, and the strength of surface convection. Similar to the clean drop dynamics, two marginally stable branches are found. Increasing the surface Peclet number results in the stabilization of the main branch while the secondary branch shifts to higher Bond numbers. The mode of breakup is also found to be strongly influenced by the strength of surface convection.

Stability of the Shape of a Viscous Drop under Buoyancy-Driven Translation in a Hele-Shaw Cell

It has been shown (N. R. Gupta, A. Nadim, H. Haj-Hariri, and A. Borhan, J. Colloid Interface Sci. 218, 338 1999) that a circular drop translating in a Hele-Shaw cell under the action of gravity is linearly stable for nonzero interfacial tension. In this paper, we use the boundary integral method to examine the nonlinear evolution of the shape of initially noncircular drops translating in a Hele-Shaw cell. For prolate initial deformations, it is found that the drop reverts to a circular shape for all finite Bond numbers considered. Initially oblate drops, on the other hand, are found to become unstable and break up if the initial shape perturbation is of sufficiently large magnitude. The critical conditions for the onset of drop breakup are examined in terms of the magnitude of the initial deformation as a function of Bond number. Two branches of marginal stability are identified and the effects of viscosity ratio and asymmetric initial perturbations on the stability diagram are discussed. Copyright 2000 Academic Press.