|
发表于 2005-7-4 11:52:33
|
显示全部楼层
大家都来说说upwind的缺点
[这个贴子最后由xtuyyy在 2005/07/04 12:23pm 第 1 次编辑]
下面引用由abcc在 2004/11/30 05:19pm 发表的内容:
xtuyyy u=c>0,if x<=0; u=c2<=0,if x>0.)无法计算。
这怎么会不可以算?
你把迎风看的太窄了, 实际上你这个问题你说不能算可能是因为你觉得它的特征在x=0附近符号变化, 这在迎风格式中都是通过 ... 忘了说明是"一阶单调迎风"差分格式.
考虑无粘Burgers方程:
$u=1,if x<=0; u=0,if x>0,-1<=x<=1,0<=t<=1.$假设边界条件为常数,假设等空间步长等时间步长.
$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delt t} + \frac{u+|u|}{2} \frac{u_i^n-u_{i-1}^n}{\Delt x} + \frac{u-|u|}{2} \frac{u_{i+1}^n-u_{i}^n}{\Delt x} = 0$
按照这个离散格式,激波不会前进。
请指教。 |
|