【讨论】fluid求解第一步总是通不过

苍山负雪

adina-f :(:(
  
:(*** ERROR *** CODE ADF2254:
Failed in solving boundary equations in boundary condition 11.
Possible causes are:
* Incorrect boundary conditions have been applied;
* Physical model is unstable;
* Solution procedure is divergent.
  
前台阶模型
进口超音速，出口自由。automatic timestep on, cfl数调到1，gmres迭代求解。
  
求解的第一步就出错了, 就是上面的那个error.  
查了一下帮助，对于supersonic flow 应该考虑1、边界条件满足状态方程；2、初值合理（又查帮助文件说要求初值设定也要满足状态方程也就是ideal gas law）；3、工质属性给定正确。
  
我设定工质属性应该没问题的，就是按照理想气体的参数给定。初值按照超音速入口边界条件给定。
  
这方面求解有什么技巧么？

BosFEA

载荷是否为逐渐变化到最终载荷？

BosFEA

如果入口采用P,v,T组合的条件，容易让载荷变化

苍山负雪

是这样设定的。 P,T 不变， V从0增大到最终值。 通过20个时间步，每个时间步间隔为1。
  
另外初值给定的应该没什么问题吧
  
不对啊不对，v从0增大到最终值，那初值设定的时候也应该是v=0才对啊。
  
改了一下初值v=0，可以计算了，但是迭代残差降不下去，最低能到0.002左右，有什么办法能降到1e-4？ 求解器有没有类似松弛因子的设定？
  
我把收敛误差调高到0.002, 也能得到结果，不过比较惨不忍睹
  
另外，很奇怪的是verfication和primer里面都没有任何关于high-speed compressible flow 的例子啊？

BosFEA

Verification中好像还有几个例题，Primer手册中我有最新的，的确是没有这方面的例子，很愁人啊
  
如果说有松弛因子，出来CFL数外，ADINA的建议是使用时间相关的材料，定义如粘度系数随时间变化，由较大值逐步减小到真实值，我没有用过，不确认对High speed compressible是否可以。

苍山负雪

有种想法不知道行不行
  
在开始的时间步时把收敛误差调高一些，让计算可以进行下去，等边界条件已经达到给定值的时候，run重启动分析，把收敛误差降低。这个时候应该可以得到比较合理的结果了吧。
  
帮助文件里说high-speed compressible flow 通过Roe格式通量分裂来对付对流项，应该说Roe的激波分辨率是不错的，就是计算量大收敛得慢了点儿。不知道具体的算法怎么处理的。

BosFEA

重启动的方法值得试一下，我感觉应该有效果；
你这个问题是稳态的吧？

aerospace

分裂项一般采用显式算法，对否？
请大侠指点！

苍山负雪

晕死了，运行了一下午都没有出来一个结果，自动时间步一直减少时间步，就是不收敛。而且算起来好像很慢的样子，15000单元而已。看了一下中间结果 基本上不堪入目，唉对adina没有信心了:(。
  
这是一个练习题目，稳态的，用来验证算法的激波捕捉能力的。要是真格的做个算例恐怕mission impossible?
  
算法的显式或隐式与分裂项没有什么关系吧，你看adina-f对流项通量分裂,算法里面不是一样有线式和隐式可以使用

aerospace

辛苦辛苦！！