Unsteady free surface wave-induced separation: Vortical structures and instabilities

Vortical structures and instability mechanisms of the unsteady free surface wave-induced separation around a surface-piercing NACA0024 foil at a Froude number of 0.37 and a Reynolds number of 1.52×10⁶ are studied using an unsteady Reynolds-averaged Navier–Stokes (URANS) code with a blended k?ε/k?ω turbulence model and a free surface tracking method. At the free surface, the separated flow reattaches to the foil surface resulting in a wall-bounded separation bubble. The mean and instantaneous flow topologies in the separation region are similar to the owl-face pattern. The initial shear-layer instability, the Karman-like instability, and the flapping instability are identified, and their scaling and physical mechanisms are studied. Validation with experimental fluid dynamics (EFD) and comparison with complementary detached-eddy simulation (DES) indicate that URANS resolves part of the organized oscillations due to the large-scale unsteady vortical structures and instabilities, thereby capturing the gross features of the unsteady separation. The URANS solutions show an initial amplitude defect of 30% for the free surface oscillations where the shear layer separates, and the defect progressively increases downstream as URANS rapidly dissipates the rolled up vortices.     

High‐order detached eddy simulation,zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth‐order discretization scheme is used and is validated using the test case of a two‐dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS–NLES approach does not suffer from this trait. For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton–Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS–NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. Copyright © 2013 John Wiley & Sons, Ltd.     

Detached eddy simulation of flow around two wall-mounted cubes in tandem

We investigate the performance of unsteady Reynolds-averaged Navier–Stokes (URANS) computation and various versions of detached eddy simulation (DES) in resolving coherent structures in turbulent flow around two cubes mounted in tandem on a flat plate at Reynolds number (Re) of 22,000 and for a thin incoming boundary layer. Calculations are carried out using four different coherent structure resolving turbulence models: (1) URANS with the Spalart–Allmaras model; (2) the standard DES [Spalart, P.R., Jou, W.H., Strelets, M., Allmaras, S.R., 1997. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In: Liu, C., Liu, Z., (Eds.), Advances in DNS/LES. Greyden Press, Columbus, OH]; (3) the Delayed DES (DDES); and (4) the DES with a low-Re modification (DES-LR) [Spalart, P., Deck, S., Shur, M., Squires, K., Strelets, M., Travin, A., 2006. A new version of detached eddy simulation, resistant to ambiguous grid densities. Theor. Comput. Fluid Dyn. 20 (3), 181–195]. The grid sensitivity of the computed solutions is examined by carrying out simulations on two successively refined grids. The computed results for all cases are compared with the experimental measurements of Martinuzzi and Havel [Martinuzzi, R., Havel, B., 2000. Turbulent flow around two interfering surface-mounted cubic obstacles in tandem arrangement. ASME J. Fluids Eng. 122, 24–31] for two different cube spacings. All turbulence models reproduce essentially identical separation of the approach thin boundary layer and yield an unsteady horseshoe vortex system consisting of multiple vortices in the leading edge region of the upstream cube. Significant discrepancies between the URANS and all DES solutions are observed, however, in other regions of interest such as the shear layers emanating from the cubes, the inter-cube gap and the downstream wake. Regardless of the grid refinement, URANS fails to capture key features of the mean flow, including the second horseshoe vortex in the upstream junction and recirculating flow on the top surface of the downstream cube for the large cube spacing, and underestimates significantly turbulence statistics in most regions of the flow for both cases. On the coarse mesh, all three DES approaches appear to yield very similar results and fail to reproduce the second horseshoe vortex. The standard DES and DDES solutions obtained on the fine meshes are essentially identical and both suffer from premature switching to unresolved DNS, due to the mis-interpretation of grid refinement as wall proximity, which leads to spurious vortices in the inter-cube region. Numerical solutions show that the low-Re modification (DES-LR) is critical prerequisite in DES on the ambiguously fine – not fine enough for full LES – mesh to prevent excessive nonlinear drop of the subgrid eddy viscosity in low cell-Re regions like in the inter-obstacle gap. Mean flow quantities and turbulence statistics obtained with DES-LR on the fine mesh are in good overall agreement with the measurements in most regions of interest for both cases.     

Hybrid URANS/LES Turbulence Simulation of Vortex Shedding Behind a Triangular Flameholder

