苍山负雪 发表于 2003-12-3 22:01:32

【讨论】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 发表于 2003-12-4 08:55:04

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

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

BosFEA 发表于 2003-12-4 09:03:46

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

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

苍山负雪 发表于 2003-12-4 10:38:12

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

是这样设定的。 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 发表于 2003-12-4 11:49:43

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

Verification中好像还有几个例题,Primer手册中我有最新的,的确是没有这方面的例子,很愁人啊

如果说有松弛因子,出来CFL数外,ADINA的建议是使用时间相关的材料,定义如粘度系数随时间变化,由较大值逐步减小到真实值,我没有用过,不确认对High speed compressible是否可以。

苍山负雪 发表于 2003-12-4 12:22:38

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

有种想法不知道行不行

在开始的时间步时把收敛误差调高一些,让计算可以进行下去,等边界条件已经达到给定值的时候,run重启动分析,把收敛误差降低。这个时候应该可以得到比较合理的结果了吧。

帮助文件里说high-speed compressible flow 通过Roe格式通量分裂来对付对流项,应该说Roe的激波分辨率是不错的,就是计算量大收敛得慢了点儿。不知道具体的算法怎么处理的。

BosFEA 发表于 2003-12-4 13:10:37

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

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

aerospace 发表于 2003-12-4 15:50:36

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

分裂项一般采用显式算法,对否?
请大侠指点!

苍山负雪 发表于 2003-12-4 18:33:54

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

晕死了,运行了一下午都没有出来一个结果,自动时间步一直减少时间步,就是不收敛。而且算起来好像很慢的样子,15000单元而已。看了一下中间结果 :( 基本上不堪入目,唉对adina没有信心了:(。

这是一个练习题目,稳态的,用来验证算法的激波捕捉能力的。要是真格的做个算例恐怕mission impossible?

算法的显式或隐式与分裂项没有什么关系吧,你看adina-f对流项通量分裂,算法里面不是一样有线式和隐式可以使用

aerospace 发表于 2003-12-5 18:56:54

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

辛苦辛苦!!
页: [1]
查看完整版本: 【讨论】fluid求解第一步总是通不过