CFD初学习
不知所云的CHDer
水力学及河流动力学专业
展开
-
一维常系数对流方程的学习——来自流沙公众号
一维对流扩散方程的学习——来自流沙公众号认识:1、时间步长越小越靠后移动;网格越小则波的形状越一致,波形失真在减小,引出Courant数:u*dt/dx<sigma;2、对numpy中的ones()、linspace()、zeros()有基础认识;3、再次练习了matplotlib中的pyplot。4、明白了一维方程编程的简单运算,基础import numpy as np# n...转载 2018-09-25 22:07:35 · 7939 阅读 · 1 评论 -
黎曼不变量——特征法的学习
一、特征法参考这个网络讲义1.苏剑林. (2017, Dec 07). 《一阶偏微分方程的特征线法 》[Blog post]. Retrieved from https://kexue.fm/archives/47182.知乎解答:https://www.zhihu.com/question/40518374二、黎曼不变量在非齐次项为零时,变量沿着特征线不变,故称为黎曼不变量:dw1/r1...原创 2018-12-06 11:43:26 · 4384 阅读 · 1 评论 -
浅水方程的的精确黎曼求解器——推导浅水方程(一)
一、浅水方程二维浅水方程可根据无粘N-S方程沿水深方向积分可得到,主要运用莱布尼茨公式。运用公式后进行合并,图中给出三维的合并:浅水方程的连续方程到此就推出了将欧拉方程(无粘N-S方程)中的压力梯度项做出假设:压强分布满足静水压强分布进一步得到简化公式:这里的s代表水位s=z将连续方程两边乘以u,动量方程乘以h:合并便可以得到u方向的动量方程:同理v方向至此,二维...原创 2018-11-16 11:29:11 · 9027 阅读 · 0 评论 -
对二维扩散方程的学习——来自流沙公众号
代码:import numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dnx = 31 # 定义x方向节点数ny = 31 # 定义y方向节点数nt = 17 # 定义时间步长nu = 0.05 # 扩散系数dx = 2 / (nx - 1)dy = 2 / (ny ...转载 2018-10-16 21:04:40 · 11924 阅读 · 1 评论 -
对二维线性对流方程的学习——来自流沙公众号
许久没接触python,又有点忘了继续学习流沙公众号对这个方程,非定常项,以及x和y方向的对流项碰巧今天看安德鲁的计算流体力学偏微分方程所讲,利用克莱姆法则和特征值法来判断方程的属性例子为一阶,在做习题碰到二阶。此为闲话同时,根据文章所写,稍加修改,将for循环和数组放在同一程序中进行对比,以温习python,发现遗忘甚多。from mpl_toolkits.mplot3d imp...转载 2018-10-16 19:14:52 · 2922 阅读 · 1 评论 -
对Laplace方程的学习——来自流沙公众号
对Laplace方程的学习代码如下:import numpy as npimport matplotlib.pyplot as pltimport timedx = 0.1dy = 0.1dx2 = dx*dxdy2 = dy*dydef laplace(u): nx,ny = u.shape #查看数组的维数 for i in range(1, nx-1): ...转载 2018-09-28 22:13:12 · 1225 阅读 · 0 评论 -
一维Burgers方程的学习——来自流沙公众号
对二维Burger方程的认识,该方程import numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dnx = 41ny = 41nt = 120c = 1dx = 2 / (nx - 1)dy = 2 / (ny - 1)sigma = 0.0009nu = ...转载 2018-10-17 22:53:34 · 2128 阅读 · 0 评论 -
一维Burgers方程的学习——来自流沙公众号
伯格斯方程(Burgers equation) 是一个模拟冲击波的传播和反射的非线性偏微分方程。伯格斯方程也称为粘性伯格斯方程,只适用于一个空间维度。转载 2018-09-27 21:07:16 · 23040 阅读 · 12 评论 -
一维扩散方程的学习——来自流沙公众号
来自流沙公众号认识:1、理论:差分离散,较为基础2、第一次接触matplotlib中ion()和ioff()。ion()开启了交互模式,ioff() 感觉像暂停图像同时,查阅其他博客得到:python可视化库matplotlib有两种显示模式:阻塞(block)模式交互(interactive)模式在Python Consol命令行中,默认是交互模式。而在python脚本中,ma...转载 2018-09-26 21:29:30 · 11408 阅读 · 5 评论 -
一维扩散方程差分格式的数值计算
摘要:采用FTCS显格式、Dufort-Frankel显格式、Lassonen隐格式、Grank-Nicolson隐格式对一维扩散方程进行数值计算,得到不同时间y方向扩散的速度。结果表明了不同差分格式的差别。关键词:扩散方程;差分格式;FORTRAN;module global implicit none real::u(1000,1000)end module prog...原创 2018-12-24 16:19:12 · 12834 阅读 · 5 评论