偏微分方程数值解法
偏微分方程数值解法
遗世独立的理想乡_
高性能计算、计算流体力学、数值分析
展开
-
计算流体力学简介(九)——拉瓦尔喷管模拟
控制方程首先根据前面一维无粘流动的控制方程u=[ρρuE],F=[ρuρu2+pu(E+p)],E=pγ−1+12ρu2 ∂u∂t+∂F∂x=0u=\left[\begin{matrix}\rho \\\rho u\\E\end{matrix}\right],F=\left[\begin{matrix}\rho u\\\rho u^2+p\\u(E+p)\end{m...原创 2020-02-03 11:35:13 · 17518 阅读 · 0 评论 -
计算流体力学简介(八)——激波管问题
问题描述原创 2020-01-24 16:11:34 · 11112 阅读 · 3 评论 -
计算流体力学简介(七)——频散关系保持(DRP)格式
色散误差前面的计算中分析了一阶格式的格式粘性,并且提到了一阶格式的格式粘性是泰勒展开的最低阶误差项,因为格式是一阶所以最低阶误差就是二阶项,表现为粘性,因此称为格式粘性。但是前面也有提到泰勒展开的误差项中偶数阶误差项表现为耗散,奇数阶表现为色散。因此其实一阶格式的最高阶误差应当分为奇数阶和偶数阶分别分析。当保留三阶误差项时即可得到一个色散项。这个色散项会使得各频率分量以不同与微分方程描述的传播速...原创 2020-01-24 16:11:03 · 2892 阅读 · 0 评论 -
计算流体力学简介(六)——有限体积法和基本无震荡(ENO)格式简介
基本原理有限体积法和通量差分在网格划分上基本相同,原理上也差不多。网格划分上,有网格点在网格的中心也有网格点在网格顶点的。不论哪种网格划分,网格边界都是通量点,每一步计算都只有网格点上有值,通量点上的值是利用插值临时计算出来用于推进的。这里我就只做网格点在网格中心的情况。对于方程∂u∂t+∂f∂x=0\frac{\partial u}{\partial t}+\frac{\partial f...原创 2020-01-24 16:10:45 · 3436 阅读 · 0 评论 -
计算流体力学简介(五)——通量差分
基本原理原创 2020-01-24 16:10:15 · 1862 阅读 · 0 评论 -
计算流体力学简介(四)——有限差分法
基本原理有限差分法的原理很简单,就是利用差商代替微分,利用泰勒展开得到微分项的差商近似,从而求出方程的近似解。记差分算子Ekf=f(x+kh)E^kf=f(x+kh)Ekf=f(x+kh),微分算子Dkf=∂kf∂xkD^kf=\frac{\partial^kf}{\partial x^k}Dkf=∂xk∂kf利用差分算子和微分算子表示泰勒展开有Ef=If+hDf+12h2D2f+⋯+1...原创 2020-01-24 16:10:04 · 6112 阅读 · 5 评论 -
计算流体力学简介(三)——有限元
一维对流方程原创 2020-01-24 16:09:45 · 2579 阅读 · 0 评论 -
计算流体力学简介(二)——伪谱法
一维对流方程∂u∂t+c∂u∂x=0\frac{\partial u}{\partial t}+c\frac{\partial u}{\partial x}=0∂t∂u+c∂x∂u=0,取c=−1c=-1c=−1,现用伪谱法求解。指数初场原创 2020-01-24 16:09:19 · 4184 阅读 · 1 评论 -
计算流体力学简介(一)——一些基本概念
偏微分方程与常微分方程偏微分方程和常微分方程的区别主要就是体现在待求解函数是一元函数还是多元函数。多元函数存在对不同自变量的偏导数,因此这种带有多元函数偏导数的方程就是偏微分方程。而对于只有一个自变量的一元函数,由它的各阶导数组成的方程就是常微分方程。描述带阻尼的简谐振动的方程md2xdt2+νdxdt+kx=0m\frac{d^2x}{dt^2}+\nu\frac{dx}{dt}+kx=...原创 2020-01-24 16:09:12 · 5028 阅读 · 0 评论