Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.     

A New Low-Reynolds-Number One-Equation Model of Turbulence

In this study, we propose a new Low-Reynolds-Number (LRN)one-equation model, which is derived from an LRN two-equation(k-ε) model. The derivation of the transport equation, in principle, is based on the assumption that the turbulent structure parameter remains constant. However, the relation for the turbulent structure parameter a ₁(=|− |/k) is modified to account for near-wall turbulence. As a result, the present one-equation model contains a term which takes the near-wall limiting behavior explicitly into account. Thus, the present model provides the correct wall-limiting behavior of turbulence in the vicinity of the wall and can be applied to the analysis of heat transfer. The validity of the present model is tested in channel flows, boundary layer flows with and without pressure gradient, plane wall jet, and flow with separation and reattachment. The calculated results showed good agreement with the direct numerical simulation (DNS) and experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Note on the Numerical Treatment of the k-epsilon Turbulence Model*

Numerical solution of the equations arising from the κ mdash; ε turbulence model has difficulties inherent to nonlinear convection-reaction-diffusion equations with strong reaction terms, resulting in that numerical schemes easily become unstable. We present a formulation that stresses on the robustness of the solution method, tackling common problems that produce instability. The main contribution concerns the loss of positivity of κ and ε, which is addressed by acting on the coefficients of the reaction and diffusion terms rather than on the turbulent variables themselves. In addition, a linearized implicit, non-iterative, treatment of the wall law is proposed.     

Turbulence Closure with a Topography-parameter-free Single Equation Model

A new topography-parameter-free turbulence closure, based on a transport equation for the pseudo eddy viscosity, is described. The model is tested against experimental data of several flow cases across the Mach number range and compared to Menter's single equation model and to the Spalart-Allmaras model. The main conclusion is that the new closure outperforms the other two models in both low and high speed flows.     

A Turbulence Model Sensitivity Study for CH4/H2 Bluff-Body Stabilized Flames

The objective of this work is to assess the performances of different turbulence models in predicting turbulent diffusion flames in conjunction with the flamelet model.The k– model, the Explicit Algebraic Stress Model (EASM) and the k– model withvaried anisotropy parameter C (LEA k– model)are first applied to the inert turbulent flow over a backward-facing step, demonstrating the quality of the turbulence models. Following this, theyare used to simulate the CH₄/H₂ bluff-body flame studied by the University of Sydney/Sandia.The numerical results are compared to experimental values of the mixture fraction, velocity field, temperature and constituent mass fractions.The comparisons show that the overall result depends on the turbulence model used, and indicate that theEASM and the LEA k– models perform better than the k– model and mimic most of the significant flow features.     

Efficiency of Scale-Similarity Model for Study of Forced Compressible Magnetohydrodynamic Turbulence

Efficiency of scale-similarity model for study of forced compressible magnetohydrodynamic turbulence is studied. The scale-similarity model has several important advantages in contrast to the eddy-viscosity subgrid closures: good reproduction of the correlation between actual and model turbulent stress tensor even when the flow is highly anisotropic, and absence of special model constants. These advantages may be very essential for study of forced magnetohydrodynamic turbulence. Numerical computations under various similarity parameters are carried out and the obtained results are analyzed by means of comparison with results of direct numerical simulation and Smagorinsky closure for magnetohydrodynamics. Linear forcing algorithm is applied to keep the characteristics of turbulence stationary in time. Influence of discrete filter shapes on the scale-similarity model is studied as well. It is shown that the scale-similarity model provides good accuracy and the results agree well with the direct numerical simulation results. The present results show that the scale-similarity model might be a useful subgrid closure for study of scale-invariance properties of forced compressible magnetohydrodynamic turbulence in the inertial range and in contrast to decaying case the scale-similarity model can serve as a stand alone subgrid model.     

A Quasi-Realizable Cubic Low-Reynolds Eddy-Viscosity Turbulence Model with a New Dissipation Rate Equation

A non-linear relationship of the Reynolds stresses in function of the strain rate and vorticity tensors, with terms up to third order, is developed. Anisotropies in the normal stresses, influence from streamline curvature or rotation of the reference frame, and swirl effects are accounted for. The relationship is linked to ak–ε model with a modified transport equation for the dissipation rate. A new low-Reynolds source term is introduced and a model parameter is written in terms of dimensionless rate-of-strain and vorticity. The model is checked on different realizability constraints. It is shown that practically all constraints are fulfilled. The model is numerically tested on a fully developed channel and pipe flow, both stationary and rotating. The plane jet–round jet anomaly is addressed. Finally, the model is applied to the flow over a backward-facing step. Results are compared with a linear low-Reynolds k–ε model and the shear stress transport model. This revised version was published online in July 2006 with corrections to the Cover Date.     

一种紊流润滑理论分析新方法—复合型紊流模式理论

