两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。     

Some developments in computational modeling of turbulent flows

In this paper, some recent developments of two turbulence closure schemes at ICOMP, NASA Lewis will be discussed. One is the Reynolds-stress algebraic equation model and the other is the Reynolds-stress transport equation model. Various model constraints required by the rapid distortion theory, the invariant theory and the realizability principle, etc. will be described in the model development. The models discussed are for high-turbulent Reynolds number flows, so that the near-wall turbulence and the low-Reynolds-number turbulence are not discussed here.     

A stability condition for turbulence model： From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale （EMMS） model for gas-solid systems, is formulated to close the dynamic constraint equa- tions of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the ＂EMMS-based turbulence model＂, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room, The numerical results show that the EMMS-hased turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas–solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the “EMMS-based turbulence model”, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

Pressure–Strain Correlation Modeling: Towards Achieving Consistency with Rapid Distortion Theory

In turbulence closure modeling, it is widely accepted that the rapid pressure–strain correlation (RPSC) model be consistent with the rapid distortion theory (RDT). It is desirable to achieve this consistency with a closure model that is computationally tractable and satisfies the requisite mathematical constraints of realizability and linearity in the appropriate variables. In this investigation, starting from a detailed modal analysis of two-dimensional mean flows, we identify important flow features to be incorporated into the model. However, the dynamical system analysis shows that the suggested physics cannot be embodied in a model with all desired computational and mathematical attributes. To resolve this conflict, we propose a slight compromise in the physical requirement and ease one of the linearity constraints leading to a “best possible” tractable model. Overall, the present work provides important insight into RPSC closure modeling challenges—arising from the interplay among physical fidelity, computational viability and mathematical constraints—and proposes avenues for future improvement.     

Modeling fluid-particle interaction in dilute-phase turbulent liquid-particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid-particle pipe flows in which the inter-stitial fluid effect on the particle fluctuating motion is significant.The impacts of employing different drag correlations and turbulence closure models to describe the fluid-particle interactions(i.e.drag force and long-range interaction)are examined at both the mean and fluctuating velocity levels.The model pre-dictions are validated using experimental data of turbulent liquid-particle flows in a vertical pipe at different particle Reynolds numbers(ReP > 400 and ReP < 400),which characterize the importance of the vortex shedding phenomenon in the fluid-phase turbulence modulation.The results indicate that(1)the fluctuating velocity level predictions at different ReP are highly sensitive to the drag correlation selec-tion and(2)different turbulence closure models must be employed to accurately describe the long-range fluid-particle interaction in each phase.In general,good agreement is found between the model predic-tions and the experimental data at both the mean and fluctuating velocity levels provided that appropriate combinations of the drag correlation and the turbulence closure model are selected depending on Rep.     

Modeling the Dynamics of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distortion theory (RDT). The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined. This revised version was published online in July 2006 with corrections to the Cover Date.     

关于植被中湍流的研究

本文从实验和模式理论两个方面介绍了近年来在植被湍流研究方面的进展.实验研究表明,植被内湍流是高度间歇性的,大尺度涡在湍流输运中起主导作用.植物枝、干、叶打碎了大尺度涡,产生的迹湍流中的小尺度涡更容易耗散成热。因此,植被内能量不按一般的能量级串理论预示的方式进行,而要发生所谓“短路”现象,使得功率谱曲线的斜率在惯性区比-2/3律更负.现有的高阶封闭模式虽然取得了很大成功,但在模拟湍流强度上高估了其大小,还需根据实验加以改进.      

Gas turbulence modulation in the pneumatic conveying of massive particles in vertical tubes

The present work is concerned with the interaction between large particles and gas phase turbulence. Gas turbulence modulation in these systems is considered to be dominated by a generation mechanism which arises due to the presence of wakes behind particles. Following a recent proposal, a closure for gas turbulence modulation accounting for the effect of wakes is employed within the context of a mathematical model for particle-laden, turbulent flows. The model accounts for particle particle and particle-wall interactions associated with larger particles based on concepts from gas kinetic theory. It is shown that due to the significant flattening of the mean gas velocity profile with the addition of particles, and the corresponding decrease in turbulent energy production, a generation mechanism must be present in order to produce gas velocity fluctuation predictions which are consistent with the experimental measurements, even in the case where the experimental results indicate a net suppression of gas phase turbulence in the presence of particles.     

简化Navier—Stokes方程的IFT理论的各向异性张力有限元法

本文首先讨论简化Navier-Stokes IFT方程组的有限元离散方式,然后对其广义解进行分析,并从而利用与之相匹配的各向异性张力单元对流函数—涡量方程进行计算。通过平板层流和台阶绕流两个算例的分析,证明这种与IFT理论相匹配的有限单元算法是成功的。     

高阶矩湍流模型研究进展及挑战

