计算流体力学读书报告

发布 2019-07-18 17:37:37 阅读 4344

5 它不像物理模型实验一开始就能给出流动现象并定性地描述,往往需要由原体观测或物理模型试验提供某些流动参数,并需要对建立的数学模型进行验证。

6 程序的编制及资料的收集、整理与正确利用,在很大程度上依赖于经验与技巧。

7 因数值处理方法等原因有可能导致计算结果的不真实,例如产生数值粘性和频散等伪物理效应。

8 cfd因涉及大量数值计算,因此,常需要较高的计算机软硬件配置。

理论分析成本最低、结果最理想、影响因素表达清楚;缺点:局限与非常简单的问题。数值方法成本较低:

数值实验、适用范围宽;缺点:可靠性差,表达困难。实验测量可靠、成本高。

将三种方法有机结合,互为补充,必然会取得相得益彰的效果。

经过四十多年的发展,cfd出现了多种数值解法。这些方法之间的主要区别在于对控制方程的离散方式。根据离散的原理不同,cfd大体上可分为三个分支:

有限差分法(finite different method,fdm); 有限元法(finite eiement method,fem);有限体积法(finite volume method,fvm)。

有限差分法是应用最早、最经典的cfd方法,它将求解域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。求出差分万程组的解,就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题的近似数值解法。

这种方法发展较早,比较成熟,较多地用于求解双曲型和抛物型问题。在此基础上发展起来的方法有pic(particle-in-cell)法、mac(marker-and-cell)法,以及南美籍华人学者陈景广提出的有限分析法(finite analytic method)等。

有限元法是20世纪80年代开始应用的—种数值解法,它吸收了有限差分法中离散处理的内核,又采用了变分计算中选择逼近函数对区域进行积分的合理方法。有限元法因求解速度较有限差分法和有限体积法慢,因此应用不是特别广泛。在有限元法的基础上,英国c a.bbrebbia等提出了边界元法和混合元法等方法。

有限体积法是将计算区域划分为一系列控制体积,将待解微分方程对每一个控制体积积分得出离散方程。有限体积法的关键是在导出离散方程过程中,需要对界面上的被求函数本身及其导数的分布作出某种形式的假定。用有限体积法导出的离散方程可以保证具有守恒特性,而且离散方程系数物理意义明确,计算量相对较小。

2023年,在其专著《numericacl heat transfer and fluid flow》中对有限体积法作了全面的阐述。此后,该方法得到了广泛应用,是目前cfd应用最广的一种方法。当然,对这种方法的研究和扩展也在不断进行,如pchow提出了适用于任意多边形非结构网格的扩展有限体积法。

有限差分法(finite difference method):用差商与代替导数;经典、成熟;数学理论基础明确;主导方法。

有限元法(finite element method):将求解区域分成若干个小的单元(element);设定待求变量在单元上的分布函数;适应性强,适用于复杂的求解区域;一度有取代有限差分法的趋势;程序技巧要求高;数学基础不如有限差分法明确。

边界单元法(boundary element method):对数学模型在边界上离散化;基于数学模型的基础解;不需要全区域求解;数学技巧要求高;通用性差;数学基础不是非常明确。

样条边界单元法(sample spectrum):改进的边界单元法;用样条插值解决边界元的基础解问题;应用范围大大拓宽;灵活性更强;缺点:通用性差、数学基础不是非常明确。

有限分析法(finite analytical method):将求解区域分成若干个子区域;给出在各个子区域上的分析解;利用边界条件耦合各个子区域上的分析解从而得到离散化方程;最大限度地引入了分析解的成分;一般可以提高求解效率和精度;数学技巧非常高;与问题的性质有关;很难形成通用程序。

数值积分变换法(numerical integration transform method):将积分变换法引入各类问题的求解;将问题进行分解,可以得到分析解的辅助问题;多个(无限多个)常微分方程;无需整体求解;数学要求高;前期准备工作量非常大;很难形成通用的求解程序。

2. 计算流体力学的步骤。

计算流体力学的求解步骤大致可以概括为如下十步:

1 给出物理模型。

2 借助基本原理/定律给出数学模型。

如质量守恒、能量守恒、动量守恒等定律。

物理模型是指把实际的问题,通过相关的物理定律概括和抽象出来并满足实际情况的物理表征。

数学模型就好理解了,就是对物理模型的数学描写。

比如n-s方程就是对粘性流体动力学的一种数学描写,值得注意的是,数学模型对物理模型的描写也要通过抽象,简化的过程。

3 确定边界条件与初始条件。

初始条件与边界条件是控制方程有确定解的前提,控制方程与相应的初始条件、边界条件的组合构成对一个物理过程完整的数学描述。 初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况。对于瞬态问题,必须给定初始条件。

对于稳态问题,不需要初始条件。 边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于任何问题,都需要给定边界条件。

对于初始条件和边界条件的处理,直接影响计算结果的精度。

4 划分计算网。

采用数值方法求解控制方程时,都是想办法将控制方程在空间区域上进行离散,然后求解得到的离散方程组。要想在空间域上离散控制方程,必须使用网格。现已发展出多种对各种区域进行离散以生成网格的方法,统称为网格生成技术。

不同的问题采用不同数值解法时,所需要的网格形式是有一定区别的,但生成网格的方法基本是一致的。目前,网格分结构网格和非结构网格两大类。简单地讲,结构网格在空间上比较规范,如对一个四边形区域,网格往往是成行成列分布的,行线和列线比较明显。

