空腔流动的动量分解及能量输运特性

空腔结构广泛应用于航空航天飞行器部件及地面交通工具中,其复杂的流声特性是相关工程设计中必须考虑的关键问题.空腔流动中的流声相互作用是空腔自持振荡的重要过程,准确识别并解耦空腔内的流体动力学模态和声模态,是深入理解空腔流声相互作用和能量转化机制的关键.通过直接求解二维Navier-Stokes方程数值模拟来流马赫数Ma=...     

Wavelet spatial scaling for educing dynamic structures in turbulent open cavity flows

In this work a methodology was developed for the selection of wavelet spatial scales to educe dynamic structures in turbulent cavity flows. The wavelet transform was applied to both the temporal signal and spatial fields to extract structures from the oscillating shear layer. The dominant frequencies were identified from the temporal transform of the shear layer oscillations, and then the corresponding wavelength was obtained using the relation U_cT=λ at each frequency. The wavelet spatial scaling was examined and a one-to-one relationship was established with respect to the wavelength. At each spatial scale, the transformed images of vorticity, velocity, and pressure fluctuations captured the vortical structures. Using this methodology, the dynamic vortical structures were extracted from the turbulent open cavity flows. Energy analysis was performed to examine the contributions of each structure.     

Vortex shedding and evolution induced by a solitary wave propagating over a submerged cylindrical structure

The shedding and evolution of the vortical structures generated by a solitary wave propagating over a submerged cylindrical structure are investigated experimentally and numerically. The cylindrical structure consists of two concentric cylinders and represents a simplified model for an offshore submerged intake structure typically used in coastal power plants. Flow visualization by dye injection is used to identify the dominant vortical structures near the structure. The flow visualization results show an unexpected flow reversal that causes shedding of secondary vortical structures. The experimental results are used to check a three-dimensional volume of fluid-large eddy simulation (VOF-LES) numerical model. The VOF-LES model is then used to further study the flow structure. A total of six dominant vortical structures generated by the wave motion are identified, followed by two more generated by the flow reversal. In summary, this paper provides the vorticity evolution for a complex fluid–structure interaction problem and a three-dimensional numerical simulation tool has also been validated, which can be extended to study more complex geometries and wave conditions.     

Flow field characteristics study of a flapping airfoil using computational fluid dynamics

The flow field of a flapping airfoil in Low Reynolds Number (LRN) flow regime is associated with complex nonlinear vortex shedding and viscous phenomena. The respective fluid dynamics of such a flow is investigated here through Computational Fluid Dynamics (CFD) based on the Finite Volume Method (FVM). The governing equations are the unsteady, incompressible two-dimensional Navier-Stokes (N-S) equations. The airfoil is a thin ellipsoidal geometry performing a modified figure-of-eight-like flapping pattern. The flow field and vortical patterns around the airfoil are examined in detail, and the effects of several unsteady flow and system parameters on the flow characteristics are explored. The investigated parameters are the amplitude of pitching oscillations, phase angle between pitching and plunging motions, mean angle of attack, Reynolds number (Re), Strouhal number (St) based on the translational amplitudes of oscillations, and the pitching axis location (x/c). It is shown that these parameters change the instantaneous force coefficients quantitatively and qualitatively. It is also observed that the strength, interaction, and convection of the vortical structures surrounding the airfoil are significantly affected by the variations of these parameters.     

Unsteady separated and reattaching turbulent flow over a two-dimensional square rib

The spatio-temporal characteristics of the separated and reattaching turbulent flow over a two-dimensional square rib were studied experimentally. Synchronized measurements of wall-pressure fluctuations and velocity fluctuations were made using a microphone array and a split-fiber film, respectively. Profiles of time-averaged streamwise velocity and wall-pressure fluctuations showed that the shear layer separated from the leading edge of the rib sweeps past the rib and directly reattaches on the bottom wall (x/H=9.75) downstream of the rib. A thin region of reverse flow was formed above the rib. The shedding large-scale vortical structures (fH/U₀=0.03) and the flapping separation bubble (fH/U₀=0.0075) could be discerned in the wall-pressure spectra. A multi-resolution analysis based on the maximum overlap discrete wavelet transform (MODWT) was performed to extract the intermittent events associated with the shedding large-scale vortical structures and the flapping separation bubble. The convective dynamics of the large-scale vortical structures were analyzed in terms of the autocorrelation of the continuous wavelet-transformed wall pressure, cross-correlation of the wall-pressure fluctuations, and the cross-correlation between the wall pressure at the time-averaged reattachment point and the streamwise velocity field. The convection speeds of the large-scale vortical structures before and after the reattachment point were U_c=0.35U₀ and 0.45U₀, respectively. The flapping motion of the separation bubble was analyzed in terms of the conditionally averaged reverse-flow intermittency near the wall region. The instantaneous reattachment point in response to the flapping motion was obtained; these findings established that the reattachment zone was a 1.2H-long region centered at x/H=9.75. The reverse-flow intermittency in one period of the flapping motion demonstrated that the thin reverse flow above the rib is influenced by the flapping motion of the separation bubble behind the rib.     

