收稿日期
基金项目 国家自然科学基金资助项目
作者简介 郭 隽 女 辽宁盘锦人 博士后 北京
非 线 性 问 题 的 MPS 无 网 格 算 法
郭 隽 陶 智
北京航空航天大学 动力系
摘 要 给出了修正的 MPSMoving Particle Semiimplicit无网格方法 指出
恰当选取权函数可避免使用虚节点来处理边界 从而解决复杂边界处虚节点难设定
的问题 以非线性 Burgers 方程为算例分析了其数值特性 结果表明 以样条函数为权
函数时不采用虚节点的结果与使用虚节点的结果一致 对一维问题 光滑权函数性能
明显优于非光滑权函数 而对二维问题 近边界处加权域不对称性的影响则不很显
著 与精确解或由传统方法所得结果的对比分析表明该方案可行 该方法物理概念清
晰 易于编程实现
关 键 词 数值模拟 非线性方程 加权 无网格方法 MPS 方法
中图分类号 O
文献标识码 A 文 章 编 号
近年来无网格方法由于只需节点信息 无需
生成单元而倍受关注
尽管其中一些方法仍
未能完全避免网格的介入 但它们确实表现出许
多优良特性 如易施行自适应技术及多尺度分析
可处理断裂和接触等常规方法难处理的问题
MPS Moving Particle Semiimplicit方法
就
是一种完全无需网格介入的方法 具有物理概念
清晰 易于编程实现等特点 实现了不可压问题的
无网格求解 但该方法在边界外需引入虚节点 限
制了复杂几何条件下的应用 为此本文对其作了
一定修正 通过优化选取加权函数 从而不再使用
虚节点 指出选用样条函数为权函数时可不引入
虚节点 解决了复杂边界虚节点难设定的问题
1 基本方程及求解方案
11 基本方程
Burgers 方程常用以检验求解 NS 方程的数值
算法的精度与稳定性 二维 Burgers 方程为
u
t
u
u
c
u
u
v
t
u
u
c
v
v
式中 u 为流体速度 u
c
为计算点速度 u
c
u
时为拉氏计算 u
c
时为欧拉计算
1 2 数值求解方法
采用显式分步法求解上述方程 分为两步
显式计算扩散项 获得速度及位置的中间
值 u
L
v
L
和 r
L
u
L
u
n
t
u
x
u
y
n
v
L
v
n
t
v
x
v
y
n
r
L
r
n
u
L
d t
显式计算对流项 获得各物理量的新值
r
n
r
n
u
n
u
L
tu
a
u
L
v
n
v
L
tu
a
v
L
式中 u
a
u
u
c
r
n
r
L
t
2 MPS 无网格算法
21 邻点交互作用模型
依据插值理论 将各点参数值与其相邻点间
通过加权函数相关联 利用配点法对微分方程进
行等效变换 进而实现数值求解 设求解区域中各
点 r
i
的物理参数值为 f
i
则任一点 r
i
的光滑值可
用其邻点的加权平均
f r
i
j
f
j
w
|
r
j
r
i
|
j
w
|
r
j
r
i
|
来表示 相当于 MLS Moving Least Square方法中
年 月
第 卷 第期
北 京 航 空 航 天 大 学 学 报
Journal of Beijing University of Aeronautics and Astronautics
January
Vol No