fluent湍流模型

2022-09-20

第一篇:fluent湍流模型

湍流怎么造句

湍流拼音【注音】: tuan liu

湍流解释

【意思】:(tuānliú)<书>流得很急的水。

湍流造句:

1、湍流将紧挨着球的曲面,从而减少足球的空气阻力。

2、我们认为一个中等或尾流结构可能存在,现在我们可以证明有大群结构位于湍流非常中心的位置。

3、虽然其结构,被称为壁结构,已经在湍流的边缘被找到,但是一个难以捉摸的中等或尾流结构至今从未被发现。

4、在汹涌的湍流中,每个人在他们内心都应该有指导他们做出决定的思想。

5、研究小组现正期找到类似的结构,如果它们存在于其它的湍流流动的案例中。

6、由于湍流而快速变化的折射在视线中会影响到光的不同颜色,这种影响也各不相同,一般会给恒星产生一种闪烁的效果。

7、这可能包括工作机械零件,涉及血液流动的医疗,和在空中,海上和公路旅行中的湍流各个方面。

8、当鲨鱼在水中游动时,水流从鳞屑的沟槽中流过有助于减少湍流,保持其流线型的泳姿。

9、《第十三个故事》情节跌宕起伏,就像湍流的河水,充满不可预知的漩涡和大浪,让读者无法逃避。

10、如果你踢球的力量足够大,使得球表面的气流形成湍流,则阻力会很小,你很可能踢成高射炮。

11、实际上,上周经历很多湍流的航班就是沿着该高压边缘。

12、然而,当气流为湍流时,边界层维持时间较长。

13、它将测定太阳磁场形成以及如何导致太阳剧烈活动,比如太阳风湍流。

14、但如果你能大力踢球使其获得一个足够快的速度,使它表面的气流形成湍流,足球将受到较小的制动力(见上图)。

15、当球在空中速度减慢时,周围的气流从湍流变为稳定的层流。

16、这架69磅重的飞行器由一位希腊奥林匹克自行车手所驱动,在靠近圣托里尼的海岸时还遭遇到了空中湍流的袭击。

17、这一新发现的湍流状态是由大量存在于一种湍拎干结构中的元素组成的,而且已经被该研究组描述为一块“打结的漩涡挂毯”。

18、球的表面流动的空气形成湍流,这使得球的阻力相对较低。

19、混沌理论先驱BenoitMandelbrot发现尼罗河每年的洪水泛滥程度符合这个性质,音乐和空气湍流中也有这个性质。

20、现在,我们确信我们所拥有的湍流经验可以帮助消费者克服困难,并能帮助商业的成功。

21、他们站在湍流的汹涌的河水中间,束手无策,天完全黑下来。他们离河岸还有25英尺之远。

22、对于大多数危险的湍流,我们花费了更多的时间来保障安全,但是只有少数情况下,这些措施才起到重大的作用。

23、从冰川包覆的山巅冲击而下的湍流携下一种具有很高价值的玉石,毛利人将这种硬质半透明的石头雕刻成为珠宝和刀刃,既是工具也可以作为武器。

24、然而,处女是简单化的,什么东西都显现在表面一目了然,天蝎却更加注重生活表象下的湍流。

25、在绵延湍流中,享受尼泊尔宁静、与世隔绝的乡间景观。

第二篇:FLUENT

FLENT ①快捷键:

缩放:中键选择要缩放区域----松开; ②

一、导入并检查网格

1、读入网格

2、网格按比例缩放

---在scaling中选择尺寸单位(建模时所用单位)---scale--- 会发现domain extent中数值变化;

Domain extent:当前模型的尺寸范围(我们需要的尺寸范围); Scale:缩放;unscale:撤销缩放(反缩放);

3、检查网格

Mesh---check:网格最小体积大于零即可; Minimum volume(m3):

4、光顺网格与交换单元面

5、网格显示

二、选择求解器及计算模型

1、选择求解器

①压力基和密度基:

分离方法(segregated)—压力基(pressure-based)--可压缩(分压):

基于压力的算法;适用于低速流体; 耦合方法(coupled)—密度基(density-based)--不可压缩:

基于密度的算法;使用与高速流体;密度可变;密度,速度 ,压强都变化,采用耦合求解;

只有密度基的时候才有积分显式/隐式选择: Solution condition---formulation--- ②设置湍流模型:

Models---viscous-laminar--- ③选择能量方程: Models---energy-off

2、设置运行环境

①打开步骤:define---operating conditions 或cell zone conditions--- operating conditions

三、定义材料

1、创建新流体:Define---materials:

<1从database中找(常用与添加液体水water liquid); <2新定义名字,自己输入参数;

2、粘性系数:material—air—properties—viscosity 温度变化不大时,一般采用southerland模型;

四、设置边界条件

1、设置流体区域边界条件

①设置计算所用材料: ---cell zone condition---edit---material name---选择所需材料---

2、设置进口边界条件 ①设置压力入口:

入口压强定义了总压和静压,即可知道马赫数;

p121(1Ma) p02Cauge total pressure 总表压(相对压强) Initial cauge pressure 初始表压---静压

3、设置出口边界条件

4、

五、求解

1、求解参数设置:--Soulution Controls ①courant number:柯朗数;C2att 2axVmin

t:步长;Vmin:最小单元体积;

通过设置柯朗数来设置计算步长;柯朗数越小,步长越小,

越容易收敛;

②显式/隐式积分:

在选择density-based的时候在solution methods里面会出现;

显式积分(explict):C<1;一般先取0.1,如果收敛,可取0.5;

隐式积分(inplict):先取0.1,接着可取0.5,1,5···

2、打开残差图:solution---monitors <1 <2设置收敛点(降低收敛值):

residuals,statistic···--- residuals-plot(双击)---equations---convergence absolute criteria---将数值调小即可;

3、初始化

4、动画设置

5、保存文件

6、开始迭代

六、后处理

1、等值线

2、矢量图

3、检查质量流量连续性:

Report---result reports---fluxes---选择面(入出口)---compute;;

