求解非稳态对流占优问题的非标准无网格Galerkin方法

用θ加权法离散时间域,并将四种稳定化方案与无网格Galerkin方法相耦合进行空间域的离散。在无网格Galerkin方法中,采用线性基和具有连续的权函数,基于移动最小二乘法构造了高阶导数连续的形函数,从而避免了有限元方法中采用线性元插值时,因忽略稳定项中二阶导数项而降低计算精度的问题。数值计算表明:本文构造的方法成功地消除了非定常对流扩散方程中对流项占优时的数值伪振荡现象,并具有计算精度高、稳定性好、算法实施简单、前后处理方便的优点。特别是所构造的MFLS方法非常适宜于求解非定常的对流扩散方程。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

采用混合场方法计算二维介质-导体复合散射体雷达散射截面

在二维介质导体复合散射体计算中,对导体部分应用积分方程时,不论电场积分方程还是磁场积分方程,去求解导体的散射问题时,在内谐振条件下将出现解答的不唯一性或不稳定性,产生谐振伪解,给计算工作带来困难.为避免谐振伪解,对导体部分采用混合场积分方程方法结合介质部分体积分方程和矩量法,对TM入射平面波情形建立了该模型求解方程,并编程求解出了该模型雷达散射截面(RCS).     

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明：如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。     

一种无数值积分的间断Galerkin有限元方法

在使用间断Galerkin有限元方法的计算过程中,需要构造相应的积分表达式作为数值求解的出发点,继而会引入体积分和面积分。对于这些积分项的值,一般需要通过数值积分的方法获得。当需要使用高阶间断Galerkin有限元方法时,数值积分计算精度的要求会相应地增加,其所需的计算量将变得很大,而数值积分的计算量又在很大程度上决定了间断Galerkin有限元方法的计算效率。针对这一问题,通过建立Lagrange插值多项式基函数和Jacobi正交多项式基函数的一定关系,构造了一种无数值积分的间断Galerkin有限元方法显式半离散格式,并对不同条件下的线性和非线性一维、二维守恒律进行了直接数值模拟,得到了理想的数值结果。该方法不再需要通过数值积分来计算每个单元的积分项,而且有效地达到了间断Galerkin有限元方法的高阶精度,其对于构造更为高效的高阶间断Galerkin有限元计算方法具有非常显著的意义。     

可渗透面对流FW-H方程伪声产生机制和抑制方法

可渗透面对流FW-H方程及其积分解在远场气动噪声预测中广泛应用，但当涡波通过可渗透面时会产生伪声传播。本文旨在清晰阐明这种伪声的产生机制，并提出有效抑制方法。耦合分析可渗透面对流FW-H方程的右端厚度源和载荷源，发现厚度源的物质导数和载荷源的散度操作可使积分面源自动过滤伪声源。伪声的产生源于对流FWH方程的求解采用了Farassat提出的分部积分公式，因部分积分项被忽略，导致伪声源自动过滤功能失效。从厚度源和载荷源中抽取含有涡波扰动的项，在对流波动方程求解过程中保留相应的物质导数和散度操作，以抑制涡波伪声传播。数值测试算例验证了伪声产生机制和抑制方法的正确性。     

对流-扩散型方程的一种简单、有效的欧拉-拉格朗日分裂格式

本文提出了一种简单而有效的数值求解一维对流-扩散型方程的欧拉-拉格朗H(E.-L.)分裂格式。无论对于线性对流-扩散方程和非线性Burgers方程,或者是对流占优或非对流占优,它都能给出较好的结果。文中给出了四个典型而又较困难的计算例子,与精确解的比较表明,数值结果是相当令人满意的。     

求解Burgers方程的特征中心型有限体积法

以Burgers方程为试验模型,提出一种新的有限体积法,空间上的离散采用中心型加权基本无振荡重构,时间上的离散采用数值积分。其中积分节点上的数值通量由特征理论回溯求解,保持了物理量沿特征方向传输的特性,计算量相对Runge-Kutta法明显减少。数值结果表明了方法的有效性和稳定性。     

背腔式缝隙天线阵列的分析与设计

