带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...     

基于网格变形的伴随方法在翼型优化中的应用

进行了基于网格变形的伴随方法在翼型气动优化设计中应用的研究。随着集成了伴随方法的流场求解器的普及和效率的不断提高,采用梯度法进行翼型气动优化设计的主要计算花费逐渐由计算梯度对流场变量的敏感度（或流场敏感度）转变为计算梯度对网格变形的敏感度（或网格敏感度）,后者的求解通常采用有限差分方法,计算花费较高。在传统伴随方法（或流场伴随方法）的基础上,引入了基于网格变形的伴随方法（或网格伴随方法）,采用网格伴随方法计算梯度,可以大幅度减少梯度计算花费,提高翼型气动优化设计的效率。     

复杂构件MDO六面体网格重构方法

从参数化模型和模型拓扑特征出发,提出了一套六面体网格重构的有限元模型建模方法.同时,使用随动点来约束网格的拓扑边界,保证网格的正交性,确保了网格重构后的网格质量;结合现有的商用软件ICEM的网格生成算法,开发了基于特征的参数化有限元(FE)建模程序.最终,针对航空发动机低压涡轮带冠叶片的叶冠参数进行的多学科优化 (MDO),实例验证了网格重构方法的有效性和可行性.      

叶片参数化有限元建模与罩量优化设计

基于网格变形技术提出了一种叶片参数化有限元建模方法,该方法具有快速实现叶片参数化有限元建模的优点,同时能够在保证叶片有限元网格拓扑形式不变的前提下,有效解决叶片优化设计过程中的网格自动更新问题.在此基础上,结合叶片罩量优化设计的灵敏度分析技术,建立了叶片罩量自动优化一般设计流程,开展了叶片实例优化设计应用研究.优化结果表明,相对于初始设计,罩量优化设计可以有效地改善应力分布,降低叶片应力水平.      

基于在轨热环境的可展开天线反射面精度调整技术

随着宇航事业的发展,网状展开天线的尺寸越来越大,对形面精度要求越来越高(特别是高频段),在轨热环境与地面调整点的热环境差别大,必须采用基于在轨热环境的网面精度调整技术以满足反射面的精度要求。就在轨热环境下的辅助牵引面式网面精度调整技术的理论分析及调整方法进行了探讨,给出了计及在轨热环境的网面形状精度调整的基本理论公式及调整步骤,并对某样机进行了调整计算。计算结果表明,文中的分析、公式、方法是正确的,可行的。     

外挂物投放问题的二维非结构动态网格技术研究

以动态网格守恒型方程为控制方程,在二维非结构网格基础上发展了弹簧式光顺法和局部网格重构法两种动态网格变形技术。结合运动轨迹计算方程,将动态网格法应用于飞机外挂物投放问题中,成功模拟了导弹自由投放的整个过程,获得了详细的包括弹体运动速度、运动轨迹和受力情况等在内的丰富的流场信息。成功的应用不仅验证了动态网格法的可行性;而且也表明,动态网格法的进一步完善与发展将成为解决外挂物投放问题的有力工具。     

网格拓扑对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.