On the capability of the curvilinear immersed boundary method in predicting near-wall turbulence of turbulent channel flows

The immersed boundary method has been widely used for simulating flows over complex geometries.However,its accuracy in predicting the statistics of near-wall turbulence has not been fully tested.In this work,we evaluate the capability of the curvilinear immersed boundary(CURVIB)method in predicting near-wall velocity and pressure fluctuations in turbulent channel flows.Simulation results show that quantities including the time-averaged streamwise velocity,the rms(root-mean-square)of velocity fluctuations,the rms of vorticity fluctuations,the shear stresses,and the correlation coefficients of u'and v"computed from the CURVIB simulations are in good agreement with those from the body-fitted simulations.More importantly,it is found that the time-averaged pressure,the rms and wavenumber-frequency spectra of pressure fluctuations computed using the CURVIB method agree well with the body-fitted results.

Enhancing heat transfer at the micro-scale using elastic turbulence

Small concentrations of a high-molecular-weight polymer have been used to create so-called ＂elastic tur- bulence＂ in a micro-scale serpentine channel geometry. It is known that the interaction of large elastic stresses created by the shearing motion within the fluid flow with streamline curvature of the serpentine geometry leads initially to a purely-elastic instability and then the generation of elastic turbulence. We show that this elastic turbulence enhances the heat transfer at the micro-scale in this geometry by up to 300% under creeping flow conditions in comparison to that achieved by the equivalent Newtonian fluid flow.

The properties of dilute debris flow and hyper-concentrated flow in different flow regimes in open channels

Particle Image Velocimetry(PIV) technique was used to test the analogues of hyperconcentrated flow and dilute debris flow in an open flume. Flow fields, velocity profiles and turbulent parameters were obtained under different conditions. Results show that the flow regime depends on coarse grain concentration. Slurry with high fine grain concentration but lacking of coarse grains behaves as a laminar flow. Dilute debris flows containing coarse grains are generally turbulent flows. Streamlines are parallel and velocity values are large in laminar flows. However, in turbulent flows the velocity diminishes in line with the intense mixing of liquid and eddies occurring. The velocity profiles of laminar flow accord with the parabolic distribution law. When the flow is in a transitional regime, velocity profiles deviate slightly from the parabolic law. Turbulent flow has an approximately uniform distribution of velocity and turbulent kinetic energy. The ratio of turbulent kinetic energy to the kinetic energy of time-averaged flow is the internal cause determining the flow regime: laminar flow(k/K<0.1); transitional flow(0.1< k/K<1); and turbulent flow(k/K>1). Turbulent kinetic energy firstly increases with increasing coarse grain concentration and then decreases owing to the suppression of turbulence by the high concentration of coarse grains. This variation is also influenced by coarse grain size and channel slope. The results contribute to the modeling of debris flow and hyperconcentrated flow.

NUMERICAL INVESTIGATION OF THREE-DIMENSIONAL VISCOUS INCOMPRESSIBLE FLOWS IN DIVERGENT CURVED CHANNELS AND TURBULENT MODEL STUDY

In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.

Numerical Research on the Fiber Suspensions in a Turbulent T-shaped Branching Channel Flow

The concentration and orientation of fiber in a turbulent T-shaped branching channel flow are investi-gated numerically. The Reynolds averaged Navier-Stokes equations together with the Reynolds stress turbulent model are solved for the mean flow field and the turbulent kinetic energy. The fluctuating velocities of the fluid are assumed as a random variable with Gaussian distribution whose variance is related to the turbulent kinetic energy. The slender-body theory is used to simulate the fiber motion based on the known mean and fluctuating velocities of the fluid. The results show that at low Reynolds number, fiber concentration is high in the flow separation regions, and fiber orientation throughout the channel is widely distributed with a slight preference of aligning along the horizontal axis. With increasing of Re, the high concentration region disappears, and fiber orientation becomes ho-mogeneous without any preferred direction. At high Reynolds number, fiber concentration increases gradually along the flow direction. The differences in the distribution of concentration and orientation between different fiber aspect ratio are evident only at low Re. Both Re and fiber aspect ratio have small effect on the variance of orientation angle.

Transport of dissolved oxygen at the sediment-water interface in the spanwise oscillating flow

The distribution and concentration of dissolved oxygen(DO)play important roles in aerobic heterotroph activities and some slow chemical reactions,and can affect the water quality,biological communities,and ecosystem functions of rivers and lakes.In this work,the transport of high Schmidt number DO at the sediment-water interface of spanwise oscillating flow is investigated.The volume-averaged Navier-Stokes(VANS)equations and Monod equation are used to describe the flow in the sediment layer and the sediment oxygen demand of microorganisms.The phase-averaged velocities and concentrations of different amplitudes and periods are studied.The dependence of DO transfer on the amplitude and period is analyzed by means of phase-average statistical quantities.It is shown that the concentration in the sediment layer is positively correlated with the turbulence intensity,and the DO concentration and penetration depth in the sediment layer increases when the period and amplitude of the oscillating flow increase.Moreover,in the presence of oscillating flow,a specific scaling relationship exists between the Sherwood number/oxygen consumption of aerobic heterotrophs and the Reynolds number.

