粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

不可压扰动方程高精度对称紧致差分数值解法及应用

从粘性不可压扰动方程一阶改型形式出发,对其实现了高精度对称紧致差分离散,就导出的扰动线性特征值问题给出了一个高铲双重迭代局部解法,以相同精度将特征值和特征函数同时得到。通过不可平面Ｐｏｉｓｅｕｉｌｌｅ流时间稳定性算例详细对比显示了算法良好的谱分辨能力和较弱的网格依赖性,并结合复矩阵广义特征隐式移ＱＺ算法获取了一个 扰动特征值谱计算结果。     

Navier-Stokes方程积分表达式的有限元解法

本文采用定常不可压Navier-Stokes（NS）方程的积分表达式方法计算粘性流动问题,用线性三角形单元划分计算区域,建立了适当的迭代求解过程。本文改进了压力计算方法,用全变量松驰加快了收敛速度,用矩阵分块计算减小了内存限制,并可以通过对矩阵特性的分析更好地控制迭代过程的稳定性。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

A policy iteration method for improving robot assembly trajectory efficiency

Bolt assembly by robots is a vital and difficult task for replacing astronauts in extra-vehicular activities (EVA), but the trajectory efficiency still needs to be improved during the wrench insertion into hex hole of bolt. In this paper, a policy iteration method based on reinforcement learning (RL) is proposed, by which the problem of trajectory efficiency improvement is constructed as an issue of RL-based objective optimization. Firstly, the projection relation between raw data and state-action space is established, and then a policy iteration initialization method is designed based on the projection to provide the initialization policy for iteration. Policy iteration based on the protective policy is applied to continuously evaluating and optimizing the action-value function of all state-action pairs till the convergence is obtained. To verify the feasibility and effectiveness of the proposed method, a noncontact demonstration experiment with human supervision is performed. Experimental results show that the initialization policy and the generated policy can be obtained by the policy iteration method in a limited number of demonstrations. A comparison between the experiments with two different assembly tolerances shows that the convergent generated policy possesses higher trajectory efficiency than the conservative one. In addition, this method can ensure safety during the training process and improve utilization efficiency of demonstration data.     

车模比例和风洞截面对轿车气动性能影响

采用三维不可压缩N-S(Navier-Stokes)方程和RNG k-ε(renormalization group k-ε)湍流模型仿真计算汽车流场特性.控制界面的物理量应用二阶迎风差分格式获得,并运用SIMPLEC压力修正法进行迭代.选用不同车模比例和风洞截面形状进行仿真计算,分析了车模比例和风洞截面形状对轿车气动性能的影响.结果表明:车模比例在一定范围内,汽车气动参数变化平缓;选用不同的车模比例在三种典型风洞模型中进行仿真计算,结果相差很大.      

基于NURBS和梯度法的飞艇形状优化

采用非一致有理样条曲线(NURBS)作为飞艇外形参数化方法、基于梯度的优化方法、飞艇体积作为约束条件、飞艇的阻力系数作为目标函数的优化流程对平流层流动条件下的飞艇外形进行了优化.阻力系数通过求解二维黏性不可压缩流场获得,阻力不但包含了形状阻力也包括了黏性阻力.优化过程以LOTTE飞艇的形状为初始形状,优化的结果表明优化后飞艇的气动性能明显改善,阻力系数有了比较明显的降低.优化的结果也同时证明了该方法适用于飞艇的优化命题.      

A boundary value problem of aerohydrodynamics in designing an axisymmetric body with jet blowout

A numerical and analytical method for constructing a shape of the axisymmetric body streamlined with jet blowout by the specified velocity distribution along its meridian section is proposed. The foundation of the method is the iteration process based on the solutions of an inverse problem in the plane case and a primal problem for the axisymmetric body. A program realizing the iteration process is set up and examples of numerical calculations are given.     

边界元的耦合界面步进法分析物体撞水响应

计算物体的撞水响应目前已有了一些专用的算法。在分析和比较这些算法的基础上,提出了一个解撞水问题的边界元耦合界面步进法,该方法不仅可充分发挥边界元计算半空间流场的优越性,更主要是在加速收敛方面有其独到之处。给出圆柱刚体和楔形刚体两个撞水算例,结果表明了该方法的可靠性和有效性。     

基于ADE-ELM的涡轴发动机建模方法

