Large injection rates on unsteady compressible boundary layer flow over cones at angles of attack

Summary Three-dimensional unsteady laminar boundary layer near the planes of symmetry of sharp cones at angles of attack subject to large rates of injection is obtained numerically by using an implicit finite difference scheme in combination with the quasi-linearization technique. Several model gases are considered with Mach numbers, wall-to-total-enthalpy ratios, and cross-flow parameters spanning the ranges of main engineering interest. A detailed study has been made of the solutions in the symmetry plane, in order to increase the understanding of the problem. Various cases are considered, when the free-stream velocity and the surface mass transfer (injection) vary arbitrarily with time. The effects of viscous dissipation and the cross-flow parameter have also been discussed.This research has been partially supported by the Research and Development Centre for Iron and Steel, Steel Authority of India Ltd. The constructive comments of Professor G. Nath and Professor A. K. Lahiri are sincerely appreciated.     

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

Inviscid,laminar and turbulent opposed flows

This paper attempts to reproduce numerically previous experimental findings with opposed flows and extends their range to quantify the effects of upstream pipes and nozzles with inviscid, laminar and turbulent flows. The choice of conservation equations, boundary conditions, algorithms for their solution, the degree of grid dependence, numerical diffusion and the validity of numerical approximations are justified with supporting calculations where necessary. The results of all calculations on the stagnation plane show maximum strain rates close to the annular exit from the nozzles and pipes for lower separations and it can be expected that corresponding reacting flows will tend to extinguish in this region with the extinction moving towards the axis. With laminar flows, the maximum strain rate increased with Reynolds number and the maximum values were generally greater than with inviscid flows and smaller than with turbulent flows. With large separations, the strain rates varied less and this explains some results with reacting flows where the extinction appeared to begin on the axis. The turbulent‐flow calculations allowed comparison of three common variants of a two‐equation first‐moment closure. They provided reasonable and useful indications of strain rates but none correctly represented the rms of velocity fluctuations on the axis and close to the stagnation plane. As expected, those designed to deal with this problem produced results in better agreement with experiment but were still imperfect. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical prediction of turbulent flow structure generated by a synthetic cross‐jet into a turbulent boundary layer

A detailed numerical study using large‐eddy simulation (LES) and unsteady Reynolds‐averaged Navier–Stokes (URANS) was undertaken to investigate physical processes that are engendered in the injection of a circular synthetic (zero‐net mass flux) jet in a zero pressure gradient turbulent boundary layer. A complementary study was carried out and was verified by comparisons with the available experimental data that were obtained at corresponding conditions with the aim of achieving an improved understanding of fluid dynamics of the studied processes. The computations were conducted by OpenFOAM C++, and the physical realism of the incoming turbulent boundary layer was secured by employing random field generation algorithm. The cavity was computed with a sinusoidal transpiration boundary condition on its floor. The results from URANS computation and LES were compared and described qualitatively and quantitatively. There is a particular interest for acquiring the turbulent structures from the present numerical data. The numerical methods can capture vortical structures including a hairpin (primary) vortex and secondary structures. However, the present computations confirmed that URANS and LES are capable of predicting current flow field with a more detailed structure presented by LES data as expected. Copyright © 2011 John Wiley & Sons, Ltd.     

Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Nature of the surface heart transfer fluctuation in a hypersonic separated turbulent flow

This paper presents the results of an experimental study of the unsteady nature of a hypersonic separated turbulent flow. The nomimal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5×10⁷/m. The separated flow was generated using finite span forward facing steps. An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make multi-channel measurements of the fluctuating surface heat trtansfer within the separated flow. Conditional sampling analysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness. The compression waves converge into a single leading shock beyond the boundary layer. The shock structure is unsteady and undergoes large-scale motion in the streamwise direction. The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave. There exists a wide band of frequency of oscillations of the shock system. Most of the frequencies are in the range of 1–3 kHz. The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave. This intermittent phenomenon is considered as the consequence of the large-scale shock system oscillations. Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed. The project supported by National Natural Science Foundation of China     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical simulations of non‐equilibrium turbulent boundary layer flowing over a bump