Mechanism of controlling turbulent channel flow with the effect of spanwise Lorentz force distribution

A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.

EXPERIMENTAL STUDY ON TURBULENT FEATURES IN THE NEGATIVE TRANSPORT REGION OF ASYMMETRIC PLANE CHANNEL FLOW

Turbulent features of streamwise and vertical components of velocity in the negative transport region of asymmetric plane channel flow have been studied experimentally in details. Experiments show that turbulent fluctuations in negative transport region are suppressed, and their probability distributions are far from Gaussian. Besides, the skewness factors attain their negative maxima at the position of the maximum mean velocity, whereas the flatness factors attain their positive maxima at the same position.

Multi-scale analysis of subgrid stress and energy dissipation in turbulent channel flow

In present study, the subgrid scale （SGS） stress and dissipation for multiscale formulation of large eddy simulation are analyzed using the data of turbulent channel flow at Ret = 180 obtained by direct numerical simulation. It is found that the small scale SGS stress is much smaller than the large scale SGS stress for all the stress components. The dominant contributor to large scale SGS stress is the cross stress between small scale and subgrid scale motions, while the cross stress between large scale and subgrid scale motions make major contributions to small scale SGS stress. The energy transfer from resolved large scales to subgrid scales is mainly caused by SGS Reynolds stress, while that between resolved small scales and subgrid scales are mainly due to the cross stress. The multiscale formulation of SGS models are evaluated a priori, and it is found that the small- small model is superior to other variants in terms of SGS dissipation.

DEVELOPMENT OF HYBRID RANS / LES CODE AND ITS VALIDATION

A computational code based on the hybrid RANS-LES approach is developed.The hybrid approach combines the delayed detached-eddy simulation ( DDES ) with an improved RANS-LES hybrid model aiming at wall modeling in LES ( WMLES ) .In the code , the convective flux is solved using the fourth-order skew-symmetric scheme so as to diminish the negative effect of numerical dissipation.The Spalart-Allmaras ( S-A ) model is applied as a subgrid scale ( SGS ) model.To validate the developed code , homogeneous isotropic turbulence and turbulent channel flow are simulated and the results are compared with experimental data and DNS results.The results of the isotropic turbulence show that the fourth-order skew-symmetric scheme is adequate enough and the model works well coupling with the convective scheme.The results of the turbulent channel flow agree well with the DNS data , the predicted velocity profiles at Reynolds number from 178to 2 700match well with the Reichardt′s law , and the organized vortical structures are well captured.

Circulation characteristics of horseshoe vortex in scour region around circular piers

This paper presents an experimental investigation of the circulation of the horseshoe vortex system within the equilibrium scour hole at a circular pier, with the data measured by an acoustic Doppler velocimeter （ADV）. Velocity vector plots and vorticity contours of the flow field on the upstream plane of symmetry （y = 0 cm） and on the planes ：e3 cm away from the plane of symmetry Cv = ~3 cm） are presented. The vorticity and circulation of the horseshoe vortices were determined using the forward difference technique and Stokes theorem, respectively. The results show that the magnitudes of circulations are similar on the planes y = 3 cm and y = -3 cm, which are less than those on the plane y = 0 cm. The circulation decreases with the increase of flow shallowness, and increases with the densimetric Froude number. It also increases with the pier Reynolds number at a constant densimetric Froude number, or at a constant flow shallowness. The relative vortex strength （dimensionless circulation） decreases with the increase of the pier Reynolds number. Some empirical equations are proposed based on the results. The predicted circulation values with these equations match the measured data, which indicates that these equations can be used to estimate the circulation in future studies.

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

We analyze the error of large-eddy simulation(LES)in wall pressure fluctuation of a turbulent channel flow.To separate different sources of the error,we conduct both direct numerical simulations(DNS)and LES,and apply an explicit filter on DNS data to obtain filtered DNS(FDNS)data.The error of LES is consequently decomposed into two parts:The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation.The second part is the difference between LES and FDNS induced by the error of LES in velocity field.By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation,it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.

Artificial neural network-based subgrid-scale models for LES of compressible turbulent channel flow

Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.

Three-Dimensional Bursting Phenomena in Meander Channel

