采用双分布函数格子Boltzmann方法的直流氩等离子体流动与传热分析

本文采用双分布函数格子Boltzmann方法(DDF-LB)模拟了直流氩电弧等离子体中的流动与传热问题。针对非局域热平衡状态下热等离子体中的原子、离子和电子三种组分,推导给出了不同组分的格子Boltzmann方程,通过耦合迭代求解得到各组分的温度和速度分布。文中采用随温度变化的松弛时间用以提高数值计算的稳定性。通过将本文的计算结果与文献结果进行比较,验证了DDF-LB方法求解等离子体中能量输运问题的稳定性和数值精度。     

尾流的双尺度简化模型及其大涡模拟建模

叶轮机中存在着复杂的湍流流动结构,如何准确地模拟湍流流动一直是工程界关注的核心问题。由于湍流在大部分区域都处于强非均衡态,根据均衡能量传输假设建立起来的RANS湍流模型和大涡模拟亚格子模型都不再适用。非均衡湍流的一个经典的模型流动是尾流,但通常意义上的尾流往往是无数不同尺度相互作用的结果,难以进行简化分析。本文构造了一个极度简化的双尺度尾流模型进行了直接数值模拟(DNS),结果显示在转捩前期就存在反向的非均衡能量传输。分别采用各种传统的RANS湍流模型和大涡模拟亚格子模型进行了数值模拟,发现各种模型都无法准确预测尾流转捩。在考虑非均衡能量传输机制后,通过流向速度梯度扭率构造的理性亚格子应力模型可以使大涡模拟模拟的转捩位置相对于标准Smagorinsky模型有所改善;不采用涡粘假设的改进速度增量(IVI)模型可以得到更加接近DNS的结果。这些结果为湍流模型进一步的改进提供了思路。     

一种改进后的亚格子特征尺度表达式

针对大涡模拟涡粘性亚格子模式中使用网格尺度为特征尺度存在的不足,提出一种改进的亚格子特征尺度表达式,利用结合Sagaut混合尺度模式获得的改进模式与Smagorinsky模式、Germano模式和Sagaut模式等进行对比研究.时间发展混合层和中性大气边界层的大涡模拟结果表明,改进的有效亚格子特征尺度能合理反映亚格子脉动的时空分布特性,同时扩展了适用范围.在对耗散程度、流动演化过程、可解湍流强度和可解雷诺应力等方面的模拟,改进模式优于基准模式,表明亚格子脉动有效特征尺度在大涡模拟中具有一定的理论基础和广泛的应用价值.     

R-T界面不稳定性及湍流混合的大涡模拟

采用大涡模拟(LES)方法研究R-T(Rayleigh Taylor)不稳定性,显示了不稳定发展的小扰动阶段、变形阶段、规则非线性阶段、不规则非线性阶段和湍流混合阶段界面演化与特征,分析了湍流混合层厚度增长的规律,同时还研究了湍流二阶相关量如动量、质量和能量通量在大尺度和亚格子尺度所占分量,验证了LES方法是处理R-T界面不稳定性和湍流混合的一种有效方法.     

低阶格式在大涡模拟计算中的适用性

采用三维Taylor-Green涡作为研究对象,利用工程中常用的低阶数值格式,研究格式本身的数值误差对大涡模拟计算的影响.结果表明:三种数值格式的数值耗散行为都与亚格子模型行为类似,即在小雷诺数下,流场比较光滑时,耗散很小,当雷诺数增加,流动转捩为湍流,流场梯度增大,耗散显著增大.对于MUSCL格式和二阶有界中心格式,在高雷诺数下,亚格子尺度模型没有明显改善计算结果,但也没有使计算结果恶化.中心格式相比其它两种格式,数值耗散最小,但是在高雷诺数湍流情况下,中心格式的数值耗散仍然主导了能量的耗散,再添加亚格子模型,计算结果反而变得稍差.对于工程中的低阶格式而言,采用中心格式计算大涡模拟是比较好的选择,而且在计算不存在稳定性问题时,采用不添加亚格子模型的隐式大涡模拟效果更好.     

大雷诺数Navier-Stokes方程的两水平亚格子模型稳定化方法

基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性.     

气膜冷却圆孔射流输运特性的大涡模拟研究

气膜冷却是燃机透平的关键冷却技术,其中射流问题又是气膜冷却的核心物理问题,本文采用大涡模拟方法研究了圆孔冷却射流的输运特性.文中首先采用流向准周期性方法实现了稳定的壁面湍流生成功能.通过对槽道湍流的大涡模拟验证,指出了在两种对流项格式下解析壁面湍流所需要的网格尺度,同时针对六种亚格子模型进行了验证.文中还发展了基于曲面...     

理性亚格子建模的实践与反思

