参考学习资料 流体输配管网 2023_Sci_Sin_Math_Survey.pdfVIP

中国科学:数学2023年第49卷第?期:125

SCIENTIASINICAMathematica

综述

广义对流扩散问题的稳定化有限元方法

吴朔男1∗

1.北京大学数学科学学院,北京100871

E-mail:snwu@

收稿日期:2023-06-30;接受日期:2023-08-30;网络出版日期:202X-XX-XX

国家自然科学基金(批准号:、北京市自然科学基金(批准号:1232007)资助项目

摘要本文介绍了稳态对流扩散问题的稳定化有限元方法.其主要难点在于当对流占优时可能出现边

界层,导致传统有限元方法在边界层内失去稳定性,从而产生剧烈震荡.在拟均匀网格下,稳定化有限

元方法可分为两类:迎风型方法和指数拟合方法.前者利用对流速度的信息在变分形式中增加稳定化

项，而后者利用边界层解的特征将指数函数引入到格式设计中.这两类方法对于设计电磁场等新型对

流扩散问题的数值方法起到重要指导作用.

关键词对流扩散方程有限元方法迎风型指数拟合电磁场对流扩散

主题分类65N12,65N15,65N30

1引言

对流扩散方程是一类常见的偏微分方程,描述了物质在流体介质中同时受到对流和扩散的影响.

其中,对流项描述了流体运动引起的物质输运,扩散项描述了物质由于浓度梯度而发生的自发扩散.对

流扩散方程在多个科学和工程领域中都具有广泛的应用,包括流体力学、天体物理学、地下水流、气

象学、半导体和反应流动等.

在有界区域Ω⊂Rd上,经典的标量对流扩散方程有如下形式:

Lu:=−ε∆u+·∇u+γu=f,在Ω内.(1.1)

这里,u表示待求的物理量(如浓度、温度等),ε0表示扩散系数,表示速度场.该问题的适定性可

以在一定的边界条件下(如u|∂Ω=0)成立.扩散项−ε∆u与对流项·∇u的相互作用导致了方程解

的特殊行为.例如,当问题是对流占优,即ε≪||时,问题(1.1)的解可能出现内部或边界层结构(在

非常小的区域解及其导数急剧变化).边界Γ:=∂Ω根据流场方向,可分为入流出流两类,即:

Γ−:={x∈Γ:·0},Γ+:={x∈Γ:·0}.

英文引用格式WuSN.Stabilizedfiniteelementdiscretizationsofgeneralconvection-diffusionproblems(inChinese).SciSin

Math,2023,49:1–25,doi:10.1360/N012017-XXXX

⃝《中国科学》杂志社

吴朔男:广义对流扩散问题的稳定化有限元方法

对于纯对流问题(即ε=0),仅需在入流边界Γ−处施加Dirichlet边界条件.因此,当ε→0+时,(1.1)

的边界条件不自洽性是导致解可能存在边界层结构的原因之一.考虑u|∂Ω=0的简单情形,其变分问

题为:求u∈H1(Ω)使得

0

a(u,v):=ε∫∇u·∇v+∫(·∇u)v+∫γuv=∫fv∀v∈H1(Ω).(1.2)

0

ΩΩΩΩ

在离散层面上,记Th为Ω的正则单纯形剖分族(简单起见,这里假设区域为多边形).对于任意

T∈Th,hT表示T的直径.对于变分问题(1.2),如果使用经典的有限元方法,并且边界层处网格尺度

hT较大,那么得到的数值解通常会在整个区域中产生大幅度的数值震荡.