Modelling of pressure-strain correlation in compressible turbulent flow

Previous studies carried out in the early 1990s conjectured that the main compressible effects could be associated with the dilatational effects of velocity fluctuation. Later, it was shown that the main compressibility effect came from the reduced pressure-strain term due to reduced pressure fluctuations. Although better understanding of the compressible turbulence is generally achieved with the increased DNS and experimental research effort, there are still some discrepancies among these recent findings. Analysis of the DNS and experimental data suggests that some of the discrepancies are apparent if the compressible effect is related to the turbulent Mach number, Mt. From the comparison of two classes of compressible flow, homogenous shear flow and inhomogeneous shear flow （mixing layer）, we found that the effect of compressibility on both classes of shear flow can be characterized in three categories corresponding to three regions of turbulent Mach numbers： the low-Mr, the moderate-Mr and high-Mr regions. In these three regions the effect of compressibility on the growth rate of the turbulent mixing layer thickness is rather different. A simple approach to the reduced pressure-strain effect may not necessarily reduce the mixing-layer growth rate, and may even cause an increase in the growth rate. The present work develops a new second-moment model for the compressible turbulence through the introduction of some blending functions of Mt to account for the compressibility effects on the flow. The model has been successfully applied to the compressible mixing layers.

Optimizing calculation of phase screen distribution with minimum condition along an inhomogeneous turbulent path

When building an experimental platform for light propagation along an inhomogeneous turbulent path,it is very essential to set up the reasonable distribution of phase screen.Based on multi-layered model of phase screen,an iterative optimization algorithm of phase screen position is given in this paper.Thereafter,the optimal position of phase screens is calculated under the Hufnagel-Valley5/7 and Hefei-day turbulence profile.The results show that the positions of phase screen calculated by the iterative algorithm can fit well with the turbulence profile rather than mechanically placed phase screens at equal distance.Compared with the uniform distribution of phase screens position,the residual phase error of the iterative algorithm decreases very significantly.The similarity degree between them is minimal when number of layers is equal to two.

Turbulence mitigation scheme based on multiple-user detection in an orbital-angular-momentum multiplexed system

Atmospheric turbulence（AT） induced crosstalk can significantly impair the performance of a free-space optical（FSO）communication link using orbital angular momentum（OAM） multiplexing.In this paper,we propose a multiple-user detection（MUD） turbulence mitigation scheme in an OAM-multiplexed FSO communication link.First,we present a MUD equivalent communication model for an OAM-multiplexed FSO communication link under AT.In the equivalent model,each input bit stream represents one user＇s information.The deformed OAM spatial modes caused by AT,instead of the pure OAM spatial modes,are used as information carriers,and the overlapping between the deformed OAM spatial modes are computed as the correlation coefficients between the users.Then,we present a turbulence mitigation scheme based on MUD idea to enhance AT tolerance of the OAM-multiplexed FSO communication link.In the proposed scheme,the crosstalk caused by AT is used as a useful component to deduce users＇ information.The numerical results show that the performance of the OAM-multiplexed communication link has greatly improved by the proposed scheme.When the turbulence strength C_n^2 is 1 × 10^（-15） m^（-2/3）,the transmission distance is 1000 m and the channel signal-to-noise ratio（SNR）is 26 dB,the bit-error-rate（BER） performance of four spatial multiplexed OAM modes l_m = ＋ 1,＋2,＋3,＋4 are all close to 10-5,and there is a 2-3 fold increase in the BER performance in comparison with those results without the proposed scheme.In addition,the proposed scheme is more effective for an OAM-multiplexed FSO communication link with a larger OAM mode topological charge interval.The proposed scheme is a promising direction for compensating the interference caused by AT in the OAM-multiplexed FSO communication link.

Effect of inertial particles with different specific heat capacities on heat transfer in particle-laden turbulent flow

