数值计算的六大方法

有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。

对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。

插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。

单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。

在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。

对于有限元方法,其基本思路和解题步骤可归纳为

02多重网格方法

多重网格方法通过在疏密不同的网格层上进行迭代,以平滑不同频率的误差分量.具有收敛速度快,精度高等优点.

多重网格法基本原理微分方程的误差分量可以分为两大类,一类是频率变化较缓慢的低频分量;另一类是频率高,摆动快的高频分量。一般的迭代方法可以迅速地将摆动误差衰减,但对那些低频分量,迭代法的效果不是很显著。

高频分量和低频分量是相对的,与网格尺度有关,在细网格上被视为低频的分量,在粗网格上可能为高频分量。

多重网格方法作为一种快速计算方法,迭代求解由偏微分方程组离散以后组成的代数方程组,其基本原理在于一定的网格最容易消除波长与网格步长相对应的误差分量。

该方法采用不同尺度的网格,不同疏密的网格消除不同波长的误差分量,首先在细网格上采用迭代法,当收敛速度变缓慢时暗示误差已经光滑,则转移到较粗的网格上消除与该层网格上相对应的较易消除的那些误差分量,这样逐层进行下去直到消除各种误差分量,再逐层返回到细网格上。

目前两层网格方法从理论上已证明是收敛的,并且其收敛速度与网格尺度无关。

多重网格法是迭代法与粗网格修正的组合,经过证明迭代法可迅速地将那些高频分量去掉,粗网格修正则可以帮助消除那些光滑了的低频分量,而对那些高频分量基本不起作用。

在多重网格计算中,需要一些媒介把细网格上的信息传递到粗网格上去,同时还需要一些媒介把粗网格上的信息传递到细网格上去。限制算子Iih(i-1)h是把细网格i-1层上的残余限制到粗网格i层上的算子,最简单的算子是平凡单射,另外还有特殊加权限制;插值算子Iih(i-1)h是把粗网格i层上的结果插值到细网格i-1层上的算子,一般采用线性插值或完全加权限制算子。

03有限差分方法

有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。

目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

04有限体积法

有限体积法(FiniteVolumeMethod)又称为控制体积法。其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。简言之,子区域法属于有限体积发的基本方法。

有限体积法的基本思路易于理解,并能得出直接的物理解释。离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒;而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程中不同的项采取不同的插值函数。

05近似求解的误差估计方法

近似求解的误差估计方法共有三大类:单元余量法,通量投射法及外推法。

单元余量法广泛地用于以FEM离散的误差估计之中,它主要是估计精确算子的余量,而不是整套控制方程的全局误差。这样就必须假定周围的单元误差并不相互耦合,误差计算采用逐节点算法进行。单元余量法的各种不同做法主要来自对单元误差方程的边界条件的不同处理办法。基于此,该方法能够有效处理局部的残余量,并能成功地用于网格优化程序。

06多尺度计算方法

近年来发展的多尺度计算方法包括均匀化方法、非均匀化多尺度方法、以及小波数值均匀化方法、多尺度有限体积法、多尺度有限元法等。

均匀化方法是一种多尺度分析的方法。该方法通过对单胞问题的求解,把细观尺度的信息映射到宏观尺度上,从而推导出宏观尺度上的均匀化等式,即可在宏观尺度上求解原问题。

均匀化方法在很多科学和工程应用中取得了巨大成功,但这种方法建立在系数细观结构周期性假设的基础上,因此应用范围受到了很大限制。

鄂维南等提出的非均匀化多尺度方法,是构造多尺度计算方法的一般框架。该方法有两个重要的组成部分:基于宏观变量的整体宏观格式和由微观模型来估计缺少的宏观数据,多尺度问题的解通过这两部分共同得到。

多尺度有限元方法是由Babuska等提出的。该法在宏观尺度上进行网格剖分,然后通过在每个单元里求解细观尺度的方程(构造线性或者振荡的边界条件)来获得基函数。从而把细观尺度的信息反应到有限元法的基函数里,使宏观尺度的解包含了细观尺度的信息。但多尺度有限元方法在构造基函数时需要较大的计算量。

