小行星轨道受木星的平均摄动研究

引进行星系统的基本作用量子η＝λＧＭ＊／ｃ≈９．２８×１０１４ｍ２ｓ－１和行星椭圆轨道的无量纲根数,推导以这种根数表示的小行星轨道受木星的平均摄动函数和三体摄动方程;再由摄动方程计算小行星轨道受木星的平均摄动;最后,给出了几颗特殊小行星轨道受木星的一阶平均摄动数值     

对流扩散方程的数值流形格式及其稳定性分析

将流形方法应用于对流扩散方程的数值求解,建立了基于标准Galerkin加权余量法的定常无源对流扩散方程的数值流形格式,采用一维定常无源对流扩散方程证明了物理覆盖的覆盖函数取完全一阶多项式的标准流形格式具有绝对的数值稳定性,并通过与一维对流扩散方程有限元解、精确解的对比,对该数值流形格式的稳定性进行了验证.同时,将基于四节点矩形有限单元覆盖系统的数值流形格式应用于二维平行管道中定常热对流扩散问题的数值分析.结果表明:在小的单元Pe(Pe<2)时,流形解的精度较有限元方法显著提高;在较大单元Pe条件下,一阶多项式覆盖函数的标准流形格式虽然绝对稳定,但假扩散作用显著,得到的数值解与真实结果存在较大的偏差.     

人造卫星相对论摄动运动方程的解

本文应用平均根数法求得了人造卫星在地球形状摄动下,相对论摄动运动方程的二阶长期摄动和一阶长周期摄动的解.     

相对论摄动运动方程的解

本文在考虑相对论效应的前提下,研究了地球形状对人造卫星运动的影响。用平均根数法求得了卫星的相对论摄动运动方程在地球形状摄动下的一阶长期摄动和一阶短周期摄动解,并将结果与非相对论摄动解进行了比较。     

定常对流扩散反应方程的指数型高阶差分格式

提出了一种数值求解一维定常对流扩散反应方程的指数型四阶差分格式．首先,由常数变易法求得模型方程的通解并应用到xi-1,xi,xi＋13点上,得到模型方程的指数型差分格式．然后,利用源项f（x）在点xi处的二阶泰勒展开,得到定常对流扩散反应方程的指数型四阶差分格式．最后,用数值算例验证了该格式的高阶精度和可靠性．     

一类变时间分数阶含源项非定常奇异摄动对流扩散方程的数值分析

考虑了变时间分数阶含源项非定常奇异摄动对流扩散方程的数值逼近问题.首先采用分段线性插值法,结合对一阶时间导数的一个二阶近似离散Coimbra变时间分数阶导数,然后用中心差分离散一阶空间分数阶导数和二阶空间分数阶导数,最后用数值例子验证了提出的数值方法,说明了数值方法的有效性.     

一类微分方程解与其平均方程解之间的关系

研究了一类微分方程解与其平均方程解之间的关系,当相应的平均方程的解为定常解时,我们所得到的原方程的一阶近似解必为周期解;而当相应平均方程的解为周期解时,原方程的一阶近似解必为调幅振动解,且在满足一定的条件下,还可以是周期解.     

有阻尼Duffing方程的二阶次谐振动及其稳定性研究

本文借助渐近法给出了有阻尼Duffing方程二阶次谐和解的解析表达式,并讨论了定常运动的稳定性。给定了一系列参数,应用Rungc-Kutta法作了数值计算,绘制了x(f)、■(t)的时间历程图以及x-■相平面图。此外,还提供了无阻尼Duffing方程的八阶和十六阶次谐振动的计算结果。所有这些结果都是十分满意的。     

Lagrange三次有限体积元法的超收敛现象

基于三角形网上求解Poisson方程的Lagrange三次有限体积元法, 给出了超收敛性的数值结果. 数值实验表明, 在三角形单元的对称点(即3边中点和3个角顶点)上, 数值解平均梯度的收敛阶约为4阶, 比按H¹模的收敛阶(O(h³))约高一阶.     

非线性椭圆奇异摄动问题的数值方法

本在非等距离网格上给出了非线性随圆异摄动问题的解，这种方法是基于积分方程的数值解法，证明了差分格式的一阶一致收敛性，并给出了数值例子。     

机翼跨音速颤振特性的频域分析

从气动力影响系数矩阵的物理意义出发,推导出以气动力时域计算为基础的跨音速气动力影响系数矩阵的计算公式.通过求解非定常Euler方程,获得机翼简谐振动下的跨音速气动力时间历程,应用最小二乘法对模态气动力系数的时间历程进行曲线拟合,得到跨音速气动力影响系数矩阵,进而在频域内进行机翼的跨音速颤振特性分析.采用该方法计算了AGRAD 445.6机翼随马赫数变化的颤振边界,所得结果与文献的实验结果吻合良好,验证了该方法的可行性和准确性.     

三维超音速和高超音速无粘非定常流动的数值模拟

对于钝头体的三维高超音速和超音速非定常绕流，给出了Ｅｕｌｅｒ方程在团结于钝头体 的非惯性坐标系中的守恒形式，通过对求解域的有限体积离散，矢通量分裂方法以及Ｗａｒｍｉｎｇ 和Ｂｅａｍ的二阶迎风格式用数值方法计算了Ｅｕｌｅｒ方程的时间精确解。给出了利用本文的非定 常流动数值计算结果及１９９２年任玉新的学位论文“气动稳定性导数的理论与数值方法”计算出 的俯仰阻尼导数的值，它与实验结果较好地符合。表明文中的非定常流动数值计算方法是可靠 的，对准确确定飞行器的动导数具有意义。     

