非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

均匀各向同性湍流的脱体涡数值模拟

采用脱体涡模拟方法对均匀各向同性湍流进行了数值模拟,并与Comte-Bellot实验测量结果进行了对比,验证了该文的脱体涡模拟方法对均匀各向同性湍流模拟的可靠性.初始速度场的生成采用Rogallo所提出的构造方法,初始湍能谱满足Von Karman波谱分布;在对流输运项的选择方面,分别采用二阶中心型格式、四阶偏斜对称型中心格式和迎风型低耗散通量分裂格式,考察它们在均匀各向同性湍流模拟中的计算精度和适用性.同时,通过改变计算域大小以及脱体涡模拟方法中的模型常数,达到能谱截断波数的改变,考察它对各向同性湍流计算的能谱以及能谱截断处小尺度涡能量积累问题的影响.      

改进的二维三阶半离散中心迎风格式

对二维三阶半离散中心迎风格式中的权函数给出了简化改进. 在保持格式精度的基础上, 改进后的权函数在二维情况下具有更加简单直接的结构而且严格非负. 该改进方法得到的格式仍然具有半离散中心迎风格式的优点, 同时保持了重构函数的非振性. 时间离散采用保持强稳定性的三阶Runge-Kutta方法, 并利用四阶Lax-Wendroff (L-W) 格式计算磁流体算例中的磁场散度. 用该修正格式计算了二维磁流体数值算例, 得到高精度无振荡的结果, 验证了此方法的有效性.      

刚性方程的一个算法

在文献[3]中构造了一个与二阶Gear方法等价的差分格式,并采用简单迭代求解。本文详细论述了差分格式在实际计算中应采取的若干措施及程序框图,最后给出了数值例子。     

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM+ (AdvectionUpstreamSplittingMethod)格式 ,推广应用于二维非结构网格中 .利用格林 高斯公式对控制体内的变量进行线性重构 ,获得空间二阶精度 ,并采用Barth型限制器以抑制数值解的振荡 .为提高计算效率 ,采用了自适应多重网格法 .最后给出了NACA0 0 1 2翼型和平板激波反射的几个Euler方程的算例 ,并进行了分析比较 .     

高超音速化学非平衡流的隐式格式迭代解法

针对高超音速化学非平衡流计算中的隐式TVD格式给出了一种迭代解法.该方法是在原有的Harten-Yee格式基础上,在每一时间步对无粘部分进行迭代运算,目的是为了加速收敛.借助对模型方程的收敛性分析,证明该迭代算法是收敛的.二维数值实验表明该方法对于高马赫数时的化学非平衡流能明显提高隐式格式的计算效率,其计算结果是满意的.     

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.     

广义拉格朗日乘子MHD方程组的两类逆风分裂格式

开发高性能的磁流体力学数值模拟方法是提高空间天气数值预报研究的一个重要方面.有限体积法的逆风分裂格式具有良好的间断捕获能力，Steger-Warming和AUSM（Advection Upstream Splitting Method）是逆风分裂格式FVS（Flux Vector Splitting）方法中具有代表性的两种格式.采用这两种格式求解具有伽利略不变性的扩展型广义拉格朗日乘子磁流体力学（EGLM-MHD）方程组，对Orszag-Tang涡流问题和三维爆炸波问题进行数值模拟，结果表明两种格式均能得到稳定精确的数值结果.与Steger-Warming格式相比，AUSM格式产生的磁场散度误差更小，计算速度更快.      

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。      

永磁电动机网格自动剖分

有限元法是一种可以用来求复杂工程问题近似解的数值方法,而求解域的网格自动剖分技术有明显的实用价值。本文根据网格自动剖分的一般原则,阐述了稀土永磁直流电动机矩形槽三角形单元的自动剖分方法及相应的几个功能子程序,并利用Auto CAD绘图软件包方便的实现了网格数据检查和图形显示。     

高可靠性航空电子设备热分析中的有限体积法

阐述了热对系统可靠性的影响,对比了数值法求解温度场分布的优缺点.针对航空电子设备热分析中同时要求满足高的分析精度和快的计算速度的问题,提出用有限体积法求解温度场分布的计算方法.根据典型航空电子设备的结构特点,建立了温度场的数学模型,给出了边界条件.在此基础上,用有限体积法对待求解方程离散化,给出了离散化过程中网格划分及边界的处理方法.结合实例进行了温度场求解,计算结果与热测量结果进行了对比,得到了较高的分析精度,验证了方法的有效性.最后对计算误差进行了分析,并提出了进一步热设计的合理化建议.      