流入和流出质量会有误差,但误差应在一定范围内,例5%,超出范围时,应降低收敛点后继续计算;

4、显示迹线图:

---graphics and animation---pathlines---set up---选择release surface---display--- 显示迹线动画:单击pulse(display右侧)即可---stop停止;

第三篇:fluent学习心得

1. 分离式求解器和耦合式求解器:都适用于从不可压到高速可压的很大范围的流动,总得来说,计算高速可压时,耦合式求解器更有优势;分离式求解器中有几个模型耦合式求解器中没有,如VOF,多项混合模型等。

2. 对于绝大多数问题,选择1st-Order Implicit就已经足够了。精度要求高时,选择2st-Order Implicit.而Explicit选项只对耦合显式求解器有效。

3. 压力都是相对压力值,相对于参考压力而言。对于不可压流动,若边界条件中不包含有压力边界条件时,用户应设置一个参考压力位置。计算时,fluent强制这一点的相对压力值为0. 4. 选择什么样的求解器后,再选择什么样的计算模型,即通知fluent是否考虑传热,流动是无粘、层流还是湍流,是否多相流,是否包含相变等。默认情况,fluent只进行流场求解,不求解能量方程。

5. 多相流模型:其中vof模型通过单独的动量方程和处理穿过区域的每一流体的容积比来模拟两种或三种不能混合的流体。

6. 能量方程:选中表示计算过程中要考虑热交换。对于一般流动,如水利工程及水力机械流场分析,可不考虑传热;气流模拟时,往往要考虑。默认状态下,fluent在能量方程中忽略粘性生成热,而耦合式求解器包含有粘性生成热。

7. 粘性模型:inviscid无粘计算;Laminar模型,层流模型;k-epsilon(2 eqn)模型,目前常用模型。

8. 材料定义:比较简单 9. 边界条件:见P210-211 10. 给定湍流参数:在计算区域的进口、出口及远场边界,需给定输运的湍流参数。Turbulence specification Method项目,意为让用户指定使用哪种模型来输入湍流参数。用户可任选其一,然后按公式计算选定的湍流参数,并作为输入。 湍流强度,湍动能k,湍动耗散率e。 11. 常用的边界条件: 压力进口:适用于可压和不可压流动,用于进口的压力一直但流量或速度未知的情况。Fluent中各种压力都是相对压力值。

速度入口:用于不可压流,如果用于可压流可能导致非物理结果。 质量进口:规定进口的质量。

压力出口:需要在出口边界处设置静压。静压只用于亚音速流动。在fluent求解时,当压力出口边界上流动反向时,就是用这组回流条件。出口回流有三种方式:垂直与边界,给定方向矢量,来自相邻单元。 出流:用于模拟求解前流速和压力未知的出口边界。适用于出流面上的流动情况由区域内外推得到,且对上游没影响。不用于可压流动,也不能与压力进口边界条件一起是用。 压力远场:只适用于可压气体流动,气体的密度通过理想气体定律来计算。

12. 设置求解控制参数:为了更好的控制求解过程,需要在求解器中进行某些设置,内容包括选择离散格式、设置欠松弛因子、初始化场变量及激活监视变量等。

Fluent允许用户对流项选择不同的离散格式。默认情况下,当是用分离式求解器时,所有方程中的对流相一阶迎风格式离散;耦合式求解时,二阶精度格式,其他仍一阶。对于2D三角形和3D四面体网格,注意要是用二阶精度格式。一般,一阶容易收敛,精度差。

欠松弛因子:为了加速收敛,在迭代10次左右后,检查残差是增加还是减小,若增大,则减小欠松弛因子的值;反之,增大它。

Pressure-velocity coupling:包含压力速度耦合方式的列表。该项只在分离式求解器中出现。可选SIMPLE、SIMPLEC、PISO。多数选择simplec,piso算法主要用于瞬态问题的模拟,特别是希望使用大的时间步长的情况。

Courant Number;设置网格的Courant数,用于控制耦合求解时的时间步长。对于耦合显示求解器,该数值不要过大,一般<2。隐式求解器,可取较大值,一般取5,有时20,甚至100,也可收敛。

13. 设置监视参数,一般残差监视。

14. 初始化流场的解:向fluent提供流场的解的初始猜测值。 15. 流畅迭代计算,稳态问题求解和非稳态问题求解。

第四篇:Fluent的并行计算

Possibilities of Parallel Calculations in Solving Gas Dynamics Problems in the CFD Environment of FLUENT Software

N. Chebotarev Research Institute of Mathematics and Mechanics, Kazan State

University, Kazan, Russia Received August 25, 2008 Abstract :The results of studying an incompressible gas flow field in a periodic element of the porous structure made up of the same radius spheres are presented; the studies were based on the solution of the Navier–Stokes equations using FLUENT software. The possibilities to accelerate the solution process with the use of parallel calculations are investigated and the calculation results under changes of pressure differential in the periodic element are given. DOI : 10.3103/S1068799809010103

Multiprocessor computers that make it possible to realize the parallel calculation algorithms are in recent use for scientific and engineering calculations. One of the fields in which application of parallel calculations must facilitate a considerable progress is the solution of three-dimensional problems of fluid mechanics. Many investigators use standard commercial CFD software that provide fast and convenient solutions of three-dimensional problems in complicated fields.

The present-day CFD packages intended to solve the Navier–Stokes equations describing flows in the arbitrary regions contain possibilities of parallel processing. The objective of this paper is to test the solutions of a three-dimensional problem of gas dynamics using FLUENT software [1] by a multiprocessor computer in the mode of parallel processing.

An incompressible gas flow in the porous structure made up of closely-arranged spheres is calculated. Structures of different sphere arrangement are widely used as models of porous media in the theory of filtration. Using the porous elements it is possible to realize the processes of filtration, phase separation, throttling, including those in aircraft engineering [2]. The hydrodynamic flows in porous structures in the domain of small Reynolds numbers are described, as a rule, under the Stokes approximation without regard for the inertia terms in the equations of fluid motion [3]. At the same time, the flow velocities in the porous media may be rather large, and the Stokes approximation will not describe a real flow pattern. In this case the solution of the complete Navier–Stokes equations should be invoked. The flow with regard for inertia terms in the equations of fluid motion in the structures of different sphere arrangement was theoretically studied in [6] and experimentally in [7].

