预估-校正LU-SGS的隐式算法

为了进一步提高Euler方程求解效率,在传统LU-SGS格式基础上提出了预估-校正LU-SGS。结合高分辨率迎风通量格式,发展了一套高效的Euler方程求解器。以NACA0012翼型和RAE2822翼型跨音速无粘流动作为算例,对比分析了预估-校正LU-SGS、传统LU-SGS和GMRES格式的计算效率。结果表明:预估-校正LU-SGS算法能显著提高求解效率,其效率接近GMRES算法,比传统LU-SGS方法提高了近三倍。     

有限元网格结点编号

在有限元分析中,求解高阶线性代数方程组时整体刚度矩阵所需存储与由网格结点编号决定的顺序有关.在基于等带宽存储的求解法与基于变带宽存储的求解法的基础上推导出它们的关系.据此,提出了有限元网格结点编号的前沿法与矩形法,并给出了这两种编号法的内存消耗与结点数量的关系.理论分析和实例表明这两种编号法能有效地减少计算机内存消耗.     

非线性弹性材料的三阶本构方程

推导了可压和不可压非线性弹性体的本构方程.应力张量,作为以单个应变张量为变量的张量值函数,用含有高阶弹性张量的张量多项式来表示.利用各种对称性来简化这些表达式,最后得到了各向同性情况下的本构方程和应变能函数.得到的表达式是完备和不可约的,满足张量函数的表示定理.     

快速大时间步长熵条件格式的分辨率研究

提出了高分辨率快速大时间步长熵条件格式的构造方法.用激波管问题对一族熵条件格式进行研究.在精度、步长、限制器方面进行了详细的数值实验,研究了同样计算量下各种格式的表现品质.从理论上保证了大时间步长格式的无振荡性质,从具体的数值实验分析中确定了大时间步长格式的分辨率问题.      

特定消谐逆变器开关角度的复形调优解法

针对特定消谐式逆变器开关角度难于求解的问题,提出了一种基于约束条件下n维极值复形调优法的开关角度求解方法,分析了目标函数和约束条件的设计方法.该算法对开关角度初值的选取只需考虑约束条件,大大增加了开关角度初值的取值范围;求解过程中不存在解发散或奇异无解问题,与迭代法解非线性超越方程组相比可靠性更高;通过约束条件的限制,使求解过程自动消除了可能出现的窄脉冲;通过目标函数和加权因子可灵活的进行开关角度优化.算例和实验结果验证了方法的正确性和有效性.      

吴方法在平面并联机构位置正解中的应用

采用吴方法对平面并联机构位置正解问题进行了研究.吴方法是一种求解非线性方程组的数学机械化方法,采用这种方法任何非线性方程组都可以在有限步内得到解决.在给出了吴方法基本原理的基础上,对本问题进行了求解,并将原始方程组转化成为一个三角化的方程组.其中单变量方程的次数为6次,说明平面并联机构可以有6个不同的位姿.最后用数值实例进行了验证,给出了计算结果.吴方法在这一问题中的应用,为求解其它机构学难题提供了新途径.     

垫块法半模型风洞试验的数值模拟研究

以M6机翼半模型为算例,采用数值方法模拟垫块法风洞试验洞壁边界层及垫块厚度对半模型试验结果的影响.控制方程采用全场三维可压、雷诺平均的Navier-Stokes方程;计算格式采用Jameson中心有限体积、多步Runge-Kutta时间步长格式进行求解;湍流模型采用Baldwin-Lomax模型.通过本文的研究,可以为跨声速风洞半模型试验垫块厚度的选取提供计算方法和设计依据.     

应用于跨声速旋翼气动声学计算的两种时域方法的比较

高速脉冲噪声(HSI)与流动非线性效应密切相关,是跨声速旋翼的主要声源之一。采用理论分析和数值实验相结合的方法,对预测HSI噪声的FfowcsWilliams-Hawkings方法(简称FW-H方法)和Kirchhoff方法进行了比较和分析。理论分析从两种方法的控制方程出发,采用广义函数方法对方程源项进行整理和比较,证实在线性区域Kirchhoff方程可视为FW-H方程的有效近似。数值实验以UH-1H悬停旋翼跨声速气动噪声计算为例,由三维Euler方程数值模拟提供近场气动数据,并在同一控制面上分别采用FW-H方法和Kirchhoff方法对远场噪声进行定量预测。计算结果表明,如果适当构造控制面,2种方法均能有效地预测跨声速旋翼非线性HSI噪声。进一步研究了控制面位置及控制面上物理参数时间导数计算精度对计算结果的影响,发现控制面位置对Kirchhoff方法计算结果具有决定性的影响,而控制面上物理参数时间导数的计算精度对FW-H方法计算结果具有重要影响。     

多守恒性的差分格式的分析和比较

球面大气浅水波方程(在忽略摩擦和外源强迫的条件下)具有总能量、总质量、总位涡、总位涡拟能和总角动量守恒等五个重要的物理守恒性.然而,将大气浅水波方程差分离散化后,常常不能保持这些守恒性.本文对新发展的两种较好的差分格式进行比较;一种可称之为修正总能量守恒格式,它可以较好地保持上述五个守恒性中的四个;而另一种可称之为准辛格式,它可以使上述五种守恒性均近似保持.文中对这两类格式作了具体分析,并作相应的数值试验和比较,结果表明这两类格式均值得推广应用.     

平衡气体效应对飞行器动态特性的影响研究

通过数值求解N-S方程,模拟升力体外形的强迫振动过程,研究了平衡气体效应对带翼飞行器俯仰阻尼导数的影响。分别在完全气体和平衡气体假定条件下,给出了h=30km,Ma∞=8,α=0°,10°,20°,30°时的俯仰阻尼导数变化的定量分析结果。