The effect of inertial particles with different specific heat on heat transfer in particle-laden turbulent channel flows is studied using the direct numerical simulation（DNS） and the Lagrangian particle tracking method. The simulation uses a two-way coupling model to consider the momentum and thermal interactions between the particles and turbulence. The study shows that the temperature fields display differences between the particle-laden flow with different specific heat particles and the particle-free flow,indicating that the particle specific heat is an important factor that affects the heat transfer process in a particle-laden flow. It is found that the heat transfer capacity of the particle-laden flow gradually increases with the increase of the particle specific heat. This is due to the positive contribution of the particle increase to the heat transfer. In addition,the Nusselt number of a particle-laden flow is compared with that of a particle-free flow.It is found that particles with a large specific heat strengthen heat transfer of turbulent flow, while those with small specific heat weaken heat transfer of turbulent flow.

“Equilibrium”and“non-equilibrium”turbulence

We review the concept of ‘‘equilibrium＇＇ in turbulence. It generally means a property of the energy spectrum, it can also be understood in terms of a scalar property, the Taylor–Kolmogorov formula relating the dissipation rate to the total energy and integral length scale. The implications of equilibrium and strong departure from equilibrium for turbulence modeling are stressed.

Response of turbulent enstrophy to sudden implementation of spanwise wall oscillation in channel flow

The response of turbulent enstrophy to a sudden implementation of spanwise wall oscillation（SWO） is studied in a turbulent channel flow via direct numerical simulation. In the beginning of the application of SWO, a significant correlation is formed between ω′yand ω′z. A transient growth of turbulent enstrophy occurs, which directly enhances turbulent dissipation and drifts the turbulent flow towards a new lower-drag condition. Afterwards, the terms related to the stretching of vorticity（ωx, ω′y, and ωz）,the inclination of ω′yby ？w/？y, the turn of z by ？v′/？z, and the horizontal shear of z by ？w′/？x are suppressed due to the presence of SWO, leading to attenuation of the turbulent enstrophy.

Effect of the Arrangement of a New-Type of Turbulator Inserts on Heat Pipe Exchanger Performances

This research tests the effect of introducing turbulators of a new type into a circular tube heat exchanger under a constant and uniform longitudinal heatflux condition.A 45 mm diameter copper tube with a length of 1,350 mm is utilized with a solid disk being inserted inside the tube,which consists of three sections,each one containing two slots.The slot is cut at a 45 degree angle toward the inner tube surface,which results in diverging theflow toward the inner hot tube surface in order to enhance the heat transfer process.Air is considered as the workingfluid with Prandtl number 0.71.The Reynolds number spans the interval from 6,000–13,500,which indicates that the consideredflow is turbulent.The heat exchanger performance is studied and analyzed in terms of average Nusselt number.The experimental results show that the Nusselt number value is directly proportional to the increase of the Reynolds number,and the number of turbulators inserts.With the use of three novel turbulators,the heat transfer was about 3.15 times higher than that in the smooth tube and the friction factor was about 1.11.

Dynamics of cavitation–structure interaction

Cavitation–structure interaction has become one of the major issues for most engineering applications. The present work reviews recent progress made toward developing experimental and numerical investigation for unsteady turbulent cavitating flow and cavitation–structure interaction. The goal of our overall efforts is to（1） summarize the progress made in the experimental and numerical modeling and approaches for unsteady cavitating flow and cavitation–structure interaction,（2） discuss the global multiphase structures for different cavitation regimes, with special emphasis on the unsteady development of cloud cavitation and corresponding cavitating flow-induced vibrations,with a high-speed visualization system and a structural vibration measurement system, as well as a simultaneous sampling system,（3） improve the understanding of the hydroelastic response in cavitating flows via combined physical and numerical analysis, with particular emphasis on the interaction between unsteady cavitation development and structural deformations. Issues including unsteady cavitating flow structures and cavitation–structure interaction mechanism are discussed.

Structural subgrid-scale modeling for large-eddy simulation: A review

Accurately modeling nonlinear interactions in turbulence is one of the key challenges for large-eddy simulation（LES） of turbulence. In this article, we review recent studies on structural subgrid scale modeling, focusing on evaluating how well these models predict the effects of small scales. The article discusses a priori and a posteriori test results. Other nonlinear models are briefly discussed, and future prospects are noted.

Unsteady bio-fluid dynamics in flying and swimming

