人工粘性

要说明白CFD计算中的人工粘性,便没有办法绕开其他几个概念,他们分别是:

截断误差

舍入误差

修正方程


先说舍入误差(它最好理解):

计算机因为数位的限制,在求解方程的每一步过程中都会对精确的结果进行一定程度的四舍五入,久而久之,积年累月,最终累计出来的误差便是舍入误差。

舍入误差外传:

舍入误差是有可能随着计算过程中每一步的迭代而逐渐增大的。(仅针对显示差分情形!显示差分情形!显示差分情形!)。也就是说在显示差分格式时,若想让求解过程稳定,即舍入误差不会不停的增大,就必须满足一定的条件(比如时间步长与空间步长不能太长等)。不同形式的方程有不同形式的条件,较为知名的限制条件为CFL条件(柯朗数),即柯朗数必须小于1,最多=1.



再说离散误差(因为它最直观):

考虑一维热传导方程

一维热传导方程典型形式(方程1)

CFD就得离散。我们对空间进行中间差分,对时间进行一阶的向前差分。得到

离散后的一维热传导方程(方程2)

先别管误差在哪儿,方程1和方程2的精确解(无舍入误差)是不一样的,两个精确解之间的误差由于离散导致,该误差称之为离散误差。离散误差+舍入误差=计算总误差。


修正方程:

考虑一维波动方程:

方程(3)

离散了它(时间一阶向前,空间一阶向后),得到

方程(4)

对方程(4)进行泰勒展开,求导,再加加减减,变换了个样子

方程(5)

方程(5)叫做修正方程。注意方程(4)和方程(5)是完全一致的。但看方程(5)和原方程(3)相比,左端项一致,右端由0变成了一长串二阶三阶导数项,多出来的这些项就是截断误差。截断误差中的偶数阶项叫数值耗散,奇数阶项称为数值色散。


人工粘性:

仔细观察修正方程中的偶数阶耗散项,是由一个系数*变量的偶数阶导数构成。与NS方程中的粘性项对比,该系数便类似于粘性系数μ。当然这个系数是数值方法引入的,人为的引入的,因此称为人工粘性。人家μ是物理粘性,是实实在在真实存在的。

人工粘性外传:

人工粘性的影响:一定程度上抹平了间断,降低了解的精度。

好处:人工粘性通常有助于提高解的稳定性。事实上经常我们还嫌弃人工粘性不够大,人为的再加点儿更多的人工粘性!得到一个不是那么准确的解总比得不到一个解要好嘛。反正掌握好一个度。

问题:人工粘性怎么加?人工粘性对计算结果的影响到底有多大?感觉像是个玄学问题,没详细研究


针对人工粘性的先进的算法:自动的在需要增加人工粘性的位置添加大小合适的人工粘性,屌不屌就问你。TVD就是这么一个概念。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,233评论 6 495
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,357评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,831评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,313评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,417评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,470评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,482评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,265评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,708评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,997评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,176评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,827评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,503评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,150评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,391评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,034评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,063评论 2 352

推荐阅读更多精彩内容

  • 1.为自己即将开设的店铺做一个基础定位,包括目标用户,产品定位,进货渠道,预计的SKU等,并给出开设的店铺地址。 ...
    3edd52f615c0阅读 203评论 0 0
  • 清晨温柔的阳光透过绿色的叶子洒向大地,年轮中曾经青涩的时光向未来讲述着美丽的传说。 用这段文字开场...
    郝吉庆阅读 194评论 0 0
  • 柠檬不是感光,柠檬皮里的柠檬精油才感光,芹菜的感光就是吃下去以后的感光。 人的皮肤生下来就是有隐形斑在基底层,那个...
    杨小凡_9d5c阅读 127评论 0 0