完全匹配层吸收条件下的交错网格弹性波模拟

LES,DNS,RANS三种模拟模型计算量比较及其原因

LES,DNS,RANS模型计算量比较 摘要：湍流流动是一种非常复杂的流动，数值模拟是研究湍流的主要手段，现有的湍流数值模拟的方法有三种：直接数值模拟（Direct Numerical Simulation: DNS），Reynolds平均方法（Reynolds Average Navier-Stokes: RANS）和大涡模拟(Large Eddy Simulation: LES)。直接数值模拟目前只限于较小Re数的湍流，其结果可以用来探索湍流的一些基本物理机理。RANS方程通过对Navier-Stokes方程进行系综平均得到描述湍流平均量的方程；LES方法通过对Navier-Stokes方程进行低通滤波得到描述湍流大尺度运动的方程，RANS和LES方法的计算量远小于DNS，目前的计算能力均可实现。 关键词：湍流；直接数值模拟；大涡模拟；雷诺平均模型 1 引言 湍流是空间上不规则和时间上无秩序的一种非线性的流体运动，这种运动表现出非常复杂的流动状态，是流体力学中有名的难题，其 性。传统计算复杂性主要表现在湍流流动的随机性、有旋性、统计[]1 流体力学中描述湍流的基础是Navier-Stokes（N-S）方程，根据N-S 方程中对湍流处理尺度的不同，湍流数值模拟方法主要分为三种：直接数值模拟（DNS）、雷诺平均方法（RANS）和大涡模拟（LES）。直接数值模拟可以获得湍流场的精确信息，是研究湍流机理的有效手段，但现有的计算资源往往难以满足对高雷诺数流动模拟的需要，从而限制了它的应用范围。雷诺平均方法可以计算高雷诺数的复杂流动，但给出的是平均运动结果，不能反映流场紊动的细节信息。大涡模拟基于湍动能传输机制，直接计算大尺度涡的运动，小尺度涡运动对大尺度涡的影响则通过建立模型体现出来，既可以得到较雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息，又比直接数值模拟节省计算量，从而得到了越来越广泛的发展和应用。

大涡模拟

