基于色散关系分析的高精度气动声学计算方法研究

来源 :湖南大学 | 被引量 : 1次 | 上传用户:xuxiaoxiu
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
气动声学问题与流体力学问题有着本质的区别,直接数值模拟气动声的产生与传播问题时,对数值算法的要求远远高于计算流体力学。为了正确捕捉波的色散性、耗散性、相位速度与群速度等传播特征,数值离散格式的色散和耗散误差必须尽可能的小,同时,在数值处理开放边界或者半开放边界问题时无明显反射。因此,发展高分辨率的数值离散格式与高性能的无反射边界条件是计算气动声学(Computational Aeroacoustics,简称CAA) CAA的两个关键问题。本文根据CAA的高精度技术需求,针对这两个关键问题,从数值格式与控制方程的色散关系(频率与波数关系)出发,研究了高精度离散格式的误差行为、高精度空间数值离散格式与高精度无反射边界条件。主要包括:高精度有限差分格式的误差分析,显式有限差分格式的保色散优化,非均匀网格的高精度紧致差分格式,适用斜流计算的线性与非线性完全耦合层(Perfectly Matched Layer,简称PML)吸收边界条件。采用国际标准CAA算例验证、测试与分析各格式及吸收边界条件的求解性能。本论文的主要研究内容如下:(1)提出了一种构造高精度色散保持显式有限差分格式的优化方法。为了获得高精度的数值离散格式,提高离散格式数值模拟长距离、长时间传播波的能力。首先,通过傅立叶与拉普拉斯变换法则将有限差分格式变换至频率与波数空间;然后,采用修正误差传播分析方法和Neumann误差分析方法,从理论上详细推导了高精度有限差分格式误差构成、误差传播以及误差累积的数学公式,为发展低耗散低色散高精度数值格式提供理论依据;在已有DRP优化思想的基础上,重点考虑群速度对格式色散性能的影响,提出了一种新的DRP优化方法。该优化方法以绝对误差最小积分值为目标函数,格式系数与积分上限为优化变量,色散、耗散及群速度误差为约束条件,采用序列二次规划方法,通过迭代计算,寻找最优解,获得高精度有限差分格式系数及积分上限值。基于该优化方法,优化了7点、9点及11点有限差分格式,并分析其色散、耗散与群速度性能,最后通过一些典型算例验证误差分析及优化方法所得的结论。(2).提出了一种构造非均匀网格高精度紧致差分格式的方法。针对非均匀网格容易导致紧致差分格式产生较大误差、引起虚假振荡等问题,首先,基于光滑的思想,通过泰勒展开分析方法,详细给出了内点与边界点紧致离散格式系数的通用求解表达式;然后,根据该理论公式构造了一个3对角6阶精度的紧致差分格式。通过傅立叶与拉普拉斯变换法则将该紧致差分格式变换至频率与波数空间,推导出该格式数值耗散与色散的理论计算公式,并以均匀网格、扰动网格、拉伸网格及突变网格4种典型网格为例,详细研究该格式对不同网格的耗散与色散特性,重点研究了扰动因子对其的影响;采用特征值方法,推导了半离散格式满足渐近稳定的条件,并研究了不同网格下该紧致离散格式的渐近稳定性。最后,通过一些典型算例测试了该格式的求解性能,数值计算结果表明所得数值解与精确解吻合,本文所提出的高精度紧致差分格式,在进行非均匀网格气动声学数值模拟时,能够抑制因网格不均匀而引起的非物理振荡,计算准确。该格式构造方便,精度高,显示了其在非均匀网格模拟中的优越性。(3)提出了一种构造二维对流线性PML吸收边界条件的七步通用方法。首先,通过傅立叶与拉普拉斯变换法则,将对流线性控制方程转换至频率与波数空间,并建立控制方程的色散关系;然后,根据该色散关系轨迹分析流场内各物理波的群速度与相位速度方向,提出了保频型与变频型两种时间与空间坐标变换关系式:基于这两类坐标变化关系式,发展了两类二维斜流线性PML吸收边界条件,并对各层吸收边界条件的控制方程进行了理论推导;采用特征值方法,分别研究了Hu的分裂型、本文建立的保频型与变频型斜流线性PML吸收边界条件的动态稳定性,重点研究了来流马赫数、波数、吸收系数及额外增加的吸收项对其稳定性的影响;最后,通过一些典型算例测试了这些斜流线性PML吸收边界条件的数值求解精度与可靠性,数值计算结果表明:PML域内各物理波具有一致的群速度与相位速度方向仅为PML稳定的必要条件而非充分条件;Hu的分裂型斜流线性PML吸收边界条件不稳定,保频型斜流线性PML吸收边界条件条件稳定,变频型斜流线性PML吸收边界条件稳定;变频型斜流线性PML吸收边界条件对声波、涡波与熵波具有良好的吸收特性,非常适合用来计算斜流CAA问题。(4)提出了三种适用斜流计算的二维非线性PML吸收边界条件。首先,通过伪时均流假设,建立了扰动量的守恒型与非守恒型控制方程;然后,采用变频时间与空间坐标变换关系,将这些扰动量控制方程变换至新的时间与空间坐标系:采用复数变换方法,分别构建了非守恒型、守恒非分裂型与守恒分裂型斜流非线性PML吸收边界条件,并推导了各层PML吸收边界条件的理论公式;通过线性假设,将非线性问题转化为线性问题,采用特征值分析方法,对建立的两种守恒型非线性PML吸收边界条件进行稳定性分析,详细研究了来流马赫数、波数、吸收系数及增加的吸收项对其稳定性的影响;最后,通过一些典型算例测试和分析了这些斜流非线性PML吸收边界条件的数值求解性能,数值计算结果表明:当角层非线性PML吸收边界条件中未额外增加吸收项时,系统特征值的最大虚部均为正数,吸收边界条件不稳定,当角层非线性PML吸收边界条件中额外增加吸收项时,系统特征值的最大虚部均为负数,吸收边界条件稳定;非守恒型斜流非线性PML吸收边界条件的吸收效果最好,分裂型PML次之,非分裂型PML反射量最大。本文建立的非守恒型、非分裂型和分裂型PML吸收边界条件均稳定,能够有效抑制边界反射,提出的构建斜流非线性PML吸收边界条件的方法切实可行。(5)提出了适用斜流计算的三维线性与非线性PML吸收边界条件。在二维线性与非线性斜流PML吸收边界条件的基础上,将其构造方法拓展至三维斜流PML吸收边界条件的理论推导过程中。将三维斜流PML计算域分为26块,7大类。根据流场内各物理波的色散关系轨迹,提出了一组适用三维斜流PML构造的变频型时间与空间坐标变换关系式。基于这组坐标变换关系,通过复数变换法,构建了三维斜流线性、三维斜流非守恒非线性、三维斜流守恒非分裂型与分裂型非线性PML吸收边界条件,详细交代了各角层PML吸收边界条件中吸收项的施加原则。此外,还重点研究了来流马赫数、波数、吸收系数及增加的吸收项对各角层PML吸收边界条件稳定性的影响。最后,通过一些典型算例测试和分析了三维斜流线性与三维斜流守恒分裂型非线性PML吸收边界条件的吸收性能,数值计算结果表明:本文建立的三维斜流线性与非线性PML吸收边界条件对声波与涡环均具有良好的吸收性能,无明显反射现象产生,计算稳定。
其他文献
目的:探析人性化护理干预方式在手术室护理里的具体应用效果。方法选取2013年10月~12月来我院入院治疗的手术患者共80例作为此次试验研究对象,将其随机均分为对照组40例和观
欧盟以碳排放交易的方式履行《京都议定书》承诺的减排要求,这种市场机制促使生产企业减少温室气体的排放,同时加大节能技术的研发。从欧盟的碳排放交易机制,我们可以清晰的
红色文化是中国人民在中国共产党的领导下,经过了新民主主义革命、社会主义革命、社会主义现代化建设和改革开放的历史阶段,形成了民族的、科学的、大众的、具有中国特色社会
目的探讨对新入院代理制护士进行岗前试岗,护士综合能力效果观察与评估。方法对2011年~2013年招聘入院的护士进行随机分组分别进入观察组和对照组,通过对护士综合能力评估表,
目的分析和探讨颅脑外伤手术治疗后尿路感染的护理体会。方法我们选取2010年6月~2014年6月颅脑外伤术后患者70例,按随机数字表的方法将其分为两组,观察组35例与对照组35例,对
基于数据包络分析法(DEA)测算小额信贷机构的经济绩效与社会绩效,在此基础上探讨其双重目标是否冲突。结果表明,在保持良好盈利的基础上小额信贷机构的经济和社会双重目标可
目的:观察特发性室性心律失常射频消融治疗的成功率、复发率及并发症的发生率。方法选择15例特发性频发室早伴发或不伴阵发性室性心动过速、药物治疗无效或不能耐受长期服药
目的:观察分析在社区居民中建立健康档案对于高血压防治的效果。方法在本社区卫生服务中心辖区内随机选取2011年~2013年的35例建立健康档案管理的高血压患者进行为期1年的观
术后疼痛是机体对手术造成的损伤所产生的一种复杂的生理反应,术后患者均有不同程度的肺功能抑制,术后疼痛使患者不敢用力呼吸导致分钟通气量不足,导致CO2潴留。根据手术部位
目的:探讨分析个性化护理模式应用于乳腺肿瘤术后护理中的临床效果。方法选择2012年12月~2013年12月90例患有乳腺肿瘤患者,由于这些患者之间没有显著性差异,因此,具有可比性,