基于变分原理的三维扩散方程离散化格式.pdfVIP

篓!登 型竺!!些些塾竺塾!釜些堡竺 鳖：鳖 [文章编号]looI-246x{2003)()4—029l—07 基于变分原理的三维扩散方程离散化格式 陈光南’， 李德元1， 万正苏2 (1北京虚用物理与“算数学研究所计算物理实验宣，北京lO∞sE；2上海夫学数学乐．上海2∞436) [摘要] 研究三维扩散方程的数值模拟．在非正规六面体网格上，崩变分原理建立扩散流形式的差分格式 将扩散流通量作为未知函数，求其泛函极值．并与原函数温度联立求解 详细推导了计算公式，并给m了数值 试验结果． [关键词] 变分原理；非正规六面体咧格；三维扩散方程；差分格式 【中圈分类号] 0241．82 [文献标识码] A 量，求泛函 0引言 (6) F(w)：IL!』dn一2IⅡdivwdo ” 在描述热传导或粒子扩散的抛物型方程的守恒 击 占 的极值，解出通量流．1f．不难推出，在每一个同定时 形式 刻满足边界条件(5)的流函数w+使泛函(6)达到极 (1) 小．因此求抛物型方程(1)满足定解条件(3)，(4)的 警“V(wadM)=o 解转化为在每一时刻先求得在满足边界条件(5)的 中，散度所作用的向量w=一*grad“表示热流或粒 子扩散流数值求解抛物型方程(1)，可以将流函数 函数类中的使泛函(6)达到极小的解，．，，再用方程 作为单独的量求解，这样就化方程(1)为方程组 情形下方程(1)的离散化格式．文[3]讨论r基于变 害=“汕， (2a) 分原理的二维热传导方程差分格式．详细推导了在 w=一ⅣgradM， (2b) 不规则四边形网格卜的计算公式，并在实际问题中 可以通过变分原理建立微分方程(1)的离散化格式 得到很好应用．本文以复杂流场为背景，在三维非正 在三维区域(z，v，z)∈n，‘O上解方程(1)并 规六面体网格上(图1)，基于变分原理提出二维扩 给卅定解条件 散方程的离散化格式，并给出相应数值结果 “(x，y．：)}，一o=“o(m，y，：)，(z，y，i)∈力， (3) 一警)伸c】。 (4) 这里an足区域力的边界，n是an的外法线方向 单位向量，u。和-厂为给定的函数，。，p表示不同类 型边界条件的参量．为简单起见，下面只讨论第二类 边界条件，即a=1，口=0的情形，这时边界条件(4) l 2 可写成形式 罔l非正规六面体嘲格 w·H f0．(5) IⅫ=厂(x，y，=，￡)l。n， Gnd l㈣ahedmn №．J