Near-wall flow in a three-dimensional boundary layer on the endwall of a 30° bend

 Turbulence measurements are reported on the three-dimensional turbulent boundary layer along the centerline of the flat endwall in a 30° bend. Profiles of mean velocities and Reynolds stresses were obtained down to y ⁺≈2 for the mean flow and y ⁺≈8 for the turbulent stresses. Mean velocity data collapsed well on a simple law-of-the-wall based on the magnitude of the resultant velocity. The turbulence intensity and turbulent shear stress magnitude both increased with increased three-dimensionality. The ratio of these two quantities, the a ₁ structure parameter, decreased in the central regions of the boundary layer and showed profile similarity for y ⁺<50. The shear stress vector angle lagged behind the velocity gradient vector angle in the outer region of the boundary layer, however there was an indication that the shear stress vector tends to lead the velocity gradient vector close to the wall. Received: 16 July 1996/Accepted: 14 July 1997     

Ignition of liquid and dust fuel layers by gaseous detonation

Abstract. Ignition of a liquid layer and dust fuel layer by a detonation wave propagating in hydrogen-oxygen and acetylene-oxygen mixtures is reported. Experiments were carried out using a shock tube equipped with optical-quality observation windows. A schlieren system and a high-speed camera were used for measurements of ignition delay. Pressure transducers provided data necessary for measurements of the detonation wave velocity and pressure variation within the front of the interacted detonation wave and fuel layer. Kerosene, nitroglycerin and PETN were used as fuels. Investigation shows that the layer of liquid fuel can be efficiently ignited by detonation wave. It was found that the ignition delay of the fuel layer depends mostly on the detonation wave velocity and sensitivity of igniting fuels, and slightly on the layer thickness. Received 12 August 2001 / Accepted 1 July 2002 Published online 4 February 2003 Correspondence to: P. Wolanski (e-mail: wolanski@itc.pw.edu.pl) An abridged version of this paper was presented at the 18th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Seattle, USA, from July 29 to August 3, 2001     

Global Equivalence for Deformable Theormoeleastic Bodies

We consider equilibria arising in a model for phase transitions which correspond to stable critical points of the constrained variational problem Here W is a double‐well potential and is a strictly convex domain. For ε small, this is closely related to the problem of partitioning Ω into two subdomains of fixed volume, where the subdomain boundaries correspond to the transitional boundary between phases. Motivated by this geometry problem, we show that in a strictly convex domain, stable critical points of the original variational problem have a connected, thin transition layer separating the two phases. This relates to work in [GM] where special geometries such as cylindrical domains were treated, and is analogous to the results in [CHo] which show that in a convex domain, stable critical points of the corresponding unconstrained problem are constant. The proof of connectivity employs tools from geometric measure theory including the co‐area formula and the isoperimetric inequality on manifolds. The thinness of the transition layer follows from a separate calculation establishing spatial decay of solutions to the pure phases. (Accepted July 15, 1996)     

The axisymmetric boundary layer beneath a Rankine-like vortex

 An experimental investigation of the three-dimensional boundary layer induced by a Rankine-like vortex with its axis normal to a stationary disk is described. The velocity field through the boundary layer was measured for Reynolds number Re (based on the tangential velocity and radius at the disk edge) ranging from 10 000 to 25 000 at various radial distances by means of a 4-beam, 2-component Laser Doppler Anemometer. Our results show that the nature of the boundary layer is affected by two factors: an inflexional instability caused by the crossflow velocity profile and a stability factor caused by the favorable pressure gradient. At lower Reynolds number, the radial pressure gradient has a very strong stabilizing effect on the boundary layer and acts to revert it to its laminar state upstream of the effusing core. At higher Re the inflexional instability caused by the crossflow velocity dominates while the stabilizing influence of the favorable pressure gradient recedes. As such, laminar reversion likely occurs closer to the effusion core. Thus, the point of laminar reversion moves closer to the effusion core as the Reynolds number is increased. Received 23 May 1996 / Accepted 29 July 1996     

Influence of roughness on attachment line boundary layer transition in hypersonic flow

 Attachment line boundary layer transition on swept cylinders is studied in a low enthalpy hypersonic wind tunnel at M _∞=7.14. Sweep angles of 60° and 70° are used and transition is detected by means of heat flux measurements. The influence on attachment line transition of single 2D-roughness elements, in the form of tripwires or slots, as well as 3D obstacles is determined and the results are analyzed with respect to Poll’s criterion. Received: 16 January 1996 / Accepted: 12 July 1996     

Convective--Diffusive Transport with Chemical Reaction in Natural Convection Flows

