分析插值误差和斜率的轨迹优化网格细化方法

提出一种基于插值误差和斜率分析的轨迹优化自适应网格细化方法,包括节点插入算法和节点删除算法。节点插入算法分析各个离散节点的控制变量的插值误差。若插值误差较大,则在该节点周围增加节点细化网格;否则,不进行细化。节点删除算法分析各个离散节点处的控制变量斜率。若某个节点的左斜率和右斜率都为零,那么删除该节点;否则,保留该节点。采用三个典型的轨迹优化算例验证了所提出的方法的有效性和特色,并且与其它几种网格细化方法进行了对比。仿真结果表明,本文方法生成的网格规模较小,需要的网格迭代次数较少,能够快速、高精度求解非光滑轨迹优化问题。     

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

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

A Parallel Adaptive 3D MHD Scheme for Modeling Coronal and Solar Wind Plasma Flows

A parallel adaptive mesh refinement (AMR) scheme is described for solving the governing equations of ideal magnetohydrodynamics (MHD) in three space dimensions. This solution algorithm makes use of modern finite-volume numerical methodology to provide a combination of high solution accuracy and computational robustness. Efficient and scalable implementations of the method have been developed for massively parallel computer architectures and high performance achieved. Numerical results are discussed for a simplified model of the initiation and evolution of coronal mass ejections (CMEs) in the inner heliosphere. The results demonstrate the potential of this numerical tool for enhancing our understanding of coronal and solar wind plasma processes. This revised version was published online in June 2006 with corrections to the Cover Date.     

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

固体火箭的鲁棒自适应耗尽关机制导方法研究

针对固体火箭的耗尽关机制导问题,提出了一种基于闭路制导和零射程线相结合的自适应耗尽关机制导方法,该方法首先利用闭路制导技术,使待增速度快速收敛;随后基于零射程线理论,设计了对关机时间变化不敏感的零射程线能量耗散管理导引律和末速精调导引律,提高了耗尽关机方法在推力、秒耗、关机时间随机摄动条件下的鲁棒性和自适应性以及制导精度,从而拓展了耗尽关机显示制导方法的工程化应用.     

火箭返回制导动力着陆段的自适应启动方法

针对重复使用火箭垂直返回着陆问题,提出了一种燃料最优的动力着陆段自适应启动方法。首先,将燃料最优启动点对应的动力着陆轨迹的推力剖面和攻角剖面描述为解析的形式,该解析形式中攻角剖面由状态量唯一确定,推力剖面仅含一个待定参数。随后,通过预测具有上述解析剖面形式的轨迹判断是否启动动力着陆。计算中引入松弛终端位置约束的策略求解推力剖面待定参数,由终端位置约束判断是否满足燃料最优启动条件。上述策略将燃料最优启动条件的判断问题简化为单一变量求解问题,实现了该问题的快速求解。仿真结果表明,该方法得到的启动点与数值优化方法得到的燃料最优启动点接近,且求解过程稳定、计算效率高。     

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.     

具有多输入和扰动作用的空空导弹滚动通道的前馈/反馈自适应控制器

中程空空导弹的参数处于激烈、大范围的时变中，采用古典控制器很难满足不同高度弹道的性能规范。为了使导弹在飞行高度和速度大范围变化时均能保持优良性能，采用数字自适应控制方法是合适的。本文提出了一种基于超稳定性理论的飞行控制系统方案，其性能明显优于不采用自适应控制器的古典系统。     

卫星ATM网络中自适应FEC情况下的有效带宽估计

针对采用自适应前向差错控制（AFEC）的卫星异步传递模式（ATM）网络，研究了两种有效带宽估计方法。由于AFEC码率的动态特性使业务实际的传输速率具有不固定性，因此估计有效带宽时要考虑AFEC的影响。首先针对一般卫星时变信道和多码率AFEC情况下的卫星ATM，应用流体近似方法分析了有效带宽的计算，然而该方法要求的条件及计算复杂度均比较高，应用中有局限性。以此方法为基础，提出了采用修正因子的估计方法。数值结果表明采用修正因子的有效带宽估计方法比较简单可靠，适用面较广。     

高阶对象的二阶模型匹配

本文证明了对于一类高阶对象，其输出可以用一系列系数连续变化的线性定常二阶系统的输出组合完全代替。