lbm波尔兹曼算法
LBM算法原理简介

LBM算法原理
• 求解瞬态问题时,计算时间少,精度高 传统方法求解的控制方程对流项是非线性的,每一 个时间步都需要迭代收敛。传统软件计算瞬态问题 时,给定最大迭代次数,计算并未收敛,最终影响 求解精度,而且计算效率低。 LBM算法求解的Lattice-Boltzmann方程对流项是 线性的,不存在这方面的问题。 采用Powerflow计算瞬态问题,精度高,耗费时间 大为减少。
Confidential doc of TS
1 T 2
为运动粘性, T 为温度
Confidential doc of TS
LBM算法原理
• 粒子分布函数 f 满足Lattice Boltzmann方程: 1 fi x ei t , t t fi x, t fi x, t f i eq x, t
Confidential doc of TS
LBM离散方法-2维问题
• D2Q9模型
Confidential doc of TS
LBM离散方法-2维问题
• 平衡分布函数 fi eq
f
eq 0
(根据文献[3])
4 3 2 1 2 u 9 2c
fi fi
eq
1 3 9 3 2 2 1 2 ei 4 ei 2 u u u 9 c 2c 2c 1 3 9 3 2 2 1 2 ei 4 ei 2 u u u 36 c 2c 2c
Confidential doc of TS
LBM算法原理
• LBM方法将连续介质看作大量位于网格节 点上的离散流体质点粒子. 粒子按碰撞和迁 移规则在网格上运动,通过对各网格流体质 点及运动特征的统计,获得流体宏观运动 规律。 • 粒子分布函数f(r,e,t)drde表示在t时刻,在 空间r处,粒子的速度在e到e+de的粒子数 密度。
格子玻尔兹曼方法

格子玻尔兹曼方法格子玻尔兹曼方法(Lattice Boltzmann Method,简称LBM)是一种基于微观粒子动力学的计算流体力学方法,它通过模拟流体微观粒子在格子空间上的运动来描述流体的宏观行为。
相比传统的有限元方法和有限差分方法,格子玻尔兹曼方法具有较好的并行性能和适应性,特别适用于多孔介质流动、复杂边界条件下的流动以及多相流等问题的模拟。
格子玻尔兹曼方法的基本思想是将流体系统离散化为一个个小的流体微团,这些微团在空间网格上运动,并通过碰撞和迁移过程来模拟流体宏观行为。
在每个时间步长内,微团在空间网格上按照一定的规则进行迁移,并在碰撞过程中遵循玻尔兹曼方程,通过碰撞和迁移过程来模拟流体的宏观行为。
通过在空间网格上迁移和碰撞的过程,可以模拟出流体的宏观运动规律,从而实现对流体流动的模拟和计算。
格子玻尔兹曼方法的优势之一是其较好的并行性能。
由于其基于网格的离散化特性,格子玻尔兹曼方法在并行计算上具有天然的优势,能够有效地利用多核、多节点的计算资源,实现对大规模流体问题的高效模拟。
这使得格子玻尔兹曼方法在计算流体力学领域得到了广泛的应用,特别是在大规模流体模拟和高性能计算方面具有很大的优势。
另外,格子玻尔兹曼方法在处理复杂边界条件和多相流问题上也具有一定的优势。
由于其基于微观粒子动力学的特性,格子玻尔兹曼方法能够比较灵活地处理复杂的边界条件,如固体边界、移动边界等,同时也能够较为方便地模拟多相流体的运动,包括气液两相流、多组分流体等,这使得格子玻尔兹曼方法在工程领域的应用具有广阔的前景。
总的来说,格子玻尔兹曼方法作为一种基于微观粒子动力学的计算流体力学方法,具有较好的并行性能和适应性,特别适用于多孔介质流动、复杂边界条件下的流动以及多相流等问题的模拟。
它在大规模流体模拟和高性能计算方面具有很大的优势,同时也能够比较灵活地处理复杂的边界条件和多相流问题,因此在工程领域具有广泛的应用前景。
格子玻尔兹曼方法的发展将为流体力学领域的研究和工程应用带来新的机遇和挑战。
格子玻尔兹曼方法及其在微通道绕流中的应用