在对主要用于高速轻载工况的常用紊理论进行简要分析和充分考察润滑流场的边界条件及内部结构的基础上，采用理论上比现有紊流润滑理论更为合理的复合型紊流模式理论，即在近壁区采用低紊流雷诺数的ｋ－ε模式，而在紊流核心区采用代数雷诺应力模式，对复杂流场的紊流润滑进行了分析，同时在计入惯性效应的情况下，推导出了一种适用于高压密封和高速重载轴承等计算，计算结果与实验数据十分吻合，验证了模型的有效性，可以应用于高压     

Effects of Spray and Turbulence Modelling on the Mixing and Combustion Characteristics of an n-heptane Spray Flame Simulated with Dynamic Adaptive Chemistry

Accurate modelling of spray combustion process is essential for efficiency improvement and emissions reduction in practical combustion engines. In this work, both unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and large eddy simulations (LES) are performed to investigate the effects of spray and turbulence modelling on the mixing and combustion characteristics of an n-heptane spray flame in a constant volume chamber at realistic conditions. The non-reacting spray process is first simulated with URANS to investigate the effects of entrainment gas-jet model on the penetration characteristics and fuel vapor distributions. It is found that the droplet motion near the nozzle has significant influence on the fuel vapor distribution, while the liquid penetration length is controlled by the evaporation process and insensitive to gas-jet model. For the case considered, both URANS with the gas-jet model and large eddy simulations can properly predict the vapor penetration. For the combustion characteristics, it is found that LES yields better predictions in the global combustion characteristics. The URANS with gas jet model yields a comparable flame length and lift-off-length (LOL) to LES, but results in a larger ignition delay time compared to the experimental data. Another focus of this work is to qualify the convergence characteristics of the dynamic adaptive chemistry (DAC) method in these transient combustion simulations, where DAC is applied to reduce the mechanism locally and on-the-fly to accelerate chemistry calculations. The instantaneous flame structures and global combustion characteristics such as ignition delay time, flame lift-off length and emissions are compared between simulations with and without DAC. For URANS, good agreements are observed both on instantaneous flame structures and global characteristics. For LES, it is shown that the errors incurred by DAC are small for scatter distributions in composition space and global combustion characteristics, while they may significantly affect instantaneous flame structures in physical space. The study reveals that for DAC application in transient simulations, global or statistic information should be used to assess the accuracy, such as manifolds in composition space, conditional quantities and global combustion characteristics. For the cases investigated, a speed-up factor of more than two is achieved by DAC with a 92-species skeletal mechanism with less than 0.2 % and 3.0 % discrepancy in ignition delay and LOL, respectively.     

Boundary Layer for the Navier-Stokes-alpha Model of Fluid Turbulence

We study boundary-layer turbulence using the Navier-Stokes-alpha model obtaining an extension of the Prandtl equations for the averaged flow in a turbulent boundary layer. In the case of a zero pressure gradient flow along a flat plate, we derive a nonlinear fifth-order ordinary differential equation, which is an extension of the Blasius equation. We study it analytically and prove the existence of a two-parameter family of solutions satisfying physical boundary conditions. Matching these parameters with the skin-friction coefficient and the Reynolds number based on momentum thickness, we get an agreement of the solutions with experimental data in the laminar and transitional boundary layers, as well as in the turbulent boundary layer for moderately large Reynolds numbers.     

主动移频式动力吸振器及其动力特性的研究

动力吸振器是振动控制中比较有效的减振装置,只要吸振器(子系统)的振动固有频率与振动物体(主系统)的振动频率相同,即可有效地消除主系统的振动。但传统动力吸振器的控制频率带宽较窄,限制了其稳定性和减振效果的提高。本文通过独特的机械设计,研制了一种可以通过调节自身的几何参数,使得其固有频率随几何参数线性变化的主动移频的新型动力吸振器,并初步设定了相应的控制方法。文中还对其动力学特性进行了理论分析和实验测试,分析了它的机理,评估了它的实际减振效果。研究结果表明该吸振器可以大范围调节自身固有频率,有效拓宽吸振频带,具有良好的减振性能和稳定性。     

正交曲线坐标系下紊流数学模型的曲率修正

考虑弯曲边界曲率效应对水流水力特性的影响,建立了正交曲线坐标系下的素流数学模型。通过计算实例说明,该数学模型能够很好地帷有复杂边界的流线弯曲水流的水力特性。     

Computation of Compressible Flows using a Two-equation Turbulence Model on Unstructured Grids

This paper presents a numerical method for solving compressible turbulent flows using a k - l turbulence model on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multi-stage Runge-Kutta time stepping scheme, while the turbulence equations are advanced using a multi-stage point-implicit scheme. The positivity of turbulence variables is achieved using a simple change of dependent variables. The developed method is used to compute a variety of turbulent flow problems. The results obtained are in good agreement with theoretical and experimental data, indicating that the present method provides a viable and robust algorithm for computing turbulent flows on unstructured meshes.     

Coherent Vortex Simulation (CVS), A Semi-Deterministic Turbulence Model Using Wavelets

