Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

一类拟线性椭圆型方程的CHOQUARD-PEKAR问题

讨论了一类拟线性椭圆型方程的CHOQUARD-PEKAR问题在无界区域中的非平凡解的存在性,对于极小问题Iλ和I∞λ,得到了对于每个λ0,存在α∈(0,λ],使得Iα和I∞α可以达到.     

平面Poisson问题的间接变量等价边界积分方程

使用变分法确立平面Laplace外问题的确切形式；在此基础上，为了避免使用无限区域上的Green公式(不成立)，通过构造紧支集上的无穷可微函数，获得满足非齐次控制微分方程的特解的重要结果。使用非解析开拓数学方法，导出满足控制微分方程的位势函数集的等价边界积分表示式。获得Poisson问题的间接变量等价边界积分方程；通过一些例子，详细讨论了传统的外问题和边界积分方程的形式，表明它们不具有普遍适用性。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

降落伞系统附加质量的研究

降落伞在减速系统中已被广泛地运用，但是其运动方程的确立，尤其是方程中的附加质量项，在中国却少有提及，为了弥补此点不足，文章对附加质量进行了较全面的阐述，从概念引入，在方程中的运用，确定方法及变化规律，到对稳定性的影响分析，对于降落伞此类轻质的，柔性的钝体，附加质量受到多种因素的影响，其规律可由试验数据得出，而由根轨迹图分析可知，附加质量对稳定下落时的降落伞系统稳定性有着重要影响。     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法，建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题，旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法，通过添加冗余约束，避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块，对基于运动嵌套网格 的空间流场进行求解，湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块，基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解，验证了网格变形模块和CFD模块的有效性，然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上，计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷，并与飞行测试数据进行了对比。计算结果表明，文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析，提高了旋翼气弹载荷的预测精度。     

基于本征方程的柔性机翼非线性颤振分析

柔性机翼在气动载荷作用下常常会产生较大的变形,颤振特性会随之发生变化,针对此问题线性理论常常难以进行合理的预测。以几何精确本征梁模型建立了机翼的运动方程,耦合ONERA-EDlin非线性气动模型,建立了柔性机翼的非线性气动弹性分析模型。利用Newton-Raphson和Backward-Differentiation-Formula(BDF)分别求解机翼的静态变形和动态响应,基于机翼平衡位置附近的线性化方程来判断系统的稳定性,进而确定颤振临界速度。通过算例验证了模型的准确性,并分析了不同刚度、后掠角、机翼安装角等参数对颤振速度的影响。     

四旋翼飞行器姿态控制建模与仿真

四旋翼飞行器是典型的欠驱动、非线性、强耦合系统。其姿态控制精度和抗干扰问题一直是研究热点之一。为了实现小型低成本四旋翼飞行器姿态准确控制,详细分析了四旋翼飞行器受力情况,利用牛顿-欧拉方程建立了四旋翼飞行器非线性动力学模型,针对四旋翼飞行器在实际飞行过程中经常会遇到阵风、气流等不确定外界干扰,设计了基于小扰动的PID控制器,并通过对俯仰通道、横滚通道、偏航通道MATLAB/Simulink仿真模型进行仿真测试和结果分析。结果表明:所设计控制算法能够满足四旋翼飞行器姿态控制要求,并具有较好的抗干扰性能。     

液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位，用于制定合理的发射流程以及紧急处置方法，采用计算流体力学(CFD)技术，对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型，为了提高Lee模型在不同压力情况下对相变过程的模拟精度，采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性，随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象，液相内会形成大的漩涡，该漩涡会使得冷热流体不断进行热交换，并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h，发生在停放4 h左右；而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低，最小值接近1 m 3/h。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。