低功率氩电弧喷射推力器中的辐射

为了研究低功率氩电弧喷射推力器中的辐射损失,采用辐射模型和非辐射模型对其工作过程进行了对比数值模拟分析。采用二阶精度无波动、无自由参数的耗散差分格式(NND格式)数值求解耦合电磁源项和辐射源项的N-S方程组,并采用隐式残差光滑法加速收敛;采用有限控制容积积分方法离散求解椭圆型偏微分电磁场方程,并采用逐点超松弛迭代方法加快收敛速度。数值模拟结果对比给出了辐射模型和非辐射模型流动分布情况,并比较了两种模型推力器推力、比冲和推进效率。研究结果表明,与非辐射模型相比,辐射模型比冲低0.137%,推进效率低1.03%,辐射对低功率电弧喷射推力器性能影响比较小。      

电弧喷射器流场模拟中电场计算的多重网格法

在建立电弧喷射器内等离子体流场电场计算模型以及基本算法基础上,研究了采用多重网格对收敛速度的改进效果。数值计算证实了此方法的有效性尤其对流场电场耦合计算是有利的。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.     

高超声速压缩拐角湍流计算的收敛性分析

应用基于shear stress transport(SST)湍流模型的隐式紧耦合方法对高超声速压缩拐角流动进行了数值计算.考察了差分格式、限制器以及湍流时间步长等7个因素对湍流计算过程稳定性和收敛性的影响.结果表明:差分格式和限制器对收敛性有很大影响,总体而言耗散性越大的差分格式和限制器收敛性越好;增加湍流时间步长对收敛性影响很小;适当降低湍流变量最低限制值有利于湍流方程的收敛;增大Courant-Friedrichs-Lewy(CFL)数和隐式内迭代次数可显著加速收敛,但均存在最优值,超过最优值则加速效果不明显.      

扩压器进气道内流及激波／边界层干扰的数值研究

应用隐式矢通量分理解差分格式结合改进的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，对扩压器跨声流动进行了数值计算，准确地捕获了激波，预测了由于强激波作用引起的湍流边界层分离现象，与实验数据比较，结果令人满意。数值实验表明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，计算结果可用于高速飞行器气动特性的研究 。     

电弧喷射推力器化学非平衡数值模拟

为了准确把握电弧喷射推力器工作过程物理机理与特征,采用化学非平衡模型对不同压缩室直径下低功率氮电弧喷射推力器工作情况进行了数值模拟。模型考虑了工作过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分,4反应的有限速率化学反应模型。采用二阶精度NND格式数值求解耦合电磁源项和化学反应源项的N S方程组,采用有限控制容积积分方法求解由麦克斯韦方程组推导出的电磁场方程。数值模拟的结果揭示了推力器内部电弧能量输入作用和高温电离气体的离解电离状况,给出了不同压缩室直径下推力器的推力、比冲和推进效率。结果分析表明,压缩室直径对推力器性能具有较大影响。     

微波等离子推力器喷管流动数值模拟

为了准确把握微波等离子推力器喷管流动的机理与特点,采用冻结流和非平衡流两种模型对其进行了对比数值模拟分析。非平衡流模型考虑了流动过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分、4反应的有限速率化学反应模型,采用二阶精度NND格式数值求解耦合化学反应源项的N S方程组。数值模拟的结果揭示了喷管内的流场结构,反映了喷管内的离解电离状况,得到了推力器的推力和比冲。分析表明,数值方法有效,计算结果合理,具有工程应用价值,能为相关研究提供参考。     

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

可压流的龙格-库塔时间步进求解法及其应用

研究了用经典四阶龙格－库塔法计算可压缩流场的可行性,并将其用到了电弧喷射推力器内部等离子体流场的数值求解中。应用情况表明,在结合了局域时间步长、隐式残值光滑加速收敛措施后,本格式能够成功地计算比较复杂的可压缩气体流场以及等离子体流场。有关计算结果揭示了气体流经电弧喷射推力器通道但无电流时形成的纯气动流场以及有电流通过时形成的等离子体流场的丰富的结构和一些重要的影响因素,为研究其过程机制提供了依据。     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

Numerical analysis of deflection control of a gas plasma jet based on magnetohydrodynamic staggered electrode configuration

To control the deflection of the gas plasma jet, a new analytical method is proposed based on the Magnetohydrodynamic(MHD) technique. Based on the typical MHD power generation model, the applied voltage is applied to the staggered electrodes, that is, a pair of electrodes on the same side wall are connected to generate an axial current in the channel. Under the action of the magnetic field perpendicular to the direction of the flow, the plasma is subjected to electromagnetic forces perpendicular...     

电弧喷射推力器电磁场模拟

为了揭示电弧喷射推力器内部电磁场特征及其与高温电离气体间的相互作用机理,建立了比较完善的电磁场模型,详细讨论了电磁场数值模拟方法及其计算效果.模型同时考虑了电场力、洛仑兹力和电子压强梯度对传导电流的贡献,由电子动力学方程和粒子微观碰撞理论推导出电流方程,由麦克斯韦方程组推导出电磁场控制方程,采用包括强隐过程(SIP)迭代法在内的十种迭代方法求解离散的电磁场方程.研究表明,SIP迭代法是适用于电弧喷射推力器电磁场模拟的快速有效的方法,模型能够准确地反映推力器内部电磁场与高温电离气体间的相互作用机理.     

利用定子槽口宽度削弱永磁电机电磁噪声研究

着重研究了永磁电机定子槽口宽度对电磁噪声的影响。通过电磁力解析计算得出电磁力与定子槽口宽度的计算关系。利用Ansoft软件对1台电动汽车用永磁电机进行电磁场计算,求解改变定子槽口宽度对电磁力波的影响,并用ANSYS软件对该永磁电机进行不同定子槽口宽度的2D声场计算。进行了试验验证,仿真计算与试验结果基本一致,可为该类电机削弱噪声设计计算提供参考。     

电磁场作用下包晶合金定向凝固过程传输模型的建立

建立了二元包晶合金在电磁场作用下定向凝固过程中多场耦合的数学模型,并根据包晶转变不同阶段的特点,提出了包晶转变中固相分数的计算办法。该模型可以计算电磁场作用下包晶合金定向凝固中的流动、传热、传质的过程,为电磁场作用下包晶合金定向凝固过程的数值计算及对实验结果的分析提供理论依据。     

辉光放电等离子体对边界层流动控制的机理研究

本文求解电位势方程得到电场分布,假定电场强度大于击穿阈值的区域为等离子体区,给定电荷密度,得到了作用于流体上的电场力。通过求解带源项的NS方程,研究了一个大气压下的均匀辉光放电等离子体对边界层流动的影响,考察了电场力做功对流动的影响。本文研究结果同文献[6]一致,即电场力总体上使边界层流动加速。另外,电场力做功对流动参数的影响可以忽略。文献[6]给出的线化电场模型和本文得到的电场相比,具有较大的差别,该差别引起了流动参数显著差别。     

平板附面层等离子体流动控制的数值模拟

通过求解电场中的拉普拉斯方程和赫姆霍兹方程的变形形式,成功地将等离子体激励对平板附面层流动的影响,以体积力向量的形式引入到NS方程之中.借助求解的电荷密度,通过将计算结果与实验数据的对比,可确定德拜长度、最大电荷密度、形状因子等可调参数的取值原则,建立起平板附面层等离子体流动控制的数值模拟方法,为将等离子体流动控制方法应用于外流及内流场中的强剪切流动控制,奠定关键性的技术支撑.