Flying and swimming in nature present sophisticated and exciting ventures in biomimetics, which seeks sustainable solutions and solves practical problems by emulating nature＇s time-tested patterns, functions, and strategies. Bio-fluids in insect and bird flight, as well as in fish swimming are highly dynamic and unsteady; however, they have been studied mostly with a focus on the phenomena associated with a body or wings moving in a steady flow. Characterized by unsteady wing flapping and body undulation, fluid-structure interactions, flexible wings and bodies, turbulent environments, and complex maneuver, bio-fluid dynamics normally have challenges associated with low Reynolds number regime and high unsteadiness in modeling and analysis of flow physics. In this article, we review and highlight recent advances in unsteady bio-fluid dynamics in terms of leading-edge vortices, passive mechanisms in flexible wings and hinges, flapping flight in unsteady environments, and micro-structured aerodynamics in flapping flight, as well as undulatory swimming, flapping-fin hydrodynamics, body–fin interac-tion, C-start and maneuvering, swimming in turbulence,collective swimming, and micro-structured hydrodynamics in swimming. We further give a perspective outlook on future challenges and tasks of several key issues of the field.

Experimental investigation on the wake interference among wind turbines sited in atmospheric boundary layer winds

We examined experimentally the effects of incoming surface wind on the turbine wake and the wake interference among upstream and downstream wind turbines sited in atmospheric boundary layer（ABL） winds. The experiment was conducted in a large-scale ABL wind tunnel with scaled wind turbine models mounted in different incoming surface winds simulating the ABL winds over typical offshore/onshore wind farms. Power outputs and dynamic loadings acting on the turbine models and the wake flow characteristics behind the turbine models were quantified. The results revealed that the incoming surface winds significantly affect the turbine wake characteristics and wake interference between the upstream and downstream turbines. The velocity deficits in the turbine wakes recover faster in the incoming surface winds with relatively high turbulence levels. Variations of the power outputs and dynamic wind loadings acting on the downstream turbines sited in the wakes of upstream turbines are correlated well with the turbine wakes characteristics. At the same downstream locations, the downstream turbines have higher power outputs and experience greater static and fatigue loadings in the inflow with relatively high turbulence level, suggesting a smaller effect of wake interference for the turbines sited in onshore wind farms.

Delayed detached eddy simulations of fighter aircraft at high angle of attack

The massively separated flows over a realistic aircraft configuration at 40？, 50？, and 60？angles of attack are studied using the delayed detached eddy simulation（DDES）.The calculations are carried out at experimental conditions corresponding to a mean aerodynamic chord-based Reynolds number of 8.93 × 10~5 and Mach number of 0.088. The influence of the grid size is investigated using two grids, 20.0×10~6cells and 31.0 × 10~6 cells. At the selected conditions, the lift,drag, and pitching moment from DDES predictions agree with the experimental data better than that from the Reynoldsaveraged Navier–Stokes. The effect of angle of attack on the flow structure over the general aircraft is also studied, and it is found that the dominated frequency associated with the vortex shedding process decreases with increasing angle of attack.

Investigation of the Flow Structure on a Flat Plate Induced by Unsteady Plasma Actuation with DNS Methods

An investigation into the flow characteristic on a flat plate induced by an unsteady plasma was conducted with the methods of direct numerical simulations（DNS）.A simplified model of dielectric barrier discharge（DBD） plasma was applied and its parameters were calibrated with the experimental results.In the simulations,effects of the actuation frequency on the flow were examined.The instantaneous flow parameters were also drawn to serve as a detailed study on the behavior when the plasma actuator was applied to the flow.The result shows that induced by the unsteady actuation,a series of vortex pairs which showed dipole formation and periodicity distribution were formed in the boundary layer.The production of these vortex pairs indicated a strong energy exchange between the main flow and the boundary layer.They moved downstream under the action of the free stream and decayed under the influence of the fluid viscosity.The distance of the neighboring vortices was found to be determined by the actuation frequency.Interaction of the neighboring vortices would be ignored when the actuation frequency was too small to make a difference.

AN INVESTIGATION ON ASM TURBULENCE MODEL IN PREDICTIING NONISOTROPIC TURBULENT FLOWS

