scholarly journals Effects of Surface Tension and Yield Stress on Mucus Plug Rupture: A Numerical Study

2020 ◽  
Vol 142 (6) ◽  
Author(s):  
Yingying Hu ◽  
Francesco Romanò ◽  
James B. Grotberg

Abstract We study the effects of surface tension and yield stress on mucus plug rupture. A three-dimensional simplified configuration is employed to simulate mucus plug rupture in a collapsed lung airway of the tenth generation. The Herschel–Bulkley model is used to take into account the non-Newtonian viscoplastic fluid properties of mucus. Results show that the maximum wall shear stress greatly changes right prior to the rupture of the mucus plug. The surface tension influences mainly the late stage of the rupture process when the plug deforms greatly and the curvature of the mucus–air interface becomes significant. High surface tension increases the wall shear stress and the time needed to rupture since it produces a resistance to the rupture, as well as strong stress and velocity gradients across the mucus–air interface. The yield stress effects are pronounced mainly at the beginning. High yield stress makes the plug take a long time to yield and slows down the whole rupture process. When the effects induced by the surface tension and yield forces are comparable, dynamical quantities strongly depend on the ratio of the two forces. The pressure difference (the only driving in the study) contributes to wall shear stress much more than yield stress and surface tension per unit length. Wall shear stress is less sensitive to the variation in yield stress than that in surface tension. In general, wall shear stress can be effectively reduced by the smaller pressure difference and surface tension.

2021 ◽  
Vol 12 ◽  
Author(s):  
Yu Wang ◽  
Junwei Wang ◽  
Jing Peng ◽  
Mingming Huo ◽  
Zhiqiang Yang ◽  
...  

Patients with heart failure (HF) or undergoing cardiogenic shock and percutaneous coronary intervention require short-term cardiac support. Short-term cardiac support using a left ventricular assist device (LVAD) alters the pressure and flows of the vasculature by enhancing perfusion and improving the hemodynamic performance for the HF patients. However, due to the position of the inflow and outflow of the LVAD, the local hemodynamics within the aorta is altered with the LVAD support. Specifically, blood velocity, wall shear stress, and pressure difference are altered within the aorta. In this study, computational fluid dynamics (CFD) was used to elucidate the effects of a short-term LVAD for hemodynamic performance in a patient-specific aorta model. The three-dimensional (3D) geometric models of a patient-specific aorta and a short-term LVAD, Impella CP, were created. Velocity, wall shear stress, and pressure difference in the patient-specific aorta model with the Impella CP assistance were calculated and compared with the baseline values of the aorta without Impella CP support. Impella CP support augmented cardiac output, blood velocity, wall shear stress, and pressure difference in the aorta. The proposed CFD study could analyze the quantitative changes in the important hemodynamic parameters while considering the effects of Impella CP, and provide a scientific basis for further predicting and assessing the effects of these hemodynamic signals on the aorta.


Cellulose ◽  
2019 ◽  
Vol 27 (1) ◽  
pp. 141-156 ◽  
Author(s):  
Tuomas Turpeinen ◽  
Ari Jäsberg ◽  
Sanna Haavisto ◽  
Johanna Liukkonen ◽  
Juha Salmela ◽  
...  

Abstract The shear rheology of two mechanically manufactured microfibrillated cellulose (MFC) suspensions was studied in a consistency range of 0.2–2.0% with a pipe rheometer combined with ultrasound velocity profiling. The MFC suspensions behaved at all consistencies as shear thinning power law fluids. Despite their significantly different particle size, the viscous behavior of the suspensions was quantitatively similar. For both suspensions, the dependence of yield stress and the consistency index on consistency was a power law with an exponent of 2.4, similar to some pulp suspensions. The dependence of flow index on consistency was also a power law, with an exponent of − 0.36. The slip flow was very strong for both MFCs and contributed up to 95% to the flow rate. When wall shear stress exceeded two times the yield stress, slip flow caused drag reduction with consistencies higher than 0.8%. When inspecting the slip velocities of both suspensions as a function of wall shear stress scaled with the yield stress, a good data collapse was obtained. The observed similarities in the shear rheology of both the MFC suspensions and the similar behavior of some pulp fiber suspensions suggests that the shear rheology of MFC suspensions might be more universal than has previously been realized.


