Finite-mode spectral model of homogeneous and isotropic Navier-stokes turbulence: a rapidly depleted energy cascade

An eddy-viscous term is added to Navier-Stokes dynamics at wave numbers k greater than the inflection point kc of the energy flux F(log(k)). The eddy viscosity is fixed so that the energy spectrum satisfies E(k) = E(kc) (k/kc)(-3) for k>kc. This resulting forcing induces a rapid depletion of the energy cascade at k>kc. It is observed numerically that the model reproduces turbulence energetics at k< or =kc and statistics of two-point velocity correlations at scales r>lambda (Taylor microscale). Compared to a direct numerical simulation of R(lambda) = 130 an equivalent run with the present model results in a gain of a factor 20 in CPU time.     

溃坝流的光滑粒子法模拟

用光滑粒子法模拟几种情况的溃坝流动,对坝外无水和有水的情况进行二维模拟,对坝外有立柱的情况进行三维模拟.流动控制方程采用雷诺平均方程模拟溃坝流动的湍流效应,采用混合长度形式的涡粘模式对控制方程进行封闭,推导其相应的光滑粒子形式的方程.模拟结果表明,数值模拟的流动特征与实验结果符合得非常好,说明发展的光滑粒子法有效.     

Shear viscosity of strongly coupled Yukawa systems on finite length scales

The Yukawa shear viscosity has been calculated using nonequilibrium molecular dynamics. Near the viscosity minimum, we find exponential decay consistent with the Navier-Stokes equation, with significant deviations on finite length scales for larger viscosity values. The viscosity is determined to be nonlocal on a scale length consistent with the correlation length, revealing the length scales necessary for obtaining transport coefficients in the hydrodynamic limit by nonequilibrium molecular dynamics methods. Our results are quasiuniversal with respect to excess entropy for excess entropies well below unity.     

环量控制翼型动态失速特性研究

本文用数值求解二维可压非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的方法研究了环量控制翼型的动态失速特性．揭示了环量控制翼型在以不同吹气动量系数、不同频率做俯仰振荡时的动态失速特征．指出吹气动量系数对环量控制翼型动态失速的特性有很大影响,同时,振荡频率的影响与传统翼型相比也有不同．进一步分析了其失速的机理．     

The energy spectrum of fully developed turbulence generated by random stirring

For fully developed turbulence in an incompressible fluid described by the Navier-Stokes equations with Gaussian random forces the relation between the energy spectrum and the stirring mechanism is investigated within a variational approach. Therein, the effect of nonlinear mode coupling is approximated by a wave number dependent eddy viscosity determined via a nonlinear integral equation for the energy spectrum. For various stirring spectra analytic approximations are compared with the solution obtained numerically with a cutoff in the integral kernel which ensures in eddy relaxing processes that the stirring forces exert strain only on scales larger than the eddy size. The results are compared with renormalization group calculations and closure approximations. Random forces injecting energy at a ratek ^–1 into the wave number banddk aroundk lead to a Kolmogorov distribution of energy. The spectrum of small-scale velocity fluctuations is shown to be universal in the sense that it remains unchanged under variations of the long wavelength stirring spectra.     

用多块多网格方法数值模拟三维粘性流动

本文给出了一个模拟三维粘性流动的数值方法．该方法来用高分辨率 ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式求解三维雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，使用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型估计湍流粘性系数，用多重网格技术加速收敛，采用多块结构化网格处理复杂的物理域．文中给出了叶轮机械多个叶片排和透平排汽缸内的全三维粘性流动的数值结果．     

大气边界层模式中随机参数的反演与不确定性分析

在大气边界层气象中湍流黏性系数是一个很重要的参数,通过直接观测往往无法得到其准确值,仅能通过间接观测获得大致范围.本文选用随机广义Ekman动力近似模式中的湍流黏性系数进行反演研究与不确定性分析.首先利用风速观测数据,并采用基于混沌多项式的集合Kalman滤波方法对系数进行反演,降低其不确定性,缩小可能取值的范围,该方法的核心思想是将集合Kalman滤波方法中求解模式不确定性传播的方法由蒙特卡罗法改为混沌多项式展开,从而避免大规模采样带来的计算资源耗费.然后进行数值实验,结果表明该方法能够有效且快速地求解出湍流黏性的后验概率分布,从而达到降低系数不确定性的目的.根据系数的先验分布计算出风速的先验分布,从而找到风速不确定性大的区域,且揭示了在不确定性大的区域内的观测数据进行系数反演可得到十分明显的效果,这对于观测点位置的选择提供了重要的指导.     