SPIV measurements around the DELFT 372 catamaran in steady drift

Secondary flow vortical patterns in arterial curvatures have the potential to affect several cardiovascular phenomena, e.g., progression of atherosclerosis by altering wall shear stresses, carotid atheromatous disease, thoracic aortic aneurysms and Marfan’s syndrome. Temporal characteristics of secondary flow structures vis-à-vis physiological (pulsatile) inflow waveform were explored by continuous wavelet transform (CWT) analysis of phase-locked, two-component, two-dimensional particle image velocimeter data. Measurements were made in a 180° curved artery test section upstream of the curvature and at the 90° cross-sectional plane. Streamwise, upstream flow rate measurements were analyzed using a one-dimensional antisymmetric wavelet. Cross-stream measurements at the 90° location of the curved artery revealed interesting multi-scale, multi-strength coherent secondary flow structures. An automated process for coherent structure detection and vortical feature quantification was applied to large ensembles of PIV data. Metrics such as the number of secondary flow structures, their sizes and strengths were generated at every discrete time instance of the physiological inflow waveform. An autonomous data post-processing method incorporating two-dimensional CWT for coherent structure detection was implemented. Loss of coherence in secondary flow structures during the systolic deceleration phase is observed in accordance with previous research. The algorithmic approach presented herein further elucidated the sensitivity and dependence of morphological changes in secondary flow structures on quasiperiodicity and magnitude of temporal gradients in physiological inflow conditions.     

A vortex axis and vortex core border grid adaptation algorithm

In the design process of hydrodynamical and aerodynamical technical applications, the numerical simulation of massively separated vortical flow is crucial for predicting, for example, lift or drag. To obtain reliable numerical results, it is mandatory to accurately predict the physical behavior of vortices. Thus, the dominant vortical flow structures have to be resolved in detail, which requires a local grid refinement and certain adaptation techniques. In this paper, a vortex flow structure adaptation algorithm is presented, which is particularly designed for local grid refinement at vortex axes positions and associated vortex core border locations. To this end, a fast and efficient vortex axis detection scheme is introduced and the algorithm for the vortex core border determination is explained. As the interaction between vortices makes the assignment of grid points to a certain vortex axis difficult, a helicity‐based vortex distinction approach in combination with a geometrical rotational sensor is developed. After describing the combined different techniques in detail, the vortex feature adaptation algorithm is applied to analytical and more realistic examples, which show that the described grid adaptation algorithm is able to enhance the grid cell resolution locally such that all significant vortical flow phenomena are resolved. Copyright © 2008 John Wiley & Sons, Ltd.     

Three-dimensional numerical study of jet-in-crossflow characteristics at low Reynolds number

A numerical study of a square jet in a cross flow is carried out at a Reynolds number of 100. The flow field and heat transfer characteristic downstream of the jet have been explored by solving three-dimensional unsteady Navier–Stokes equations and energy equation using higher order spatial and temporal discretization. The projection of vortical structure on a plane is seen to give the component of vortex normal to the plane. Four combinations of velocity profile namely (1) uniform crossflow and uniform jet, (2) laminar boundary layer crossflow and uniform jet, (3) uniform crossflow and parabolic jet profile, and (4) laminar boundary layer crossflow and parabolic jet are compared at same phase to see their effect on the flow field and heat transfer characteristic. All the four cases are seen to exhibit unsteadiness but the jet with parabolic profile is seen to show stronger unsteadiness. The instantaneous vortical structures of all the cases at the same phase show that the structures are more complex for the jet with parabolic velocity profile. The temperature field is seen to be correlated with the vortical structures. Comparison of the time averaged flow field reveals that the jet penetration is the highest for the jet having parabolic profile and boundary layer crossflow. The adiabatic effectiveness is observed to be more for the jet with uniform velocity profile and uniform crossflow and was least for the jet with parabolic velocity profile and boundary layer crossflow.     

