问题描述:
扩散公式:
离散成差分公式:
同样我不写这个差分的程序了。这个流体力学的差分相对于刚刚写的弹性力学的差分,感觉简单许多。
板1为上面的,板2为下面的。两板间距1cm,划分为11个节点。即空间步长为0.1。时间步长选取为0.001,在0到0.1秒内 板1的速度斜坡上升。0.1秒后稳定为10cm/s
这是层流,与板1接触的流体速度与板1相同。同理与板2接触的流体,速度为0. 在上表中,填入初始条件与边界条件。则可以通过差分公式,计算其余的节点值。
举例:
计算z=0.9处的流体在0.001s时的速度 即f(9,1)。 我懒得打公式也不怎么会用markdown,所以就手画了。
同理,其余节点也可以算出来。
结果: