基于LU-SGS迭代的DGM隐式方法研究 被引量：3

1

作者 郝海兵 张强 +1 位作者 杨永 梁益华 《西北工业大学学报》 EI CAS CSCD 北大核心 2014年第3期346-350,共5页

考虑到LU-SGS迭代法已经在基于非结构网格的有限体积法中得到了成功应用,文章借鉴其思想,将其推广到高精度间断Galerkin有限元隐式格式求解中来,并对其性能进行了研究。为了避免隐式算法中对大型稀疏矩阵求逆,采用LU-SGS迭代法,只需要... 考虑到LU-SGS迭代法已经在基于非结构网格的有限体积法中得到了成功应用,文章借鉴其思想,将其推广到高精度间断Galerkin有限元隐式格式求解中来,并对其性能进行了研究。为了避免隐式算法中对大型稀疏矩阵求逆,采用LU-SGS迭代法,只需要在每步时间推进中沿网格号从前到后和从后到前2次扫描计算即可,并且还能有效降低内存需求。通过对NACA0012翼型和ONERA M6机翼跨声速无粘流动进行数值模拟,计算结果表明:与TVD-RKDG显式时间格式相比,隐式格式所需的迭代步数和CPU时间均得到了很大程度上的减少,并且精度保持不变。 展开更多

关键词 间断Galerkin有限元 lu—sgs 欧拉方程 非结构网格 时间隐式方法

在线阅读 下载PDF 职称材料

隐式时间积分方法的拟动力实验 被引量：6

2

作者 邱法维 《世界地震工程》 CSCD 1995年第3期44-48,共5页

本文介绍了采用隐式时间积分方法实现的拟动力实验。目前拟动力实验中所用的时间积分方法是显式条件稳定的,所以时间步长Δt的选择受到试件刚度和自由度的限制,对于刚度大自由度很多的试件需采用很小的Δt,而Δt太小将造成实验累积误差... 本文介绍了采用隐式时间积分方法实现的拟动力实验。目前拟动力实验中所用的时间积分方法是显式条件稳定的,所以时间步长Δt的选择受到试件刚度和自由度的限制,对于刚度大自由度很多的试件需采用很小的Δt,而Δt太小将造成实验累积误差增大,实验结果失真。隐式时间积分方法是无条件稳定的,Δt的选择不受试件特性的限制,可以比显式算法的稳定极限大很多,从而拓宽了拟动力的应用范围。 展开更多

关键词 隐式时间 积分方法 抗震结构 拟动力实验

在线阅读 下载PDF 职称材料

结构动力响应分析的三阶显隐式时程积分方法 被引量：4

3

作者 陈元昌 张邦基 +1 位作者 张农 郑敏毅 《应用力学学报》 CAS CSCD 北大核心 2016年第2期195-200,366,共6页

基于泰勒级数展开式提出了一种用于结构动力响应分析的高精度时程积分方法,该方法假设t时刻的速度和加速度由t-Δt时刻、t时刻、t+Δt时刻的速度和加速度加权表示,并可根据求解需要调节权值,将积分算法构造成隐式格式或显式格式。通过... 基于泰勒级数展开式提出了一种用于结构动力响应分析的高精度时程积分方法,该方法假设t时刻的速度和加速度由t-Δt时刻、t时刻、t+Δt时刻的速度和加速度加权表示,并可根据求解需要调节权值,将积分算法构造成隐式格式或显式格式。通过理论分析和数值算例,计算讨论了该算法的稳定性和精度,确定了最佳的权值和允许的时间步长。结果表明:本文算法最高具有三阶精度,且具有振幅衰减率低、周期延长率极小等优点。最后结合一个铁道工程实例,表明本文算法适用于大型非线性动态响应的精确快速求解。 展开更多

关键词 动力响应分析 数值积分方法 显隐式 稳定性和精度

在线阅读 下载PDF 职称材料

利用弹簧质点模型和隐式方法的布料模拟研究 被引量：4

4

作者 孔川 罗大庸 《计算机工程与应用》 CSCD 北大核心 2011年第9期225-227,232,共4页