4.6.3大涡模拟LSE 大涡模拟LES 基本思想是：湍流运动是湍流运动是由许多大小不同尺度的涡旋组成，大尺度的涡旋对平均流动影响比较大，各种变量的湍流扩散、热量、质量、动量和能量的交换以及雷诺应力的产生都是通过大尺度涡旋来实现的，而小尺度涡旋主要对耗散起作用，通过耗散脉动来影响各种变量。不同的流场形状和边界条件对大涡旋有较大影响，使它具有明显的各向不均匀性。而小涡旋近似于各向同性，受边界条件的影响小，有较大的共性，因而建立通用的模型比较容易。据此，把湍流中大涡旋(大尺度量)和小涡旋(小尺度量)分开处理，大涡旋通过N-S 方程直接求解，小涡旋通过亚格子尺度模型，建立与大涡旋的关系对其进行模拟，而大小涡旋是通过滤波函数来区分开的。对于大涡旋，LES 方法得到的是其真实结构状态，而对小涡旋虽然采用了亚格子模型，但由于小涡旋具有各向同性的特点，在采用适当的亚格子模式的情况下，LES 结果的准确度很高。 大涡模拟LES 有四个一般的步骤： ①定义一个过滤操作，使速度分解u(x,t)为过滤后的成分(),u x t 和亚网格尺度成分u ’(x,t)，这里要特别指出：过滤操作和Reynolds 分解是两个不同的概念，亚网格尺度SGS 成分u ’(x,t)与Reynolds 分解后的速度脉动值是两个不同的量。过滤后的三维的时间相关的成分()t x u ,表示大尺度的涡旋运动； ②由N-S 方程推导过滤后的速度场进化方程，该方程为一个标准形式，其中包含SGS 应力张量； ③封闭亚网格尺度SGS 应力张量，可采用最简单的涡黏性模型； ④数值求解模化方程，从而获得大尺度流动结构物理量。 （1）过滤操作 LES 方法和一般模式理论不同之处在于对N-S 方程第一步的处理过程不一样。一般模式理论方法是对变量取平均值，LES 方法是通过滤波操作，将变量分成大尺度量和小尺度量。对任一流动变量(),u x t 划分为大尺度量(,)u x t 和小尺度量(),u x t '（亚格尺度）： (,)(,)(,)u x t u x t u x t '=+ 其中大尺度量是通过滤波获得:，过滤操作定义为： ()?-=dr t r x u x r G t x u ),(),(, （4.78） 式中积分遍及整个流动区域，(,)G r x 是空间滤波函数，它决定于小尺度运动的尺寸和结构。 滤波器G 要满足正规化条件 ?=1),(dr x r G （4.79） 亚网格尺度SGS 成分定义为 ),(),(),('t x u t x u t x u -= （4.80） 与Reynolds 分解不同的是，),(t x u 为一个随机的场分布，且 0),('≠t x u

边界条件的设置

第二章：边界条件 这一章主要介绍使用边界条件的基本知识。边界条件能够使你能够控制物体之间平面、表面或交界面处的特性。边界条件对理解麦克斯韦方程是非常重要的同时也是求解麦克斯韦方程的基础。 §2.1 为什么边界条件很重要 用Ansoft HFSS求解的波动方程是由微分形式的麦克斯韦方程推导出来的。在这些场矢量和它们的导数是都单值、有界而且沿空间连续分布的假设下，这些表达式才可以使用。在边界和场源处，场是不连续的，场的导数变得没有意义。因此，边界条件确定了跨越不连续边界处场的性质。 作为一个 Ansoft HSS 用户你必须时刻都意识到由边界条件确定场的假设。由于边界条件对场有制约作用的假设，我们可以确定对仿真哪些边界条件是合适的。对边界条件的不恰当使用将导致矛盾的结果。 当边界条件被正确使用时，边界条件能够成功地用于简化模型的复杂性。事实上，Ansoft HFSS 能够自动地使用边界条件来简化模型的复杂性。对于无源RF 器件来说，Ansoft HFSS 可以被认为是一个虚拟的原型世界。与边界为无限空间的真实世界不同，虚拟原型世界被做成有限的。为了获得这个有限空间，Ansoft HSS使用了背景或包围几何模型的外部边界条件。 模型的复杂性通常直接与求解问题所需的时间和计算机硬件资源直接联系。在任何可以提高计算机的硬件资源性能的时候，提高计算机资源的性能对计算都是有利的。 §2.2 一般边界条件 有三种类型的边界条件。第一种边界条件的头两个是多数使用者有责任确定的边界或确保它们被正确的定义。材料边界条件对用户是非常明确的。 1、激励源 波端口（外部） 集中端口（内部） 2、表面近似 对称面 理想电或磁表面 辐射表面 背景或外部表面 3、材料特性 两种介质之间的边界 具有有限电导的导体 §2.3 背景如何影响结构 背景边界:所谓背景是指几何模型周围没有被任何物体占据的空间。任何和背景有关联的物体表面将被自动地定义为理想的电边界（Perfect E）并且命名为外部（outer）边界条件。你可以把你的几何结构想象为外面有一层很薄而且是理想导体的材料。 有耗边界:如果有必要，你可以改变暴露于背景材料的表面性质，使其性质与

基于交错网格有限差分弹性波正演模拟及波场特征分析

基于交错网格有限差分弹性波正演模拟及波场特征分析 【摘要】为研究和认识多种储层中弹性波的波场特征，以利于多波地震资料解释，高精度数值模拟是有效的方法之一。本文在弹性波方程基础上，采用高阶交错网格有限差分技术模拟地震波在各向同性介质和各向异性介质中的传播，可得到不同类型介质的弹性波场。同时，文中也分析了各向异性系数对多波波场特征的影响。通过对高精度数值模拟得到的波场快照对比研究表明，该方法可有效获得高精度弹性波正演结果，为研究各种复杂介质中弹性波的波场特征和传播规律奠定了基础。 【关键词】多波多分量波场特征各向异性弹性波正演 1 引言 随着油气田勘探技术的不断发展[1][2]，人们对地震资料的认识也不断加深，纵波地震资料在含油气的显示上存在一些不确定性，单一纵波资料解释的多解性问题尤为突出。在地震勘探领域中，过去一直把各向同性弹性体理论作为研究地下介质的前提，但是在实际地层中普遍存在各向异性，地下介质的各向异性（如周期薄互层引起的各向异性、以及裂隙引起的各向异性）产生的弹性波场与各向同性介质产生的弹性波场存在着不可忽略的差异。由此，多波地震勘探作为油储地球物理的主要方法之一应运而生。在多波资料解释过程中，要求搞清楚储层的岩性与多波的波场特征之间的关系，因此，多波波场数值模拟技术显得非常重要。高精度数值模拟技术是联系地震、地质、测井以及油藏工程的纽带，其作用主要体现在提高人们对各种复杂介质中地震波传播规律的认知，并可为新技术、新方法提供试验数据，以满足方法技术研究的需要，同时也可以检验解释结果的正确性。 弹性波波动方程高精度数值模拟可以得到全波场信息，包含了地震波的动力学和运动学特点，为准确描述地震波场特征和波的传播规律奠定基础，本文在弹性波方程基础上，采用高阶交错网格有限差分技术模拟地震波在各向同性介质和各向异性介质中的传播，比较地震波在各向同性介质和各向异性介质中的波场响应异同，并分析了各向异性系数对多波波场特征的影响，这对研究各种复杂介质中弹性波的波场特征和传播规律有着重要的意义。 2 高精度正演模拟方法原理 多波波场正演模拟以弹性波方程为基础，其核心是研究高精度数值模拟技术，实现高精度正演。 2.1 各向同性介质中的弹性波方程 弹性波理论主要是依据物体受力与形变之间的关系，通过应力和应变分析、弹性体运动平衡关系分析，进而建立弹性波波动方程[3]：

第三章-数值模拟理论与方法

第三章 数值模拟理论与方法 §3.1 流体力学的基本方程 流体运动所遵循的规律是由物理学三大守恒定律规定的，即质量守恒定律，动量守恒定律和能量守恒定律[44]。 （一）连续方程 0)(=?+??v t ρρ （3.1） 式中 ρ－流体密度 u －流体速度分量 （二）动量方程（x 方向） 对于不可压流体（即0=?v ） x p f v u v x u x ??-+??=??+??ργρρρ)()()( （3.2） 式中 γ－运动粘性系数 p －压力 对于可压缩流体 ()()()()()x p f v x u u v x u x ??-+???+??=????ργργρρρ 31 （3.3） 式中等号后前两项是粘性力 y ，z 方向上的动量方程可类似推出。 （三）能量方程 ()()()v q T k e v e t ερρ++???=??+?? （3.4） 其中 T C e v = 式中等号左边第一项是瞬变项，第二项是对流项，等号右边第一项是扩散项，第二、三项是源项。 所以，流体力学基本方程组为： ()0=?+??v t ρρ

()x p f u u v f t u x ??-+??=??+??ργρ)( ()()y p f v v v f t v y ??-+??=??+??ργρ （3.5） ()()w p f w w v f t w w ??-+??=??+??ρλρ ()()v q e c k e v f e t v ερ++??? ? ????=??+?? §3.2 紊流模式理论概况 §3.2.1 基本方程 在自然界中，真实的流体都具有粘性。粘性流体存在两种不同的运动方式和流态，即层流和紊流。而在自然界和工农业生产中所遇见的流体流动大部分都是紊流。 三维的N-S 方程是目前描述粘性流体运动较为理想的模型，其优点一是应用范围广，在空气、水流、传热等方面均用N-S 方程描述；二是对于有分离、旋涡等情况的复杂三维流动更为适用。 三维直角坐标下的N-S 方程[45]，[46]，即不可压缩粘性流体的动量方程式为： ?????????????+??+??+??-=??+??+??+??-=??+??+??+??-=)()()(222222222222222222z w y w x w z p F Dt Dw z v y v x v y p F Dt Dv z u y u x u x p F Dt Du z y x μρρμρρ μρρ （3.6） 不可压缩流体的连续性方程为： （3.7） 式（3.6）和（3.7）共有四个未知数（u 、v 、w 、p ）和四个方程，加上边界条件，从理论上来讲其解是存在的。但是，要直接求解复杂而详细的粘性流体运动是十分复杂和困难的。其原因是：直接求解N-S 方程要求求解从反映消散运动的最小涡漩尺度到反映大尺度涡体的所有流动尺度，因而只有对简单情况下才有理论解。 0=??+??+??z w y v x u

各类边界条件fluent

Fluent技巧 边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域 2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变 (如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤： 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者，1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.，选择边界条件区域将会打开，并且你可以指定适当的边界条件

大涡模拟的FLUENT算例2D

Tutorial:Modeling Aeroacoustics for a Helmholtz Resonator Using the Direct Method(CAA) Introduction The purpose of this tutorial is to provide guidelines and recommendations for the basic setup and solution procedure for a typical aeroacoustic application using computational aeroacoustic(CAA)method. In this tutorial you will learn how to: ?Model a Helmholtz resonator. ?Use the transient k-epsilon model and the large eddy simulation(LES)model for aeroacoustic application. ?Set up,run,and perform postprocessing in FLUENT. Prerequisites This tutorial assumes that you are familiar with the user interface,basic setup and solution procedures in FLUENT.This tutorial does not cover mechanics of using acoustics model,but focuses on setting up the problem for Helmholtz-Resonator and solving it.It also assumes that you have basic understanding of aeroacoustic physics. If you have not used FLUENT before,it would be helpful to?rst review FLUENT6.3User’s Guide and FLUENT6.3Tutorial Guide. Problem Description A Helmholtz resonator consists of a cavity in a rigid structure that communicates through a narrow neck or slit to the outside air.The frequency of resonance is determined by the mass of air in the neck resonating in conjunction with the compliance of the air in the cavity. The physics behind the Helmholtz resonator is similar to wind noise applications like sun roof bu?eting. We assume that out of the two cavities that are present,smaller one is the resonator.The motion of the?uid takes place because of the inlet velocity of27.78m/s(100km/h).The ?ow separates into a highly unsteady motion from the opening to the small cavity.This unsteady motion leads to a pressure?uctuations.Two monitor points(Point-1and Point-2) act as microphone points to record the generated sound.The acoustic signal is calculated within FLUENT.The?ow exits the domain through the pressure outlet.

FLUENT中各种边界条件的适用范围

FLUENT中各种边界条件的适用范围 速度入口边界条件：用于定义流动入口边界的速度和标量。 压力入口边界条件：用来定义流动入口边界的总压和其它标量。 质量流动入口边界条件：用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流，因为当密度是常数时，速度入口边界条件就确定了质量流条件。压力出口边界条件：用于定义流动出口的静压（在回流中还包括其它的标量）。当出现回流时，使用压力出口边界条件来代替质量出口条件常常有更好的收敛速度。 压力远场边界条件：用于模拟无穷远处的自由可压缩流动，该流动的自由流马赫数以及静态条件已知。这一边界类型只用于可压缩流。 质量出口边界条件：用于在解决流动问题之前，所模拟的流动出口的流速和压力的详细情况还未知的情况。在流动出口是完全发展的时候这一条件是适合的，这是因为质量出口边界条件假定出了压力之外的所有流动变量正法向梯度为零。不适合于可压缩流动。 进风口边界条件：用于模拟具有指定的损失系数、流动方向以及周围（入口）环境总压和总温的进风口。 进气扇边界条件：用于模拟外部进气扇，它具有指定的压力跳跃、流动方向以及周围（进口）总压和总温。 通风口边界条件：用于模拟通风口，它具有指定的损失系数以及周围环境（排放处）的静压和静温。 排气扇边界条件：用于模拟外部排气扇，它具有指定的压力跳跃以及周围环境（排放处）的静压。 速度入口边界条件：速度入口边界条件用于定义流动速度以及流动入口的流动属性相关标量。这一边界条件适用于不可压缩流，如果用于可压缩流它会导致非物理结果，这是因为它允许驻点条件浮动。应该注意不要让速度入口靠近固体妨碍物，因为这会导致流动入口驻点属性具有太高的非一致性。 压力入口边界条件：压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压缩流，也可以用于不可压缩流。压力入口边界条件可用于压力已知但是流动速度和/或速率未知的情况。这一情况可用于很多实际问题，比如浮力驱动的流动。压力入口边界条件也可用来定义外部或无约束流的自由边界。 质量流动入口边界条件：用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流，因为当密度是常数时，速度入口边界条件就确定了质量流条件。当要求达到的是质量和能量流速而不是流入的总压时，通常就会使用质量入口边界条件。调节入口总压可能会导致解的收敛速度较慢，所以如果压力入口边界条件和质量入口条件都可以接受，应该选择压力入口边界条件。 压力出口边界条件：压力出口边界条件需要在出口边界处指定静（gauge）压。静压值的指定只用于亚声速流动。如果当地流动变为超声速，就不再使用指定压力了，此时压力要从内部流动中推断。所有其它的流

大涡模拟的fluent算例

Introduction:This tutorial demonstrates how to model the2D turbu-lent?ow across a circular cylinder using LES(Large Eddy Simula-tion),and compute?ow-induced noise(aero-noise)using FLUENT’s acoustics model. In this tutorial you will learn how to: ?Perform2D Large Eddy Simulation(LES) ?Set parameters for an aero-noise calculation ?Save surface pressure data for an aero-noise calculation ?Calculate aero-noise quantities ?Postprocess an aero-noise solution Prerequisites:This tutorial assumes that you are familiar with the menu structure in FLUENT,and that you have solved or read Tu-torial1.Some steps in the setup and solution procedure will not be shown explicitly. Problem Description:The problem considers turbulent air?ow over a2D circular cylinder at a free stream velocity U of69.19m/s. The cylinder diameter D is1.9cm.The Reynolds number based on the?ow parameters is about90000.The computational do-main(Figure3.0.1)extends5D upstream and20D downstream of the cylinder,and5D on both sides of it.If the computational domain is not taken wide enough on the downstream side,so that no reversed?ow occurs,the accuracy of the aero-noise prediction may be a?ected.The rule of thumb is to take at least20D on the downstream side of the obstacle. c Fluent Inc.June20,20023-1

fluent边界条件(一)

边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分，设定边界条件必须小心谨慎。 边界条件的分类：进出口边界条件：压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇；壁面、repeating, and pole boundaries:壁面，对称，周期，轴；内部单元区域：流体、固体(多孔是一种流动区域类型) ；内部表面边界：风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面，这意味着它们没有有限厚度，并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件，并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍，模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型，并且打开其他的面板以设定每一区域的边界条件参数 菜单：Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前，必须检查所有边界区域的区域类型，如有必要就作适当的修改。比方说：如果你的网格是压力入口，但是你想要使用速度入口，你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下：: 1.在区域下拉列表中选定所要修改的区域

2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时，点击确认 确认改变之后，区域类型将会改变，名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 ！注意：这个方法不能用于改变周期性类型，因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是，只能在图一中每一个类别中改变边界类型(注意：双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中，边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件，请遵循下面的步骤： 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者，1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.，选择边界条件区域将会打开，并且你可以指定适当的边界条件 在图像显示方面选择边界区域 在边界条件中不论你合适需要选择区域，你都能用鼠标在图形窗口选择适当的区域。如果你是第一次设定问题这一功能尤其有用，如果你有两个或者更多的具有相同类型的区域而且你想要确定区域的标号（也就是画出哪一区域是哪个）这一功能也很有用。要使用该功能请按下述步骤做： 1.用网格显示面板显示网格。 2.用鼠标指针（默认是鼠标右键——参阅控制鼠标键函数以改变鼠标键的功能）在图形窗口中点击边界区域。在图形显示中选择的区域将会自动被选入在边界条件面板中的区域列表中，它的名字和编号也会自动在控制窗口中显示改变边界条件名字 每一边界的名字是它的类型加标号数（比如pressure-inlet-7）。在某些情况下你可能想要对边界区域分配更多的描述名。如果你有两个压力入口区域，比方说，你可能想重名名它们

燃气轮机模型燃烧室的大涡模拟

燃气轮机模型燃烧室的大涡模拟 徐宝鹏1，曾佑杰1，马宏宇2，赵凯岚2，金戈2 （1.大连理工大学能源与动力学院，辽宁大连116024；2.中航工业沈阳发动机设计研究所，沈阳110015） 摘要：燃烧室内的燃油雾化、蒸发以及和空气进行混合过程对燃烧过程有重要影响。提出1种基于大涡模拟的数学模型来模拟燃烧室内燃料喷射、蒸发和混合过程。被空间滤波掉的亚网格尺度涡对大尺度涡的影响由求单方程SGS 湍流模型进行模拟。采用拉格朗日法和蒙特卡洛技术对流场中的喷雾粒子进行采样跟踪，采样喷雾粒子在流场中作为点源项与气相进行质量、动量和能量的双向耦合。提出1个基于SGS 湍流动能的双向耦合模型来模拟SGS 脉动速度对喷雾粒子运动的影响以及喷雾相对SGS 湍流动能的影响。通过对1个同轴模型燃烧室中的喷雾蒸发及混合过程的大涡模拟，将预测结果和试验值进行了比较，预测值和试验值吻合良好，验证了模型的可靠性。 关键词：燃烧室；燃气轮机；大涡模拟；双向耦合；燃油雾化 中图分类号：V211.3文献标识码：A doi ：10.13477/https://www.wendangku.net/doc/fe2260114.html,ki.aeroengine.2014.03.003 Large Eddy Simulation of a Gas Turbine Model Combustor XU Bao-peng 1,ZENG You-jie 1,MA Hong-yu 2,ZHAO Kai-lan 2,JIN Ge 2 (1.School of Energy and Power Engineering,Dalian University of Technology,Liaoning Dalian 116024,China;2.AVIC Shenyang Engine Design and Research Institute,Shenyang 110015,China ) Abstract:Fuel atomization,evaporation and mixing with air in gas turbine combustors are vital to the subsequent combustion process.Numerical formulation based on large eddy simulation is proposed to model fuel injection,evaporation and mixing in a gas turbine combustor.The proposed model adopts a one-equation subgrid scale turbulent model to handle the effect of the filtered subgrid scale eddies on the solved large scale eddies.Spray droplets are tracked using both Lagrangian method and Mento Carlo technique,and the sampled spray particles are regarded as point sources to conduct two-way couplings of mass,momentum and energy.The two-way coupling model based on SGS turbulent kinetic energy is used to model the mutual influences between SGS fluctuating velocity and the movement of spray droplets.The proposed models are validated against a large simulation of a co-axial model combustor and the predictions are compared to the experimental data.Good agreements are obtained,which demonstrate the reliability of the proposed models.Key words:combustor ;gas turbine;large eddy simulation;two-way coupling;atomization 航空发动机Aeroengine 第40卷第3期 Vol.40No.3 Jun.2014 收稿日期：2013-08-14基金项目：燃气轮机重大项目联合培育基金（2011LH006）资助 作者简介：徐宝鹏（1969）,男，博士，教授，主要研究方向为计算流体力学、两相流和燃烧学；E-mail:xbp624@https://www.wendangku.net/doc/fe2260114.html, 。引用格式：0引言 试验研究表明，燃油在燃气轮机燃烧室内的雾 化、蒸发及与空气进行混合的过程对直喷或预混燃烧 室中的燃烧过程起至关重要的作用[1-3]。Fric [3]通过试 验研究发现，燃料混合在空间上的不均匀性和时间上 的脉动性对氮氧化物的排放量有显著影响；此外,燃 烧室内的流动状态以及燃油喷雾粒子的大小和速度 分布对燃烧过程同样有重要影响。光学测试和数值模拟是目前研究燃烧室内燃油混合及燃烧过程的2种主要方法。光学测试已被成功应用于燃油喷射及混合过程的研究中[4-6]，但使用数值模拟方法对燃烧前喷雾的分布情况进行理论研究的文献较少，且现有的数值模拟工作大都采用基于雷诺时间平均的RANS 方法。以往的数值研究表明，RANS 方法无法准确预测燃烧室内的回流流动，且其稳态特性也不能准确预测燃料的混合及燃烧过程[7]。大涡模

声波方程有限差分正演

题目：使用Ricker 子波，刚性边界条件，并且初值为零，在均匀各向同性介质条件下，利用交错网格法求解一阶二维声波方程数值解。 解： 一阶二维声波方程： 22222221z P x P t P c ??+??=?? (1) 将其分解为： 21P c t P x P z x z x z V V x z V t V t ????=+????????=???????=???? (2) 对分解后的声波方程进行离散，可得到： 1 12211,-1,,,122[]N n n n n m i m j i m j xi j xi j m t V V c P P h + -+---=?=+-∑ 1 1 221 1,1,,,122 []N n n n n m i j m i j m zi j zi j m t V V c P P h +-++---=?=+-∑ 111121 2222,,m 1,,,,11 []N n n n n n n i j i j m xi j xi m j zi j m zi j m m tc P P c V V V V h +++++++-+--=?=+-+-∑ h z x =?=? 针对公式(1)，使用二阶中心差商公式： 2P(,,1)2(,,)(,,1)i j n P i j n P i j n t +-+-?222(1,,)2(,,)(1,,)(,1,)2(,,)(,1,)P i j n P i j n P i j n x c P i j n P i j n P i j n z +-+-??+?????=??+-+-??????? (3) 变形： P(,,1)=2(,,)(,,1)i j n P i j n P i j n +--

大涡模拟简单介绍

《粘性流体力学》小论文 题目：浅谈大涡模拟 学生姓名：丁普贤 学生学号：103911018 完成时间：2010/12/16

浅谈大涡模拟 丁普贤 （中南大学,能源科学与工程学院，湖南省长沙市，410083） 摘要：湍流流动是一种非常复杂的流动，数值模拟是研究湍流的主要手段，现有的湍流数值模拟的方法有三种：直接数值模拟、大涡模拟和雷诺平均模型。本文主要是介绍大涡模拟，大涡模拟的思路是：直接数值模拟大尺度紊流运动，而利用亚格子模型模拟小尺度紊流运动对大尺度紊流运动的影响。大涡模拟在计算时间和计算费用方面是优于直接数值模拟的，在信息完整性方面优于雷诺平均模型。本文还介绍了对N-S方程过滤的过滤函数和一些广泛使用的亚格子模型，最后简单对一些大涡模拟的应用进行了阐述。 关键词：计算流体力学；湍流；大涡模拟；亚格子模型

A simple study of Large Eddy Simulation DING Puxian (Central South University, School of Energy Science and Power Engineering, Changsha, Hunan, 410083) Abstract：Turbulent flow is a very complex flow, and numerical simulation is the main means to study it. There are three numerical simulation methods: direct numerical simulation, large eddy simulation,Reynolds averaged Navier-Stokes method. Large eddy simulation (LES) is mainly introduced in this paper. The main idea of LES is that large eddies are resolved directly and the effect of the small eddies on the large eddies is modeled by subgrid scale model. Large eddy simulation calculation in computing time and cost is superior to direct numerical simulation, and obtain more information than Reynolds averaged Navier-Stokes method. The Navier-Stokes equations filtering filter function and some extensive use of the subgrid scale model are simply discussed in this paper. Finally, some simple applications of large eddy simulation are told. Key words：computational fluid dynamics; turbulence; large eddy simulation; subgrid scale model

自由边界条件下克拉尼图形的模拟

课程设计 题目：自由边界条件下克拉尼图形的模拟 学院：物理工程学院 姓名：文美乐 专业：应用物理1203班 学号：121410020322

t=0:pi/20:4*pi; %0到4pi，每隔Pi/20取一个数 x=0:0.1*pi:4*pi; %0到4pi，每次增加0.1pi m=moviein(40); %定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面for i=1:40; y=sin(x+t(i)); plot(x,y ,'o'); m(:,i)=getframe; %依次将40帧画面存入矩阵m end movie(m,10,20) %将40帧画面演示10次，每秒20帧画面 图形：

t=0:pi/20:4*pi; %0到4pi，每隔Pi/20取一个数 x=0:0.1*pi:4*pi; %0到4pi，每次增加0.1pi m=moviein(40); %定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面for i=1:40 ; y1=sin(x+t(i)); y2=sin(x-t(i)); y=y1+y2;%两列行波叠加 plot(x,y1,x,y2,x,y); m(:,i)=getframe; %依次将40帧画面存入矩阵m end movie(m,10,20) %将40帧画面演示10次，每秒20帧画面 图形：

t=[0:pi/20:4*pi];%0到4*pi，每隔pi/20取一个数 x=[0:0.1*pi:4*pi];%0到4*pi,每次增加0.1*pi y=x; [X,Y]=meshgrid(x,y);%定义数组x,y， m=moviein(40); %定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面for i=1:40 ; Z=sin(X+t(i))+sin(X-t(i));%定义z函数 mesh(X,Y,Z);%画三维图 m(:,i)=getframe; %依次将40帧画面存入矩阵m end movie(m,10,20) %将40帧画面演示10次，每秒20帧画面 图形：

Midas各种边界条件比较

Midas各种边界条件比较 Midas的提供的边界条件非常多，而且各有用途，初学Midas的朋友们都想看看到底不同边界条件之间有什么区别，下面在Midas帮助文件选取下来的，只是作一个比较，各种边界条件的具体使用参照MIDAS帮助文件。 1.定义一般弹性支承类型 SDx-SDy 整体坐标系X轴方向和Y轴方向(或已定义的节点局部坐标系x方向和y方向)的相关弹性支承刚度。 注 一般弹性支承通常用于反映桩的支承刚度，结构分析时可以考虑与各个自由度有关的桩支承刚度。 在典型的建筑结构中，分析模型不包括桩基础。而是假定在基础底面或桩帽处存在弹性边界。 下面的通用刚度给出了桩单元的实际刚度。对斜桩，用节点局部坐标轴计算斜向的刚度。 2.一般弹性支承 分配定义的一般弹性支撑类型，或输入节点通用刚度矩阵(6×6)。其中包括选定的节点在整体坐标系或节点局部坐标系内各自由度之间相关的刚度，也可以替换或删除先前定义的弹性支承刚度 SDxSDySDzSRxSRySRz 注 在一般弹性支承类型对话框中，上述6个弹性支承刚度值只表示6 x 6阶刚度矩阵中的6个对角线刚度值。实际分配给节点的刚度值为6 x 6阶刚度。 3.面弹性支承

输入平面或实体单元单位支承面上的弹簧刚度形成弹性支承。并可同时形成弹性连接的单元。 该功能主要用于在基础或地下结构分析中考虑地基的弹性支承条件。 弹性连接长度：弹性连接单元的长度。该数据对分析结果没有影响，只是为在分析中定义一个内部矢量。 只受拉，只受压：选中选项指定弹性连接为只受拉或只受压单元。 4.弹性连接 形成或删除弹性连接。由用户定义弹性连接及其弹性连接的两个节点。 SDxSDySDzSRxSRySRz。 5.一般连接特性值 建立、修改或删除非线性连接的特性值。一般连接功能应用于建立减隔振装置、只受拉/受压单元、塑性铰、弹性支撑等模型。一般连接可利用弹簧的特性，赋予线性或非线性的特性。 一般连接的作用类型分为单元类型和内力类型。 单元类型一般连接在进行分析过程中，用更新单元刚度矩阵直接反映单元的非线性。 内力类型的一般连接不更新单元刚度矩阵，而是根据非线性的特性计算出来的内力置换成外部荷载，间接的考虑非线性。单元类型的一般连接提供的类型有弹簧、线性阻尼器、弹簧和线性阻尼器3种类型的连接单元。 内力类型的一般连接提供的类型有粘弹性消能器(Viscoelastic Damper)、间隙(Gap)、钩(Hook)、滞后系统(Hysteretic System)、铅芯橡胶支承隔震装置(Lead Rubber