测试双曲方程的 SUPG 方法和其他方法

计算科学 有限元 稳定 双曲-pde
2021-12-14 01:46:01

我对整合简单方程很感兴趣

ϕt+uϕ=0
在流入边界处具有狄利克雷边界条件(un<0)。域是一个简单的单位平方。我假装使用 SUPG,我想使用制造解决方案的方法来验证该方法。我想出了一个简单的解决方案

ϕ={cos(5xt) for y>0.50 else

平流场u=(1.0/5.0,0). 这使得方程的源项等于零,我适当地调整了 Dirichlet BC

仅使用拉格朗日元素离散该方程,不添加任何稳定方法,并与 Crank-Nicholson 及时积分会产生没有振荡的良好解,这是我期望在具有不连续函数的双曲方程上使用拉格朗日元素时所期望的。为什么会这样?是否有其他制造的解决方案会真正触发振荡?

0个回答
没有发现任何回复~