MATLAB在CFD中的应用研究.pdfVIP

MATLAB在CFD中的应用 金凤云 刘小兵 陈冬冬 曹军 (西华大学能源与环境学院，成都610039) 摘要 本文通过作者对MATLAB在CFD中的应用实践，详细描述了MATI。AB软件在 CFD中的应用前景，并附以算例及作者的一些编程技巧。 关键词CFD，MATLAB，偏微分方程，数值计算，可视化 1 前言 Fluid CFD(ComputationalDynamics)，即计算流体力学，是由流体力学理论、计算机技术 和数值方法等交叉产生的一门应用基础学科，它与实验流体力学及分析流体力学共同构成了 现代流体力学的一个完整的体系。CFD主要研究内容可分为科学计算和结果的可视化。科 学计算的主要目标是针对流体的流动特性，建立相应的数学模型、采用合适的算法求解流体流 动的偏微分方程，以便对流体的流动进行仿真，准确模拟出流体运动的状态和结构。而计算结 果的可视化是理解仿真信息的关键所在，其主要内容包括网格生成、计算结果的图形显示和可 视化分析，即我们常说的前处理和后处理过程。在这两个研究内容中，计算机技术都起着至关 重要的作用。 如何选择适当的程序语言进行科学计算和结果可视化就成为我们CFD工作者的一项重 要工作。在众多的数学软件中，MATLAB以其优秀的数值计算能力和卓越的数据可视化能 力脱颖而出，成为在CFD研究领域公认的优秀的数学应用软件之一。 算软件，它是一个以矩阵为基本单位的高效数值计算语言，是一个集科学计算、图像处理、声音 处理于一体的高度集成系统。MATI。AB语言被称为第四代计算机的语言，比起在CFD领域 应用广泛的c和FORTRAN等传统语言，其明显的优点是简洁且功能强大，特别是具有强大 的编辑图形界面的能力。并具有语言简单自然，扩充能力强，编程效率高等诸多优点。 本文将讨论MATLAB在CFD中主要应用，并附以两个算例加以说明。 2 MATLAB中的偏微分方程工具箱在CFD中的应用 CFD中的科学计算的核心内容就是对描述流体流动的偏微分方程组进行求解，而MAT— LAB中的偏微分方程工具箱(PDEToolbox)提供了研究和求解空间二维偏微分方程问题的 一个强大而又灵活实用的环境，只要使用用户界面或M文件，画出所需求解区域，输入方程类 型和有关系数，就可以显示解的图形或输出解的数值。 2．1 PDE Toolbox的主要功能 2．1．1 设置PDE(偏微分方程)的定解问题 一般可以直接利用图形用户界面(GUI)来实现对定解区域、边界条件以及方程的形式及 系数的设定。对于标准问题，一般过程为三个步骤：①利用Draw模式，使用CSG(几何结构实 体模型)对话框画几何区域。②利用Boundary模式在所绘几何图形的各个边界段上给出边 界条件。③利用PDE模式来确定方程的类型、系数c，a，f和d。对于不同材料还可以设置不 同的系数，以反映其材料特性。 对于非标准的问题，可以利用PDEToolbox的函数，或者用FEM(有限元)方法求解更为 复杂的问题。 2．1．2 利用FEM求解PDE 即进行网格生成、方程离散、J并求出数值解，其主要模式有：①Mesh模式：利用此模式来 生成网格，并自动控制网格参数。②Solve模式：对于特征值问题，能够求出给定区间的特征 值，对于抛物线型和双曲型方程，设置边值条件后能求出给定时刻的解，对于椭圆型方程还能 求非线性和自适应解。 2．1．3 解的可视化 tor等命令作图。对于抛物型和双曲型方程，还可以根据需要生成解的动画。 2．2算例分析 以抛物型偏微分方程中的热传导方程为例，它描述的是平面的热传输现象，或者是由轴对 称三维问题经过降维后的热传输问题。其流动控制方程的一般形式为： 心瓦8T—V·(kvT)一Q-t-h(T,。一T) 其中： 丁：温度 ID：密度 C：比热容 k：传热系数 Q：内热源热量 h：对流换热表面传热系数 L：环境温度 本算例为不同介质系数的热传输问题，其边界满足齐次Dirichlet边界条件。讨论区域为 一半径为4的圆形，其导热系数为15，密度为l，比热容为0．3。在此圆形区域中心有一正方 形区域，其边长为2，热源为10，导热系数为15，密度为5，比热容为0．1。其计算