基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

单斜面膨胀喷管几何参数对流场和性能的影响

基于CFD技术,通过求解R eynolds平均N-S方程附以标准k-ε双方程湍流模型,数值研究了单斜面膨胀喷管(SERN)主要几何参数对其内部流场和性能的影响,计算结果表明,SERN扩张段上、下斜面长度和角度对喷管流场结构和性能有明显的影响,随着上斜面长度的增加,喷管从欠膨胀状态过渡到过膨胀状态,在完全膨胀状态下,推力系数最大,存在某一最佳的上斜面角度,在该角度下,喷管推力系数最大,同样地,对于下斜面也存在最佳斜面长度和角度。      

基于全解耦奇偶方程的动态系统执行器故障检测与识别

给出了基于全解耦奇偶方程的动态系统执行器故障的检测与识别方法。讨论了全解耦奇偶向量的产生方法,给出了全解耦奇偶向量的存在条件,并结合卡尔曼滤波方法得到了故障模型参数的估计方法。最后给出了仿真实例。结果表明,全解耦奇偶方程方法能估计出动态系统执行器故障模型。     

降落伞稳定下降阶段流场的数值模拟

随着计算流体力学 (CFD)和结构有限元分析法的兴起 ,精确描绘降落伞力学特性成为可能。文中采用有限体积法求解不可压流的N -S方程 ,采用Spalart-Allmaras(SA)模型作为湍流模型 ,以模拟降落伞稳定下降阶段时的流场特性。流场的计算结果与实验数据基本一致 ,这表明该方法对降落伞的分析是有效并且可行的     

基于多变量特征模型的机械手自适应控制

机械手系统动力学模型所设计出来的控制器一般具有比较复杂的形式。为了解决机械手系统控制器复杂性问题,推导出机械手系统的线性多变量时变二阶离散差分方程,得到了其方程系数所满足的特征关系,并由分析可知当采样时间间隔足够小时,该方程系数具有慢时变特性。基于建立的差分方程及其系数的特征关系,给出了带有遗忘因子的最小二乘辨识方法和多变量黄金分割控制律相结合的且形式简单的自适应控制算法。利用系数的特征关系,能使需要辨识的参数个数减少。此算法具有只需要关节的位置信息而不需要速度信息的优点。数学仿真结果表明该方法对机械手控制的有效性和鲁棒性。     

守恒双曲型方程组的一类高效差分格式

本文将文［１］所构造的差分格式推广到方程组情形。数值实验结果表明,格式具有高效分辨激波的能力。     

二维抛物化稳定性方程计算

针对二维抛物化稳定性方程(PSE)在法向上采用正交函数展开法进行了计算.在求解PSE的初始条件以及L0的中性曲线中应用了Chebyshev多项式展开.计算结果表明可以获得非常精确的解,证实了在PSE中应用Chebyshev多项式的有效性.     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

非线性边界条件的Burgers－KdV方程的静态解

对几类非线性边值条件，讨论了Ｂｕｒｇｅｒｓ－ＫｄＶ方程的边值问题的静态解的存在性和解的个数．借助不动点定理，证明了该边值问题至少有一个，至多有有限个静态解．     

The use of Gaussian equations of motions of a satellite for local gravity anomaly recovery

The orbital elements of a low Earth orbiting satellite and their velocities can be used for local determination of gravity anomaly. The important issue is to find direct relations among the anomalies and these parameters. Here, a primary theoretical study is presented for this purpose. The Gaussian equations of motion of a satellite are used to develop integral formulas for recovering the gravity anomalies. The behaviour of kernels of the integrals are investigated for a two-month simulated orbit similar to that of the Gravity field and steady-state ocean circulation explorer (GOCE) mission over Fennoscandia. Numerical investigations show that the integral formulas have neither isotropic nor well-behaved kernels. In such a case, gravity anomaly recovery is not successful due to large spatial truncation error of the integral formulas. Reformulation of the problem by combining the orbital elements and their velocities leads to an integral with a well-behaved kernel which is suitable for our purpose. Also based on these combinations some general relations among the orbital elements and their velocities are obtained which can be used for validation of orbital parameters and their velocities.