Główne pojęcia
本文提出了一种基于多尺度有限元法和隐显式方案的算法,用于求解具有多尺度特性的非平稳Stokes-Darcy模型。该算法首先在Darcy区域并行计算多尺度基函数,然后基于这些多尺度基函数采用隐显式方案求解Stokes-Darcy方程。该算法可以在相对较粗的网格上求解问题,大幅降低了计算成本,同时相比于标准有限元方法在相同网格尺度下具有更高的精度。
Streszczenie
本文研究了具有多尺度特性的非平稳Stokes-Darcy模型的数值求解方法。主要内容如下:
-
引入了Stokes-Darcy耦合模型,其中Darcy区域的渗透系数呈现多尺度特性。
-
提出了一种基于多尺度有限元法和隐显式方案的算法(MsFEM-ImEx)来求解该模型。该算法分两步进行:
- 首先在Darcy区域并行计算多尺度基函数。
- 然后基于这些多尺度基函数,采用隐显式方案求解Stokes-Darcy方程。
-
该算法可以在相对较粗的网格上求解问题,大幅降低了计算成本,同时相比于标准有限元方法在相同网格尺度下具有更高的精度。
-
在假设渗透系数是周期性函数且与时间无关的条件下,证明了该算法的稳定性和收敛性。
-
通过三个数值实验验证了该算法的合理性和有效性,数值结果与理论分析一致。
Statystyki
∥uN
h ∥2
0 + ∥ϕN
h ∥2
0 +
N−1
X
n=0
(∥un+1
h
− un
h∥2
0 + ∥ϕn+1
h
− ϕn
h∥2
0) + C1ν∆t
2
∥uN
h ∥2
1 + λmin∆t
2S0
∥ϕN
h ∥2
1
≤ C(T)( C2
p
3C1ν ∆t
N−1
X
n=0
∥fn+1
f
∥2
0 +
‹
C2
p
S0λmin
∆t
N−1
X
n=0
∥f n+1
p
∥2
0
∥u0∥2
0 + ∥ϕ0∥2
0 + C1ν∆t
2
∥u0∥2
1 + λmin∆t
2S0
∥ϕ0∥2