In the spirit of Ha Minh's semi-deterministic model, we propose a new method for computing fully-developed turbulent flows, called Coherent Vortex Simulation (CVS). It is based on the observation that turbulent flows contain both an organized part, the coherent vortices, and a random part, the incoherent background flow. The separation into coherent and incoherent contributions is done using the wavelet coefficients of the vorticity field and the Biot–Savart kernel to reconstruct the coherent and incoherent velocity fields. The evolution of the coherent part is computed using a wavelet basis, adapted at each time step to resolve the regions of strong gradients, while the incoherent part is discarded during the flow evolution, which models turbulent dissipation. The CVS method is similar to LES, but it uses nonlinear multiscale band-pass filters, which depend on the instantaneous flow realization, while LES uses linear low-pass filters, which do not adapt to the flow evolution. As example, we apply the CVS method to compute a time developing two-dimensional mixing layer and a wavelet forced two-dimensional homogeneous isotropic flow. We also demonstrate how walls or obstacles can be taken into account using penalization and compute a two-dimensional flow past an array of cylinders. Finally, we perform the same segmentation into coherent and incoherent components in a three-dimensional homogeneous isotropic turbulent flow. We show that the coherent components correspond to vortex tubes, which exhibit non-Gaussian statistics and long-range correlation, with the same k ^−5/3power-law energy spectrum as the total flow. In contrast, the incoherent components correspond to an homogeneous random background flow which does not contain organized structures and presents an energy equipartition together with a Gaussian PDF of velocity. This justifies their elimination during the CVS computation to model turbulent dissipation. This revised version was published online in August 2006 with corrections to the Cover Date.     

直线振动筛相似模拟动应力研究

对27m2大型直线振动筛的相似试验模型进行了动应力测试,得到了试验模型筛筛体梁、板关键点处的最大动应力应变值。根据量纲分析法和方程分析法推导的工作状态参数相似关系,通过对原型筛与缩小的相似试验模型筛的试验结果对比分析,验证了相似模型与原型试验数据的一致性及可靠性,证实了利用缩小的相似试验模型进行动态测试及分析的有效性。本文结果为提高直线振动筛的使用寿命及结构的优化设计提供了理论依据,也为相似试验模型的发展研究提供了参考。     

高层建筑横风向湍流脉动风压动力反应计算

根据湍流脉动压力的产生机理,从湍流理论的基本方程出发,根据Taylor关于湍流的"冻结"假定,导出了湍流脉动风压谱密度函数的解析计算公式.在分析过程中考虑了压力方程源项中全部"紊动-剪切"项的影响.若取用合适的湍流积分尺度,则由此公式得出的横风向脉动风压谱密度函数值与足尺观测数据相吻合,因此对以前的研究成果有一定的改进.由于在接近结构第一阶自振频率时,谱函数值仍处于较高的水平,因此根据湍流脉动风压谱密度函数计算得出的结构横风向风振动力反应位移值与加速度值均高于由日本规范等公式中规定的漩涡脱落扰力引起的反应值,而且推测随结构高度的增加这一影响也有增加的趋势.建议对此问题进行更深入的研究并在高层建筑结构设计时考虑这一差别的影响.     

Mathematical Model of the Unsteady Motion of a Shaft in a Hydrodynamic Plain Bearing

An asymptotic solution is obtained that describes the unsteady motion of a shaft in a cylindrical plain bearing with hydrodynamic lubrication in the case of a constant external load. Oscillatory modes of transition to a steady-state position of the shaft for various values of the external load are considered. The characteristic time of velocity relaxation to the quasiequilibrium values determined from the inertialess approximation equations is obtained. Oscillation frequencies and amplitudes, shaft paths, and oscillation decay times are determined. The effect of a thin elastic liner on the characteristics of the transient process is explored.     

Normal Momentum Transfer Study by a Dynamic Technique

A novel dynamic technique for normal momentum accommodation study in a gas–solid body system is developed. The method can also be used for pressure measurements and adsorption studies. The experimental setup is described. The data obtained are interpreted with allowance for geometric parameters. The normal momentum accommodation coefficients are estimated on the basis of measurement results in systems hydrogen (nitrogen) gas mixture–singlecrystal silicon.     

A Dynamic Two-Phase Pore-Scale Model of Imbibition

We present a dynamic pore-scale network model of imbibition, capable of calculating residual oil saturation for any given capillary number, viscosity ratio, contact angle, and aspect ratio. Our goal is not to predict the outcome of core floods, but rather to perform a sensitivity analysis of the above-mentioned parameters, except from the viscosity ratio. We find that contact angle, aspect ratio, and capillary number all have a significant influence on the competition between piston-like advance, leading to high recovery, and snap-off, causing oil entrapment. Due to significant CPU-time requirements we did not incorporate long-range correlations among pore and throat sizes in our network, but were limited to small-range correlations. Consequently, the gradual suppression of snap-off occurs within one order of magnitude of the capillary number. At capillary numbers around 10⁸ - 10⁷ snap-off has been entirely inhibited, in agreement with results obtained by Blunt (1997) who used a quasi-static model. For higher aspect ratios, the effect of rate and contact angle is more pronounced.