2020 ◽  
Vol 76 (1) ◽  
pp. 9-26
Author(s):  
Saeed Bahrami ◽  
Mahmood Norouzi

Cardiovascular disease is now under the influence of several factors that encourage researchers to investigate the flow of these vessels. Oscillation influences the blood circulation in the volume of red blood cells (RBC) strongly. Therefore, in this study, its effects have been considered on hemodynamic parameters in the elastic wall and coronary bifurcation. In this study, a 3D geometry of non-Newtonian and pulsatile blood circulation is considered in the left coronary artery bifurcation. The Casson model with various hematocrits is analyzed in elastic and rigid walls. The wall shear stress (WSS) cannot show the stenosis artery alone, therefore, the oscillatory shear index (OSI) is represented as a hemodynamic parameter of WSS individually of time. The results are determined using two-way fluid-structure interaction (FSI) coupling method using an arbitrary Lagrangian-Eulerian method. The most prominent difference in velocity happened in the bifurcation and at hematocrit 30 with yield stress 6.59E-04 Pa. The backflow and vortex flow in the LCx branch grown with increasing shear rates. The likelihood of plaque generation at the ending of the LM branch is observed in hematocrits 10 and 20, while the WSS magnitude is normal in the hematocrit 60 with the greatest yield stress in the bifurcation. The shear stress among the rigid and elastic models is the highest at the ending of the LM branch. The wall shear stress magnitude among the models decreased at most of 24.49% by dividing the flow. Time-independent results for models showed that there is the highest value of OSI at the bifurcation, which then quickly dropped.


Author(s):  
Khaled J. Hammad

Wall-bounded separating and reattaching flows are encountered in biological applications dealing with blood flows through arteries and prosthetic devices. Separated and reattached flow regions have been associated in the past with the most common arterial disease, atherosclerosis. Previous studies suggest that local wall shear stress (WSS) patterns affect the location and progression rate of atherosclerotic lesions. A parametric study is performed to investigate the influence of hemorheology on the wall shear stress distribution in a separated and reattached flow region. Recent hemorheological studies quantified and emphasized the yield stress and shear-thinning non-Newtonian characteristics of unadulterated human blood. Numerical solutions to the governing equations that account for yield stress and shear-thinning rheological effects are obtained. A low WSS region is observed around the flow reattachment point while a peak WSS always exists close to the vortex center. The yield shear-thinning hemorheological model always results in the highest observed peak WSS. The yield stress impact on WSS distribution is most pronounced in the case of severe restrictions to the flow.


Water ◽  
2018 ◽  
Vol 11 (1) ◽  
pp. 19 ◽  
Author(s):  
Dhruv Mehta ◽  
Adithya Thota Radhakrishnan ◽  
Jules van Lier ◽  
Francois Clemens

