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.     

Integral bounds on the strain energy for the traction problem in finite elasticity

For the traction boundary value problem in finite elasticity, a bound is obtained for the total strain energy in terms of the L₂ integral norms of the surface tractions and body forces, under the assumptions that the unstressed body occupies a convex domain and the displacement gradients are sufficiently small.This is an extension of known results in linear (infinitesimal) elasticity into finite elasticity.     

A bound on the strain energy for the traction problem in finite elasticity with localized non-zero surface data

For the boundary value problem in finite elasticity in which nonzero tractions are given on a connected subdomain of the boundary, the rest of the boundary is stress-free, and there are no body forces, a bound is obtained for the strain energy in terms of the L ₂ integral norm of the surface tractions with the constant involved depending only upon and the material constants.The result is obtained in the context of finite elasticity under the assumptions that the unstressed body occupies a convex domain and the displacement gradients are sufficiently small. In the context of the linear theory, the same result is obtained without these assumptions.

---

Zusammenfassung Wir betrachten ein Randwertproblem in der nichtlinearen Elastizität in dem nur ein zusammenhängendes Teilgebiet der Randfläche belastet ist, sonst aber keine Randbelastung oder Körperkräfte vorhanden sind. Eine Schranke für die Verzerrungsenergie wird mittels der L ₂ Integralnorme der Randbelastung hergeleitet, wobei die auftretende Konstante nur von dem Teilgebiet und von den Eigenschaften des Materials abhängig ist.Das Ergebnis gilt für die nichtlinearen Elastizität, unter den Vorraussetzungen dass das unbelastete Material ein konvexes Gebiet besetzt und dass die Verschiebungsgradiente hinreichend klein sind. Das gleiche Ergebnis gilt in der linearen Elastizität ohne diese Vorraussetzungen.

---

---

The second author was a visitor at Georgia Institute of Technology, School of Mathematics, at the time that the revised version was prepared.     

An integral bound for the strain energy in nonlinear elasticity in terms of the boundary displacements

For the displacement boundary value problem in nonlinear elastostatics with zero body force, an integral bound for the strain energy is obtained in terms of theL ₂-norms of the given boundary displacements and their tangential derivatives (assumed sufficiently small). The constants involved depend upon the strain energy density function and upon the geometry of the domain.     

Heat transfer through partitioned enclosures

A closed-form model for the computation of heat transfer rates through the rectangular-partitioned enclosures is investigated. The rectangular-partitioned enclosures may contain solid or gas with or without a constant and uniformly distributed heat generation. The conduction in the enclosures is considered as two-dimensional, whereas one-dimensional heat transfer through the fin-type partition is assumed. Dimensionless heat flux plots are parametrically studied by varying the aspect ratio (L/H) of the enclosure, the ratio of thermal conductivities of the enclosure to the fin-type partition (k _a /k _f ), the Biot number (β _a =h _a L/k _a ), and the reduced partition thickness (t ^*/L). It is demonstrated through an example problem that there is a large error in using one-dimensional analysis, particularly at lower values of k _a /k _f , and β _a .     

UNSTEADY LIFT AND SOUND PRODUCED BY AN AIRFOIL IN A TURBULENT BOUNDARY LAYER

An analysis is made of the unsteady lift exerted on a stationary rigid body immersed in an incompressible, plane-wall turbulent boundary layer. The lift is expressed as a surface integral over the body involving theupwash velocity induced by the “free” vorticity Ω (found by taking explicit account of the interaction of the body with the flow and excluding the bound vorticity) and a harmonic function X₂that depends only on the shape of the body. The upwash velocity is the free-field velocity given in terms of Ω by the Biot–Savart formula, augmented by the velocity field of a conventional distribution of image vortices in the wall. The function X₂can be interpreted as the velocity potential of flow past the body, produced by motion of the wall at unit speed towards the body. Detailed predictions are made of the lift on a slender airfoil placed in the outer region of the boundary-layer. When the airfoil chord is large compared to the boundary-layer thickness, vortex shedding into the wake causes the magnitude of the net upwash velocity near the trailing edge to be small. The main contributions to the surface integral are then from the nose region, where the upwash velocity may be estimated independently of the fluctuations near the trailing edge. Analytical results for a thin plate airfoil of chord 2a at distance h from the wall show that the lift increases as a/h increases; it is ultimately independent of a and scales with the ratio of h to the hydrodynamic wavelength. Application is made to determine the sound generated by the airfoil in a weakly compressible boundary layer flow over a finite elastic plate.     