Large‐eddy simulation (LES) and Reynolds‐averaged Navier–Stokes simulation (RANS) with different turbulence models (including the standard k?ε, the standard k?ω, the shear stress transport k?ω (SST k?ω), and Spalart–Allmaras (S–A) turbulence models) have been employed to compute the turbulent flow of a two‐dimensional turbulent boundary layer over an unswept bump. The predictions of the simulations were compared with available experimental measurements in the literature. The comparisons of the LES and the SST k?ω model including the mean flow and turbulence stresses are in satisfied agreements with the available measurements. Although the flow experiences a strong adverse pressure gradient along the rear surface, the boundary layer is unique in that intermittent detachment occurring near the wall. The numerical results indicate that the boundary layer is not followed by mean‐flow separation or incipient separation as shown from the numerical results. The resolved turbulent shear stress is in a reasonable agreement with the experimental data, though the computational result of LES shows that its peak is overpredicted near the trailing edge of the bump, while the other used turbulence models, except the standard k?ε, underpredicts it. Analysis of the numerical results from LES confirms the experimental data, in which the existence of internal layers over the bump surface upstream of the summit and along the downstream flat plate. It also demonstrates that the quasi‐step increase in skin friction is due to perturbations in pressure gradient. The surface curvature enhances the near‐wall shear production of turbulent stresses, and is responsible for the formation of the internal layers. The aim of the present work is to examine the response and prediction capability of LES with the dynamic eddy viscosity model as a sub‐grid scale to the complex turbulence structure with the presence of streamline curvature generated by a bumpy surface. Aiming to reduce the computational costs with focus on the mean behavior of the non‐equilibrium turbulent boundary layer of flow over the bump surface, the present investigation also explains the best capability of one of the used RANS turbulence models to capture the driving mechanism for the surprisingly rapid return to equilibrium over the trailing flat plate found in the measurements. Copyright © 2010 John Wiley & Sons, Ltd.     

Construction of a method for calculating a plane boundary layer in weak solutions of polymers with laminar,transitional, and turbulent flow zones

The investigation of the effect of small polymer additives on the characteristics of the flow of a viscous liquid is, at the present time, one of the most promising approaches to the lowering of the friction resistance. One interesting question in this connection is the study of the effect of small polymer additives on the characteristics of the transitional region of flow in a boundary layer, as well as on the value of the friction resistance with the presence of laminar, transitional, and turbulent sections in the boundary layer. The article sets forth a possible method for calculation of a plane boundary layer and the friction resistance for the case of the motion of a body in weak polymer solutions with a constant concentration, taking account of the change in the flow conditions in the layer and based on the use of integral relationships. Questions connected with the development of a boundary at a body, with the feeding of a polymer in it, as well as with the effect of degradation or destruction of the polymer in the solution, are not discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–48, May–June, 1977.     

The delay of transition of a laminar boundary layer over compliant wall

In a previous paper^[1], a method has been developed to study the stability characteristics of laminar boundary layers over compliant walls. In this paper, the effect of double layered compliant wall and Kramer type compliant wall on delaying the transition is investigated, and it is shown that there does exist the possibility to delay the transition by applying such compliant walls. The project supported by the National Natural Science Foundation of China.     

Instability of the separated shear layer in flow past a cylinder: Forced excitation

The receptivity of the separated shear layer for Re = 300 flow past a cylinder is investigated by forced excitation via an unsteady inflow. In order to isolate the shear layer instability, a numerical experiment is set up that suppresses the primary wake instability. Computations are carried out for one half of the cylinder, in two dimensions. The flow past half a cylinder with steady inflow is found to be stable for Re = 300. However, an inlet flow with pulsatile perturbations, of amplitude 1% of the mean, results in the excitation of the shear layer mode. The frequency of the perturbation of the inlet flow determines the frequency associated with the shear layer vortices. For a certain range of forced frequencies the recirculation region undergoes a low‐frequency longitudinal contraction and expansion. An attempt is made to relate this instability to a global mode of the wake determined from a linear stability analysis. Interestingly, this phenomenon disappears when the outflow boundary of the computational domain is shifted sufficiently downstream. This study demonstrates the need of carefully investigating the effect of the location of outflow boundaries if the computational results indicate the presence of low‐frequency fluctuations. The effect of Re and amplitude of unsteadiness at the inlet are also presented. All computations have been carried out using a stabilized finite element formulation of the incompressible flow equations. Copyright © 2007 John Wiley & Sons, Ltd.     

Control of cross flow in the three-dimensional boundary layer using space-periodic body force action

The possibility of attenuation of the cross flow in the three-dimensional incompressible laminar boundary layer on a sideslipping wing under the action of body force sources simulating the time-average forces generated by a surface electric discharge is estimated. The effect of the distance between the sources and the sideslip angle of the wing on the cross flow velocity is investigated for the source intensity observed experimentally.     

Residual‐based variational multiscale methods for laminar,transitional and turbulent variable‐density flow at low Mach number

In the present study, residual‐based variational multiscale methods are developed for and applied to variable‐density flow at low Mach number. In particular, two different formulations are considered in this study: a standard stabilized formulation featuring SUPG/PSG/grad‐div terms and a complete residual‐based variational multiscale formulation additionally containing cross‐ and Reynolds‐stress terms as well as subgrid‐scale velocity terms in the energy‐conservation equation. The proposed methods are tested for various laminar flow test cases as well as a test case at laminar via transitional to turbulent flow stages. Stable and accurate results are obtained for all numerical examples. Substantial differences in the results between the two approaches do not become notable until a high temperature gradient is applied and the flow reaches a turbulent flow stage. The more pronounced influence of adding subgrid‐scale velocity terms to the energy‐conservation equation on the results than adding analogous terms to the momentum‐conservation equation in this situation appears particularly noteworthy. Copyright © 2009 John Wiley & Sons, Ltd.     

