基于有限分析法的平面铅垂壁羽流的相似解

应用计及浮力对湍动能及其耗散率的影响的k-epsilon湍流模式,结合有限分析法对密度差引起的平面铅垂紊动壁羽流进行了数值分析．在均匀环境条件下控制该类流动的连续性方程、流动方向的动量方程、浓度扩散方程、湍动能及其耗散率方程存在相似解．考虑到浮力通量守恒条件,应用有限分析法给出了铅直壁羽流的速度、相对密度差、湍动能及耗散率的分布,进而给出了各物理量最大值沿主流方向变化的关系式．湍流Schmidt数为1.0时的计算结果与实验资料吻合较好,表明应用有限分析法分析铅垂平面羽流是有效的,即在分析壁羽流时浮力对湍动能及其耗散率的影响应该予以考虑．     

明渠温差异重流中湍浮力回流的数学模型

本文应用湍浮力回流k-s双方程模型及SIMPLE算法,对浅水明渠温差异重流的产生及其界面掺混规律进行了分析和数值计算。计算给出了速度、温度、湍动能、湍动能耗散率和涡粘性系数等的分布,所得到的速度、温度分布与实验测定值符合较好,表明这一预测方法是行之有效的.     

湍流速度和耗散率联合几率密度函数的输运方程

根据湍流动能耗散率满足对数正态分布的假定,给出耗散率的随机模型,并对速度的Langevin模型进行了修正,从而建立起完全封闭的湍流速度和耗散率联合几率密度函数的输运方程.为检验模型的可靠性,用Monte Carlo方法对输运方程进行了数值求解,给出衰变湍流和均匀剪切湍流的满意结果.     

风浪发展对海面阻力系数的影响

考虑海气界面附近大气边界层内不同风浪发展程度下的海气动量交换 .建立有效波加表面粗糙度的海气界面模型 ,提出了表面粗糙度的计算公式 .采用κ-ε模式以计算湍流流场 .计算所得海面阻力系数CD 随风速及波龄变化的规律与实测结果符合良好;所得流场结构包括风速廓线、湍流动能、湍能耗散及湍流粘性系数等都较合理     

壁湍流多尺度湍涡结构推广的自相似标度律

在水槽中测量了中等雷诺数下平板湍流边界层中的瞬时流向速度的时间序列,验证了Benzi提出的推广的自相似标度律,用子波变换将壁湍流脉动速度分解为多尺度湍涡结构的速度,研究了每一个尺度的湍涡速度结构函数的推广的自相似标度律。主要结论如下:湍流的统计性质是自相似的,这不仅适用于充分发展湍流,而且适用于中等雷诺数和低雷诺数湍流,而且具有相同的标度指数;推广的自相似标度律的适用的尺度范围远远大于惯性子区的范围,可以一直延伸至耗散区的尺度范围;推广的自相似标度律不仅适用于均匀各向同性湍流,也适用于剪切湍流如边界层湍流。     

近壁区Reynolds应力及湍能耗散率输运特性的理论模型 *

提出在湍流边界层的近壁区采用三维不稳定波来描述湍流相干结构 ,然后根据理论模型对Reynolds应力及湍能耗散率的输运特性进行系统的计算和定量的分析 ,并计算了平均速度分布 .所得理论结果与直接数值模拟结果符合较好 ,表明本文方法正确地体现了湍流边界层近壁区的物理实质 .这不仅有益于对湍流机制的了解 ,而且可能为湍流的近壁模型化开辟一条新的途径 .     

底部加热肥皂泡上准二维湍流的数值模拟北大核心CSCD

底部加热的肥皂泡是一种全新的二维热对流模型,在实验中已发现肥皂泡上的岛涡运动规律与飓风轨迹规律一致.然而,肥皂泡的曲面特征对其准二维流场的数值模拟以及数据分析造成了较多困难.针对肥皂泡球面几何特征,该文介绍了其直接数值模拟(DNS)方法,及其流场空间波数谱、湍流通量和结构函数的计算分析方法.开展了Ra=3×10^(7),3×10^(9),3×10^(11)的数值计算,并获得了相应的波数谱、通量和湍流结构函数.计算结果表明,肥皂泡上速度的小尺度脉动特征满足Bo59的理论标度律,通过湍动能与拟涡能通量特征,发现在该准二维湍流场中存在湍流能量双级串现象.且随着Rayleigh数的增加,大尺度结构湍能量减小,更小尺度湍流结构能量增加.     

