Three-dimensional Flow Structures Research Articles

A genuinely three-dimensional Riemann solver based on self-similar variables in curvilinear coordinates

A genuinely three-dimensional Riemann solver named MuSIC-3DC (Multidimensional, Self-similar, strongly-Interacting, Consistent, Three-dimensional, Curvilinear coordinates) in curvilinear coordinates is proposed to simulate three-dimensional complex engineering problems. Following Balsara's idea, a three-dimensional Riemann problem at each grid vertex is studied. In order to accurately capture contact discontinuities, the linear variations are assumed in the strongly-interacting zone. Also, the self-similar variables are adopted to simplify the derivation according to the self-similar evolution of the Riemann problem, and the strongly-interacting state could be determined by combining the contributions of all the two-dimensional Riemann problems at the surfaces bounded the strongly-interacting state. Besides, the one-dimensional Riemann problem at the center of each cell interface is solved by the one-dimensional HLLE scheme. After these, the numerical flux at each cell interface is obtained by assembling the three-dimensional fluxes at the vertices and the one-dimensional flux at the center of the cell interface. Several numerical test cases are simulated to validate the proposed solver. The spherical blast wave and three-dimensional Riemann problems indicate that the MuSIC-3DC scheme is capable of capturing three-dimensional self-similar flow structures accurately and suppressing the mesh imprinting phenomenon effectively. The transonic case demonstrates the MuSIC-3DC scheme's high resolution in capturing shock waves. In the hypersonic cases, the MuSIC-3DC scheme is more accurate in predicting stagnation and reattachment heating loads. Moreover, the MuSIC-3DC scheme is capable of obtaining more physical surface heating distribution by considering the information traveling both normal and transverse to the cell interfaces.

Multi-Scale Flow Estimation Within Stereoscopic PIV Processing Based On Deep Learning

Classical methods of particle-image velocimetry (PIV) based on state-of-the-art cross-correlation algorithms constitute a spatial filtering operation that leads to a loss of spatial resolution of the estimated velocity fields. This is associated with numerous limitations, e.g., when it comes to aligning high-fidelity numerical simulations, since the detection of small turbulent scales is often precluded. An optical flow neural network adapted for PIV processing, called RAFT-PIV, can overcome the drawback of reduced spatial resolution while maintaining the reliability and robustness of the classical approach. The present work aims to extend the application of the aforementioned deep learning based method by integrating it within a stereoscopic PIV processing workflow. Experimental data from an optically accessible single-cylinder combustion engine the flow field of which is characterized by complex, three-dimensional flow structures and a broad range turbulent scales are processed using RAFT-PIV. The 2D-3C vector fields obtained from stereoscopic reconstruction are compared against a classical cross-correlation based PIV evaluation, providing insight into instantaneous quantities and turbulent structures. A qualitative comparison of the instantaneous in-plane and out-of-plane velocity components reveals that the reconstructed flow obtained from the deep learning based estimation includes a wide range of resolved turbulent length scales. It outperforms the classical method. Furthermore, the findings of the present study demonstrate that the pixel accuracy has an effect on both physical quantities and turbulence characteristics. In particular, an increase in turbulent kinetic energy is observed across the engine cycles of the investigated engine flow states. Moreover, a higher fraction of three-dimensional isotropic turbulent structures is derived from RAFT-PIV compared to the results from the cross-correlation based technique.

Flow dynamics and turbulent coherent structures around sediment reduction plates of a sewer system