PROBLEM STATEMENT

We consider an incompressible gas flow in the three-dimensional periodic element of the porous structure made up of closely-arranged spheres of the same diameter d , the centers of which are in the nodes of the ordered grid (Fig. 1a). The porosity of the structure under consideration determined as the ratio of the space occupied by the medium to the total volume is equal to 0.26. Taking into account symmetry and periodicity of the flow, we will separate in the space between spheres the least element of the region occupied by air (Fig. 1b). In connection with a difficulty in dividing the calculation domain, small cylindrical areas are excluded in the vicinity of points at which the element spheres are in contact.

Fig. 1. Scheme of sphere arrangement (a) and a periodic element in the air space between spheres

(b).

The gas flow velocities inside the porous structures are so small that it is possible to neglect gas compressibility and adopt a model of incompressible fluid. The laminar flow of the incompressible gas is described by the stationary Navier–Stokes equations:

where u are the gas velocity vector and its Cartesian components; p is the pressure; μ and

ρ are the dynamic viscosity coefficient and air density. At the end bounds of the periodic element we lay down the conditions of periodicity

where L = d is the periodic element length along the flow (along the y axis); at the lateral faces we lay down the conditions of symmetry. The pressure at the end element bounds is described by the formula

where Δp is the pressure differential in the element limits. The conditions of symmetry are taken not only at the lateral faces but also at the upper and lower faces. On the spherical surfaces the conditions of adhesion are specified.

System of equations (1)–(2) is solved with the aid of the SIMPLE algorithm in the finite volume method in FLUENT software environment (FLUENT 6.3.26 version). For the calculation domain the irregular tetrahedral grid division is used (Fig. 2).

Fig. 2. Division of the periodic element into finite volumes.

ANALYSIS OF CALCULATION RESULTS IN THE PARALLEL MODE The calculations were carried out on the computational cluster of Kazan State University consisting of eight servers. Each server includes two AMD Opteron 224 processors with the clock rate 1.6 GHz and 2 GB of main memory. The servers operate under the control of the Ubuntu 7.10 version of the Linux operating system. The communication between the servers is based on the Gigabit Enternet technology. In the calculations the HP Message Passing Interface library (HP-MPI) delivered together with the FLUENT program is used. At the moment of experiments four servers were accessible for operation. To analyze the efficiency of parallel processing of the numerical solutions of the Navier–Stokesequations in the FLUENT environment, the calculations were performed with three variants of the grid division of the solution domain 116895, 307946, and 510889 finite volumes (variants A , B , C ). The number of iterations in all cases was taken to be equal to 760 resulting in solution convergence to10. All computation experiments were performed in the package mode.

One of the basic moments in solving problem with the use of FLUENT software in the parallel processing mode is the division of the initial domain into subdomains. In this case, each computation unit, that is, a processor is responsible for its subdomain. In dividing into subdomains FLUENT software uses the method of bisection, that is, when it is necessary to divide into four subdomains, the initial domain is first divided into two and then recursively the daughter subdomains are divided into two. If it is necessary to divide into three subdomains, the initial domain is divided into two subdomains so that one subdomain is twice as large as the other and then the larger subdomain is divided into two. FLUENT software incorporates several algorithms of bisection, and the efficiency of each algorithm depends on the problem geometry.

We studied the acceleration factor that is determined as the ratio of the calculation time t

1 by one processor to the calculation time tn by n processors ka .To provide the experiment purity, the acceleration factor was calculated four times for each case, and the calculation time obtained in each case was somewhat different. The minimal time for four calculations was chosen for the analysis. For the variants A , B , C without paralleling it amounted to 747, 2234, and 3600 s, respectively.

The dependence of ka on the number of n processors being used is given in Fig. 3. If the number of processors is small ( n < 5), the acceleration factor is about the same for all variants of the grid division and is close to the number n . As n increases, the parameter ka becomes much less than the number of processors and becomes different for different variants of the grid division. As a whole, the acceleration factor behavior with changes in the number of processors is in conformity with the theoretical concepts, and k a tends to the final limiting value as n grows. For the variant A with a lesser number of finite volumes, the efficiency of paralleling is lower than for the variants B and C .

Fig. 3. Dependence of the acceleration factor on the number of processors.

For a more detailed analysis of the calculation time with the different number of processors, we study the time t of data exchange between the processors and the factor tn of unloaded state [8] that represents the share of exchange time in the total calculation time. The less is ku , the higher is the efficiency of paralleling. The time of data exchange is determined by a share of boundary finite volumes between subdomains in the total number of finite volumes. Tables 1 and 2 present a share of boundary cells between subdomains and the values of the factor of unloaded state. It is seen that as the number of processors grows, the share of boundary volumes between subdomains increases resulting in the growth of the relative exchange time tu , that is, the factor of unloaded state. In this case, it is clear that the variants B and C are close to each other in both variants. In the variant A , the ratio of the exchange time to the calculation time increases much faster.

The flow pattern in the element under study is mainly determined by the value of pressure differential along the fluid flow. Denoting by v0 the average velocity in the inlet cross-section, we will introduce the Reynolds number Re:

At low pressure differential (small Re) a symmetric flow that is periodic in all directions is formed. The typical vector field of velocities for this range is shown in Fig. 4a. As the Reynolds number increases (pressure differential grows), the inertia effects become apparent (Fig. 4b). The inertia flow regime is characterized by the presence of complex three-dimensional vortex structures in the region between the spheres being streamlined.

At low velocities the gas flow in the porous structure is described by the Darcy law:

where v is the filtration velocity (the volumetric air flowrate per unit of time through unit area in the porous medium); kd is the permeability factor. The Darcy law establishes proportionality of the filtration velocity to the pressure gradient. In the inertia flow regime the Darcy law is violated, and to express the dependence of the filtration velocity on the pressure differential, the well-known Forhheimer formula is widely used [9]:

