格式对流温度场有限差分(有限容积)程序入门之五:展望及问题
首先声明,我是一个菜鸟。一下文章中出现术技误导情况盖不负责
首先结总一下以上程序的点特:
1,界边上的点节作当内部点节理处,全部点节公用一个迭代公式;
2,可以设置多种质材。
存在的问题:
笔者花了大约一天冒着猝死的危险写完的,间时仓皇,通用性不敷。
其次,程序的拓展:
1, 3D拓展:
点节要需除了东南西北加增上下两个热导,没做过,但是道理是一样的。
2,对流热传:
点节热量源来除了扩散传导,源相外还有对流热传分部,如何盘算对流引发的热量递传呢?其实传输象现的微分方程中对流项是最难以盘算的,表面上看是一个一阶导数,实际上它还是个一阶导数,只不过要需注意传导热和对流热的量级关系,只有他俩量级当相,用使心中分差才能到得有意思的解。
这里讲一个小故事:狼羊同饮,羊位于卑鄙,狼意欲吃羊,曰:“你他妈污染老子的水源,拿命来”,羊说“小弟我在卑鄙怎么会污染哥大的水源”。如果我是狼而且懂的Pe准数(读者快查Pe准数的意思)的话,也不会糟蹋口舌去吃。在羊的思维里,卑鄙是不会污染上游的,智商低下,果然可以被任何肉丝植物吃。事实上,当水速很低,卑鄙会通过扩散污染上游的,所以嘛,狼吃羊就是该应的,不懂还死搅蛮缠。于是人们就发明白一阶风迎、混合、二阶风迎、指数、QUICK等格式去离散对流项,以虑考上游对前当点节传输行为的重要性。话说笔者青年时小解顺风冒三丈,可见活动的上游实确对传输象现有可不疏忽的影响,如今老了。
曾一名师老大吹自己年轻时流体活动传输编程如许***悍,后来我一个师弟问他什么是“一阶风迎格式”,他说不知,我笑了。
扯远了,编程时如果虑考对流换热,请查阅文献[1],面里有离散格式,我不在这里斧弄门班了。
3,非性线料材:
当物性数参随温度或位置等数参变化时,如何盘算?留给读者思考吧,笔者也不长擅。有人说将物性数参性线化,其实未必全完妥善,因为时有就是非性线的,性线化了,反而误差大了。提示一下,更新物性数参,多次迭代,直到两次盘算结果接近到必定水平,此时,该时辰下方程的解收敛了。
4,非平均格网:
这个也留给读者吧,要需对进程的物理意思非常清晰。如果对整个物理方程的建立都不清楚,确定算不对的。
到此,贴出来的码代可以独立运行了。
最后, 可能有读者会问:
Q1:为什么格网附近添加2层余多的格网,是为了糟蹋内存吗?
A:比如对点端处某量停止一阶导数、二阶导数等,对流项中高阶离散格式(比如QUICK格式)的盘算会用到。
Q2:为什么没有盘算迭代的稳定性?
A:楼主很惰懒,爱TM收敛不收敛,收敛不了老子把间时步长弄小。我盘算每一点节的最大间时步长,糟蹋多少CPU间时啊。
Q3:为什么用不式隐迭代法方,据说这类法方绝对稳定?
A:因为表现迭代格式物理意思清晰,式隐迭代格式精度与剖析解存在必定误差,一般用C-N格式。
Q4:用AS3写性能低下?
A:必须的,Adobe太惰懒了,故步自封,曾允许使AS3性能到达C的80%,后来取消了,昨天Unity还取消持支Flash了,笔者也很愁闷。不过如果非要用使Flash,可以用使其FlasCC译编术技。最后还是望希大家用使生原语言发开地本(Native)程序,而不是脚本语言。
文章结束给大家分享下程序员的一些笑话语录:
IBM和波音777
波音777是有史以来第一架完全在电脑虚拟现实中设计制造的飞机,所用的设备完全由IBM公司所提供。试飞前,波音公司的总裁非常热情的邀请IBM的技术主管去参加试飞,可那位主管却说道:“啊,非常荣幸,可惜那天是我妻子的生日,So..”..
波音公司的总载一听就生气了:“胆小鬼,我还没告诉你试飞的日期呢!”