Numerical investigation of flow structures around a sphere

A numerical investigation of flow around a sphere is performed and compared with previous studies. Here, a second-order accurate, finite volume method is used in order to predict the instantaneous and time-averaged flow characteristics using large eddy simulation (LES) on the multi-block grid system. Namely, the objectives of this article are: (i) the presentation of flow structures in the wake region downstream of the sphere with a wide variety of flow properties such as the distribution of velocity vectors, patterns of streamlines, Reynolds stress correlations, root mean square of velocity components and other time-averaged flow data in order to reveal the vortical flow structures in detail and (ii) to demonstrate the abilities of computational methods in simulation of vortical flow data. Finally, it has been concluded that there are good agreements between the experimental results and numerical predictions.     

Eduction of vortical structures by pressure measurements in noncircular jets

Vortical structures in the noncircular jets excited at the interaction mode were educed by measuring fluctuating static pressure, and their characteristics were discussed in relation to the mixing mechanism of the noncircular jets. The contours of phase-average fluctuating pressure show clearly the vortical structures of the noncircular jets, which compare reasonably with those observed by the flow visualization experiments. The evolution of the vortical structures is characterized by partial merging, stretching and splitting to smaller vortices. The effects of the noncircular vortical structures on mixing were also discussed based on the quantitative measurements of the velocity fields, the results suggesting that the interactions of vortical structures in the noncircular jets are very effective to enhance mixing.     

Vortical structures and pressure pulsations in draft tube of a Francis-99 turbine at part load: RANS and hybrid RANS/LES analysis

Recognizing the limitations of the conventional linear-eddy-viscosity (LEVM) Reynolds-averaged Navier–Stokes (RANS) models to reproduce complex three-dimensional unsteady flows in hydraulic machinery, we performed a comparative assessment of a second-moment (Re-stress model, RSM) RANS closure and a hybrid RANS/LES method in capturing the flow and vortical structures in the draft tube of a Francis hydroturbine at off-design conditions. Considered is a case of part load (PL) at a flow rate of only 35% of the best efficiency point (BEP) characterised by multiple unsteady vortex systems. Despite some remaining uncertainties in generating the inflow conditions, both approaches reproduced reasonably well the measured mean velocity and the rms of its fluctuations, as well as the pressure spectrum with peaks detecting the precessing vortex core. In contrast to the common LEVMs, the Re-stress closure showed sufficient receptivity to intrinsic unsteadiness and reproduced well the overall flow and vortical patterns as well as the associated pressure pulsations in accord with the experiments. The hybrid RANS/LES method gave similar predictions as the RSM, but resolving a wider range of scales, which however, showed no significant effect on the dynamics of the dominant processing vortex core and the pressure pulsations.     

Effects of the jet-to-crossflow momentum ratio on a sonic jet into a supersonic crossflow

Numerical investigation of a transverse sonic jet injected into a supersonic crossflow was carried out using large-eddy simulation for a free-stream Mach number M = 1.6 and a Reynolds number Re = 1.38 × 10⁵ based on the jet diameter. Effects of the jet-to-crossflow momentum ratio on various fundamental mechanisms dictating the intricate flow phenomena, including flow structures, turbulent characters and frequency behaviors, have been studied. The complex flow structures and the relevant flow features are discussed to exhibit the evolution of shock structures, vortical structures and jet shear layers. The strength of the bow shock increases and the sizes of the barrel shock and Mach disk also increase with increasing momentum ratio. Turbulent characters are clarified to be closely related to the flow structures. The jet penetration increases with the increase of the momentum ratio. Moreover, the dominant frequencies of the flow structures are obtained using spectral analysis. The results obtained in this letter provide physical insight in understanding the mechanisms relevant to this complex flow.     

