Averaged Flow Field Research Articles

A continuously scanning spatiotemporal averaging method for obtaining volumetric mean flow measurements with stereoscopic PIV

In this study, a variation of the stacked stereoscopic PIV technique is proposed to perform fully volumetric (3-dimensions, 3-components) measurements of average flow fields within a single experiment through the usage of an automated traversing system that continuously scans the SPIV light sheet over a linear path. The simultaneous measurement of the traverse location and the laser Q-switch pulse enables the automated assignment of instantaneous PIV fields to known physical coordinates, enabling spatiotemporal averaging in post-processing to obtain volumetric measurements of a flow field. This method provides a trade-off between spatial resolution of the volume measurements and statistical convergence of the spatiotemporal averages, enabling volumetric measurements under challenging experimental conditions where only stereoscopic PIV is viable. A comparison with the more traditional temporal averaging method and planar PIV is presented to demonstrate the capabilities and limitations of this technique in realistic, challenging experimental setups. It is found that the spatiotemporal averaging convergence behavior differs slightly from the traditional temporal averaging for the wake of a bluff body model, however relative errors lower than two standard deviations can still be attained. Thus, this technique presents a viable alternative for rapid 3D reconstruction of averaged flow fields that can provide invaluable insight of various flow topologies.

Multi-layer flow field mapping in a small-scale shallow water reservoir by coastal acoustic tomography

Long-term, high-precision and spatiotemporal observation of current structure in small-scale waters is important for underwater research and water resources management. In this paper, a reciprocal sound transmission experiment was conducted using three Coastal Acoustic tomography (CAT) systems within 300 m spacing in Huangcai reservoir, Hunan, China from 14 to 16, Sep 2020. The multi-arrival reciprocal paths between each pair of stations were identified with the help of ray simulation. 3-D CAT data which penetrated the full depth of water were obtained. The flow field along vertical slice were mapped and showed consistency with Acoustic Doppler Current Profiler (ADCP) and Temperature Depth (TD) data. A standard deviation of 0.73 m/s at the bottom layer were observed, indicating strong bottom flows in the reservoir. Furthermore, to visualize the strong bottom flow, the horizontal flow fields in the bottom layer were mapped by 3-D inversion technique. The feasibility of the technique was validated by three measures: 1) comparison between range average current and horizontal flow field; 2) Uncertainty analysis; 3) Check net volume transport. This work proves the CAT combined with 3-D inversion technique has the potentials to provide real-time 3-D current data for water researches.

Investigation on frequency influence on the transverse pulsed jet in a supersonic crossflow

AbstractThe pulsed jet is a novel and effective active mixing enhancement approach. For the transverse pulsed jet in the supersonic crossflow, the frequency influence is investigated using the three-dimensional Reynolds-averaged Navier–Stokes (RANS) equations coupled with the SST k-ω turbulence model. The averaged flow field properties of the pulsed jet are better than those of the steady jet when considering mixing efficiency and jet penetration depth, especially for the case with the pulsed frequency being 50kHz. The flow field structures of the pulsed jet are connected with the time, with periodic wave structures generating in the flow field and moving downstream. The size of the wave structures and its distance are related to the frequency, namely the size and flow distance are relatively small at 50kHz, and it takes some time for the pulsed jet to establish its influence in the full flow field. At low frequencies, the flow field produces large fluctuations, and this may be detrimental to the stable operation of the engine.

Mechanistic study on modification of convective internal boundary layer by spanwise surface motion

Modification of a convective internal boundary layer (IBL) by spanwise motion of a warm surface is investigated by imposing different surface moving speeds in the present study. Our analysis shows that the spanwise surface motion reduces the Reynolds shear stress right after the increase in the surface temperature in the convective IBL. The maximum decreasing rate of the Reynolds shear stress is found to be approximately 75% at the largest moving speed of the warm surface considered in the manuscript. Due to the reduction of the Reynolds shear stress, the vertical momentum transport is fundamentally altered, and the mean flow accelerates immediately after the increase in the surface temperature. By scrutinizing the instantaneous and conditional averaged flow fields as well as the pre-multiplied energy spectra, we have attributed the reduction of the Reynolds shear stress to the suppression of the near-surface velocity streaks and quasi-streamwise vortices, and the delayed growth of the convective structures, such as thermal plumes. Our investigation suggests that the developments of the convective IBL can be influenced by a strong spanwise motion of the warm surface, which should be taken into consideration in the prediction model for practical applications.