Fig. 4. Vector field of velocities at the periodic element faces at x = 0 (a, c) and x = d /2 (b, d) for

Re = 0.05 (a, b) and Re = 140 (c, d).

The permeability factor kd and the coefficient β are usually found by approximation of experimental results. At the same time the numerical solution of the equations for gas flow in the periodic element of the porous structure can be used for their definition. Let us write the equation for momentum variation when gas passes through the periodic element in the form [6]:

where n is the unit vector that is exterior with respect to the boundary element surfaceAf ; τ is the vector of viscous stresses. The boundary surface of the element consists of the boundary “fluid–fluid”Aff, the boundary “fluid–solid” Afs:Af=Aff+Afs

Let us rewrite Eq. (8) with regard for symmetry in the form:

where A is the area of the inlet cross-section in the periodic element. Taking into account that in connection with flow periodicity the last integral is zero, we will write Eq. (9) in the dimensionless form:

where we introduced the dimensionless parameter λ= d2Δp/ρvL. The gas velocity distribution in the porous element found from the solution of the Navier–Stokes equations makes it possible to calculate the forces acting on the spherical surfaces in the element limits, that is, the integrals in the right-hand part of Eq. (10). The dimensionless force of resistance λRe includes the resistance that is due to the normal fp and viscous f τ stresses on the spheres. Figure 5 presents the dependence of Re λ on the Reynolds number obtained from the calculations of the gas flow field at different values of pressure differential (calculated points), and the curve found by linear approximation of the calculated values Re λ at Re 4 > . The linear approximation obtained corresponds to the Forchheimer formula and is in a good agreement with the calculated data thus confirming the feasibility of the formula mentioned for description of flow characteristics in the porous media at the large Reynolds numbers. Using formulas (7) and (10), we will obtain from the calculated data the

permeability factor kd/d2=6.1 × 10-4.

Fig. 5. Dependence of the dimensionless resistance force on the Reynolds number.

Thus, using FLUENT software in the parallel calculations mode by a multi-processor computer, we solved a problem of the incompressible gas flow in the periodic element of the porous structure made up of the closely-arranged spheres. It is shown that the calculation time is reduced as the number of processors grows, and the efficiency of paralleling is higher with a larger number of finite volumes in dividing the calculation domain. The results of studying a flow field under variation of the pressure differential in the periodic element are presented. When the pressure differential increases (the large Reynolds numbers), the inertia flow regime that is characterized by a complicated pattern of vortex structures is formed in the porous element.

第五篇:fluent使用总结(本站推荐)

3.1计算流体力学基础与FLUENT软件介绍 3.1.1计算流体力学基础

计算流体力学(Computational Fluid Dynamics,简称CFD)是利用数值方法通过计算机求解描述流体运动的数学方程,揭示流体运动的物理规律,研究定常流体运动的空间物理特性和非定常流体运动的时空物理特征的学科[}ss}。其基本思想可以归纳为:把原来在时间域和空间域上连续的物理量的场,如速度场和压力场,用一系列有限个离散点上的变量值的集合来代替,通过一定的原则和方式建立起关十这些离散点上场变量之间的关系的代数方程组,然后求解代数方程组获得场变量的近似值[f=}}l

计算流体力学可以看作是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制下对流动的数值仿真。通过这种数值仿真,可以得到流场内各个位置上的基本物理量(如速度、压力、温度和浓度等)的分布以及这些物理量随时间的变化规律。

还可计算出相关的其它物理量,如旋转式流体机械的转矩、水力损失和效率等。此外,与CAD联合还可进行结构优化设计等。

过去,流体力学的研究主要有实验研究和理论分析两种方法。实验研究主要以实验为研究手段,得到的结果真实可信,是理论分析和数值计算的基础,其重要性不容低估。然}fu实验往往受到模型尺寸、流场扰动和测量精度等的限制,有时可能难以通过实验的方法得到理想的结果。此外,实验往往经费投入较大、人力和物力耗费较大及周期较长;理论分析方法通常是利用简化的流动模型假设,给出所研究问题的解析解或简化方程。然}fu随着时代的发展,这些方法已不能很好地满足复杂非线性流体运动规律的研究。理论分析方法的优点是所得结果具有普遍适用性,各种影响因素清晰可见,是指导试验研究和验证新的数值计算方法的理论基础。但是,它往往要求对计算对象进行抽象和简化,才有可能得出理论解。}fU对十非线性情况,只有少数流动才能得到解析结果。

计算流体力学方法很好地克服了前面两种方法的弱点,与传统的理论分析方法、实验研究方法一同组成了研究流体流动问题的完整体系。计算流体力学的发展,先后经历

2 FLUENT软件介绍

FLUENT软件是由美国FLUENT公司开发的著名的CFD计算分析软件,在航空、航天、透平机械、汽车、船舶、机械、化工、石化、计算机、半导体、能源、医学等领域得到了广泛的应用。能够解决流动、传热、化学反应、燃烧、多相流、旋涡流动等问题。

FLUENT软件研究的流动模型包括了定常和非定常流动,层流(包括各种非牛顿流模型),紊流(包括最先进的紊流模型),不可压缩和可压缩流动,传热和化学反应等。FLUENT软件设计基于“CFD计算机软件群的概念”,针对每一种流动的物理问题的特点,采用适合于它的数值解法在计算速度、稳定性和精度等各方面达到最佳。不同领域的计算软件组合起来,成为CFD软件群,从而高效率地解决各个领域的复杂流动的计算问题,在各软件之间可以方便地进行数值交换,采用统一的前后处理工具,省去了科研工作者在计算方法、编程、前后处理等方面投入的重复、低效的劳动,而可以将主要精力用十物理问题本身的探索上。

流体有限体积法(Finite Volume Method,简称FVM)是目前计算流体动力学领域内应用最普遍的一种对偏微分方程组的离散方法。FLUENT软件就是采用C语言编写的基于非结构化网格和有限体积法的通用CFD求解器,它推出了多种优化的物理模型,如定常和非定常流动;层流(包括各种非牛顿流模型);紊流(包括最先进的紊流模型);不可压缩和可压缩流动;传热;化学反应等。对每一种物理问题的流动特点,有适合它的数值解法,用户可对显式或隐式差分格式进行选择,以期在计算速度、稳定性和精度等方面达到最佳。

