三维凸域中候选卫星点集的确定算法

局部网格生成是无缝有限元并行计算的关键问题,然而,不恰当的局部网格生成算法会产生不一致网格.借助Delaunay三角剖分的唯一性,指出了三维凸域中Delaunay网格局部生成的不一致性,分析其产生的原因,建立避免局部网格不一致性的关键算法--搜索候选卫星点集的优化探索球算法,为网格并行生成做好数据准备.并从理论上证明了任意节点的优化探索球域包含其所有卫星点,从而,保证全域网格没有不一致性.通过算例分析了算法的正确性和局部性.     

基于切割单元直角坐标网格的DSMC方法及优化技术

对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。     

航天器燃耗最优轨道直接/间接混合法延拓求解

针对转移时间和始末状态固定的航天器燃耗最优轨道的求解,给出了一种延拓方法:以双脉冲轨道为初值,首先求解全程推进轨道,然后逐步增加推力幅值,应用直接/间接混合法依次求解所有推力幅值下的、满足包括开关函数在内的所有必要条件的转移轨道,包括连续和脉冲推力轨道。通过基于开关函数曲线变化趋势的开关序列预设方法,以及基于已有优化结果的延拓步长自适应方案,实现了延拓方法的自动运行。为实现该延拓方法,给出了适用于改进春分点根数模型的脉冲最优转移轨道主矢量必要条件,推导了无推力轨道段改进春分点根数协态变量状态转移矩阵。通过3个算例对延拓求解会遇到的不同情况进行了具体说明。延拓方法可以看作现有直接/间接混合法的进一步完善与拓展,延拓过程和结果有助于对燃耗最优轨道与系统参数之间的关联获得更为深刻的认识。     

无网格方法关于导数计算的程序验证

为了提高无网格方法对空间导数的数值计算能力,提出并论证了新的适用于直角坐标系、柱坐标系和球坐标系的空间导数无网格算法。基于数值解与理论解的对比分析,对数值误差和收敛性进行了后验评估。评估结果表明:该算法对函数和导数的估计均为2阶精度,支撑域尺度对数值误差的大小有一定影响,但不影响数值计算的精度等级。在离散尺度h为0.01的条件下,对所选函数及其导数的数值计算相对误差不大于0.65%。      

立式捏合机桨叶结构参数对混合釜流场影响的仿真分析

为揭示立式捏合机桨叶几何结构参数与混合釜流场之间的关系,以1 L两桨立式捏合机为研究对象,采用计算流体力学ANSYS Fluent 14.5软件对立式捏合机进行了三维模型仿真研究,分析了桨叶结构参数(桨桨间隙、桨叶螺旋角)的改变对混合物料挤压应力和剪切应力的影响。结果表明,增大桨桨间隙可减小桨叶对混合物料的挤压应力和剪切应力,空心桨螺旋角分别为40°、45°、50°,桨桨间隙由1.0 mm增加到3.0 mm,挤压应力、剪切应力分别平均减小了82.4%、57.23%;增大桨叶螺旋角可增大桨叶对混合物料的挤压应力和剪切应力,空心桨间隙分别为1.0、2.0 mm,桨叶螺旋角由35°增加到55°,挤压应力、剪切应力分别平均增加了92.8%、55.7%。     

一类非线性空间分布系统自适应容错控制

研究了一类由偏微分方程组(Partial differential equations,PDE)描述的非线性空间分布系统的自适应容错控制(Fault-tolerant control,FTC)问题。首先,采用模态分解方法将PDE系统表示为一个有限维慢子系统与一个无限维快子系统相耦合的型式;然后,基于慢子系统模型及小增益定理设计了自适应FTC律,使闭环系统在所有容许的未知非线性动态以及执行器卡死故障的影响下都能渐近稳定。     

混合式CRP吊舱推进器水动力性能数值模拟及试验

建立混合式CRP吊舱推进器数值模型并置于数值空泡水洞内,结合RANS方程和SST k-ω湍流模型,运用滑移网格方法对混合式CRP吊舱推进器在均匀流场中水动力性能进行了数值模拟,得到了敞水性能预报结果,并分析了混合式CRP推进器的非定常力;同时,利用空泡水洞、吊舱动力仪及长轴动力仪对混合式CRP吊舱推进器敞水性能进行了试验研究,通过3组重复性试验表明了试验方法的稳定性;并将3组试验平均值与将数值模拟所得的敞水性能结果进行了比较,彼此吻合较好。结果表明,本文所用试验方法和数值计算方法对吊舱推进器模型敞水特征的预报具有较高的可信度,能达到工程应用要求的要求;和国外学者在空泡水洞内得到的不附带吊舱相同对转桨敞水试验数据进行比较,分析了吊舱对前后桨敞水性能的影响。结果表明,吊舱对前桨的敞水性能基本没有影响,但对后桨的敞水性能影响显著。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。