计算流体力学课程作业(共15页).docx
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《计算流体力学课程作业(共15页).docx》由会员分享,可在线阅读,更多相关《计算流体力学课程作业(共15页).docx(15页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上计算流体力学大作业有限差分法解Poisson方程五点格式解区域内Poisson方程摘要:本文结合计算流体力学课上所学知识,采用数值解法中的有限差分法求解Poisson方程(偏微分方程中椭圆型方程的一种),并用其五点格式采用高斯塞德尔(Gauss-Seidel)迭代求解。并比较了数值近似解与真实解,以及不同步长情况下误差的大小,得到了一定的结论。关键词:Poisson方程 有限差分法 五点格式一、计算流体流体力学的特点计算流体力学中许多问题求解最终都会变成偏微分方程的求解,而在数学上,除了几种极少数情况外,要求出它们精确解是很难的。计算机技术的发展使得这一难题的一很好地
2、解决。二、偏微分方程的种类2.1、 椭圆型偏微分方程椭圆型偏微分方程的一般形式为 其中:若,为的梯度,则其定义为 散度的定义为 这样,可以更明确地表示为 若为常数,则进一步化简为 其中,又称为Laplace算子。这样椭圆型偏微分方程可以简单地写为 2.2、抛物型偏微分方程抛物型偏微分方程的一般形式为 根据上面叙述,若为常数,则该方程可以更简单地写为 2.3、双曲型偏微分方程双曲型偏微分方程的一般形式为 若为常数,则可以将该方程简化为 三类方程的直接的区别在于对的导数的阶次。若对没有求导,可以理解为其值为常数,故称为椭圆型的。若取对时间的一阶导数,则与对的二阶导数直接构成了抛物线关系,故称为抛物
3、型偏微分方程。若取对时间的二阶导数,称其为双曲型偏微分方程。三、Poisson方程:泊松方程是数学中一个常见于静电学、机械工程和理论物理的偏微分方程。是从法国数学家、几何学家及物理学家泊松而得名的。泊松方程一般可写为:=f在这里 代表的是拉普拉斯算符(也就是哈密顿算符的平方),而 f 和 可以是在流形上的实数或复数值的方程。拉普拉斯方程: 因此泊松方程通常写成: 在三维直角坐标系,可以写成 如果没有f, 这个方程就会变成拉普拉斯方程=0.泊松方程可以用格林函数来求解;如何利用格林函数来解泊松方程可以参考screened Poisson equation。现在有很多种数值解法,数学上,泊松方程属
4、于椭圆型方程。四、Poisson方程的解法泊松首先在无引力源的情况下得到泊松方程,=0(即拉普拉斯方程);当考虑引力场时,有=f(f为引力场的质量分布)。后推广至电场磁场,以及热场分布。该方程通常用格林函数法求解,也可以分离变量法,特征线法求解。由于Poisson方程难以求得其解析解,计算机技术发展之后,数值解法成为工程实际中应用最广泛的求Poisson方程解的方法,常见的数值解法有:里兹(Ritz)法,加权余量法,有限差分法,有限元法,边界元法及有限体积法等。4.1、里兹(Ritz)法瑞利-里兹法(也称里兹法)是通过泛函驻值条件求未知函数的一种近似方法,是英国的瑞利于1877年在声学理论一书
5、中首先采用,后由瑞士的W.里兹于1908年作为一个有效方法提出。这一方法在许多力学、物理学、 量子化学问题中得到应用。同时它也是广泛应用于应用数学和机械工程领域的经典数值方法,它可以用来计算结构的低阶自然频率。它是直接变分法的一种,以最小势能原理为理论基础。通过选择一个试函数来逼近问题的精确解,将试函数代入某个科学问题的泛函中, 然后对泛函求驻值,以确定试函数中的待定参数,从而获得问题的近似解。4.2、加权余量法加权余量法(Weighted residual approach),又称加权残量法,加权残余法。当 n 有限时,定解方程存在偏差(余量)。取权函数,强迫余量在某种平均意义上为零。采用使
6、余量的加权积分为零的等效积分的“弱”形式来求得微分方程近似解的方法称为加权余量法。加权余量法在固体力学中,是求解线性、非线性微分方程的一种有效方法,它是基于等效积分形式的近似方法,也是通用的数值计算方法有限元法、边界元法、无网格法都是加权余量法的特殊情况,由于这三种方法各有其特点,所以都各自发展为一种独立的方法,加权余量法最早是用于流体力学,传热等科学领域,后在固体力学中得到了更大的发展。权函数的选择加权余量法是求解微分方程近似解的一种有效方法显然,任何独立的完全函数都可用来作为权函数,加权余量法可分为内部法、边界法和混合法,在内部法中,又可分为:配置法(以笛拉克函数作为权函数),子域法,最小
7、二乘法,矩量法,伽辽金法等。4.3、有限差分法有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替, 这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似, 积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组 , 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题,即所谓的有限差分
8、法。有限差分法求解偏微分方程的步骤如下:1、区域离散化,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格;2、近似替代,即采用有限差分公式替代每一个格点的导数;3、逼近求解。换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。4.4、有限元法有限元法(finite element method)是一种高效能、常用的计算方法。有限元法在早期是以变分原理为基础发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中(这类场与泛函的极值问题有着紧密的联系)。自从1969年以来,某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最
9、小二乘法等同样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系。基本思想:由解给定的泊松方程化为求解泛函的极值问题。将连续的求解域离散为一组单元的组合体,用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数,近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。有限元法的一般步骤:步骤1:剖分:将待解区域进行分割,离散成有限个元素的集合元素(单元)的形状原则上是任意的二维问题一般采用三角形单元或矩形单元,三维空间可采用四面体或多面体等每个单元的顶
10、点称为节点(或结点)步骤2:单元分析:进行分片插值,即将分割单元中任意点的未知函数用该分割单元中形状函数及离散网格点上的函数值展开,即建立一个线性插值函数步骤3:求解近似变分方程用有限个单元将连续体离散化,通过对有限个单元作分片插值求解各种力学、物理问题的一种数值方法。有限元法把连续体离散成有限个单元:杆系结构的单元是每一个杆件;连续体的单元是各种形状(如三角形、四边形、六面体等)的单元体。每个单元的场函数是只包含有限个待定节点参量的简单场函数,这些单元场函数的集合就能近似代表整个连续体的场函数。根据能量方程或加权残量方程可建立有限个待定参量的代数方程组,求解此离散方程组就得到有限元法的数值解
11、。有限元法已被用于求解线性和非线性问题,并建立了各种有限元模型,如协调、不协调、混合、杂交、拟协调元等。有限元法十分有效、通用性强、应用广泛,已有许多大型或专用程序系统供工程设计使用。结合计算机辅助设计技术,有限元法也被用于计算机辅助制造中。4.5、边界元法边界元法(boundary element method)是一种继有限元法之后发展起来的一种新数值方法,与有限元法在连续体域内划分单元的基本思想不同,边界元法是只在定义域的边界上划分单元,用满足控制方程的函数去逼近边界条件。所以边界元法与有限元相比,具有单元个数少,数据准备简单等优点。但用边界元法解非线性问题时,遇到同非线性项相对应的区域积
12、分,这种积分在奇异点附近有强烈的奇异性,使求解遇到困难。边界元法以定义在边界上的边界积分方程为控制方程,通过对边界分元插值离散,化为代数方程组求解。它与基于偏微分方程的区域解法相比,由于降低了问题的维数,而显著降低了自由度数,边界的离散也比区域的离散方便得多,可用较简单的单元准确地模拟边界形状,最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数 ,而具有解析与数值相结合的特点,通常具有较高的精度。特别是对于边界变量变化梯度较大的问题 ,如应力集中问题 ,或边界变量出现奇异性的裂纹问题,边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算 流体力学 课程 作业 15
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内