基于层流动能湍流模型的数值模拟方法

从层流动能出发提出数值模拟原则;综合考虑自然、旁路和分离流转捩的因素构建实用的层流动能湍流模型,结合预处理和基本求解技术发展出适于转捩流动的数值模拟方法和程序.针对预处理技术,以Weiss-Smith矩阵为基础,考虑湍流粘性的影响;针对基本离散格式和边界条件,结合模型方程进行对角占优强化等特殊处理.最后通过平板边界层和典型翼型,特别是低雷诺数翼型的数值模拟,验证数值方法的有效性和鲁棒性.算例表明本文的方法能够为求解更复杂的流动提供参考.     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

Large eddy simulation of a shocked gas cylinder instability induced turbulence

The Navier-Stokes equations for compressible fluid are solved with the operator splitting technique and LES (large eddy simulation) with the Smagorinsky model. A computational code MVFT (multi-viscosity-fluid and turbulence) is developed to study hydrodynamic instability and the induced turbulent mixing for multi compressible fluid. In order to validate the code MVFT,the LANL's shock tube experiment of shocked SF6 gas cylinder is simulated with the initial state of SF6 gas cylinder described by dissipative ...     

On the Dynamics of Navier-Stokes and Euler Equations

This is a detailed study on certain dynamics of Navier-Stokes and Euler equations via a combination of analysis and numerics. We focus upon two main aspects: (a) zero viscosity limit of the spectra of linear Navier-Stokes operator, (b) heteroclinics conjecture for Euler equation, its numerical verification, Melnikov integral, and simulation and control of chaos. Due to the difficulty of the problem for the full Navier-Stokes and Euler equations, we also propose and study two simpler models of them.     

On the structure of the velocity field under the free spheroidal surface of a viscous liquid drop oscillating in an electrostatic field

An asymptotic analytical solution to an initial boundary-value problem considering (i) the time evolution of the capillary oscillation amplitude as applied to a viscous spheroidal liquid drop placed in a uniform electrostatic field and (ii) the liquid flow velocity field inside the drop is found. The problem is solved in an approximation that is linear in two small parameters: the dimensionless oscillation amplitude and the dimensionless field-induced constant deformation of the equilibrium (spherical) shape of the drop. Terms proportional to the product of the small parameters are retained. In this approximation, interaction between oscillation modes is revealed. It is shown that the intensity of the eddy component of the oscillation-related velocity field depends on the liquid viscosity and the external uniform electrostatic field strength. The intensity of the eddy component decays rapidly with distance from the free surface. The depth to which the eddy flow (which is caused by periodical flows on the free surface) penetrates into the drop is a nonmonotonic function of the polar angle and increases with dimensionless viscosity and field strength.     

Convergence of the Vanishing Viscosity Approximation for Superpositions of Confined Eddies

A confined eddy is a circularly symmetric flow with vorticity of compact support and zero net circulation. Confined eddies with disjoint supports can be superimposed to generate stationary weak solutions of the two-dimensional incompressible inviscid Euler equations. In this work, we consider the unique weak solution of the two-dimensional incompressible Navier-Stokes equations having a disjoint superposition of very singular confined eddies as the initial datum. We prove the convergence of these weak solutions back to the initial configuration, as the Reynolds number goes to infinity. This implies that the stationary superposition of confined eddies with disjoint supports is the unique physically correct weak solution of the two-dimensional incompressible Euler equations.     

Stochastic Lagrangian Particle Approach to Fractal Navier-Stokes Equations

In this article we study the fractal Navier-Stokes equations by using the stochastic Lagrangian particle path approach in Constantin and Iyer (Comm Pure Appl Math LXI:330–345, 2008). More precisely, a stochastic representation for the fractal Navier-Stokes equations is given in terms of stochastic differential equations driven by Lévy processes. Based on this representation, a self-contained proof for the existence of a local unique solution for the fractal Navier-Stokes equation with initial data in \mathbb W^1,p{{\mathbb W}^{1,p}} is provided, and in the case of two dimensions or large viscosity, the existence of global solutions is also obtained. In order to obtain the global existence in any dimensions for large viscosity, the gradient estimates for Lévy processes with time dependent and discontinuous drifts are proved.     

