Experimental-computational study of fibrous particle transport and deposition in a bifurcating lung model

Experiments carried out using a lung model with a single horizontal bifurcation under different steady inhalation conditions explored the orientation of depositing carbon fibers, and particle deposition fractions. The orientations of deposited fibers were obtained from micrographs. Specifically, the effects of the sedimentation parameter (γ), fiber length, and flow rate on orientations were analyzed. Our results indicate that gravitational effect on deposition cannot be neglected for 0.0228 < γ < 0.247. The absolute orientation angle of depositing fibers decreased linearly with increasing γ for values 0.0228 < γ < 0.15. Correspondence between Stokes numbers and γ suggests these characteristics can be used to estimate fiber deposition in the lower airways. Computer simulations with sphere-equivalent diameter models for the fibers explored deposition efficiency vs. Stokes number. Using the volume-equivalent diameter model, our experimental data for the horizontal bifurcation were replicated. Results for particle deposition using a lung model with a vertical bifurcation indicate that body position also affects deposition.     

Unsteady gravity-driven slender rivulets of a power-law fluid

Unsteady gravity-driven flow of a thin slender rivulet of a non-Newtonian power-law fluid on a plane inclined at an angle α to the horizontal is considered. Unsteady similarity solutions are obtained for both converging sessile rivulets (when 0 < α < π/2) in the case x < 0 with t < 0, and diverging pendent rivulets (when π/2 < α < π) in the case x > 0 with t > 0, where x denotes a coordinate measured down the plane and t denotes time. Numerical and asymptotic methods are used to show that for each value of the power-law index N there are two physically realisable solutions, with cross-sectional profiles that are ‘single-humped’ and ‘double-humped’, respectively. Each solution predicts that at any time t the rivulet widens or narrows according to |x | ^{(2N+1)/2(N+1)} and thickens or thins according to |x | ^N/(N+1) as it flows down the plane; moreover, at any station x, it widens or narrows according to |t | ^?N/2(N+1) and thickens or thins according to |t | ^?N/(N+1). The length of a truncated rivulet of fixed volume is found to behave according to |t | ^N/(2N+1).     

Evaluation of subgrid-scale models in large-eddy simulation of flow past a two-dimensional block

Large-eddy simulations of flow past a two-dimensional (2D) block were performed to evaluate four subgrid-scale (SGS) models: (i) the traditional Smagorinsky model, (ii) the Lagrangian dynamic model, (iii) the Lagrangian scale-dependent dynamic model, and (iv) the modulated gradient model. An immersed boundary method was employed to simulate the 2D block boundaries on a uniform Cartesian grid. The sensitivity of the simulation results to grid refinement was investigated by using four different grid resolutions. The velocity streamlines and the vertical profiles of the mean velocities and variances were compared with experimental results. The modulated gradient model shows the best overall agreement with the experimental results among the four SGS models. In particular, the flow recirculation, the reattachment position and the vertical profiles are accurately reproduced with a relative coarse grid resolution of (N_x × N_y × N_z=) 160 × 40 × 160 (n_x × n_z = 13 × 16 covering the block). Besides the modulated gradient model, the Lagrangian scale-dependent dynamic model is also able to give reasonable prediction of the flow statistics with some discrepancies compared with the experimental results. Relatively poor performance by the Lagrangian dynamic model and the Smagorinsky model is observed, with simulated recirculating patterns that differ from the measured ones. Analysis of the turbulence kinetic energy (TKE) budget in this flow shows evidence of a strong production of TKE in the shear layer that forms as the flow is deflected around the block.     

LES of turbulent flow in a concentric annulus with rotating outer wall

