声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

声振论坛 展示 基础理论 流体力学 查看内容

[系列]Python计算CFD问题

2016-4-5 14:45| 发布者: aspen| 查看: 11711| 评论: 0|原作者: funi|来自: lorenabarba.com

摘要: Python计算CFD问题1:一维线性对流问题 #-*-coding=utf-8 -*- import numpy as np import matplotlib.pyplot as plt import time,sys nx=41 #空间节点数 dx=2.0/(nx-1) #网格间距 nt=25 #时间步数 dt ...
[系列]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


最新评论

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-4-28 07:09 , Processed in 0.077281 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

返回顶部