基于时空最优低维动力系统的多尺度可压缩湍流数值模拟方法（英文）

按照周培源教授关于研究湍流数值模拟建模时必须分析和求解脉动速度场的思想,该研究基于第一性原理,系统地建立了基于时空低维最优动力系统的多尺度可压缩湍流数值模拟方法(LMS方法),并将其应用于多次冲击Richtmyer-Meshkov问题的数值模拟中,首次得到了可压缩湍流的中尺度流场和不同于DNS近似解的湍流近似解.数值结果表明,LMS方法可以用较少的网格获得更精确的湍流近似解.首先解决了研究中遇到的几个问题,为LMS方法的构建铺平了道路.这些问题是:基于湍流的物理特性,提出了湍流大、中、小尺度分解的新概念;找到了box滤波空间相关性的计算方法;指出了湍流建模理论中长期存在的逻辑错误,提出了多尺度湍流模型的概念;讨论了湍流封闭问题的本质和关键,给出了克服湍流封闭问题的数值方法.采用box滤波方法/空间网格平均方法且在大尺度网格的意义下,LMS方法的本质是一种将RANS、LES、DES和DNS等湍流数值模拟方法统一的全新湍流数值模拟方法.需要指出的是,LMS方法也可以作为湍流模型研究的辅助工具,以检验SGS尺度方程/脉动方程中各项所对应的湍流模型是否正确.     

二维剪切流的粘性-无粘湍流干扰理论

对二维不可压缩近壁剪切湍流,本文提出一个粘性-无粘湍流干扰理论.主要内容有:从分子粘性考虑出发确定干扰湍流的流动结构及其物理尺度,导出空间为小尺度的局部流动结构随顺流距离的演变规律,导出支配干扰湍流流动的简化Reyno-lds(SR)方程和扩散抛物化K-ε方程.该SR方程是作者早先提出的简化Navier-Stokes(SNS)方程的湍流形式,它的重要性质是“简化运算”和时间Reynolds平均运算的顺序可以交换.关于最大湍流剪应力、本理论计算值与实验测量值很好相符.经典湍流边界层理论、Clauser平衡湍流边界层以及湍流分离Triple-deck理论均是本文理论的特例.证实了顺流方向长度尺度随干扰增强而显著减小的实验结论.     

反映强流动曲率效应的非线性湍流模型

首先定性地分析了流线曲率效应对流场湍流结构的影响,然后以U型槽道流为典型算例,对多种湍流模型进行了评估．评估的模型包括:线性涡粘性模型,二阶和三阶非线性涡粘性模型,二阶显式代数应力模型和Reynolds应力模型．评估结果表明,性能良好的三阶非线性涡粘性模型,如黄于宁等人发展的HM模型以及CLS模型,可以较好地描述流线的曲率效应对湍流结构的影响,如凸曲率作用下内壁附近湍流强度的衰减和凹曲率作用下外壁附近湍流的增强,并且较好地确定了管道下游的分离点位置和分离泡长度,其预测的结果和实验符合较好,与Reynolds力模型的结果十分接近,因此可以较好地应用于具有曲率效应的工程湍流的计算．     

分层剪切湍流大涡模拟的一种动力学亚格子尺度模型

基于Yoshizawa涡粘性模型提出了一种适用于热分层剪切湍流大涡模拟的动力学亚格子尺度 (SGS)模型 ,包括SGS湍流应力和湍流热通量模型 ,计算验证了该模型的正确性 .进一步采用大涡模拟方法和该动力学SGS模型研究了稳定分层和不稳定分层槽道湍流的湍流特性 ,以及湍流主要变量关联量的变化规律 .同时 ,根据计算结果 ,预测了热分层槽道湍流的临界Richardson数 ,与理论分析结果基本相符 .     

一种各向异性多重尺度的湍流模型

本文给出了一种适用于复杂湍流流动计算的各向异性、多重尺度的湍流模型(MS/ASM).这种模型对雷诺应力进行直接的模拟,并可模拟湍流流动的多重尺度影响.对自由应力流动、旋转流动和回流的湍流流动的计算表明,它比常用的单重尺度的k-ε模型有明显的改进.由于计算机工作量增加得不多,所以它在工程计算中,具有广泛的应用前景.     