Fully-developed turbulent flow in a concentric annulus, r₁/r₂ = 0.5, Re_h = 12,500, with the outer wall rotating at a range of rotation rates N = U_θ,wall/U_b from 0.5 up to 4 is studied by large-eddy simulations. The focus is on the effects of moderate to very high rotation rates on the mean flow, turbulence statistics and eddy structure. For N up to ∼2, an increase in the rotation rate dampens progressively the turbulence near the rotating outer wall, while affecting only mildly the inner-wall region. At higher rotation rates this trend is reversed: for N = 2.8 close to the inner wall turbulence is dramatically reduced while the outer wall region remains turbulent with discernible helical vortices as the dominant turbulent structure. The turbulence parameters and eddy structures differ significantly for N = 2 and 2.8. This switch is attributed to the centrifuged turbulence (generated near the inner wall) prevailing over the axial inertial force as well as over the counteracting laminarizing effects of the rotating outer wall. At still higher rotation, N = 4, the flow gets laminarized but with distinct spiralling vortices akin to the Taylor–Couette rolls found between the two counter-rotating cylinders without axial flow, which is the limiting case when N approaches to infinity. The ratio of the centrifugal to axial inertial forces, Ta/Re² ∝ N² (where Ta is the Taylor number) is considered as a possible criterion for defining the conditions for the above regime change.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

Neutrality in the case of N-phase elliptical inclusions with internal uniform hydrostatic stresses

A novel method is proposed to design neutral N-phase (N ? 3) elliptical inclusions with internal uniform hydrostatic stresses. We focus on the study of the internal and external stress states of an N-phase elliptical inclusion which is bonded to an infinite matrix through (N ? 2) interphase layers. The interfaces of the N-phase elliptical inclusion are (N ? 1) confocal ellipses. The design of the resulting overall composite material consists of four stages: (i) an inner perfectly bonded interphase/inclusion interface which is necessary to make the internal uniform stress state hydrostatic; (ii) outer imperfect interphase layers properly designed to make the coated inclusion harmonic (i.e., the uniform mean stress of the original field within the matrix is unperturbed); (iii) the aspect ratio of the elliptic inclusion uniquely chosen for a given material and thickness parameters to make the resulting coated inclusion neutral (i.e., the prescribed uniform stress field in the matrix remains undisturbed); and finally (iv) the derivation of a simple condition relating the remote uniform stresses and the thickness parameters of the (N ? 2) interphase layers for given material parameters which lead to internal uniform hydrostatic stresses. We note that another interesting feature of the present results is that the mean stress is found to be constant within each interphase layer, and the hoop stress in the innermost interphase layer is uniform along the entire interphase/inclusion interface.     

Extrudate swell of linear and branched polyethylenes: ALE simulations and comparison with experiments

Extrudate swell is a common phenomenon observed in the polymer extrusion industry. Accurate prediction of the dimensions of an extrudate is important for appropriate design of dies for profile extrusion applications. Prediction of extrudate swell has been challenging due to (i) difficulties associated with accurate representation of the constitutive behavior of polymer melts, and (ii) difficulties associated with the simulation of free surfaces, which requires special techniques in the traditionally used Eulerian framework. In a previous work we had argued that an Arbitrary Lagrangian Eulerian (ALE) based finite element formulation may have advantages in simulating free surface deformations such as in extrudate swell. In the present work we reinforce this argument by comparing our ALE simulations with experimental data on the extrudate swell of commercial grades of linear polyethylene (LLDPE) and branched polyethylene (LDPE). Rheological behavior of the polymers was characterized in shear and uniaxial extensional deformations, and the data was modeled using either the Phan–Thien Tanner (PTT) model or the eXtended Pom–Pom (XPP) model. Additionally, flow birefringence and pressure drop measurements were done using a 10:1 contraction–expansion (CE) slit geometry in a MultiPass Rheometer. Simulated pressure drop and contours of the principal stress difference were compared with experimental data and were found to match well. This provided an independent test for the accuracy of the ALE code and the constitutive equations for simulating a processing-like flow. The polymers were extruded from long (L/D = 30) and short (L/D = 10) capillaries dies at 190 °C. ALE simulations were performed for the same extrusion conditions and the simulated extrudate swell showed good agreement with the experimental data.     

Finite-length effects on dynamical behavior of rod-like particles in wall-bounded turbulent flow