A study of laminar natural convection flow over a semi-infinite vertical plate at constant species concentration is examined. The plate is maintained at a given concentration of some chemical species while convection is induced by diffusion into and chemical reaction with the ambient fluid. In the absence of chemical reaction, a similarity transform is possible. When chemical reaction occurs, perturbation expansions about an additional similarity variable dependent on reaction rate must be employed. Two fundamental parameters of the problem are the Schmidt number, Sc, and the reaction order, n. Results are presented for the Schmidt number ranging from 0.01 to 10000 and reaction order up to 5. In the presence of a chemical reaction, the diffusion and velocity domains expand out from the plate. This results in a larger, less distinct convection layer. Received 21 July 1998 and accepted 24 June 1999     

Comparison of kinetic and continuum approaches for simulation of shock wave/boundary layer interaction

The present paper, which is a collaboration between three different research groups, analyzes the efficiency of various numerical approaches to describe the complex problem of shock wave/boundary layer interaction. Computations were carried out based on a kinetic approach (Direct Simulation Monte Carlo method) and on two continuum approaches (Navier-Stokes equations and quasigasdynamic equations), which are validated by comparison with experimental results obtained in the R5Ch blowdown Hypersonic Wind Tunnel in ONERA. The influence of the slip boundary conditions for two continuum approaches are also studied. The results obtained by all models display the good prediction of the main structure of the flow and the levels of the flux coefficients are very close to those measured. The implementation of the slip boundary condition for the continuum approaches improves the agreement with the experimental data. Received 12 July 2001 / Accepted 24 May 2002 /Published online 4 December 2002 Correspondence to: D. Zeitoun (e-mail: David.Zeitoun@polytech.univ-mrs.fr) An abridged version of this paper was presented at the 23rd Int. Symposium on Shock Waves at Fort Worth, Texas, from July 22 to 27, 2001     

Electrical measurement of sediment layer thickness under suspension flows

 This paper reports a new technique to measure the thickness of a layer of deposited sediment as a function of time, independent of the flow conditions or presence of suspended sediment above the layer. Small electrodes on the bottom and a reference electrode in the fluid above were used to measure the resistance of the layer with a small AC current and a bridge circuit. Using a multiplexer and an Analog-to-Digital converter the growth of the layer can be accurately monitored at many locations on the tank bottom. In a trial experiment the sedimentation under a stagnant column of a monodisperse suspension was examined. The results show that changes in the sediment layer thickness of less than 0.3% can be measured for layers up to 0.2 g/cm². Received: 8 February 1998/Accepted: 19 July 1998     

Laminated shells with debonding between laminas in temperature field

A heat-conduction problem is formulated for laminated plates and shells with a heat-conducting layer and debonding between laminas. The approach consists in analyzing how the layer thickness changes in the process of debonding of laminas and deformation of plates and shells. The three-dimensional thermoelastic and heat-conduction equations are expanded into polynomial Legendre series in thickness. The first-order, Timoshenko’s, and Kirchhoff-Love equations are examined. A numerical example of laminated shells with a heat-conducting layer is considered Published in Prikladnaya Mekhanika, Vol. 42, No. 7, pp. 135–141, July 2006.     

Effects of Mach number on turbulent separation behaviours induced by blunt fin

 An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M _∞=7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin. Received: 21 July 1996/Accepted: 4 February 1998     

Unsteady free convection on a vertical cylinder with variable heat and mass flux

The unsteady natural convection boundary layer flow over a semi-infinite vertical cylinder is considered with combined buoyancy force effects, for the situation in which the surface temperature T ^′ _w(x) and C ^′ _w(x) are subjected to the power-law surface heat and mass flux as K(T ^′/r) = −ax ⁿ and D(C ^′/r) = −bx ^m . The governing equations are solved by an implicit finite difference scheme of Crank-Nicolson method. Numerical results are obtained for different values of Prandtl number, Schmidt number ‘n’ and ‘m’. The velocity, temperature and concentration profiles, local and average skin-friction, Nusselt and Sherwood numbers are shown graphically. The local Nusselt and Sherwood number of the present study are compared with the available result and a good agreement is found to exist. Received on 7 July 1998     

Flow past a cylinder: shear layer instability and drag crisis