低散射目标谐振区散射特性分析

用有限面元模拟导体表面，简化了任意三维金属物体的散射计算问题，用矩量法对标准球和仿球体的单站散射截面进行了计算，分析了雷达目标关于谐振区的散射特性，阐述了低频区将成为“隐身”目标的频域易控制窗口的观点。     

Space debris cumulative flux considering the Interval Distance-based method

The availability of engineering models to estimate the risk from space debris is essential for space missions. According to current research, cumulative flux calculation is mostly carried out based on the equal-width interval discretization. The method discretizes the volume around the Earth into cells defined in earth centered inertial coordinates. The resulting debris flux onto a target object is shown to depend on the chosen size of the cells. To avoid a discretization error, this must be accounted for. In order to present reliable flux predictions for space mission, the algorithm improvement is an ongoing topic for the related research field. The aim of this study was to examine the discretization error during the cumulative flux determination process. Both the effect of interval step length and the orbital boundary are under investigation. Several typical orbits are selected as examples here and the 2018/01/03 TLE (Two Line Element) data published by the US Space Surveillance Network is used as the debris background in this paper. Furthermore, the Interval Distance-Based method for Discretization (IDD) is adopted in this paper. A position-centered flux determination method is introduced based on the IDD method. According to the example analysis, the IDD used in the flux calculation process provides results which are less affected by the interval step-size setup; and the orbital boundary has no effect on the calculation process. In other words, the discretization error is significantly reduced. The position-centered method provided a possible suggestion for the improvement of space debris environment models.     

有限体积法中面积分离散格式的精度分析

讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差,同时与中心格式和顶点格式进行了比较.格式精度分析采用泰勒级数展开的方法.分析结果表明:顶点中心格式和顶点格式类似,总能达到至少是一阶的精度,并且可以通过加密网格的办法减少误差,而且顶点中心格式更适合于分区计算.通过圆弧亚音速绕流和NACA0012翼型绕流流场的数值计算,验证了分析的结论.     

2D约束Delaunay剖分生成表面模型的表面网格

提出一种生成表面模型的表面非结构化网格的方法.该方法仅假定表面模型的每个面片在取定适当的投影平面时可看成是一个单值函数.对表面模型的每个面片,首先在其相应的投影面上进行二维约束Delaunay剖分,然后对网格顶点进行插值.由于采用联动剖分的策略使得单独生成的各个表面网格在共同边界处匹配.与通常的生成表面网格的网格前沿法相比,本方法无需预先对域的边界进行离散,边界的离散体现在约束Delaunay剖分中恢复限定线段的边界细分过程中,减少了用户干预.通过合理地指定2D约束Delaunay三角化时网格单元尺寸分布函数,可以有效控制表面网格对表面模型的逼近精度和实现自适应的表面网格.实验结果证明了本算法的有效性.     

高精度法矢下切割面自适应的凸曲面射线寻迹

针对电大尺寸目标难以精确解析表达带来的一致性几何绕射算法应用难问题,提出了基于三角网格、适用于任意凸曲面的射线寻迹（TM-tracing）算法。应用工程中较易获取的三角网格及其协议,设计了一种满足快速多边搜索条件的网状数据存储链表;提出了满足寻迹要求的高精度法矢求解算法;采用切割面自适应调整的弧形拟合寻迹方法实现了爬行波寻迹算法;结合一致性几何绕射理论（UTD）实现了暗区场值求解算法。任意网格曲面射线寻迹结果表明:本文提出的寻迹算法适用于包括球、柱和锥在内的任意光滑凸曲面,寻迹偏差小于1.61%,寻迹速度为2.8 s,具有一定的工程应用价值。      

高频随机响应分析中的一种有限元网格细化标准

根据波动理论,满足随机响应工程计算精度要求的有限元网格大小与待求高频随机响应问题的频率范围上限有关。提出了一种新的有限元网格细化标准———三分之一波动细化标准。以简支板高频加速度响应分析为例,通过NASTRAN分析结果与解析解的对比验证了其有效性。