CFD Numerical Simulation of Wind Field and Vehicle Aerodynamic Characteristics on Truss Bridge Deck under Crosswind

Due to the complexity and changeability of the wind field in deep-cut gorges, the vehicles on the bridge deck are easily affected by a strong crosswind. Thus, to accurately evaluate the wind field characteristics of a suspension bridge deck, the wind profiles of wind speed and angle of attack (AoA) and the vehicle aerodynamic parameters were investigated by the CFD. The results show that the shape of wind speed profile is mainly controlled by the AoA but less affected by the Reynolds number. The main girder's shielding effect can accelerate the local wind field, and the closer to the windward lane, the less the interference; thus, a suitable location of measuring points to represent the incoming flow is found. Furthermore, the equivalent wind speed based on the equivalent side force is generally larger than the value based on the rolling moment, and the responding value is greatly affected by the AoA. In addition, the vortex is a time-dependent phenomenon, the averaged flow field produced less force in the wake, but the high local wind speed variations may affect the traffic unfavorably. The results provide an essential reference significance for studying the local wind field characteristics of the bridge deck.

Effect of expansion on the wall heat flux in a supersonic turbulent boundary layer

Direct numerical simulation of a spatially developing supersonic turbulent boundary layer at a Mach number of 2.25 and a friction Reynolds number of Reτ = 769 subjected to an expansion corner with a deflection angle of 12° is performed to investigate the effect of expansion on the characteristics of the wall heat flux (WHF). The effect of expansion on the statistical and structural properties of the fluctuating WHF is analyzed systematically in terms of probability density function, frequency spectra, and space-time correlations. Normalization using the local root mean square value yields good collapse of the probability density function curves. Unlike with wall pressure frequency spectra, it is found that expansion has little influence on the low-frequency components of the WHF spectrum. The correlation results show that the main effect of expansion is to increase the characteristic length scales and convection velocity of the WHF fluctuation in the post-expansion region. Furthermore, a direct comparison between the conditionally averaged flow fields and those presented in the authors' previous work [Tong et al., Phys. Fluids 34, 015127 (2022)] is performed to uncover the effect of expansion on the flow structures associated with extreme positive and negative WHF fluctuation events. We highlight that the extreme positive event emerges below a small hot spot under the action of a strong Q4 event, whereas the extreme negative event is relatively insensitive to expansion and still occurs between a pair of strong oblique vortices. In addition, we decompose the mean WHF into seven physics-informed contributions and quantify the effect of expansion on the dominating components with the aid of the bidimensional empirical mode decomposition method. The scale-decomposed results demonstrate quantitatively that expansion decreases the contribution of the large-scale structures in the outer region but the small-scale structures in the near-wall region contribute heavily to the mean WHF generation in the downstream region.

Numerical Simulation of Flow and Heat Transfer of a Discontinuous Single Started Helically Ribbed Pipe

In the present study, the turbulent flow field and the heat transfer in a single started helically ribbed pipe with a discontinuous rib are investigated. A large-eddy simulation (LES) technique is applied in a pipe section with cyclic boundary conditions. The aim of this study is to explain and further analyze the findings from the heat transfer measurements at such complex structures with the help of detailed flow simulations. The simulation results are validated with measurements at a Reynolds number of Re = 21,100 and a Prandtl number of Pr = 7 with water as fluid. The comparison clearly shows that the current method delivers accurate results concerning average flow field, turbulence quantities and local heat transfer. The results demonstrate that the applied method is capable of correctly simulating flows with heat transfer in complex three-dimensional structures. The overall heat transfer performance of the helically ribbed pipe with a discontinuous rib is compared to a smooth pipe and a continuous rib configuration. The impact of the interruption of the rib structure on pressure drop and heat transfer are analyzed in detail.