Characteristics of wall pressure fluctuations in separated flows over a backward-facing step: Part II. Unsteady wavelet analysis

 Time-dependent characteristics of wall pressure fluctuations in separated and reattaching flows over a backward-facing step were investigated by means of the continuous wavelet transform. Emphasis was placed on the combination of time-localized analyses of the wavelet transform and multi-point measurements of pressure fluctuations. Synchronized wavelet maps revealed the evolutionary behavior of pressure fluctuations and gave further insight into the modulated nature of large-scale vortical structures. It was found that there exist two modes of shed vortices: one is the global oscillation and the other is the vortex convection. The two alternating modes are synchronized with the flapping frequency component of pressure fluctuations. The flapping motion gives rise to the difference in pressure spectra, indicating more intensive pressure activity during the shrinking period of the recirculation region.     

Correlation between vortex structures and unsteady loads for flapping motion in hover

During the past decade, efforts were made to develop a new generation of unmanned aircrafts, qualified as Micro-Air Vehicles. The particularity of these systems resides in their maximum dimension limited to 15 cm, which, in terms of aerodynamics, corresponds to low Reynolds number flows (Re ≈ 10² to 10⁴). At low Reynolds number, the concept of flapping wings seems to be an interesting alternative to the conventional fixed and rotary wings. Despite the fact that this concept may lead to enhanced lift forces and efficiency ratios, it allows hovering coupled with a low-noise generation. Previous studies (Dickinson et al. in Science 284:1954–1960, 1999) revealed that the flow engendered by flapping wings is highly vortical and unsteady, inducing significant temporal variations of the loads experienced by the airfoil. In order to enhance the aerodynamic performance of such flapping wings, it is essential to give further insight into the loads generating mechanisms by correlating the spatial and temporal evolution of the vortical structures together with the time-dependent lift and drag. In this paper, Time Resolved Particle Image Velocimetry is used as a basis to evaluate both unsteady forces and vortical structures generated by an airfoil undergoing complex motion (i.e. asymmetric flapping flight), through the momentum equation approach and a multidimensional wavelet-like vortex parameterization method, respectively. The momentum equation approach relies on the integration of flow variables inside and around a control volume surrounding the airfoil (Noca et al. in J Fluids Struct 11:345–350, 1997; Unal et al. in J Fluids Struct 11:965–971, 1997). Besides the direct link performed between the flow behavior and the force mechanisms, the load characterization is here non-intrusive and specifically convenient for flapping flight studies thanks to its low Reynolds flows’ sensitivity and adaptability to moving bodies. Results are supported by a vortex parameterization which evaluates the circulation of the multiple vortices generated in such complex flows. The temporal evolution of the loads matches the flow behavior and hence reveals the preponderant inertial force component and that due to vortical structures.     

Three-dimensional numerical simulation of a vortex ring impacting a bump

A vortex ring impinging on a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re = 4 × 10⁴ based on the initial translation speed and diameter of the vortex ring. The effects of bump height on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Furthermore, the mechanism of flow transition from laminar to turbulent state has been revealed based on analysis of turbulent kinetic energy.     

Assessment of the organization of a turbulent separated and reattaching flow by measuring wall pressure fluctuations

The effect of local forcing on the organization of a turbulent separated and reattaching flow was assessed by measuring wall pressure fluctuations. Multi-arrayed microphones were installed on the surface to measure the simultaneous spatial and temporal wall pressure fluctuations. Local forcing at the separation edge was applied to the separated flow over a backward-facing step through a thin slit. The organization of the separated and reattaching flow was found to be greatest at the effective forcing frequency. The flow structure was diagnosed by analyzing several characteristics of the wall pressure fluctuations: the wall pressure fluctuation coefficients, wall pressure spectrum, wavenumber-frequency spectrum, coherence, cross-correlation, and multi-resolution autocorrelations of pressure fluctuations using the maximum overlap discrete wavelet transform and continuous wavelet transform. Features indicative of the amalgamation of vortices under the local forcing were observed; this amalgamation process accounted for the observed reduction of the reattachment length. Examination of the wall pressure fluctuations revealed that introduction of local forcing enhanced flapping motion as well as the streamwise and spanwise dispersions of vortical structures.     

Vortical patterns in turbulent flow downstream a 90° curved pipe at high Womersley numbers

