高精度紧致差分方法及其应用研究

来源 :上海大学 | 被引量 : 0次 | 上传用户:bluebird2000
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
正确数值模拟多尺度复杂流场是计算流体力学(CFD)中最为重要,也是最具有挑战性的课题之一.近年研究表明,在提高数值模拟可靠性和有效性方面,高精度和高分辨率要求已成为技术发展中的一个决定性因素.紧致差分格式因其网格基架点少、精度高、边界易处理且能达到与谱方法相近的分辨率等优点,日益受到人们的重视.本文基于在应用中提出的问题和未来发展的需求,针对粘性不可压缩复杂流动的数值模拟和高阶紧致(high-order compact,HOC)差分方法,研究了对流扩散问题的指数型高阶紧致(ItOC exponential type,EHOC)差分方法和“导数型”紧致差分逼近的高阶组合格式,并建立了求解原始变量和涡量流函数形式的粘性不可压缩Navier-Stokes(N-S)方程、Navier-Stokes/Boussinesq(N-S/Boussinesq)方程等的高精度差分格式与算法,直接数值模拟了驱动方腔流、自然对流换热、双扩散对流系统等复杂流动问题,得到了若干具有理论和实际意义的成果. 全文共分九章,要点和主要工作如下: 第一章为绪论,概述了CFD数值方法,其中以粘性不可压缩流动问题的计算方法和计算技术为重点;特别介绍了高精度紧致型差分逼近方法和格式的研究进展;同时还概述了本文研究的主要工作内容. 第二章的主要目标是建立定常对流扩散模型问题的具有较高分辨率的无振荡EHOC方法.从建立一维常系数对流扩散方程的指数型紧致格式入手,利用余项修正技术、降维法和符号紧致算子技术等,分别提出了—维变系数、二维常系数/变系数对流扩散模型问题的EHOC格式;利用Fourier分析法,对一维常系数问题的离散逼近式进行定量分析,讨论和分析了离散式的行为特性,并与其它高阶紧致格式相比较,得出了EHOC格式的色散误差和耗散误差最小;数值结果表明,本文的EHOC方法不仅具有高分辨率和高精度,而且是稳健的. 在第三章中,运用符号紧致算子技巧,提出了非定常对流扩散方程具有系数矩阵严格对角占优的指数型高阶紧致ADI(EHOC ADI)格式.以建立的定常一维常系数/变系数对流扩散方程的指数型四阶紧致格式为基础,分别提出了非定常二维常系数/变系数对流扩散方程的EHOC ADI方法;以常系数模型方程为例,从理论上分析了所建立格式的精度和稳定性,给出了系数矩阵严格对角占优的证明;数值实验结果表明,本文提出的:EHOC ADI格式不仅精度和求解效率高,而且适用于对大梯度或边界层问题的求解,优于其它求解非定常对流扩散问题的HOC方法. 第四章的主要目的是引入组合紧致迎风(CCU)差分方法,为建立流动与传热问题的高精度差分算法打下理论基础.利用解析特解方法构造紧致型差分逼近式,推导出可任意组合的三点或五点紧致型差分逼近式的基本格式;在此基础上,提出了三阶、四阶三点组合紧致迎风(CCU34)格式和四阶、五阶五点组合紧致迎风(CCu45)格式;采用Fourier分析法对组合紧致迎风格式的数值解的精度及行为进行了分析,并通过数值算例对组合紧致迎风格式的有效性作了进一步的验证. 在第五章中,运用在第二章和第三章为对流扩散模型问题建立的EHOC格式,提出了涡量流函数形式的不可压缩N-S方程和N-S/Boussinesq方程的指数型高阶紧致型差分逼近和算法,给出了细致的算法描述;数值求解了驱动方腔流和方腔自然对流问题,验证了本文方法的可靠性和有效性;在此基础上对带后台阶封闭腔内自然对流换热问题进行了直接数值模拟,重点研究了不同Rayleigh数下腔体台阶下方宽度和腔体台阶上方高度变化对腔体内能量传递的影响. 第六章提出了求解粘性不可压缩N-S方程的高精度差分投影算法.基于投影法和第四章提出的组合迎风紧致格式,在交错网格与非交错网格系统上分别建立了数值求解原始变量形式的二维非定常、粘性不可压缩N-S方程的高精度紧致投影算法,并相应地提出了求解压力梯度的显式四阶紧致格式和压力Poisson方程的四阶紧致格式;对两个带有周期边界条件的流动问题分别在交错网格和非交错网格上进行了数值试验,验证了方法的精度,考核了算法的有效性;在此基础上,对中高Reynolds数下的平面驱动方腔流进行数值模拟,确定腔内流动第一次分叉的临界Reynolds数的范围,描述和分析了Reynolds数等于10000时方腔流动周期性演化的过程,并确定了周期解的周期. 第七章建立了求解原始变量形式的粘性不可压N-S/Boussinesq方程的四阶五阶组合紧致迎风(CCU45)差分投影算法,并直接数值模拟了倾斜闭腔内自然对流换热问题.采用投影法和四阶五阶组合紧致迎风(CCU45)格式,构造了求解二维非定常不可压缩N-S/B方程的高阶紧致差分投影算法;对方腔自然对流换热问题进行了数值模拟,将所得计算结果与精细网格上的“标准解”比较,验证了本文建立的高精度紧致差分投影算法的有效性和可靠性;对重力场作用下的与倾斜角度有关的闭腔内自然对流换热问题进行了直接数值模拟研究,探讨了腔体倾斜角度、Rayleigh数、初值条件、边界条件和腔体尺寸对腔内流体流动状态的转变及热量传递的变化规律的影响;研究了高Rayleigh数情况下的闭腔内的流场、温度变化规律,描述了不同倾斜角度下腔体内流动出现周期性变化的过程. 第八章针对双扩散对流系统的数值模拟,基于第四章提出的三阶四阶组合紧致迎风(CCU34)格式,构造了温度方程、浓度方程和涡量方程的高精度紧致格式,建立了涡量流函数形式的N-S/Boussinesq方程的高精度紧致数值方法和算法,给出了细致的算法描述;直接数值模拟了高宽比为2:1闭腔内的双扩散对流系统,将所得结果与文献中的结果进行对比,验证了所提出算法的正确性;探讨浮力比和腔体尺寸对腔内流场、温度场和浓度场的影响,直接数值模拟研究了腔体高宽比为4:1的双扩散对流系统,确定流场发生周期振荡的温度浓度浮力比的临界值范围,描绘了部分浮力比下流场发生周期振荡的全过程,研究分析了流场产生周期振荡的机理及性质. 第九章为本文研究工作的主要贡献、结论和展望.
其他文献
永磁行星齿轮是一种新型的磁力耦合传动装置,它由永磁太阳轮、永磁行星轮、永磁外齿圈以及行星架组成。借助于三个端口,即太阳轮端口、外齿圈端口和行星架端口,永磁行星齿轮可实
近年来电力市场不断发展和完善,电力体制改革不断深化,厂网分开、试行竞价上网等市场化运作方式对发电企业提出了新的要求。电厂除了重视关口表电量外,内部的发电量、厂用电损耗的考核也越来越受到重视。目前电厂现有的生产管理信息系统MIS(Management Information System)虽然具有强大的功能,但是在某些具体的管理应用上不能满足更深层次的要求,例如市场化运行后,不但要计量有功电能,还要
异步电动机负载是电力系统中最重要的负载,它所带的机械负载有不同的特性,有的转矩与转速无关,有的则随转速而变化。有变化比较缓慢的一般性负载,也有在短时间内发生急剧变化的冲击性负载,因此,特性十分复杂,要进行准确的模拟,是一个十分复杂和困难的课题。传统的模拟方法是采用物理模拟的方法,在进行模拟时,首先应考虑异步电动机本身参数和特性的模拟,其次是阻力机械特性的模拟,总的来说,物理模拟是采用相似原理制造出
在特种电源领域,逆变电源已逐渐成为其发展重点。而随着当今社会对节能、可靠性要求的提高,软开关技术已在逆变电源中占据了不可取代的地位。本文创造性的提出了一种焊割电源的理念,将移相全桥软开关技术应用其中,并对其进行了系统的研究。文章对零电压开关(ZVS)和零电压零电流开关(ZVZCS)两类软开关进行了比较分析,并对ZVZCS软开关的工作原理进行了详细分析,总结了几种典型的实现方案,最终采用由饱和电感实
在电力系统中,无功补偿与无功平衡对于提高全网经济效益和改善供电质量至关重要。因此,新型无功补偿装置的研制与应用显然已成为电力系统中的重要研究课题。近年来新型无功补偿装置发展迅速,种类也日益增多。电网中现有的无功补偿装置主要包括同步调相机、晶闸管投切电容器(TSC)、晶闸管控制电抗器(TCR)、静止无功发生器(STATCOM)等。同步调相机由于其响应速度慢、运行维护困难等诸多缺点,目前已很少使用。晶
协同神经网络是Haken提出的一种模式识别方法,它是一个自上而下的模式识别过程。本文在分析经典.Haken协同神经网络的主要关键技术问题的基础上,提出了一种改进的协同神经网络新
电动执行器是常见的执行元件,具有定位点多、定位精度高等特点,在诸多领域广泛应用。为了节能、减小运行成本,用户在选择执行器会考虑实际运行工况下执行器的能耗。电动执行