二阶STC格式及对跨音速湍流流动的数值计算.pdf 6页VIP

中田工程热物理学会 编号 赫靠气钴羹力学学术会议 二阶S!g垫塞及其对鬯童鎏谴遒萋壁的爨堡过蔓1 枥吉民徐建中黄修乾 (中田科学藏工程搏暂理研究所．北京10∞∞) ● ■■吼格式悬一种采■积分骺式物理守恒棒的计算方法，^有守恒性好．精度高的的特点．目曾求■ E岫方程的一卅和二阶sTc格式已垤建立起来．乜求算N．s方程的s1℃格式尚在发晨之中．为了推广slc ● 格式在采解齄性蠢动时的应用，奉文在曲域坐标摹下构造出了柬■二健H．s方程的二阶sTc格式．井用谈格 ● 式对辟音蕾■童蠢动闩矗进行了计算． ● 1 引言 空间一时问守恒Csl℃)格式ll川是在cE／sE方法Hl的基础上建立起来的一种求解积 分形式的钧理守恒律的计算方法．与其它守恒型有限差分格式及有限体积方法相比．它 具有优良的局部和全局守恒性能．文献【I，2】中已经建立起了求解Elller方程的STc格式， 并对多个流磅问曩进行了计算，计算表明s’rc格式是一种具有较高精度的计算方法，能 够根有效地处理古有激波、接触间断等复杂流动现象的流动问题．与s1℃格式用于求解 EIller方程的发展相比．用以求解N．s方程的s1℃格式还不成熟，尚在继续发展之中． 为了推广STc格式在求解粘性流动时的应用．本文用二阶TayIor展开式来表示流动变量 在解元中的分布，建立起了求解曲线坐标系下的二维N．s方程的二阶sT℃格式。用建立 起来的格式对跨音速滴流流动问题进行了计算． 2二维N．S方程的二阶STC格式 在曲线坐标系下，二维可压缩N．S方程的无量纲微分形式为： 堑§生：上f堕+鲴 (1) 西必卸鼬。L必卸J 一 其中 ■ ‘ ’ 、 1 11 ● J j j J j 垂=号．童=号占+号，．声=号E+号，，豆=号E+等E．丘=号E+号E∞】’ j’ j j’ 上式中的Q、E、F、E和只的表达式可见文献(3】，J为： ．，=号，T1，一号，q， (3) 令z1={，屯=11和x，=f为三维Euclid空间毛中的坐标。在毛空间中，根据GauSs 公式，得到方程(1)对应的的积分形式为： ‘本文得到田家重点基础研究发展规划项目课题(No．G1999022305)和国家自然科学基金重大项目(No． 59呻54“)的经费资助． 柏 护’出=u (4) 缸r， 其中，(1)双矿)为E，空间中任一空间一时间区域y的边界表面，(2)西=dⅫ，db为 烈矿)上的微元表面的面积．H为指向微元表面外侧的单位法向量，(3)日为E，空间中 的通量密度向量： ● 肚l豇击豆芦击觚．j ‘5) 通过方程(4)，空间坐标和时间坐标被紧密地结合起来，同样地处理。当“矿)只是 一个网格单元的边界表面时．(4)式给出流通量的局部守恒关系．当s(y)包含许多网格 f 单元乃至整个求解域时．(4)式则给出流通量的整体守恒关系。方程(4)是离散求解的对象。 在E，空间中，网格点在({，11)平面上的布置如图l所示．空间步长取为：△芎=Aq=l。 0 随着时间的推进．以·代表的网格点和以O代表的阿格点是交错排列的。 1r 、 1 ， 1 丫 ÷ I ’ f： ： i i ’ 、 l i’。’ i。 圈I网格点在(毛，T1)平面上的布置 解元疆(，．)和守恒元C已(P)的构造分 别如图2和3所示。在解元中．假定流动变量 是连续的，不出现流动间断，而在守恒元中可 阻有流动间断的存在．这样．在解元中可以将 流动变量按T’aylor展开式展开来表示其分布． 知 STc格式的精度与所用的T可10r展开式的阶 羲有关，用高阶的展开式可以建立出具有高精 ‘ 度的格式，以下采用二阶展开式来构造格式。 图3守恒元CE(尸’) 在解元舾(，’)内，以P’U，．i}，，h址)为基点，将任意点({，nf)处的互、豆和户按二阶 T对lor展开式展开为： 垂({，唯，)，=釜，+(磊)，砖+(磊)，卸+(磊)，6r+÷(蟊)，6E2+÷(豆。)，却2 +{(磊，)，2+(磊。)，．砖卸+(毛)，髦6f+(磊，)，矶6f 富({，1，，)，=富，+(丘)，．鸶+(言。)，砷+(E)，+{(毒E)，5E2+；(营。。)，矶2 +{(爵)，6f2+(丘。)，．6号姻+(露)，毽5f+(瓦)，两6r 41 威％，n，f)，=戽．十(最)，砖+皈)，如+晖)，+÷(元)，6{2+{(露)，．两2 ‘ ‘ (8) +；(只)，甜2+(气)，砖6n+(元)一6秘+(瓦)，卸6r 式中，砖={一{，·矾=11一”P，5r=r一0· 在解元船(，’)上，任意点({，1，f)处的通量密度向量口为： ● 聪)，=【凤，咖·一击豆㈣，r)，一，户(w)，一击丘㈣，『)，，亟(㈨’J(9) 由于流动在解元驻(P．)中是连续的，所以微分方程(1)在解元中成立。将式(6卜(8)1 代入方程(1)中，在网格点P处得到如下关系： ， (10) ● (磊)，=一(巨)，一(E)，+击((豆)t】，+击【(丘)r】一 (磊沁一缸卜(最．¨击㈣氍¨奇‘丘)“】， (11) (12) (毛)一=．(茸．)，一(瓦)r+i}【(官-)nb+i}【(置)nnb (13) (区)，。-(毛)r(‘)，+击f(言r)口】，+击【(丘)w】， 记脚；抛y蚀，Ⅳ=a可勉，庸=§，村+{，Ⅳ，霄=q，村+n，Ⅳ·则由此可推得言、 户的各阶导数与童的相关导数之间的关系，通过联台式(10H13)分析知，若(Q)，．、 (g)，、(岛)，．、(Q。)，、(如)，及(如)，的值已知·则其它的量就可根据它们求出来。 所以这里将(西，、(磊)，、(磊)，、(磊)，，(瓦)，和(嚣)，作为不同的未知量同时 进行求解，这是与传统的数值方法不同的地方． 将方程(4)应用到守恒元CE(P)中．这时S(矿)为守恒元CE(P)的边界表面·则有： 俨协=o (14) 在守恒元CF(，’)的各边界表面上，利用各相应解元上的式(7H8)，对上式积分得： ● 式中， Ⅲ， 百=船，地+亟(一*o]，+垂(一*0]。+耍(抖0)。 云=吱。毒书．一瓦。，；，钧，一言(吁÷，钧。+童(o，一；，第。 ct·， F=文扣书。+d一；，o，訇。一户(一；，o，书，一d扣书。 c㈣ (19) 瓦；壶呱)．+曦)r+(西)e+哦)。+(瓦)一+瓯)r+(氐)e+(∞。】 (21) 巨=或(o’；。等)。+茜(o，一；，等)。一丘(0，；，警)。一置(0，一；，等)。 (22) 丘=晟t{0．争。+丘(一；，o，缸一丘({，0，瓤一只(一扣争。 (23) ● (磊)，和匝．)，．的计算采用加权平均的方法m司，即： ●● (耍)矿卟磊)；，匣)；，伍】I浅)，=咋磊)硝磊)；，q】 (24) 式中．(磊)；．、(区)；．、(磊)；和砭)；的计算式分别为： 亡 c磊肛一陬o，o，钧。+币to，訇。一z西] czs， ，醢)；H．=硝0，0，铜+爿o。0，譬)一2西 (26) 、 ‘，， 、 ‘／o c磊，；一陬o，o，钧．+域o，o，钧。吨玉] c：，， (磊臻=攻0Io，笃。+攻0．o。匀。一珥 (28) 这里用加权平均方法韵目的是为了抑制通常出现在漉动阔断附近的数值解的撮蒋， 使格式艟够处理吉有间断的流动． 二阶导致(垂毫)，、(瓦)，．及(垂。)，悬用点』、c、￡和G处的垂E、磊的值通过 中心差分来计算的．即： (菇)，：堕坚丛尝幽 (29) (k)，；堕k鲥掣址丝丝 (30) ● 1) ● (瓦)矿堕k型掣垒生二塑世 (3 3数值计算 下面用上述sTC格式对绕表面带有圆弧形鼓包的圆柱的跨音速特流流动流动问题 进行计算。流动结构和圆弧形鼓包的几何形状如图4所示。圆柱的直径是15．2cm，圆柱 的前端到圆弧形鼓包的前缘的距离是61cm．鼓包的高度是1．9cm．弦长是20．3cm。 对马赫数为M。=0．875、0．925两种来流情况进行了计算．阻鼓包的弦长作为特征 长度，来流的单位长度的雷诺数是RⅣm=13-6×106／m。计算中将坐标原点放在鼓包的起 点·计算域为：在流动方向从_6lcm到81．2啪，垂直于流动方向为从圆柱体壁面到距壁 面32cm。初始条件为：在初始时刻整个流场是均匀的，把自由来流曲参数指定为整个流 ．。。．。。．．。。．．．．．．．L 场的参数，参变量的导数均指定为零。边界条件为：左边界为自由流边界条件．给定总 温：总压和来马赫数：在右边解指定背压i上边界采用自由流边界条件i下边界是固体 壁面，采用无滑移边界条件，壁面温度保持恒定。湍流模型采用B．L代数模型。网格点 数为201×81．d呓=n9，a=l。由于计算机条件的限制，壁面附近的网格没有取得太 密．距离壁面的第一个网格点没有布置在牯性底层，而是布置在对数律层，它与壁面间 的距离约为平均网格间距的千分之四，换算得到的y+约为66。 图5是计算出的马赫敷等值线，图丘是壁面 ● 处的压力分布．在两种来流情况下，流动都在鼓 包上产生激波，在激波后边界层出现分离。在 肘一。=0上75时，激波的位置在x，c兰06l处，边界 ： 层分离出现在j尼兰0．70处，计算的最大马赫数为 村．。=131．在札。=0．925时，激波的位置在 ● r，c兰¨3处，边界层分离出现在z，c皇0．72处，最 大马赫救为肘一=138，流动中有超音速舌形区的 存在．在札。=0盯5时的实验结果是”l：激波的 位置在x，c兰％3处．流动分离出现在x，o兰n70 处，最大马赫教为虬。=l_，2．箍着来流马赫数 豳7马赫数等值线州 的增大，澈波的弯曲程度也增大了．图7是文献【91对^t。=0．92S时的流动计算得到的马 赫教等值线．从各方面的比较可以看到．本文的格式计算得到的结果还是比较好的。 4结论 本文为了推广STc格式在求解粘性流动时的应用．采用二阶T．aylor展开式来表示流 动变量在解元中的分布，构造出了求解曲线坐标系下的二维N-S方程的二阶ST℃格式。 用建立起来的格式对绕表面带有圆弧形鼓包的圆柱的跨音速湍流流动问题进行了计算。 计算结果表明，本文的二阶sTc格式具有较高的精度，瞻够有效地处理含有激波、边界 层分离等流动现象的流动问题． 参 考 文 献 矗 …刘湃涛，镣建中求■Eukr方程的空问一时闻守恒格式工程热物理学报，1997，18【3)：2舛～299 ● H．T J zxu·叫x con矧vm啪scheme Lju s呻oe．Timc in№n叭boI帅柚c山Yilinc盯c∞rdin批cljmo啦 【2】 QH啪B J．0f^a口n-叫苗。19蛆．11(3)：165～176 【3l黄售乾，按建中．用s1℃格式末解二雄激菠一边界层相互柞用砑置．工程垄}I物理学报，2000．2I(1)：23～28 №vic卜s咄酋勖un‘哪J．咖口IlLH舭，l钟5，119：2蜘监4 Y Y．a·w-1．ds EI削啪t¨d w帅‘c sp∞e-Thc帅鹫r”don 【却x cch¨‘州曲R瑚lub帅Ed盯so～∞B^蝴on SdI咖日曲螂nMe【Ⅲ^I^^P■畔r96·d764。19％ r YL0f．柚d E Yu，X P C fq^H抽也s．C．‰‘Stw嵋‘C J郴；帅H啦R啪I州∞Ge呻眦峥 M山lidi咻吐i帅Il Ehle吡枷Sduuon日ementMtl|lod SoIuⅡ加0fC呷oHⅧi帅LI懈by山eSp丑ce-TIn坼C帅鲫n呻n ^l从hD甘∞-10帖．19鲫 YChow11lc soIu‘ionElemenlMeth。d．AHew sC 4啊C 【，l a瑚％xY．w蛐B S畔}Th帖Con，叮v帅帅日em咖●IId c。蛳m师日h惦J c鲫puPIlys】999．1始 HjIh-R酬岫∞鲥∞u№睁^l业讨jm酬j。枷P啪dj胛衙酬vjnl 棒136 ¨ ——事■—__—飘{r B·叫o PI叫啪向r‘姗出