Experiments were conducted in a U-shaped open-channel flume with the intention of investigating the bursting phenomena in the meander channel. The experimental results of the secondary flow fields and the Reynolds shear stress distributions show that the velocity and velocity fluctuation in the transverse direction are not negligible. Moreover, the bursting process is investigated using the three-dimensional quadrant analysis, which is more accurate than using the traditional two-dimensional quadrant analysis for the meandering channel. It is obtained from the experimental results that the internal group of events occurs more frequently than the external group, particularly the internal ejection and internal sweep events. In addition, the transition probabilities of the movements, which are defined as the changes of events from the current situation to the next situation in a time series, show that the stable organizations of events are the most possible movements, whereas the cross organizations of events have the least possible movements.

Generalized Block Markov Superposition Transmission over Free-Space Optical Links

In free-space optical(FSO) communications, the performance of the communication systems is severely degraded by atmospheric turbulence. Channel coding and diversity techniques are commonly used to combat channel fading induced by atmospheric turbulence. In this paper, we present the generalized block Markov superposition transmission(GBMST) of repetition codes to improve time diversity. In the GBMST scheme, information sub-blocks are transmitted in the block Markov superposition manner, with possibly different transmission memories. Based on analyzing an equivalent system, a lower bound on the bit-error-rate(BER) of the proposed scheme is presented. Simulation results demonstrate that, under a wide range of turbulence conditions, the proposed scheme improves diversity gain with only a slight reduction of transmission rate. In particular, with encoding memory sequence(0, 0, 8) and transmission rate 1/3, a diversity order of eleven is achieved under moderate turbulence conditions. Numerical results also show that, the GBMST systems with appropriate settings can approach the derived lower bound, implying that full diversity is achievable.

Three dimensional flow motions in the viscous sublayer

We employ novel digital Fresnel reflection holography to capture the 3D flows within the viscous sublayerof a smooth-wall turbulent channel flow at Reτ=400.The measurements reveal unsteady and diverse flow patterns in the sublayer including nearly uniform high and low speed flows and strong small-scale(onthe order of viscous wall units)spanwise meandering motions.The probability density functions(PDFs)ofwall shear stresses show a clear discrepancy in high stress range with those from direct numerical simu-lation(DNS),which is attributed to the unresolved streamwise and spanwise motions by DNS.Moreover,the PDF of Lagrangian particle accelerations yields a stretched exponential shape like that in homogenousisotropic turbulence,indicating strong intermittency in the sublayer.We find a significant fraction of highaccelerations is associated with the small-scale meandering motions.Our study helps explain the effectof sublayer-scale roughness on reducing drag and flow separation reported in the literature.

Independent component analysis of streamwise velocity fluctuations in turbulent channel flows

Independent component analysis(ICA)is used to study the multiscale localised modes of streamwise velocity fluctuations in turbulent channel flows.ICA aims to decompose signals into independent modes,which may induce spatially localised objects.The height and size are defined to quantify the spatial position and extension of these ICA modes,respectively.In contrast to spatially extended proper orthogonal decomposition(POD)modes,ICA modes are typically localised in space,and the energy of some modes is distributed across the near-wall region.The sizes of ICA modes are multiscale and are approximately proportional to their heights.ICA modes can also help to reconstruct the statistics of turbulence,particularly the third-order moment of velocity fluctuations,which is related to the strongest Reynolds shear-stressproducing events.The results reported in this paper indicate that the ICA method may connect statistical descriptions and structural descriptions of turbulence.

A Novel Trellis Coded Overlapping Amplitude and Pulse Position Modulation Scheme for Gamma-Gamma Channel Free-Space Optical Communication

A novel trellis coded-4×8 overlapping amplitude and pulse position modulation(TC-4×8AOPPM) scheme is proposed to enhance bit error rate(BER) performance of free-space optical communication(FSO) system. In addition, an uncoded AOPPM referential scheme is also designed. The schemes manage to decrease BER by designing gamma-gamma(GG) channel applicable decoding and demodulation methods. Simulation results of 8, 16 and 64-state TC-4×8AOPPM show 2.5-3.3 dB SNR gain against traditional TC-4×8AOPPM scheme respectively. Thus significant BER performance improvement is achieved and the reliability of the FSO system is also enhanced.

Transient response of enstrophy transport to opposition control in turbulent channel flow

The transient response of the turbulent enstrophy transport to opposition control in the turbulent channel flow is studied with the aid of direct numerical simulation. It is found that the streamwise enstrophy and the spanwise enstrophy are suppressed by the attenuation of the stretching terms at first, while the vertical enstrophy is reduced by inhibiting the tilt of the mean shear. In the initial period of the control, the streamwise enstrophy evolves much slower than the other two components. The vertical vorticity component exhibits a rapid monotonic decrease and also plays an important role in the attenuation of the other two components.

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation （DNS） of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodi- cally well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.