Flow past a circular cylinder for Re=100 to 10⁷ is studied numerically by solving the unsteady incompressible two‐dimensional Navier–Stokes equations via a stabilized finite element formulation. It is well known that beyond Re ～ 200 the flow develops significant three‐dimensional features. Therefore, two‐dimensional computations are expected to fall well short of predicting the flow accurately at high Re. It is fairly well accepted that the shear layer instability is primarily a two‐dimensional phenomenon. The frequency of the shear layer vortices, from the present computations, agree quite well with the Re^0.67 variation observed by other researchers from experimental measurements. The main objective of this paper is to investigate a possible relationship between the drag crisis (sudden loss of drag at Re ～ 2 × 10⁵) and the instability of the separated shear layer. As Re is increased the transition point of shear layer, beyond which it is unstable, moves upstream. At the critical Reynolds number the transition point is located very close to the point of flow separation. As a result, the shear layer eddies cause mixing of the flow in the boundary layer. This energizes the boundary layer and leads to its reattachment. The delay in flow separation is associated with narrowing of wake, increase in Reynolds shear stress near the shoulder of the cylinder and a significant reduction in the drag and base suction coefficients. The spatial and temporal power spectra for the kinetic energy of the Re=10⁶ flow are computed. As in two‐dimensional isotropic turbulence, E(k) varies as k^?5/3 for wavenumbers higher than energy injection scale and as k^?3 for lower wavenumbers. The present computations suggest that the shear layer vortices play a major role in the transition of boundary layer from laminar to turbulent state. Copyright © 2004 John Wiley & Sons, Ltd.     

Investigation of Anisotropy-Resolving Turbulence Models by Reference to Highly-Resolved LES Data for Separated Flow

