SIP-AMR-CESE模式对CR2055背景太阳风的模拟

使用三维太阳行星际自适应网格守恒元解元太阳风模型（SIP-AMR-CESE MHD）,模拟从太阳表面到地球轨道附近的太阳风.该模型使用六片网格技术,同时利用PARAMESH软件包实现网格自适应.在该模型的基础上,通过增加广义拉格朗日乘子（GLM）磁场散度误差消去方法,完善网格加密放粗判据,微调加速加热形式等方法,使模拟结果与观测更好地符合.另外,通过控制不同时刻的计算区域,显著提高了模型的计算效率.在此基础上,给出了模型改进后模拟得到的CR2055太阳风稳态解与观测的对比分析.      

特定消谐逆变器开关角度的复形调优解法

针对特定消谐式逆变器开关角度难于求解的问题,提出了一种基于约束条件下n维极值复形调优法的开关角度求解方法,分析了目标函数和约束条件的设计方法.该算法对开关角度初值的选取只需考虑约束条件,大大增加了开关角度初值的取值范围;求解过程中不存在解发散或奇异无解问题,与迭代法解非线性超越方程组相比可靠性更高;通过约束条件的限制,使求解过程自动消除了可能出现的窄脉冲;通过目标函数和加权因子可灵活的进行开关角度优化.算例和实验结果验证了方法的正确性和有效性.      

强粘性层流理论及在粘性流计算中的应用

对高雷诺数流动计算，为了解决部分粘性项甚至全部惯性粘性项落入误差与流动物理尺度的关系问题，本文提出强粘性流动理论。强粘性流中至少有一个粘性项与惯性项同量阶，理论包含了物理尺度各向极限、经典边界层和多层边界层理论为其特例，给出了从经典边界层向物理尺度各向相同极限演化的尺度规律和粘性惯性诸项变化的量阶关系，阐明了粘性与惯性力强相互作用将在剪切层的法向以流向同时“激发”涉尺度结构。对粘性流计算，利用强粘性剪切流尺度律重新标度NS格式的修正微分方程，给出临界网格尺度与流动物理尺度和差分格式精度的关系，得到部分粘性项落入误差和计算结果为非物理数值粘性解的二个判据。并以流场中的边界层、驻点和分离点领域计算为例说明理论的应用，对强粘性剪切流计算、证实部分粘性项甚至全部粘性惯性项落入误差的问题值得重视。     

航空发动机状态监控亚定方程组有效数字算法

缺少状态监控的信息量是基于数学模型的航空发动机状态监控的主要困难。航空发动机状态监控亚定方程组缺少成熟可靠的求解算法。提出一种求解航空发动机状态监控亚定方程组的有效数字算法。此算法按实际工程应用所需的有效数字位数将解空间离散化,再加上亚定方程组本身的物理背景约束,然后用求解整数规划问题的方法将亚定方程组求解出来。仿真结果表明此算法是有效的,且2或3位有效数字的仿真结果曲线拟合较好,能较好地满足实际工程应用的需要。     

二维可压缩Euler方程组的SUPG有限元数值解

本文在文献[1]的基础上,构造了二维非定常可压缩Euler方程组的SUPG变分方程组。文中对由文献[1]提出的预估——校正算法作了分析,并用该算法对变分方程组进行解算,所做算例分别为:无粘性激波在固壁上的反射、零攻角10％抛物翼型亚、跨声速绕流、带攻角NACA 0012翼型跨声速绕流,取得了较好的数值结果。     

两冕流间CME事件数值模拟的改进

通过偶极子场和六极子场适当叠加，改进猜解磁场，使猜解磁场在太阳南北极符号相反，然后采用理想磁流体力学方程组(MHD)，由猜解磁场与太阳风流动相互作用计算出稳态自洽解，得到定性上与观测比较接近的具有两个冕流的背景结构．在两个冕流间采用具有同心圆磁场位形的触发模型触发CME事件，研究CME的日冕传播特征．模拟结果表明，CME被约束在两冕流间传播，CME闭磁场位形和磁云横截面磁场位形相似，可以解释1AU处观测磁云的部分特征；在CME附近，存在压力和Lorentz力起主要作用的区域，这可以为分析1AU处CME事件的观测数据提供帮助．     

使用GAO-YONG湍流方程组 对顶盖驱动方腔流的计算

运用SIMPLER方法、QUICK格式及交错网格技术求解GAO-YONG湍流方程组,对二维顶盖驱动方腔流进行数值模拟,得到方腔内的涡量分布和方腔中心水平剖线上的速度变化曲线,并与DNS(Directly Numerical Simulation)结果进行定性对比,有较好的一致性.结果表明,GAO-YONG湍流方程组中与湍流多尺度现象相对应的机械能方程对此湍流方程组能够合理地模拟出顶盖驱动方腔流中复杂的湍流粘性分布以及湍流能量逆转现象起到了关键作用.研究进一步证明了GAO-YONG湍流方程组能够得到复杂流动的真实粘性场.     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

三维磁电热弹材料裂纹问题的超奇异积分方法

应用位移基本解,磁电热弹材料三维裂纹问题被转化为求解一组以裂纹表面位移间断为未知函数的超奇异积分方程问题.进而,采用主部分析法得到了裂纹尖端的应力奇异指数及奇性应力场的解析表达式.然后,通过将裂纹表面位移间断未知函数表达为位移间断基本密度函数,使用有限部积分方法对超奇异积分方程组建立了数值方法.最后,通过典型算例计算,得到广义应力强度因子的变化规律.