第二章－数值计算2008.08.pptVIP

第二章 物理现象的数学描写 主要包含以下内容： §2.1 控制微分方程 §2.2 坐标的性质 目的：熟悉控制传热、流动及其它有关过程规律的微分方程的形式和意义（这些方程形式上的一致性是构成一个通用解法的基础），并讨论方程中一些变量的性质。 §2.1 控制微分方程 §2.1-1 微分方程的意义 1. 意义：将要讨论或曾经讨论过的各个微分方程各自代表着一定的守恒原理，如连续性方程、动量方程、能量方程。方程本身代表着影响因变量的各个因素之间存在着的某种平衡。 2. 因变量：通常具有“比”的性质，即以单位质量为基础，如质量分量（一定容积内所含组分 l 的质量与相同容积中所包含的混合物的总质量之比kg/kg）、速度（单位质量的动量）以及比焓 。 温度例外 从通用的控制方程出发建立离散方程，使所建立的离散方程既可适用于不同的坐标系，又能包括不同类型的物理问题； 把整个程序分成不变的通用部分及与具体问题有关的用户部分两块。不变部分在该程序的适用范围内不必作任何改动，而用户部分则与所计算问题的初始条件、边界条件及物性参数有关； 采用分块式结构来写程序，即把整个程序分成一些相互间有相当独立性的段或块，每一块有一定职能，并以其输入输出与其它块相联系。 整个程序应有良好的结构，逻辑清楚，按顺序一块一块地执行。不必为节省很少的计算时间而采取难懂的算法； 离散方程的建立与代数方程的求解应属于两个不同的块，以便于进行差分格式与求解方法的改动。 3. 利用坐标的性质来区分微分方程的类型 抛物型的微分方程表示一种单向状态，而椭圆或双曲型则是表示双向状态的。如果能把物理状态描述成对一给定坐标表示的抛物型或是椭圆型的，将是非常有意义的。 抛物型方程描写物理学中一类步进问题，这类问题中因变量与时间有关，或问题中有类似于时间的变量，故称初值问题。非稳态导热问题—抛物型问题。实质上是时间坐标上的抛物型，空间坐标上的椭圆型。 偏微分方程有3种类型：抛物型、椭圆型和双曲型。 稳态导热问题对所有的坐标都是椭圆型的；对于边界层类型的流动与换热问题，因为略去了主流方向的扩散作用，所以以主流方向坐标表示的是抛物型的，而以横向的坐标表示时则为椭圆型的。 椭圆型方程描写物理学中一类稳态问题，这类问题的变量与时间无关，需要在空间的一个闭区域内求解。 二维边界层流动方程： 结论：至少存在一个单向坐标的状态是抛物型状态；否则就是椭圆型状态。所以边界层型的流动与换热属于抛物型问题。 双曲型：具有一种单向的特性，但这种特性并不是沿坐标方向的，而是沿着一些称之为特征线的特殊线的方向。所以，在此把双曲问题当成是椭圆问题这一总类的一个组成部分（即全部是双向坐标）。 所研究的问题中有一个空间坐标是单向的，就称这种流动或换热为边界层型的问题；如果所有的坐标都是双向的，就称之为回流流动。抛物型与椭圆型是从数学的角度来命名的，而边界层型与回流型则是物理意义的称谓。 4. 边界层流动与回流流动 * 太 原 理 工 大 学 * / 36 Thermal 3. 方程项的意义：微分方程的各项代表着以单位容积为 基础的效果。 密度 4. 几个基本概念 流量密度：单位时间、单位面积流过的质量, 其单位 是kg/(s·m2)，实质上相当于ρu的单位。 （ ρu A） 因变量?的流量密度： ρu ? ，如比焓的流量密度 ρuh? kg/(s·m2) · J/ kg? J/(s·m2) ?W/m2 设想用 符号J 代表一个典型因变量? 的流量密度 单位体积流出的净流量的表示： 控制容积上的流量平衡 通过该面整个面积上流出的净 流量： 同理可写出y、z方向的净流量… ∴ 单位体积流出的净流量为： dx dy dz 因变量?的变化速率： 对于变量比焓： 两者具有相同的量纲。 一个微分方程是这样一些项的组合，每一项代表一个以单位容积为基础的效应，而所有项合在一起反映着某种平衡或守恒。 接下来的任务：通过几个标准的微分方程，找出一个通用的形式 §2.1-2 化学组分的守恒 1.组分方程 令 ml 代表一种组分l 的质量分量，当存在速度场u 时， ml 的守恒方程为： 2.方程各项的物理意义 单位体积内化学组分 l 的质量变化率; 组分l 的对流流量密度 组分l 的扩散流量密度 单位体积 组分 l 的净流量 质量分量: 一定容积内所含组分 l 的质量与相同容积中所包含的混合物的总质量之比。 单位 kg/kg 对流流量密度是由一般化流场ρu 所携带的流量密度; 扩散流量密度由ml 的梯度引起，用菲克扩散定律表示 A B 组分A、B的相互扩散 菲克定律：组分A及组分B 被分隔在容器的两侧，抽去隔板时，由于浓度梯度的存在，组