|
发表于 2004-9-16 22:58:00
|
显示全部楼层
疑问:非结构网格单元节点编号的问题。
是酱子的。解方程时候,方程是不分先后的。你可以认为瞬态问题是稳态问题的叠加。但不要在什么学术会议上面或者答辩上这么说的,因为这么说其实不科学。另外求解时候要注意隐性和显性格式。
对于每一个时间步长deltaT, 后面一个时间T(n+1)道理可以这么说是和前面一个时间T(n)有关系,但离散格式又分显性和隐性,显形格式比较好处理,所有T(n+1)时间的参数比如速度,都是完全由上一个时间算出的,比较简单,直接累加计算就可以了。隐性格式相对麻烦,因为T(n+1)的参数(比如速度)还和本时间的状况有关,方程是线性方程组,要编程高斯消去或者高斯迭代解的。
大家都喜欢用显性格式,但隐性格式有时候是必须的,因为显性格式很多时候不收敛并且让人觉的水平较低。我上CFD课时候,老师明写不许用显性格式,虽然显形也能做,呵呵
另外时间步长也很重要,太大了不收敛,太小了浪费计算资源。如果英文够好,推荐看下面这本书,不用全看,看头几章,CFD的经典教材。我觉的比林林总总国内的CFD书强多了。
非结构化2维还好说,3维就麻烦了。我的理解是,经典的CFD计算都是以结构化网格为基础的,对于2维来说就是,一个计算节点,它上下左右4个方向各有一个点。但如果是非结构化的,比如3角网格,我还真不知道怎么离散合适。
"Computational Fluid Mechanics and Heat Transfer"
by John C. Tannehil, Dale A. Anderson, and Richard H.
Pletcher
|
|