大气湍流模型_.pdf
激光(Laser：Light Amplification by Stimulated Emission of Radiation)是光受激辐射放大的简称。原子受光子照射时不仅发生受激辐射同时还发生受激吸收。这两种过程是矛盾的。通常情况下吸收过程是主要的，受激辐射过程是次要的. 如果能够通过某种方法破坏粒子数的热平衡分布，受激辐射能量将大于吸收能量，受激过程将胜于吸收过程。这里E[v2]是单位为(m/s)2的速度平方平均值,离开地面高度h的单位为km,h的范围5~20km。近期实验证实了这个模型的近似合理性,并指出在儿百米以卜的近地面范围内,在白天Cn(h)≈C2()h43,在仅晚C2(h)≈C2()h23,另外Hal!发现8在草地覆盖面上空,Cn(h)≈C2()h3, Davidson发现海面上空C(h)≈C2()h23,1m≤h≤10m。一般m言,在近地面处C2的典型值从102m23(对于强湍流)到108m23(对于弱湍流)的范围内折射率功率谱密度折射率的随机起伏n1(r)主要是由温度空间分布中的随机微观结构而引起的,这种微观结构的起源则在于地球表面不同区域被太阳不同加热而引起的极大尺度的温度非均匀性,这种大尺度的温度非均匀性进而又引起大尺度的折射率非均匀性,它们最后被湍流风和对流冲碎,使非均匀性的尺度变得越来越小。湍流的大小范围通常从几毫米到几米,分别用内长度l和外长度L表示通常把大气折射率的非均匀性称为湍流“旋涡”,可以把他们想像成一些空气包,每个空气包鄞有一个特征的折射率。均匀湍流的功率谱密度Φ,(K)可以看成是尺度为L=2x/k,L=2x/和L,=2x/k的旋涡的相对丰度的一种量度。在各向同性湍流的情况下,Φ(K)仪是波数k的函数,k通过L=2x/k与旋涡大小L相联系。在 Kolmogorov关于湍流理论的终典工作的基础,普遍认为功率谱密度Φn(K)包括三个不同的区。对于很小的k=2n/L0(很大规模的尺寸)的区域叫输入区,在这个区域内谱的形状取决于特定的湍流是如何发生的,而且它通常是各向异性的。在这个区域被理论不能预言Φ(K)的数学形式。但当k大于某一临界波数λ时,Φ(K)的形状由制约着大湍流旋涡破碎为小旋涡的物理定律来决定。当k大于k0时,k0≈2x/,就进入了谱的惯性子区间。这里的Φn的形式可以由已确立的制约湍流的物理定律描述。由Kσ oImogoro湍流理论,Φ,为D(K)=0033C2k-1/32.5)当k达到了另一个临界值k的形式再次改变,这个区域叫耗散区,在这个区域里能量的耗散超过了动能。因此能量很小。所以,当k>k时,Φ很快下降。这里kn≈2π/l。 Tatarski川如下模型来概括k>k时Φn的快速下降:29C1994-2010ChinaAcademicJournalElcctronicPublishinghOusc.Allrightsrescrved.http://www.cnki.ncton(k)=0.033C, exp(k/km)(2.6)若选取kn=5.92/l,并且k>kn,上式是一个合理的近似。由式(25)和(2.6)所表示的谱在原点均偶不可积的极点,为了克服这种模型的缺点,常采用一种称为Ⅴ on karman谱的形式。这时谱近似地表示为0.033CΦ(k)≈(k2+k2)16 exp(- /km)Andrews提出了一个的近似谱8Φ(0.033 CR exp(-k2/k2)7/61+aQ2((k2+k2)16(28)k1k1其屮,a1=1.802,a2=0.254,k,-3.3/。注意,在a1=a2=0和作k=kn代换后Andrews模式简化为 Von karman谱;当k=l=0时,上式退化为(25)式。激光在大气中的传输方程假定大气的磁导率为常数,介电常量是随空间变化的。对于单色电磁波在地球大气中的传播,这时 Maxwell方程取以下形式l9V·H=0V×E=jboH(2.9)EOoEV·(EE)=0式(2.9)中,E是电场;H是磁场;a是角频率;而V矢量的分量为(O/ax,0/Cy,O/az)。把V×运算应用于式(29)的第二个方程,并将第二方程代入第三个方程,并考虑第四个方程,我们得到V·E+AHEE+V(E·Vlns)=0(2.10)这里代表以e为底的对数。波传播的局域速度即某点上的速度是(e)12,它也等于ch,式中,c是自由空间中的光速,而n是同一点上的局部折射率,因此u8=n/(2.11)由于和c是常数,有VIne=2vIn n(2.12)将式(2.11)和(2.12)代入式(210),得到C1994-2010ChinaAcademicJournalElcctronicPublishinghOusc.Allrightsrescrved.http://www.cnki.nct
- 2020-12-04下载
- 积分：1
fluent udf中文帮助
fluent udf中文帮助本章简要地介绍了用户自定义函数(UDF)及其在 Fluent 中的用法。在 1.1 到 1.6 节中我们会介绍一下什么是 UDF；如何使用 UDF，以及为什么要使用 UDF，Makefile udfmakefile udf2DEF INE PROPERTYDEFINE DRIFT DI AM基础输运方程单元()面,区域()和线(操作求解器数据运仃输运方程求解器建立在有限容积法的基础上,这种方法将计算域离散为有限数目的控制体或是单元。网格单元是中基本的计算单元,这些单元的守恒特性必须保证。也就是说普通输运方程,例如质量,动量,能量方程的积分形式可以应用到锊个单元:/pad+,p·dA=中 TODa+/,SwAunsteadyeonv画 ctiondiffusiongeneration此处,是描述普通输运数量的变量()根据所求解的输运方程它可取不同的值。下面是在输运方程中可求解的的子集E守恒与否需要知道通过单元边界的通量。因此,需计算出单元和面上的属性值(单元()面,区域()和线(单元和单元面被组合为一些区域(),这些区域规定了计算域(例如,入口,出口,壁面)的物理组成)。当用户使用的时,用户的叮调用流体区域或是边界区域的计算变量(需要获得适当的变量,比如说是区域参考)和单元,以便标定各个单元。区域)是一群单元或单元面的集合,它可以由模型和区域的物理特征(比如入∏,出∏,壁面,流伓区域)来标定。例如,·些被指定为面域()的单元面可以被指定为类型,由此,速度也就可指定了。线()是据结构的内部名称,可被用来指定一个区域结构可作为数据储存器米使用,这些数据对于它所表示的单元和面来说是公用的(操作多数的任务需要在一个线的所有单元和面上重复执行。比如,定义一个自定义轮廓函数()则会对一个面线上)的所有单元和面进行循环。为了用户方使向用户提供了一些循环宏工具(来执行对单元,面,节点()和线()的重复操作例如,单元循环宏()可以对给定单元线上的所有单元进行循环操作()。而面循环宏()则可调用所有给定面线()的面。提供的循环工具请见在某些情况下需要对某个变量操作,而这个变量恰恰又不能直接被当作变量来传递调用。比如,如果用户使用宏来定义,求解器将不会向它传递指针。这种情况下,用户函数需要用提供的宏来调用线指针(见求解器数据通过用户界面将函数(它已被编译和连接)连接到求解器上可实现调用求解器变量。旦和求解器正确连接,无论何时,函数都可调用求解器数据。这些数据将会被作为用户变量自动地传递给注意,所有的求解器变量,不管是求解器传递给的,还是传递给求解器的,都使用单位。运行将会在侦定时刻被调用。但是,也可对它们进行异步执行,使用宏,还可在需要时()执行。详情请见解释和编译的比较编译和的构建方式一样。脚木被用来调用编译器来构建一个当地目标代码库()。目标代码库包含高级语言源代码的机器语言翻译。代码库在运行时由“动态加载”()过程连接到上。连接后,与共亨库的联系()将会被保存在用户的文件中这样,当以后再读入文件时,此编译库将会与自动连接。这些库是针对计算机的体系结构和·定版本的使用的。所以,当更新,或计算机操作系统改变,或是在不同类型的机器上运行时,这些库必须重新构建而解释则是在运行吋,直接从语言源代码编译和装载()。在运行中,源代码被编译为中介的独立于物理结构的使用预处理程序的机器代码(被调用时,机器代码由內部仿真器(),或注释器)执行。注释器不具备标准编译器的所有功能;它不支持语言的某些原理。所以,在使用吋,有语言限制(见)。例如,不能够通过废弃结构()来获得数据。要获得数据结构,必须使用由提供的预定义宏。另一个例子是不能识别指针数组。这些功能必须由来执行。编译后,用户的函数名称和内容将会被储存在文件中。函数将会在读入文件时被自动编译。独立于物理结构的代码的外层()叮能会导致执行错误(),但却可使共亨不同的物理结构,操作系统,和版本。如果运行速度较慢,不用被调节就可以编译代码的形式(运行中的和请见选择或是时,注意以下内容:对其它平台是便捷的(可作为()来运行不需编译器。慢需要较多的代码。在使用语言上有限制。不能与编译系统或用户库()连接。只能使用预定义宏来获得结构中的数据。(见和)比运行快在使用语言上不存在限制可用任何编译器编译。能调用以其他语言编写的函数机器物理结构需要用户建立的每个版本的共享库(如果包含有注释器()不能处理得语言元素,则不能作为)运行总的来说,当决定使用那种类型的时使用作为简单的函数使用作为复杂的函数,这些函数对有较大要求(例如每次运行时,在每个单元上均须调用的属性需要使用编译库一个例子编辑代码,并且在用户的模型中有效使用它,须遵循以卜七个基木步骤定义用户模型。编制语言源代码。运行读入,并设置文件。编译或注释()语言源代码。在中激活开始计算。分析计算结果,并与期望值比较。在开始解决问题前,用户必须使用定义希望解决的问题()。例如,加入用户希望使用来定义一个用户化的边界条件()。用户首先需要定义一系列数学方程来描述这个条件。接下来用户需要将这些数学方程(概念设计)用语言写成一个函数)。用户可用文本编辑器米完成这一步。以为后缀名米把这个文件保存在工作路径下。写完语言函数后,用户即可运行并且读入或设置文件()。对C语言源代码进行注释,编译,和调试),并在中激活用户函数()。最后,运行计算(),分析结果并与期望值比较。()。根据用户对结果的分析,可将上述整个过程重复几次。具体如下。定义用户模型生成和使用的第一步是定义用户的模型方程。如图所示的涡轮叶片。模拟叶片冑围的流玚使用了非结构化网格。计算域由底端的周期性边界()延伸到顶端的相同部分(),速度入口在左边,压力出口在右边。Irvine ua〔15512405危e893nde5)文中对入口速度为常数分布和抛物线分布的流场进行了比较。分段线性的分布可由边界场选项得到(),而多项式分布则只能使用用户自定义函数得到。进口速度为常数()的结果如图1.7.2和1.7.3所示。当流动沿着涡轮叶片进行时,初始速度场被改变了。了4101265105 tces. 89E.11电+01Turbine vane〔1551cl,2405他e,的93nde引假定现在要设涡轮叶片入∏速度不是常数值,其分布如下2un=20-2000745士变量在入口中心处为0.0,在入口上部和下部则分别为而入口中心处的速度为边界上为0。用户可用描述这一分布,并将它应运到模型中来解决这类问题。编制语言源代码。选定方程定义后,用户可用任意文本编辑器来书写C语言代码。以扩展名保存源代码文件保存到工作路径下。关于的书写请参考下面是一个怎样在中应用方程的例子。的功能由主要的宏()米定义。此处,宏用来表示下面的代码旨在给求解器捉供边界的轮廓信息。书中将在以后部分讨论其它的宏宏的第一个变量用来定义速度入口面板中的函数。名称可任意指定。在给定的边界区域上的所有单元面()上将会使用函数的这个方程。当用户在用户界面选定作为边界条件时,将会自动定义线()。下标由应用程序自动定义。中被用来形成对边界区域上所有单元面的循环()。对于每个面,面的质心()的坐标可由宏来获得。抛物线方程中用到了坐标,速度值通过宏来返回给面。宏和宏都是提供的宏。详情请见5读入,并设置文件建立后,用户开始设置在工作路径下启动读入(或设置)文件(如果文件以前设置过,请确认它是否被保存在了工作路径编详或注释()语言源代码这部分将例中的源代码作为来编译。注意,这个例子不可应用于的并行网络()。完整的编译和连接请见确认的文件(如果以前设置过)和语言源代码在工作路径卜。用面板编译例如, ud fexamp leInterpreted UDFSSource file rlameample+CCPP Command nameStack sizeH Display Assembly ListingH Use Contributed CPPCompilecospHelp在下键入语言源代码文件例如, udfexample.c。如果用户源代码不在目前工作路径下,则在编译,需在面板中间如文件完整的路径。在栏里,选择预处理器。省设置为如果用户函数的局部变量数目大于将会导致堆栈溢出。这种情況下,应将设置为比局部变量大的数
- 2020-11-04下载
- 积分：1
|
联系客服（QQ：3324461878）