在FLUENT 5.0之后的版本中,都采用GAMBIT的专用前处理软件。GAMBIT软件是面向CFD的专业前处理器软件,它包含全面的几何建模能力,也可以从主流的CAD/CAE软件导入几何体和网格,GAMBIT强大的布尔运算能力为建立复杂的几何模型提供的极大的方便。GAMBIT功能强大的网格划分工具,可以划分出包含边界层等CFD特殊要求的高质量的网格。GAMBIT中专有的网格划分算法可以保证在较为复杂的几何区域直接划分出高质量的六面体网格。GAMBIT中的TGRID方法可以在极其复杂的几何区域中划分出与相邻区域网格连续的完全非结构化的网格,GAMBIT网格划分方法的选择完全是智能化的,在选择一个几何区域后GAMBIT会自动选择最合适的网格划分算法,使网格划分过程变得极为容易。

通用CFD软件包,用来模拟从不可压缩到高度可压缩范围内的复杂流动。由于采用了多种求解方法和多重网格加速收敛技术,因而FLUENT能达到最佳的收敛速度和求解精度。灵活的非结构化网格和基于解的自适应网格技术及成熟的物理模型,使FLUENT在转捩与湍流、传热与相变、化学反应与燃烧、多相流、旋转机械、动/变形网格、噪声、材料加工、燃料电池等方面有广泛应用。

FLUENT软件具有以下特点:

☆ FLUENT软件采用基于完全非结构化网格的有限体积法,而且具有基于网格节点和网格单元的梯度算法; ☆ 定常/非定常流动模拟,而且新增快速非定常模拟功能;

☆ FLUENT软件中的动/变形网格技术主要解决边界运动的问题,用户只需指定初始网格和运动壁面的边界条件,余下的网格变化完全由解算器自动生成。网格变形方式有三种:弹簧压缩式、动态铺层式以及局部网格重生式。其局部网格重生式是FLUENT所独有的,而且用途广泛,可用于非结构网格、变形较大问题以及物体运动规律事先不知道而完全由流动所产生的力所决定的问题;

☆ FLUENT软件具有强大的网格支持能力,支持界面不连续的网格、混合网格、动/变形网格以及滑动网格等。值得强调的是,FLUENT软件还拥有多种基于解的网格的自适应、动态自适应技术以及动网格与网格动态自适应相结合的技术;

☆ FLUENT软件包含三种算法:非耦合隐式算法、耦合显式算法、耦合隐式算法,是商用软件中最多的; ☆ FLUENT软件包含丰富而先进的物理模型,使得用户能够精确地模拟无粘流、层流、湍流。湍流模型包含Spalart-Allmaras模型、k-ω模型组、k-ε模型组、雷诺应力模型(RSM)组、大涡模拟模型(LES)组以及最新的分离涡模拟(DES)和V2F模型等。另外用户还可以定制或添加自己的湍流模型; ☆ 适用于牛顿流体、非牛顿流体;

☆ 含有强制/自然/混合对流的热传导,固体/流体的热传导、辐射; ☆ 化学组份的混合/反应;

☆ 自由表面流模型,欧拉多相流模型,混合多相流模型,颗粒相模型,空穴两相流模型,湿蒸汽模型; ☆ 融化溶化/凝固;蒸发/冷凝相变模型; ☆ 离散相的拉格朗日跟踪计算;

☆ 非均质渗透性、惯性阻抗、固体热传导,多孔介质模型(考虑多孔介质压力突变); ☆ 风扇,散热器,以热交换器为对象的集中参数模型; ☆ 惯性或非惯性坐标系,复数基准坐标系及滑移网格; ☆ 动静翼相互作用模型化后的接续界面;

☆ 基于精细流场解算的预测流体噪声的声学模型; ☆ 质量、动量、热、化学组份的体积源项; ☆ 丰富的物性参数的数据库;

☆ 磁流体模块主要模拟电磁场和导电流体之间的相互作用问题;

☆ 连续纤维模块主要模拟纤维和气体流动之间的动量、质量以及热的交换问题;

☆ 高效率的并行计算功能,提供多种自动/手动分区算法;内置MPI并行机制大幅度提高并行效率。另外,FLUENT特有动态负载平衡功能,确保全局高效并行计算;

☆ FLUENT软件提供了友好的用户界面,并为用户提供了二次开发接口(UDF); ☆ FLUENT软件采用C/C++语言编写,从而大大提高了对计算机内存的利用率。

1. Fh lent 软件的结构组成 ()1 前 处 理 器 :9幽bit 科u en t软 件包的前处理器是galllbit,galnbit具有前处 理器建模及网格划分的功能,是进行数值模拟计算前处理器 的首选。但是,gambit适合于简单模型的建立,对于复杂模 型,可以采用Pr‘ug等软件进行建模,复杂模型建模完成 后,可以导入ganlbit软件再进行网1各划分。网格划分完成 后保存dbs文件和愉出msh文件。

前处 理 阶 段需耍用户进行如下操作: 定义 计 算 域、绘制简化物理模型 对计 算 域 进行网格划分

定义 域 边 界单元的边界条件 定义 流 体 的属性参数 (2 )求 解 器 : nuent

FLUENT简介

fluent是用于计算流体流动和传热问题的程序。它提供的非结构网格生成程序,对相对复杂的几何结构网格生成非常有效。可以生成的网格包括二维的三角形和四边形网格;三维的四面体 和六面体及混合网格。fluent很能够根据计算的结果调整网格,这种网格自适应能力对于精确求解 有较大梯度的流场有很实际的作用。由于网格自适应和调整只是在需要加密的流动区域里实施,而非整个流场,因此可以节约计算时间。

一、程序的结构

fluent程序软件包由以下几个部分组成:

(1)GAMBIT—用于建立几何结构和网格的生成。