This article follows from a previous study by the authors on the computational fluid dynamics-based analysis of Herschel–Bulkley fluids in a pipe-bounded turbulent flow. The study aims to propose a numerical method that could support engineering processes involving the design and implementation of a waste water transport system, for concentrated domestic slurry. Concentrated domestic slurry results from the reduction in the amount of water used in domestic activities (and also the separation of black and grey water). This primarily saves water and also increases the concentration of nutrients and biomass in the slurry, facilitating efficient recovery. Experiments revealed that upon concentration, domestic slurry flows as a non-Newtonian fluid of the Herschel–Bulkley type. An analytical solution for the laminar transport of such a fluid is available in literature. However, a similar solution for the turbulent transport of a Herschel–Bulkley fluid is unavailable, which prompted the development of an appropriate wall function to aid the analysis of such flows. The wall function (called ψ 1 hereafter) was developed using Launder and Spalding’s standard wall function as a guide and was validated against a range of experimental test-cases, with positive results. ψ 1 is assessed for its sensitivity to rheological parameters, namely the yield stress, the fluid consistency index and the behaviour index and their impact on the accuracy with which ψ 1 can correctly quantify the pressure loss through a pipe. This is done while simulating the flow of concentrated domestic slurry using the Reynolds-Averaged Navier–Stokes (RANS) approach for turbulent flows. This serves to establish an operational envelope in terms of the rheological parameters and the average flow velocity within which ψ 1 is a must for accuracy. One observes that, regardless of the fluid behaviour index, ψ 1 is necessary to ensure accuracy with RANS models only in flow regimes where the wall shear stress is comparable to the yield stress within an order of magnitude. This is also the regime within which the concentrated slurry analysed as part of this research flows, making ψ 1 a requirement. In addition, when the wall shear stress exceeds the yield stress by more than one order (either due to an inherent lower yield stress or a high flow velocity), the regular Newtonian wall function proposed by Launder and Spalding is sufficient for an accurate estimate of the pressure loss, owing to the relative reduction in non-Newtonian viscosity as compared to the turbulent viscosity.


2011 ◽  
Vol 705 ◽  
pp. 258-279 ◽  
Author(s):  
Parsa Zamankhan ◽  
Brian T. Helenbrook ◽  
Shuichi Takayama ◽  
James B. Grotberg

AbstractWe study numerically the steady creeping motion of Bingham liquid plugs in two-dimensional channels as a model of mucus behaviour during airway reopening in pulmonary airways. In addition to flow analysis related to propagation of the plug, the stress distribution on the wall is studied for better understanding of potential airway epithelial cell injury mechanisms. The yield stress behaviour of the fluid was implemented through a regularized constitutive equation. The capillary number, $\mathit{Ca}$, and the Bingham number, $\mathit{Bn}$, which is the ratio of the yield stress to a characteristic viscous stress, varied over the ranges 0.025–0.1 and 0–1.5, respectively. For the range of parameters studied, it was found that, while the yield stress reduces the magnitude of the shearing along the wall, it can magnify the amplitude of the wall shear stress gradient significantly, and also it can elevate the magnitude of the wall shear stress and wall pressure gradient up to 30 % and 15 %, respectively. Therefore, the motion of mucus plugs can be more damaging to the airway epithelial cells due to the yield stress properties of mucus. The yield stress also modifies the profile of the plug where the amplitude of the capillary waves at the leading meniscus decreases with increase in $\mathit{Bn}$. Other findings are that: the thickness of the static film increases with increasing $\mathit{Bn}$; the driving pressure difference increases linearly with $\mathit{Bn}$; and increasing $\mathit{Bn}$ extends any wall stagnation point beneath the leading meniscus to an unyielded line segment beneath the leading meniscus. With an increase in $\mathit{Bn}$, the unyielded areas appear and grow in the adjacent wall film as well as the core region of the plug between the two menisci. The plug length, ${L}_{P} $, mostly modifies the topology of the yield surfaces. It was found that the unyielded area in the core region between the two menisci grows as the plug length decreases. The very short Bingham plug behaves like a solid lamella. In all computed liquid plugs moving steadily, the von Mises stress attains its maximum value near the interface of the leading meniscus in the transition region. For Bingham plugs moving very slowly, $\mathit{Ca}\ensuremath{\rightarrow} 0$, the driving pressure is non-zero.


1971 ◽  
Vol 93 (2) ◽  
pp. 102-108 ◽  
Author(s):  
L. S. Galowin ◽  
M. J. Desantis