首先运用弹簧-质点模型建立布料的面模型,然后对质点进行力的分解以及受力分析并优化。提出逼近的隐式数值积分方法模拟质点的运动轨迹,这解决了显式数值积分方法的不稳定性和小时间间隔的缺点和其他隐式方法计算量大的缺点,这也是实现... 首先运用弹簧-质点模型建立布料的面模型,然后对质点进行力的分解以及受力分析并优化。提出逼近的隐式数值积分方法模拟质点的运动轨迹,这解决了显式数值积分方法的不稳定性和小时间间隔的缺点和其他隐式方法计算量大的缺点,这也是实现基于物理模型的布料仿真的关键技术。针对具体碰撞对象采用简单的包围盒方法进行碰撞检测,和利用二分法进行碰撞的处理,大大增加了碰撞处理的逼真效果。实验证明其模拟方法具有稳定性和实用性。 展开更多

关键词 弹簧-质点模型 隐式数值积分 包围盒方法 碰撞检测

在线阅读 下载PDF 职称材料

一种数值阻尼耗散可控的结构动力方程积分方法

5

作者 刘伟 童小龙 金蓉 《应用数学和力学》 CSCD 北大核心 2024年第7期922-935,共14页

数值耗散是数值积分方法的重要特性,直接影响数值仿真结果的准确性.对于含有虚假高频振动的动力系统,数值耗散能够改善数值仿真结果,但是对于具有真实高频振动的动力系统,数值耗散则会导致计算结果失真.该研究针对结构动力系统的求解,... 数值耗散是数值积分方法的重要特性,直接影响数值仿真结果的准确性.对于含有虚假高频振动的动力系统,数值耗散能够改善数值仿真结果,但是对于具有真实高频振动的动力系统,数值耗散则会导致计算结果失真.该研究针对结构动力系统的求解,提出了一种数值耗散可控的两子步隐式数值积分方法.基于理论推导,详细介绍了新积分方法的谱半径、稳定性、振幅衰减和周期延长等数值特性.新隐式积分方法通过算法参数α能够对高频虚假振动数值耗散完全可控,相应的耗散比例为1-α,其中-1≤α≤1.通过单自由度动力系统、高频虚假振动系统和多自由度非线性弹簧质量系统三个典型算例,分别证明了新隐式积分方法在计算精确性、高频数值耗散和非线性求解能力方面的优越性. 展开更多

关键词 结构动力方程 隐式积分方法 数值阻尼 放大矩阵 谱半径

在线阅读 下载PDF 职称材料

基于二维非结构网格的GMRES隐式算法 被引量：5

6

作者 李春娜 叶正寅 王刚 《西北工业大学学报》 EI CAS CSCD 北大核心 2007年第5期630-635,共6页

将广义极小残差GMRES(Generalized Minimum RESidual)隐式算法应用到二维非结构网格上,并结合LU-SGS(Lower Upper-Symmetric Gauss-Seidel)方法对所求解方程组的残值向量进行预处理,发展了一套高效、可靠的二维Euler方程的求解器。NACA0... 将广义极小残差GMRES(Generalized Minimum RESidual)隐式算法应用到二维非结构网格上,并结合LU-SGS(Lower Upper-Symmetric Gauss-Seidel)方法对所求解方程组的残值向量进行预处理,发展了一套高效、可靠的二维Euler方程的求解器。NACA0012翼型和某四段翼型的2个算例,表明该隐式算法的计算效率要比传统的四步Runge-Kutta显式算法高出几十倍,与LU-SGS隐式算法的效率相比,该算法的效率高出近1个量级。应用了重启型的GMRES算法,并对2种构造系数Jacobian矩阵的方法进行了比较。 展开更多

关键词 广义极小残差隐式算法 非结构网格 lu—sgs算法 重启型的GMRES算法

在线阅读 下载PDF 职称材料

柔性多体系统动力学典型数值积分方法的比较研究 被引量：3

7

作者 郭晛 章定国 《南京理工大学学报》 EI CAS CSCD 北大核心 2016年第6期726-733,共8页

