scholarly journals Effects of surfactant on propagation and rupture of a liquid plug in a tube

2019 ◽  
Vol 872 ◽  
pp. 407-437 ◽  
Author(s):  
M. Muradoglu ◽  
F. Romanò ◽  
H. Fujioka ◽  
J. B. Grotberg

Surfactant-laden liquid plug propagation and rupture occurring in lower lung airways are studied computationally using a front-tracking method. The plug is driven by an applied constant pressure in a rigid axisymmetric tube whose inner surface is coated by a thin liquid film. The evolution equations of the interfacial and bulk surfactant concentrations coupled with the incompressible Navier–Stokes equations are solved in the front-tracking framework. The numerical method is first validated for a surfactant-free case and the results are found to be in good agreement with the earlier simulations of Fujioka et al. (Phys. Fluids, vol. 20, 2008, 062104) and Hassan et al. (Intl J. Numer. Meth. Fluids, vol. 67, 2011, pp. 1373–1392). Then extensive simulations are performed to investigate the effects of surfactant on the mechanical stresses that could be injurious to epithelial cells, such as pressure and shear stress. It is found that the liquid plug ruptures violently to induce large pressure and shear stress on airway walls and even a tiny amount of surfactant significantly reduces the pressure and shear stress and thus improves cell survivability. However, addition of surfactant also delays the plug rupture and thus airway reopening.

Author(s):  
Payam Sharifi ◽  
Asghar Esmaeeli

Most of the studies conducted so far on EHD-driven instability of superimposed fluids have been concerned with liquid layers of modest depths. In many applications, however, the liquid layers can be very thin. Since the dynamics in thin films is generally governed by lubrication equations rather than full Navier-Stokes equations, it is expected that the interface dynamics will be quite different from that of the liquids with modest depths. The objective of this study is to explore the effect of initial liquid thickness on the dynamics of the phase boundary. To do this end, we perform Direct Numerical Simulations (DNS) using a front tracking/finite difference scheme, in conjunction with Taylor’s leaky dielectric model. For the physical parameters used here, it is shown that for sufficiently thick liquid layers, the interface instability leads to formation of liquid columns that merge together to form a big column. However, for thin layers, the interactions between the columns are weaker and lead to a short and a longer column that are connected by a thin liquid film.


1994 ◽  
Vol 271 ◽  
pp. 1-16 ◽  
Author(s):  
Peter Y. Huang ◽  
Jimmy Feng ◽  
Daniel D. Joseph

We do a direct two-dimensional finite-elment simulation of the Navier–Stokes equations and compute the forces which turn an ellipse settling in a vertical channel of viscous fluid in a regime in which the ellipse oscillates under the action of vortex shedding. Turning this way and that is induced by large and unequal values of negative pressure at the rear separation points which are here identified with the two points on the back face where the shear stress vanishes. The main restoring mechanism which turns the broadside of the ellipse perpendicular to the fall is the high pressure at the ‘stagnation point’ on the front face, as in potential flow, which is here identified with the one point on the front face where the shear stress vanishes.


2015 ◽  
Vol 1105 ◽  
pp. 105-109 ◽  
Author(s):  
Varvara Yu. Gordeeva ◽  
Andrey V. Lyushnin

Evaporation of a thin layer of a polar liquid (water) having a free surface and located on a solid substrate is investigated. A surfactant is solved in the liquid film. The surface tension is a linear function of the surface concentration of the surfactant. The surface energy of the solid-liquid interface is a nonmonotonic function of the layer thickness and is the sum of the Van der Waals interaction and the specific interaction of the double electric layer on the interface. The effect of the solvable surfactant on the dynamics of the propagation of the evaporation front in the thin liquid film is analyzed in the long-wave approximation in the system of Navier-Stokes equations.


Author(s):  
Man-Woong Heo ◽  
Tae-Wan Seo ◽  
Chung-Suk Lee ◽  
Kwang-Yong Kim

This paper presents a parametric study to investigate the aerodynamic and aeroacoustic characteristics of a side channel regenerative blower. Flow analysis in the side channel blower was carried out by solving three-dimensional steady and unsteady Reynolds-averaged Navier-Stokes equations with the shear stress transport turbulence closure. Aeroacoustic analysis was conducted by solving the variational formulation of Lighthill’s analogy on the basis of the aerodynamic sources extracted from the unsteady flow analysis. The height and width of the blade and the angle between inlet and outlet ports were selected as three geometric parameters, and their effects on the aerodynamic and aeroacoustic performances of the blower have been investigated. The results showed that the aerodynamic and aeroacoustic performances were enhanced by decreasing height and width of blade. It was found that angle between inlet and outlet ports significantly influences the aerodynamic and aeroacoustic performances of the blower due to the stripper leakage flow.


Author(s):  
D.V. Singh ◽  
R. Sinhasan ◽  
S.P. Tayal

Additives are extensively used in the commercial lubricants to improve their specific qualities. These lubricants are therefore non-Newtonian and their nonlinear relations between shear stress and shear strain rate are generally represented by cubic shear stress laws. The Navier-Stokes equations and the continuity equation in clindrical coordinates, representing the flow-field in the clearance space of each lobe of the three-lobe hydrodynamic journal bearings having Newtonian fluids, are solved by the finie element method using Galerkin’s technique. The solution for non-Newtonian lubricants is obtained by an iteration technique modifying the viscosity term in each iteration. The static performance characteristics have been obtained for both Newtonian and the non-Newtonian lubricants. The load capacity and friction of the bearing decrease with increase in the nonlinearity of the lubricant whereas the end flow is relatively unaffected.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Filipe S. Pereira ◽  
Luís Eça ◽  
Guilherme Vaz