A theoretical investigation was conducted to obtain velocity, pressure, and shear stress distributions for incompressible, steady, fully developed, laminar flow through a cylinder with a uniformly porous wall. Ejection/injection at the walls results from the pressure difference across the porous wall. Fluid flow phenomena in porous tubes and ducts have previously been investigated with the velocity prescribed as the boundary condition at the wall. An accurate wall condition must account for the variable wall velocity being dependent upon the pressure difference across the wall, the properties of the fluid, the thickness and the permeability of the structure. An integral momentum technique was employed to reduce the axisymmetric Navier-Stokes equations in cylindrical coordinates to a nonlinear, second-order ordinary differential equation with appropriate boundary conditions. The velocity condition at the wall was established for the ejection/injection at the surface resulting from the pressure difference across the porous wall derived from Darcy’s law. Numerical solutions were obtained for a range of axial flow Reynolds numbers, wall permeabilities, and initial pressure difference across the porous wall. The calculated static pressure variation in the axial flow direction, the velocity components, and the wall shear stress are presented. For the case of fluid ejection, the results of the analysis show that the wall shear stress and static pressure decrease in the axial flow direction. The rates of decrease are functions of the wall porosity, initial pressure gradient across the wall, and inlet flow Reynolds number. The present analysis treats the realistic problem of flow adjustment to the condition where zero pressure differential across the porous wall occurs (the normal wall velocity vanishes). Previous models are based upon the assumptions of constant radial velocity at the wall and/or prescribed wall shear stress without taking into account the pressure drop through the wall. Such assumptions imply that a variable pressure exists external to the pipe, or that the pipe has walls of variable permeability and thickness rather than the hypothesized condition that the pipe has uniformly porous walls. For one set of boundary conditions it is shown that the outflow through the walls completely discharges the entering flow. As a result no far downstream axial flow occurs. Such effects were not previously discussed by other investigators. For other sets of boundary conditions reductions in centerline velocity and shear stress occur.


Author(s):  
Khaled J. Hammad

The influence of inflow conditions and human blood rheology on the wall shear stress distribution in a confined separated and reattached flow region is investigated. The governing mass and momentum conservation equations along with the Herschel-Bulkley rheological model are solved numerically using a finite-difference scheme. A parametric study is performed to reveal the influence of uniform and fully-developed inflow velocity profiles on the wall shear stress (WSS) characteristics using hemorheological models that account for the yield stress and shear-thinning non-Newtonian characteristics of human blood. The highest WSS or WSSmax, is always observed inside the flow separation region at a location corresponding to that of the corner vortex center. Uniform inflow results in higher WSSmax values in comparison with fully-developed inflow for moderate upstream flow restrictions. The opposite trend is observed for severe flow restrictions. Uniform inflow always results in smaller flow separation regions and WSSmax values at locations closer to the flow restriction plane. The yield shear-thinning hemorheological model always results in the highest observed peak WSS. The yield stress impact on WSS distribution is most pronounced in the case of severe restrictions to the flow.


2018 ◽  
Vol 24 (1) ◽  
pp. 19-24
Author(s):  
Sertaç Çadırcı ◽  
Ufuk Demir ◽  
Semra Zuhal Birol ◽  
Levent Trabzon ◽  
Hasan Güneş

2018 ◽  
Vol 141 (1) ◽  
Author(s):  
Kofi Freeman K. Adane ◽  
Martin Agelin-Chaab

In this study, a qualitative assessment of transitional velocity engineering models for predicting non-Newtonian slurry flows in a horizontal pipe was performed using data from a wide range of pipe diameters (25–268 mm). In addition, the gamma theta transition model was used to compute selected flow conditions. These models were used to predict transitional velocities in large pipe diameters (up to 420 mm) for slurries. In general, it was observed that most of the current engineering models predict transitional velocities conservatively. Based on the gamma theta transition model results, for large Hedström numbers (He ≳ 105), other methods should be used to predict transitional velocities if a change in the pipe diameter (scale-up) results in an order of magnitude increase in the He value. It was also found that the gamma theta transition model predicted a laminar flow condition in the fully developed region for flow conditions with a small plug region (low-yield stress-to-wall shear stress ratio), which is contrary to what has been observed in some experiments. This is attributed to the local fluid rheological parameters values, which might be different from those reported. However, the gamma theta transition model results are in good agreement with the experimental data for flow conditions that have a large plug region (high-yield stress-to-wall shear stress ratio).


Sign in / Sign up

Export Citation Format

Share Document