Active control of turbulence for drag reduction based on the detection of near-wall streamwise vortices by wall information

The spatial relations between the measurable wall quantities （streamwise shear stress τwx, spanwise shear stress τwz, and pressure fluctuations Pw） and the near-wall streamwise vortices （NWSV） are investigated via direct numerical simulation （DNS） databases of fully developed turbulent channel flow at a low Reynolds number. In the stan- dard turbulent channel flow, the results show that all the wall measurable variables are closely associated with the NWSV. But after applying a stochastic interference, the relation based on τwx breaks down while the correlations based on Pw and τwz are still robust. Hence, two wall flow quantities based on Pw and τwz are proposed to detect the NWSV. As an appli- cation, two new control schemes are developed to suppress the near-wall vortical structures using the actuation of wall blowing/suction and obtain 16 % and 11% drag reduction, respectively.

Large eddy simulation of high-Reynolds-number atmospheric boundary layer flow with improved near-wall correction

It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.

PROBABILITY DISTRIBUTION FUNCTION OF NEAR-WALL TURBULENT VELOCITY FLUCTUATIONS

By large eddy simulation （LES）, turbulent databases of channel flows at different Reynolds numbers were established. Then, the probability distribution functions of the streamwise and wall-normal velocity fluctuations were obtained and compared with the corresponding normal distributions. By hypothesis test, the deviation from the normal distribution was analyzed quantitatively. The skewness and flatness factors were also calculated. And the variations of these two factors in the viscous sublayer, buffer layer and log-law layer were discussed. Still illustrated were the relations between the probability distribution functions and the burst events-sweep of high-speed fluids and ejection of low-speed fluidsIin the viscous sub-layer, buffer layer and loglaw layer. Finally the variations of the probability distribution functions with Reynolds number were examined.

Effect of Magnetic Mirror on the Asymmetry of the Radial Profile of Near-Wall Conductivity in Hall Thrusters

Considering the actual magnetic field configuration in a Hall thruster, the effect of magnetic mirror on the radial profile of near-wall conductivity （NWC） is studied in this paper. The plasma electron dynamic process is described by the test particle method. The Monte Carlo scheme is used to solve this model. The radial profile of electron mobility is obtained and the role of magnetic mirror in NWC is analysed both theoretically and numerically. The numerical results show that the electron mobility peak due to NWC is inversely proportional to the magnetic mirror ratio and the asymmetry of electron mobility along the radial direction gets greater when the magnetic mirror is considered. This effect indicates that apart from the disparity in the magnetic field strength, the difference in the magnetic mirror ratio near the inner and outer walls would actually augment the asymmetry of the radial profile of NWC in Hall thrusters.

The Effect of the Gap Ratio on the Flow and Heat Transfer over a Bluff Body in Near-Wall Conditions

In order to study the effect of different gap ratios on the thermofluid-dynamic field around a bluff body located in proximity to a heated wall,a series of experiments and numerical simulations have been conducted.The former were carried out using an open circulating water tank experimental platform and a single cylinder and square column as geometrical models(their characteristic length being D).The latter were based on the well-known SIMPLE algorithm for incompressible flow.The results show that the gap ratio is an important factor affecting the wake characteristics of near-wall bluff bodies.When the gap ratio is small,the influence of the wall on the bluff body wake is large.With an increase in the gap extension,periodic vortex shedding is enabled and heat transfer is strengthened accordingly;in addition,the vortex shedding period is larger for the square column.The square column displays hysteresis compared with the cylinder at the same gap ratio(the critical gap ratio of cylinder is 0.2~0.4,while that of square column is 0.4
0.6).

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.

MATHEMATICAL MODELING OF NEAR-WALL FLOWS OF TWO- PHASE MIXTURE WITH EVAPORATING DROPLETS

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream.

SECOND-MOMENT CLOSURE FOR MODELLING THE NEAR-WALL TURBULENCE IN 3D MEAN FLOWS

A second-moment closure for the near-wall turbulence is proposed. The limiting behaviour of this closure near a wall is consistent with that of the exact Reynolds-stress transport equations, and it converts asymptotically into a high- Reynolds-number closure remote from the wall. The closure is applied to a pressure- driven 3D transient channel flow. The predicted results are in fair agreement with the DNS data.

SPH Modelling of the Vortex-Induced Vibration of A Near-Wall Cylinder

The frequency-locked phenomenon commonly occurs in the vortex-induced vibration(VIV)of bluff bodies.Numerical simulation of this lock-in behavior is challenging,especially when the structure is positioned in close proximity to a solid boundary.To establish a robust simulator,an enhanced smoothed particle hydrodynamic(SPH)model is developed.The SPH model incorporates a particle shifting algorithm and a pressure correction algorithm to prevent cavity formation in the structure's wake area.A damping zone is also established near the outlet boundary to dissipate the vortices that shed from the structure.Additionally,GPU parallel technology is implemented to enhance the SPH model's computational efficiency.To validate the mo del,the predicted results are compared with the available refere nce data for flow past both stationary and oscillating cylinders.The verified SPH model is then employed to comparatively investigate the motion re sponse,lift characteristic,and vortex shedding mode of cylinders with and without accounting for the effect of boundary layers.Numerical analyses demonstrate that the developed SPH model is a proficient tool for efficiently simulating the vibration of near-wall bluff bodies at low Reynolds number.

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.