导电薄板的磁弹性组合共振分析

基于Mexwell方程,给出了导电薄板的非线性磁弹性振动方程、电动力学方程和电磁力表达式.在此基础上,研究了横向磁场中梁式导电薄板的磁弹性组合共振问题,应用Galerkin法导出了相应的非线性振动微分方程组.利用多尺度法进行求解,得到了系统稳态运动下的幅频响应方程,分析了组合共振激发的条件.根据Liapunov近似稳定性理论,对稳态解的稳定性进行了分析,得到了稳定性的判定条件.通过数值计算,给出了一、二阶模态下共振振幅随调谐参数、激励幅值和磁场强度的变化规律曲线图,以及系统振动的时程响应图、相图、Poincare映射图和频谱图,进一步分析了电磁、机械等参量对解的稳定性及分岔特性的影响,并讨论了系统的倍周期和概周期等复杂动力学行为.     

应用于非惯性系湍流模拟的扩展内禀旋转张量

研究扩展内禀旋转张量在非惯性系湍流模拟中的作用,特别是对代数Reynolds应力湍流模式(如非线性K-ε模式)的重要性．为此,采用几个近年来发展的非线性K-ε湍流模式模拟旋转坐标系下均匀剪切湍流,并且和大涡模拟的结果进行比较．计算结果和分析表明,需要发展更先进的非线性K-ε模式从而更好地描述非惯性系下的复杂湍流．     

二阶动态亚格子尺度应力模型

提出了一个基于亚格子尺度应力与速度梯度张量之间关系的二阶动态模型.然后利用在高雷诺数的流场直接数值解的结果对此二阶模型进行检验.直接数值解的流场包括均匀各向同性强迫湍流,衰减湍流以及均匀旋转湍流.数值检验结果发现与一阶动态模型相比,二阶模型的相关系数提高.     

黏弹性屈曲梁非线性内共振稳态周期响应

研究了内共振下简支边界屈曲黏弹性梁受迫振动稳态周期幅频响应.考虑Kelvin黏弹性本构关系,并通过对非平凡平衡位形做坐标变换,建立屈曲梁横向振动的非线性偏微分-积分模型.基于对控制方程的Galerkin截断,得到多维非线性常微分方程组.在前两阶模态内共振存在的条件下,运用多尺度法分析截断后的控制方程,利用可解性条件消除长期项,获得一阶主共振下的幅值与相角方程.通过数值算例以展示系统稳态幅频响应关系以及失稳区域,从而聚焦系统共振中存在的非线性现象,如跳跃现象、滞后现象,并讨论了双跳跃现象随轴向荷载的演化.通过直接数值方法处理截断方程,数值验证近似解析解,计算结果表明多尺度法具有较高精度.     

面内平动黏弹性板非线性振动的内-外联合共振

研究了同时存在受迫共振和1∶3内共振时的面内平动黏弹性板的横向非线性振动问题．板的黏弹性材料用Kelvin本构关系描述．基于系统的运动方程和四边简支的边界条件,对偏微分方程应用直接多尺度法建立了联合共振时的可解性条件．应用Routh-Hurvitz判据对系统幅频响应的稳定性进行了判别．给出了黏弹性系数、面内平动速度和激励幅值3个参数对幅频响应的影响．最后,应用微分求积数值方法验证了近似解析方法的结论.     

A rheological model for fluid flows with arbitrary Reynolds numbers