In this study a detached eddy simulation (DES) model, which belongs to the group of hybrid URANS/LES turbulence models, is used for the simulation of vortex shedding behind a triangular obstacle. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS shear-stress transport (SST) model is used. In the other regions with higher grid resolution a LES model, which uses a transport equation for the turbulent subgrid energy, is applied. The DES model is first investigated for two standard test cases, namely decaying homogeneous isotropic turbulence and the backward facing step, respectively. For the decaying homogeneous isotropic turbulence test case the evolution of the energy spectra in wavenumber space for different times are studied for both the DES and a Smagorinsky type LES model. Different grid resolutions are analyzed with a special emphasis on the modeling constant connecting the filter length scale to the grid size. The results are compared to experimental data. The backward facing step test case is used to study the model behavior for a case with a transition region between a RANS modeling approach close to the wall and LES based modeling in the intense shear flow region. The final application is the simulation of the vortex shedding behind a triangular obstacle. First, the influence of the inlet condition formulation is studied in detail as they can have a significant influence especially for LES based models. Detailed comparisons between simulation and experiment for the flow structure past the obstacle and statistical quantities such as the shedding frequency are shown. Finally the additional temporal and spatial information provided by the DES model is used to show the predicted anisotropy of turbulence.     

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.     

Evaluation of linear and nonlinear convection schemes on multidimensional non‐orthogonal grids with applications to KVLCC2 tanker

Almost all evaluations of convection schemes reported in the literature are conducted using simple problems on uniform orthogonal grids; thus, having limited contribution when solving industrial computational fluid dynamics (CFD), where the grids are usually non‐orthogonal with distortions. Herein, several convection schemes are assessed in uniform and distorted non‐orthogonal grids with emphasis on industrial applications. Linear and nonlinear (TVD) convection schemes are assessed on analytical benchmarks in both uniform and distorted grids. To evaluate the performance of the schemes, four error metrics are used: dissipation, phase and L₁ errors, and the schemes' effective order of accuracy. Qualitative and quantitative deterioration of these error metrics as a function of the grid distortion metrics are investigated, and rigorous verifications are performed. Recommendations for effective use of the convection schemes based on the range of grid aspect ratio (AR), expansion ratio (ER) and skewness (Q) are included. A ship hydrodynamics case is studied, involving a Reynolds averaged Navier–Stokes simulation of a bare‐hull KVLCC2 tanker using linear and nonlinear convection schemes coupled with isotropic and anisotropic Reynolds‐stress (ARS) turbulence models using CFDShip‐Iowa v4. Predictions of local velocities and turbulent quantities from the midships to the nominal wake plane are compared with experimental fluid dynamics (EFD), and rigorous verification and validation analyses for integral forces and moments are performed for 0° and 12° drift angles. Best predictions are observed when coupling a second‐order TVD scheme with the anisotropic turbulence model. Further improvements are observed in terms of prediction of the vortical structures for 30° drift when using TVD2S‐ARS coupled with DES. Copyright © 2009 John Wiley & Sons, Ltd.     

Turbulence modelling of problem aerospace flows

Unsteady Reynolds averaged Navier–Stokes (URANS) and detached eddy simulation (DES) related approaches are considered for high angle of attack NACA0012 airfoil, wing–flap, generic tilt‐rotor airfoil and double‐delta geometry flows. These are all found to be problem flows for URANS models. For DES fifth‐order upwinding is found too dissipative and the use of, for high speed flows, instability prone centred differencing essential. An existing hybrid ILES–RANS modelling approach, intended for flexible geometry, relatively high numerical dissipation codes is tested along with differential wall distance algorithms. The former gives promising results. The standard turbulence modelling approaches are found to give perhaps a surprising results variation. Results suggest that for the problem flows, the explicit algebraic stress and Menter shear stress transport (SST) URANS models are more accurate than the economical Spalart–Allmaras (SA). However, the explicit algebraic stress model (EASM) in its k–ε form is impractically expensive to converge. Here, SA predictions lack a rotation correction term and this is likely to improve these results. Copyright © 2005 John Wiley & Sons, Ltd.     

URANS prediction of flow fluctuations in rod bundle with split-type spacer grid

