Three-Dimensional Direct Numerical Simulation of Flow past A Near-Wall Circular Cylinder:Combined Effects of Gap Ratio and Boundary Layer Thickness on Flow Profiles and Pressure Distribution

Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Numerical Modeling of Mass Transfer in the Interaction between River Biofilm and a Turbulent Boundary Layer

In this article dedicated to the modeling of vertical mass transfers between the biofilm and the bulk flow, we have, in the first instance, presented the methodology used, followed by the presentation of various results obtained through analyses conducted on velocity fields, different fluxes, and overall transfer coefficients. Due to numerical constraints (resolution of relevant spatial scales), we have restricted the analysis to low Schmidt numbers (Sc=0.1, Sc=1, and Sc=10) and a single roughness Reynolds number (Re*=150). The analysis of instantaneous concentration fields from various simulations revealed logarithmic concentration profiles above the canopy. In this zone, the concentration is relatively homogeneous for longer times. The analysis of results also showed that the contribution of molecular diffusion to the total flux depends on the Schmidt number. This contribution is negligible for Schmidt numbers Sc≥0.1, but nearly balances the turbulent flux for Sc=0.1. In the canopy, the local Sherwood number, given by the ratio of the total flux (within or above the canopy) to the molecular diffusion flux at the wall, also depends on the Schmidt number and varies significantly between the canopy and the region above. The exchange velocity, a purely hydrodynamic parameter, is independent of the Schmidt number and is on the order of 10% of in the present case. This study also reveals that nutrient absorption by organisms near the wall depends on the Schmidt number. Such absorption is facilitated by lower Schmidt numbers.

DIRECT NUMERICAL SIMULATION OF VIBRATIONINDUCED T-S WAVE ON BOUNDARY LAYER WALL

The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.

Direct Numerical Simulation of Particle Dispersion in Gas-Solid Compressible Turbulent Jets

A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.

Direct numerical simulation of viscoelastic-fluid-based nanofluid turbulent channel flow with heat transfer

Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.

AN ANALYSIS OF SUBGRID-RESOLVED SCALE INTERACTIONS WITH USE OF RESULTS FROM DIRECT NUMERICAL SIMULATIONS

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Direct numerical simulation study of the interaction between the polymer effect and velocity gradient tensor in decaying homogeneous isotropic turbulence

Direct numerical simulation of decaying homogeneous isotropic turbulence （DHIT） of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase ＂strain generation resistance＂ and ＂vorticity generation resistance＂, i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.

Validation of the RANS-SOM Combustion Model Using Direct Numerical Simulation of Incompressible Turbulent Reacting Flows

The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

Numerical simulation on the NACA0018 airfoil self-noise generation

Airfoil self-noise is a common phenomenon for many engineering applications. Aiming to study the underlying mechanism of airfoil self-noise at low Mach number and moderate Reynolds number flow, a numerical investigation is presented on noise generation by flow past NACA0018 airfoil. Based on a high-order accurate numerical method, both the near-field hydrodynamics and the far-field acoustics are computed simultaneously by performing direct numerical simulation. The mean flow properties agree well with the experimental measurements. The characteristics of aerodynamic noise are investigated at various angles of attack. The obtained results show that inclining the airfoil could enlarge turbulent intensity and produce larger scale of vortices. The sound radiation is mainly towards the upper and lower directions of the airfoil surface. At higher angle of attack, the tonal noise tends to disappear and the noise spectrum displays broad-band features.

Direct numerical simulation of matrix diffusion across fracture/matrix interface

Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relatively large diffusion coefficients and small fracture spacings, due to an intrinsic assumption of an equilibrium concentration profile in the matrix blocks. Motivated and required by practical applications, we propose a direct numerical simulation （DNS） approach without any empirical assumptions. A three-step Lagrangian algorithm was developed and validated to directly track the particle dynamics across the fracture/matrix interface, where particle＇s diffusive displacement across the discontinuity is controlled by an analytical, one-side reflection probability. Numerical experiments show that the DNS approach is especially efficient for small diffusion coefficients and large fracture spacings, alleviating limitations of the classical modeling approach.