In the management of urban drainage networks, great interest has been generated in the removal of sediments from sewer systems. The unsteady three-dimensional (3D) flow and turbulent coherent structures surrounding sediment reduction plates in a sewer system are investigated by means of the detached-eddy simulation (DES). Particular emphasis is given to detailing the instantaneous velocity and vorticity fields within the grooves, along with an examination of the three-dimensional, long-term, average flow structure at a Reynolds number of approximately 105. Velocity vectors demonstrate continuous flapping of the flow on the groove wall, periodically interacting with ejections of positive and negative vorticity originating from the grooves. The interaction between the three-dimensional groove flow and the shear flow leads to the downstream transport of patches of positive and negative vorticity, which significantly influence sediment transport. The high-velocity shear flows and strong vortices generated in undulating topography, as identified by the Q-criteria, are the key factors contributing to the efficient sediment reduction capabilities of the sediment reduction plates. The sediment reduction plates with partially enclosed structures exhibit low sedimentation rates in grooves on the plate, a broader acceleration region, and a lesser impact on the flow capacity. The results improve the understanding of the hydrodynamics and turbulent coherent structures surrounding the sediment reduction plates while elucidating the driving factors behind the enhancement of sediment scouring and suspension capacities. These results indicate that the redesign of the plates as partially enclosed structures contributes to further improving their sediment reduction performance.

Experimental and numerical investigation of the effects of porous bleed systems on a model scramjet inlet under high backpressure conditions

In this study, we conducted experimental and numerical investigations to assess the effects of backpressure and configurations of porous bleeds on a scramjet inlet's performance at Mach 5. Experiments in a high-enthalpy wind tunnel, using a backpressure controller and a porous bleed with 1 mm holes and 10 mm lateral spacing, were paired with Reynolds-averaged Navier-Stokes simulations to explore the mitigation of flow separation and the prevention of inlet unstart under various backpressure conditions. The simulations provided insights into how backpressure levels and porous bleed geometries impact internal flow structures, as well as the critical backpressure ratio at which inlet unstart occurs. Additional simulations varied the bleed system's hole diameter and spacing, identifying configurations that optimize aerodynamic performance by enhancing total pressure recovery and Mach number at the isolator exit, while reducing bleed mass flow rates. The study also quantified the impact of these configurations on engine thrust, determining that certain bleed system designs significantly improve thrust by efficiently suppressing the interaction between the core flow and corner recirculation zones. This study examined the internal three-dimensional flow structure characteristics under high back pressure and provided insights into porous bleed configuration capable of efficiently suppressing inlet unstart.

Combining Computational Fluid Dynamics and Experimental Data to Understand Fish Schooling Behavior.

Understanding the flow physics behind fish schooling poses significant challenges due to the difficulties in directly measuring hydrodynamic performance and the three-dimensional, chaotic, and complex flow structures generated by collective moving organisms. Numerous previous simulations and experiments have utilized computational, mechanical, or robotic models to represent live fish. And existing studies of live fish schools have contributed significantly to dissecting the complexities of fish schooling. But the scarcity of combined approaches that include both computational and experimental studies, ideally of the same fish schools, has limited our ability to understand the physical factors that are involved in fish collective behavior. This underscores the necessity of developing new approaches to working directly with live fish schools. An integrated method that combines experiments on live fish schools with computational fluid dynamic (CFD) simulations represents an innovative method of studying the hydrodynamics of fish schooling. CFD techniques can deliver accurate performance measurements and high-fidelity flow characteristics for comprehensive analysis. Concurrently, experimental approaches can capture the precise locomotor kinematics of fish and offer additional flow information through particle image velocimetry (PIV) measurements, potentially enhancing the accuracy and efficiency of CFD studies via advanced data assimilation techniques. The flow patterns observed in PIV experiments with fish schools and the complex hydrodynamic interactions revealed by integrated analyses highlight the complexity of fish schooling, prompting a reevaluation of the classic Weihs model of school dynamics. The synergy between CFD models and experimental data grants us comprehensive insights into the flow dynamics of fish schools, facilitating the evaluation of their functional significance and enabling comparative studies of schooling behavior. In addition, we consider the challenges in developing integrated analytical methods and suggest promising directions for future research.

Three-dimensional flow structure in a confluence-bifurcation unit