The response of strain gages to longitudinally sweeping strain pulses

The current state-of-the-art for estimating the maximum rise time of a strain gage which is subjected to an axially sweeping strain pulse ist _rg≤0.8L/c+0.5μsec wheret _rg is the 10/90 rise time of the strain gage,^* L is the gage length andc is the strain-pulse velocity. This paper shows that the effect of the 0.8L/c term can be significantly reduced by utilizing an analytical compensation technique. In addition, it is shown that the 0.5 μsec additive constant can be reduced to approximately 0.1 μsec by reinterpreting data published by a previous author.     

Prediction of micro-bubble dissolution characteristics in water and seawater

This paper is concerned with the prediction of micro-bubble dissolution characteristics in water and seawater when microbubbles are generated by a Sadatomi-type micro-bubble generator (2003) with a spherical body in a flowing liquid tube. In the experiments, in order to know the effects of the salinity on the characteristics, tap water and an artificial seawater with different salt concentrations of 1 and 3 wt% were used as the test liquids. Parameters measured were the Sauter mean diameter of bubbles, d_BS, the void fraction, α, the rising velocity of bubbles, u_G, the interfacial area concentration, a, the volumetric mass transfer coefficient, K_La, and the liquid-side mass transfer coefficient, K_L. In the analysis, for predicting α, K_La and K_L, some correlations in the literatures were tested against the present data. Furthermore, in order to improve the predictability, new correlations were developed based on the present data. The prediction of K_La with the new correlation agreed well with Nishino et al.’s [T. Nishino, K. Terasaka, M. Ishida, Application for several micro-bubble generators for gas absorber, in: Proceedings of the Annual Meeting of the Japanese Society for Multiphase Flow, 2006, pp. 276–277 (in Japanese)] and Li and Tsuge’s [P. Li, H. Tsuge, Water treatment by induced air flotation using microbubbles, Journal of Chemical Engineering of Japan 39 (2006) 896–903; P. Li, H. Tsuge, Ozone transfer in a new gas-induced contactor with microbubbles, Journal of Chemical Engineering of Japan 39 (2006) 1213–1220] data for different aeration systems using several different micro-bubble generators.     

Turbulence measurements in the bubbly flow region of hydraulic jumps

A hydraulic jump is characterized by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air–water flow measurements were performed in a large-size facility using two types of phase-detection intrusive probes: i.e. single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent integral time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of Chanson [H. Chanson, Bubbly flow structure in hydraulic jump, Eur. J. Mech. B/Fluids 26 (3) (2007) 367–384], providing the turbulent scales of the eddy structures advecting the air bubbles in the developing shear layer. The length scale L_xz is an integral air–water turbulence length scale which characterized the transverse size of the large vortical structures advecting the air bubbles. The experimental data showed that the dimensionless integral turbulent length scale L_xz/d₁ was closely related to the inflow depth: i.e. L_xz/d₁ = 0.2–0.8, with L_xz increasing towards the free-surface.     

Integral representations of constitutive equations

This paper is concerned with constitutive equations in which the stress at an instant of time,t, say, is a functional of the deformation gradient history up to and includingt. It is assumed that the deformation gradient histories are of bounded variation and are continuous, or piece-wise continuous with a countable number of salti. It is shown that if the stress functional isn times Fréchet differentiable in the sense of the supremum norm, it can be represented by a series of multiple integrals in the deformation gradient history, with an error or ordern + 1 in the displacement gradients. The relation of this type of integral representation to that of Green and Rivlin is discussed.     

