求解对流问题和对流占优扩散问题的值域离散网格方法

来源 :中国科学技术大学 | 被引量 : 0次 | 上传用户:marker1900
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对流问题和对流占优扩散问题存在于诸多自然界与工程实际问题中。这类流动问题的主要特点是:流场内会存在物理量在很小空间尺度内发生剧烈变化的锋面,甚至激波;而在远离锋面或者激波的区域内物理量变化相对平缓。在求解这种存在多尺度结构的流场时,传统数值方法往往会遭遇很多困难:使用常规有限差分格式进行求解时可能会导致物理上不合理的数值震荡,甚至引起稳定性问题;工程上常用的一阶格式往往会在锋面或者激波处产生较大的数值耗散。为减小数值耗散,提高求解分辨率,锋面或者激波处的网格尺度必须足够细分。若采用全场均匀的细网格,则会耗费大量的计算资源,因此采用自适应网格是一种现实可行的方法。然而,传统的自适应网格法往往基于空间离散网格,通常需要复杂的自适应操作以及大量的网格重构过程,这在一定程度上增加了计算复杂度。因此,如何利用有限的计算资源准确模拟锋面或者激波的形态及其附近流场的变化仍然面临巨大的挑战。本文注意到传统数值方法往往基于空间离散网格,在采用空间离散网格求解对流和对流占优问题时遇到困难的一种解释是:空间离散网格使用有限小的网格单元描述零宽度的激波或者接近零宽度的锋面。鉴于此,本文提出了一种新的离散网格系统——值域离散网格。顾名思义,与空间离散网格将流动问题的位域空间离散成若干网格节点不同,值域离散网格是通过将流动问题中物理量的值域空间离散化而得到的求解网格。本文发现,与均匀的空间离散网格相比,值域离散网格能够锐利地描述处于任意位置的间断,这对于求解存在激波的问题是重要的;另外,值域离散网格节点在锋面或激波处聚集,而在物理量变化较为平缓的区域内分布较为稀疏,这意味着,该网格具有与生俱来的自适应特性,有利于在消耗较少计算资源的前提下获得较高的求解精度。与传统自适应网格相比,值域离散网格具有简洁的数据结构;并且,相比于传统自适应网格法须在每个计算步进行网格重构以达到自适应的目的,值域离散网格的自适应过程是随着流动问题的求解而自动完成的。因此,如果将值域离散网格应用到对流问题和对流占优扩散问题的求解中,则有望获得一种计算速度快、计算精度高的数值方法。基于值域离散网格的特点,在求解流动问题的过程中,并非求解空间离散点上的物理量值的变化,而是追踪具有已知离散物理量值的网格节点的运动。然而,在实施值域离散网格法的过程中需要解决一些复杂的数学物理问题。首先,本文将值域离散网格应用到一维空间中可能会产生激波的对流问题的求解中。由于值域离散网格所描述的激波强度只能是离散的,如果采用Rankine-Hugonoit激波关系计算离散激波的运动速度,那么可能会导致守恒性问题。为了解决这一问题,本文通过守恒关系推导出了修正的Rankine-Hugonoit公式。分析发现,利用该修正公式不但在求解过程中满足了守恒性要求,而且激波位置的求解能够达到二阶精度。另外,相邻网格节点可能会在求解过程中互相发生位置逾越而导致物理上不合理的多值性解。为了解决这一问题,本文设计了时间步长的选取规则,并引入了基于熵条件和守恒性的后处理步。如此,不但多值性问题得以避免,从而可以模拟出激波的产生、演化和衍灭等复杂过程;更重要的是,该方法的时间步长不受CFL条件的限制。其次,本文提出了求解一维对流扩散问题的值域离散网格法。由于扩散项的存在,网格节点的推进速度不仅与当地物理量值有关,而且还与其周围物理量分布形态有关。如何确定网格节点的运动速度是将值域离散网格应用到求解对流扩散问题的关键。本文基于值域离散网格建立了在值域空间内固定,而在位域空间内动态变化的控制体;通过控制容积积分法得到了流动方程在值域离散网格上的守恒型离散格式。本文通过在极值点处设定的网格节点及其控制体,并利用抛物线线型逼近给出了离散格式的调整,从而将值域离散网格法推广到非单调问题的求解中。对于初始值中存在的常值分布段的流动问题,本文将常值分布的两端视作移动边界条件,并给出了移动边界的高精度处理方法。本文将值域离散网格法推广到二维两相渗流问题的求解中。根据Buckley-Leverett理论,两相渗流问题的压力方程为类拉普拉斯方程,而饱和度方程为非线性双曲型方程。鉴于有限分析格式能够高效求解类拉普拉斯方程,而值域离散网格法能够零耗散地求解双曲型方程,本文将二者相结合,得到了一种求解速度快、精度高的顺序求解方法。本文将该方法应用到两相驱替问题的研究中,通过对比不同流度比下两相界面形态,探讨了粘性指进现象和网格取向效应的内在联系。综上所述,本文针对采用值域离散网格法求解对流问题和对流占优扩散问题时所遇到的数学物理问题进行了初步的研究。设计并初步实现了求解这些流动问题的值域离散网格法,通过数值算例证实了值域离散网格法的可靠性和有效性,并利用所提数值方法研究了两相渗流问题中有趣的驱替现象。本文所提方法有望为对流问题和对流占优问题的研究提供一种快速的、精确的数值工具。
其他文献
700 t双体起重工作船可将修理的船舶尾梢抬起,以便进行水下工程检修和保养。其船体结构采用双体和局部连接桥结构。两个片体在尾部采用连接桥连接而成,连接桥在航行环境中要
以酿酒过程中产生的酒糟为原料,制备得到氨化酒糟吸附剂,用于去除电镀废水中的六价铬〔Cr(Ⅵ)〕,利用FTIR、EDS、BET和SEM等对吸附剂组成和结构进行分析和表征,在不同的pH及
<正> 为了提高家蚕选育工作的效果,必须提高选择差数(被选择的蚕茧后代的指标比平均值群体的优越)。根据多年试验结果,我们研究和介绍提高选择差数的方法子生产,包括多次交配
期刊
行政审批制度改革作为政府转变职能的突破口,在我国政府行政体制改革过程中发挥了重要作用。十八届三中全会之后,我国的行政审批制度改革进入全面深化改革阶段,行政审批制度
对于如何培养孩子们的创造能力,主要是以孩子们的兴趣为主,激励孩子们经常发散思维,大胆想象,不要让孩子们形成思维定势,经常举办美术活动,丰富孩子们的知识,让孩子们自己思
有效教学理念的核心就是教学中的效益问题。教师如何在一段时间的教学之后使学生获得具体的进步或发展? 本文根据自己的教学实践,结合情意原理和序进原理探讨化学教学中的有
主要阐述了油轮单壳改双壳工程中的重点和难点一变形控制,介绍了在工程实施过程中采取分期分阶段施工的方法,合理安排现场工艺孔的开设,并加强施工过程中的现场监控和测量,使整体
根据2000—2018年中国大气污染治理技术专利数据,采用专利分析法,从产出趋势、领域分布、空间分布、申请人等方面研究中国大气环保产业技术创新能力。结果表明:中国大气环保
加强生物实验教学培养学生实验能力山东费县第一中学(273400)王志鹏生物学是以实验为基础的学科,实验教学是生物教学的一个重要组成部分。通过实验教学,让学生掌握一定的实验技能,学习