Enhanced understanding of flow structure in braided rivers is essential for river regulation, flood control, and infrastructure safety across the river. It has been revealed that the basic morphological element of braided rivers is confluence-bifurcation units. However, flow structure in these units has so far remained poorly understood with previous studies having focused mainly on single confluences/bifurcations. Here, the flow structure in a laboratory-scale confluence-bifurcation unit is numerically investigated based on the FLOW-3D® software platform. Two discharges are considered, with the central bars submerged or exposed respectively when the discharge is high or low. The results show that flow convergence and divergence in the confluence-bifurcation unit are relatively weak when the central bars are submerged. Based on comparisons with a single confluence/bifurcation, it is found that the effects of the upstream central bar on the flow structure in the confluence-bifurcation unit reign over those of the downstream central bar. Concurrently, the high-velocity zone in the confluence-bifurcation unit is less concentrated than that in a single confluence while being more concentrated than that observed in a single bifurcation. The present work unravels the flow structure in a confluence-bifurcation unit and provides a unique basis for further investigating morphodynamics in braided rivers. Highlights 3D flow structure in a confluence-bifurcation unit (CBU) is numerically investigated. Flow convergence/divergence in CBU is relatively weak when central bars are submerged. Effects of the upstream central bar on CBU flow reign over those of the downstream central bar. The high-velocity zone is less/more concentrated in the CBU than in a single confluence/bifurcation when the central bars are exposed.

Corner Vortex Structures Within Recirculation Zones of Confined Bluff-Body Flows

The vortical structures of recirculation zones in turbulent nonreacting and premixed reacting flows around confined equilateral-triangle bluff bodies are investigated using large-eddy simulations. A three-dimensional flow structure is observed within the recirculation zones of reacting and nonreacting cases; the structure includes both the “spanwise” vortices that extend the length of the recirculation zone and “corner vortex structures” that are situated adjacent to the bluff-body trailing edge and spanwise walls. The corner vortex structures enhance the mixing and residence times of fluid inside the recirculation zone. Fluid is circulated between the spanwise vortices and corner vortex structures. Corner vortex structures always appear downstream (relative to the local flow) of the boundary-layer separation location on the spanwise walls within the recirculation zone, and they are absent when boundary-layer formation on the spanwise walls is prevented. Laminar boundary-layer theory predicts the boundary-layer separation locations within ∼12% of the calculated values from the large-eddy simulations in all but one case. These results suggest that 1) boundary-layer separation on the spanwise walls is a necessary condition for the formation of corner vortex structures, and 2) boundary-layer separation on the spanwise walls is caused by an adverse pressure gradient in the recirculation zone.

An Investigation of Contaminant Transport and Retention from Storage Zone in Meandering Channels

Contaminant trapping by recirculation zones occurring at the apex of natural meandering channels induces a long tail in the contaminant cloud, thereby complicating the prediction of mixing behaviors. Thus, the understanding of the interaction between solute trapping and recirculating flow is important for responding to and mitigating water pollution accidents. In this research, the EFDC model was employed to reproduce three-dimensional flow structures of recirculating flow at the channel apex and investigate the influence on contaminant mixing. To investigate the contaminant transport characteristics from the storage zone in meandering channels, simulations were conducted using various discharge values to assess the impact of storage zone development on the concentration–time curves. The analysis of the relationship between the storage zone size and mixing behaviors indicates that an increase in discharge could result in a shorter tail and larger longitudinal dispersion even with the larger storage zone size. On the other hand, the enlarged recirculation zone size contributes to reducing transverse dispersion, evidenced by flatter dosage curves under lower flow rate conditions. These findings suggest that the increase in longitudinal dispersion with a larger flow rate is primarily caused by the reduction in transverse dispersion resulting from the formation of the recirculation zone.

Flow three-dimensionality of wavy elliptic cylinder: vortex shedding bifurcation