为合理选用积分方法,对柔性多体系统动力学中8种典型的数值积分方法的性能进行了比较研究。以中心刚体-柔性悬臂梁系统为研究对象,运用第二类拉格朗日方程建立系统高次耦合动力学模型。采用8种典型的数值积分方法对方程进行求解,比较了... 为合理选用积分方法,对柔性多体系统动力学中8种典型的数值积分方法的性能进行了比较研究。以中心刚体-柔性悬臂梁系统为研究对象,运用第二类拉格朗日方程建立系统高次耦合动力学模型。采用8种典型的数值积分方法对方程进行求解,比较了计算效率、数值精度等。结果表明:显式方法较隐式方法更依赖于时间步长的选取;在同等时间步长下,隐式方法的计算效率低于显式方法,隐式方法可以通过放大时间步长提高计算效率;自启动自动变步长且自动变阶的吉尔(Gear)法计算效率高且其更合适于计算广义质量阵为常元素阵的动力学方程;希尔伯特-修斯-泰勒(HHT)法、广义-α法可以通过放大时间步长提高计算效率,但计算精度降低。 展开更多

关键词 柔性多体系统 动力学 数值积分 中心刚体-柔性悬臂梁系统 第二类拉格朗日方程 显式方法 隐式方法 时间步长 广义质量阵

在线阅读 下载PDF 职称材料

一种三维饱和土-基础-结构动力相互作用分析方法 被引量：14

8

作者 陈少林 赵宇昕 《力学学报》 EI CSCD 北大核心 2016年第6期1362-1371,共10页

地震波入射时土与结构动力相互作用分析是地震工程领域的重要问题.由于问题的复杂性,以往的研究大多考虑地基土为干土情形.而实际工程中,土体中经常充满孔隙水,土层往往是含水层或部分含水层.孔隙水对土层的地震反应影响较大,进而影响... 地震波入射时土与结构动力相互作用分析是地震工程领域的重要问题.由于问题的复杂性,以往的研究大多考虑地基土为干土情形.而实际工程中,土体中经常充满孔隙水,土层往往是含水层或部分含水层.孔隙水对土层的地震反应影响较大,进而影响支撑于其上的基础和上部结构的响应.因此,有必要考虑饱和土-基础-结构动力相互作用问题.基于此,土体采用Biot模型,利用集中质量显式有限元方法并结合多次透射人工边界进行模拟;结构经有限元离散后,采用Newmark隐式时步积分方法进行分析,可通过ANSYS等有限元软件实现;基础假定为刚性,采用显式时步积分进行求解;土体和结构(及基础)可分别采用不同的时间步距;通过FORTRAN程序实现了三维饱和土-基础-结构系统在地震作用下的整体分析.从饱和多孔介质动力微分方程出发可知,干土是饱和土的特殊情形,通过将流体体积模量及孔隙率置为零,可将饱和土退化到干土,从而将干土统一到饱和土的计算框架中,通过土-结构相互作用算例对此进行了验证,进一步实现了干土与饱和土混合情形时的土-结构动力相互作用分析,使得问题的模拟更贴近实际工程(如考虑地下水位情形).通过算例对比了饱和土地基、干土地基、干土覆盖饱和土地基(考虑地下水位)情形时,土-结构相互作用对基础和结构响应的影响,结果表明地下水位对基础和结构响应的影响较大. 展开更多

关键词 土-结构动力相互作用 饱和多孔介质 显-隐式积分方法 人工边界条件 地震响应分析

在线阅读 下载PDF 职称材料

一种改进稳定性的时域电场积分方程法 被引量：10

9

作者 唐煜 顾长青 《微波学报》 CSCD 北大核心 2008年第4期31-35,共5页

研究了时域电场积分方程隐式MOT算法中时间基函数对其稳定性的影响,由理论分析和数值结果得出了选用高频分量较少的时间基函数不利于隐式MOT算法稳定性的结论。在此基础上,将一种多项式时间基函数与新的时间平均方法相结合,提出了一种... 研究了时域电场积分方程隐式MOT算法中时间基函数对其稳定性的影响,由理论分析和数值结果得出了选用高频分量较少的时间基函数不利于隐式MOT算法稳定性的结论。在此基础上,将一种多项式时间基函数与新的时间平均方法相结合,提出了一种改进稳定性的隐式MOT算法,计算实例表明比先五步后三步平均法更加精确,比以前的隐式MOT算法更加稳定。 展开更多