A cellular scale numerical study of the effect of mechanical properties of erythrocytes on the near-wall motion of platelets

The effect of mechanical properties of erythrocytes on the near-wall motion of platelets was numerically studied with the immersed boundary method. Cells were modeled as viscous-fluid-filled capsules surrounded by hyper-elastic membranes with negligible thickness. The numerical results show that with the increase of hematocrit, the near-wall approaching of platelets is enhanced, with which platelets exhibit larger deformation and orientation angle of its near-wall tank-treading motion, and the lateral force pushing platelets to the wall is increased with larger fluctuation amplitude. Meanwhile the near-wall approaching is reduced by increasing the stiffness of erythrocytes.

Jetting of a near-wall cavitation bubble induced by another tandem bubble

Double bubbles near a rigid wall surface collapse to produce a significant jet impact,with potential applications in surface cleaning and ultrasonic lithotripsy.However,the dynamic behaviors of near-wall bubbles remain unexplored.In this study,we investigate the jetting of a near-wall bubble induced by another tandem bubble.We define two dimensionless standoff distances,γ_(1),γ_(2),to represent the distances from the center of the near-wall bubble to the rigid wall and the center of controlling bubble to the center of the near-wall bubble,respectively.Our observations reveal three distinct jetting regimes for the near-wall bubble:transferred jetting,double jetting,and directed jetting.To further investigate the jetting mechanism,numerical simulations are conducted using the compressibleInterFoam solver in the open-source framework of OpenFOAM.A detailed analysis shows that the transferred jet flow is caused by the pinch-off resulting from the axial contraction velocity at the lower end of the near-wall bubble being greater than the vertical contraction velocity,leading to a maximum jet velocity of 682.58 m/s.In the case of double jetting,intense stretching between the controlling bubble and the wall leads to a pinch-off and a double jetting with a maximum velocity of 1096.29 m/s.The directed jet flow is caused by the downward movement of the high-pressure region generated by the premature collapse of the controlling bubble,with the maximum jet velocity reaching 444.62 m/s.

A theoretical model for Reynolds-stress and dissipation-rate budgets in near-wall region

A 3-D wave model for the turbulent coherent structures in near-wall region is proposed. The transport nature of the Reynolds stresses and dissipation rate of the turbulence kinetic energy are shown via computation based on the theoretical model. The mean velocity profile is also computed by using the same theoretical model. The theoretical results are in good agreement with those found from DNS, indicating that the theoretical model proposed can correctly describe the physical mechanism of turbulence in near wall region and it thus possibly opens a new way for turbulence modeling in this region.

Extension of Near-Wall Domain Decomposition to Modeling Flows with Laminar-Turbulent Transition

The near-wall domain decomposition method(NDD)has proved to be very efficient for modeling near-wall fully turbulent flows.In this paper the NDD is extended to non-equilibrium regimeswith laminar-turbulent transition(LTT)for the first time.The LTT is identified with the use of the e^(N)-method which is applied to both incompressible and compressible flows.TheNDD ismodified to take into account LTT in an efficientway.In addition,implementation of the intermittency expands the capabilities of NDD to model non-equilibrium turbulent flows with transition.Performance of the modified NDD approach is demonstrated on various test problems of subsonic and supersonic flows past a flat plate,a supersonic flow over a compression corner and a planar shock wave impinging on a turbulent boundary layer.The results of modeling with and without decomposition are compared in terms of wall friction and show good agreement with each other while NDD significantly reducing computational resources needed.It turns out that the NDD can reduce the computational time as much as three times while retaining practically the same accuracy of prediction.

气体引射对临近空间升力体飞行器气动力热影响

针对临近空间高速飞行器的降热减阻需求,采用三维数值模拟手段研究了临近空间典型升力体飞行器气动力/热特性,分析了壁面气体质量引射对气动力/热特性的影响规律和作用机理。结果表明:壁面质量引射会略微增大激波脱体距离,对激波层和速度边界层有增厚作用,降低了边界层内的速度梯度,使得表面剪切力减小,从而具有较好的减阻效果。对于典型的升力体构型,在高度为60~70 km时,壁面质量引射减阻效率均不低于47.5%;在高度为80 km时,壁面质量引射减阻效率不低于10.2%;且壁面质量引射减阻效果随攻角增大而增强,随飞行高度升高而减弱。同时,壁面引射排出的低温气体对飞行器具有隔热保护作用,因此壁面质量引射能显著降低迎风面的热流,具有较好的降热效果。

近壁旋转圆柱流场特性数值模拟分析