The three-dimensional flow of a wavy elliptic cylinder and its wake are investigated for wavelength (λ) to hydraulic diameter (Dm) ratio λ/Dm = 3.43, 4.58, and 6.01. The study aims to gain insights into flow separation, flow dynamics, and three-dimensional flow topologies for the selected wavelengths. Three-dimensional flow structures are visualized through time-averaged streamlines with critical point theory. The results for λ/Dm = 3.43 show spiral flow directed from the node to the saddle, counter-rotating vortices, and a wavy vortex structure. On the contrary, the flows for λ/Dm = 4.58 and 6.01 respectively undergo bifurcation in stable recirculation and vortex shedding on a half-span of the cylinder wavelength. The bifurcation results in a spiral flow toward the saddle and the other spiral flow toward the node. The bifurcations of recirculation and vortex shedding are revealed for the first time. The findings contribute to an improved understanding of the insight into intricate flow physics around wavy elliptic cylinders.

Tip effects on three-dimensional flow structures over low-aspect-ratio plates: mechanisms of spanwise fluid transport

Three-dimensional flows over low-aspect-ratio rectangular flat plates ( ${A{\kern-4pt}R} = 1.00$ – $1.50$ ) are investigated using tomographic and planar particle image velocimetry techniques. The chord-based Reynolds number is $5400$ , and the angle of attack is fixed at $6^\circ$ . This study reveals for the first time the interplay between spanwise fluid transport and downwash, both originating from the tip effects. Spanwise fluid transport promotes the formation and subsequent coherent development of leading-edge vortices, whereas downwash stabilizes the flow. Specifically, two mechanisms related to spanwise fluid transport are revealed. First, the spanwise fluid transport enhances the intensity of the reversed flow, promoting the shear layer roll-up and vortex shedding. Second, the near-wall spanwise flow interacts with the shed C-shape vortices, thereby strengthening the vortex heads. In particular, through these interactions, spanwise fluid transport can sustain the coherence of the C-shape vortices until the vortex heads split in a regular fashion. Consequently, the C-shape vortices are transformed into novel Þ-shape vortices for the plates of ${A{\kern-4pt}R} \leq 1.25$ , which supplements the previously discovered transformation from C-shape to M-shape vortices for larger ${A{\kern-4pt}R}$ plates. Downstream of this novel vortex-splitting transformation, two fundamental processes contribute to the formation of hairpin vortices. The above comprehensive understanding of complete vortex evolution routine provides valuable insights into the tip effects on the formation of three-dimensional flows over low- ${A{\kern-4pt}R}$ plates.

Revisiting the surface-mounted cube: An updated perspective of the near wake and near-wall flow field

The flow around a surface-mounted cube has been investigated for decades, yet the fundamentals of the mean flow have not been updated to reflect the latest advances regarding the flow around surface-mounted finite-height square prisms in general. One of the main gaps is the flow field very near the cube walls, especially the sides, and its relationship to the near wake. To investigate these features, large-eddy simulations of the flow around a surface-mounted cube were carried out at a Reynolds number Re =1×104. Two boundary layers were considered: a thin and laminar boundary layer and a thick and turbulent one, to provide an overview of the flow around the cube for contrasting boundary layers. The major flow structures in the mean wake were the horseshoe vortex, the arch vortex and the dipole structures, with other regions of vorticity also present. The time-averaged flow fields presented similar flow features for both boundary layers, but the thicker and turbulent one caused the horseshoe vortex to be located closer to the cube, the wake to become narrower and shorter, the coherent structures and downwash to weaken, the pressure to decrease around the sides and top of the cube, the drag force coefficient to decrease, the normal force coefficient to increase, and the trailing edge vortices to weaken. Intermittent flow reattachment on the top and side faces of the cube, as well as corner vortices, were found for both boundary layers, while the side and free end vortices were found exclusively with the thicker boundary layer, yielding a headband vortex. The three-dimensional flow structures and features were related to the near-wall flow field on the cube and ground plane surfaces, giving a complete and updated description of the mean flow characteristics.

Flow-induced oscillations of pitching swept wings: stability boundary, vortex dynamics and force partitioning