(2)FLUENT—用于进行流体模拟计算的求解器。

(3)prePDF—用于模拟PDF燃烧过程。

(4)TGrid—用于从现有的边界网格生成体网格。

(5)Filter(Translator)—转换其他程序生成的网格,用于FLUENT计算。

利用FLUENT软件进行流体的流动和传热计算的模拟计算的流程一般是,首先利用GAMBIT进行流动区 域几何形状的构建、定义边界类型和生成网格,然后将GAMBIT中的网格文件输出用于FLUENT求解器计算的格式,在FLUENT 中读取所输出的文件并设置条件对流动区域进行求解计算,最后对计算的结果进行后处理。

二、FLUENT 程序可以求解的问题

FLUENT 可以求解计算二维和三维问题,在计算过程中,网格可以自适应调整。fluent软件的应用范围非常广泛,主要范围如下:

(1)用非结构自适应网格模拟2D或者3D流场,它所使用的非结构网格主要有三角形/五边形、四边形/五边形,或者混合网格,其中混合网格有棱柱形和金字塔形。(一致网格和悬挂节点网格都可以)

(2)不可压或可压流动

(3)定常状态或者过渡分析

(4)无粘,层流和湍流

(5)牛顿流或者非牛顿流

(6)对流热传导,包括自然对流和强迫对流

(7)耦合热传导和对流

(8)辐射热传导模型

(9)惯性(静止)坐标系非惯性(旋转)坐标系模型

(10)多重运动参考框架,包括滑动网格界面和rotor/stator interaction modeling的混合界面

(11)化学组分混合和反应,包括燃烧子模型和表面沉积反应模型

(12)热,质量,动量,湍流和化学组分的控制体源

(13)粒子,液滴和气泡的离散相的拉格朗日轨迹的计算,包括了和连续相的耦合

(14)多孔流动

(15)一维风扇/热交换模型

(16)两相流,包括气穴现象

(17)复杂外形的自由表面流动

三、FLUENT程序求解问题的步骤

利用FLUENT求解问题的步骤如下:

(1) 确定几何形状生成计算网格(用GAMBIT,也可以读取其他指定程序生成的网格)。

(2)输入并检查网格。

(3)选择求解器(2D或3D)

(4)选择求解的方程(层流或是湍流、化学组分或化学反应、传热模型等),确定其他需要的模型

(5)确定流体的材料的物性

(6)确定边界的类型及其边界条件(前者在GAMBIT中确定,但在FLUENT中可以修改,后者在FLUENT中实现)

(7)条件计算的控制参数

(8)流场的初始化

(9)求解计算

(10)判断收敛

(11)保存结果并进行后处理

四、关于FLUENT求解器的说明

在打开后会出现如下对话框,对话框中各个项代表的意义是表示求解器的精度。

(1)FLUNT2D—表示二维单精度求解器;

(2)FLUENT3D—表示三维单精度求解器;

(3)FLUENT2ddp—表示二维双精度求解器; (4)FLUENT3ddp—表示三维双精度求解器。

五、FLUENT求解方法的选择

FLUENT中所涉及的求解方法有非耦合求解(segregated)、耦合隐式求解(coupled implicit)和耦合显示求解(coupled explicit)。

非耦合求解方法主要用于不可压缩或低马赫数压缩性流体的流动.耦合求解方法则可以用在高速可压缩流体。fluent默认设置为非耦合求解,但对于高速可压 流动,或需要考虑体积力的流动,求解问题时网格要比较密,建议采用耦合隐式求解方法求解能量和动量方程,可较快地得到收敛解。缺点是需求的内存比较大,大 约是非耦合求解迭代时间的1.5-2.0倍。如果必须要耦合求解,但是机器的内存不够的条件下,可以考虑用耦合显示解法器求解问题。该解法也耦合了动量、 能量及组分方程,但是内存却比隐式求解方法小。缺点是收敛的时间比较长。

而且fluent5.5以前的版本(包括5。5),其物理模型,(比如粘性流体的几个模型)都是预先设定的,所以,对于那些做探索性或者检验新方法而进行的模拟,就不适合用。

同时gambit做网格,对于粘性流体,特别是计算湍流尺度,或者做热流计算来说其网格精度一般是不可能满足 的,除非是很小的计算区域。所以,用fluent做的比较复杂一点的流场(除了经典的几个基本流场)其计算所得热流,湍流,以及用雷诺应力模拟的粘性都不 可能是准确的,这在物理上和计算方法已经给fluent判了死刑,有时候看到很多这样讨论的文章,觉得大家应该从物理和力学的本质上考虑问题。

但是,fluent往往能计算出量级差不多的结果,曾经做了一个复杂的飞行器热流计算,高超音速流场,得到的 壁面热流,居然在量级上是吻合的,但是,从计算热流需要的壁面网格精度来判断,gambit所做的网格比起壁面网格所满足的尺寸的要大了至少2个数量级, 到现在还不明白fluent是怎么搞的。

综上,如果对付老板的一些工程项目,可以用fluent对付过去,但是如果真的做论文,或者需要发表文章,除非是做一些技术性工作,比如优化计算一般用fluent是不适合的。

fluent做力的计算是很不错的,做流场结构的计算,即使得出一些涡,也不是流场本身性质的反应,做低 速流场计算,fluent的优势在于收敛速度快,但是低速流场计算,其大多数的着眼点在于对流场结构的探索,所以计算得到的结果就要好好斟酌一下了,高速 流场的模拟中,一般着眼点在于气动力的结果,压力分布以及激波的捕捉,这些fluent做的很不错。。

对于运用fluent来求解问题,首先要对本身求解的物理模型有充分的了解,只有在这个基础上,才能够选择出正确的,计算模型以及相应的边界条件。

对于fluent计算的方法,确实是采用的有限体积法,不过对基于非结构网格的5.X,觉得其采用的应该 是同位网格而不是交错网格,因为非结构网格情况下,交错网格的方法处理起来比同位网格方法要复杂很多。一般见到的非结构网格下FVM(有限体积法)多半还 是采用的同位网格而非交错网格,这个问题还可以进一步探讨。对于非结构网格而言,目前能够做到的离散精度也只能是二阶精度了,再高精度目前还没法做到,或 者说还没有做到很实用。