Noise simulations of flap devices for wind turbine rotors

AbstractThe rotor of a large diameter wind turbine experiences more substantial and more dynamic loads due to the fluctuating and heterogeneous wind field. The project SmartBlades 2.0 investigated rotor blade design concepts that alleviate aerodynamic loading using active and passive mechanisms. The present work evaluates the acoustics of the two load alleviating concepts separately, an inboard slat and an outboard flap, using the Fast Random Particle Mesh/Fast Multipole Code for Acoustic Shielding (FRPM/FMCAS) numerical prediction toolchain developed at DLR with input from the averaged flow field from RANS. The numerical tools produce a comparable flap side‐edge noise spectrum with that of the measurement conducted in the Acoustic Wind Tunnel Braunschweig (AWB). The validated FRPM/FMCAS was then used to analyze the self‐noise from a slat at the inboard section of a rotor blade with a 44.45 m radius and compared with that from the outboard trailing edge. Furthermore, the rotational effect of the rotor was included in the post‐processing to emulate the noise observed at ground level. The findings show an increase in the slat's overall sound pressure level and a maximum radiation upwind of the wind turbine for the case with the largest wind speed that represents the off‐design condition. In operational conditions, the slat adds at most 2 dB to the overall sound pressure level. The toolchain evaluates wind turbine noise with conventional or unconventional blade design, and the problem can be scaled up for a full‐scale analysis. As such, the tools presented can be used to design low‐noise wind turbines efficiently.

Currents shaped the genetic structure of Tridentiger barbatus (Günther 1861) along the Chinese coast

The modern coast of China is characterized by complicated hydrology, which might be a potential factor shaping genetic structure of marine organisms. Tridentiger barbatus is a widespread benthic goby in coastal China whose range expansion might depend on larval dispersal in its spawning season. We explored the genetic structure of T. barbatus and its interrelationships with coastal currents, planktonic larval duration, and Pleistocene climate fluctuations. Genealogies of Cytb, Rh, and RAG1 reveal two lineages in T. barbatus, namely the southern and northern lineages. The southern lineage is distributed through the South China Sea to the South Yellow Sea but is absent from the North Yellow Sea and Bohai Sea. The HYCOM average flow field indicates an eddy between the North and South Yellow Seas, which might cause larval retention and limit dispersal of the southern lineage into the North Yellow Sea and the Bohai Sea. The northern lineage is distributed through the South Yellow Sea to the Bohai Sea but is absent from the South China Sea and East China Sea. The distribution of the northern lineage fits the eastward deflection of the Yellow Sea Coastal Current at the boundary of the East China Sea and South Yellow Sea, which might limit southward larval dispersal of the northern lineage into the East China Sea. We also uncovered two Cytb haplotype groups within the southern lineage, namely groups A and B, whose distributions might be shaped by the northward Taiwan Strait Current. The coincidence between the genetic structure of T. barbatus and the coastal currents and eddies strongly implies significant roles of larval dispersal and current transportation in shaping the genetic structure of T. barbatus. In addition, significant isolation by distance is found in the southern lineage, which coincides with its short planktonic larval duration (20–30 days). The inference on divergence time suggests that the observed lineages and haplotype groups in T. barbatus originated in the Pleistocene glaciations. Therefore, the coastal currents, planktonic larval duration, and historical glacial isolation are proposed as important factors shaping the genetic structure of T. barbatus.

Unsteady Analysis of a Pulsating Alternate Flow Pattern in a Radial Vaned Diffuser