The importance of the turbulence closure to the modeling accuracy of the partially-averaged Navier–Stokes equations (PANS) is investigated in prediction of the flow around a circular cylinder at Reynolds number of 3900. A series of PANS calculations at various degrees of physical resolution is conducted using three Reynolds-averaged Navier–Stokes equations (RANS)-based closures: the standard, shear-stress transport (SST), and turbulent/nonturbulent (TNT) k–ω models. The latter is proposed in this work. The results illustrate the dependence of PANS on the closure. At coarse physical resolutions, a narrower range of scales is resolved so that the influence of the closure on the simulations accuracy increases significantly. Among all closures, PANS–TNT achieves the lowest comparison errors. The reduced sensitivity of this closure to freestream turbulence quantities and the absence of auxiliary functions from its governing equations are certainly contributing to this result. It is demonstrated that the use of partial turbulence quantities in such auxiliary functions calibrated for total turbulent (RANS) quantities affects their behavior. On the other hand, the successive increase of physical resolution reduces the relevance of the closure, causing the convergence of the three models toward the same solution. This outcome is achieved once the physical resolution and closure guarantee the precise replication of the spatial development of the key coherent structures of the flow.


1985 ◽  
Vol 154 ◽  
pp. 357-375 ◽  
Author(s):  
J. A. C. Humphrey ◽  
H. Iacovides ◽  
B. E. Launder

The paper reports numerical solutions to a semi-elliptic truncation of the Navier–Stokes equations for the case of developing laminar flow in circular-sectioned bends over a range of Dean numbers. The ratios of bend radius to pipe radius are 7:1 and 20:1, corresponding with the configurations examined experimentally by Talbot and his co-workers in recent years. The semi-elliptic treatment facilitates a much finer grid than has been possible in earlier studies. Numerical accuracy has been further improved by assuming radial equilibrium over a thin sublayer immediately adjacent to the wall and by re-formulating the boundary conditions at the pipe centre.Streamwise velocity profiles at Dean numbers of 183 and 565 are in excellent agreement with laser-Doppler measurements by Agrawal, Talbot & Gong (1978). Good, albeit less complete, accord is found with the secondary velocities, though the differences that exist may be mainly due to the difficulty of making these measurements. The paper provides new information on the behaviour of the streamwise shear stress around the inner line of symmetry. Upstream of the point of minimum shear stress, our numerical predictions display a progressive shift towards the result of Stewartson, Cebici & Chang (1980) as the Dean number is successively raised. Downstream of the minimum, however, in contrast with the monotonic approach to an asymptotic level reported by Stewartson, the numerical solutions display a damped oscillatory behaviour reminiscent of those from Hawthorne's (1951) inviscid-flow calculations. The amplitude of the oscillation grows as the Dean number is raised.


2014 ◽  
Vol 348 ◽  
pp. 88-93 ◽  
Author(s):  
S. Aggoune ◽  
El Hachemi Amara

We consider in the present work the fusion laser cutting of stainless steel sheets under a nitrogen laminar gas jet. The molten metal is treated as a laminar and steady viscous incompressible fluid. The mathematical model describing our problem is set in terms of Navier-Stokes equations, solved numerically using the finite differences method, where the effect of the gas jet velocity on the molten boundary layer is considered. The generated shear stress occurring on the gas-liquid interface and its contribution in the momentum is carried out, and it is found that when the skin friction and the shear stress decrease, the thickness and the velocity at the edge of the molten boundary layer increase along the kerf surface. The layer thickness reduces when the assisting gas velocity is increased.


2021 ◽  
Vol 24 (2) ◽  
pp. first
Author(s):  
Truong V. Vu ◽  
Vinh T. Nguyen ◽  
Phan H. Nguyen ◽  
Nang X. Ho ◽  
Binh D. Pham ◽  
...  

Introduction: Compound fluid filaments appear in many applications, e.g., drug delivery and processing or microfluidic systems. This paper focuses on the numerical simulation of an incompressible, immiscible, and Newtonian fluid for the contraction process of a fluid compound filament by solving the Navier-Stokes equations. The front-tracking method is used to solve this problem, which uses connected segments (Lagrangian grid) that move on a fixed grid (Eulerian grid) to represent the interface between the liquids. Methods: The interface points are advected by the velocity interpolated from those of the fixed grid using the area weighting function. The coordinates of the interface points are used to construct the indicators specifying the different fluids and compute the interfacial tension force. Results: The simulation results show that under the effects of the interfacial tension, the capsuleshaped filament can transform into a spherical compound droplet (i.e., non-breakup) or can break up into smaller spherical compound and simple droplets (i.e., breakup). When the density ratio of the outer to middle fluids increases, the filament changes from non-breakup to breakup upon contraction. Conclusion: Increasing the density ratio enhances the breakup of the compound filament during contraction. The breakup is also promoted by increasing the initial length of the filament.


2016 ◽  
Vol 38 (3) ◽  
pp. 193-204
Author(s):  
Vu Van Truong

In this paper, direct numerical simulations are presented for solidification with the effects of density difference between the solid and liquid phases. A front-tracking method is used. The solidification front, i.e. the solid-liquid interface separating solid and liquid, is represented by connected elements that move on a rectangular and stationary grid. The Navier-Stokes equations are solved by a projection method on the entire domain including the solid phase. An indicator function reconstructed from the front information is used to set the velocities in the solid phase to zero, and thus to enforce the no-slip condition at the interface. The method is validated through comparisons with exact solutions for one- and two-dimensional problems. The method is then used to simulate the solidification processes with the effects of volume change due to density difference


Sign in / Sign up

Export Citation Format

Share Document