TVD有限体积法在喷流计算中的应用

该文根据ＷａｎｇＪ．Ｃ．Ｔ．的二阶ＴＶＤ有限体积法求解守恒型Ｅｕｌｅｒ方程，对原文中的数值格式进行了部分修改。计算网格单元是非均匀的任意四边形。文中计算了楔形风洞中的斜激波反射问题和无粘喷流近场流动结构问题。数值计算结果表明，该文方法具有激波分辨率高、计算稳定、伪振荡小等优点。同时该文方法可以处理任意复杂边界的二维及轴对称流动问题的计算，也可以将该文的有限体积法思想推广到三维具有复杂边界的流动的计     

SF6压气式灭弧室冷态气流场与介质恢复强度计算

应用Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ法对可压缩气体非稳态控制方程组进行了数值求解，计算了ＳＦ６压气式断路器灭弧室在不同压比，不同喷口形状下的流场分布。通过对计算结果的分析，讨论了在不同情况下灭弧室对开断小电容性电流后的过电压的承受能力。     

二维机翼非定常气动反命题的欧拉解法

提出在映象坐标下求解二维非定常欧拉方程获取二维振荡机翼的气动反命题解,给出了非定常气动反问题的边界条件和求解方法．该方法克服了以往基于势函数方程解反命题的局限,统一了亚音、跨音与超音流反问题的解法,并避免了非定常尾涡面的直接求解与反问题未知边界的不确定性．     

二维水库水流泥沙数值模拟

根据水库的水沙运动特点及库区地形条件,考虑悬移质和推移质运动,建立了二维水流泥沙数学模型。网格单元用拉普拉斯方程和代数方法生成,模型方程用有限元法离散求解。利用水槽试验资料和水库实测资料进行模型验证和稳定性分析,结果表明：计算的断面流速分布、流场和河床变形与实测结果基本吻合。     

基于动态嵌套网格方法的摆动翼型粘性绕流数值模拟

该文根据摆动翼型的流场特性运用动态嵌套网格技术将其划分为不同的子区,在不同的子区中使用不同的流动控制方程,以提高计算效率.对翼型生成一个随其摆动的相对较小较密的贴体网格,在它上面采用Navier Stokes方程;对整个流场生成了一个大的相对较稀的静止背景网格,在它上面采用Euler方程.两个网格间的流场信息交换通过动态嵌套网格技术来实现.通过与单块翼型网格的计算过程和计算结果的对比分析,表明该方法在不损失计算精度的前提下计算非定常问题时具有相对较高的计算效率.     

翼型绕流不同频率展向振荡电磁力减阻控制实验研究

研究结果表明展向振荡电磁力可控制湍流边界层,电磁力的振荡频率对湍流的控制效果有影响,但并未讨论电磁力振荡频率对控制效果的影响机理。实验研究了不同频率展向振荡电磁力控制翼型绕流的减阻效果及其影响机理。实验在转动的水槽中进行,在翼型的背风面包覆展向振荡电磁力激活板,并将其浸入水槽中,利用应变传感器测量翼型的阻力,基于意法半导体公司生产的微处理器开发电磁力控制器,用于控制电磁力的方向和振荡频率。研究结果表明展向振荡电磁力对翼型绕流具有减阻效果,对比分析了不同频率的展向振荡电磁力的减阻效果,发现电磁力的振荡频率为20 Hz时减阻效果较优,减阻效率可达到18%;展向振荡电磁力可减小翼型阻力的振动幅值,具有减震功能;当电磁力的振动频率与阻力曲线内小波动频率相近时,电磁力的减阻减震的效果最佳。     

振荡管内非定常流动的数值模拟

建立了振荡管内非定常流动与传热的数学模型 ,分别采用Godunov格式和控制容积法对气动方程及传热方程进行离散 ,得到了管内压力波形、气流速度及温度随时间的变化历程 .所得管内压力波在波形和相位上与实验结果吻合较好 ;排气初期的气流速度最大且温度最低 ;管内出现强反射激波 .提出的求解Riemann问题迭代初值可使迭代快速、稳定收敛     

大曲率结构化表面湍流流场特性对比研究

湍流模型选择对于小尺寸大曲率圆管内软性磨粒流湍流流场模拟计算具有重要影响.针对此问题,提出使用标准k-ε模型与考虑两相磨粒流平均流动中的旋转及旋流流动情况的重整化群(RNG)k-ε模型进行仿真对比研究的方法.基于颗粒动力学理论的欧拉模型,对软性磨粒流的控制方程进行求解.通过使用两种湍流模型,对小尺寸大曲率圆管内湍流流动进行仿真计算,得到了两湍流模型下软性磨粒流流场的速度、压力、湍动能等数值模拟结果.仿真结果对比分析表明:标准k-ε模型模拟得到的速度值与实验值存在一定误差,而RNG k-ε模型仿真的速度值误差相对较小;标准k-ε模型对湍动能的估值计算比RNG k-ε湍流模型高,难以计算有滞止点的流动过程;RNG k-ε湍流模型能更有效地模拟有大曲率带分离的湍流流动,从而证明该模型更适合小尺寸大曲率圆管内软性磨粒流湍流流动的研究.