壁面加热湍流耗散率标度指数测量的实验研究

对中等雷诺数下壁面常温和壁面加热的平板湍流边界层中速度和温度粗粒化的耗散率结构函数标度指数进行了实验测量.用热线风速仪测量了风洞中壁面常温和加热的平板湍流边界层中不同法向位置的流向速度分量和温度的时间序列信号,研究了由于湍流边界层近壁区域相干结构的存在而导致的非各向同性、非均匀性对湍流耗散率结构函数标度指数的影响,研究发现,中等雷诺数下壁面加热的边界条件和剪切湍流的平均速度梯度对速度和温度耗散率结构函数的标度指数没有影响,均匀各向同性湍流的耗散率结构函数标度指数的层次结构模型对壁面加热平板湍流边界层的速度和温度耗散率结构函数的标度指数也是适用的.     

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.     

An assessment of turbulence models applied to the simulation of a two-dimensional submerged jet

This paper is concerned with the investigation of the performance of different turbulence models in the numerical prediction of transient flow caused by a confined submerged jet. Several widely used models, i.e., the standard k–ε, RNG k–ε, low Reynolds number k–ε models and the differential Reynolds stress model, as included in CFD codes, were compared with each other for a two-dimensional, incompressible, turbulent jet flow and with reported experimental data. A flapping oscillation was predicted regardless of the model used. A chosen Strouhal (St) number definition brought the fundamental frequencies from both the experiments and computations into close proximity. However, different turbulence models have exhibited quite different behaviours in terms of the frequency and regularity of the oscillation and in terms of the scale and duration of the vortices generated. All versions of the k–ε model yielded regular oscillations, which agree with experimental observations. On the other hand, the Reynolds stress (RS) model produced a complex pattern but a slower dissipation of vortices. In addition, some aspects of gridding and inflow representation are also discussed.     

Investigations of shear free turbulent diffusion in a rotating frame

We reconsider the problem of shear free turbulent diffusion in a rotating frame, rotating about x₁. Shear free turbulence is generated at a vibrating grid in the x₂–x₃ plane and diffuses away from the grid in x₁ direction. An important property of this flow case is that there is no mean flow‐velocity. With the help of Lie‐group methods Reynolds‐stress transport models can be analyzed for this kind of flow in a rotating frame. From the analysis it can be found, that the turbulent diffusion only influences a finite domain. Implicating this solution in the model equations shows that even fully nonlinear Reynolds‐stress transport models (non‐linear in the Reynolds‐stresses for the pressure‐strain model) are insensitive to rotation for this type of flow. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Coupled maps and analytic function spaces

We consider ergodic properties of weakly coupled analytic and expanding circle maps. For weak enough coupling a natural ergodic measure exists and exhibits exponent decay of time correlations. The marginal densities of the natural measure are analytic. A spatial decay of correlations (e.g. polynomial) in these densities may arise from a similar spatial decay of the couplings. The space of couplings and observables is a Banach algebra of analytic functions of infinitely many variables. This algebra acts upon a Banach module of complex measures with analytic marginal densities. A Perron-Frobenius type operator acts on the Banach module and we apply a re-summation technique to derive uniform bounds for this operator. Explicit bounds are calculated for some examples.     

Direct numerical simulation of compressible isotropic turbulence

Direct numerical simulation (DNS) of decaying compressible isotropic turbulence at turbulence Mach numbers of M_t = 0.2-0.7 and Taylor Reynolds numbers of 72 and 153 is performed by using the 7th order upwind-biased difference and 8th order center difference schemes. Results show that proper upwind-biased difference schemes can release the limit of“start-up” problem to Mach numbers. Compressibility effects on the statistics of turbulent flow as well as the mechanics of shocklets in compressible turbulence are also studied, and the conclusion is drawn that high Mach number leads to more dissipation. Scaling laws in compressible turbulence are also analyzed. Evidence is obtained that scaling laws and extended self similarity (ESS) hold in the compressible turbulent flow in spite of the presence of shocklets, and compressibility has little effect on scaling exponents.     

Comparison of turbulence models in simulating swirling pipe flows

