十二步解N-S方程之第四步

2019独角兽企业重金招聘Python工程师标准>>>

十二步解N-S方程之第四步

通过前面的学习，我们已经知道如何去解线性和非线性的一维对流方程，而且也知道CFL是如何影响数值求解的。在这里我们考虑一维扩散方程，如下所示：

import numpy as np
import matplotlib.pyplot as plt
#
nx=41
dx=2.0/(nx-1)
nt=20
nu=0.3
sigma=0.2
dt=sigma*dx**2/nu
#
u=np.ones(nx)
u[0.5/dx:1+1/dx]=2.0
#
u0=np.ones(nx)
u0=u.copy()
#
un=np.ones(nx)
#
for n in range(nt):
un=u.copy()
for i in range(1,nx-1):
u[i]=un[i]+nu*dt/dx**2*(un[i+1]-2*un[i]+un[i-1])
plt.figure()
xc=np.linspace(0,2,nx)
plt.plot(xc,u,linewidth=2.0,label="converged")
plt.plot(xc,u0,linewidth=2.0,label="initial")
plt.legend(loc='upper right')
plt.savefig("F:\PYTHON\python_cfd\step4.png",dpi=150)
plt.show()

运行结果如下图所示：

转载于:https://my.oschina.net/cfdvalidation/blog/309807

十二步解N-S方程之第四步相关推荐

五步移相法matlab程序_四步移相法提取相位（相位解包裹）基础问题
最近在自学四步移相法提取相位,有个初级疑问: 有四幅移相图片a.b.c.d如下:完整解包裹代码如下: 程序运行后,主要是figure(7)和figure(8)没看太明白,;从figure(5)和fig ...
[算法课]全面翻新计划！第十二周全解
文章目录上课部分回顾: 例1字典序最小问题 PS: 颜老板版本更新版吐槽思路 AC代码安全版非安全,极简版思考:任意一串字符,删除其中的重复字符! 题目描述 AC代码哈希表哈希表 ...
从零开始学USB（二十二、USB接口HID类设备（四）_报表描述符Local类）
本地项目标签定义控件的特征. 这些项目不会转移到下一个主要项目. 如果Main项定义了多个控件,则可能会有几个类似的Local项标签. 例如,一个Input项可能有几个与之关联的Usage标签,每个控 ...
Flink实战（八十五）：flink-sql使用（十二）Flink 与 hive 结合使用（四）Hive Read Write
0 简介 Using the HiveCatalog and Flink's connector to Hive, Flink can read and write from Hive data as ...
HyperLynx（三十二）高速串行总线仿真（四）
高速串行总线仿真(四) 在上节的基础上 1.通过导出到LineSim验证一个串行通道 2.快速眼图仿真 3.高速串行总线设计注意事项 1.通过导出到LineSim验证一个串行通道接下来,将导出网络到 ...
Python之基础详解（十二）：进阶，用pyecharts绘制世界地图、玫瑰图、雷达图
我之前发过一篇博文基本讲解了基本绘制地图的方法. 我这里打算直接调用api来创建实时动态地图,并且分析数据绘制了玫瑰图,雷达图. 步骤配套资源下载: 第一步:配置环境第二步:爬取数据第三步:创建 ...
今天，「神舟十二号」发射！这次载人航天有什么不一样？
来源|爱范儿文|周宇上个月,天问一号成功登陆火星,并释放出我国首辆火星车「祝融」,在一次发射中出色地完成了围绕火星飞行.降落.巡航的任务,堪称人类航天史上的壮举. 而今天,我们将再进一步. 就在今 ...
大二上学期软件工程概论学习进度表（第十二周）
第十二周所花时间(包括上课) 四节课有4小时,自己编程3小时, 写作业2小时,实训8小时.共17小时. 代码量(行) 约400行博客量(篇) 2篇了解到的知识点实现用户注册转载于:https ...
开源的调色板软件：焰火十二卷（Rickrack）
一.前言 Adobe Color 是一个配色网站,可以在线生成不同风格的色彩主题,受到设计师的喜爱.然而,Adobe Color 服务器在国外,导致连接速度很慢,甚至会连接不上.并且,Adobe Co ...

十二步解N-S方程之第四步

十二步解N-S方程之第四步

十二步解N-S方程之第四步相关推荐

最新文章

热门文章