Direct numerical simulation of elastic turbulence and its mixing-enhancement effect in a straight channel flow

In this paper,we present a direct numerical simulation（DNS） of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number（Re = 1） in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.

Numerical Simulation of Wake Deflection Control around NACA0012 Airfoil Using Active Morphing Flaps

This study demonstrates an active flow control for deflecting a direction of wake vortex structures behind a NACA0012 airfoil using an active morphing flap. Two-dimensional direct numerical simulations are performed for flows at the chord Reynolds number of 10,000, and the vortex pattern in the controlled and noncontrolled wakes as well as the effect of an actuation frequency on the control ability are rigorously investigated. It is found that there is an optimum actuation-frequency regime at around F + = 2.00 which is normalized by the chord length and freestream velocity. The wake vortex pattern of the well-controlled case is classified as the 2P wake pattern according to the Williamson’s categorization [1] [2], where the forced oscillation frequency corresponds to the natural vortex shedding frequency without control. The present classification of wake vortex patterns and finding of the optimum frequency regime in the wake deflection control can lead to a more robust design suitable for vortex-induced-vibration (VIV) related engineering systems.

Characteristics and generation of elastic turbulence in a three-dimensional parallel plate channel using direct numerical simulation

Direct numerical simulations（DNSs） of purely elastic turbulence in rectilinear shear flows in a three-dimensional（3D） parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures＇ behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures＇ elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.

A particle-resolved heat-particle-fluid coupling model by DEM-IMB-LBM

Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.

Boundary-layer disturbances subjected to free-stream turbulence and simulation on bypass transition

The phenomena associated with the entrainment of free-stream turbulence （FST） into boundary-layer flows are relevant for a number of subjects. It has been be- lieved that the continuous spectra of the Orr-Sommerfeld （O-S）/Squire equations describe the entrainment process, and thus they are used to specify the inlet condition in simulation of bypass transition. However, Dong and Wu （Dong, M. and Wu, X. On continuous spectra of the Orr-Sommerfeld/Squire equations and entrainment of free-stream vortical disturbances. Journal of Fluid Mechanics, 732, 616-659 （2013）） pointed out that continuous spectra exhibit several non-physical features due to neglecting the non-parallelism. They further proposed a large-Reynolds-number asymptotic approach, and showed that the non-parallelism is a leading-order effect even for the short-wavelength disturbance, for which the response concentrates in the edge layer. In this paper, the asymptotic solution is verified numerically by studying its evolution in incompressible boundary layers. It is found that the numerical results can be accurately predicted by the asymptotic solution, implying that the latter is adequate for moderate Reynolds numbers. By introducing a series of such solutions as the inflow perturbations, the bypass transition is investigated via the direct numerical simulation （DNS）. The transition processes, including the evolution of streaks, the amplification of secondary-instability modes, and the emergence of turbulent spots, agree with the experimental observations.

NUMERICAL INVESTIGATION OF EVOLUTION OF DISTURBANCES IN SUPERSONIC SHARP CONE BOUNDARY LAYERS

The spatial evolution of 2-D disturbances in supersonic sharp cone boundary layers was investigated by direct numerical simulation （DNS） in high order compact difference scheme. The results suggested that, although the normal velocity in the sharp cone boundary layer was not small, the evolution of amplitude and phase for small amplitude disturbances would be well in accordance with the results obtained by the linear stability theory （LST） which supposes the flow was parallel. The evolution of some finite amplitude disturbances was also investigated, and the characteristic of the evolution was shown. Shocklets were also found when the amplitude of disturbances increased over some value.

Numerical study on evolution of subharmonic varicose low-speed streaks in turbulent channel flow

The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose （SV） mode, is studied by a direct numerical simulation （DNS） method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous （SS） low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing ＋~. When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.