实现边界元并行分布计算的一种新方法

 给出了实现边界元并行分布计算的一种新方法。它将边界元离散节点分为与计算网络所拥有的处理机台数相对应的若干组, 以此实现任务分配, 独立地计算形成每组节点对应的离散边界元积分方程, 最后利用迭代方法行求解分布的边界元方程组。以弹性问题为例, 在Inmos T800 Transputer 系统上实现了边界元计算。验证了本文方法的可行性与有效性。     

双材料界面裂纹应力强度因子的无网格分析

应用改进的无网格方法对双材料界面裂纹进行了分析,积分区域采用基于节点的Voronoi图,通过在裂纹尖端加密节点提高位移场的计算精度,采用位移外插法确定出应力强度因子,避免了应力场的奇异性和振荡性。本文对由环氧树脂和铝 环氧树脂组成的双材料进行分析,对比了均匀离散、均匀细化和局部细化三种不同模型,通过与光弹性实验结果比较可知,采用裂纹尖端局部细化的模型计算量低并且计算效果好,其计算结果与实验结果基本吻合。     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之二:简化模型及应用

提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier-Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架,对压气机不稳定边界的判断通过计算和检验特征值的方式来实现,给出了有限差分法和基于Gauss-Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性.     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之一:理论模型及算法分析

本文提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier－Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架。对压气机不稳定边界的判断通过计算和检验特征值的方式来实现。本文给出了有限差分法和基于Gauss－Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性。     

角度-空间间断谱元法的遮蔽条件下二维热辐射特性研究

航空发动机燃烧室温度分布对燃烧室结构设计、燃烧过程检测均具有重要意义。燃烧室结构封闭且复杂、存在大面积的辐射遮挡,这给计算其内部由热辐射引起的温度变化分析带来了一定困难。本文提出采用角度-空间间断谱元法求解遮蔽条件下封闭腔内热辐射问题。在求解过程中,通过黎曼迎风方案对空间单元边界上的辐射信息进行分割,实现单元间信息的有效传递。通过与文献结果进行对比发现：角度上采用间断谱元法可有效抑制辐射遮蔽带来的数值振荡,相同角度离散单元数下,可以获得比离散坐标法更光滑的计算结果。通过分析网格数和基函数阶数对数值结果的影响,发现该方法具有hp收敛特性。此外,控制角度计算节点相同,在与离散坐标法相当的计算精度下,消耗更少的计算时间。     

敏感度法估算转子危险转速区

提供一种估算含多层叠片阻尼器转子危险转速区的“敏感度”分析法.采用传递系数法,计算刚度区间各离散点对应的临界转速,通过分析临界转速相对支承刚度变化的敏感程度,确定转子危险转速区.首先,划分转子计算节点,计算各节点物性参数值;其次,取定支承刚度存在的范围,按二分法离散刚度区间,利用传递系数法,循环计算每个离散刚度点对应的临界转速和振型,直到找出振型突变的两个相邻刚度点;接着,做出临界转速随刚度变化的规律曲线,根据曲线找出临界转速变化剧烈时对应的子刚度区间,该区间与上步找出的刚度点行成的区间应是相符的;然后,若希望得到更高精度的结果,可以在子刚度区间内,重复步骤二,直到得到满意结果.该方法易于实现模块化操作,避免在敏感支承刚度区出现危险转速.     

多局部节点异步抗差航迹关联算法

为解决航迹异步与系统误差并存情况下的多局部节点航迹关联问题,提出一种基于区间序列离散度的多局部节点异步抗差航迹关联算法。定义区间型数据集的离散信息度量,给出系统误差下航迹序列区间化方法,通过累次积分计算离散度,结合多维分配进行关联判定。针对多局部节点上报目标不完全一致现象,设置零号航迹管理关联质量。与传统算法相比,无需时域配准,可在系统误差下对异步航迹直接关联。仿真结果表明,算法能在局部节点上报目标不完全一致场景下实现有效关联,且正确关联率随局部节点数目的增加或目标密集程度的增大而提高。     

基于SPH的无网格点快速生成算法及翼型绕流研究

提出了一套基于光滑粒子流体动力学(SPH)方法的高效无网格点生成算法。首先，采用Catmull-Rom三次样条对边界曲线进行参数化建模，由笛卡尔点生成背景节点系统并通过改进加速射线法判断背景节点位置。然后，对复杂二维固体边界和二维圆柱入水两个算例分别进行布点测试，结果表明所提算法能够大幅度减少射线法的判断次数，极大提升了离散布点的计算效率。最后，针对Re=10³时的NACA0012翼型绕流问题，通过所提布点算法对其离散布点后，进一步采用SPH方法进行了数值模拟和气动特性分析。仿真结果表明，计算结果与文献结果吻合较好，进一步验证了布点算法的合理性和有效性。     

