求解对流扩散方程和涡量流函数方程的改进谱方法

来源 :山东大学 | 被引量 : 0次 | 上传用户:yan2541023
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对流扩散方程和涡量流函数方程是流体力学中常用的两类方程,对流扩散方程还广泛应用于传热学、化学和环境科学等领域,涡量流函数方程在研究粘性不可压缩流体的流动时起着重要作用。在实际问题中由于求解区域和边界条件的复杂性,很难求得这两类方程的解析解,数值求解便成为一种重要的方法,因此构造高效和稳定的数值方法有着重要的应用价值。基于谱方法的思想和积分方程理论,本文提出了一种改进的谱方法,用于数值求解对流扩散方程和涡量流函数方程。在应用改进谱方法求解稳态对流扩散方程时,对第一类齐次边界条件,可通过变量代换消去对流项中的一阶导数项。然后引入关于格林函数的辅助边值问题,把格林函数和未知函数用正弦级数展开,并求得格林函数表达式。再根据格林函数的性质,把对流扩散方程化为积分方程。最后根据正弦级数的正交性,把含有空间变量和级数系数的积分方程简化为只含有级数系数的积分方程组,求解该方程组就能求得原方程的级数有限和形式的近似解。应用改进谱方法求解了一维、二维稳态对流扩散问题的算例,计算结果表明该方法与有限体积法和QUICK格式相比,有着更高的精度并且不需要网格划分。与积分方程法相比,在保证同样精度的前提下,简化了推导和编程流程。对于非稳态对流扩散问题,方程中不仅有空间变量,还出现了时间变量。可用和前面介绍的稳态情形一样的方法离散空间变量,把关于空间和时间变量的积分方程简化为只关于时间变量的积分方程组。用四阶龙格-库塔法求解这一积分方程组,得到用级数有限和形式表示的原问题的近似解。最后用一维、二维非稳态对流扩散问题的算例检验了改进的谱方法,讨论了 Peclet数、截断项数、时间步长、计算时间以及不同源项对计算误差的影响。对于非线性涡量流函数方程,先引入关于格林函数的辅助初边值问题,并用格林函数的性质将涡量流函数方程转化为积分方程。然后对于周期边界条件和第一类齐次边界条件下的未知函数和格林函数,分别用傅里叶级数和正弦级数展开。再利用级数的正交性将积分方程转化成只含时间变量的积分方程组,用四阶龙格-库塔法求解该方程组,得到原方程的级数有限和形式的近似解。最后用不同边界条件下的算例检验了改进谱方法,发现其计算结果和解析解吻合得很好。
其他文献
近年来,随着东部地区基础设施建设的逐渐完善,我国地下隧道工程逐渐向中西部地形陡峭、构造发育区转移,随之而来的是更高的施工难度和更大的灾害风险。为进行隧道灾害的超前判识,就要进行超前地质预报工作,在众多超前地质预报方法中,超前钻探是隧道工程探明前方地质情况,进行超前地质预报的重要手段。然而对于超前钻探获取的数据,往往不能进行完全的利用,造成数据浪费。因此急需建立一种基于钻进测试的不良地质快速识别方法
学位
我国已是海底隧道建设大国,目前在建和规划了一批世界级的海底隧道工程。断层破碎带围岩松散、力学强度低、开挖后自稳能力不足,在开挖扰动及海水渗流作用下,海底隧道断层破碎带施工过程中围岩塌方失稳风险大。研究海底隧道断层破碎带塌方机理,阐明围岩塌方失稳灾变演化机制对于指导灾害防控具有重要意义。本文针对海底隧道断层破碎带围岩塌方失稳机制及灾变演化模式,综合采用理论分析、模型试验、有限元数值模拟等手段,开展了
学位
复合地层是隧道修建过程中经常遇到的特殊地质类型,在施工过程中极易产生围岩大变形、拱架突出、衬砌裂损、塌方等工程事故,造成工期延误、增加施工成本甚至会造成人员伤亡。因此,研究复合地层隧道开挖过程中围岩变形特征,并对其智能预测对于确保施工安全,保证施工工期尤为重要。本文依托新成昆铁路特尔莫隧道,采用理论分析、室内试验、数值模拟和机器学习方法,对特尔莫隧道非均匀变形特征及量化方法、变形影响因素和变形等级
学位
随着我国高速公路及铁路建设规模的不断扩大,桥梁基础的承载力和沉降控制等问题日益得到关注。钻孔灌注桩因其具有承载力高和沉降小等优点而被广泛应用于桥梁基础工程。为进一步提高桩基承载能力、减少桩基沉降、降低工程造价,研发了喷扩锥台压灌桩。由于喷扩锥台压灌桩是近年逐渐兴起的一种新型桩基础,在工程实践中多借鉴现有规范进行经验设计,缺少相应的理论指导和设计依据,同时,双锥台扩径体和肋板的设置使得喷扩锥台压灌桩
学位
因为地质、地形和天气等各种客观因素及其在建筑设计、建造和运营中的不同主观要求影响,在隧道建造和运营使用的过程中,产生了如衬砌出现裂纹、衬砌渗漏水和围岩空洞等各种病害。隧道在建成后的长时间运营中往往会由于病害的出现而陷入质量危机的状况,所以对于病害的隧道加以处治是十分必要的。因此本文将对在病害前提条件下,衬砌结构的稳定性及其处治对策等加以深研探讨,主要的研究内容包括:(1)对围岩空洞和FRP的研究现
学位
在全球范围内,多年冻土区和季节性冻土区分布范围非常广。而我国的寒区主要分布在东北、西北地区及青藏高原区。特别是我国西北地区,早晚温差过大,导致岩土体内部水分出现冻融循环现象,反复的冻胀力使岩体结构发生形变,而这种岩体变形不能完全恢复,反复循环进而加剧了岩石内部的胀缩、损伤和开裂。然而寒区岩土工程往往不仅要面临冻融循环问题,岩体节理内部不同类型的填充物也会对岩体稳定性有很大影响,由于填充物一般在岩体
学位
长期以来,水力压裂技术已成为油气开发中的核心技术之一,如何通过水力压裂技术进行高效开发,对于油气生产非常重要。缝洞型油气储层中存在大量分布随机的既有裂缝和洞室,这些既有裂缝和洞室对水力裂缝的扩展影响复杂。现有的研究表明,既有裂缝对水力裂缝的扩展行为有着重要影响,国内外学者对水力裂缝与既有裂缝相互作用进行了大量的试验分析与机理研究,并取得了一定的成果。然而,既有裂缝后存在的洞室对水力裂缝穿透既有裂缝
学位
介电材料在外力作用下产生的非均匀变形可以引起电响应,这种现象称为挠曲电效应。挠曲电效应存在于所有的介电材料中且具有尺寸依赖性。已有的研究表明,通过减小器件尺寸和提高应变梯度均可以明显提高器件的力电耦合效应。对于薄膜类二维材料,很容易产生较大的曲率,即应变梯度。那么,对于介电薄膜来说,由于应变梯度和极化之间的耦合作用,很容易产生很强的挠曲电效应。介电薄膜的这一特性,使之成为挠曲电应用的一种重要载体。
学位
黄泛区粉土是黄河流域特有的一种土体,其土颗粒磨圆度高,颗粒之间有较多的孔隙,缺乏粘细颗粒填充,土体不易被压实,尤其在毛细作用下,土体的水稳定性极差。利用这种粉土建造的路基在季冻环境下极容易因水分迁移发生冻胀融沉现象,对道路的安全造成严重破坏,如何控制土体内部含水量的变化是一件亟待解决的难题。本文用高吸水树脂(SAP)作为改良材料,将其掺入土体中,对改良后的土体进行了相关性质的研究,并利用室内冻融循
学位
随着高速铁路建设的发展,出现了越来越多穿越软岩地层的高铁隧道,为了保证高速列车的正常运行,隧道结构动力学问题亟待解决。一方面,当列车速度提高时,列车对隧底结构的振动冲击作用增强,影响隧底结构的稳定性。另一方面,当隧道基底为软岩时,列车振动荷载与地下水的耦合作用下将其碾磨成泥浆,影响隧底结构的基础条件。因此,开展列车荷载作用下炭质页岩高铁隧道动力响应特征及稳定性控制研究,具有很强的现实意义。论文依托
学位