采用等效磁流分析方法,将波导激励的单背腔缝隙单元阵场求解问题离散,得到缝隙表面的边界积分方程。应用矩量法求解该边界积分方程,得到单背腔缝隙单元阵的场结构。再应用阵列天线理论仿真计算了一种易于共形的缝隙天线阵的方向图,与实际测量结果进行了对比,取得了较好的效果。     

粒子群混合算法在变导向器面积计算中的应用

求解非线性方程组经典方法具有严格的局部收敛性;粒子群等进化算法解决了全局收敛的问题,但计算效率偏低,存在最优解不稳定的问题。结合经典Newton-Raphson法的超线性收敛速度和粒子群算法全局收敛能力的粒子群混合算法具备2类算法的优点。在迭代初期采用粒子群算法获得的近似全局解作为Newton-Raphson算法的初始值,以确定高精度的解。利用粒子群混合算法在发动机变导向器面积的大偏离计算中获得了较好的收敛效果,解决了常规Newton-Raphson法不收敛的问题。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究，将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元，采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元，更好地表述了真实物面特征；物面边界相邻体网格单元相应构造为高阶体网格单元，其余体网格单元采用八节点六面体单元，以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解，验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外，由于采用了隐式计算方法，仅需较少的时间步就能迭代收敛。     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

混合初边值问题解的渐近性质及其差分解

 本文讨论了当t→∞时拟线性方程组初边值解的一些渐近性质。在假设非定常问题的解唯一、定常问题的解存在时,对同一边界条件,证明了当t→∞时非定常方程组的解趋于定常问题的解,讨论了用差分方法求解这种问题时,边界点上一些量的精确算法;为了消除在物体尖角附近出现的局部非物理解,得到正确解,我们还给出了一个简单的耗散型物面边界格式;最后给出了一些尖角凹形体的算例。     

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

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

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

高超声速飞行器飞行动力学特性不确定分析

 高超声速飞行器分析模型存在较大不确定性,其给出的动力学特性与真实值之间存在偏差,因此研究飞行动力学特性分析结果的可信度水平对控制系统设计具有重要意义。针对典型的乘波体构型高超声速飞行器,在建立气动/结构/推进相互耦合的动力学模型基础上,利用非概率区间来描述模型中的不确定参数,并将复特征根不确定范围求解问题转化为频率和阻尼比两个实数的不确定区间分别求解,给出了动力学模态特征根、频率及阻尼比的不确定边界。分别采用直接蒙特卡罗(DMC)模拟方法、基于泰勒展开的区间分析方法(TIAM)和基于多项式逼近的区间分析方法(CIAM)对高超声速飞行器飞行动力学不确定性进行了研究。结果表明:CIAM计算时间适中,且给出的边界更为准确、安全,适合在控制系统设计和验证过程中使用。     

基于泰勒展开法的转子系统动力特性区间分析方法

对转子系统动力特性运用一种非概率-区间分析方法进行分析.基于区间数学和1阶泰勒展开理论的区间分析方法将非确定参数支承刚度和连接结构刚度视为区间向量,运用泰勒展开法建立了转子系统固有频率的公式.区间分析方法降低了传统的概率分析方法对不确定参数信息的过分要求,为解决含有非确定参数的转子系统动力特性问题提供了一个途径.应用区间泰勒展开法和概率方法对数值算例进行了分析,并比较了结果.当参数非确定性小于20%时,计算得到的转子系统固有频率区间上下界与真实值区间上下界误差小于2.2%.建立了转子系统动力特性试验装置,试验结果验证了方法的有效性.      

模拟MKDV方程的格子Boltzmann方法

目前,格子Boltzmann方法已被广泛应用于模拟各种非线性物理方程。文中用D1Q2模型给出MKDV方程的带修正项的BGK型格子Boltzmann方法。通过对演化方程的泰勒展开并应用多尺度技术恢复了宏观方程。数值模拟表明文中所建立的模型是有效的。     

跨声速伽辽金有限元法所用上风技术的改进

伽辽金有限元法是解跨声速绕流的有效方法,对超临界流动,为了自动捕捉到激波,通常采用了上风技术,但导致了气流通过激波时质量守恒条件得不到满足。本文通过Bateman变分原理和加罚泛函找到了强制满足质量守恒约束条件的附加项,大大改善了用伽辽金有限元法解超临界绕流的计算结果。