【讨论】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、工质属性给定正确。
我设定工质属性应该没问题的,就是按照理想气体的参数给定。初值按照超音速入口边界条件给定。
这方面求解有什么技巧么?
回复: 【讨论】fluid求解第一步总是通不过
载荷是否为逐渐变化到最终载荷?回复: 【讨论】fluid求解第一步总是通不过
如果入口采用P,v,T组合的条件,容易让载荷变化回复: 【讨论】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 的例子啊?
回复: 【讨论】fluid求解第一步总是通不过
Verification中好像还有几个例题,Primer手册中我有最新的,的确是没有这方面的例子,很愁人啊如果说有松弛因子,出来CFL数外,ADINA的建议是使用时间相关的材料,定义如粘度系数随时间变化,由较大值逐步减小到真实值,我没有用过,不确认对High speed compressible是否可以。
回复: 【讨论】fluid求解第一步总是通不过
有种想法不知道行不行在开始的时间步时把收敛误差调高一些,让计算可以进行下去,等边界条件已经达到给定值的时候,run重启动分析,把收敛误差降低。这个时候应该可以得到比较合理的结果了吧。
帮助文件里说high-speed compressible flow 通过Roe格式通量分裂来对付对流项,应该说Roe的激波分辨率是不错的,就是计算量大收敛得慢了点儿。不知道具体的算法怎么处理的。
回复: 【讨论】fluid求解第一步总是通不过
重启动的方法值得试一下,我感觉应该有效果;你这个问题是稳态的吧?
回复: 【讨论】fluid求解第一步总是通不过
分裂项一般采用显式算法,对否?请大侠指点!
回复: 【讨论】fluid求解第一步总是通不过
晕死了,运行了一下午都没有出来一个结果,自动时间步一直减少时间步,就是不收敛。而且算起来好像很慢的样子,15000单元而已。看了一下中间结果 :( 基本上不堪入目,唉对adina没有信心了:(。这是一个练习题目,稳态的,用来验证算法的激波捕捉能力的。要是真格的做个算例恐怕mission impossible?
算法的显式或隐式与分裂项没有什么关系吧,你看adina-f对流项通量分裂,算法里面不是一样有线式和隐式可以使用
回复: 【讨论】fluid求解第一步总是通不过
辛苦辛苦!!
页:
[1]