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.

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.

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.

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.

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.

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.

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.

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.

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.

Lagrangian time scales and its relationship to Eulerian equivalents in turbulent channel flow

Lagrangian and Eulerian time scales were obtained from the direct numerical simulation of turbulent channel flow at two Reynolds numbers based on the friction velocity and channel half-height, Rer= 80, 100. The Lagrangian integral time scales and time microscales were compared to their Eulerian equivalents. It is found that the ratio of Lagrangian to TL Eulerian integral time scales is given by TE/TiE= 1 ＋ 0.1y＋ for y＋ ≤ 10, and that the ratios between the Lagrangian to theEulerian time microscales are almost the same irrespective of the components. Those increase with y＋ are approximated by ≈ 2.75 - 1.75 exp （-v＋/a） . These results also show that these expressions are independent of the Reynolds number.

A Novel Non-Equidifferent Optical APPM Mapping Scheme for Strong Turbulent Atmospheric Channel

An optical Amplitude and Pulse Position Modulation(APPM) mapping scheme for strong turbulent atmospheric channel is proposed to optimize Bit Error Rate(BER) performance.In this scheme,a nonequidifferent amplitude series is designed based on quantitative BER analysis of the specific A×M APPM demapping procedures containing time slot selection and amplitude decision in selected time slot,which are different from traditional ones.Simulation results of 4×4,4×8 and 4×16 APPM show 4,3.4 and 6.9 d B SNR gain against traditional APPM scheme respectively.Thus significant BER performance improvement is achieved which helps to enhance reliability of freespace optical communication systems.

Formation of Gel-Like Shear-Induced Structure with Dosing of Dilute Surfactant Solution and Its Effect on Turbulent Channel Flow

A shear-induced structure (SIS) is formed under appropriate concentration and shear conditions in a surfactant micellar solution. In this study, we performed experiments of surfactant solution dosing in a fully developed two-dimensional turbulent channel flow from a sintered metallic wire mesh plate attached to a side wall. We investigated the behavior of the solution under the elongation during its passing through the wire mesh and under the strong shear due to the channel flow. It was confirmed that the dosed solution containing a laser dye was visualized by a laser sheet, and the accumulated gel from the wire mesh formed a layer and developed with time. Consequently, on dosing the dilute surfactant solution from the wire mesh, a gel-like SIS layer was formed, which majorly covered the wire mesh plate. The gel-like SIS layer on the wire mesh plate acted as a sticky solid and restricted the flow in the channel. This layer continued to grow while dosing, owing to which the pressure drop of the channel flow significantly increased. The gel-like SIS layer grew rapidly even in the turbulent flow and reached the equilibrium thickness. After the termination of the dosing, the gel layer collapsed gradually. In addition, the thickness of the gel-like SIS layer (indicating the strength indirectly) strongly depended on the surfactant concentration and the elongation rate in the wire mesh.

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.

Hybrid URANS/LES Method of Flow Fields in Axial-flow Compressor Rotor Rotor

Accurate and efficient prediction of the aerodynamic performance and flow details of axial-flow com-pressors is of great engineering application value for the aerodynamic design and flow control of axial-flow compres-sors.In this work,a delayed detached eddy simulation method is developed and applied to numerically simulate the tur-bulent channel flow and the aerodynamic performance of NASA Rotor 35.Several acceleration techniques including parallel implementation are also used to speed up the iteration convergence.The mean velocity distribution and Reyn-olds stress distribution in the boundary layer of turbulent channel flow and the aerodynamic performance curve of NASA Rotor 35 are predicted.The good agreement between the present delayed detached eddy simulation results and the available direct numerical simulation results or experimental data confirms the effectiveness of the developed meth-od in the accurate and efficient prediction of complex flow in turbomachinery.

Investigation on the vortical structures by the Liutex method in turbulent channels at Re_(τ)=180 with scalloped and triangular riblet control

Riblets are a series of small protrusions formed along the flow direction,which have been extensively studied as a passive turbulent drag reduction technique.Experiments and numerical simulations have shown that well-designed riblets can significantly reduce drag in turbulent flows,making them highly promising and valuable for various applications.In this study,we focus on a scalloped riblet,which is designed by smoothly connecting two third-order polynomials,and thus the sharpness of the tip and the curvature of the valley can be well defined.We conduct direct numerical simulations of turbulent channel with smooth plate,scalloped riblet-mounted and triangular riblet-mounted walls.Width in wall units of W^(+)=20 and height-width ratio ofγ=0.5 are selected for both riblet cases.Compared with the smooth plate case,the scalloped riblet case achieves an 8.68%drag reduction,while the triangular riblet case achieves a 4.79%drag reduction.The obtained drag reduction rate of the triangular riblet is consistent with previous experiments and simulations,and the results indicate that the scalloped riblet is more effective in reducing drag and deserves further investigation.We compare turbulent statistics of the scalloped riblet case with those of the triangular riblet case.The mean velocity profiles of riblets are similar,but both the Reynolds shear stress and second-order statistics of velocity fluctuations and Liutex are significantly reduced in the scalloped riblets controlled turbulent channel,indicating that the scalloped riblet can more effectively suppress the spanwise and wall-normal turbulent intensity near the wall.We also compare the pre-multiplied spectra of streamwise velocity and streamwise Liutex component for the three cases to investigate the energy distribution and characteristics of Liutex distribution.The Liutex vortex identification method is also utilized to analyze the instantaneous flow field,which provides insights into the flow field and could be beneficial for the further optimization of riblet.