亚格子建模是大涡模拟中的核心问题之一.理性亚格子建模的初衷是为尽量去除建模过程中的经验和唯象的成分,让建模过程中的每一步都在物理和数学上有迹可循.根据理性亚格子建模的理论已经发展出一系列不依赖经验常数和唯象模型的亚格子模型,但在实际问题中的效果却仍有不足,从而引发了近年来对亚格子建模思想本身的反思.本文回顾已有的研究和新的进展,目标是对未来的大涡模拟给出理性且自洽的理论指导.     

不同亚格子模式在后台阶湍流流动大涡模拟中的应用

本文用大涡模拟方法研究了湍流后台阶流场中的大涡演变过程,并在此基础上研究了目前大涡模拟中比较常用的六种亚格子模式。在相同的流动几何参数以及计算条件下,给出了不同亚格子模式下湍流流动瞬时压力场以及流场中瞬时粘性大小的分布,尤其给出了流场瞬时演变的大尺度涡结构。将六种亚格子模式从回流区长度、计算时间、计算结果的准确性和流场稳定性等不同角度进行了比较。为选取合适的亚格子模式深入研究湍流耗散机理奠定了基础。     

等离子弧焊接熔池演变过程的格子Boltzmann模拟与验证

本文应用基于分子动理论的格子Boltzmann方法,建立了描述定点等离子弧焊接熔池动态演变过程的二维数理模型,对相变过程的传热与流动现象开展模拟;根据焊接过程能量分布特点改进等离子弧的组合式热源模型,采用total-enthalpy模型求解温度、速度分布及追踪相界面。研究结果表明,模拟的熔合线形状与实验焊缝吻合,格子Boltzmann模拟得到的计算精度及计算效率均优于基于连续流体假设的有限容积法,验证了格子Boltzmann方法用于等离子弧焊接模拟的可行性和优越性;熔池中出现两个方向相反的环流,流动对焊缝形状的作用不容忽略;熔池的流动方式影响了温度场、速度场及二者协同度,直接影响固相线上的热量传递,促进了焊缝中部凸起的形成。     

Numerical simulation of instantaneous flow structure of swirling and non-swirling coaxial-jet particle-laden turbulence flows

A subgrid scale two-phase second-order-moment (SGS-SOM) model based on the two-fluid continuum approach is presented for the analysis of the instantaneous flow structures of swirling and non-swirling coaxial-jet particle-laden turbulence flows. Since the interaction between the two-phase subgrid scale stresses and the anisotropy of two-phase subgrid scale stresses is fully considered, it is superior to the conventional subgrid scale model on the basis of single gas phase or together with their similar forms for the particle phase for not taken these characters thoroughly into account. The swirling numbers s=0.47 and s=0 of coaxial-jet particle-laden turbulence flows (measured by M. Sommerfeld, H.H. Qiu, Detailed measurements in a swirling particulate two-phase flow by a phase Doppler anemometer, Int. J. Heat Fluid Flow 12 (1991) 20-28) are numerically simulated by large eddy simulation using this model, together with a Reynolds-averaged Navier-Stokes model using the unified second-order-moment two-phase turbulence model (RANS-USM). The instantaneous results show that the multiple recirculating gas flow structure is similar to that of single-phase swirling flows; but the particle flow structure contains less vortices. Both SGS-SOM and RANS-USM predicted that the two-phase time-averaged velocities and the root-mean-square fluctuation velocities are validated and are in good agreement with the experimental results. It is seen that for the two-phase time-averaged velocities both the models give almost the same results, hence the RANS-USM modeling is validated by large eddy simulation. For the two-phase root-mean-square fluctuation velocities the SGS-SOM results are obviously better than the RANS-USM results.     

可压缩湍流研究进展

高马赫数可压缩湍流的运动是一个多尺度多过程的物理现象。采用了多过程分解的方法,将可压缩湍流分解为剪切和胀压过程,分析其统计行为和动力学行为。发展了一种新的紧致差分和WENO格式相结合的混合型数值格式,准确模拟了可压缩湍流场;研究了其多尺度多过程行为和对粒子的输运影响;研究了激波结构对湍流场的影响;在高雷诺数可压缩湍流中,证明存在惯性区,其中流运动和压力做功引起的动能流通量都是常数;证明可压缩湍流中存在从大尺度到小尺度的动能级串过程;证明动能流通量的剪切部分和胀压部分在惯性区都为常数;分析亚格子应力项和亚格子质量流动项对动能级串的影响。     

Level set simulations of turbulent thermonuclear deflagration in degenerate carbon and oxygen

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localized determination of subgrid scale closure parameters is essential. We formulate a semi-localized model based on the dynamical equation for the subgrid scale turbulence energy k _sgs. The turbulent flame speed s _t is of the order √2k _sgs. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, s _t∝√2k _sgs, are investigated.     

格子Boltzmann亚格子模型的研究

为了将格子Boltzmann法应用于大雷诺数流动的模拟,本文将Smagorinsky亚格子模型和LBGK模型相结合,并对该亚格子LBM模型进行了研究。利用该亚格子LBM模型,对二维顶盖驱动流进行了模拟,得到了若干大雷诺数下流线图和方腔中心线上无量纲速度分布。计算结果与基准解进行比较,两者相互吻合。     