格子玻尔兹曼方法及其在微通道绕流中的应用格子玻尔兹曼方法(LBM)是一种基于格子模型和玻尔兹曼方程的流体力学仿真方法。
相比于传统的Navier-Stokes方程求解方法,LBM具有更强的并行计算能力和数值稳定性,因此在微通道绕流等流体力学问题的数值模拟中得到了广泛的应用。
LBM的基本原理是将流体划分为一系列的格子点,每个格子点上有一个分布函数,该函数描述了在该点上的流体微粒的速度和密度。
通过在每个时间步中更新这些分布函数,可以计算出流体的速度场和密度分布。
在LBM中,流体微粒只在一个离散的速度集合中进行碰撞和弛豫过程,在碰撞过程中,微粒的速度和密度会根据玻尔兹曼方程进行更新,而在弛豫过程中,微粒的速度和密度会收敛到平衡态。
通过迭代求解所有格子点上的分布函数,可以得到整个流体域的速度场和密度分布。
LBM在微通道绕流中的应用主要包括两个方面:流动行为的模拟和微通道的设计优化。
在流动行为的模拟方面,LBM可以用来研究不同条件下微通道中流体的流动行为。
通过调节微通道尺寸、入口边界条件和流体的性质等参数,可以模拟和分析在不同流速和黏度条件下微通道中的流动行为。
例如,可以研究微通道中的速度分布、压力损失、剪切层和流动不稳定等现象,从而为微通道的设计提供理论依据。
在微通道的设计优化方面,LBM可以用来研究和优化微通道的几何形状和结构。
通过在LBM中引入边界条件和障碍物,可以模拟不同形状和结构的微通道,并通过优化算法来寻找最佳的流体设计。
例如,可以通过改变微通道的形状、尺寸和结构来获得更好的流体传热效果或流体混合效果。
此外,LBM还可以用来研究和优化微通道的表面润湿性,从而实现更好的流体控制和微流控操作。
总之,格子玻尔兹曼方法(LBM)是一种在微通道绕流中广泛应用的流体力学仿真方法。
它通过离散化流体微粒的速度和密度,在网格上更新分布函数来模拟流体的速度场和密度分布。
LBM不仅可以用来模拟微通道中的流动行为,还可以用来研究和优化微通道的设计和操作。
浸润边界法和格子玻尔兹曼方法

浸润边界法(immersed boundary method)是一种非边界贴合方法,在流固耦合问题中得到了广泛的应用。
该方法需要两套网格,流场使用固定直角坐标网格来求解,而浸入边界则用拉格朗日点来标识,物体与流场的作用通过两套网格之间的信息交互传递来完成。
传统的浸入边界法主要分为直接力法、惩罚力法以及动量变化法。
格子玻尔兹曼方法(Lattice Boltzmann Method,简称LBM)是一种CFD算法,可求解流动、传热等常见CFD问题。
LBM基于格子玻尔兹曼方程(LBE),从介观尺度描述了流体运动。
LBE的通用表达形式为:式中,左边为迁移项(streaming term),右边为碰撞项(collision term),fi 为粒子分布函数。
对粒子分布函数进行积分处理,可得流体密度、宏观流体速度、流体压力等宏观物理量。
lbm 高超声速计算