In centrifugal compressors, Mild Surge (MS) leads to unstable operation. Previous experimental work on a centrifugal compressor designed and built by Safran Helicopter Engines (SafranHE) showed that MS corresponds to the pulsation of an alternate stall pattern at the Helmholtz frequency of the test rig on two channels in the radial diffuser. The present contribution experimentally investigates the impact of the Inlet Guide Vane (IGV) stagger angle on this alternate flow and numerically studies the topology of this pulsating alternate flow. The experimental investigation is performed with unsteady pressure sensors, and shows that the IGV stagger angle only impacts the pulsation frequency of the alternate flow pattern. This change is explained by the dependence of the Helmholtz frequency on the compressor inlet section. The topological analysis of the average flow field, computed from wall-resolved Unsteady Reynolds-Averaged Navier–Stokes (URANS) simulations, demonstrates that the saddle point (major critical point) in the corner hub/suction side of the stalled blade migrates upstream while staying in the corner if the mass flow rate decreases. One main blade over two is stalled on both sides because the flow originating from this corner separation circumvents the trailing edge and migrates upstream along the pressure side. In the simulation, the pulsation of the alternate stall is coupled with the reflection of acoustic waves on the inlet and outlet planes, regarded as an environmental effect.

Stereoscopic PIV Measurements Of The Human Thermal Plume

A human thermal plume contributes significantly towards indoor fluid flows. Such flows impact the ventilation system designs like mixing and displacement ventilation. Further, the human thermal plume influences the dispersion of air constituents to and away from the breathing zone. This is also relevant in the transport of respiratory droplets which results in spread of respiratory diseases including COVID-19. One of the general scenarios in many ventilated spaces is the sitting condition which is the focus here. Experiments were performed with a thermal manikin and comparisons were made with a real human. The thermal plume is visualised by Stereoscopic PIV with the visualisation plane located 30 cm above the head parallel to ground. The results presented here show that the thermal plume from a manikin without breathing function have many similarities in the average flow field above the head with the field induced by a real person. The peak velocity observed with the real person is slightly lower than the one of the thermal manikin and the structures of the thermal plumes resemble to some extent in both cases. The real person with breathing or talking revealed a bi-lobe like structure with the secondary lobe due to the breath. The thermal manikin also showed a second lobe which possibly is due to the thermal plume from the legs. Furthermore, the effect due to breathing, speaking, and movement of the real person are analysed. The integral values do not vary with respect to the different situations except for the movement case. The characterizations of these phenomena are important because they cannot be simulated by the manikin, while they affect the extent and strength of the thermal convection flow.

FLOW Estimation and Rose Superposition (FLOWERS): an integral approach to engineering wake models

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each discrete wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions and speeds. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed 2 orders of magnitude faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a relatively smooth solution space for the gradient-based optimizer to excel in comparison to the existing weighted-averaging power calculation.

Comparison of Flow Behavior in Saccular Aneurysm Models Using Proper Orthogonal Decomposition

Aneurysms are abnormal ballooning of a blood vessel. Previous studies have shown presence of complex flow structures in aneurysms. The objective of this study was to quantify the flow features observed in two selected saccular aneurysm geometries over a range of inflow conditions using Proper Orthogonal Decomposition (POD). For this purpose, two rigid-wall saccular aneurysm models geometries were used (i.e., the bottleneck factor of 1 and 1.6), and the inflow conditions were varied using a peak Reynolds number (Rep) from 50 and 270 and Womersley number (α) from 2 and 5. The velocity flow field data for the studied aneurysm geometries were acquired using Particle Image Velocimetry (PIV). The average flow field from the PIV measurement showed that the model geometry and Rep have more significant impact on the average flow field than the variations in α. The POD results showed that the method was able to quantify the flow field characteristics between the two model geometries. The mode shapes obtained showed different spatial structures for each inflow scenarios and models. The POD energy results showed that more than 80% of the fluctuating kinetic energy were captured within five POD modes for BF=1.0 flow scenarios, while they were captured within ten modes for BF=1.6. The time varying coefficient results showed the complex interplay of POD modes at different inflow scenarios, highlighting important modes at different phases of the flow cycle. The low-order reconstruction results showed that the vortical structure either proceeded outward or stayed within the aneurysm, and this behavior was highly dependent on α, Rep, and model geometry that were not evident in average PIV results.

Time-averaged flow characteristics of pulsed swirling coaxial jets with annular blockage