An approach to constructing a phenomenological model of motion for a viscous fluid alternative to the Prandtl mixing-length hypothesis is suggested. The approach makes it possible to describe the motion of a fluid independently of the regime realized in a given region of the flow. On the basis of this approach, a differential one-parameter model for the flow of a viscous fluid applicable to any regimes of motion, called a uniform laminar-turbulent model, is constructed. For this purpose, the field of a scalar turbulence measure is introduced, which equals the ratio of the Reynolds stress to the total stress in the case of a simple shear flow. This makes it possible to write new expressions for turbulent viscosity. The influence of the turbulence measure field on the flow is taken into account by using an additional transport differential equation. The model is applicable to both compressible and incompressible fluids and makes it possible to obtain solutions in quadratures for steady simple shear flows. Various forms of the system of equations of motion and boundary conditions are given. Original Russian Text ￠ V.A. Pavlovskii, D.V. Nikushchenko, 2009, published in Vestnik Sankt-Peterburgskogo Universiteta. Seriya 1. Matematika, Mekhanika, Astronomiya, 2009, No. 1, pp. 104–112.     

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

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

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

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

One method for solving systems of nonlinear partial differential equations

A method for reducing systems of partial differential equations to corresponding systems of ordinary differential equations is proposed. A system of equations describing two-dimensional, cylindrical, and spherical flows of a polytropic gas; a system of dimensionless Stokes equations for the dynamics of a viscous incompressible fluid; a system of Maxwell’s equations for vacuum; and a system of gas dynamics equations in cylindrical coordinates are studied. It is shown how this approach can be used for solving certain problems (shockless compression, turbulence, etc.).     

Numerical simulation of unsteady flows with transient regimes

Self-oscillatory flows in aerodynamics and astrophysics are studied. The two-dimensional compressible gas equations are solved using the implicit Runge-Kutta scheme of the third order with respect to the inviscid terms and of the second order with respect to the viscous terms. An algebraic Cebeci-Smith turbulence model is used. Weakly unsteady and strongly unsteady flow regimes are observed. The former occur in a supersonic flow past a cylinder with a front projection and in the heliosphere. Such flows became stable when the turbulent diffusion is taken into account. The latter flows occur when a supersonic jet meets an obstacle and when such a jet penetrates a cavity. In these flows, the amplitude of oscillations slightly decreases when the turbulent diffusion is taken into account.     

Numerical simulation of airfoil vibrations induced by turbulent flow

The subject of this paper is the numerical simulation of the interaction between two-dimensional incompressible viscous flow and a vibrating airfoil. A solid elastically supported airfoil with two degrees of freedom, which can rotate around the elastic axis and oscillate in the vertical direction, is considered. The numerical simulation consists of the stabilized finite element solution of the Reynolds averaged Navier–Stokes equations with algebraic models of turbulence, coupled with the system of ordinary differential equations describing the airfoil motion. Since the computational domain is time dependent and the grid is moving, the Arbitrary Lagrangian–Eulerian (ALE) method is used. The developed method was applied to the simulation of flow-induced airfoil vibrations.     

两相流基本方程

本文按照严格的办法,从纳维斯笃克斯方程出发,经过雷诺平均得到了两相流动的雷诺方程组,从更大的程度上包含了一些物理量的湍流脉动乘积的平均项。本文还讨论了两相之间的相互作用力,分散相的压力和压力以外的应力以及两相之间能量交换项的表达式。     

Multiscale behavior and fractional kinetics from the data of solar wind–magnetosphere coupling

Multiscale phenomena are ubiquitous in nature as well as in laboratories. A broad range of interacting space and time scales determines the dynamics of many systems which are inherently multiscale. In many systems multiscale phenomena are not only prominent, but also they often play the dominant role. In the solar wind–magnetosphere interaction, multiscale features coexist along with the global or coherent features. Underlying these phenomena are the mathematical and theoretical approaches such as phase transitions, turbulence, self-organization, fractional kinetics, percolation, etc. The fractional kinetic equations provide a suitable mathematical framework for multiscale behavior. In the fractional kinetic equations the multiscale nature is described through fractional derivatives and the solutions of these equations yield infinite moments, showing strong multiscale behavior. Using a Lévy flights approach, we analyze the correlated data of the solar wind–magnetosphere coupling. Based on this analysis a model of the multiscale features is proposed and compared with the solutions of diffusion-type equations. The equation with fractional spatial derivative shows strong multiscale behavior with infinite moments. On the other hand, the equation with space dependent diffusion coefficients yield finite moments, indicating Gaussian type solutions and absence of long tails typically associated with multiscale behavior.