Scale-wise coherent vorticity extraction for conditional statistical modeling of homogeneous isotropic two-dimensional turbulence

Classical statistical theories of turbulence have shown their limitations, in that they cannot predict much more than the energy spectrum in an idealized setting of statistical homogeneity and stationarity. We explore the applicability of a conditional statistical modeling approach: can we sort out what part of the information should be kept, and what part should be modeled statistically, or, in other words, “dissipated”? Our mathematical framework is the initial value problem for the two-dimensional (2D) Euler equations, which we approximate numerically by solving the 2D Navier-Stokes equations in the vanishing viscosity limit. In order to obtain a good approximation of the inviscid dynamics, we use a spectral method and a resolution going up to 8192². We introduce a macroscopic concept of dissipation, relying on a split of the flow between coherent and incoherent contributions: the coherent flow is constructed from the large wavelet coefficients of the vorticity field, and the incoherent flow from the small ones. In previous work, a unique threshold was applied to all wavelet coefficients, while here we also consider the effect of a scale by scale thresholding algorithm, called scale-wise coherent vorticity extraction. We study the statistical properties of the coherent and incoherent vorticity fields, and the transfers of enstrophy between them, and then use these results to propose, within a maximum entropy framework, a simple model for the incoherent vorticity. In the framework of this model, we show that the flow velocity can be predicted accurately in the L² norm for about 10 eddy turnover times.     

Zero Viscosity Limit for Analytic Solutions of the Navier-Stokes Equation on a Half-Space.¶ II. Construction of the Navier-Stokes Solution

This is the second of two papers on the zero-viscosity limit for the incompressible Navier-Stokes equations in a half-space in either 2D or 3D. Under the assumption of analytic initial data, we construct solutions of Navier-Stokes for a short time which is independent of the viscosity. The Navier-Stokes solution is constructed through a composite asymptotic expansion involving the solutions of the Euler and Prandtl equations, which were constructed in the first paper, plus an error term. This shows that the Navier-Stokes solution goes to an Euler solution outside a boundary layer and to a solution of the Prandtl equations within the boundary layer. The error term is written as a sum of first order Euler and Prandtl corrections plus a further error term. The equation for the error term is weakly nonlinear; its linear part is the time dependent Stokes equation. This error equation is solved by inversion of the Stokes equation, through expressing the solution as a regular (Euler-like) part plus a boundary layer (Prandtl-like) part. The main technical tool in this analysis is the Abstract Cauchy-Kowalewski Theorem. Received: 5 September 1996 / Accepted: 14 July 1997     

双大涡模拟方法模拟提升管内气固两相流动

基于Smagorinsky涡黏模型以及颗粒动理学理论,建立了气固两相流双大涡模拟模型。考虑大涡模拟中过滤尺度的影响,给出颗粒相亚格子压力和热传导系数计算模型。考虑颗粒聚团对两相作用的影响,给出了考虑颗粒聚团作用的气固两相多尺度曳力系数模型。数值模拟了提升管内气固两相流动特性,合理地预测出了提升管内气固两相环-核流动结构。模拟结果与Knowlton等实测结果相吻合。     

二维卡门涡街的格子Boltzmann仿真

格子气自动机和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新的方法。本文中，我们介绍了Boltzmann方法，解决了格子气方法中的缺点，通过选择适当平衡分布及其参数，导出了Navier－Stokers方程，并得到了声速和粘性系数。最后在微机上模拟了在无限长平板流动问题及绕单分离板的流动问题，得到了卡门涡街。结果表该模型有格子气方法及其它的数值方法所没有的优点，计算更精确、更直观、更有效。     

Remarks about the Inviscid Limit of the Navier–Stokes System

In this paper we prove two results about the inviscid limit of the Navier-Stokes system. The first one concerns the convergence in H ^s of a sequence of solutions to the Navier-Stokes system when the viscosity goes to zero and the initial data is in H ^s . The second result deals with the best rate of convergence for vortex patch initial data in 2 and 3 dimensions. We present here a simple proof which also works in the 3D case. The 3D case is new.