Direct Evaluation of Cohesive Law in Mode I of Pinus pinaster by Digital Image Correlation

The direct identification of the cohesive law in pure mode I of Pinus pinaster is addressed in this work. The approach couples the double cantilever beam (DCB) test with digital image correlation (DIC). Wooden beam specimens loaded in the radial-longitudinal (RL) fracture propagation system are used. The strain energy release rate in mode I (G _I) is uniquely determined from the load–displacement curve by means of the compliance-based beam method (CBBM). This method relies on the concept of equivalent elastic crack length (a _eq) and therefore does not require the monitoring of crack propagation during test. DIC measurements are processed with two different purposes. Firstly, the physical evidence of a _eq is discussed with regard to actual estimation of the crack length based on post-processing full-field displacement measurements. Secondly, the crack tip opening displacement in mode I (w _I) is determined from the displacements near the initial crack tip. The cohesive law in mode I (σ _I???w _I) is then identified by numerical differentiation of the G _I???w _I relationship. The methodology and accuracy on this reconstruction are addressed. Moreover, the proposed procedure is validated by finite element analyses including cohesive zone modelling. It is concluded that the proposed data reduction scheme is adequate for assessing the cohesive law in pure mode I of P. pinaster.     

Size-dependent rigidities of nanosized torsional elements

A theory for the prediction of the size dependence of torsional rigidities of nanosized structural elements is developed. It is shown that, to a very good approximation, the torsional rigidity (D) of a nanosized bar differs from the prediction of standard continuum mechanics (D_c) as (D−D_c)/D_c=Ah₀/a where A is a non-dimensional constant, a is the size scale of the cross-section of the bar and h₀ is a material length equal to the ratio of the surface elastic constant to the bulk elastic constant. The theory developed is compared with direct atomistic calculations (“numerical experiment”) of the torsional rigidity bars made of several FCC metals modeled using the embedded atom method. Very good agreement is obtained between theory and simulation. The framework presented here can aid the development of design methodologies for nanoscale structural elements without the need for full scale atomistic simulations.     

Pressure potential formulation of 2-D stokes problems in multiply-connected domains

In general Stokes problems, no boundary conditions exist for the pressure. But pressure is an L²(Ω) function and can uniquely be represented as the divergence of a precisely defined vector field. In the 2-D case, this vector field can in turn be represented as the sum of a gradient (of a pressure-potential) and the curl of a second scalar potential. The latter potential is entirely determined by the first one. A variational equation is obtained for such pressure potential class, which exists and is uniquely characterized. This variational problem is well-posed. Finite element approximations can easily be realized and ensure high convergence rates for the L²(Ω) norm of the pressure.     

New type of equation for the relation between viscosity and particle content in suspensions

Many disperse systems show a typical non-Newtonian flow at relatively high concentrations of the disperse particle. However, two Newtonian viscosities η_∞ and η₀ can be, respectively, determined at high and low rates of shear. Expect for very low particle content, η_∞/η_s is proportional to exp(mϕ), where η_s is a medium viscosity, m a constant which might reflect the particle-particle interaction and ϕ the volume fraction. In considering this relationship, a new type of equation which describes the relation between the zero shear relative viscosity η_r0( ≡ η₀/η_s) of the disperse system and ϕ is proposed as follows. ln η_r0 = A(p)ϕ + am³ϕ², where A(p), the Einstein-Simha constant, is a function of the axial ratio p of dispersing particles, and a is a constant (⋍ 0.03) which depends slightly on the particle shape.The equation has been compared with the experimental results obtained for several disperse systems. A number of disperse systems of spherical particles are described well by the choice A(p) = 2.5 and a = 0.027, and a system of rod-like particles with p = 50 by the choice A(p) = 215.6 and a = 0.033. m for rod-like particles is larger than that for spherical particles.     

Optimizing strength and toughness of nanofiber-reinforced CMCs