[目的]为探讨近壁旋转圆柱尾流及流体力特性,对典型间隙比下旋转圆柱绕流进行研究。[方法]对雷诺数Re=200下3种典型间隙比(G/D=0.2,0.8,1.4)的旋转圆柱绕流展开数值模拟,对比不同间隙比和转速比下的圆柱尾流及流体力特性。[结果]结果显示:当G/D=0.2时,圆柱表面脱涡会受到显著抑制,圆柱表面升阻力无波动;当G/D=0.8和1.4且转速比较低时,会发生“尾流涡”脱落现象,其结构与2S模式相似,升阻力系数呈正弦周期性波动,振幅较小;当正旋转速较大时,圆柱表面无漩涡脱落,形成稳定的D模式尾流(随转速比增大由D+模式变为D-模式),“尾流涡层”与“壁面涡层”发生分离,“壁面涡”呈现多周期性脱落现象,升阻力系数呈多周期波动,振幅显著增大;当反旋转速较大时,圆柱表面被一层正涡量的涡层包裹,漩涡脱落受到显著抑制,升阻力无波动。[结论]所得结论可为高效流动控制技术发展提供参考。

近壁面深水爆炸气泡射流演化特性的研究

近壁面气泡载荷是水中战斗部对障碍物毁伤作用的重要载荷之一。以近壁面深水爆炸气泡脉动及射流载荷特征问题为中心,开展了模拟300 m水深环境下的爆炸实验。通过高速摄影仪得到了气泡射流的演化过程;利用Autodyn轴对称模型对气泡射流的演化过程进行了计算,分析了比例距离对气泡脉动周期、射流演化形成时刻及射流强度特征的影响;探讨了壁面近场范围内的压力特征,总结得到射流载荷的演化规律。可为深水环境下水中战斗部对障碍物的毁伤效应研究提供理论参考。

近场爆炸下波纹双钢板混凝土组合墙板的损伤破坏及抗爆性能

相对于传统的钢筋混凝土墙板和平面双钢板-混凝土组合墙板(profiled double-skin composite wall,PDSCW),波纹双钢板-混凝土组合墙板(concrete-infilled double steel corrugatedplate wall,CDSCW)具有更高的轴向抗压承载力、更大的侧向抗弯刚度以及良好的抗冲击和抗震性能,在船舶和军事领域有广阔的应用前景。制作2种CDSCW试件,通过近场爆炸试验对比分析了2种试件的损伤模式及动态响应;采用ANSYS/LS-DYNA软件建立了CDSCW和PDSCW的有限元模型,研究了近场爆炸下2种混凝土组合墙板的损伤机理和爆炸响应,并与试验结果进行对比分析;分析了混凝土厚度、钢板厚度、药量对CDSCW抗爆性能的影响。结果表明:近场爆炸作用下,相较于PDSCW,相同混凝土方量和尺寸(长、宽)的CDSCW具有更大的抗弯刚度、更强的抗变形能力以及更优的抗爆性能;增加波纹深度能有效提高CDSCW的抗爆性能,可为抗爆构件设计和相关工程研究提供参考。

基于浸没边界法的壁面凸起对仿生鱼游动影响的研究

流体流经壁面凸起后在壁面凸起的后方产生周期性脱落的尾涡,尾涡内流场复杂。当尾涡继续流动流经游动中的仿生鱼时,仿生鱼的水动力系数也随之发生剧烈变化。基于多重直接力浸没边界法,结合狄拉克函数,通过求解不可压缩Navier-Stokes(N-S)无量纲方程,对低雷诺数情况下的仿生鱼游动过程进行无量纲化数值模拟研究。为了验证算法和自主开发程序的准确性,对圆柱绕流和NACA0012外形鱼体游动流场进行了验证算例计算,计算结果与文献结果吻合较好,验证了算法和程序的有效性。在此基础上,研究固壁半圆形凸起对仿生鱼游动水动力系数的影响,计算结果表明:当仿生鱼距壁面较近时(d=0.25L),仿生鱼受到远离壁面方向的侧向力和顺时针方向力矩;当仿生鱼距壁面较远时(d=1.25L),仿生鱼受到指向壁面的侧向力和逆时针方向力矩。总体而言,壁面凸起使得仿生鱼更容易在与壁面有一定距离的高度区间内游动,壁面凸起对仿生鱼游动起到了一定减阻作用。

波动鳍极地两栖机器人近壁区鳍面运动仿真研究

随着北极海冰消融,北极的战略地位日益凸显,而两栖机器人在极地极端环境中具有广阔的应用前景。本文提出一种波动鳍极地两栖机器人概念样机,并建立波动鳍结构及运动学模型。针对近冰底观测任务,将波动鳍分为水平摆动和垂直摆动2种运动模式,通过数值仿真方法探究了波动鳍在近壁区不同运动模式下受力、压强以及涡量场的变化规律。结果表明,水平摆动的波动鳍离壁面越近,升力越大,且存在临界幅值,幅值高于临界值时,推力随之减小;幅值低于临界值时,推力随之增加。垂直摆动的波动鳍离壁面越近,推力越小,负升力越大。波动鳍推力和升力的变化原因是鳍面与壁面之间空间的变化,改变了流体流动结构,从而改变了波动鳍周围的压力场和涡量场。