A theoretical mode for the instability of turbulent boundary layer over compliant surface

A theoretical model for the instability of turbulent boundary layer over compliant surfaces is described. The investigation of instability is carried out from a time-asymptotic space-time perspective that classifies instabilities as either convective or absolute. Results are compared against experimental observations of surface waves on elastic and viscoelastic compliant layers.     

Numerical simulation of an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder

This article presents a numerical investigation of turbulent flow in an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder. The governing equations were discretized by the finite-volume method and SIMPLER method was applied to solve the equations on a staggered grid. The turbulent flow was numerically simulated using the standard k–ε, Abe–Kondoh–Nagano (AKN) and Shear Stress Transport (SST) turbulence models. The comparisons made between numerical results and experimental measurements showed that the SST model is superior to other models in the present calculation.Computations were performed for three different Reynolds numbers of 6000, 10 000 and 20 000 based on the cylinder diameter. To our knowledge, this study represents the first numerical investigation of the present flow configuration. The computational results were validated with the available experimental data of reattachment length, mean velocity distribution and wall static pressure coefficient in the turbulent blunt circular cylinder flows. Further, other characteristics of the flow, such as turbulent kinetic energy, pressure, streamlines, and the velocity vectors are discussed.The results show that the main characteristics of the turbulence flow in the separation region, such as reattachment length or velocity profiles, are nearly independent of the Reynolds number. The obtained results showed that a secondary separation bubble may appear in the main separation bubble near the leading edge. Furthermore, it was found that the turbulent kinetic energy has a large effect on the formation of the secondary bubble.     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq     

High-fidelity numerical simulation of the flow field around a NACA-0012 aerofoil from the laminar separation bubble to a full stall

In the present work, large eddy simulations of the flow field around a NACA-0012 aerofoil near stall conditions are performed at a Reynolds number of 5 × 10⁴, Mach number of 0.4, and at various angles of attack. The results show the following: at relatively low angles of attack, the bubble is present and intact; at moderate angles of attack, the laminar separation bubble bursts and generates a global low-frequency flow oscillation; and at relatively high angles of attack, the laminar separation bubble becomes an open bubble that leads the aerofoil into a full stall. Time histories of the aerodynamic coefficients showed that the low-frequency oscillation phenomenon and its associated physics are indeed captured in the simulations. The aerodynamic coefficients compared to previous and recent experimental data with acceptable accuracy. Spectral analysis identified a dominant low-frequency mode featuring the periodic separation and reattachment of the flow field. At angles of attack α ≤ 9.3°, the low-frequency mode featured bubble shedding rather than bubble bursting and reformation. The underlying mechanism behind the quasi-periodic self-sustained low-frequency flow oscillation is discussed in detail.     

Open and closed-loop experiments to identify the separated flow dynamics of a thick turbulent boundary layer

Open and closed-loop flow control experiments were performed on the transient attachment and separation mechanisms of a thick turbulent boundary layer (TBL). Without actuation, the TBL is subjected to an adverse pressure gradient and separates downstream of a sharp variation in the wall geometry. Departing from a given geometry and steady operations of vortex generator actuators, the control objective was to attach the flow in the separated region with a minimum of injected fluid using adaptation of the closed-loop control. The large scale of the facility (i.e., δ = 20 cm upstream of separation) induces large time scales and large Reynolds numbers of the flow to be controlled. It is found to consequently induce large time scales of the separation/attachment mechanisms, making the dynamic closed-loop implementation easier. Open-loop tests were first performed to extract the adequate input/output variables for closed-loop implementations. The chosen input variable was the Duty Cycle, DC, which enables sending of a control action at least 10 times faster than the time scales of the attachment/separation process. The chosen output variable was the voltage signal from a hot-film probe located on the flap which characterizes the degree of separation. In open loop, both the large scale (i.e., large time scales) of the present facility (Carlier and Stanislas in J Fluid Mech 535(36):143–188, 2005) and the well-defined excitation (Braud and Dyment in Phys Fluids 24:047102, 2012) help to extract the different time scales involved and to identify the whole system (actuators, baseline flow and sensor). Three Reynolds numbers based on the momentum thickness of the boundary layer near the actuators and upstream of separation were investigated (Re _θ  = 7,500, 10,500 and 12,600) through variation of the free-stream velocity (U _∞ = 5, 8, 10 m/s). These three systems were found to behave like first-order linear systems, with coefficients that need to be adapted depending on the Reynolds number. From Re _θ  = 7,500 to Re _θ  = 12, 600, the time scale and static gain of the linear system needed to be almost doubled. A simple controller (Proportional-Integral) was implemented in closed-loop configuration, improving the reactivity of the system. Robustness was tested by varying the free-stream velocity. Closed-loop control based on a fixed reference was unsuccessful as it failed to account for the effect of the Reynolds number. This was successfully overcome by tracking a given state of the flow using a simple P controller to adapt the reference according to variations of Re. The P controller, acting on the DC variable, compensates the corresponding variations of VR (ratio between the free-stream and the jet exit velocity).