Combined Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) measurements have been performed in dilute suspensions of rod-like particles in wall turbulence. PIV results for the turbulence field in the water table flow apparatus compared favorably with data from Direct Numerical Simulations (DNS) of channel flow turbulence and the universality of near-wall turbulence justified comparisons with DNS of fiber-laden channel flow. In order to examine any shape effects on the dynamical behavior of elongated particles in wall-bounded turbulent flow, fibers with three different lengths but the same diameter were used. In the logarithmic part of the wall-layer, the translational fiber velocity was practically unaffected by the fiber length l. In the buffer layer, however, the fiber dynamics turned out to be severely constrained by the distance z to the wall. The short fibers accumulated preferentially in low-speed areas and adhered to the local fluid speed. The longer fibers (l/z > 1) exhibited a bi-modal probability distribution for the fiber velocity, which reflected an almost equal likelihood for a long fiber to reside in an ejection or in a sweep. It was also observed that in the buffer region, high-speed long fibers were almost randomly oriented whereas for all size cases the slowly moving fibers preferentially oriented in the streamwise direction. These phenomena have not been observed in DNS studies of fiber suspension flows and suggested l/z to be an essential parameter in a new generation of wall-collision models to be used in numerical studies.     

Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

Experimental study on the effect of drag reducing polymer on flow patterns and drag reduction in a horizontal oil–water flow

In this study, a HMW anionic co-polymer of 40:60 wt/wt NaAMPS/acrylamide was used as a drag reducing polymer (DRP) for oil–water flow in a horizontal 25.4 mm ID acrylic pipe. The effect of polymer concentration in the master solution and after injection in the main water stream, oil and water velocities, and pipe length on drag reduction (DR) was investigated. The injected polymer had a noticeable effect on flow patterns and their transitions. Stratified and dual continuous flows extended to higher superficial oil velocities while annular flow changed to dual continuous flow. The results showed that as low as 2 ppm polymer concentration was sufficient to create a significant drag reduction across the pipe. DR was found to increase with polymer concentration increased and reached maximum plateau value at around 10 ppm. The results showed that the drag reduction effect tends to increase as superficial water velocity increased and eventually reached a plateau at U_sw of around 1.3 m/s. At U_sw > 1.0 m/s, the drag reduction decreased as U_so increased while at lower water velocities, drag reduction is fluctuating with respect to U_so. A maximum DR of about 60% was achieved at U_so = 0.14 m/s while only 45% was obtained at U_so = 0.52 m/s. The effectiveness of the DRP was found to be independent of the polymer concentration in the master solution and to some extent pipe length. The friction factor correlation proposed by Al-Sarkhi et al. (2011) for horizontal flow of oil–water using DRPs was found to underpredict the present experimental pressure gradient data.     

Large-eddy simulation of thermal fatigue in a mixing tee

Temperature fluctuations occur due to thermal mixing of hot and cold streams in the T-junctions of the piping system in nuclear power plants, which may cause thermal fatigue of piping system. In this paper, three-dimensional, unsteady numerical simulations of coolant temperature fluctuations at a mixing T-junction of equal diameter pipes were performed using the large eddy simulation (LES) turbulent model. The experiments used in this paper to benchmark the simulations were performed by Hitachi Ltd. The calculated normalized mean temperatures and fluctuating temperatures are in good agreement with the measurements. The influence of the time-step ranging from 100 Hz to 1000 Hz on the numerical simulation results was explored. The simulation results indicate that all the results with different frequencies agree well with the experimental data. Finally, the attenuation of fluctuation of fluid temperature was also investigated. It is found that, drastic fluctuation occurs within the range of less than L/D = 4.0; the fluctuation of fluid temperature does not always attenuate from the pipe center to the wall due to the continuous generation of vortexes. At the top wall, the position of L/D = 1.5 has a minimum normalized mean temperature and a peak value of root-mean square temperature, whereas at the bottom wall, the position having the same characteristics is L/D = 2.0.     

Image decomposition method for the analysis of a mixed dislocation in a general multilayer