We study experimentally the aeroelastic instability boundaries and three-dimensional vortex dynamics of pitching swept wings, with the sweep angle ranging from 0 $^\circ$ to 25 $^\circ$ . The structural dynamics of the wings are simulated using a cyber-physical control system. With a constant flow speed, a prescribed high inertia and a small structural damping, we show that the system undergoes a subcritical Hopf bifurcation to large-amplitude limit-cycle oscillations (LCOs) for all the sweep angles. The onset of LCOs depends largely on the static characteristics of the wing. The saddle-node point is found to change non-monotonically with the sweep angle, which we attribute to the non-monotonic power transfer between the ambient fluid and the elastic mount. An optimal sweep angle is observed to enhance the power extraction performance and thus promote LCOs and destabilize the aeroelastic system. The frequency response of the system reveals a structural-hydrodynamic oscillation mode for wings with relatively high sweep angles. Force, moment and three-dimensional flow structures measured using multi-layer stereoscopic particle image velocimetry are analysed to explain the differences in power extraction for different swept wings. Finally, we employ a physics-based force and moment partitioning method to correlate quantitatively the three-dimensional vortex dynamics with the resultant unsteady aerodynamic moment.

Numerical Simulation of 3D Flow Structure and Turbulence Characteristics near Permeable Spur Dike in Channels with Varying Sinuosities

Owing to the different degrees of bending in rivers in nature, it is difficult to conduct experiments in situ. In this study, the renormalization group (RNG) k-ε turbulence model in ANSYS Fluent was used to analyze the three-dimensional flow structure and turbulence characteristics near a spur dike and to evaluate the variation trend of flow in rivers with different degrees of bending. The results show that in channels with different curvatures, the vortex appears between the spur dikes and is disturbed by the permeable hole, and the backflow area moves downstream. The strength of secondary flow (SSF) fluctuates greatly in the vicinity of the spur dike and the downstream region, and the peak value appears 3.22 m (21.5 times L) away from the inlet of the bend. The SSF increases as the bend curvature increases. The SSF displays similar variation trends in the three kinds of bends. The peak value of normalized turbulent kinetic energy (NTKE) appears 3.14 m away from the entrance of the bend, the NTKE is the largest in the 45° bend and the smallest in the 180° bend, and it decreases only at distances of 3.25–4.19 m away from the entrance of the bend as the bend curvature increases.

Three-dimensional flow structures in turbulent Rayleigh–Bénard convection at low Prandtl number Pr = 0.03

In this paper we report on an experimental study focusing on the manifestation and dynamics of the large-scale circulation (LSC) in turbulent liquid metal convection. The experiments are performed inside a cylinder of aspect ratio $\varGamma = 0.5$ filled with the ternary alloy GaInSn, which has a Prandtl number of $Pr = 0.03$ . The large-scale flow structures are classified and characterized at Rayleigh numbers of ${Ra} = 9.33 \times 10^6, 5.31 \times 10^7$ and $6.02 \times 10^8$ by means of the contactless inductive flow tomography which enables the full reconstruction of the three-dimensional (3-D) flow structures in the entire convection cell. This is complemented with the multi-thermal-probe method for capturing the azimuthal temperature variation induced by the LSC at the sidewall. We use proper orthogonal decomposition (POD) to identify the dominating modes of the turbulent convection. The analysis reveals that a single-roll structure of the LSC alternates in short succession with double-roll structures or a three-roll structure. This is accompanied by dramatic fluctuations of the Reynolds number, whose instantaneous values can deviate by more than 50 % from the time-average value. No coherent oscillations are observed, whereas a correlation analysis indicates a residual contribution of the torsion and sloshing modes. Results of the POD analysis suggest a stabilization of the single-roll LSC with increasing $Ra$ at the expense of flow structures with multiple rolls. Moreover, the relative lifetime of all identified flow states, measured in units of free-fall times, increases with rising $Ra$ .

Drag reduction effect of distributed roughness on the transitional flow state using direct numerical simulation

