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

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

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

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

机翼非平行边界层稳定性研究

从Navier-Stocks(N-S)方程导出曲线坐标系下的抛物化稳定性方程(Parabolicstabilityequation,PSE),研究机翼非平行的可压缩边界层稳定性问题。发展了求解PSE的高效数值方法:引进法向变换,使得在临界层与壁面之间的扰动量变化最快的区域有更多法向网格点;采用包含边界邻域在内的完全四阶精度的法向差分格式,这对方程精确离散至关重要;以及全局法和局部法相结合的数值方法及其新的迭代公式,能大大加速收敛并得到更精确的特征值。算例分析研究了扰动增长因子和形状函数等演化曲线。     

采用流固耦合方法数值分析减压器动态响应

采用流固耦合数值模拟技术研究了减压器工作过程的动态特性。流体控制方程为三维基于动态网格描述的Eu ler方程,采用非结构动网格有限体积格式求解;固体结构的运动方程由拉格朗日方程导出,采用纽马克方法求解。提出了可变边条零厚度虚拟挡板技术模拟了活门由闭到开的瞬态过程。计算得到了上游不同开启压力以及不同开启持续时间减压器弹性元件的动态响应特性,结果表明,虚拟挡板技术模拟物体间由零距离开始分离的实际过程是成功的,且对最小网格尺度没有限制;将减压器上游高压的开启持续时间延长能够避免阀芯开度的振动。     

探针式光纤探头的有效检测深度

研究了探针式光纤探头的有效检测深度。采用扩散方程描述光在强散射生物组织中的传播规律,运用纽曼边界条件,使用灵敏度矩阵对有效检测深度进行了研究。灵敏度矩阵和光源与检测器的空间位置分布密切相关。光源与检测器对的灵敏度矩阵定义为：光源位置和检测器位置分别作为光源的光场分布的乘积。设计了6组针对不同参数的模型,参数主要包括光源、检测器光纤芯径、光源、检测器光纤中心距、光纤探头所处的深度、生物组织的吸收系数和散射系数。对探针式光纤探头在各组模型中的有效探测深度进行了仿真,重点研究了六种参数对检测深度的影响。所得结论可以推广到分析光纤检测器的空间分辨率和多光纤检测器的研究中。     

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

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

边界条件对三维空腔流动振荡的影响

基于全隐式无分裂方法,数值求解三维非定常雷诺平均的N av ier-S tokes方程,紊流模型采用修正的B-L代数模型。运用此方法,本文数值模拟了三维开式空腔的超声速、跨声速、亚声速流动。研究了空腔口处剪切层的非定常运动现象,着重探讨了入口边界条件对三维空腔流自持性振荡的影响。数值计算结果表明,与层流速度型相比,入口条件为紊流速度型和特征边界型得到的结果能够吻合。     

广义系统静态输出反馈镇定的一个充要条件

研究广义系统静态输出反馈镇定的问题．基于广义李雅普诺夫方程,给出广义系统通过输出反馈镇定的一个充要条件,该结果是正常系统输出反馈镇定相应结果的推广。另外给出了求取镇定控制器参数的一个算法,并通过数字算例说明该算法是有效的。     

基于带波动算子非线性Schr(o)dinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schr(o)dinger方程进行了数值分析,提出了一个含参数的二阶守恒差分格式,根据参数选取的差异,该格式既可隐式计算也可显式计算.对初值条件进行了中心差分离散,使其具有二阶精度,从而与守恒格式的精度一致.利用矩阵理论证明了差分解的存在惟一性,并利用一个重要的不等式在先验估计的基础上,运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解.数值实验表明该格式具有较高的计算效率.     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。