通过差分计算一维扩散

问题描述:

通过差分计算一维扩散_第1张图片
扩散公式:
离散成差分公式:

同样我不写这个差分的程序了。这个流体力学的差分相对于刚刚写的弹性力学的差分,感觉简单许多。

板1为上面的,板2为下面的。两板间距1cm,划分为11个节点。即空间步长为0.1。时间步长选取为0.001,在0到0.1秒内 板1的速度斜坡上升。0.1秒后稳定为10cm/s
通过差分计算一维扩散_第2张图片
这是层流,与板1接触的流体速度与板1相同。同理与板2接触的流体,速度为0.   在上表中,填入初始条件与边界条件。则可以通过差分公式,计算其余的节点值。
举例:
计算z=0.9处的流体在0.001s时的速度 即f(9,1)。
通过差分计算一维扩散_第3张图片  我懒得打公式也不怎么会用markdown,所以就手画了。
同理,其余节点也可以算出来。

结果:
通过差分计算一维扩散_第4张图片
通过差分计算一维扩散_第5张图片

你可能感兴趣的:(数值分析)