Nanofibers used in current ceramic matrix composites (CMCs) are usually wavy and of finite length. Here, a shear-lag model for predicting the tensile strength and work of fracture of a composite containing a “single matrix crack”, as a function of all the material parameters including constant plus Coulomb interfacial friction, is presented for a CMC containing wavy, finite-length nanofibers having a statistical distribution of strengths. Literature results are recovered for straight infinite fibers, with strength and toughness depending on a characteristic strength σ_c and a characteristic length δ_c. For nanofibers of finite length L, radius of curvature R, and with Coulomb friction coefficient μ, the strength and toughness are found to depend only on σ_c, δ_c, and two new dimensionless parameters μδ_c/R and L/δ_c. Parametric results for a typical carbon nanotube CMC show an optimal region of morphology (L and R) that maximizes composite strength and toughness, exceeding the properties of composites with straight (R=∞) and/or long (L=∞) reinforcements. Therefore, while factors such as the nanofiber strength distribution and the nanofiber–matrix interface sliding resistance may not be easily controlled, the tuning, via processing, of purely geometrical properties of the nanofibers (L and R) alone can aid in maximizing composite properties. These results have important application in hybrid CMCs such as “fuzzy fiber” CMCs, where micron-scale fibers are covered with a forest of nanofibers such that the nanofiber array can bridge longitudinal cracks and thus improve delamination resistance.     

A stress analytical solution for Mode III crack within modified gradient elasticity

The strain gradient exists near a crack tip may significantly influence the near-tip stress field. In this paper, the strain gradient and the internal length scales are introduced into the basic equations of mode III crack by the modified gradient elasticity (MGE). By using a complex function approach, the analytical solution of stress fields for mode III crack problem is derived within MGE. When the internal length scales vanish, the stress fields can be simplified to the stress fields of classical linear elastic fracture mechanics. The results show that the singularity of the shear stress is made up of two parts, r^−1/2 part and r^−3/2 part, and the sign of the stress σ_yz changes. With the increase of l_x, the peak value of σ_yz decrease and its location moves farther from the fracture vertex. The influence of strain gradient for mode III crack problem cannot be ignored.     

Numerical study on coalescence of two pre-existing coplanar flaws in rock

Crack propagation and coalescence processes are the fundamental mechanisms leading to progressive failure processes in rock masses, in which parallel non-persistent rock joints are commonly involved. The coalescence behavior of the latter, which are represented as pre-existing coplanar flaws (cracks), is numerically investigated in the present study. By using AUTODYN as the numerical tool, the present study systematically simulates the coalescence of two pre-existing coplanar flaws in rock under compression. The cumulative damage failure criterion is adopted in the numerical models to simulate the cumulative damage process in the crack initiation and propagation. The crack types (shear or tensile) are identified by analyzing the mechanics information associated with the crack initiation and propagation processes. The simulation results, which are generally in a good accordance with physical experimental results, indicate that the ligament length and the flaw inclination angle have a great influence on the coalescence pattern. The coalescence pattern is relatively simple for the flaw arrangements with a short ligament length, which becomes more complicated for those with a long ligament length. The coalescence trajectory is composed of shear cracks only when the flaw inclination angle is small (such as β ? 30°). When the pre-existing flaws are steep (such as β ? 75°), the coalescence trajectory is composed of tensile cracks as well as shear cracks. When the inclination angle is close to the failure angle of the corresponding intact rock material, and the ligament length is not long (such as L ? 2a), the direct shear coalescence is the more favorable coalescence pattern. In the special case that the two pre-existing flaws are vertical, the model will have a direct tensile coalescence pattern when the ligament length is short (L ? a), while the coalescence between the two inner flaw tips is not easy to achieve if the ligament length is long (L ? 2a).     

Transient Behavior of Turbulent Scalar Transport in Premixed Flames