Turbulent flow in a rod bundle with split-type spacer grid has been studied using Unsteady Reynolds-Averaged Navier–Stokes (URANS) approach. In the previous studies of turbulent flow in rod bundles URANS (as well as steady-state RANS) simulations predicted mean velocity profiles fairly well. However, they severely underpredicted velocity fluctuations, which is investigated in the present study. Our simulations were performed with the Shear Stress Transport (SST) turbulence model and automatic wall-treatment using OpenFOAM, an open-source CFD code. Results of URANS simulations are compared with the measurements of the MATiS-H experiment, which was performed at Korean Atomic Energy Research Institute (KAERI) in 2011–2012.The URANS predictions of velocity fluctuations have been improved by appropriately summing up fluctuations resolved by the basic URANS model and non-resolved fluctuations, which were modelled with the turbulence model. This treatment of turbulent fluctuations, which are directly measured in high-quality experiments, allows more detailed evaluation of various URANS turbulence models. It was found out that the best agreement is achieved when resolved and modelled fluctuations are assumed to be uncorrelated, which indicates that the large-scale structures in this particular flow are distinct in the spectral space from the rest of turbulence. Turbulent flow in the MATiS-H experiment was reproduced by numerous authors using different approaches and our results are among the most accurate.     

Unsteady CFD simulations of a pump in part load conditions using scale-adaptive simulation

The scope of this work is to demonstrate the applicability of an eddy resolving turbulence model in a turbomachinery configuration. The model combines the Large Eddy Simulation (LES) and the Reynolds Averaged Navier Stokes (RANS) approach. The point of interest of the present investigation is the unsteady rotating stall phenomenon occurring at low part load conditions. Since RANS turbulence models often fail to predict separation correctly, a LES like model is expected to give superior results. In this investigation the scale-adaptive simulation (SAS) model is used. This model avoids the grid dependence appearing in the Detached Eddy Simulation (DES) modelling strategy. The simulations are validated with transient measurement data. The present results demonstrate, that both models are able to predict the major stall frequency at part load. Results are similar for URANS and SAS, with advantages in predicting minor stall frequencies for the turbulence resolving model.     

Unsteady RANS and detached-eddy simulations of flow around a circular cylinder in ground effect

Unsteady Reynolds-averaged Navier–Stokes (URANS) simulations and detached-eddy simulations (DES) were performed of flow around a circular cylinder placed near and parallel to a moving ground, on which substantially no boundary layer developed to interfere with the cylinder. The results were compared with experiments previously reported by the authors to examine how accurately the URANS and DES can predict the cessation of von Kármán-type vortex shedding and the attendant critical drag reduction of the cylinder in ground effect. The DES, which were performed in a three-dimensional domain with spanwise periodicity imposed, correctly captured the cessation of the vortex shedding, whereas both two- and three-dimensional URANS also predicted it but at a much smaller gap-to-diameter ratio compared with the experiments. The wake structures of the cylinder predicted by the DES were in good agreement with the experiments in both large- and small-gap regimes, and also in the intermediate-gap regime, where the DES captured the intermittence of the vortex shedding in the near-wake region. Based on the results obtained, further discussions are also given to the reason why the von Kármán-type vortices in the URANS solutions incorrectly ‘survived’ until the cylinder came much closer to the ground.     

In-cylinder flows of a motored four-stroke engine with flat-crown and slightly concave-crown pistons

The temporal and spatial evolution processes of the in-cylinder flow structures and turbulence intensities in the symmetry and offset planes of a motored four-valve, four-stroke engine during the intake and compression strokes are diagnosed by using a particle image velocimeter. Two pistons of different crown shapes (flat-crown and slightly concave-crown pistons) are studied. The inception, establishment, and evolution of the tumbling vortical flow structures during the intake and compression strokes are clearly depicted. Quantitative strengths of the rotating vortical flow motions are presented by a dimensionless parameter, the tumble ratio, which can represent the mean angular velocity of the vortices in the target plane. The turbulence intensity of the in-cylinder flow is also calculated by using the measured time-varying velocity data. The results show that the flat-crown piston induces higher bulk-averaged tumble ratio and turbulence intensity than the slightly concave-crown piston does because the tumble ratio and turbulence generated by the flat-crown piston in the offset planes during the compression stroke are particularly large. The engine with the flat-crown piston also presents larger torque and power outputs and lower hydrocarbon emission than that with the slightly concave-crown piston. This might be caused by the enhanced combustion in the engine cylinder due to the stronger tumble ratio and turbulence intensity.     

