[系列]Python计算CFD问题<4>:一维扩散方程
与对流方程不同,扩散方程是2阶的偏微分方程。一维扩散方程可用下面的方程进行描述: [attach]69946[/attach] 该方程可以离散为: [attach]69945[/attach] 写成迭代的形式为: [attach]69947[/attach] 采用与前例相同的初始条件,及0.5<=x<=1时,u=1,其他位置u=2。扩散系数mu=0.3。 程序代码为: [attach]69948[/attach] 利用函数diffision可查看各时刻速度分布(时间步数分别为0、50、100、150、200、250)。 [attach]69949[/attach] [attach]69950[/attach] [attach]69951[/attach] [attach]69952[/attach] [attach]69953[/attach] [attach]69954[/attach] 放到一起看起来可能更直观一些: [attach]69944[/attach] 【来源:hhttp://nbviewer.ipython.org/gith ... ons/04_Step_3.ipynb】 |
GMT+8, 2024-11-24 22:08 , Processed in 0.043931 second(s), 17 queries , Gzip On.
Powered by Discuz! X3.4
Copyright © 2001-2021, Tencent Cloud.