非线性耗散方程的Bcklund变换及对称约化

应用齐次平衡法思想,求出非线性耗散项的Burgers-Fisher方程Bcklund变换、若干对称约化和相似解.     

新型降压槽结构设计及流动特性分析

针对三元复合驱分注过程中低视粘度损失和高压降不能同时实现的问题,设计出了一种新型降压槽,运用CFD软件对其环空流道的流动特性进行数值模拟,研究了新型降压槽内的压力、速度、湍动能、视粘度和应变率分布情况.同时将其与环形降压槽进行对比,结果表明,相同长度的降压槽,新型降压槽的压差更大,视粘度损失更小.     

三角形穿孔翅片对自然对流传热的强化作用

研究一种设置水平矩形翅片对自然对流传热的强化作用,翅片内含三角形穿孔,三角形的底边平行并朝向翅片顶端.比较了这种多孔翅片和同等实心翅片的热耗散率.考虑的参数包括翅片及其穿孔的几何尺寸和热性能.讨论了翅片穿孔后引起的传热强化和翅片重量的减轻.结果表明,在某些三角形穿孔和穿孔间距的值域内,穿孔翅片改进了同等情况实心翅片的热耗散.当翅片导热率及其厚度增大时,穿孔翅片的传热也强化.     

Approximating the larger eddies in fluid motion V: Kinetic energy balance of scale similarity models

We first review a classical scale-similarity model used to simulate the motion of large eddies in a turbulent flow. The kinetic energy balance of this model is very unclear in theory. Experiments with it often have reported that an additional Smagorinski type subgridscale term is needed. This term is not benign; it can alter significantly the predicted long term dynamics of the large eddies. However, we also show that the principal of scale-similarity (introduced in 1980 by Bardina, Ferziger and Reynolds) can also give rise to other scale similarity models which have the correct kinetic energy balance.     

HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER-MESHKOV INSTABILITY TO LATE TIMES:DYNAMICS,DEPENDENCE ON INITIAL CONDITIONS,AND COMPARISONS TO EXPERIMENTAL DATA

The dynamics of the reshocked multi-mode Richtmyer-Meshkov instability is investigated using 513×257 2three-dimensional ninth-order weighted essentially nonoscillatory shock-capturing simulations.A two-mode initial perturbation with superposed random noise is used to model the Mach 1.5 air/SF6 Vetter-Sturtevant shock tube experiment. The mass fraction and enstrophy isosurfaces,and density cross-sections are utilized to show the detailed flow structure before,during,and after reshock.It is shown that the mixing layer growth agrees well with the experimentally measured growth rate before and after reshock.The post-reshock growth rate is also in good agreement with the prediction of the Mikaelian model.A parametric study of the sensitivity of the layer growth to the choice of amplitudes of the short and long wavelength initial interfacial perturbation is also presented.Finally,the amplification effects of reshock are quantified using the evolution of the turbulent kinetic energy and turbulent enstrophy spectra,as well as the evolution of the baroclinic enstrophy production,buoyancy production,and shear production terms in the enstrophy and turbulent kinetic transport equations.     

Investigation of physiological pulsatile flow in a model arterial stenosis using large-eddy and direct numerical simulations

Physiological pulsatile flow in a 3D model of arterial stenosis is investigated by using large eddy simulation (LES) technique. The computational domain chosen is a simple channel with a biological type stenosis formed eccentrically on the top wall. The physiological pulsation is generated at the inlet using the first harmonic of the Fourier series of pressure pulse. In LES, the large scale flows are resolved fully while the unresolved subgrid scale (SGS) motions are modelled using a localized dynamic model. Due to the narrowing of artery the pulsatile flow becomes transition-to-turbulent in the downstream region of the stenosis, where a high level of turbulent fluctuations is achieved, and some detailed information about the nature of these fluctuations are revealed through the investigation of the turbulent energy spectra. Transition-to-turbulent of the pulsatile flow in the post stenosis is examined through the various numerical results such as velocity, streamlines, velocity vectors, vortices, wall pressure and shear stresses, turbulent kinetic energy, and pressure gradient. A comparison of the LES results with the coarse DNS are given for the Reynolds number of 2000 in terms of the mean pressure, wall shear stress as well as the turbulent characteristics. The results show that the shear stress at the upper wall is low just prior to the centre of the stenosis, while it is maximum in the throat of the stenosis. But, at the immediate post stenotic region, the wall shear stress takes the oscillating form which is quite harmful to the blood cells and vessels. In addition, the pressure drops at the throat of the stenosis where the re-circulated flow region is created due to the adverse pressure gradient. The maximum turbulent kinetic energy is located at the post stenosis with the presence of the inertial sub-range region of slope −5/3.     

