常见的计算方法 偏微分方程数值解.ppt

第五章 偏微分方程数值解 Numerical Methods for Partial Differential Equations 5.1 偏微分方程简介 5.2 离散化公式 5.3 几种常见偏微分方程的离散化计算 5.4吸附床传热传质模型中偏微分方程求解 5. 1 偏微分方程简介 偏微分方程 如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 在化工或化学动态模拟方程中，常常有一个自变量是时间，其它的自变量为空间位置。如果只考虑一维空间，则只有两个自变量；如果考虑两维空间，则有3个自变量。 许多化工过程均是通过对偏微分方程的求解进行工艺参数的确定或数值模拟。 5.1 偏微分方程简介 偏微分方程的分类 5.1 偏微分方程简介 数学上的分类： 椭圆方程 Elliptic 抛物线方程 Parabolic 双曲线方程 Hyperbolic 物理实际问题的归类： 波动方程(双曲型）一维弦振动模型： 热传导方程（抛物线型）一维线性热传导方程 拉普拉斯方程（椭圆型）稳态静电场或稳态温度分布场） 5.1 微分方程的求解思路 求微分方程数值解的一般步骤： Step1区域剖分：首先按一定规则将整个定义域分成若干小块 Step2微分方程离散：构造离散点或片的函数值递推公式或方程 Step3初始、边界条件离散：根据递推公式，将初值或边界值离散化，补充方程，启动递推运算 Step4 数值解计算：求解离散系统问题 微分方程的定解问题 离散系统的求解问题 5.2 离散化公式 将自变量在时间和空间上以一定的间隔进行离散化，则应变量就变成了这些离散变量的函数。 一阶偏导的离散化公式 一般采用欧拉公式表示 有时为了保证系统和稳定性， 对时间的差分往往采用向后公式 5.2 离散化公式 对于二阶偏导，我们可以通过对泰勒展开式处理技术得到下面离散化计算公式： 5.2 离散化公式推导 将uk+1在uk处按二阶泰勒式展开： 将uk-1在uk处按二阶泰勒式展开： 二式相加得： 5.3几种常见偏微分方程的离散化计算 5.3.1 波动方程求解 对于初值问题，是已知t=0时，u与 依赖于x的函数形式，求解不同位置，不同时刻的u值。而 u是定义在 的二元函数，即上半平面的函数。 对于混合问题除初值外，还有边值。是已知初值及x=0及x=l 时u依赖于t的函数，求解不同位置x，不同时刻的u值。此时u是定义在 的带形区域上的二元函数。 5.3.1 波动方程求解 5.3.1 波动方程求解 5.3.2 一维流动热传导方程 与波动方程的情形类似，用差商近似代替偏商，可以得到一维流动传热传导方程的混合问题的差分方程，以其解作为流动传热传导方程的近似解。 5.3.2 一维流动热传导方程 将上式进行处理得到： 该式是显式格式。只要保证式中各项系数大于零，一般情况下是稳定的，可以获得稳定的解。 分析上式可以发现，当为了提高数值精度取适当小的Δx 时，最有可能小于零的系数是 uin的系数，若要保证此项系数大于零，此时Δt必须相应地更小,会导致计算量将大大增加，这是显式格式的缺点，为了克服此缺点，下面提出一种隐式格式： 偏微分方程在 点上进行离散化，且对时间的偏微分采用向后欧拉公式得到原偏微分方程的离散化公式： 5.3.2 一维流动热传导方程 5.3.2 一维流动热传导方程 5.3.3 稳态导热/扩散方程 5.3.3稳态导热/扩散方程求解 离散化公式： 取 ，经化简得： 外节点（边界节点）和内节点 求解方法 划分网格 建立节点离散方程 迭代求解（或解稀疏方程组） 5.3.3稳态导热/扩散方程求解 5.3.3 稳态导热/扩散方程求解 5.3.3 稳态导热/扩散方程求解 Microsoft Excel 迭代计算公式中的循环引用 在“工具”菜单上，单击“选项”，再单击“重新计算”选项卡。 选中“迭代”复选框。 若要设置 Microsoft Excel 进行重新计算的最大次数，请在“最多迭代次数”框中键入迭代次数。迭代次数越高，Excel 用于计算工作表的时间越多。 若要设置两次迭代结果之间可以接受的最大误差，请在“最大误差”框中键入所需的数值。数值越小，结果越精确，Excel 用于计算工作表的时间也越多。 5.4 吸附床传热传质模型中偏微分方程求解实例 5.4.1 基本设定及假设 5.4.2　流体传热模型的建立 4.