提出了基于自适应微分进化-极端学习机(ADE-ELM)求解平衡方程的高精度涡轴发动机实时部件级模型建立方法.基于牛顿-拉夫逊(N-R)迭代模型,以迭代计算前模型平衡方程残差为输入,迭代收敛后平衡方程猜值修正量为输出,训练极端学习机,并采用自适应微分进化(ADE)算法优化极端学习机(ELM)参数,提高猜值修正量映射精度.ADE算法中采用sigmoid型自适应缩放因子,提高了微分进化算法的寻优能力.在涡轴发动机不同飞行状态下的测试结果表明,以N-R迭代算法模型为基准,基于ADE-ELM的发动机模型,最大建模误差约为一次通过算法的1/3,运算耗时约为一次通过算法的1/3,验证了算法的有效性.      

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

非定常流动计算的混合时间推进方法研究

结合ADI方法的二阶时间精度特性和高斯-赛德尔迭代法的良好收敛性特性,提出了一种ADI-高斯-赛德尔混合迭代格式,应用该格式,对NACA0012翼型的跨声速俯仰振荡进行了EULER方程求解,空间离散格式采用ROE’S FDS格式,计算结果与实验及成熟的LUSGS双时间步迭代方法计算结果进行了分析比较。计算结果表明,该格式精度能达到LUSGS双时间步迭代方法完全收敛的二阶精度,且该格式精度对内迭代次数的依赖性较低。分析表明,本文提出的方法具有较高的精度和较好的收敛性。     

高级曲面设计功能的算法原理及实现

论述了国产化软件BSURF-GI交互式曲面造型系统中高级曲面设计功能的算法原理及其实现途径。提出了一种能沿曲面与曲面交线和过渡曲面边界线进行曲面裁剪的算法,指出了该系统今后进一步研究的方向及与实体造型系统集成的途径。为解决曲面与曲面求交算法中自动给出Newton-Raphson迭代法初值的优选难问题,本文还提出了一种有效的自动优选初值的方法——面积映射法。最后给出了一些高级曲面功能图例。     

一种不需要迭代的发动机辅助变量建模方法

现代航空发动机建模计算时常用牛顿迭代法求解非线性方程组,但由于发动机模型非线性太强且十分复杂而不易在全包线收敛。本文针对这一不足之处提出了一种不需要迭代就可以求解此非线性方程组的方法:引入航空发动机容积动力学中的变量及方程作为辅助变量及辅助方程,从而使整个非线性方程组闭合而可以直接求解。本文将此方法成功应用于某涡扇发动机,取得了良好的效果,验证了本方法的可行性。      

面向突发任务的空间站任务重规划方法

针对空间站在轨运营出现突发任务的情况,提出一种基于启发式规则的任务重规划方法,满足了方案重规划的快速响应需求。根据任务执行的连续性特点和冲突状态,建立了空间站突发任务规划领域模型。考虑重规划过程中任务包含活动间复杂约束关系传播的影响,提出了时间回溯迭代冲突化解策略,同时依据任务执行时间间隔,提出针对间隔插空的时间冗余启发式规则。基于时间回溯迭代冲突化解策略和时间冗余启发式规则,对原任务执行计划进行实时重规划,实现了突发任务的快速响应。应用算例分析表明,提出的重规划方法可以成功地满足空间站突发任务规划需求,实现实时更新空间站在轨任务执行详单的目的。     

可压流双方程紊流模型及其应用

 本文提出了可压流紊流模型。所用的方法类似于K-KL双方程模型,但考虑了气流压缩性的影响。方程中的常数已由实验确定。用本模型计算了管内超音速混合流动,其理论计算曲线与实验数据相当符合,证明该模型较不可压紊流模型更为适用。     

一维ADI-FDTD方法的数值色散分析

给出了交替变向隐式时域有限差分法(ADI-FDTD)的一维差分格式，利用增长矩阵推导了数值色散关系式。将其与过去研究中给出的色散关系式进行了比较。利用牛顿迭代法对其进行了全面的分析，发现了时间步长的选取与网格、色散误差的大小有一定的联系。     

固体推进剂药柱松弛模量随机粘弹性有限元分析

为了研究性能参数的分布对药柱结构完整性的影响，以不可压缩粘弹增量有限元和摄动法为基础，采用局部平均方法对随机场进行离散，发展了一种松弛模量随机粘弹有限元方法，讨论了模量的随机性对结构分析的影响。对药柱进行随机模拟，程序实施简单，计算效率和精度较高。     

共轴异向旋转筒间不可压流的临界主不稳定性

基于长隙比无穷大和Tollmien-Schilchting 行谐波假设，通过导出坐标系下不可压一阶改型边疆特征值问题，实现了异向旋转筒流动临界主不稳定性的线性分析，在稳定性方程离散中引入了一类高精度对称紧致差分格式，在特征值搜索中发展了一个高效率复矩阵双重反迭代算法。双临界点分析揭示了几何参数、、流动参数、流动参数、扰动参数相互间的变化规律。理论结果同实验测量数据基本一致。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。