Highly Resolved LES and URANS of Turbulent Buoyancy-Driven Flow Within Inclined Differentially-Heated Enclosures

In the present study, both LES and unsteady RANS computations are presented, for turbulent natural convection of air inside differentially-heated rectangular tilted cavities using a finite volume code (Code_Saturne), for an aspect ratio of H/L?=?28.6 and Rayleigh number of 0.86×10⁶. Attention is focused on two angles of inclination: 15° to the horizontal with hot lower and cold upper wall, the 15° unstable case, and the mirror image of this case where the angle is the same but with a hot upper and cold lower wall, the 15° stable case. In accordance with recent experimental data, the LES computations for both the stable and unstable tilted cavities returned three-dimensional time-averaged flow fields. In the case of the unstably stratified enclosure, the flow is highly unsteady with coherent turbulent structures in the core of the enclosure. Time-averaged temperature, velocity and resolved turbulence intensities resulting from LES computations show close agreement to measured data. Subsequent comparisons of different URANS schemes with LES are used in order to explore to what extent these models are able to reproduce the large-scale unsteady flow structures. All URANS schemes have been found to be able to reproduce the 3-D unsteady flow features present in the 15° unstable cavity. However, the low-Reynolds-number model tested, as well as requiring a high resolution near-wall grid, also needed a finer grid in the core region than the high-Reynolds-number models, thus making it computationally very expensive. Flow within the 15° stable cavity also shows some 3-D features, although it is significantly less unsteady, and the URANS models tested here have been less successful in reproducing this flow pattern. The overall heat transfer is presented here for both differentially heated enclosures.     

Turbulence modelling of the flow past a pitching NACA0012 airfoil at and Reynolds numbers

This paper provides a study of the NACA0012 dynamic stall at Reynolds numbers 10⁵ and 10⁶ by means of two- and three-dimensional numerical simulations. The turbulence effect on the dynamic stall is studied by statistical modelling. The results are compared with experiments concerning each test case. Standard URANS turbulence modelling have shown a quite dissipative character that attenuates the instabilities and the vortex structures related to the dynamic stall. The URANS approach Organised Eddy Simulation (OES) has shown an improved behaviour at the high Reynolds number range. Emphasis is given to the physical analysis of the three-dimensional dynamic stall structure, for which there exist few numerical results in the literature, as far as the Reynolds number range is concerned. This study has shown that the downstroke phases of the pitching motion are subjected to strong three-dimensional turbulence effects along the span, whereas the flow is practically two-dimensional during the upstroke motion.     

Non-zonal detached eddy simulation coupled with a steady RANS solver in the wall region

Xiao and Jenny (2012) proposed an interesting hybrid LES/RANS method in which they use two solvers and solve the RANS and LES equations in the entire computational domain. In the present work this method is simplified and used as a hybrid RANS-LES method, a wall-modeled LES. The two solvers are employed in the entire domain. Near the walls, the flow is governed by the steady RANS solver; drift terms are added to the DES equations to ensure that the time-averaged DES fields agree with the steady RANS field. Away from the walls, the flow is governed by the DES solver; in this region, the RANS field is set to the time-averaged LES field. The disadvantage of traditional DES models is that the RANS models in the near-wall region – which originally were developed and tuned for steady RANS – are used as URANS models where a large part of the turbulence is resolved. In the present method – where steady RANS is used in the near-wall region – the RANS turbulence models are used in a context for which they were developed. In standard DES methods, the near-wall accuracy can be degraded by the unsteady agitation coming from the LES region. It may in the present method be worth while to use an accurate, advanced RANS model. The EARSM model is used in the steady RANS solver. The new method is called NZ S-DES . It is found to substantially improve the predicting capability of the standard DES. A great advantage of the new model is that it is insensitive to the location of the RANS-LES interface.     

Numerical and experimental investigation on droplet dynamics and dispersion of a jet engine injector

