Chebyshev-Gauss-Lobatto节点的一个注记

利用Lagrange插值基函数和Chebyshev多项式的性质,推导以Chebyshev-Gauss-Lobatto点为插值点构造的插值基函数的一阶、二阶微分矩阵的显示格式,并由插值点的性质得出两者之间的关系.通过对具有解析解的一维对流扩散方程进行数值求解,验证了一阶、二阶微分矩阵显式格式的正确性.数值结果表明:由微分矩阵显式格式可以方便地构造配置点谱方法中的拟谱算子,利用其求解微分方程,在较少的网格点时,即可得到快速收敛的高精度的数值结果.研究工作对配置点谱方法的应用具有一定的理论指导意义.     

计算矩形直通道内紊流流动的紊流代数应力模型的研究

应用全三维紊流计算方法，采用由Ｌａｕｎｄｅｒ、Ｒｅｅｃｅ及Ｒｏｄｉ提出的紊流代数应力模型（ＬＲＲＡＳＭ），对矩形直通道内的复杂紊流流动进行了计算，提出了用无量纲距离ｙ＊改进ＬＲＲＡＳＭ中考虑固壁效应的计算方法，并进行了多种不同形式的对比计算．用改进的考虑固壁效应计算方法所计算的主流速度、紊流量均与实验值符合良好     

谱元方法求解波动方程及影响其数值精度的相关因素

为探讨波动方程的高精度数值模拟,采用Chebyshev谱元方法结合隐式Newmark时间积分方法求解波动方程.求解一个具体算例验证了数值方法的可行性,讨论了时间步长、Newmark因子以及计算区域的网格剖分方式对数值精度的影响.结果表明:和差分法相比,谱元方法求解波动方程具有所用网格节点少,数值精度高的特点;数值误差随时间步长减小而减小;在满足稳定性要求的前提下,数值误差随着Newmark因子的减小而减小;当总网格节点数相同时,不同的网格剖分方式所得数值误差不同.所述方法和结论可用于模拟声波在空气中的传播.     

Chebyshev谱元法求解含吸收边界的二维均匀稳定流场的声传播

从群速度的角度推导了包含均匀稳定来流的二维波动方程的1阶吸收边界条件,基于Che-byshev谱元法提出了二维均匀稳定来流波动方程的求解方法.在空间上采用谱元方法,在时间上采用隐式Newmark积分法,从而获得了波动方程的离散形式.经具体算例验证表明:与1阶Clay-ton-Engquist-Majda吸收边界条件相比,所推导的吸收边界条件能更有效地削弱边界上的数值反射,避免解的失真,求解方法在空间上具有谱精度,在时间上达到了2阶精度.     

低速气流中二元叶片颤振数值模拟与Hopf分岔分析

针对低速来流,建立了一种低速气流中的二元叶片颤振数值模拟、运动稳定性分析的方法.以速度势函数形式的Laplace方程、非定常Bernoulli方程作为流动控制方程,结合二元叶片沉浮、扭转结构动力模型,通过迭代的方式实现了叶片颤振过程的数值模拟.同时,从非线性动力学角度对计算结果进行了分析,对颤振产生的机理进行了研究.对不同来流速度下某二元叶片的自激振动进行了模拟,结果表明:来流速度是一个影响叶片运动稳定性的重要因素,来流速度增加而导致的叶片颤振是以来流速度为分岔参数的Hopf分岔产生的结果.     

谱元方法求解波动方程时显式与隐式差分方法的比较

分别将显式中心差分和隐式Newmark差分方法与Chebyshev谱元方法相结合,探讨了当采用谱元方法求解气动噪声问题时,这2种时间差分方法对数值精度和计算稳定性的影响.通过模拟噪声扰动在自由空间的传播过程,比较了2种时间差分方法的数值误差,研究了显式中心差分方法的稳定性条件.由此得出以下结论:当显式中心差分方法稳定时,2种差分方法均可得到有效的数值解,在同一时刻的同一网格节点上,Newmark方法的数值误差约为显式方法的2倍左右;当计算采用的时间步长大于显式中心差分方法的临界步长时,显式方法发散,Newmark方法的数值结果仍和解析解保持一致.     

谱元方法求解含吸收边界的声学问题

推导并实现了带吸收边界条件的二维波动方程的切比雪夫谱元解法．该解法引入一阶Clayton-Engquist—Majda吸收边界条件，在空间上利用谱元方法，在时间上利用中心差分的积分方法得到波动方程的离散形式，并给出具体算例进行了验证．结果表明：该解法在空间上具有谱精度，在时间上达到二阶精度；与第一类边界条件相比，吸收边界有效地削弱了边界上的数值反射，避免了解的失真；使用中心差分的时间积分方法同隐式积分方法相比，适合波的传播特性，避免了矩阵求逆运算，并且占用的计算机内存小．     

谱元方法求解对流扩散方程及其稳定性分析

探讨了一维对流扩散方程的一种高精度数值解法,该解法在空间上采用了Chebyshev谱元方法,在时间上结合了半隐式Adams方法。通过数值算例验证了解法的可行性,利用特征分析法得到了对流扩散方程谱元求解时不同离散形式的稳定性条件,并对数值求解的稳定性进行了预测。通过时间步长、网格划分、对流项和黏性项插值阶数的影响研究表明:耦合Chebyshev谱元方法和半隐式Adams方法在求解对流扩散方程时能够获得高精度的数值解;时间离散时Adams方法的黏性项采用一阶插值形式、对流项采用二阶插值形式,在未增加计算量的同时能够获得较大的稳定区域和较高的计算精度。     

用切比雪夫谱元方法求解均匀流场中的声传播问题

利用谱元方法中的无穷光滑插值函数的高阶精度特点,结合隐式时间推进算法的稳定性,推导并实现了低马赫数均匀流场中声波动方程的切比雪夫谱元解法,进而得到了流场影响下的声传播问题的数值解.该解法对均匀流场中的声传播问题在空间上进行谱元离散,在边界上引入Clay-ton-Engquist-Majda吸收边界条件,在时间上利用隐式Newmark积分方法推进求解.算例与解析解的对比验证表明:该解法在空间上可以实现高阶精度,在时间上达到2阶精度;使用的隐式New-mark时间积分方法稳定性好,计算工作量相对较小;当数值解达到稳态传播时和解析解吻合得非常好.随着计算条件的飞速发展,加密网格并采用更高阶的切比雪夫谱逼近可以进一步提高精度,以适应计算气动声学的精度要求,另外可尝试采用更高精度的吸收边界条件以改善边界反射对计算声场的干扰.