Non-isotropy is evident for flows with large streamline curvature. Algebraic stress model can fairly simulate non-isotropy of turbulence. Some flow situations were investigated with a well developed ASM model. The calcu- lated results show that the distributions of Reynolds stresses in different directions are obviously different. The non-isotropy of the flow is evident. The effect of streamline curvature on the flow feature of submerged water jet restricted by solid boundaries were studied, and its application in the design of energy dissipators of free trajectory jet type was disscused.

Two-phase turbulence models for simulating dense gas-particle flows

The two-fluid model is widely adopted in simulations of dense gas-particle flows in engineering facili- ties. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas-particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas-particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-O and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle-particle collision using the transport equation model for the two-phase velocity correlation.

Double-averaging analysis of turbulent kinetic energy fluxes and budget based on large-eddy simulation

The turbulent flow over a channel bed roughened by three layers of closely packed spheres with a Reynolds number of Re= 15 000 is investigated using the large eddy simulation（LES） and the double-averaging（DA） method. The DA velocity is compared with the results of the corresponding laboratory experiments to validate the LES results. The existence of the types of vortex structures is demonstrated by the Q-criterion above the permeable bed. The turbulent kinetic energy（TKE） fluxes and budget are quantified and discussed. The results show that the TKE fluxes are directed downward and downstream near the virtual bed level. In the TKE budget, the form-induced diffusion rate is significant in the vicinity of the crest bed level, and the TKE production rate and the dissipation rate attain their peaks at the crest bed level and decrease sharply below it.

Vector properties of a tunable random electromagnetic beam in non-Kolmogrov turbulence

Analytical formulas for a class of tunable random electromagnetic beams propagating in a turbulent atmosphere through a complex optical system are derived with the help of a tensor method. One finds that the far field intensity distribution is tunable by modulating the source correlation structure function. The on-axis spectral degree of polarization monotonically increases to the same value for different values of order M in free space while it returns to the initial value after propagating a sufficient distance in turbulence. Furthermore, it is revealed that the state of polarization is closely determined by the initial correlation structure rather than by the turbulence parameters.

Large-eddy simulation: Past, present and the future

Large-eddy simulation（LES） was originally proposed for simulating atmospheric flows in the 1960 s and has become one of the most promising and successful methodology for simulating turbulent flows with the improvement of computing power. It is now feasible to simulate complex engineering flows using LES. However, apart from the computing power, significant challenges still remain for LES to reach a level of maturity that brings this approach to the mainstream of engineering and industrial computations. This paper will describe briefly LES formalism first, present a quick glance at its history, review its current state focusing mainly on its applications in transitional flows and gas turbine combustor flows, discuss some major modelling and numerical challenges/issues that we are facing now and in the near future, and finish with the concluding remarks.

Prediction of oil-water flow patterns, radial distribution of volume fraction, pressure and velocity during separated flows in horizontal pipe

Flow of two immiscible fluids gives rise to variety of flow patterns,which influence transportation process.In this work,we present detailed analysis on the prediction of flow pattern maps and radial distribution of volume fraction,pressure and velocity of a pair of immiscible liquids through a horizontal pipeline by computational fluid dynamics（CFD） simulation using ANSYS FLUENT 6.3.Moderately viscous oil and water have been taken as the fluid pair for study.Volume of fluid（VOF） method has been employed to predict various flow patterns by assuming unsteady flow,immiscible liquid pair,constant liquid properties,and co-axial flow.From the grid independent study,we have selected 47 037 number of quadrilateral mesh elements for the entire geometry.Simulation successfully predicts almost all the flow patterns（viz.,plug,slug,stratified wavy,stratified mixed and annular）,except dispersion of oil in water and dispersion of water in oil.The simulated results are validated with experimental results of oil volume fraction and flow pattern map.Radial distribution of volume fraction,pressure and velocity profiles describe the nature of the stratified wavy,stratified mixed and annular flow pattern.These profiles help to developing the phenomenological correlations of interfacial characteristics in two-phase flow.

LARGE-EDDY AND DETACHED-EDDY SIMULATIONS OF THE SEPARATED FLOW AROUND A CIRCULAR CYLINDER

The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes （URANS）. The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale （SGS） model in LES, the DES-based Spalart-Allmaras （S-A） and κ-ω Shear-Stress- Transport （SST） models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.