格式对流温度场有限差分(有限容积)程序入门之五:展望及问题

首先声明,我是一个菜鸟。一下文章中出现术技误导情况盖不负责

      首先结总一下以上程序的点特:

    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..”..
  波音公司的总载一听就生气了:“胆小鬼,我还没告诉你试飞的日期呢!”

你可能感兴趣的:(格式对流温度场有限差分(有限容积)程序入门之五:展望及问题)