The present experimental work focuses on highly pulsatile, i.e. inertia dominated, turbulent flow downstream a curved pipe and aims at investigating the vortical characteristics of such a flow. The flow parameters (Dean and Womersley number) investigated are of the same order as those met in the internal combustion engine environment. The technique employed is time-resolved stereoscopic particle image velocimetry at different cross-sections downstream the pipe bend. These measurements allow the large-scale structures that are formed to be analyzed by means of proper orthogonal decomposition. The flow field changes drastically during a pulsatile cycle, varying from a uniform flow direction across the pipe section from the inside to the outside of the bend to vortical patterns consisting of two counter-rotating cells. This study characterizes and describes pulsatile curved pipe flow at Womersley numbers much higher than previously reported in the literature. Furthermore, the oscillatory behaviour of the Dean cells for the steady flow – the so-called ‘swirl switching’ – is investigated for different downstream stations from the bend exit and it is shown that this motion does not appear in the immediate vicinity of the bend, but only further downstream.     

Direct numerical simulation of a turbulent Couette-Poiseuille flow,part 2: Large- and very-large-scale motions

A direct numerical simulation dataset of a fully developed turbulent Couette-Poiseuille flow is analyzed to investigate the spatial organization of streamwise velocity-fluctuating u-structures on large and very large scales. Instantaneous and statistical flow fields show that negative-u structures with a small scale on a stationary bottom wall grow throughout the centerline due to the continuous positive mean shear, and they penetrate to the opposite moving wall. The development of an initial vortical structure related to negative-u structures on the bottom wall into a large-scale hairpin vortex packet with new hairpin vortices, which are created upstream and close to the wall, is consistent with the auto-generation process in a Poiseuille flow (Zhou et al., J. Fluid Mech., vol. 387, 1999, pp. 353–396). Although the initial vortical structure associated with positive-u structures on the top wall also grows toward the bottom wall, the spatial development of the structure is less coherent with weak strength due to the reduced mean shear near the top wall, resulting in less turbulent energy on the top wall. The continuous growth of the structures from a wall to the opposite wall explains the enhanced wall-normal transport of the streamwise turbulent kinetic energy near the centerline. Finally, an inspection of the time-evolving instantaneous fields and conditional averaged flow fields for the streamwise growth of a very long structure near the centerline exhibits that a streamwise concatenation of adjacent large-scale u-structures creates a very-large-scale structure near the channel centerline.     

Large-scale structures of turbulent flows over an open cavity

Low Mach number turbulent flows over an open cavity were studied to investigate the quantitative characteristics of large-scale vortical structures responsible for self-sustained oscillations. Wind tunnel experiments with particle image velocimetry (PIV) were conducted in the range of the ratio of cavity length (L) to depth (D), 1<L/D<4, when the incoming boundary layer is turbulent at Re_θ=830 and 1810. Self-sustained oscillation modes were classified by varying the conditions of L/D and Re_θ. The oscillation modes were consistent with the number of vortical structures existing between the leading and trailing edges of the cavity. Proper orthogonal decomposition (POD) was employed to the spatial distributions of vertical velocity correlations on the lip line of cavity geometry. By examining the conditionally averaged distributions of the correlation coefficients of POD, the spatial characteristics of large-scale vortical structures for self-sustained oscillations were examined.     

Vortex identification methods based on temporal signal-processing of time-resolved PIV data

The lack of a universally accepted mathematical definition of a vortex structure has led to a considerable number of Eulerian criteria to identify coherent structures. Most are derived from the instantaneous local velocity gradient tensor and its derivatives and require appropriate thresholds to extract the boundaries of the structures. Notwithstanding their great potential for studying coherent structures, most criteria are not frame-independent and they lack a clear physical meaning. The Lyapunov exponent, a popular tool in dynamical system theory, appears as a promising alternative. This Lagrangian criterion does not suffer from the drawbacks of the Eulerian criteria and is constructed on a simple physical interpretation that includes information on the history of the flow. However, since the computation of the Lyapunov exponent involves the knowledge of fluid particle trajectories, experimental applications are currently restricted to laminar flows and two-dimensional turbulence, provided that velocity fields are time-resolved. In this work, we explore temporal post-treatment methods to extract vortical structures developing in a flow through a smooth axisymmetric constriction. Data from planar time-resolved Particle image velocimetry, measuring two or three components of the velocity vectors, are transformed via the Taylor hypothesis to quasi-instantaneous three-dimensional velocity field and are interpreted in terms of the discrete wavelet decomposition, the finite-time Lyapunov exponent, and the linear stochastic estimation. It appears that these methods can concurrently provide very rich and complementary scalar fields representing the effects of the vortical structures and their interactions in the flow.