2012年5月10日星期四
CFD and adina
Today finally understand clear about CFD how it works.
there is actually not such thing call pressure-pressure condition, it sounds i am too silly, as there is a redundant condition, you speficy the velocity condition, at the same time, you are also manage to settle the pressure , as in cfd, both side are equalized. one you give magnitdue, another one it will automatically calculate for you. i finally can go into CFD.
i also manage to solve adina problem
给版主做个补充:overlap不单指流场网格和结构场发生重叠。
由于adina在FSI求解时,采用的是ALE列式,简单些说即流场和结构场交接区域的网格是可变形的,当发生“强耦合”运动时,或者网格初始质量不太好(比如网格的初始夹角太小、塌陷率过大等),一旦网格运动,很容易导致网格变形后出现“雅克比阈值”接近0甚至小于0,产生overlap。
各个CFD程序内置的overlap标准也不同,原因在于网格出现压缩或者拉伸后回直接影响积分精度,因此,理论上,这个阈值设置高一些,精度保证就高一些,但出现overlap的概率就大一些。
要解决overlap,根本上只能从网格初始质量寻求解决办法。当然,adina内置的一些辅助措施,如设置leader follower对于缓减overlap也是很有效的,甚至有些时候会让“貌似不可算”(主要指初始网格质量很差,能保证开算,但很容易就overlap的模型)的网格顺利计算完毕。
有兴趣的可以试试看,在稳态流场计算时,只要模型初始网格雅克比阈值大于0,是不会遇到overlap的,原因在于稳态流场中一般用纯欧拉格式的有限体积法,网格是不会变形的,差别只是收敛时间的长短或者难易。
简单些说,雅克比矩阵就是一个多维函数在某点位置对应的各个维度上的偏导数矩阵(该导数矩阵使得改点对应母函数具有线性最优逼近特征),如果偏导数矩阵小于0 ,表示函数在改点存在反向,否则保持原有方向。对于一个单元来说,如果在节点位置出现反向则意味着什么呢?你可以思考一下,再进一步,对一维网格、二维网格、三维网格分别意味着什么呢
http://www.forumadina.com/forum.php?mod=viewthread&tid=209&extra=page%3D1
http://forumadina.com/forum.php?mod=viewthread&tid=274&highlight
overlapp的判断标准问题
订阅:
博文评论 (Atom)
没有评论:
发表评论