lbm 高超声速计算
LBM(Lattice Boltzmann Method)是一种基于微观粒子动力学的流体动力学模拟方法,它可以用于模拟高超声速流动。
在高超声速流动中,流体的速度远远超过声速,因此需要考虑诸如激波、脱离层等复杂的流动现象。
LBM作为一种基于格子的方法,可以模拟这些复杂的流动现象。
要进行高超声速流动的LBM模拟,首先需要选择适当的离散速度模型和格子类型,以及相应的边界条件。
对于高超声速流动,通常会选择D3Q27格子模型,它包含27个离散速度方向,能够更好地描述流体的运动。
在进行高超声速流动的LBM模拟时,需要考虑流体的压力、密度、温度等物理量的耦合,以及化学反应等因素。
此外,还需要考虑流体与固体或流体与流体的相互作用,以及可能存在的激波、脱离层等现象对流动的影响。
在实际计算过程中,需要考虑模拟的精度和计算的稳定性,选择合适的时间步长和网格分辨率,以确保模拟结果的准确性和可靠性。
同时,还需要考虑并行计算的方法,以提高计算效率。
总之,高超声速流动的LBM计算涉及到多个方面的物理和数值计算问题,需要综合考虑流体动力学、热力学、化学反应以及数值模拟等知识,以及计算机科学和并行计算技术,才能进行全面、准确的模拟。
波尔兹曼算法lbm单位换算

1.对于无量纲系统考虑无量纲系统:假定格子数为N,时间步数为N t,离散空间和时间尺度δx和δt分别给出如下:δx=1/N , δt=1/N t.此时,u p和u l的关系给出如下:u p=δxδtu l , C p=δxδtC l.同时,我们有νp=1/Re=δx2δtνl.那么,格子速度和黏性决定如下:u l=δtδxu p, νl=δtδx21Re.由于无量纲速度为1, 所以我们有u p=1 , u l=δtδx.对于,真实物理系统如法炮制,只需根据物理u p,重新定义νp,然后得到u l.2.真实物理系统网友问题物理参数为:P in=101325Pa P out=101225Pa, 压强差为100Pa Lx=25mm Ly=1mm,求出特征速度=13.9m/s, Re=920 空气粘度为1.5X(10^-5),划分为250X10问题:如何确定格子参数?<1>. 确定特征尺度,特征速度和Re:对于压力驱动的流动,对于此问题(Posieuille flow),如果为层流,由压力差可以确定特这速度U(y)=ΔP p2η((L/2)2−y2),y=±L/2(L=Ly)由此可以计算出入口的平均速度U mean,p.U mean,p=∫L/2−L/2U(y)d y/L.<2>. 首先确定格子解析度假定格子解析度为N ,那么物理空间网格尺度为δx=L/N. 根据δx可以确定其他方向的格子解析度。
<3>. 压力差转换. 在已知了ρp(物理密度),压力差的转换关系可以下面来定义ΔP pρp U2mean,p/2=ΔP lρl U2mean,l/2<4>. 确定物理时间步长.假定物理时间步长为δt,(C p和C l为物理声速和格子声速),那么物理速度和格子速度的关系为u p=δxδtu l.如果也考虑格子声速和物理声速的关系,也需要下列关系成立C p=δxδtC l.由Re可得物理黏性和格子黏性之间的关系,νp=U mean.p LRe=δx2δtνl.如果关心Ma数或者弱可压缩性,可以根据物理Ma数确定格子速度如下:u l=u p C p C l.否则,你可以自己定义低Ma对应u l减小可压缩性误差. 如此, 时间步长可以确定δt=δxu l u p.实际物理单位:L(长度),U(速度),C(声速) ,ν(黏度),T(时间)对应的格子单位:L l,U l,C s(格子速度),νl,T l对应关系:L r=LL l,U r=UU l=CC s,νr=ννl,T r=TT l因为Re相同,所有:Re=ULν=U l L lνl(1),对于一个已知的格子模型,已知量有L l(等于N*δx,δx一般取1,所以L l等于N),νl(可由式子ν=C s(τρ−0.5)获得),ν(可查表获得)且有(1)式可知ννl=ULU l L l=U r L r (2)U r=UU l=CC s (3)该式子的结果是已知的,因为实际声速可查表获得,而格子声速也是固定的。
格子玻尔兹曼方法