The elastic solutions for a mixed dislocation in a general multilayer with N dissimilar anisotropic layers are obtained via a generalized image decomposition method. The original problem is decomposed into N homogeneous subproblems with strategically placed continuously distributed image (virtual) dislocations which satisfy the consistency conditions for degenerate N − M (M < N) layer problems. The image dislocations are used to satisfy the interface or free surface conditions, and represent the unknowns of the problem. The resulting singular Cauchy integral equations are transformed into non-singular Fredholm integral equations of the second kind using certain H- and I-integral transforms. The Fredholm integral equations are then solved via the classical Nyström method. The general decomposition and the elimination of all singular integrals yield an exact formulation of the problem; the approximation arises only in the Nyström method. The dislocation mixity and the number of layers dissimilar in thickness and elastic anisotropy can be handled without difficulty, constrained only by the number of linear algebraic equations in the Nyström method for large N. For the numerical study, image forces on a dislocation in two- and three-layer systems are calculated. The accuracy of the results is verified by checking the boundary conditions and by comparison with previous results. The dependence of the image force on the dislocation position and mixity, and on the layer thicknesses and elastic anisotropies, is also illustrated via numerical investigations.     

The degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary

This paper investigates the degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary. Inside the boundary, there are many voids with arbitrary configurations. The problem is analyzed with a relevant homogenous BIE (boundary integral equation). It is assumed that all the inner void boundary tractions are equal to zero, and tractions on the outer circular boundary are constant. Therefore, all the integrations in BIE are performed on the outer circular boundary only. By using the relation z * conjg(z) = a * a, or conjg(z) = a * a/z on the circular boundary with radius a, all integrals can be reduced to an integral for complex variable and they can be integrated in closed form. The degenerate scale a = 1 is found in the Laplace equation and in plane elasticity regardless of the void configuration.     

Flow evolution of a turbulent submerged two-dimensional rectangular free jet of air. Average Particle Image Velocimetry (PIV) visualizations and measurements

The paper presents average flow visualizations and measurements, obtained with the Particle Image Velocimetry (PIV) technique, of a submerged rectangular free jet of air in the range of Reynolds numbers from Re = 35,300 to Re = 2200, where the Reynolds number is defined according to the hydraulic diameter of a rectangular slot of height H. According to the literature, just after the exit of the jet there is a zone of flow, called zone of flow establishment, containing the region of mixing fluid, at the border with the stagnant fluid, and the potential core, where velocity on the centerline maintains a value almost equal to the exit one. After this zone is present the zone of established flow or fully developed region. The goal of the paper is to show, with average PIV visualizations and measurements, that, before the zone of flow establishment is present a region of flow, never mentioned by the literature and called undisturbed region of flow, with a length, L_U, which decreases with the increase of the Reynolds number. The main characteristics of the undisturbed region is the fact that the velocity profile maintains almost equal to the exit one, and can also be identified by a constant height of the average PIV visualizations, with length, L_CH, or by a constant turbulence on the centerline, with length L_CT. The average PIV velocity and turbulence measurements are compared to those performed with the Hot Film Anemometry (HFA) technique. The average PIV visualizations show that the region of constant height has a length L_CH which increases from L_CH = H at Re = 35,300 to L_CH = 4–5H at Re = 2200. The PIV measurements on the centerline of the jet show that turbulence remains constant at the level of the exit for a length, L_CT, which increases from L_CT = H at Re = 35,300 to L_CT = 4–5H at Re = 2200. The PIV measurements show that velocity remains constant at the exit level for a length, L_U, which increases from L_U = H at Re = 35,300 to L_U = 6H at Re = 2200 and is called undisturbed region of flow. In turbulent flow the length L_U is almost equal to the lengths of the regions of constant height, L_CH, and constant turbulence, L_CT. In laminar flow, Re = 2200, the length of the undisturbed region of flow, L_U, is greater than the lengths of the regions of constant height and turbulence, L_CT = L_CH = 4–5H. The average PIV and HFA velocity measurements confirm that the length of potential core, L_P, increases from L_P = 4–5H at Re = 35,300 to L_P = 7–8H at Re = 2200, and are compared to the previous experimental and theoretical results of the literature in the zone of mixing fluid and in the fully developed region with a good agreement.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Comprehensive experimental investigation of the hydrodynamics of large-scale, 3D,oscillating bubble plumes