THE END
1.求加数的方法求加数的方法 我记得有一次,朋友小明来问我数学题。他拿着作业本,一脸苦恼地说:“我在做加法运算,知道和与其中一个加数,可就是想不起来怎么求另一个加数。” 我看着他笑了笑,开始给他讲我的理解。我心里想着,这可是数学里很基础的概念呢,得让他完全明白。 我对小明说:“你看啊,加法其实就是把几个数https://wenku.baidu.com/view/73f68ecf8d9951e79b89680203d8ce2f0166656f.html
2.在线解决多元多次方程组,轻松提高数学效率(在线解多元多次方程组在数学领域,多元多次方程组是一类常见的问题,它涉及到多个变量和多个方程的复杂计算。对于学生和研究人员来说,解决这类问题往往需要耗费大量的时间和精力。幸运的是,随着科技的进步,我们现在可以利用在线工具来简化这一过程。 传统的解方程组方法需要通过手工计算,这不仅容易出错,而且效率低下。而https://ai.zaixianjisuan.com/daishu/article-zai-xian-jie-jue-duo-yuan-duo-ci-fang-cheng-zu-qing-song-ti-gao-shu-xue-xiao-lv.html
3.在线求签求事业的方法在线求签求事业的方法 求签正确步骤,求签有什么讲究 求签是一种古老的传统习俗,人们通过求签来寻求神灵的指引和启示。然而,求签并不是一种简单的行为,其中包含着一些正确的步骤和讲究。在进行求签之前,我们需要了解这些步骤和讲究,以确保求签的准确性和有效性。https://www.mtrjc.com/wen/264524.html
4.HTTP请求方法对照表在线HTTP请求方法查看 HTTP Request Method共计15种 序号方法描述 1GET请求指定的页面信息,并返回实体主体。 2HEAD类似于get请求,只不过返回的响应中没有具体的内容,用于获取报头 3POST向指定资源提交数据进行处理请求(例如提交表单或者上传文件)。数据被包含在请求体中。POST请求可能会导致新的资源的建立和/或已有资http://www.kjson.com/files/httpmethod/
5.千万别开启苹果手机的旁白功能,我现在已经后悔,在线求关闭方法千万别开启苹果手机的旁白功能,我现在已经后悔,在线求关闭方法#数码##数码那点事儿# Video Player is loading. 00:00/00:00 Loaded: 0% 视频加载失败,请查看其他精彩视频 特别声明:以上文章内容仅代表作者本人观点,不代表新浪网观点或立场。如有关于作品内容、版权或其它问题请于作品发表后的30日内与新浪网http://k.sina.com.cn/article_6248139533_m1746b0b0d033018v40.html
6.千万别开启苹果手机的旁白功能,我现在已经后悔,在线求关闭方法82 评论 千万别开启苹果手机的旁白功能,我现在已经后悔,在线求关闭方法 发现更多热门视频 打开APP https://k.sina.cn/article_6248139533_m1746b0b0d033018v40.html
7.爆笑两姐妹:在线求爆炒土豆的n种方法套路!高清在线观看PP视频为您提供爆笑两姐妹:在线求爆炒土豆的n种方法套路!高清视频在线观看,爆笑两姐妹:在线求爆炒土豆的n种方法套路!主要内容:https://v.pptv.com/show/tChS1ziblFVO2NJw.html
8.乔欣脚被蚊子咬后肿成“猪蹄”在线发文求消肿方法近日,乔欣在绿洲晒出被蚊子咬了以后肿成“猪蹄”的脚。https://ent.huanqiu.com/article/9CaKrnKrTTq
9.gitclone和TortoiseGit的Pull报错,在线求解决方法,折腾很久了请问各位大神,我的配置好了linux 的git库,windows克隆后就报以下错误,用TortoiseGit 的Pull 也报错,在线求等解决方法,谢谢!!! git.exe pull -v --no-rebase --progress "origin" Could not chdir to home directory /home/git: No such查看原文Win10 Git 推/拉仓库代码无权限问题 问题如下: git.https://www.pianshen.com/article/15731312635/
10.impala集成kudu问题,求解决方法。急,在线等!【现象】[root@xardc2 ~]# impala-shell Starting Impala Shell without Kerberos authentication Connectedhttps://bbs.csdn.net/topics/392281536
11.家里LED灯常见三个故障及解决方法,简单维修不求人家用LED灯常见的3个故障及解决方法,家装电工必备知识https://www.56.com/u25/v_MTU0NTQ3Mjg2.html
12.民间求孕快又灵的方法有哪些?民间有很多求孕方法,其中快又灵的有枕头下放红枣桂圆,寓意早生贵子;放金铲子在床下,寓意顺利产子;还有在门前种石榴树,寓意多子多福。另外,念经文也可以求子,一般建议念《观世音普门品》或者是《药师琉璃光如来本愿功德经》,只要心诚就会显灵。不过这些方法都只是民间流传的罢了,没有科学证实,所以可信度一般并不是https://shiguan.myzx.cn/21509.html
13.求简单化学方法鉴别1,1与求简单化学方法鉴别1,1-二甲基环丙烷和环己烷,在线等!相关的问答 「(3,5-二甲基苯基)硼酸」有哪些品牌? 「3,5-二甲基苯基硼酸」172975-69-8 「(3,5-二甲基苯基)硼酸」172975-69-8 「3,5-二甲基苯基乙酸」42288-46-0 有优质的「3,5-二甲基苯基硫脲」97480-60-9供货厂家吗?希望提供详细报价。 有https://www.chem960.com/ask/q-4ed772f63c4f4ed8bced858c9b0b658c
14.利用函数方法求参数的取值范围.doc全文免费在线预览全文 利用函数方法求参数的取值范围用函数方法求参数的取值范围既是高中数学中常见的问题,也是历年高考的热点和重点,这类问题中一般都含有多个变量,需要一定的解题技巧,关键是转化化归,通过等价转化把问题顺利解决.下面介绍四种方法供大家参考. 一、分离参数当所求的参数能够与其它变量完全分离出来,并且分离后https://max.book118.com/html/2015/0327/13863912.shtm
15.在线等没头了没头了没头了没头了没头了!!!在线跪求解決方法求救在线等穿不穿装也是这样图片:在线跪求解決方法[/ 求神人解救!!!https://bbs.3dmgame.com/thread-4346839-1-919.html
16.《重塑证券交易心理把握市场脉搏的方法和技术》pdf下载在线阅读《重塑证券交易心理》(BrettN.Steenbarger)电子书网盘下载免费在线阅读资源链接:链接:密码:fxhj 书名https://www.1633.com/ask/289645.html
17.智能语音交互在线语音合成接口智能语音交互API文档必须订购在线语音合成服务,才能调用该接口。 该接口同步返回结果,等待引擎结果最长超时时间为1分钟。即该接口最长1分钟会返回结果。 每次合成的文本字符数最大为300字符;1个中文、1个英文、1个标点符号各算1个字符。 二. URI URI构造方法如下 请求方法 URI POST https://api-wuxi-1.cmecloud.cn:8443/api/lihttps://ecloud.10086.cn/op-help-center/doc/article/36974