We analyzed the details of Tollmien–Schlichting wave transitions on a sandy rough surface through direct numerical simulations by resolving each roughness level using the volume penalization and zonal methods. Although wall roughness has generally been discussed in terms of increased drag in fully developed turbulent flows, we found that the sandy roughness on a flat plate suppresses the growth of both the turbulent kinetic energy and friction drag coefficient during the laminar-turbulent transitional flow state. The roughness size is on the microscale order, assuming similarity based on the Reynolds number and considering the case of the surface of a transonic airplane. Therefore, we call such three-dimensional and complex-shaped roughness with this effect distributed microroughness (DMR). On the sandy rough surface, the Tollmien–Schlichting wave collapses into a three-dimensional flow structure in the initial stage; however, the turbulent kinetic energy remains small. Generally, the turbulent kinetic energy of the subsequent three-dimensional turbulent component increases in correlation with the size of the primary T-S vortex. On the rough surface, the T-S roll vortex collapses before it maximally grows due to the spanwise disturbance caused by the three-dimensionality of the sandy roughness shape, resulting in the three-dimensional turbulence associated with the roll remaining small. Near the wall surface, the velocity fluctuates due to the distributed roughness, and the statistical peaks of the kinetic energy budget tend to be away from a wall. Although the details of the mechanism by which changes in the energy balance cause a decrease in the productive turbulent energy have not been elucidated, we will clarify the changes due to roughness parameters in a concurrent study.

A Combined Delayed Detached Eddy Simulation and Linearized Navier–Stokes Equation Study on the Generation and Reduction of Aerodynamic Noises Inside Steam Turbine Control Valve With Acoustic Liner

Abstract This study was aimed at numerically investigating the source, generation mechanism, and strategy for reducing aerodynamic noises inside a steam turbine control valve. A delayed detached eddy simulation was performed to extract the three-dimensional unsteady turbulent flow structures formed within the serpentine flow passage of the turbine valve. Acoustic analogies, spatial Fourier transform, and spectral proper orthogonal decomposition on the delayed detached eddy simulation-simulated flow data were complementarily combined to clarify the generation mechanism of tonal and broadband aerodynamic noises. The results showed that broadband noises were produced by wall-attached jet flow and turbulent mixing flow between the annular wall jets and central reverse flow. High-intensity tonal noises were generated by the excitation of multi-order natural acoustic modes of the bell-shaped valve spindle. The intensive acoustic pressure pulsations concentrated inside the bell jar and propagated along the diffuser to the downstream turbine chamber. A novel ring acoustic liner was designed using the acoustic impedance model to reduce the valve noises without sacrificing the flow performance. The noise reduction effectiveness was evaluated by solving the linearized Navier–Stokes equations in the frequency domain.

In silico model to assess thrombosis risk of TAVR with hemodynamic predictors using fluid structure interaction

Abstract The promising results of transcatheter aortic valve replacement (TAVR) over the past two decades indicate an expansion of the patient cohort toward patients with intermediate or low surgical risk. Since some complications of TAVR have already been minimized, subclinical leaflet thrombosis (SLT) has gained importance in recent years. SLT is manifested by a thrombotic layer on the prosthetic leaflets that gradually reduces leaflet motion. The resulting decrease in functionality of the TAVR causes a need for re-intervention. The origin of SLT and approaches to prevent SLT are still unexplored. For this reason, we have developed an in silicomodel that can be used during the design development process of TAVR devices to estimate the thrombosis risk of the implant. Based on passive scalar transport, hemodynamic metrics are used to quantify platelet activation and aggregation which are associated with the formation of thrombosis. In conjunction with a numerical simulation model considering the fluid-structure interaction between the blood mimicking fluid and the TAVR implanted in an aortic root, the thrombosis risk can be modeled. The simulation model can be used to calculate the three-dimensional flow structures within the native sinus and neo-sinus and also provides the ability to derive metrics to assess the risk of thrombosis. We demonstrated that this in silicomodel is a time-effective tool to assess thrombosis risk in TAVR product development.