In the case of turbine combustors operating with liquid fuel the combustion process is governed by the liquid fuel atomization and its dispersion in the combustion chamber. By highly unsteady flow field conditions the transient interaction between the liquid and the gaseous phase is of interest, because it results in a temporal variation of air–fuel ratio which leads to a fluctuating temperature distribution. The objective of this research was the investigation of transient flow field phenomena (e.g. large coherent structures) on droplet dynamics and dispersion of an isothermal flow (of inert water droplets) as a necessary first step towards a full analysis of spray combustion in real-life devices. The advanced injector system for lean jet engine combustors PERM (Partial Evaporated Rapid Mixing) was applied, generating a dilute polydispersed spray in a swirled flow field. Experiments were performed using Phase Doppler Anemometry (PDA) and a patternator to determine the droplet polydispersity, concentration maps, and velocity profiles in the flow. An important finding is the effect of large-scale coherent structures due mainly to the precessing of the vortex core (PVC) of the swirling air jet on the particle dispersion patterns. The experimental results then serve as reference data to assess the accuracy of the Eulerian–Lagrangian computations using a Large Eddy Simulation (LES), a Unsteady Reynolds-Average Navier–Stokes Simulation (URANS) and two simplified (steady-state) simulations. There, a simplified droplet injection model was used and the required boundary conditions of injected droplet sizes were obtained from measurements. Important transient effects of deterministic droplet separation observed during experiments, could be perfectly replicated with this injection model. It is convincingly shown, through extensive computations, that the resolution of instantaneous vortical structures is indeed crucial; hence the LES, or a reasonably-well resolved URANS are preferred over the steady-state solutions with additional, stochastic-type, turbulent dispersion models.     

Large-eddy simulation of counter-rotating Taylor–Couette flow: The effects of angular velocity and eccentricity

In the present work, turbulent flow in the annulus of a counter-rotating Taylor-Couette (CRTC) system is studied using large-eddy simulation. The numerical methodology employed is validated, for both the mean and second-order statistics, with the direct numerical simulation (DNS) data available in the literature, for a range of Reynolds numbers from 500 to 4000. Thereafter, turbulent flow occurring in this system at Reynolds numbers of 8000 and 16000 are studied, and the results obtained are analyzed using mean and second-order statistics, vortical structures, velocity vector plots and power energy spectra. Further, the spatio-temporal variation of azimuthal velocity, extracted near the inner cylinder, shows the existence of herringbone like patterns similar to that observed in the previous studies. The effect of eccentricity of the inner cylinder with respect to the outer cylinder is studied, on the turbulent flow in the CRTC system, for two different eccentricity ratios of 0.2 and 0.5 and for two different Reynolds numbers of 1500 and 4000. The results of the eccentric CRTC are analyzed using contours of pressure, mean and second-order statistics, velocity vectors, vortical structures, and turbulence anisotropy maps. It is observed from the eccentric CRTC simulations that the smaller-gap region seems to contain higher amplitude fluctuations and more vortical structures when compared with the larger-gap region. The mean turbulent kinetic energy contours do not change qualitatively with the Reynolds number, however, quantitatively a higher turbulent kinetic energy is observed in the higher Reynolds number case of 4000.     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

The Vortical Structures in the Rear Separation and Wake Produced by a Supersonic Micro-Ramp

The vortical structures in the rear separation and wake region produced by a micro-ramp that immersed in a supersonic turbulent boundary layer are investigated. The small scale separation close to the trailing edge was revealed and this confirms the previous experimental observation. Between the reverse region and surrounding fast moving flow, a three-dimensional shear layer was formed, and vortices are generated. By using vortex line method, the spiral points were illustrated as the cross-sections of the Ω-shaped vortices that follow the shape of the separation. The vortical structure was analogous to that in the wake region, where similar Ω-shaped vortex which follows the deficit region caused by the micro-ramp can be observed. Finally, the revealed flow topology was conceived new and beneficial to the studying of wall bounded turbulence which involves similar vortical structures but in a smaller scale, compared with the vortical pattern in the current micro-ramp wake.     

Numerical simulation of vortical and coherent structures in compressible jet flows

The compressible flows of plane free jets and jets of the intake-stroke of a rectangular piston-engine model are investigated by numerical simulations. The observed vortical structures appear to be the well-known coherent structures of turbulent shear layers. The simulated structures are compared to experimental data by means of density fields and turbulent statistics taken from different authors. The computed flow depends on physical as well as on numerical parameters. The good agreement with the experimental data is obtained by direct simulation without any turbulence model.