Numerical modeling on inlet aperture effects on flow pattern in primary settling tanks

Inlets should be designed to dissipate the kinetic energy or velocity head of the mixed liquor and to prevent short-circuiting, mitigate the effects of density currents, and minimize blanket disturbances. Flow in primary settling tank is simulated by means of computational fluid dynamics. The fluid is assumed incompressible and non-buoyant. A two-dimensional computational and one phase fluid dynamics model was built to simulate the flow properties in the settling tank including the velocity profiles, the flow separation area and kinetic energy. In this study, the RNG turbulent model was solved with the Navier–Stokes equations. In order to evaluate hydraulic influences on the velocity profile, separation length and kinetic energy, three different of opening positions and two and three aperture in inlets were simulated. The flow model uses to apply a fixed-grid of cells that are all rectangular faces; the fluid moves through the grid and free surfaces are tracked with the volume-of-fluid (VOF) technique. Effects of numbers and locations of inlet apertures on the flow field are presented and the results show the positions of inlet apertures are affected on the flow pattern in the settling basin and increasing the numbers of slots can reduce kinetic energy in the inlet zone and produce uniform flow.     

Fast spectral solutions of the double-gyre problem in a turbulent flow regime

Several semi-analytical models are considered for a double-gyre problem in a turbulent flow regime for which a reference fully numerical eddy-resolving solution is obtained. The semi-analytical models correspond to solving the depth-averaged Navier–Stokes equations using the spectral Galerkin approach. The robustness of the linear and Smagorinsky eddy-viscosity models for turbulent diffusion approximation is investigated. To capture essential properties of the double-gyre configuration, such as the integral kinetic energy, the integral angular momentum, and the jet mean-flow distribution, an improved semi-analytical model is suggested that is inspired by the idea of scale decomposition between the jet and the surrounding flow.     

固液两相流中的K-ε双方程湍流模式*

建立了固液两相流中更一般的K-ε双方程湍流模式。模化了固相和液相的连续方程、动量方程及K方程和ε方程。该湍流模型考虑了固液两相间速度的滑移,颗粒间的作用及相间作用。使用本文所建立的湍流模型,数值预测了一管湍流中的沙水混合流动,其预测结果与实验结果比较一致。     

Study of spectral characteristics of a homogeneous turbulent flow

A spectral representation of kinetic energy for a vortex cascade of instability in a compressible inviscid shear flow is considered, and the Rayleigh-Taylor instability is studied. A comparative analysis is given to the spectral decompositions of kinetic energy for both problems. The classical Kolmogorov −5/3 power law is proved to hold for developed turbulent flows.     

Hybrid LES-RANS Modeling of Complex Turbulent Flows

The paper describes a state-of-the-art hybrid LES-URANS method for the simulation of complex internal and external turbulent flows. Relying on a unified LES-URANS approach with a soft interface the methodology is designed for wall-bounded non-equilibrium flows. The unsteady Reynolds-averaged Navier-Stokes (URANS) mode within the hybrid approach is taken into account by an explicit algebraic Reynolds stress model (EARSM), which guarantees an appropriate representation of the anisotropic near-wall turbulence. All non-closed terms in the transport equation of the turbulent kinetic energy are modeled by enhanced formulations using the EARSM (production and diffusion term) and the splitting of the dissipation rate into a homogeneous and an inhomogeneous contribution. The former is expressed analytically by a Taylor series expansion of the homogeneous lateral Taylor microscale in the vicinity of the wall guaranteeing the correct asymptotic behavior. The latter is incorporated into the diffusion term. The interface location between the large-eddy simulation (LES) mode and the URANS mode is determined automatically on-the-fly based on the modeled turbulent kinetic energy and the distance to the wall. For transitional (external) flows an additional dynamic transition criterion is applied which determines the laminar and the turbulent flow regimes and is combined with the existing interface criterion. An internal flow over a periodic arrangement of hills and an external flow past a SD7003 airfoil with a laminar separation bubble are taken into account for a detailed evaluation of the method. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.