A random 2D walk of a submerged free-floating disc in a convective layer

Free motions of a thermally insulating disc submerged in a liquid layer heated from below are studied experimentally. The disc is fixed at a certain height from the bottom of the convective cell, but can move freely in horizontal directions. Blocking of vertical motion and heat flux leads to the appearance of large-scale vortices, which in turn, drags the disc. The interplay of the large-scale circulation and the disc produces complex two-dimensional disc motions. The disc dynamics is mainly determined by large-scale modes: toroidal vortex with descending fluid in the centre and large-scale vortices along each coordinate, existing on the background of intense, but relatively short-lived small-scale convective cells. The wandering disc forms a cold spot above itself, which affects large-scale circulation but has little effect on cell convection at scales on the order of the layer thickness, which provide the main heat transport and dominate the upper part of the layer. Despite the essentially three-dimensional flow structure, which is fundamentally different from that of the convective flow in an elongated tank with one-dimensional disc motions, the travels along each coordinate look similar to one-dimensional motions.

Control of Cowl Shock/Boundary Layer Interaction in Supersonic Inlet Based on Dynamic Vortex Generator

A shock wave/boundary layer interaction (SWBLI) is a common phenomenon in supersonic inlet flow, which can significantly degrade the aerodynamic performance of the inlet by inducing boundary layer separation. To address this issue, in this paper, we propose the use of a dynamic vortex generator to control the SWBLI in a typical supersonic inlet. The unsteady simulation method based on dynamic grid technology was employed to verify the effectiveness of the proposed method of control and investigate its mechanism. The results showed that, in a duct of finite width at the inlet, the SWBLI generated complex three-dimensional (3D) flow structures with remarkable swirling properties. At the same time, vortex pairs were generated close to the side wall as a result of its presence, and this led to the intensification of transverse flow and, in turn, the formation of a complex 3D structure of the flow of the separation bubble. The dynamic vortex generator induced oscillations of variable intensity in the vortex system in the supersonic boundary layer that enhanced the mixing between the boundary layer flow and the mainstream. Meanwhile, the unique effects of “extrusion” and “suction” in the oscillation process continued to charge the airflow, and the distribution of velocity in the boundary layer significantly improved. As the oscillation frequency of the vortex generator increased, its charging effect on low-velocity flow in the boundary layer increased, and its control effect on the flow field of the SWBLI became more pronounced. The proposed method of control reduced the length of the separation bubble by 31.76% and increased the total pressure recovery coefficient at the inlet by 6.4% compared to the values in the absence of control.

Three-dimensional flow structures and scalar mixing of a sand dune-inspired jet in crossflow

This paper performed particle image velocimetry (PIV) measurements and delayed detached-eddy simulations (DDES) to comprehensively investigate the flow and mixing behaviors of a sand dune(SD)-inspired JICF (jet in crossflow). PIV measurements were conducted in planes parallel to the mid-plane (Z/D = 0, 0.5, and 1), and three velocity ratios (VR = 0.4, 0.8, and 1.2) were studied. The time-averaged and turbulent flow features behind the sand dune-inspired JICF is strongly three-dimensional: the jet tends to move from the mid-plane to the outer planes and can stay attached well to the wall even at Z/D = 1. Then DDES method was applied to reveal detailed three-dimensional flow structures and mixing characteristics of the SD-inspired JICF, and a cylindrical JICF was also simulated as a benchmark for comparison. The horseshoe vortex, arch-shaped shear layer vortices, and streamwise rollers are mainly responsible for the formation of the anti-counter rotating vortex pair (anti-CRVP) observed in the mean flow field of SD configuration. Analysis of concentration distributions indicates that scalar mixing between the jet and the mainstream is highly correlated to the flow features and vortex structures, and the mixing strength of the sand dune-inspired JICF is generally weaker than that of the cylindrical JICF.