The predictive properties of several non-linear eddy-viscosity models are investigated by reference to highly-resolved LES data obtained by the authors for an internal flow featuring massive separation from a curved surface. The test geometry is a periodic segment of a channel constricted by two-dimensional (2D) `hills' on the lower wall. The mean-flow Reynolds number is 21560. Periodic boundary conditions are applied in the streamwise and spanwise directions. This makes the statistical properties of the simulated flow genuinely 2D and independent from boundary conditions, except at the walls. The simulation was performed on a high-quality, 5M-node grid. The focus of the study is on the exploitation of the LES data for the mean-flow, Reynolds stresses and macro-length-scale. Model solutions are first compared with the LES data, and selected models are then subjected to a-priori studies designed to elucidate the role of specific model fragments in the non-linear stress-strain/vorticity relation and their contribution to observed defects in the mean-flow and turbulence fields. The role of the equation governing the length-scale, via different surrogate variables, is also investigated. It is shown that, while most non-linear models overestimate the separation region, due mainly to model defects that result in insufficient shear stress in the separated shear layer, model forms can be derived which provide a satisfactory representation of the flow. One such model is identified. This combines a particular quadratic constitutive relation with a wall-anisotropy term, a high-normal-strain correction and a new form of the equation for the specific dissipation ω = ∈/k. This revised version was published online in July 2006 with corrections to the Cover Date.     

Errors in hot-wire X-probe measurements induced by unsteady velocity gradients

 Errors in hot-wire X-probe measurements due to unsteady velocity gradients are investigated by a comparison of hot-wire and laser Doppler velocimetry (LDV) measurements. The studied flow case is a laminar boundary layer subjected to high levels of free-stream turbulence, and the hot-wire data shows a local maximum in the wall-normal fluctuation velocity inside the boundary layer. The observed maximum is in agreement with existing hot-wire data, but in conflict with the present LDV measurements as well as existing results from numerical simulations. An explanation for the measurement error is suggested in the paper. Received: 16 October 2000 / Accepted: 23 July 2001     

From Two-Point Closures of Isotropic Turbulence to LES of Shear Flows

We first recall the EDQNM two-point closure approach of three-dimensional isotropic turbulence. It allows in particular prediction of the infrared kinetic-energy dynamics (with ak ₄ backscatter) and the associated time-decay law of kinetic-energy, useful in particular for one-point closure modelling. Afterwards, we show how the spectral eddy viscosity concept may be used for large-eddy simulations: we introduce the plateau-peak model and the spectral-dynamic models. They are applied to decaying isotropic turbulence, and allow recovery of the EDQNM infrared energy dynamics. Anew infrared k ₂ law for the pressure spectrum, predicted by the closure, is also well verified. Assuming that subgrid scales are not too far from isotropy, the spectral-dynamic model is applied to the channel flow at h ₊= 390, with statistics in very good agreement with DNS, while reducing considerably the computational time. We study with the aid of DNS and LES the case of the channel rotating about an axis of spanwise direction. The calculations allow to recover the universal linear behaviour of the mean velocity profile, with a local Rossby number equal to −1. We present also LES (using the Grenoble Filtered Structure-Function Model), of a turbulent boundary layer passing over a cavity. Finally, we make some remarks on the future of LES for industrial applications. This revised version was published online in July 2006 with corrections to the Cover Date.     

An internal crack problem for an infinite elastic layer

The elastostatic plane problem of an infinite elastic layer with an internal crack is considered. The elastic layer is subjected to two different loadings, (a) the elastic layer is loaded by a symmetric transverse pair of compressive concentrated forces P/2, (b) it is loaded by a symmetric transverse pair of tensile concentrated forces P/2. The crack is opened by an uniform internal pressure p ₀ along its surface and located halfway between and parallel to the surfaces of the elastic layer. It is assumed that the effect of the gravity force is neglected. Using an appropriate integral transform technique, the mixed boundary value problem is reduced to a singular integral equation. The singular integral equation is solved numerically by making use of an appropriate Gauss–Chebyshev integration formula and the stress-intensity factors and the crack opening displacements are determined according to two different loading cases for various dimensionless quantities.     

Anisotropy measurements in the boundary layer over a flat plate with suction

Laser Doppler velocity measurements are carried out in a turbulent boundary layer subjected to concentrated wall suction (through a porous strip). The measurements are taken over a longitudinal distance of 9× the incoming boundary layer thickness ahead of the suction strip. The mean and rms velocity profiles are affected substantially by suction. Two-point measurements show that the streamwise and wall-normal autocorrelations of the streamwise velocity are reduced by suction. It is found that suction alters the redistribution of the turbulent kinetic energy k between its components. Relative to the no-suction case, the longitudinal Reynolds stress contributes more to k than the other two normal Reynolds stresses; in the outer region, its contribution is reduced which suggests structural changes in the boundary layer. This is observed in the anisotropy of the Reynolds stresses, which depart from the non-disturbed boundary layer. With suction, the anisotropy level in the near-wall region appears to be stronger than that of the undisturbed layer. It is argued that the mean shear induced by suction on the flow is responsible for the alteration of the anisotropy. The variation of the anisotropy of the layer will make the development of a turbulence model quite difficult for the flow behind suction. In that respect, a turbulence model will need to reproduce well the effects of suction on the boundary layer, if the model is to capture the effect of suction on the anisotropy of the Reynolds stresses.     

COUPLING BETWEEN LAMINAR FILM CONDENSATION AND NATURAL CONVECTION ON OPPOSITE SIDES OF A VERTICAL PLATE

Theoretical analyses which incorporate one-dimensional heat conduction along a plate and transverse heat conduction approximations are presented to predict the net heat transfer between laminar film condensation of a saturated vapour on one side of a vertical plate and boundary layer natural convection on the other side. It is assumed that countercurrent boundary layer flows are formed on the two sides. The governing boundary layer equations of this problem and their corresponding boundary conditions are all cast into dimensionless forms by using a non-similarity transformation. Thus the resulting system of equations can be solved by using the local non-similarity method for the boundary layer equations and a finite difference method for the heat conduction equation of the plate. The plate temperature and the heat flux through the plate are repetitively determined until the solutions for each side of the plate match. The predicted results show that the effect of Pr_c is not negligible for larger values of A* (thermal resistance ratio between natural convecti on side and condensing film side) and the approximation of transverse heat conduction overpredicts the plate temperature for lower values of R_t (thermal resistance ratio between plate and condensing film). However, no significant differences are observed between the two different approximations for higher values of R_t. © by 1997 John Wiley & Sons, Ltd.     

Reacting, Circular Mixing Layers in Transition to Turbulence

The evolution of a reacting, circular mixing layer - a model of round-jet flow - in its transition to turbulence was studied by direct numerical simulation. An economical Fourier pseudospectral method was combined with the third-order Adams-Bashforth scheme to integrate Navier-Stokes and scalar transport equations. The Reynolds number based on initial mixing-layer diameter and velocity difference was 1600. The initially thin mixing layer encloses a cylindrical core of fuel that mixes and reacts with the surrounding oxidizer. Both fast and finite-rate reactions were examined. The stages in transition are characterized by roll-up of the mixing layer into a sequence of vortex rings, pairing of adjacent rings, azimuthal instability, and breakdown to a disordered (turbulent) state. Reaction surfaces in the fast reaction limit become extended, folded and pinched off at various times corresponding to the dynamics of the vortices observed in the simulations. When the equivalence ratio is O(1) or smaller,the progress of reaction is determined by the dynamics of vortex rings. For larger ratios there is a qualitative difference: Initially, the flame is located well outside the rings and is relatively unaffected. Following breakdown to turbulence, there is a steep increase in flame surface area resulting in a noticeable change in fuel consumption rate. At smaller reaction rates (small Damkohler numbers), the reaction zones are diffuse and fill the vortical (mixed) regions. Product accumulates in and its presence raises the temperature of vortex cores, but reaction rates remain low due to low reactant concentrations. Reaction rates are highest in the braids between vortex rings where scalar dissipation rates and compressive strain rates show the highest values. This revised version was published online in July 2006 with corrections to the Cover Date.