关键词 时域电场积分方程 隐式时间步进算法 时间基函数 平均方法

在线阅读 下载PDF 职称材料

高精度高分辨率格式模拟离心叶轮内湍流流场 被引量：2

10

作者 谭大治 袁新 《汽轮机技术》 北大核心 2003年第3期139-141,共3页

应用LU-SGS-GE隐式格式和改良型高阶MUSCLTVD格式,通过求解全三维可压缩Reynolds平均的Navier-Stokes方程和低Reynolds数q-ω双方程湍流模型,对后弯角30°的Krain离心叶轮在设计工况下通道内复杂的三维流动进行了计算。计算结果与... 应用LU-SGS-GE隐式格式和改良型高阶MUSCLTVD格式,通过求解全三维可压缩Reynolds平均的Navier-Stokes方程和低Reynolds数q-ω双方程湍流模型,对后弯角30°的Krain离心叶轮在设计工况下通道内复杂的三维流动进行了计算。计算结果与实验数据符合良好,说明该方法适用于大逆压梯度流动的数值模拟,可以用来深入研究离心压气机内各种复杂的三维流动现象。 展开更多

关键词 数值模拟 q—w山湍流模型 高分辨率格式 MUSCL TVD格式 lu—sgs—GE隐式格式

在线阅读 下载PDF 职称材料

基于边界值方法的微分动力系统数值计算方法 被引量：1

11

作者 汪芳宗 潘明帅 杨萌 《计算力学学报》 CSCD 北大核心 2017年第6期718-724,共7页

对高维非线性初值问题,微分求积法在每一步的积分过程中需要求解一个更高维的非线性方程组,因而计算量巨大。基于微分求积法与边界值方法两者之间的关系,可以将广义向后差分方法和扩展的隐式梯形积分方法看作是经典微分求积法的稀疏表... 对高维非线性初值问题,微分求积法在每一步的积分过程中需要求解一个更高维的非线性方程组,因而计算量巨大。基于微分求积法与边界值方法两者之间的关系,可以将广义向后差分方法和扩展的隐式梯形积分方法看作是经典微分求积法的稀疏表达形式。将广义向后差分方法以及扩展的隐式梯形积分方法这两类边界值方法应用于微分动力系统的数值计算,提出了一类新的数值计算方法。理论分析及算例结果表明,对高维非线性微分初值问题的数值计算,本文方法相对于经典的微分求积法具有更高的计算效率。 展开更多

关键词 动力系统 边界值方法 微分求积法 广义向后差分方法 扩展的隐式梯形积分方法

在线阅读 下载PDF 职称材料

三维弹性动力学的改进的无单元Galerkin方法 被引量：5

12

作者 魏庆节 彭妙娟 《应用力学学报》 CAS CSCD 北大核心 2019年第1期22-27,250,共7页

基于改进的移动最小二乘法建立三维弹性动力学问题的形函数,结合三维弹性动力学的Galerkin积分弱形式,采用罚函数法施加位移边界条件,并引入隐式时间积分,建立了三维弹性动力学的改进的无单元Galerkin方法。该方法由于引入了改进的移动... 基于改进的移动最小二乘法建立三维弹性动力学问题的形函数,结合三维弹性动力学的Galerkin积分弱形式,采用罚函数法施加位移边界条件,并引入隐式时间积分,建立了三维弹性动力学的改进的无单元Galerkin方法。该方法由于引入了改进的移动最小二乘法,避免了病态或奇异方程,在保证计算精度的同时提高了传统的无单元Galerkin方法的计算效率。最后通过数值算例对收敛性进行了分析,并证明了该方法比传统的无单元Galerkin方法计算效率提高了15%。 展开更多

关键词 改进的移动最小二乘法 改进的无单元Galerkin方法 隐式时间积分 弹性动力学

在线阅读 下载PDF 职称材料

动能BGK格式模拟高雷诺数定常流的收敛性研究

13

作者 蒋瑾 钱跃竑 《空气动力学学报》 EI CSCD 北大核心 2012年第4期464-470,共7页

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有... 为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。 展开更多

关键词 动能BGK格式 当地时间步长 隐式lu—sgs方法 多重网格 收敛性 定常流