而对非结构网格在空间分布上没有明显的行线和列线。

5 建立离散方程。

对于在求解域内所建立的偏微分方程,理论上是有真解(或称精确解或解析解)的。但由于所处理的问题自身的复杂性,一般很难获得方程的真解。因此,就需要通过数值方法把计算域内有限数量位置(网格节点或网格中心点)上的因变量值当作基本未知量来处理,从而建立一组关于这些未知量的代数方程组,然后通过求解代数方程组来得到这些节点值,而计算域内其他位置上的值则根据节点位置上的值来确定。

由于所引入的应变量在节点之间的分布假设及推导离散化方程的方法不同,就形成了有限差分法、有限元法、有限元体积法等不同类型的离散化方法。

6 离散初始条件和边界条件。

前面所给定的初始条件和边界条件是连续性的,如在静止壁面上速度为0,现在需要针对所生成的网格,将连续型的初始条件和边界条件转化为特定节点上的值,如静止壁面上共有90个节点,则这些节点上的速度值应均设为0。这样,连同在各节点处所建立的离散的控制方程,才能对方程组进行求解。

在商用cfd软件中,往往在前处理阶段完成了网格划分后,直接在边界上指定初始条件和边界条件,然后由前处理软件自动将这些初始条件和边界条件按离散的方式分配到相应的节点上去。

7 给定求解控制参数。

在离散空间上建立了离散化的代数方程组,并施加离散化的初始条件和边界条件后,还需要给定流体的物理参数和紊流模型的经验系数等。此外,还要给定迭代计算的控制精度、瞬态问题的时间步长和输出频率等。

在cfd的理论中,这些参数并不值得去**和研究,但在实际计算时,它们对计算的精度和效率有着重要的影响。

8 求解离散方程。

在进行了上述设置后,生成了具有定解条件的代数方程组。对于这些方程组,数学上已有相应的解法,如线性方程组可采用guass消去法或guass-seidel迭代法求解,而对非线性方程组,可采用newton-raphson方法。在商用cfd软件中,往往提供多种不同的解法,以适应不同类型的问题。

这部分内容,属于求解器设置的范畴。

9 判断解的收敛性。

对于稳态问题的解,或是瞬态问题在某个特定时间步上的解;往往要通过多次迭代才能得到。有时,因网格形式或网格大小、对流项的离散插值格式等原因,可能导致解的发散。对于瞬态问题,若采用显式格式进行时间域上的积分,当时间步长过大时,也可能造成解的振荡或发散。

因此,在迭代过程中,要对解的收敛性随时进行监视,并在系统达到指定精度后,结束迭代过程。这部分内容属于经验性的,需要针对不同情况进行分析。

10 显示和输出计算结果。

线值图:在二维或三维空间上,将横坐标取为空间长度或时间历程,将纵坐标取为某一物理量,然后用光滑曲线或曲面在坐标系内绘制出某一物理量沿空间或时间的变化情况。

矢量图:直接给出二维或三维空间里矢量(如速度)的方向及大小,一般用不同颜色和长度的箭头表示速度矢量。矢量图可以比较容易地让用户发现其中存在的旋涡区。

计算流体力学的步骤。

3. 模型方程。

模型方程既要简单又要反映流体特征。

1 n-s方程等。

流体在静止时虽不能承受切应力,但在运动时,对相邻两层流体间的相对运动,即相对滑动速度却是有抵抗的,这种抵抗力称为粘性应力。流体所具有的这种抵抗两层流体间相对滑动速度,或普遍说来抵抗变形的性质,称为粘性。

n-s方程:

输运方程:对流方程:

2 流体热传导及扩散。

除了粘性外,流体还有热传导(heat transfer)及扩散(diffusion)等性质。当流体中存在着温度差时,温度高的地方将向温度低的地方传送热量,这种现象称为热传导。同样地,当流体混合物中存在着组元的浓度差时,浓度高的地方将向浓度低的地方输送该组入的物质,这种现象称为扩散。

乡土中国读书报告

第九章 无讼 点明了中国在从乡土社会蜕变的过程中法治秩序与礼治秩序的矛盾,现行的司法制度在乡下发生了很多的 第十章 无为政治 论述了社会冲突中的 横暴权力 与社会合作中的 同意权力 但在农业性的乡土社会中 无为 的政治是最现实也是最理想的。第十一章 长老统治 指明在中国传统乡土社会的权力结构中,与 ...

读书报告格式与规范

关于读书报告的格式和规范。读书报告格式 1 引言。简要交代本文撰写的缘起,开宗明义点明话题。2 内容概述。概括性地归纳书的主要内容,包括作者陈述的基本观点。3 思路辨析。根据阅历,清理本书思想脉络,并与其他相关著述做比较分析。4 引申性 结合文本现象或个案对总体问题进行引申性把握。5 评价与结语。对...

新数值分析读书报告

第四章数值积分和数值微分。考虑了函数逼近之后,我们现在要考虑的就是对于一些较麻烦的函数怎么求。它的某些运算。例如函数的积分运算,如果函数的性质不好的话,我们就不能求出它的原函数,从而不能用牛顿莱布尼茨公式计算积分函数的值,那么我们就要找一个运算来逼近它,我们介绍的方法是公式 f x dx akf x...