An experimental study was conducted to investigate the mixing and time-averaged turbulence characteristics of pulsed swirling coaxial jets with annular blockage. The study was conducted using a custom-made device by pulsing the central jet using acoustic excitation. A high-speed camera was used to identify the flow structures. The mixing characteristics and the turbulence properties of the flow were studied through a hot-wire anemometer and high-speed particle image velocimetry (PIV) techniques. Flow visualization, mean velocity profiles, vorticity and Reynolds stress fields, and Lagrangian time and length scales were used to investigate the turbulence and mixing characteristics of the flow. The flow was excited under two distinct conditions. At the lower pulsation intensity, mixing was shown to improve near the exit because of the entrainment induced during the growth of the periodic vortex rings. At the higher pulsation intensity, mixing was shown to improve further downstream due to increased jet spreading caused by the breakup of the vortex rings. Investigation of the ensemble average flow field showed that the structure of the unexcited flow was altered because of the excitation. At the higher pulsation intensity, a stagnation point was formed on the centerline which was not observed in the unexcited flow and excitation at the lower pulsation intensity. Within a single period, the position of the stagnation point moved cyclically along the centerline. In comparison with the naturally evolving jets, excitation at the higher pulsation intensity has reduced the mean axial velocity allowing better mixing of the jets in the downstream area.

Numerical study of a pair of spheres in an oscillating box filled with viscous fluid

When two identical spherical particles are submerged in an oscillating fluid, they align themselves perpendicularly to the direction of the flow, leaving a small gap between them. The formation of this structure is attributed to a non-zero residual flow known as steady streaming. We performed direct numerical simulations of a fully resolved, oscillating flow in which the pair of particles is modeled using an immersed boundary method. In equilibrium, the particles oscillate both parallel and perpendicularly to the oscillating flow. Two scaling regimes for the particle dynamics, which are related to changes in the (average) flow field, are observed.

Coherent structures decomposition of the flow field in Francis turbine runner under different working conditions

The dynamic mode decomposition (DMD) method is adopted to do flow field decomposition in the runner of Francis turbine. The coherent structures corresponding to main frequency components in runner passage under different working conditions are extracted, and the relative error of the reconstructed flow field and the predicted flow field are analyzed. The results show that the DMD method successfully decomposes the position and shape of the coherent structure corresponding to a single frequency in the unsteady flow, and its temporal and spatial evolution. Under part load (PL), best efficiency point (BEP), and high load (HL) working conditions, the coherent structures with high frequency induced by rotor-stator interaction (RSI) are the most stable, which located at the entrance of the runner passage. The low-frequency coherent structures mostly appear inside the runner passage, which is mainly caused by the short blade wake and flow separation under off-design conditions. The averaged flow field contains the most flow energy (larger than 90%), and unsteady coherent structures consume or transport less energy. The flow field in runner passage predicted by the DMD method is consistent with the computational flow field, which relative error is stable at 10−2. The study of coherent structures in runner passage of Francis turbine is beneficial to achieve a better design of runner and more stable operation of power plants.

Supersonic turbulent channel flows over spanwise-oriented grooves

In the present study, we perform direct numerical simulations to study the influences of the spanwise-oriented grooves, which are emulated by the reasonably designed “relaxed” boundary conditions, on the kinetic and thermodynamic statistics in a supersonic turbulent channel flow at the Mach number of 1.5 and Reynolds number of 3000. The phase averaged flow fields show that the relaxed boundary induces compressive and expansive waves that travel across the whole channel and are reflected by the upper wall. These waves are isentropic in the average sense except in the viscous sublayer. In the near-wall region, vortices and streaks that constitute the self-sustaining cycles are less populated and less meandering, while in the outer region, especially near the channel center, the velocity divergence is as strong as the vorticity. The temperature, density, and pressure fluctuations are enhanced by these waves. The correlations between the velocity and temperature are altered, due to the counter effects caused by the vortical motions and isentropic waves.

Wall shear stress and wall heat flux in a supersonic turbulent boundary layer

