没有合适的资源?快使用搜索试试~ 我知道了~
首页LU-SGS迭代法提升DGM隐式求解效率:降低CPU时间和迭代步数
本文档深入探讨了"基于LU-SGS迭代的DGM隐式方法研究"这一主题,发表于2014年6月的《西北工业大学学报》。LU-SGS迭代法在非结构网格的有限体积法中表现出色,论文作者将其理念引入到高精度间断Galerkin(DGM)有限元的隐式求解策略中。这种改进旨在避免隐式算法中对大型稀疏矩阵求逆带来的挑战,显著提升计算效率。 采用LU-SGS迭代方法,仅需在每个时间步推进过程中进行两次网格扫描,一次从前向后,一次从后向前,大大减少了计算复杂度和内存需求。这种方法特别适合处理大规模数值模拟,如跨声速无粘流场的数值模拟,研究对象包括NACA0012翼型和ONERA M6机翼。 对比TVD-RKDG显式时间格式,隐式格式显示出显著的优势,不仅迭代步数大幅度减少,而且所需的CPU时间也有了明显缩减。尽管引入了隐式处理,但精度保持在同一高水平,证明了这种方法的有效性和实用性。 关键词:间断Galerkin有限元、LU-SGS迭代、欧拉方程、非结构网格、时间隐式方法。这篇论文的学术价值在于它填补了国内在DGM隐式方法领域的研究空白,展示了国内学者在这一领域的逐渐成熟和发展,特别是在处理复杂工程问题上的创新应用。对于从事数值计算和流体力学的工程师和技术人员来说,这篇论文提供了宝贵的经验和参考依据。
资源详情
资源推荐
%&’(
年
W
月
第
3%
卷
第
3
期
西北工业大学学报
!0"D#?1 0E +0DLBS$OL$D# V016L$5B#45?1 >#4C$DO4L6
!"#$
/012 3%
%&’(
+02 3
收
稿日期
: %&’3G’&G%%
基金项目
:
国家自然科学基金
( -&7,W&’,)
与国家
“.W3”
计划
( %&’%HH&’H3&()
资助
作者简介
:
郝海兵
( ’7.’—)
,
中国航空计算技术研究所工程师
、
博
士
,
主要从事高精度数值算法研究
。
基
于
R>G*=*
迭代的
:=I
隐式方法研究
郝
海兵
’
,
张
强
%
,
杨
永
%
,
梁
益华
’
’2
中国航空计算技术研究所
,
陕
西 西安
,’&&W.
%2
西北工业大学 翼型叶栅空气动力学国防科技重点实验室
,
陕西 西安
,’&&,%
(
)
摘 要
:
考
虑到
R>G*=*
迭代法已经在基于非结构网格的有限体积法中得到了成功应用
,
文章借鉴其
思想
,
将其推广到高精度间断
=?1$DT4#
有限元隐式格式求解中来
,
并对其性能进行了研究
。
为了避免
隐式算法中对大型稀疏矩阵求逆
,
采用
R>G*=*
迭代法
,
只需要在每步时间推进中沿网格号从前到后
和从后到前
%
次扫描计算即可
,
并且还能有效降低内存需求
。
通过对
+H8H&&’%
翼型和
9+;XH IW
机翼跨声速无粘流动进行数值模拟
,
计算结果表明
:
与
F/:GXZ:=
显式时间格式相比
,
隐式格式所需
的迭代步数和
8V>
时间均得到了很大程度上的减少
,
并且精度保持不变
。
关 键 词
:
间断
=?1$DT4#
有限元
; R>G*=*;
欧拉方程
;
非结构网格
;
时间隐式方法
中图分类号
: /%’’
文献标志码
: H
文章编号
: ’&&&G%,-.( %&’() &3G&3(WG&-
间断
=?1$DT4#
方法
[’]
( :=I)
由于具备易于实
现
高精度以及求解间断问题等优点而受到广泛关
注
。
该方法在任意单元内部通过提高逼近多项式的
阶次实现高精度
;
空间离散时允许单元边界上存在
间断
,
具有灵活处理复杂区域边值问题的能力
;
在单
元边界处可以采用有限体积法中的思路
,
基于
X4$G
N?##
问题来构造数值通量函数
,
实现逆风格式
,
从
而更有利于求解间断问题
。
目前
,
国外对该方法已
经展开了近十多年的研究
,
并取得了较大的进
展
[%]
,
而
国内相关研究起步较晚
,
在近几年才逐渐
受到 重 视
。
国内主要进展包括邱建贤等
[3]
将
\;+9
格式作为限制
器引入到
:=I
中
;
张来平
等
[(]
提出了静动态混合
重构的
:= ]Q/
混合格式
等
。
当前
:=I
主要应用于空间离散
,
而时间方向
推进多沿用了
*B"
等
[%]
构
造的
F/:GXZ:=
方法
。
因此
,
当求解更高精度或大规模计算问题时
,
稳定性
条件变得越来越苛刻
,
从而导致当地时间步长取值
越来越小
,
收敛率会急剧下降
,
计算非常耗时
。
针对
该问题
,
目前主要采用
K
型多重网格方法
[%]
和
隐式
时间离散
%
种加速收敛技术来提高计算效率
。
本文
主要对隐式格式离散进行研究
,K
型多重网格方法
见参考文献
[%,- ]。
目前
,
国外很多学者对间断
=?1$DT4#
有限元隐式算法进行了多种试探性的研
究
,
不过目前还没有形成成熟稳定的算法
。
主要进
展包括如下
: A ?OO4
和
X$P?6
首次将
=IX;*
隐式算
法应 用 于
:=
求解二维可压缩
+G*
方 程
[W]
;
U?DLN?##
等
人使用
=IX;*G+$SL0#
算法求解二维
;"1$D
和
+G*
方程
[,]
等
。
本
文基于非结构网格
,
针对定常
;"1$D
方程组
,
发展相应的
:=I
隐式算法
。
考虑到
R>G*=*
迭代
法不需要存储和处理大型矩阵
,
并且非常适合于非
结构网格的无序性
,
具有显著的计算效率
,
并且在有
限体积法中得到了广泛应用
。
本文借鉴其思想
,
将
该方法推广到高精度
:=I
隐式格式中
,
构造适合
于
:=I
隐式 求解的
R>G*=*
迭代 格式
。
通过对
+H8H&&’%
翼型和
9+;XH IW
机翼跨声速无粘流
动进行数值模拟验证其加速收敛性能
。
1
控
制方程
考
虑非定常
;"1$D
方程在直角坐标系的守恒
形式
:
∂
U( x,t )
∂
t
+
∇
·F( U( x,t)
)
=
& ( ’)
下载后可阅读完整内容,剩余4页未读,立即下载
weixin_38631225
- 粉丝: 5
- 资源: 908
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 最优条件下三次B样条小波边缘检测算子研究
- 深入解析:wav文件格式结构
- JIRA系统配置指南:代理与SSL设置
- 入门必备:电阻电容识别全解析
- U盘制作启动盘:详细教程解决无光驱装系统难题
- Eclipse快捷键大全:提升开发效率的必备秘籍
- C++ Primer Plus中文版:深入学习C++编程必备
- Eclipse常用快捷键汇总与操作指南
- JavaScript作用域解析与面向对象基础
- 软通动力Java笔试题解析
- 自定义标签配置与使用指南
- Android Intent深度解析:组件通信与广播机制
- 增强MyEclipse代码提示功能设置教程
- x86下VMware环境中Openwrt编译与LuCI集成指南
- S3C2440A嵌入式终端电源管理系统设计探讨
- Intel DTCP-IP技术在数字家庭中的内容保护
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功