用泰勒级数求解非线性代数和微分方程组

提出一种求解非线性代数方程和非线性常微分方程的新方法。所论方程组被转换成增量形式，而解被表示成泰勒级数。在代数方程（多项式形式）的情况下，此算法归结为求解一系列线性方程组，而且在每一增量步中系数矩阵是不变的；在常微分方程组的情况下，此算法归结为一组递归算式，并不要求解方程组。此方法固有的自动走步功能，可保证得到收敛的解，从而大大节省计算时间。     

基于广义正则化最小二乘的天基测向初定轨

传统天基测向初定轨的不足,主要是由于观测数据存在系统误差和观测方程组的系数矩阵病态或不可逆。文章建立观测方程的半参数回归模型,提出基于补偿最小二乘估计和岭估计的广义正则化最小二乘估计,推导了估计公式,并证明了相关统计性质。引入选主列Givens-QR分解算法,提高观测方程求解效率和数值计算稳定性。仿真结果表明:该方法应用于天基测向初定轨可行,可以提高定轨精度和解算成功率。     

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式，具有十分优异的分辨率特性。但在求解多维方程组时，HCW-R格式需要求解块状三对角方程组，因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P，构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式，其可沿迎风方向、由边界处向内推进求解，避免了处理三对角或块状三对角方程组，从而其计算代价与显式格式无异。虽然就分辨率而言，HCW-E稍逊于HCW-R，但前者的计算效率要显著高于后者。因此，当花费相同的计算代价，HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

柔性杆子杆模型参数的遗传算法优化求解

研究了应用实数型的遗传算法解决柔性机器人子杆法建模中,子杆模型参数的优化求解问题.子杆模型参数形成高度非线性的方程组,它的求解是子杆法能否实现的关键.普通的牛顿迭代法及其变形,要求初始向量足够靠近解,一些所谓大范围收敛的算法又有方程组的Jacobian矩阵非奇异等要求,在此都难以应用.本文基于遗传算法的解法可求得模型参数的最优解,实例仿真显示了这一方法的可行性和高效率.      

三站点备件供应保障关系建模与分析

就建模层面而言备件相关的现有研究在考虑横向保障时存在一定的不足。本文选取三站点组成的备件供应保障系统作为研究对象,通过引入站点间的备件供应保障关系实现了将站点间备件纵向保障和横向保障的综合考虑。运用定量化的概率值表征站点间的备件需求关系,基于库存平衡理论建立了应得备件数量、初始库存量、可用库存量和备件短缺数之间的关系表达式。在建立单站点库存平衡方程的基础上,又依据三站点间的供应保障关系联立得到了三站点备件延期交货量方程组。对方程组解的情况进行了分析,在确定其解唯一性的基础上基于压缩映射原理给出了方程组的迭代求解方法。对于简单示例进行了计算,验证了模型的正确性和求解方法的有效性。     

吴方法在平面并联机构位置正解中的应用

采用吴方法对平面并联机构位置正解问题进行了研究.吴方法是一种求解非线性方程组的数学机械化方法,采用这种方法任何非线性方程组都可以在有限步内得到解决.在给出了吴方法基本原理的基础上,对本问题进行了求解,并将原始方程组转化成为一个三角化的方程组.其中单变量方程的次数为6次,说明平面并联机构可以有6个不同的位姿.最后用数值实例进行了验证,给出了计算结果.吴方法在这一问题中的应用,为求解其它机构学难题提供了新途径.     

子午面太阳风结构及其对日冕抛射物运动的影响

本文提出了一种求解二维理想磁流体力学方程组的组合格式-网格方法。该方法在保证一定精度的前提下稳定性好,计算量小。文中使用该方法算得了子午面太阳风的盔状流动解。以该解作为初态,研究了盔状冕旒位形对日冕抛射物运动的影响,发现该位形使得抛射物向低纬赤道方向偏转,较好地解释了观测结果。     

粘性流体换热问题用直角坐标网格的贴体解法

采用直角坐标网格但通过在边界附近保留不规则控制体的非结构化直角坐标网格方法求解粘性不可压流体流动和换热问题.在边界附近切割单元的求解采用常规的方法处理,而不采用插值方法处理,并在同位网格上对控制方程组进行离散.因而该算法是贴体的,并可适用于解决复杂边界流动换热问题.通过对同轴环行空间和二维倾斜空腔内的对流换热问题的计算对该数值方法在传热问题中的应用进行验证.对比结果表明,计算结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法可用于传热问题.      

函数传递法求解悬链线问题

悬链线问题是一类经典又多变的力学问题，其曲线构形指导着桥梁等工程应用的结构设计。为了求得结构特征或载荷特征不同的悬链线的变形，提出了一种基于传递矩阵法思想的、通用性较好的求解方法。首先，提炼出悬链力学模型，将悬链顺序划分成若干简单单元，结合单元力平衡和本构-几何关系解得单元特征函数方程组，顺序嵌套单元特征函数方程组获得整体特征函数方程组；然后，使用离散Newton迭代法对该非线性方程组进行求解，获得悬链的受力和变形；最后，算例验证了结果与解析解的一致性。函数传递法对具有复杂结构特征和载荷特征的悬链线问题有很好的适用性，对求解其他可划分为若干首尾相接结构单元的结构系统的广义变形也适用。      

基于非Favre过滤的大涡模拟

