一种求解对流扩散问题的高精度无网格方法.pdfVIP

中国王程热物理学会， 传热传质学 学术会议论文 编号：073150 一种求解对流扩散问题的高精度无网格方法 吴学红李增耀马良栋 陶文铨 两安交通大学能源与动力工程学院，西安，710049 1’el：029—82669l 06． E—mail： wqtao@mailxjtu．edu．cn 摘要：无网格数值汁算方法是近’f．几年发展起来的一种新的数值计算方法，该方法主要是在计算区域内 布置一些离散的点，并利用节点的权函数近似地表示其影响域内的物理量，进而形成与节点物理量相 是～种商精度的数值计算方法，。能有效的计算商Pe数的流动问题。 关键词：无网格法，Petro．Galerkin方法，数值模拟 1引言 近十几年来，有限元法与有限容积法成为解决工程问题的主要计算方法。在应用这 些方法的过程中，遇到的困难之一是对复杂几何形状的网格剖分，网格质量的好坏直接 影响到求解问题的精度，这样就对网格的质量提出了很高的要求。鉴于此问题，近些年 来国内外许多学者提出了无网格方法I卜31。 最初，无网格法主要用来求解力学问题。在近JL年，一些作者用它来求解传与流动 方法需要用背景网格进行积分，因此它不是一种真正的无网格法；而光滑粒子动力学方 法则对于非均匀布置盼节点不满足一致性条件，且在边界处的精度较低。1998年Atluri 和zhu吲提出了无网格局部petrov．GaJerkjn(meshIess 它是用移动最小二乘法建立场函数的近似，在局部子域上进行积分，是一种真正的无网 并把它与有限容积的高阶格式进行比较，说明它是一种高精度的计算方法。 2Simth．Hutton问题 甜=2y(1一x2)，v=一2x(1～y2) 边界条件： 在一l≤工≤O，y=0： 乙(x)=1+tanh【a(1+2x)】 ． 在o≤xl，y：0： 娶：q：o 鲫 ～ 其它边界上： r=0 其中口称为陡度系数，其值越大，进口阶梯形温度分布越陡。其流场与进口温度场分布 如图1所示： y X 图lSimth．Hunon问题的流场与进口温度(口=100) 用二维、稳态无源项的对流．扩散方程来求解这一问题，引入Pe数作为参数，’该 基金项}：：：l：国家自然科学綦金重点项目(No的资助。 1155 方程可表示如’F： a丁 1 a271 哟瓦2瓦碍 3MLPG方法的实施 无网格局部Petrov—Galerkin(MLPG) 方法是在节点X处的局部予域Q。 (Q，∈Q)上进行积分，而且此局部子域 可以是任意的形状(如右图2)。为了对式 (1)在点X的局部子域上进行积分，j4j 权函数w与式(1)相乘，其积分方程如 ‘ 下： 心署一去器w皿一㈤ 图2移动最小二乘近似的试探函数的定义区域及 节点的支持域(其中r。为第二类边界条件，r，为 对上式第二项进行分步积分，并用 第一类边界条件) Gauss定理可得： Jc吩詈w+去善考，拯，一砖器，w卯=。 c3， 对于式(3)，第二类边界条件可以直接代入上式进行积分，但它不能直接处理第一 类边界条件；本文应用转换法处理第一类边界条件，详细描述参见文献[1]。 为了得到点X在子域Q，上离散方程(它包含N个节点)，