瞬态热应力问题的组合杂交有限元方法

求解瞬态热应力问题的杂交元方法面临Ladyshenskaya-Babuska-Brezzi(LBB)条件的检查,其稳定性难以保证,从而提出瞬态热应力问题的组合杂交有限元方法。建立瞬态热应力问题的基于区域分解的组合杂交变分原理及其有限元离散;通过数值算例验证求解瞬态热应力问题的组合杂交元的数值性能。结果表明:在大宽厚比网格下,八节点六面体组合杂交元(CHH(0-1))可以获得与二十节点六面体元(B20)精度相当的位移计算结果和更好的应力计算结果。     

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

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

Fast approximate maximum likelihood period estimation from incomplete timing data

 To estimate the period of a periodic point process from noisy and incomplete observations, the classical periodogram algorithm is modified. The original periodogram algorithm yields an estimate by performing grid search of the peak of a spectrum, which is equivalent to the periodogram of the periodic point process, thus its performance is found to be sensitive to the chosen grid spacing. This paper derives a novel grid spacing formula, after finding a lower bound of the width of the spectral mainlobe. By employing this formula, the proposed new estimator can determine an appropriate grid spacing adaptively, and is able to yield approximate maximum likelihood estimate (MLE) with a computational complexity of O(n²). Experimental results prove that the proposed estimator can achieve better trade-off between statistical accuracy and complexity, as compared to existing methods. Simulations also show that the derived grid spacing formula is also applicable to other estimators that operate similarly by grid search.     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

一种适用于飞行器复杂外形的表面与空间网格生成方法

本文通过建立反映机身横截面外形特征的控制点概念，及在控制点间分配网格点数方式，给出了一种非常实用、交互式的表面网格生成方法，所生成的网格失真度小，物体保形好，生成适用于CFD计算需要的表面网络花费的时间少，适用性强。基于无限插值理论，引入本文提出的网格正交控制、物面法向量控制、加权平均光顺措施，给出了一种改进的空间网格生成方法，有效克服了传统无限插值网格生成方法在复杂外形网格生成方面的缺陷，有效改善了网格的法向疏密性、贴体性及周向网格的均匀性，对真实飞行器部件网格生成非常迅速。以某飞行器复杂机身表面与空间网格生成，及无人机翼身组合体空间网格生成为例，检验了本文网格生成方法的实用性和有效性。     

具有元外节点的相容元素

使用有限元素法做应力分析时,经常需要在不同的部位采用疏密程度不同的网格。在应用所谓“选择加密网格”进行最佳网格分析时尤其如此。这时,在粗细网格区域的分界处出现了位移不协调。本文建议一种利用元外节点的内插方案,从而得到完全相容的元素。文中叙述了这类元素形状函数的构造方法,适用性和优缺点。     

二维扩压器和进气道内正交贴体曲线坐标的快速生成法

贴体曲线坐标系的形成是用有限差分法或有限元素法在诸如扩压器和进气道这类任意形状的域内数值求解各类场时,所必需要做的第一步工作。 通常,对这种贴体曲线坐标系有若干要求。首先,所形成的曲线坐标系必须是“贴体”的,以便于精确地施加边界条件。其次,曲线坐标系最好是正交的,以使得计算平面内的控制方程变得较为简单（不出现交叉导数项）,同时也使得便于应用双层零方程紊流模型。第三,所生成的网格系统的网格距要能随意控制,以便在域内计算变量变化     

数值网格生成技术中方程非线性系数对网格正交性和疏密分布的影响

对椭圆型方程中非线性系数对网格正交性和空间分布的影响进行了理论分析和数值验证。结果表明 :方程中混合偏导数项系数具有控制网格正交性的作用;其它两个二次偏导数项系数具有调节网格空间分布的作用,且其调节的强弱和作用与其比值有直接关系。利用方程系数的这些特性,提出一种 Numman边界条件的给定方法,较好地实现了网格正交和疏密分布的有效控制。     

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。     

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

结冰计算中翼面网格重构研究

对结冰数值计算中的网格生成进行了深入研究,针对结冰后几何边界的复杂性及边界依赖于时间步长变化对网格生成的影响展开了分析.网格生成是结冰计算的基础,文中基于求解椭圆形偏微分方程(PDEs)生成了未结冰时翼面周围的背景网格,采用Hilgenstock源项构造方法,合理控制了结冰边界周围的网格间距和正交性.模型表面结冰后,为...     

两种椭圆型方程求源项方法在喷管内流场网格生成中的应用

用椭圆型方程生成二维网格，发展两种Hilgenstock求源项方法，实现对网格与边界间距及与边界正交的控制。