通用壁面边界条件.docx
- 文档编号:26260095
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:9
- 大小:191.59KB
通用壁面边界条件.docx
《通用壁面边界条件.docx》由会员分享,可在线阅读,更多相关《通用壁面边界条件.docx(9页珍藏版)》请在冰豆网上搜索。
通用壁面边界条件
通用壁面边界条件
1、简介—目的
目前,科学商业软件CFD的主要缺陷之一是由于壁面边界条件有效设定而引起的。
通常情况下,用户可在两种标志方法之间进行选择。
第一种:
高雷诺数Re边界条件—假定对于近壁处网格使用壁面对数律算法;第二种:
低雷诺数Re边界条件—可以解决整个边界层包括网格质量很好的内层。
不幸地,两种方法都需要依赖于流体流动特定的网格特性。
结果,用户在创建网格时需要预先了解流动的各种细节。
为了克服这一问题,在这里提出了通用壁面边界条件这一基本构想。
这种新方法是以高雷诺数Re边界条件
D
为基础,通过与低雷诺数边界条件区域混合而延伸来的,这一区域通常看作是当地网格分辨率的距离小到接近壁面的距离。
所提出的混合方法允许在没有任何用户输入的情况下来适应它的当地的流态。
因此,新的通用壁面边界条件没有近壁网格分辨率的具体要求,它适用于二维流动也适用于三维。
作为新方法的承诺是独立网格简化了网格生成过程,这一方法也可能是工业流领域中对CFD应用的一个重要的改进。
Grotjans和Menter博士基于壁面函数形式提供了类似的方法,这种壁面函数形式提供了通用网格,但是,它总是与咼雷诺数形式一样产生相同的缺点。
新方法的实施和测试在两个简单的二
维流动中进行:
一种是Kim的槽道湍流,一种是在大升力条件下的ONERA-A的翼型绕
流。
此外,这种方法已经应用于其它翼型绕流,比如:
Devenport和Simpson的wall-mountedwing、
2、方法
新方法提供了近壁网格的平均动量和湍流特性条件。
介绍高雷诺数边界条件之间的混合函数时,基于到壁面的无量纲距离:
■U■
丫high=~E-及Y|ow=U
为了得到丫+表达式,对e"/用泰勒级数展开可得:
-f
1I:
二u
=Ue-VU
EII
对小参数来说,这种扩展的收敛速度很快然而,与eU项相比,对大参数而言收敛速
度就非常的小了。
正确的方法是去掉九阶项之后各项。
1)动量方程边界条件
解动量方程时需要壁面切应力,为了确定这一特性,引入一种迭代算法,这种算法要求对每一个点进行逐个迭代求解。
在下一
次迭代中壁面切应力或切应力速度必须被使用。
除了混合速度U+外,所有湍流量都用这种相似的方式进行处理。
2)湍流特性边界条件
对于低雷诺数和高雷诺数两种形式来说,相同黎曼边界条件下湍动能都是可以被使用的,此条件仍然为:
=0
icn丿壁面
此外,湍动能的产生受到近壁处控制体
积的影响,米用标号P表示为:
fE2、
pPp=
Ikn丿p减少到湍流长度尺度边界条件是第二湍流数量的边界条件,这一边界条件被应用与耗散速率;方程或湍流频率•方程中。
Wolfshtein制定了近似的近壁的长度尺度公式:
1
2
Ci
公式是以湍流雷诺数Re为基础。
参数
=0.2不同于Wolfshtein给出的值(=0.263),这是为了确保从壁面函数到低雷诺数形式相融合的正确性。
到壁面的湍流长度尺度的渐进式的演化研究得出了湍流长度尺度的表达式。
通过用泰勒级数展开的线性截断表达的指数函数可以在近壁处控制体积内进行分析。
对于湍流频率的边界条件可以用类似的方法进行导出。
3)涡粘性't边界条件
用同样的方法可以给出Spalart-Allmaras一方程模型的通用边界条件。
结合这种新的
边界条件,建议使用Spalart-Allmaras模型的Edwards修正,这是由于在近壁处区域可以提高它的稳定性。
就Y400和涡粘性非常大情况下而言,问题的产生受到了一定的限制。
更多的详细信息在可以参看文献[9]和
[13]。
3、范例
1)槽道湍流流动
对于一个基本的确定的新的边界条件在简单的槽道湍流流动的应用。
在这种情况下,在近壁网格求解时,改变特定网格处的流入雷诺数,则计算结果也会轻微变动。
大的雷诺数将导致无量纲距离的增加。
将这种新的混合边界条件应用到槽道流动中,可以预测雷诺数的变化范围内的速度的精准变化。
图1展示了用Wilcox的-模型求得的计算结果。
近壁层从丫=1,
Re=3640到丫十>10,Re=56000时仍然具有
相同的速度
“口
」』.■・LH'.j
110
图1由混合壁面边界条件和Wilcox的-模型获得
的不同雷诺数时槽道湍流流动的无量纲速度
00
100
nn
DOJq_3
ID11011D2ID
0.2
n
200
图2由混合壁面边界条件和Wilcox的-模型获
得的在不同雷诺数下的槽道湍流流动的湍动能和切
应力
在图2中,速度、湍动能和切应力都是
由不同雷诺数和不同网格独立求得的。
结果证明,至少对这种简单流动,应用新的通用的边界条件,网格的影响完全消失了。
尽管如此,通过与Kim、Moin和Moser三人的直接数值求解(DNS)比较,的不足之处还没有得到改善。
'-;模型和Spalart-Allmaras一方程模型能够提供相似的图形。
这些图形
和图2中绘制的基本相同。
2)高升力翼型流动
在以下确定的试验测试中,在其关键位置的流动中,新边界条件和数值网格对计算结果的影响需要另外研究。
在Re=2106,
Ma=0.15和攻角为13.3时的绕ONERA-A翼型流动时,其特征为在翼型尾迹区域带有湍流再附着点的层流分离泡和湍流边界层分离。
由于这种试验测试模拟对湍流网格的细节是非常敏感的,因此对新边界条件的测试考虑这些是合适的。
在Spalart-Allmaras
模型集中的研究表明,对于这种流动产生的结果是令人满意的。
转捩点定在翼型上表面X0.12处和翼型下表面X/C=0.30处,根据以往的研究,
在第一个控制体到壁面的最大的距离
Y1.3处,网格数为22000时可以被接受的。
这些网格结果与修正网格获得的结果进行比较,用无量纲距离绘制出的图形如图3所
示,在丫+=2到丫+=20的范围内有相同的网格数。
图3ONERAA翼型在高雷诺数网格(HG)下的近壁处丫•分布
Experiment
Low-Re(LG)
Low-Re(HG)
Hybrid(HG)
Uni.High-Re(HG)
CL
1.52
1.545
1.483
1.606
1.645
CD
0.024
0.0197
0.0169
0.0146
0.0140
表1Re=2.106,=13.3处ONERAA翼型
的升力和阻力系数的测量值与计算值
对于不同的边界条件得到的两种网格计算结果如表1所示,在普通边界条件和低雷诺数边界条件下的LG网格结果完全一致,因此,仅仅在表1中展示了后者。
与质量高的LG网格结果相比,不同边界条件下的HG网格的压力分布与表面摩擦系数如图4和图
5所示。
表面摩擦系数的结果与试验采集的数据有很好的吻合程度以及合适的尺度。
正如预测压力的期望值导致粗网格HG
大幅度的分离。
与标准壁面函数形式相比,通用边界(HYB展示了很好的结果,但是它在前缘区超出了预期的结果。
根据具体的丫分布,通用边界条件和高雷诺数Re边界条件两者之间的不同仅仅出现在X/C=0.4的下游区域。
对于HG网格采用低雷诺数边界条件所获得的结果需谨慎处理。
速度剖面和边界层厚度似乎优越其它结果,但在转捩区的表面摩擦系数低于预测值,且对整个流场而言,表面摩擦系数非常接近流动分离,这在前缘区域产生很坏的结果。
这种新的混合边界条件可以提供更加合理的结果,至少在特定的
流动区域可以提高计算结果的稳定性
图4ONERAA翼型压力系数
-2.0
0.00.20.40.60.81.0
X/C
图5ONERAA翼型表面摩擦系数
参考文献
1.Cazalbou,J.B,Spalart,.RP.andBradshaw,.P"Onthebehavioroftwo-equationmodelsattheedgeofaturbulentregion".Phys.Fluids,Vol.6,No5,May1994,pp.1797-1804.
2.Edwards,J.R."Comparisonofeddyviscosity-transportturbulencemodelsforthreedimensional,shock-separatedflowfields".AIAAJournalVol.34,No4,1996,pp.756-763.
3.Devenport,W.J.andSimpson,R.L."Time-dependentandtime-averagedturbulencestructurenearthenoseofawing-bodyjunction".J.FluidMech.Vol.210,1990,pp.23-55.
4.Grotjans,H.andMenter,F."WallfunctionsforgeneralCFDcodes".InProc.of4thEuropeanComp.FluidDynamicsconference,Athens,Greece,1998,pp.1112-1117.
5.Haase,W.,Bradsma,F.,Elsholz,E.,Leschziner,.A.andSchwamborn,D.(Hrsg.)."EUROVAL?
EuropeanInitiativeonValidationofCFDCodesBand42",inNotesonNumericalFluidMechanics.Vieweg,Braunschweig/Wiesbaden.1993.
6.Kim,J.,Moin,P.andMoser,R."TurbulencestatisticsofafullydevelopedchannelflowatlowReynoldsnumber".J.Fluid.Mech.Vol.177,1987,pp.133-166.
7.Lubcke,H.M.
"EntwicklungexpliziterDarstellungenzweiterstatistischerMomentezurnumerischenSimulationturbulenterStr?
mungen.Disseration,TechnischeUniversit?
tBerlin,2001.
8.Catalano,P.,Dol,H.,Elsholz,E.,Gleyzes,C.,O'Brien,P.,Schatz,M.,Tourette,L.andWallin,S."RANSsimulationsofthehighly3Dflowoverasweptwingmountedinawindtunnel".CEASAerospaceAerodynamicsResearchConference,Cambridge,UK,10-13June,2002.
9.Rung,T."FormulierunguniversellerWandbedingungenfur
Transportgleichungsturbulenzmodelle"InstitutsberichtNr.02/99,HFI,1999.
10.Spalart,P.R.andAllmaras,S.R."Aone-equationturbulencemodelforaerodynamicflows".AIAAPaper92-0439,AmericanInst.ofAeronauticsandAstronautics,1992.
11.Wallin,S.andJohansson,A."AnexplicitalgebraicReynoldsstressmodelforincompressibleandcompressibleflows".J.FluidMech.Vol.403,2000,pp.89-132.
12.Wilcox,D.C."Reassessmentofthescaledeterminingequationforadvancedturbulencemodels".AIAAJournal,Vol.26,No11,1988,pp.1299-1310.
13.Rung,T.,Lubcke,H.M.andTh.i"eUlen,iFversalwall-boundaryconditionsforturbulence-transportmodels".Zeitschriftf黵angewandteMathematikundMechanik,2000,pp.1756-1758,
14.Wolfshtein,M.W."Thevelocityandtemperaturedistributioninone-dimensionalflowwithturbulenceaugmentationandpressuregradient".Int.F.HeatMassTransfer,Vol.12,1969,pp.301-318.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 通用 边界条件