高阶多矩算法及其在全球模式中应用.pdfVIP

第十五届全国流体力学数值方法研讨会 北京香山201[年8月23-25日 目 目 膏阶多矩算法及其在全球黼的应用 Ⅷ㈣自{ⅡRm☆㈣暑毳^女％№，100Ⅲ IIM：肚1虻V*想，m$”究了约束矩的t择#目地自自度的定义．提mT一种新的≈m$ 矩掉让。＆碇T《“^《求解模型，＆用TvBm制#上除TM断R域的非物Ⅲ振荡．目H保jjf_r 光mK域的模拟精＆n“Ⅱ方球日格r构建r高阶牟球多《拽水披幢型．为±球环流模式垃《 r新的数位框架。 关●I■·高阶算沾．多如＆多矩∞柬．有R体${法．《＆^程，±球横“ 基于多矩思想发展的多矩约束型有限体积方法(Multi—momefltcon$【r￡incdnⅡ11e volame method．MCVmethod)¨1为发胜任意高阶精度的计算格式提供了统框架。奉 研究进一步时约束锋的选择和埔地F|由度的定义进行了研究。通过在H算单元内部定 义新的约束矩，避免了鞭掉法单纯依籁单元边界岛阶导数矩约束的局限性。该方法一 方【自】减少，单兀边界愤息交换。娃著提高r此类数值模直在进行火规模并{I运算时的 通讯效半，另 方解决了拒阐格不连续处可能m现的计算小稳定问题。研究r心地自 由度定义位置对H掉稳定性厦}_算精度的影响。相比等距定义方式，按GuassLobaao 点酣置计算荦，o内韶一m度厦约来条件捉商丁葬沽竹性能。本研究柑建了j至卉阶 MCV格式。通过数值试鸵椅验丁尊法收敛性(如表l所小)。 表l}Ⅻ∞求窄弁m体${方法收敛性梭骑 —丽F—r—～—1；———一一—]r]i—一 (aepl——再i——面五■—1而r—1；aF一面；i—ir面 ————可—了菇磊i——■——亍：百丽酥丁_两醑ir一_ fil 4 2 ． ． 08ll“ 2惦73“ 21064e-5 · — f30015 7∞∞-7 7j518e-7 73983e-7 ． ． 6 26l弱“ 2617口e-8 2猖150-8 4 l 398 3∞ 121 3155em 131聃邮 l3175e-64帅 f6帅1 5 22347e^505 22377自-85∞ 22347e_85∞ 6 4嘴1l}10600 t1149e-105鲫 41504^1D6W ———i—面if；r—ji…一nii—Ⅱr1■目—1jr 241 {s船5缸8＆卿 8街％时 ^鲫 ＆船目}8 3船 5 6 501 602120d050l 50l f12001 9188810 6。28缸10 6 644d3e-125蛆 64弱仳125鲫 64T0缸12600 本研咒笈展了求解欧拉疗程的多矩约束型有限体积模犁。使蹦TvB限制器[2]保证 了数值葬法住求解问断问题时十m现1F物理振荡现象．『d时可以保证光滑K域的H掉 精度。利用多矩约胤柯限体积方法11_算sod问题及激波／湍流相n作用问题的鲒粜如图 1 2所示。TVB限制器的缺点足存在需人Tr预的参数．且随问趟小ld咖变化。时r求解 川断问题．F步的研究方向足使用WENO思想【3】构造更高性能的限制器。