- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
JFNK方法迭代过程与物理约束.pdf
第29卷第5期 计 算 物 理 V01.29,No.5
CHINESE OF
2012年9月 JOURNALCOMPUTATIONALPHYSICS Sep.,2012
文章编号:1001二246X(2012)05-0654-07
JFNK方法迭代过程与物理约束
安恒斌, 莫则尧
(北京应用物理与计算数学研究所计算物理实验室,北京100094)
摘要:对Jacobian—free
JFNK方法进行改进.改进后的JFNK方法迭代过程中的迭代序列总是满足物理约束,克服了迭代过程中可能-出现
的非物理现象.采用改进之后的算法求解二维三温能量方程,可以保证在迭代过程中不会出现负温度问题,使
JFNK方法的健壮性得到提高.
关键词:JFNK;Newton—Krylov方法;迭代方法;物理约束;二维三温能量方程
中图分类号:0241.7 文献标识码:A
U 引吾
采用各类非线性迭代方法求解实际问题时,常常需要考虑迭代过程中的物理约束.例如,迭代求解能量
方程时,迭代过程中需要保持能量(温度)非负.采用Jacobian—free
方程时,在迭代过程中需要评估非线性代数方程的残差.由三温能量方程离散之后得到的代数方程的表达式
中涉及到温度的开方运算,因此在迭代过程中评估非线性残差时,如果温度出现负值,则迭代会发生中断.
为保证偏微分方程离散得到的代数方程组具有满足物理约束的解,需要从三个层面加以考虑.首先,物
理模型(即偏微分方程自身加上合适的边界条件)要保证偏微分方程具有满足物理约束的解.其次,离散格
式需要保证由偏微分方程离散得到的代数方程组具有满足物理约束的解.针对二维变形网格上扩散方程的
九点离散格式¨‘51以及几何奇异点的特殊离散方法M1,均需要具有保正性.再次,对于离散得到的代数方程
组,需要采用合适的迭代方法获得满足物理约束的解.只有以上三方面相互配合,才能共同保证得到满足物
理约束的解.本文中,我们关心的是第三个层面的问题.
Jacobi—free
法¨”。.本文对于JFNK方法的迭代过程进行分析,特别是关注迭代过程中可能不满足物理约束的关键步
骤.通过吸收物理量约束的信息,对迭代过程进行改进,使迭代过程总是保持物理约束.
非线性方程组的一般形式
F(x)=0, (1)
其中F:尺“_R“为非线性映射.从数学上看,物理约束其实就是对于映射F的定义域进行限制.设映射F的
定义域为够cR“,由实际模型离散得到的代数方程组
F(X)=0,F:够c (2)
R“一尺“.
这里,够中的每个点均满足物理约束.因此,在迭代过程中,迭代点以及每次评估非线性残差时对应的自变
量均落在够内.本文中,我们假定定义域够为R8空间中的超立方体
够=n(z‘,tti).
收稿日期:2011—10—26;修回日期:2012一01—20
(2010AA012301,2010AA012303)资助项目
作者简介:安恒斌(1974一),男,宁夏彭阳,副研究员,博士,主要从事大规模非线性方程组迭代解法及并行算法研究,
E—mail:811一hengbin@iapcm.ac.ca
第5期 安恒斌等:JFNK方法迭代过程与物理约束 655
此时,工∈够也可记为LzU,其中L=(1i’,U=(11,‘)∈R“.
定义泛函
0 F(工)02:,,:够cR8-+R
,(工)=丁1
.厂(z)在JFNK算法的全局收敛性设计和分析中具有重要作用‘7“
您可能关注的文档
最近下载
- 2023-2024学年上海市长宁区高中数学人教A版选修二第四章-数列章节测试-19-含解析.pdf VIP
- 大学英语六级词汇汇总.pdf VIP
- 乡愁英文版.pptx
- 2023-2024学年江苏省南通市高中数学人教A版选修二第四章-数列章节测试-19-含解析.pdf
- 2022年黑龙江省大庆市中考道德与法治试卷(word精校版).docx
- JTT329-2010 公路桥梁预应力钢绞线用锚具、夹具和连接器.docx VIP
- 水电站大坝主体混凝土施工组织设计.doc
- 2023-2024学年北京市西城区高中数学人教A版选修二第四章-数列章节测试-9-含解析.pdf
- 20秋新疆电大一体化 社会保障学 (50627)章节测验1-24答案.pdf
- 2023年全国乙卷理科综合高考试卷(原卷+答案).pdf
文档评论(0)