对于gambit做网格,确实不是十分的理想,不过这个也不能怪罪gambit,因为非结构网格的生成方法,本 身

在理论上就有一些瑕疵(姑且这样说吧,不能说是错误,呵呵)所以对于一些十分复杂,而且特殊的流场,可能最终生成的网格会很不理想,这个时候多半需要采 取一些其它的迂回的方法,例如将复杂区域分区,分成一些简单的区域,然后在简单区域里面生成网格,最后再组合,而不是将整个复杂区域教给gambit让其 一次生成网格。有时在软件做不到的地方,就需要人想法补上了。

对于壁面网格的问题,gambit中提供了生成边界层网格的方法,恩,不知道是否这个功能也同样不能满足所需。gambit中边界层网格只是在壁面法向进行特别的处理。对于壁面切向方向则是和边界层外网格尺度相当的。

对于fluent的适用范围,本身fluent是一个比较成熟的商业软 件,换句话说,其适用的数值方法,多半也是目前相对比较成熟的方法之一。因此用fluent来做工程项目确实是很适合的,因为它相对效率较高,而且实际上 fluent中有一些对特殊问题的简化处理其目的也是直接针对工程运用的。因此如果是完全的基于fluent做流场分析,然后做论文,这样是不行的。需要 强调的是,fluent仅仅是一种CFD的工具,一个相对好用的工具。

对于fluent做高速可压流动问题,由于有限体积法本身对于求解有间断(激波)的流动问题就存在一定的误差的,有限体积法实际上应该更加的适合于不可压流动问题,因为这个方法本身 的特点就保证了通量的守恒,对于不可压流动,那就是保证了整个流场的质量守恒。对于算激波的问题似乎还是得要实用一些高精度格式,例如 NND,TVD,时空守恒格式等。顺便问stipulation一个问题,在算钝头体(导弹)小攻角来流夸音速流动问题时,在计算中是否有激波的振荡现

对于旋转机械的流动问题,fluent中提供了几种方法,一种是就是很简单用坐标变换的概念化旋转为静止,然后 添加一个惯性力。一种是所谓的多参考坐标系方法,还有就是混合面方法,最后是滑移网格方法。第一种方法自不用说,理论上是精确的,后面三钟方法中, fluent中以滑移网格方法计算的准确度最好,前面两种方法都有很强的工程背景并且是在此基础上简化而来的。但这些方法的运用都有一些前提条件。

fluent公司还有另外的一个工具,MixSim是针对搅拌混合问题的专用CFD软件内置了专用前处理器,可迅速建立搅拌器和混合器的网格及计算模型。

:

解决问题的步骤

确定所解决问题的特征之后,你需要以下几个基本的步骤来解决问题: 1.创建网格. 2.运行合适的解算器:2D、3D、2DDP、3DDP。 3.输入网格 4.检查网格 5.选择解的格式

6.选择需要解的基本方程:层流还是湍流(无粘)、化学组分还是化学反应、热传导模型等 7.确定所需要的附加模型:风扇,热交换,多孔介质等。 8..指定材料物理性质 8.指定边界条件 9.调节解的控制参数 10.初始化流场 11.计算解 12.检查结果 13.保存结果

14.必要的话,细化网格,改变数值和物理模型。

想起CFD,人们总会想起FLUENT,丰富的物理模型使其应用广泛,从机翼空气流动到熔炉燃烧,从鼓泡塔到玻璃制造,从血液流动到半导体生产,从洁净室到污水处理工厂的设计,另外软件强大的模拟能力还扩展了在旋转机械,气动噪声,内燃机和多相流系统等领域的应用。今天,全球数以千计的公司得益于FLUENT的这一工程设计与分析软件,它在多物理场方面的模拟能力使其应用范围非常广泛,是目前功能最全的CFD软件。

FLUENT因其用户界面友好,算法健壮,新用户容易上手等优点一直在用户中有着良好的口碑。长期以来,功能强大的模块,易用性和专业的技术支持所有这些因素使得FLUENT受到企业的青睐。

网格技术,数值技术,并行计算

计算网格是任何CFD计算的核心,它通常把计算域划分为几千甚至几百万个单元,在单元上计算并存储求解变量,FLUENT使用非结构化网格技术,这就意味着可以有各种各样的网格单元:二维的四边形和三角形单元,三维的四面体核心单元、六面体核心单元、棱柱和多面体单元。这些网格可以使用FLUENT的前处理软件GAMBIT自动生成,也可以选择在ICEM CFD工具中生成。

在目前的CFD市场, FLUENT以其在非结构网格的基础上提供丰富物理模型而著称,久经考验的数值算法和鲁棒性极好的求解器保证了计算结果的精度,新的NITA算法大大减少了求解瞬态问题的所需时间,成熟的并行计算能力适用于NT,Linux或Unix平台,而且既适用单机的多处理器又适用网络联接的多台机器。动态加载平衡功能自动监测并分析并行性能,通过调整各处理器间的网格分配平衡各CPU的计算负载。

湍流和噪声模型

FLUENT的湍流模型一直处于商业CFD软件的前沿,它提供的丰富的湍流模型中有经常使用到的湍流模型、针对强旋流和各相异性流的雷诺应力模型等,随着计算机能力的显著提高,FLUENT已经将大涡模拟(LES)纳入其标准模块,并且开发了更加高效的分离涡模型(DES),FLUENT提供的壁面函数和加强壁面处理的方法可以很好地处理壁面附近的流动问题。

气动声学在很多工业领域中倍受关注,模拟起来却相当困难,如今,使用FLUENT可以有多种方法计算由非稳态压力脉动引起的噪音,瞬态大涡模拟(LES)预测的表面压力可以使用FLUENT内嵌的快速傅立叶变换(FFT)工具转换成频谱。Fflow-Williams&Hawkings声学模型可以用于模拟从非流线型实体到旋转风机叶片等各式各样的噪声源的传播,宽带噪声源模型允许在稳态结果的基础上进行模拟,这是一个快速评估设计是否需要改进的非常实用的工具。

