利用并行八叉树划分技术开展复杂构型航天器放气污染分析研究

放气污染对航天器造成的有害影响需要定量评估。放气产物分子的输运过程可以采用Monte-Carlo方法模拟,并用射线表示放气分子输运轨迹。如此一来,即转换为射线追踪问题。由于航天器构型复杂,射线追踪的计算很耗时,所以有必要研究其加速算法。文章将复杂构型航天器表面用三角形非结构网格表示,利用八叉树划分技术开展分子运动轨迹与表面相交的加速模拟计算。同时,在基于共享内存多核计算平台,利用OpenMPAPI软件实现了计算的并行化。该方法在对空间站的实例计算中显著提高了运算速度,并对空间站的污染情况进行了初步评价分析。     

二维跨音速Euler方程分区并行算法

对二维跨音速流动的Euler方程分区算法、并行算法以及多区计算的有效内边界耦合条件进行了探讨,应用Van Leer矢通量分裂方法和一维数组方式,研究了多种区域分解数目的分区计算效率.并行计算中采用"先进先出"的同步控制等待机制,采用纯结点并行编程方式进行了单区、二区和四区并行计算对比,分析了影响并行效率和通讯过载比的主要因素.     

共轴刚性旋翼直升机旋翼控制相位角问题分析

利用共轴刚性旋翼直升机飞行动力学模型,以XH-59A共轴刚性旋翼直升机为研究对象,分析了旋翼控制相位角对纵向配平特性、需用功率以及上、下旋翼桨毂弯矩的影响。基于分析结果,提出了一种针对共轴刚性旋翼直升机的旋翼控制相位角的配置方法。该配置方法以降低直升机需用功率为目标,并保证上、下旋翼桨毂弯矩和配平特性满足要求。通过该方法能使XH-59A直升机在0~80 m/s的飞行速度范围内满足上、下旋翼最大桨毂弯矩和纵向操纵限幅的要求,并且能最多降低8%的直升机需用功率,为共轴刚性旋翼直升机的设计提供了参考依据。     

轴向磁通磁场调制型电机的解析分析

提出了一种新型二维解析分析方法用于计算轴向磁通磁场调制型电机的气隙磁场分布和电机特性。对电机每个子域列出拉普拉斯方程或泊松方程并结合边界条件,求解出每个子域的矢量磁位,进而得到气隙磁场分布、空载反电势、电磁转矩和轴向磁拉力的解析表达式。为了评估解析方法的计算精度,对二维解析分析结果和三维有限元仿真结果进行了比较,结果表明二维解析计算结果的计算误差在10.3%以内。由于解析分析在保证计算精度的前提下大大缩短了计算时间,因此该方法适用于轴向磁通磁场调制型电机初始设计阶段电机参数的确定。     

流场计算程序并行化方法的研究

利用“流场分区”的方法把流场计算程序并行化,并成功地在Transputer并行计算机上运行,并行效率可达到80％。介绍了流场划分、并行编程以及提高并行效率的方法。     

温度比对二维受限垂直射流穿透的影响

用数值方法研究了温度比对受限垂直射流穿透深度的影响。发现进口温度比是独立于动量比以外的重要影响参数。在一定的动量比下,冷射流垂直射入高温来流中,受热膨胀加速,动量较等温射流增加,从而穿透深度增加;热射流垂直穿入冷主流中,受冷却而减速,动量减少,导致穿透深度减小。      

振动翼型跨音速非定常Euler方程数值解方法

通过坐标变换使变化的物面边界条件及无穷远边界条件变为固定的有限范围内的边界条件。因此可在不增加计算机时的情况下,边界条件可以精确地满足。由于主控方程是Euler方程,因此,凡适合于Euler方程解的问题（跨音速、大扰动、非定常、高简化频率）均可应用此法。同时,在计算过程中采用贴体、正交网格与可控网格技术获得最佳网格分布。算例计算结果是令人满意的。     

基于空间算子代数的空间多体系统动力学递推计算

应用空间算子代数方法对空间多体系统进行动力学建模、分析与仿真。空间算子代 数 作为基于李群及旋量方法的算子理论,采用估计理论中卡尔曼滤波及平滑方法,实现了具有 明确物理含义的O(n)计算效率的动力学实时递推求解算法。首先使用空间算子代数理论 对多体系统动力学进行递推形式的描述、建模及分析,并通过与估计理论对比提出了质量矩 阵求逆的高效递推算法,而后基于面向对象技术编制了软件,并通过漂浮基双臂机器人典型 算例的求解,与商业动力学软件Simpack进行对比,验证了软件及算法的正确。结果表明, 空间算子代数方法概念清晰、推导方便,计算效率及精度能够满足工程需求,可用于空间多 体系统动力学实时仿真及控制系统设计。〖JP〗     

耦合高效配平策略的旋翼气动特性分析方法

建立了一个耦合高效配平策略和计算流体力学（CFD）技术的旋翼气动特性分析方法.采用CFD求解器进行雅克比矩阵和旋翼流场的计算,建立了一个包含直接配平法的旋翼气动特性分析平台;引入差量法思想对配平策略进行简化,发展了一种差量配平法,解决了直接配平法中雅克比矩阵计算效率低下的问题,显著缩短了计算时间.旋翼流场求解中控制方程采用Navier-Stokes方程,空间离散格式选取2阶迎风Roe格式,时间推进采用隐式lower-upper symmetric Gauss-Seidel (LU-SGS)格式,湍流模型为Baldwin-Lomax (B-L)模型.分别采用建立的基于直接配平方法和差量配平方法的气动特性分析平台对AH-1G和Helishape 7A直升机的旋翼进行了配平计算和研究,对比分析了两种配平方法的计算效率和精度及背景网格的网格量对配平结果的影响.结果表明:与直接配平法相比,耦合差量配平法的气动特性分析方法计算圈数更少,可节省至少60%的计算时间;使用疏化的背景网格虽然不利于捕捉流场细节,但是可以在兼顾配平结果准确性的同时,降低配平计算代价,提高旋翼配平和流场计算的效率.      

飞行器气动/结构多学科延迟耦合伴随系统数值研究

基于自主研发的大规模并行结构化网格雷诺平均Navier-Stokes（RANS）求解器PMB3D以及流固耦合代码FSC3D建立了飞行器静气动弹性数值模拟技术,进一步基于并行化伴随方程求解器PADJ3D,开展了Navier-Stokes方程与结构动力学方程耦合离散伴随的推导、构造。对各个学科伴随变量进行延迟处理,进行耦合伴随系统的解耦,学科之间的影响通过右端强迫项的形式在方程中体现,通过松耦合形式进行各个学科方程右端项数据传递,各个学科的伴随方程一定程度上能够相对独立,进一步实现了基于LDLT方法的结构伴随方程的高效求解;对典型客机柔性机翼进行梯度信息求解,并与考虑气动弹性影响的差分结果进行对比分析。数值模拟结果表明,延迟耦合伴随形式更有利于保持原有程序结构形式以及程序模块化,梯度计算精度完全满足优化设计需要,能够为柔性机翼飞行器气动/结构多学科优化设计提供研究基础与技术平台。