将直接过滤的可压缩Navier-Stokes(NS)方程组作为可压缩湍流大涡模拟方程组.方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,得到可压缩张量扩散模型.与传统的采用质量加权过滤(Favre过滤)的可压缩大涡模拟方程组相比,由于这种方法没有密度信息的损失,所以能反映真实的密度脉动信息;方程所求量为流场物理量的过滤值,物理意义明确,使大涡模拟结果与直接数值模拟或实验结果更具有可比性.用这种大涡模拟方程组对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验、Baldwin-Lomax(BL)模型及Smagorinsky亚格子模型的结果进行了对比.      

半刚性尾涡在冲压空气涡轮性能计算中的应用

与Goldstein假设一致的刚性尾涡是螺旋桨升力线(面)计算中的常用模型,适于轻负荷工况,但在较大负荷工况时则有一定的误差.与Theodorsen假设一致的半刚性尾涡将Goldstein尾涡推广至较大负荷工况.对Theodorsen尾涡模型应用于冲压空气涡轮升力线计算进行了探讨,对一个模型冲压空气涡轮作了计算,给出了涡距、诱导速度、环量及叶素风能利用系数,其分布趋势合理,风能利用系数与实验值符合的较好.     

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.     

垂直下降状态下的旋翼三维流场数值模拟

应用激励盘模型,采用计算流体力学(CFD, Computational Fluid Dynamics)方法结合动量叶素理论,计算了旋翼在垂直下降时的流场特征与气动性能.计算过程中,旋翼载荷被描绘成沿桨盘径向分布、与当地流动参数相关的压力源项.在直角坐标系中,运用有限体积方法直接求解了定常不可压N-S(Navier-Stokes)方程,湍流模型为S-A(Spalart-Allmaras)一方程模型,重点使用诱导速度经验公式计算出了压力源项.旋翼流场模拟结果和性能预计结果同实验测量数据吻合较好,表明这种CFD方法可以有效地分析旋翼下降飞行时的空气流动特性,为进一步的涡环状态研究提供了参考.     

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.      

旋翼涡尾流与下洗流场的计算方法

 以重点考察旋翼尖涡的运动为目标来构想一种悬停和前飞旋翼涡尾流的理论计算方法.首先从预定的广义尾流入手,并引入尖涡涡核作用的半经验修正,对旋翼的自由尾流进行计算分析.然后在环量收敛与尾迹收敛的条件下,计算拉力以满足拉力收敛准则.之后,作为本方法的工程应用,文中给出了计算的旋翼沿火箭发射线的下洗速度分布并作了简单分析.从理论上说,文中所发展的方法可以为进一步研究尖涡的生成模式与涡核结构所应用,同样也可以预测旋翼的涡诱导流场.     

后掠翼身干扰区流动显示

在水槽中利用激光片光源及荧光素纳染色液显示方法,研究了圆柱、机翼与平板交接区及后掠圆柱、后掠机翼变迎角情况下的干扰流场结构、特性及参数影响规律.实验结果表明,除Re数之外,模型迎角、后掠角等参数对干扰区马蹄涡特性有很大影响.研究发现后掠圆柱及后掠机翼在一定条件下存在一类既不同于马蹄涡也不同于卡门涡的空间稳定发展的三维旋涡系即背涡.背涡在一定的迎角及后掠角条件下发生破裂;机翼背涡随迎角变化具有不同于圆柱背涡的特点.讨论了干扰背涡产生的机制及其与马蹄涡的相互关系.此外还给出了对于干扰流场典型截面的PIV测量结果.     

一种战斗机外形涡流场Euler方程数值模拟

 用三步显式格式时间推进求解物理空间曲线网格上有限体积离散的Euler方程,数值模拟一种战斗机外形的涡流场.计算的机翼表面压力分布与实验符合较好,用Euler方程捕获到机翼前缘分离涡(主涡)及翼面上的二次分离涡.Euler方程解中出现的二次涡可能是逆压梯度、尖锐边缘和人工粘性共同作用的结果,它的流谱与实验定性符合.     

悬停状态旋翼固定尾迹分析中三角级数的应用

为了在固定尾迹分析中更好地描述桨叶附着涡环量分布,结合叶素理论,采用一个三角级数表达附着涡环量的连续分布,并应用于悬停状态下单旋翼和共轴双旋翼的气动特性分析,通过小量计算即获得连续光滑的附着涡环量分布,将其应用到推导的公式中计算出连续分布的诱导速度,进而分析了共轴双旋翼上下旋翼气动干扰及扭矩平衡条件.计算结果与实验结果进行了比较分析,验证了该方法的有效性,并得出结论:下旋翼总距变化对上旋翼诱导作用影响微小,但对自身诱导作用明显;当上下旋翼扭矩平衡时,下旋翼总距比上旋翼总距约大1°.     

背景太阳风数值模拟的熵守恒格式

背景太阳风研究是根据行星际扰动传播情况预测空间天气状况的基础,磁流体（MHD）模拟是背景太阳风研究的重要手段.采用一种新的数值计算方式,利用Ideal GLM-MHD将计算过程中产生的磁场散度以c_h的速度向计算区域外传播,从而消去磁场散度;重构部分使用受约束的最小二乘法,将磁场散度作为约束条件添加到重构中,进一步对重构后的磁场梯度进行修正;通量计算采用满足热力学第二定律的熵守恒格式,该格式能够确保在计算过程中熵不增,保证数值稳定.研究结果表明,该方法应用于太阳风数值模拟的求解得到了更加稳定的结果.      

用于表面温度测量的导热反问题的误差研究

讨论了用于温度测量的非稳态一维线性导热反问题。通过引入正则参数 ,这种不适定问题得到了稳定的解。在给出非级数形式的分析解后 ,计算误差可以通过数值解与准确解的比较确定。同时讨论了影响计算误差的各种因素。