在线阅读 下载PDF 职称材料

扩展的Euler法求解奇异摄动Volterra型多滞量积分微分系统

14

作者 赵然 王磊 《应用数学》 CSCD 北大核心 2006年第S1期90-92,共3页

本文考虑了用隐式欧拉方法算法求解奇异摄动Volterra型积分微分延迟系统.在本文中给出了隐式Euler方法在一类特殊奇异摄动问题中仿真的算法.文章末尾我们给出了数值试验证明了算法的有效性.

关键词 奇异摄动 延迟 Volterra型延迟积分微分方程 隐式Euler方法 数值实验

在线阅读 下载PDF 职称材料

一种基于经典帕德逼近的传输线瞬态响应数值计算方法 被引量：3

15

作者 付倩倩 陈胜泉 汪赳羚 《电力科学与技术学报》 CAS 北大核心 2019年第4期184-188,共5页

针对均匀有损传输线系统瞬态响应的数值计算,提出一种基于2阶或4阶经典帕德逼近的数值计算方法。在利用有限差分对电报方程的空间坐标进行离散化的基础上,通过2阶或4阶经典帕德逼近方法求解离散后的一阶常系数线性微分方程,从而得到各... 针对均匀有损传输线系统瞬态响应的数值计算,提出一种基于2阶或4阶经典帕德逼近的数值计算方法。在利用有限差分对电报方程的空间坐标进行离散化的基础上,通过2阶或4阶经典帕德逼近方法求解离散后的一阶常系数线性微分方程,从而得到各空间离散点处的时域数值解。数值仿真算例表明:2阶经典帕德逼近方法与隐式梯形积分法以及4阶经典帕德逼近方法与4级4阶显式Runge-Kutta方法的仿真结果完全吻合,从而间接地验证该文方法的高精度和有效性。 展开更多

关键词 传输线系统 经典帕德逼近 电报方程 隐式梯形积分法 RUNGE-KUTTA方法

在线阅读 下载PDF 职称材料

一种求解抛物化Navier-Stokes方程的空间推进算法 被引量：4

16

作者 陈兵 徐旭 蔡国飙 《力学学报》 EI CSCD 北大核心 2008年第2期162-170,共9页

讨论了抛物化NS方程(parabolized Navier-Stokes equations,PNS)的数学性质,对比分析多种处理流向压力梯度的方法的优缺点.以此为基础,成功地将LU-SGS隐式时间积分方法推广到PNS方程的流向空间积分上,发展了基于PNS方程的有限体积单次... 讨论了抛物化NS方程(parabolized Navier-Stokes equations,PNS)的数学性质,对比分析多种处理流向压力梯度的方法的优缺点.以此为基础,成功地将LU-SGS隐式时间积分方法推广到PNS方程的流向空间积分上,发展了基于PNS方程的有限体积单次扫描空间推进算法(single-sweep parabolized Navier-Stokes algorithm,SSPNS).在该算法中,横向无黏数值通量和黏性通量分别采用混合型迎风格式和中心格式求解.用SSPNS算法计算了4个典型流场,包括超声速平板流、15°楔板压缩高超声速流、带攻角的高超声速锥形流和侧压式高超声速进气道流动.SSPNS计算结果与NASA UPS程序数值结果、文献提供的实验数据及理论分析结果符合得很好.对比研究表明,SSPNS法与传统时间迭代法相比,二者计算精度相当,而SSPNS计算速度快1～2个量级,存储量至少低1个量级. 展开更多

关键词 抛物化NS方程 空间推进算法 lu—sgs隐式积分方法 混合型迎风格式 超声速/高超声速流动

在线阅读 下载PDF 职称材料

用空间推进算法模拟高超声速进气道流场 被引量：3

17

作者 陈兵 徐旭 蔡国飙 《宇航学报》 EI CAS CSCD 北大核心 2006年第6期1319-1323,共5页