Investigation of anomalous very fast decay regimes in homogeneous isotropic turbulence

The emergence of anomalous fast decay regimes in homogeneous isotropic turbulence (HIT) decay is investigated via both theoretical analysis and eddy-damped quasi-normal Markovian simulations. The work provides new insight about a fundamental issue playing a role in HIT decay, namely the influence of non-standard shapes of the energy spectrum, in particular in the large energetic scale region. A detailed analysis of the kinetic energy spectrum E(k) and the non-linear energy transfer T(k) shows that anomalous decay regimes are associated with the relaxation of initial energy spectra which exhibit a bump at energetic scales. This feature induces an increase in the energy cascade rate, toward solutions with a smooth shape at the spectrum peak. Present results match observations reported in wind-tunnel experiments dealing with turbulence decay in the wake of grids and bluff bodies, including scaling laws for the dissipation parameter C_?. They also indicate that the ratio between the initial eddy turnover time and the advection time determines of how fast anomalous regimes relax toward classical turbulence free-decay. This parameter should be used for consistent data comparison and it opens perspectives for the control of multiscale effects in industrial applications.     

Particle subgrid scale modelling in large-eddy simulations of particle-laden turbulence

This study is concerned with particle subgrid scale (SGS) modelling in large-eddy simulations (LESs) of particle-laden turbulence. Although many particle-laden LES studies have neglected the effect of the SGS on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM) and the stochastic models of Fukagata et al. (Dynamics of Brownian particles in a turbulent channel flow, Heat Mass Transf. 40 (2004), 715–726) Shotorban and Mashayek (A stochastic model for particle motion in large-eddy simulation, J. Turbul. 7 (2006), 1–13) and Berrouk et al. (Stochastic modelling of inertial particle dispersion by subgrid motion for LES of high Reynolds number pipe flow, J. Turbul. 8 (2007), pp. 1–20) are analysed. The particle SGS models are assessed using both a priori and a posteriori simulations of inertial particles in a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re_λ = 74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient SGS energy, but show a large range of results dependent on the Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers, but are unable to predict preferential concentration.     

A triple level finite element method for large eddy simulations

Since large eddy simulation (LES) was introduced by Smagorinsky in 1963, it has been improved with various thoughts from many researchers. Unfortunately, despite that, the filtered approach that is widely used at present still suffers because quite empirical factors are used to determine non-closured subgrid Reynolds stresses. Based on a new definition of LES and multiscale finite element concepts, this work presents an attempt to remove such factors. Using direct sum decomposition of the solution space, we devised a hierarchical multi-level formulation of the Navier–Stokes equations for turbulence. The base-level, bearing the information of large eddies, is calculated by the conventional finite element method. The finer levels are for small scale eddies of turbulence. We address the solution methods for the small scale movements. In particular, a spectral element method is introduced for the finer level solutions. Thus large eddies and small eddies to some extent may be accurately obtained. The introduced approach offers not only access to calculate turbulence in complex geometries because of the nature of finite element method but also an effective tool for multiscale physical problems with turbulence, such as reaction flows. It is worth noting that the approach introduced here is similar to the implicit LES in finite volume and finite difference methods.     

New equation for the polarization fourier components of the small-scale velocity of an incompressible fluid in anisotropic turbulence

A new equation for the small-scale polarization Fourier components of the incompressible fluid velocity in the case of anisotropic turbulence is suggested. The principal invariant of the strain rate tensor for the large-scale velocity is found. This invariant is of most significance for the subgrid simulation of fully developed turbulence.     

三维均匀各向同性两相湍流的直接模拟

本文对三维气粒两相均匀各向同性湍流进行了直接数值模拟。气相控制方程组采用分布投影方法进行求解,微分方程采用六紧致阶差分格式和快速Fourier变换结合求解;计算颗粒场时,采用Lagrangian方法。由该方法得到的能谱和各统计量与由谱方法得到的对应值进行了比较,吻合十分理想,对不同Stocks数颗粒在流场内的瞬态分布也进行了初步模拟,并观察到局部富集现象,证明该方法是进行两相湍流直接数值模拟行之有效的方法。     

Compressible turbulence: the cascade and its locality

We prove that interscale transfer of kinetic energy in compressible turbulence is dominated by local interactions. In particular, our results preclude direct transfer of kinetic energy from large-scales to dissipation scales, such as into shocks, in high Reynolds number turbulence as is commonly believed. Our assumptions on the scaling of structure functions are weak and enjoy compelling empirical support. Under a stronger assumption on pressure dilatation cospectrum, we show that mean kinetic and internal energy budgets statistically decouple beyond a transitional conversion range. Our analysis establishes the existence of an ensuing inertial range over which mean subgrid scale kinetic energy flux becomes constant, independent of scale. Over this inertial range, mean kinetic energy cascades locally and in a conservative fashion despite not being an invariant.