Swirling flow is a common phenomenon in engineering applications. A numerical study of the swirling flow inside a straight pipe was carried out in the present work with the aid of the commercial CFD code fluent. Two-dimensional simulations were performed, and two turbulence models were used, namely, the RNG k–ε model and the Reynolds stress model. Results at various swirl numbers were obtained and compared with available experimental data to determine if the numerical method is valid when modeling swirling flows. It has been shown that the RNG k–ε model is in better agreement with experimental velocity profiles for low swirl, while the Reynolds stress model becomes more appropriate as the swirl is increased. However, both turbulence models predict an unrealistic decay of the turbulence quantities for the flows considered here, indicating the inadequacy of such models in simulating developing pipe flows with swirl.     

Langevin equation of a fluid particle in wall-induced turbulence

We derive the Langevin equation describing the stochastic process of fluid particle motion in wall-induced turbulence (turbulent flow in pipes, channels, and boundary layers including the atmospheric surface layer). The analysis is based on the asymptotic behavior at a large Reynolds number. We use the Lagrangian Kolmogorov theory, recently derived asymptotic expressions for the spatial distribution of turbulent energy dissipation, and also newly derived reciprocity relations analogous to the Onsager relations supplemented with recent measurement results. The long-time limit of the derived Langevin equation yields the diffusion equation for admixture dispersion in wall-induced turbulence.     

On the decay exponent of isotropic turbulence

It has long been observed that after a short initial transient period of time the decay of the velocity .uctuations of high Reynolds number homogeneous isotropic turbulence follows closely the algebraic law ∼ t^–n. From experiments and DNS data it is noticed that the numerical value for the exponent is n ≈ 1.2. In the last decade values close to the latter one have been adopted to calibrate turbulence models such as the k‐ε‐model. It will be shown that if such an algebraic decay exists it corresponds to the invariant (similarity) solution of the equations of fluid dynamics under the group of scaling. Most important it is shown that if such solutions exist the decay exponents are fixed due to certain invariants. For the Navier‐Stokes equation the invariant is a constant Reynolds number and it follows n = 1. For the Euler equation the value is prescribed to owing to the conservation of energy with finite initial energy. It is interesting to note that the experimental and DNS data indicate towards a decay induced by the Euler equations rather than the one conformal to the Navier‐Stokes equations.     

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

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

Theoretical Considerations about Near-Wall Turbulence and Resulting Flow Control Schemes

Based on analytical considerations about how near-wall turbulence needs to be modified in order to reduce the momentum loss towards solid walls and to yield lower energy losses, a mechanism of turbulent drag reduction is proposed. This mechanism suggests that drag reducing flow control at high Reynolds numbers should be designed to minimize the turbulent dissipation rate. A previously published approach on how a reduction of the turbulent dissipation in the near-wall region can be achieved is analyzed further. The obtained results provide some new insight on the parameters that need to be considered when designing flow control schemes for skin friction drag reduction. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Aerodynamics of an isolated slot-burner from a tangentially-fired boiler

The aerodynamic development of fully turbulent isothermal jets issuing from rectangular slot-burners was modelled by obtaining a solution to the Reynolds averaged Navier–Stokes equations. A finite-volume method was used with the standard k–ε, RNG k–ε and Reynolds stress turbulence models. The slot-burners were based on physical models, which were designed to be representative of typical burner geometries found in tangentially-fired coal boilers. Two cases were investigated, in which jets from three vertically stacked rectangular nozzles discharged at 90° and then 60° to the wall containing the burner. The nozzle angle had little effect on jet centreline velocity decay, with the 60° nozzle showing a marginally higher rate of decay. The jets from the 60° nozzles were found to deviate slightly from their geometric axis slightly due to internal pressure redistribution in the flow at the nozzles. The simulations were validated against the physical models and were found to reproduce the flow field of the jets accurately with the Reynolds stress model producing the best results.     

Φ—伪压缩型映象的具误差的lshikawa和Mann迭代程序的收敛性问题

引入Φ－伪压缩型映象的概论,并研究了此类映象的具误差的Ｉｓｈｉｋａｗａ和Ｍａｎｎ迭代程序的收敛性问题。本文结果改进和发展了许多人的最新结果。     

各向异性湍浮力模型的研究

通过对隐式代数应力模型应力结构的分析,认为湍浮力流动的应力由应变和浮力作用引起,从而利用应力输运的线性近似和平衡近似提出了非线性各向异性湍浮力模型·　该模型避免了数值奇异,其可靠性通过计算和实验的对比得到了验证·     

Numerical solution of unsteady incompressible Navier-Stokes equations using high-order method of lines