Transition from gradient to countergradient scalar transport in a statistically planar, one-dimensional, developing, premixed turbulent flame is studied both theoretically and numerically. A simple criterion of the transition referred to is derived from the balance equation for the combustion progress variable, with the criterion highlighting an important role played by flame development. A balance equation for the difference in velocities $\bar{u}_b$ and $\bar{u}_u$ conditioned on burned and unburned mixture, respectively, is numerically integrated. Both analytical and computed results show that; (1) The flux $\overline{\rho u'' c''}$ is gradient during an early stage of flame development followed by transition to countergradient scalar transport at certain instant t _tr . (2) The transition time is increased when turbulence length scale L is increased or when the laminar flame speed S _L and/or the density ratio are decreased. (3) The transition time normalized using the turbulence time scale is increased by u??. Moreover, the numerical simulations have shown that the transition time is increased by u?? if a ratio of u??/S _L is not large. This dependence of t _tr on u?? is substantially affected by (i) the mean pressure gradient induced within the flame due to heat release and (ii) by the damping effect of combustion on the growth rate of mean flame brush thickness. The reasonable qualitative agreement between the computed trends and available experimental and DNS data, as well as the agreement between the computed trends and the present theoretical results, lends further support to the conditioned balance equation used in the present work.     

Lift force acting on single bubbles in linear shear flows

Lift coefficients, C_L, of single bubbles in linear shear flows are measured to investigate effects of the bubble shape, the liquid velocity gradient and the fluid property on C_L. The range of the Morton number, M, tested is from logM = − 6.6 to − 3.2. The shapes of bubbles are spherical and ellipsoidal. A correlation of bubble aspect ratio for single bubbles in infinite stagnant liquids proposed in our previous study can give good evaluations for bubbles in the linear shear flows. The C_L of spherical bubbles at low bubble Reynolds numbers, Re, depend on the dimensionless shear rate Sr and Re and decrease with increasing Re. These characteristics agree with the Legendre-Magnaudet correlation. The use of a single dimensionless group such as Re, the Eötvös number, the Weber number and the Capillary number cannot correlate C_L of non-spherical bubbles. The trend of the critical Re for the reversal of the sign of C_L is the same as that for the onset of oscillation of bubble motion, which supports the mechanism proposed by Adoua et al., at least within the range of −6.6 ≤ logM ≤ −3.2. An experimental database of C_L is provided for validation of available C_L models and CFD.     

POD analysis of a finite-length cylinder near wake

The near wake of a wall-mounted finite-length square cylinder with an aspect ratio of 7 is investigated based on the proper orthogonal decomposition (POD) of the PIV data measured in three spanwise planes, i.e., z/d = 6, 3.5 and 1.0, near the cylinder free end, mid-span and fixed end (wall), respectively. The Reynolds number based on free-stream velocity (U _∞) and cylinder width (d) is 9,300. A two-dimensional (2D) square cylinder wake is also measured and analyzed at the same Reynolds number for the purpose of comparison. The structures of various POD modes show marked differences between the two flows. While the coefficients, a ₁ and a ₂, of the POD modes 1 and 2 occur within an annular area centered at a ₁ = a ₂ = 0 in the 2D wake, their counterparts are scattered all over the entire circular plane at z/d = 1.0 and 3.5 of the finite-length cylinder wake. Flow at z/d = 6 is dominated by POD mode 1, which corresponds to symmetrical vortex shedding and accounts for 54.0 % of the total turbulent kinetic energy (TKE). On the other hand, the POD modes 1 and 2, corresponding to anti-symmetrical vortex shedding, are predominant, accounting for about 45.0 % of the total TKE, at z/d = 3.5 and 1. It has been found that the flow structure may be qualitatively and quantitatively characterized by the POD coefficients. For example, at z/d = 6, a larger a ₁ corresponds to a smaller length of flow reversal zone and a stronger downwash flow. At z/d = 3.5 and 1, two typical flow modes can be identified from a ₁ and a ₂. While large a ₁ and/or a ₂ correspond to anti-symmetrical vortex shedding, as in a 2D cylinder wake, small a ₁ and a ₂ lead to symmetrical vortex shedding. Any values between the large and small a ₁ and/or a ₂ correspond to the flow structure between these two typical flow modes. As such, the probability of occurrence of a flow structure may be determined from the distribution of the POD coefficients.