动态和移动网格

内燃机、阀门、弹体投放和火箭发射都是包含有运动部件的例子,FLUENT提供的动网格模型满足这些具有挑战性的应用需求。它提供几种网格重构方案,根据需要用于同一模型中的不同运动部件,仅需要定义初始网格和边界运动。动网格与FLUENT提供的其他模型如雾化模型、燃烧模型、多相流模型、自由表面预测模型和可压缩流模型相兼容。搅拌槽、泵、涡轮机械中的周期性运动可以使用FLUENT中的动网格模型(moving mesh)进行模拟,滑移网格和多参考坐标系模型被证实非常可靠,并和其他相关模型如LES模型、化学反应模型和多相流等有很好的兼容性。

传热、相变、辐射模型

许多流体流动伴随传热现象,FLUENT提供一系列应用广泛的对流、热传导及辐射模型。对于热辐射,P1和Rossland模型适用于介质光学厚度较大的环境,基于角系数的surface to surface模型适用于介质不参与辐射的情况,DO模型(Discrete ordinates)适用于包括玻璃的任何介质。DTRM模型(Discrete ray tracing module)也同样适用。太阳辐射模型使用光线追踪算法,包含了一个光照计算器,它允许光照和阴影面积的可视化,这使得气候控制的模拟更加有意义。

其他与传热紧密相关的汽蚀模型、可压缩流体模型、热交换器模型、壳导热模型、真实气体模型、和湿蒸汽模型。相变模型可以追踪分析流体的融化和凝固。离散相模型(DPM)可用于液滴和湿粒子的蒸发及煤的液化。易懂的附加源项和完备的热边界条件使得FLUENT的传热模型成为满足各种模拟需要的成熟可靠的工具。 化学反应模型

化学反应模型,尤其是湍流状态下的化学反应模型在FLUENT软件中自其诞生以来一直占着很重要的地位,多年来,FLUENT强大的化学反应模拟能力帮助工程师完成了对各种复杂燃烧过程的模拟。涡耗散概念、PDF转换以及有限速率化学模型已经加入到FLUENT的主要模型中 :涡耗散模型、均衡混合颗粒模型,小火焰模型以及模拟大量气体燃烧,煤燃烧、液体燃料燃烧的预混合模型。预测NOx生成的模型也被广泛的应用与定制。

许多工业应用中涉及发生在固体表面的化学反应,FLUENT表面反应模型可以用来分析气体和表面组分之间的化学反应及不同表面组分之间的化学反应,以确保表面沉积和蚀刻现象被准确预测。对催化转化、气体重整、污染物控制装置及半导体制造等的模拟都受益于这一技术。

FLUENT的化学反应模型可以和大涡模拟(DES)及分离涡(DES)湍流模型联合使用,这些非稳态湍流模型耦合到化学反应模型中,才有可能预测火焰稳定性及燃尽特性。 多相流模型

多相流混合物广泛应用于工业中,FLUENT软件是在多相流建模方面的领导者,其丰富的模拟能力可以帮助工程师洞察设备内那些难以探测的现象,Eulerian多相流模型通过分

别求解各相的流动方程的方法分析相互渗透的各种流体或各相流体,对于颗粒相流体采用特殊的物理模型进行模拟。很多情况下,占用资源较少的的混合模型也用来模拟颗粒相与非颗粒相的混合。FLUENT可用来模拟三相混合流(液、颗粒、气),如泥浆气泡柱和喷淋床的模拟。可以模拟相间传热和相间传质的流动,使得对均相及非均相的模拟成为可能。

FLUENT标准模块中还包括许多其他的多相流模型,对于其他的一些多相流流动,如喷雾干燥器、煤粉高炉、液体燃料喷雾,可以使用离散相模型(DPM)。射入的粒子,泡沫及液滴与背景流之间进行发生热、质量及动量的交换。

VOF模型(Volume of Fluid)可以用于对界面的预测比较感兴趣的自由表面流动,如海浪。汽蚀模型已被证实可以很好的应用到水翼艇、泵及燃料喷雾器的模拟。沸腾现象可以很容易地通过用户自定义函数实现。 前处理和后处理

FLUENT提供专门的工具用来生成几何模型及网格创建。GAMBIT允许用户使用基本的几何构建工具创建几何,它也可用来导入CAD文件,然后修正几何以便于CFD分析,为了方便灵活的生成网格,FLUENT还提供了TGrid,这是一种采用最新技术的体网格生成工具。这两款软件都具有自动划分网格及通过边界层技术、非均匀网格尺寸函数及六面体为核心的网格技术快速生成混合网格的功能。对于涡轮机械,可以使用G/Turbo,熟悉的术语及参数化的模板可以帮助用户快速的完成几何的创建及网格的划分。

FLUENT的后处理可以生成有实际意义的图片、动画、报告,这使得CFD的结果非常容易地被转换成工程师和其他人员可以理解的图形,表面渲染、迹线追踪仅是该工具的几个特征却使FLUENT的后处理功能独树一帜。FLUENT的数据结果还可以导入到第三方的图形处理软件或者CAE软件进行进一步的分析。 定制工具

用户自定义函数在用户定制FLUENT时很受欢迎。功能强大的资料库和大量的指南提供了全方位的技术支持。FLUENT的全球咨询网络可以提供或帮助创建任何类型装备设施的平台,比如旋风分离器、汽车HVAC系统和熔炉。另外,一些附加应用模块,比如质子交换膜(PEM)、固体氧化物燃料电池、磁流体、连续光纤拉制等模块已经投入使用。

FLUENT自豪的是能持续满足广大行业客户的应用需求。客户能够得到业内最有经验的流体工程师的技术支持,以他们丰富的专业技能作为依靠。联系您当地的FLUENT分支机构,看看FLUENT能为您的工程项目提供何种帮助吧。

上一篇:整脏治乱方案及制度下一篇:争做四有老师演讲稿

本站热搜