A high-order method of lines is devised for solving the unsteady incompressible Navier-Stokes equations in the vorticity-stream function formulation. The vorticity transport equation is solved by the eight- or tenth-order method of lines and the Poisson equation for the stream function is solved by a high-order multigrid method. The numerical results of the two-dimensional (2D) homogeneous isotropic turbulence and the turbulent mixing layer are presented. In the homogeneous isotropic turbulence with tenth order of spatial accuracy, the power law of the inertial energy spectrum at the climax stage coincides with the predictions by Batchelor, Leith and Kraichnan. In the turbulent mixing layer with eight order of spatial accuracy, the vortex pairing are reproduced and the coherent structure of the Reynolds stress at the pairing is noticed.     

Viscoelastic flow in a curved channel: A similarity solution for the Oldroyd-B fluid

A similarity solution is used to analyse the flow of the Oldroyd fluid B, which includes the Newtonian and Maxwell fluids, in a curved channel modelled by the narrow annular region between two circular concentric cylinders of large radius. The solution is exact, including inertial forces. It is found that the non-Netonian kinematics are very similar to the Newtonian ones, although some stress components can become very large. At high Reynolds number a boundary layer is developed at the inner cylinder. The structure of this boundary layer is asymptotically analysed for the Newtonian fluid. Non-Newtonian stress boundary layers are also developed at the inner cylinder at large Reynolds numbers.     

Convergence acceleration by varying time‐step size using Bi‐CGSTAB method for turbulent flow computation

A design of varying step size approach both in time span and spatial coordinate systems to achieve fast convergence is demonstrated in this study. This method is based on the concept of minimization of residuals by the Bi‐CGSTAB algorithm, so that the convergence can be enforced by varying the time‐step size. The numerical results show that the time‐step size determined by the proposed method improves the convergence rate for turbulent computations using advanced turbulence models in low Reynolds‐number form, and the degree of improvement increases with the degree of the complexity of the turbulence models. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 454–474, 2001.     

Turbulence modeling from a new perspective

In this work we establish a link between a Reynolds averaged turbulence modeling methodology, containing interactions up to the second order correlations between the velocity fluctuations at various scales, and a multi-objective optimization problem with the constraints expressed in terms of equality and inequality, imposed by the given boundary conditions and the positive semi-definiteness of the Reynolds stress tensor, etc. The information unavailability and uncertainty associated with the boundary conditions for the fluctuation correlations of various orders is delineated, and the information from the Navier–Stokes equations is utilized to the extent allowed by the available input data necessary for simulations; turbulence from the perspective of systems simulation is explored and some objective functions are proposed. Finally, the challenges faced by the formulation and the issues yet to be resolved are discussed.     

Effect of constriction height on flow separation in a two-dimensional channel

We studied numerically the effect of the constriction height on viscous flow separation past a two-dimensional channel with locally symmetric constrictions. A numerically stable scheme in primitive variables (velocity and pressure) for the solution of two-dimensional incompressible time-dependent Navier–Stokes equations is employed using finite-difference approximation in staggered grid. The wall shear stresses at different heights of the constriction are computed and presented graphically. It is noticed that the maximum stress and the length of the recirculating region associated with two shear layers of the constriction increase with the increase of the area reduction of the constriction. The critical Reynolds number for symmetry breaking bifurcation for the 50%, 60% and 70% area reduction are obtained numerically. The flow field separates after the symmetry breaking bifurcation and the symmetry of the flow depends on the Reynolds number and the height of the constriction.     

2D numerical simulation of submerged hydraulic jumps

Hydraulic jumps are usually used to dissipate energy in hydraulic engineering. In this paper, the turbulent submerged hydraulic jumps are simulated by solving the unsteady Reynolds averaged Navier–Stokes equations along with the continuity equation and the standard k–? equations for turbulence modeling. The Lagrangian moving grid method is employed for the simulation of the free surface. In the developed model, kinematic free-surface boundary condition is solved simultaneously with the momentum and continuity equations, so that the water elevation can be obtained along with velocity and pressure fields as part of the solution. Computational results are presented for Froude numbers ranging from 3.2 to 8.2 and submergence factors ranging from 0.24 to 0.85. Comparisons with experimental measurements show that numerical model can simulate the velocity field, variation of free surface, maximum velocity, Reynolds shear and normal stresses at various stations with reasonable accuracy.