03_控制方程的离散化方法.pdfVIP

Elective Course for Graduate Students 3.1 概述 什么是离散化 CFD计算之前，要对计算区域进行离散化，包括空间 3 控制方程的离散化方法 离散和时间离散(瞬态问题) 。 空间离散是对空间上连续的计算域进行划分，把它 划分成许多子域，并确定每个区域中的节点，从而生 成网格。然后将控制方程在网格上离散，即将偏微分 孙晓颖 格式的控制方程转化为各个节点上的代数方程组。 Harbin Institute of Technology 1 离散化的目的 3.2 常用的离散化方法 在计算域上建立偏微分方程理论上有真解，但由于 有限差分法(Finite Difference Method)—点近似 问题的复杂性，很难获得真解，需通过数值方法把计 基本思想：将求解域划分为差分网格，用有限个网格 算域内有限数量的因变量作为未知量处理，建立一组 节点代替连续的求解域，然后将偏微分方程(控制方程) 关于未知量的代数方程。 导数项用差商来代替，推导出含有离散点上有限个未 偏微分方程定解问题的数值解法分成两阶段： 知数的差分方程组，求解方程组的近似解。 1) 用网格线将连续的计算域划分为有限离散点(网格 FDM法是一种直接将微分问题变为代数问题的近似 节点)集，并选择适当路径将微分方程及其定解条件转 数值解法，数学概念直观，表达简单，是数值解法中 化为网格节点上相应的代数方程组，即建立离散方程 最古老也是最经典的方法。 组。 多用于求解双曲型和抛物型问题；对边界条件复杂、 2) 在计算机上求解离散方程组，得到节点上的解 尤其是椭圆型问题不如有限元法或有限体积法方便。 2 3 差分格式 有限差发法的出发点是守恒型方程的微分形式，求 1) 从格式的精度来划分，有一阶格式、二阶格式和 解域用网格覆盖，在每个网格节点上，通过将偏导数 高阶格式。 近似为节点函数值的代数方程式来近似原偏微分方 2) 从差分的空间形式来考虑，可分为中心格式和逆 程。最终在每一个节点上都有一个代数方程，方程的 风格式。 未知数是中心节点以及相邻节点上的变量值。