弹道导弹中段机动变轨模型研究

鉴于弹道导弹在中段机动变轨会导致飞行轨迹和落点的改变,对弹道导弹中段机动变轨进行了研究.针对弹道导弹中段飞行,引入小火箭的推力作用,建立弹道导弹中段变轨运动模型,运用四阶Runge-Kutta积分法求解其飞行轨迹,并将两种时机的机动变轨结果与自由飞行情况进行了比较,变轨时的中段结束点位置改变较为明显.仿真结果表明了该模型的有效性.该变轨模型可为目标的跟踪研究提供依据.     

Characterizing Strong Stability Preserving Additive Runge-Kutta Methods

Space discretization of some time dependent partial differential equations give rise to ordinary differential equations containing additive terms with different stiffness properties. In these situations, additive Runge-Kutta (additive RK) methods are used. For additive RK methods the curve of absolute monotonicity gives stepsize restrictions for monotonicity. Necessary conditions for nontrivial curves of absolute monotonicity are the nonnegativity of the additive RK coefficients and some inequalities on some incidence matrices. In this paper we characterize strong stability preserving additive Runge-Kutta methods giving some order barriers and structural properties. Research supported by the Ministerio de Educación y Ciencia, Project MTM2005-03894.     

An efficient solver for the RANS equations and a one-equation turbulence model

A three-stage Runge-Kutta (RK) scheme with multigrid and an implicit preconditioner has been shown to be an effective solver for the fluid dynamic equations. Using the algebraic turbulence model of Baldwin and Lomax, this scheme has been used to solve the compressible Reynolds-averaged Navier–Stokes (RANS) equations for transonic and low-speed flows. In this paper we focus on the convergence of the RK/Implicit scheme when the effects of turbulence are represented by the one-equation model of Spalart and Allmaras. With the present scheme the RANS equations and the partial differential equation of the turbulence model are solved in a loosely coupled manner. This approach allows the convergence behavior of each system to be examined. Point symmetric Gauss-Seidel supplemented with local line relaxation is used to approximate the inverse of the implicit operator of the RANS solver. To solve the turbulence equation we consider three alternative methods: diagonally dominant alternating direction implicit (DDADI), symmetric line Gauss-Seidel (SLGS), and a two-stage RK scheme with implicit preconditioning. Computational results are presented for airfoil flows, and comparisons are made with experimental data. We demonstrate that the two-dimensional RANS equations and a transport-type equation for turbulence modeling can be efficiently solved with an indirectly coupled algorithm that uses RK/Implicit schemes.     

Initial value methods for the primary two point boundary value problems in modeling viscoelastic flows

Initial value solvers through shooting methods will be used to compute solutions to the primary two-point boundary value problem arising in modeling viscoelastic flow. We will show that the classical Runge-Kutta methods will have at least h ²-order of convergence due to the presence of the singularity. Comparison with the work of others will also be presented through some numerical results and examples.     

时间最优开关控制的非线性规划方法

提出单输入非线性时间最优开关控制的非线性规划算法．针对不同的开关控制划分未知时间段，以时间段为变量建立与最优控制等价的非线性规划模型．构造每个未知时间段的等分龙格库塔格式，该格式不但不增加未知变量，而且可构造出不含导数的非线性规划算法．利用龙格库塔格式的收敛性和非线性规划的一阶最优性条件证明该方法的收敛性．最后以实例验证该算法．     

Note on filon-type integration for higher order exponential time differencing methods in stiff systems

The Filon-type quadrature is efficient for highly oscillatory functions — Fourier transforms. Based on Cox and Matthews’ ETD schemes, the higher order single step exponential time differencing schemes are presented based on the Filon-type integration and the A-stability of the two-order Adams-Bashforth exponential time differencing scheme is considered. The effectiveness and accuracy of the schemes is tested. Foundation item: Projects(02JJY2006, 03JJY2001) supported by Natural Science Foundation of Hunan Province; project supported by JSPS Fellowship Research Program     

非线性时滞控制系统Runge-Kutta方法的IS稳定性

研究求解非线性时滞控制系统的Runge-Kutta方法的IS稳定性,首先给出非线性时滞控制系统理论解IS稳定的充分条件,在约束网格下,将(k,l)为代数稳定的Runge-Kutta方法求解该系统,获得了方法IS稳定的条件.     

改进的油浸式电力变压器绕组热点温度计算和在线监测及运行分析

针对变压器绕组热点温度的计算模型,采用变步长的Runge-Kutta方法对计算模型求解,实时计算变压器的顶层油温和绕组最热点温度。建立针对变压器的实时在线监测系统,实现变压器运行情况的监测分析,采用先进的计算方法对变压器热点温度计算模型求解,提高计算精度;并对变压器短期过负载能力进行科学合理预测。     

龙格库塔方法在求解瞬态中子扩散方程中的应用

本工作开发了NGFMN-K程序求解瞬态中子扩散方程。空间离散采用第2类边界条件节块格林函数方法,时间离散分别选取向后欧拉格式和四阶精度对角线隐式龙格库塔(DIRK)格式。对DIRK格式采用嵌入三阶精度格式估计截断误差实现时间步长自调节。数值验证结果表明,两种格式的计算结果与参考程序结果符合很好,对于剧烈瞬变情况,DIRK格式较向后欧拉格式更为精确、高效。     

用经典R-K法求解扩散方程

扩散方程,通常是带有初值-边值的时倚偏微分方程.一般数值求解这一类偏微分方程是将偏微分方程转化为初值问题的常微分方程,再对常微分方程进行求解.目前常用差分法、有限体积法或有限元法来求解方程,但在处理实际问题时,这些方法往往因为解的精度不高而使计算结果不合理.用经典R-K方法求解扩散方程,可以明显提高解的精度,通过对实际算例计算结果的比较,该方法解的精度几乎与解析解的精度相同.