格子玻尔兹曼方法格子玻尔兹曼方法(Lattice Boltzmann Method,简称LBM)是一种基于微观粒子动力学的计算流体力学方法,它是由Lattice Gas Automata(LGA)经过演化和发展而来的。
LBM是一种离散的方法,它通过在空间网格上模拟分子碰撞和传输过程来描述流体的宏观运动。
与传统的有限差分法、有限体积法相比,LBM具有计算效率高、并行性好、适应复杂边界条件等优点,因此在流体力学领域得到了广泛的应用。
LBM的基本思想是将流体系统离散化,将连续的流体宏观运动转化为离散的微观碰撞和传输过程。
在LBM中,流体被看作是由大量微观粒子组成的,这些微观粒子在空间网格上按照一定的规则进行碰撞和传输。
通过对微观粒子的运动状态进行统计,可以得到流体的宏观性质,如密度、速度等。
LBM的核心是格子玻尔兹曼方程(Lattice Boltzmann Equation,简称LBE),它描述了微观粒子在空间网格上的运动规律。
在LBM中,流体的宏观性质由分布函数来描述,分布函数是表示在某一时刻某一空间点上流体微观粒子的分布情况。
在每个时间步内,分布函数按照一定的规则进行碰撞和传输,通过迭代计算可以得到流体在空间网格上的演化过程。
LBM的计算过程可以并行化,因此在计算效率上具有明显的优势。
LBM的另一个优点是它对复杂边界条件的处理能力强。
由于LBM是基于离散网格的方法,因此可以比较容易地处理复杂的边界条件,如曲面边界、移动边界等。
这使得LBM在模拟复杂流体系统时具有一定的优势。
除此之外,LBM还有一些其他的优点,如对多相流、多孔介质流动等复杂流体现象的模拟能力强,对于非稳态流动和湍流流动的模拟也有一定的优势。
总之,格子玻尔兹曼方法作为一种新兴的计算流体力学方法,具有诸多优点,逐渐得到了流体力学领域的广泛关注和应用。
随着计算机硬件性能的不断提升,LBM的应用前景将更加广阔,相信它会在流体力学领域发挥越来越重要的作用。
格子玻尔兹曼方法顶盖驱动流