We report the characteristics of wall shear stress (WSS) and wall heat flux (WHF) from direct numerical simulation (DNS) of a spatially developing zero-pressure-gradient supersonic turbulent boundary layer at a free-stream Mach number M∞ = 2.25 and a Reynolds number Reτ = 769 with a cold-wall thermal condition (a ratio of wall temperature to recovery temperature Tw/Tr = 0.75). A comparative analysis is performed on statistical data, including fluctuation intensity, probability density function, frequency spectra, and space–time correlation. The root mean square fluctuations of the WHF exhibit a logarithmic dependence on Reτ similar to that for the WSS, the main difference being a larger constant. Unlike the WSS, the probability density function of the WHF does not follow a lognormal distribution. The results suggest that the WHF contains more energy in the higher frequencies and propagates downstream faster than the WSS. A detailed conditional analysis comparing the flow structures responsible for extreme positive and negative fluctuation events of the WSS and WHF is performed for the first time, to the best of our knowledge. The conditioned results for the WSS exhibit closer structural similarities with the incompressible DNS analysis documented by Pan and Kwon [“Extremely high wall-shear stress events in a turbulent boundary layer,” J. Phys.: Conf. Ser. 1001, 012004 (2018)] and Guerrero et al. [“Extreme wall shear stress events in turbulent pipe flows: Spatial characteristics of coherent motions,” J. Fluid Mech. 904, A18 (2020)]. Importantly, the conditionally averaged flow fields of the WHF exhibit a different mechanism, where the extreme positive and negative events are generated by a characteristic two-layer structure of temperature fluctuations under the action of a strong Q4 event or a pair of strong oblique vortices. Nevertheless, we use the bi-dimensional empirical decomposition method to split the fluctuating velocity and temperature structures into four different modes with specific spanwise length scales, and we quantify their influence on the mean WSS and WHF generation. It is shown that the mean WSS is mainly related to small-scale structures in the near-wall region, whereas the mean WHF is associated with the combined action of near-wall small-scale structures and large-scale structures in the logarithmic and outer regions.

Parameterised drag model for the underlying surface roughness of buildings in urban wind environment simulation

Knowledge of the flow field in the urban canopy layer (UCL) is of great significance in guiding the layout of buildings in blocks and mitigating local air pollution. Large eddy simulation (LES) with relatively high precision can simulate the flow field in the UCL. Considering the limitations of computational power, building resistance source parameterisation is an effective method of solving the average blocking effect of buildings on airflow and reducing the number of grids. Belcher's drag model has been widely used for calculating the average flow field in the UCL. However, the model does not consider the influence of the horizontal recirculation region and the variation of the canopy-drag length scale in the height direction. Therefore, a parameterised drag model for the resistance source of building(s) was improved by considering the horizontal recirculation region and effective roughness density. The simulation results suggest that the horizontal recirculation area of an isolated building is mainly affected by the spanwise width of the building at a certain height and presents a quadratic function distribution in the height direction, with a maximum at approximately z/H = 0.45. The validation of the isolated building and building arrays with different layouts demonstrated that the parameterised drag model could reveal the average flow field more accurately than Belcher's drag model.

Large eddy simulation on the flow characteristics of an argon thermal plasma jet

Large eddy simulations based on the CFD software OpenFOAM have been used to study the effect of Reynolds number and turbulence intensity on the flow and mixing characteristics of an argon thermal plasma jet. Detailed analysis was carried out with respect to four aspects: the average flow field, the instantaneous flow field, turbulence statistical characteristics and the self-similarity. It was shown that for the argon thermal plasma jet with low Reynolds number, increasing the turbulence intensity will increase the turbulent transport mechanism in the mixing layer rather than in the jet axis, leading to the faster development of turbulence. The effect of the turbulent transport mechanism increases with increasing Reynolds number. However, the characteristics of flow and mixing are not affected by turbulence intensity for high Reynolds number situations. It was also found that the mean axial velocity and mean temperature in the axis of the turbulent thermal plasma jet satisfy the self-similarity aspects downstream. In addition, decay constant K is 1.25, which is much smaller than that (5.7–6.1) of the turbulent cold gas jet and has nothing to do with the Reynolds number or turbulence intensity in the jet inlet.