大家都来说说upwind的缺点

a9002

迎风产生的数值粘性太大了，当然流场稳定性相对也就比较好了，所以商业化软件喜欢采用它，但要看到细致的涡阶好像迎风格式搞不定啊。

whgoodboy

什么啊。不要把迎风说得这么牛，我在做sonic point glitch时，迎风有着不可克服的困难。

newcfd

[这个贴子最后由newcfd在 2005/05/26 04:42am 第 2 次编辑]

Probelms:
  1. slow convergence
  2. wrong solutions by flows with low Mach numbers (Ma < 0.01),   
  3. have hard time to get a solution by supersonic flows
    ( Ma > 2.0 ) with strong shocks on unstructured mesh.
Solutions:
  1. implicit (normally a must )
  2. preconditioning
  3. make a nice mesh across shocks and tune limiters  
Advantages:
  1. high resolution
  2. use fewer mesh points in boundary layer
  3. have more physical meanings
  

whgoodboy

一个很重要的缺陷：sonic point glitch.几乎所有的迎风型格式在音速点都会有问题

xtuyyy

[这个贴子最后由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方程:
&#36;u=1,if x<=0; u=0,if x>0,-1<=x<=1,0<=t<=1.&#36;假设边界条件为常数，假设等空间步长等时间步长.
&#36;\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&#36;
按照这个离散格式，激波不会前进。
请指教。

laofang

传统迎风TVD都存在以下弱点既在间断点或者极值点降阶，出现所谓的极值“修剪”现象.ENO和WENO的提出解决了这个问题。但是文献中来验证格式精度的算例基本是纯量疏运方程或者长波情况。对于短波（初始解变化率较大）、非线性方程，并不见得取得多好的效果。我的经验是波越短，极值衰减越严重。
以上是我个人学习所见，那位高手如有兴趣，请指教或者讨论。
qq：154309922

hztang

For whgoodboy:
H.Z. Tang, On the Sonic Point Glitch, J. Comput. Phys., 202 (2): 507-532 JAN 20 2005  

ymcnuaa

陶老师的《数值传热学》上面讲得很详细的

whgoodboy

for hztang:
  您是汤老师吧？我九所的

ypbing

一阶的逆风格式数值耗散太大，解纯对流问题隐含了数值粘性项，比较烦人,尤其在处理间断问题的时候，耗散将特征给均匀化了。

lu2000

下面引用由abcc在 2004/11/30 05:19pm 发表的内容：
xtuyyy u=c>0,if x<=0; u=c2<=0,if x>0.）无法计算。
这怎么会不可以算?
你把迎风看的太窄了, 实际上你这个问题你说不能算可能是因为你觉得它的特征在x=0附近符号变化, 这在迎风格式中都是通过　...