在有限体积法框架下，采用空间推进算法SSPNS（Single-Sweep Parabolized Navier-Stokes Algorithm）求解抛物化NS方程（即PNS方程），在流向采用LU—SGS隐式积分，而横向无粘和粘性通量则分别采用AUSM系列格式和中心格式计算。用该方法... 在有限体积法框架下，采用空间推进算法SSPNS（Single-Sweep Parabolized Navier-Stokes Algorithm）求解抛物化NS方程（即PNS方程），在流向采用LU—SGS隐式积分，而横向无粘和粘性通量则分别采用AUSM系列格式和中心格式计算。用该方法对1个二维高超声速进气道和2个三维高超声速进气道流场进行了数值模拟，得到的流场波系结构、壁面压力及传热系数分布与文献中相关数值解和实验数据基本一致，表明SSPNS法能够准确地模拟超燃冲发动机进气道内的高超声速流动。对比研究表明，SSPNS法与求解FNS（Full Navier＼lStokes Equations）方程的传统时间迭代法相比，二者计算精度相当，而SSPNS计算速度快1～2个量级，存储量至少低1个量级。本文的研究为CFD在超燃冲压发动机部件及一体化优化设计中的集成，以及大型高超声速工程流动的高效计算，打下了良好的基础。 展开更多

关键词 抛物化NS方程 空间推进算法 高超声速进气道 lu—sgs隐式积分 AUSM系列格式

在线阅读 下载PDF 职称材料

基于物理模型的水波动画 被引量：4

18

作者 秦培煜 陈传波 +2 位作者 解洪涛 夏晖 吕泽华 《计算机科学》 CSCD 北大核心 2007年第10期238-240,共3页

近年来,自然景物的模拟一直是计算机图形学中最具挑战性的问题之一。关于山、水、云、烟、火焰等自然景物的模拟,在计算机游戏、影视、广告等各种领域中有着广泛的用途。作为自然景物模拟的重要内容,对流水、波浪的模拟正日益引起人们... 近年来,自然景物的模拟一直是计算机图形学中最具挑战性的问题之一。关于山、水、云、烟、火焰等自然景物的模拟,在计算机游戏、影视、广告等各种领域中有着广泛的用途。作为自然景物模拟的重要内容,对流水、波浪的模拟正日益引起人们的关注。本文基于二维浅水波方程模型,给出了一种实现水波动画的数值模拟方法。和以往的方法不同的是,本文既没有去特意构造具体的波形函数,也没有去求解复杂的Navier-Stokes方程,而是基于二维浅水波方程模型,采用隐式半拉格朗日积分方法进行求解,在保证稳定性的同时,可以允许大时间步长。实践证明,用该算法模拟水波,效果比较真实,而且在普通的PC平台上即可满足一般动画的实时需求。 展开更多

关键词 水波动画 浅水波方程 半拉格朗日方法 隐式积分

在线阅读 下载PDF 职称材料

拱坝横缝影响及有效抗震措施的研究 被引量：13

19

作者 林皋 胡志强 《世界地震工程》 CSCD 2004年第3期1-8,共8页

大量研究结果和某些拱坝的地震震害表明,横缝对拱坝的地震响应有很大的影响。通过采用非光滑方程组方法以及考虑碰撞时刻动量、动能守恒来模拟横缝所引起的动接触问题,同时为了提高计算效率,采用隐-显式积分方法对坝-基系统的动力平衡... 大量研究结果和某些拱坝的地震震害表明,横缝对拱坝的地震响应有很大的影响。通过采用非光滑方程组方法以及考虑碰撞时刻动量、动能守恒来模拟横缝所引起的动接触问题,同时为了提高计算效率,采用隐-显式积分方法对坝-基系统的动力平衡方程进行求解。针对在拱坝中上部配筋这一抗震措施,也作了探讨。通过对小湾拱坝的分析,为高拱坝工程抗震措施的选择提供技术依据。 展开更多

关键词 拱坝抗震 横缝 动接触 非光滑方程组方法 隐-显式积分

在线阅读 下载PDF 职称材料

超声速进气道可压及不可压流动数值模拟 被引量：2

20

作者 陈兵 徐旭 蔡国飙 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2006年第1期31-35,45,共6页

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmet-ric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstrea... 采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmet-ric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致. 展开更多

关键词 超声速进气道 预处理方法 lu-sgs隐式时间积分 AUSM%PluS%(P)格式

在线阅读 下载PDF 职称材料