An extensive study of the most important hydrodynamic characteristics of fairly large-scale bubble plumes was conducted using several measurement techniques and a variety of tools to analyze the data. Particle image velocimetry (PIV), double-tip optical probes (OP) and photographic techniques were extensively applied to measure bubble and liquid velocities, void-fraction and bubble sizes. PIV measurements in a vertical plane crossing the centre of the injector provided the instantaneous velocity fields for both phases, as well as hydrodynamic parameters, such as the movement of the axis of the plume and its instantaneous shape. Statistical studies were performed using image processing to determine the distribution of the apparent instantaneous plume diameter and centreline position. An important finding was that there is little instantaneous spreading of the bubble plume core; the spreading of the time-averaged plume width (as measured from the time-averaged void-fraction and time-averaged liquid velocity fields) is largely due to plume meandering and oscillations. The liquid-phase stress tensor distributions obtained from the instantaneous velocity data indicate that, for the continuous phase, these stresses scale linearly with the local void-fraction in the range of 0.5% < α < 2.5%. The bubbles were found to be ellipsoidal, with shape factor e ≈ 0.5.     

Investigation of a continuous adjoint-based optimization procedure for aeroacoustic control of plane jets

A control optimization technique using the continuous adjoint of the compressible Navier–Stokes equations is implemented for aeroacoustic optimization of plane jet flows. The purpose of the adjoint equations is to provide sensitivity information, which is afterwards used in a gradient-based minimization of a prescribed cost functional, designed to describe the far-field sound pressure level (SPL). The objective of the present paper is to demonstrate the ability to reduce the sound in the near far-field of plane jets. Furthermore, as the continuous adjoint approach can become inaccurate, due to inconsistencies between the continuous and the discretized system, the accuracy of the continuous adjoint approach is investigated. The considered cases exhibit a nozzle exit Reynolds number of Re_jet = ρu_jetD/μ = 2000 and a Mach number of M_jet = 0.9, performed using two-dimensional direct numerical simulation and three-dimensional large-eddy simulation, respectively. A comparison of the obtained gradient via adjoint and finite differences is presented and it is shown, that in order to obtain reliable gradient directions, the length of the optimization time needs to be restricted. Furthermore, a receding horizon optimization for the two-dimensional plane jet simulation is used to obtain a sound reduction over much longer time intervals. The influence of different formulations of the viscosity in the adjoint equations is finally investigated.     

Wake structures and vortex-induced vibrations of a long flexible cylinder—Part 1: Dynamic response

Results showing the dynamic response of a vertical long flexible cylinder vibrating at low mode numbers are presented in this paper. The model had an external diameter of 16 mm and a total length of 1.5 m giving an aspect ratio of about 94, with Reynolds numbers between 1200 and 12 000. Only the lower 40% of its length was exposed to the water current in the flume and applied top tensions varied from 15 to 110 N giving fundamental natural frequencies in the range from 3.0 to 7.1 Hz. Reduced velocities based on the fundamental natural frequency up to 16 were reached. The mass ratio was 1.8 and the combined mass–damping parameter about 0.05. Cross-flow and in-line amplitudes, x–y trajectories and phase synchronisation, dominant frequencies and modal amplitudes are reported. Cross-flow amplitudes up to 0.7 diameters and in-line amplitudes over 0.2 were observed with dominant frequencies given by a Strouhal number of 0.16.     

Criteria for the appearance of recirculating and non-stationary regimes behind a cylinder in a viscoplastic fluid

The appearance of a recirculation zone and the formation of non-stationary vortices behind a cylinder in the unconfined flow of a Herschel–Bulkley fluid have been studied by numerical simulation. The Herschel–Bulkley constitutive equation was regularised by using the Papanastasiou model. Special attention was paid to determining the numerical parameters and comparing them to existing results. The influence of the Oldroyd number and power-law index on flow morphology and, in particular, on the unyielded zones was studied over a wide spectrum (0 ≤ Od ≤ 10) and (0.3 ≤ n ≤ 1.8). It was seen that the greater the Oldroyd number, the greater the critical Reynolds numbers and Strouhal number for the two flow regimes. The influence of the power-law index is more complex.