高阶矩模型是湍流模式理论研究中的难点和前沿.自周培源先生首次建立一般湍流的雷诺应力输运方程起,为了更精确的预测复杂流动,人们从未间断过对高阶矩模型的研究.尤其进入新世纪以来,随着计算机硬件水平的飞跃和高精度数值算法的突破,湍流模拟方法正由RANS向LES转变.而无论对于RANS框架、LES框架还是两者混合,高阶矩模式都...     

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state,and the flow is frequently unsteady and turbulent.To assess the state-of-the-art of computational capabilities for unsteady cavitating flows,different cavitation and turbulence model combinations are conducted.The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics.The kε turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively.We have also blended these alternative cavitation and turbulence treatments,to illustrate that the eddy viscosity near the closure region can significantly influence the capture of detached cavity.From the experimental validations regarding the force analysis,frequency,and the cavity visualization,no single model combination performs best in all aspects.Furthermore,the implications of parameters contained in different cavitation models are investigated.The phase change process is more pronounced around the detached cavity,which is better illus-trated by the interfacial dynamics model.Our study provides insight to aid further modeling development.     

基于颗粒湍流和颗粒碰撞相互作用规律的湍流模型的研究

颗粒湍流和颗粒碰撞的相互作用规律是两相流动中的核心问题。用颗粒湍流模型和颗粒碰撞的动力论模型叠加的方法在研究两相湍流流动方面取得了一定的成效,但是还有待改进。本文基于颗粒湍流形成大尺度脉动和颗粒间碰撞引起小尺度脉动的概念,从双流体模型出发,建立了两相流动的双尺度kp-pε两相湍流模型。利用该模型对下行床和突扩室内的气固...     

涡度拟能传输惯性区二维湍流标量场方差谱

本文应用非平衡统计力学封闭方法,给出涡度拟能传输惯性区二维湍流标量场方差谱的完整表达式,数值计算该表达式中的比例系数B。由于传输过程的非局部性,B依赖于表征波数变化范围的局部化因子,不再是普适常数。     

On the viscosity and thermal conduction of fluids with multivalued internal energy

This work is concerned with an extension of the classical compressible Euler model of fluid dynamics in which the fluid internal energy is a measure-valued quantity. This model can be derived from the hydrodynamic limit of a kinetic model involving a specific class of collision operators. In the present paper, we investigate diffusive corrections of this fluid dynamical model derived from a Chapman–Enskog expansion of the kinetic model, in the case where the collision time depends on the particle energy in the fluid frame. We show that the closure relations for the stress tensor and heat flux vector differ from their expression in the usual Navier–Stokes model. We argue why such a feature could be used as a tool towards an understanding of fluid turbulence from kinetic theory.     

On the Interaction of Turbulent Shear Layers with Harmonic Perturbations

The problem of coherent perturbations in a turbulent shear layer is considered for the purpose of developing a mathematical model based on a triple decomposition that extracts the coherent components of random fluctuations. The governing equations for the mean and the coherent parts of flow are derived, assuming the eddy-viscosity equivalence for the random part of flow, and solved by iterations to provide a coupled solution of the problem as a whole. Calculations agree well with experimental data in the upstream part of the layer where the mean–coherent flow interaction is the most important. In this region, the interaction changes the mean flow velocity distribution in such a manner that the neutral stability curve is shifted upstream relative to its position in the undisturbed layer and the perturbation intensity decreases further downstream. Experiments show that the coherent waves suppress the turbulent Reynolds stress production downstream of this region, but the model fails to predict the layer spreading correctly probably due to an inadequate turbulence closure of the mean flow. For the case of a turbulent mixing layer, we suggest a new closure relation that takes into account this coherent-random interaction.     

小波变换在湍流数值研究中的应用

小波变换具有时空双局部性特点,恰好适应了湍流特性。本文主要阐述了小波在湍流数值计算中的两大研究进展:一个是利用连续小波(高斯小波)可使L ap lace算子降阶的特性来求解N-S方程;另一个是正交小波与有限元法相结合的方法——相干涡模拟。指出它们的优缺点及其存在的问题,并对小波在湍流计算中的应用前景作了展望。     

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.     

High-order strand grid methods for shock turbulence interaction

ABSTRACT

In this work, we examine the flux correction method for three-dimensional transonic turbulent flows on strand grids. Building upon previous work, we treat flux derivatives along strands with high-order summation-by-parts operators and penalty-based boundary conditions. A finite-volume like limiting strategy is implemented in the flux correction algorithm in order to sharply capture shocks. To achieve turbulence closure in the Reynolds-Averaged Navier–Stokes equations, a robust version of the Spalart–Allmaras turbulence model is employed that accommodates negative values of the turbulence working variable. Validation studies are considered which demonstrate the flux correction method achieves a high degree of accuracy for turbulent shock interaction flows.