格子玻尔兹曼方法顶盖驱动流格子玻尔兹曼方法(LBM)是一种近年来在流体力学模拟中被广泛应用的数值模拟方法。
该方法可以在复杂的几何和边界条件上精确地模拟各种流体现象。
其中的顶盖驱动流模拟,是一种受到广泛关注的研究领域。
1.什么是顶盖驱动流?顶盖驱动流是指由上端盖板施加外力所形成的流场,它是在一个密闭的正方形边界中进行模拟。
流场通常是由重力和顶部盖板的驱动力组成的。
对于使用LBM计算的顶盖驱动流,由于模型简单,计算效率高,因此已取得了广泛的应用。
2.LBM是如何模拟顶盖驱动流的?LBM的基本原理是通过在大量的离散速度上进行策略性模拟,以模拟物质运动。
与传统的流体动力学方法不同,LBM使用离散的速度和密度来描述流体的运动,因此它可以非常方便的用于复杂的几何和边界条件下的流体模拟。
在顶盖驱动流模拟中,LBM将二维正方形边界分为许多离散化的小单元格,并在每个单元格上施加离散速度和压力值来模拟流体行为。
3.顶盖驱动流模拟的应用领域是什么?顶盖驱动流模拟在许多领域具有广泛的应用,包括地质、生物学、工程学和环境科学等领域。
在地质学中,它被用于模拟岩石的岩石圈运动并对流体流动的地质效应进行分析。
在工程学中,它被用于模拟汽车空气动力学和水力学作用以及结构物的振动和熱传导等现象。
4.顶盖驱动流模拟存在的挑战是什么?尽管顶盖驱动流模拟为模拟流体行为提供了强有力的工具,但仍然存在一些挑战。
例如,该方法需要高度离散化的速度空间和网格结构,这可能会导致计算效率低下和计算成本高昂。
此外,顶盖驱动流模拟还需要对物理设置和计算参数进行大量的调整和测试,才能使模拟结果更加准确和可靠。
总之,格子玻尔兹曼方法顶盖驱动流是一种新兴的数值模拟方法,在流体力学和其他领域的应用越来越广泛。
将来,随着计算机硬件和软件的不断发展,顶盖驱动流模拟将进一步提高计算精度和计算效率,为工程学、生物学和环境科学等领域能提供更准确的解决方案。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
波尔兹曼方法基本原理格子Boltzmann 方法是使用简单的微观模型来模拟流体的宏观行为的一种新的方法。
格子Boltzmann 方法是建立在微观粒子运动论基础上的数值计算方法。
其求解过程一般需要通过编程来实现!一般来说研究流体的行为有两种方法:一种是从宏观的角度出发,假设流体连续分布于整个流场,注入密度、速度、压力等物理量均是时间可空间的足够光滑的函数。
另一种是从微观的角度,从非平衡统计力学的观点出发,假设流体是由大量的微观的例子组成,这些例子遵守力学定律,同时服从统计定律,运用统计的方法来讨论流体的宏观性质。
然而流体是由大量的粒子组成的,当我们从宏观的角度研究流体行为的时候,并没有涉及到单个粒子的行为。
通常我们所感兴趣的事代表某个点的宏观量,例如密度、速度、压力。
根据连续性假设我们可以推导出N-S 方程,并且利用数学上的微积分知识来求解,然而由于N-S 方程是高度非线性化的偏微分方程,仅仅一些具有简单变界或者比较严格物理闲着的现象才能够得到理论分析界,如果从微观的角度了研究单个粒子的真是行为,对于一个包含大量例子的系统来说粒子的运动方程往往是得不到解的。
统计学可以考虑整个系统所有的状态以及处理这个状态的概率来解决这些困难,对于稀薄气体所得到的就是Boltzmann 方程,但是得到的方程还不够,我们还要借助于统计方法得到流体的宏观性质,这就要求解Boltzmann 方程,然而Boltzmann 方程是一非线性微分方程,一般情况下严格求解也是非常困难的。
格子气方法是近年来发展起来的模拟流体力学以及其他系统的比较新的方法,格子气自动机模拟流场,就是将流体及其存在的时间和空间完全离散,给出离散的流体粒子之间相互作用以及迁移的规则。
流体只存在于空间网格上,用一系列布尔变量,.....,2,1)(,(b i t x n i =来描述在时刻t 位于x 处节点的每一个速度方向是否有粒子存在,其中b 表示每一个节点的速度方向的数目,粒子在每一个时间步长的演化包括两部分:()a 迁移,粒子沿它的速度方向向距离最近的节点运动;()b 碰撞,当不同的粒子同时到达某个节点时,按照一定的碰撞规则发生碰撞并改变运动的方向,格子气模型具有两重意义:()a 尽可能建立一个简单的模型是指能够用来模拟一个有大量粒子组成的系统;()b 反映粒子真实碰撞的本质,这样经过长时间我们可以获得流体的宏观特性。
粒子的演化过程能够用来模拟宏观的流体过程是基于下列事实,即流体的宏观特性是系统内大量粒子整体行为的结果。
分子之间的相互作用可以改变流体的传输特性,比如粘度,但是并不改变宏观方程的基本形式。
格子气的HPP 模型与FPH 模型HPP 模型将流体存在的空间划分为间距为单位长度的正方形网格,将流体想象成许多有质量没有体积的微小粒子组成,在同一时刻同一网格节点上,每一个速度方向最多允许存在一个粒子,每个粒子可以向四个方向的其中之一运动,并且遵守以下碰撞准则:当且仅当只有两个粒子沿相反方向达到某节点时(对头碰撞),它们沿另外的两个方向离开该节点,其他情形则直接穿透,PHP 模型则是将流场划分为间距为单位长度的正三角网格,并且增加了相应的碰撞准则。
格子气的微观方程为简单起见,以HPP 模型为例,用()x ,t n i 代表在时刻t 位置x 处的节点上第i 个方向的粒子数,则整个布尔场的更新可以写成()()()()231312,1++++++-Λ-ΛΛ-Λ-ΛΛ-Λ=++i i i i i i i i i i i n n n n n n n n n e x t n ν写成代数式即:()()i i i x t n e x t Ω+=++,,1n ii Ω碰撞因子从微观方程出发利用多尺度分析以及Chapmen_Enskog 展开可以得到N-S 方程。
格子波尔兹曼模型由于格子气方法的节点用一系列的0或者1来表征每个速度方向上是否存在粒子所以对于格子气的运算只涉及个速度方向上是否存在粒子,对于格子气的运算只涉及衡分布为Fermi-Dirac 分布,从而导致非伽利略不变性,既N-S 方程对流项前面有一个依赖于密度的系数(应该为1),压力依赖于宏观速度(通常的压力只依赖于密度和温度,既p=p(ρ,T)),随即统计噪声比较大(因为变量为整形)等。
为了克服上述缺陷,McNamara 与Zanetti 提出直接由Boltz-mann 方程去替换格子气自动机,简单的用单粒子分布函数f 去代替布尔变量n ,他们求Ωi(f)的方法是根据碰撞规则表求出Ωi(n),然后用f 代替n,这样波尔兹曼方程就变为:()()()f x t f e x t i i Ω+=++,,1f i把原来的整数运算变成了实数运算(f 是一个0~1之间的实数),解决了随机噪声问题。
后来,Higuera 、Jimenez 和Succi 引用线性碰撞算子。
两者的平衡分布仍然是Fermi_Di-rac 分布,没有解决其他两个问题。
后来在1991~1992年间陈十一以及钱跃宏提出了基于单一松驰时间模型来简化碰撞函数,提出了各向同性,满足伽利略不变性和使压力与速度无关的平衡分布,使问题得到解决,通常将这个模型称为LB-BGK 模型。
即将()f i Ω用()eq i i f f -τ1代替,其中eqi f 局部平衡分布函数,τ称为驰豫时间,反映非平衡态趋向平衡态的快慢。
于是LB-BGK 模型为:()()()eq i i i i f f x t f e x t -+=++τ1,,1f i 对于两维的情况一般的采取九点格子模型,其中每个节点上允许一个静止粒子存在,加上与其相邻得有8个节点记为D2Q9型。
平衡分布函数采取钱跃宏提出的: ⎪⎭⎫ ⎝⎛-=223194u f eq i ρ ()()⎥⎦⎤⎢⎣⎡-*+*+=-225291913191u u e u e f i i eq ρ ()()⎥⎦⎤⎢⎣⎡-*+*+=-2296919131361u u e u e f i i eq ρ 宏观参数为:∑=i f ρi i e f m ∑=ρ()23u e f RT i ii -=∑ρ边界条件的处理边界问题是流体力学数值计算当中相当复杂的问题,传统的数值计算方法对于边界的处理非常的困难。
格子波尔兹曼方法的边界条件处理则相对简单的多。
下面给出几种比较典型的边界问题的处理方法:刚体固壁边界条件在粘性流体力学中,边界速度为零,即切向速度和法向速度都零,因此我们可以采取下面得处理方式:在t=n-时刻以速度u 入射进来的粒子度u 入射进来的粒子度u 入射进来的粒子自由滑移边界条件流体力学中常见的一种边界条件是开口边界,如矩形空腔流等等。
在这类边界情况下,切向速度为常数,而法向速度为零。
对于这样的边界条件一般采取下面的处理方式:即把边界看成一个镜面,将入射进来的粒子沿和入射方向相对称的方向折回这样的话,假设速度为()v s u u u ,=)的粒子入射进来,以()v s u u u -=,反射出去,则平均起来看,恰好有s s u n u =*以及0=*v n u 其中v s v s n n u u 分别为粒子切向速度。
粒子法向速度,边界的单位切向速度,边界的单位法向速度边界处理方法1反弹格式反弹格式主要用于无滑移壁面。
所谓反弹格式,是指如果一个流体节点上的粒子沿粒子运动方向流动一步达到边界节点,则该粒子沿原方向反弹回原流体节点。
LBM 方程可以演化为碰撞和流动两个过程。
2插值外推格式插值外推格式是根据Taylor 展开,利用内部流体点的分布函数,根据阶数似方法得到边界上的分布函数。
非平衡态外推方法非平衡态外推方法是根据宏观物理量与分布函数确定未知量的方法。
包括压力边界和速度边界。
格子波尔兹曼方法在各个领域的应用1用于障碍绕流尾流模拟2 Benard 热对流现象的模拟3 空腔流的模拟4 在两相流方面的应用5 激波模拟综上所述,我们可以知道格子波尔兹曼算法是一种数值计算方法,它可以利用matlab编程来解决相关的数学方程的求解问题。
Ansys Flotran CFD分析简介Ansys10.0程序中的FLOTRAN CFD分析模块的功能,是一个分析二维记三维流体流体场的工具,,使用ansys1.0中用于FLOTRAN CFD分析的FLUID 141及142单元可以解决如下问题:作用于气动翼型上的升力和阻力。
超声速喷管中的流场。
弯管中流体的复杂的流体三维流动。
FLOTRAN CFD分析的种类层流活湍流。
传热或者绝热。
可压缩或不可压缩。
牛顿流体或非牛顿流体。
多组分析传输。
FLOTRAN CFD分析的主要步骤:1.确定问题的区域必须确定所分析问题的明确范围,将问题的边界设置摘田间已知的地方,如果不能知道精确的边界条件而必须假设是,不要将分析的边界条件设在靠近感兴趣的地方,也不要将边界条件设在求解变量变化梯度较大的地方。
2.确定流体的状态用户需要估计流体的特征,流体的特征是流体性质、几何边界以及流场速度幅值的函数。
FLOTRAN 能求解的流体包括气体和流体,其性质可随温度而发生显著的变化。
FLOTRAN 中的气体只能是理想气体。
再大多数的情况下,近似认为流体性质是常数,既不随温度而变化,都可以得到精确的解。
通常用雷诺数来判断流体是层流还是湍流。
通常用马赫数来判断流体是否可压缩。
流场中任意一点的马赫数是该点的速度与该点的声速的比值。
当马赫数大于0.3是考虑用可压缩算法求解。
3.生成有限元网格用户必须事确定流场中哪个地方流体的梯度变化较大,在这些地方网格必须做适当的调整。
4.施加边界条件可在网格划分之前或之后对模型施加边界条件,此时要将模型所有的边界条件都考虑进去,如果有某个相关变量的条件没有加上去,则改变量沿边界法相值得梯度将被假设为0。
5.设置FLOTRAN CFD分析参数流体性质等项目的设置,是与所求解的流体问题的类型是相关的。
6求解通过在观察求解过程中的相关变量的改变率,可以建树求解的收敛性以及稳定性。
这些变量包括速度、压力、温度、动能(ENKE自由度)和动能耗散率(ENDS自由度)等等湍流量以及有效粘性(EVIS).7检查结果以及后处理综上所述我们可以知道,ansys是可以用来解决流体运动的及其相关参数的软件,通过建模设置单元类型,模型相关的流体和热的相关特征参数,合理划分网格设置求解参数,可以得到相应的结果!。