用任意拉格朗日--欧拉方法求解流体力学方程组

从积分形式的二维Lagrange流体力学方程组出发，用有限体积格式进行计算，针对不规则的四边形网格，在用次级网格的方法确定网眼内物理量的梯度的基础上，提出一种体积加权的方法来构造网眼内的线性插值多项式。将新网格进行了细划，先确定新网格小网眼中心点的密度，将小网眼的质量直接计算出来，再计算新网格的质量，然后确定相应的密度，从而实现高精度重映。最后用ALE方法进行数值模拟，得到了预期的效果。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

计算参数对数值修正载荷识别算法的影响

针对在载荷识别计算中经常遇到的累积误差问题,提出了一种在每个时间步长内迭代修正的 载荷识别算法。首先利用拟静态算法得到载荷初值,再使用数值迭代算法对其进行修正计算,仿真结果表明,该修正算法可以有效地减小由于累积误差导致的发散,得到收敛的识别结果。针对上述算法,本文以多输入多输出简支梁为模型,分别分析了区间放大系数、区间分割系数和精确度指标3个计算参数对于识别结果的影响。计算结果显示,参数的选择对算法的效率和精度影响很大,不当的参数甚至可能引起识别结果严重发散,所以选择合适的计算参数对于数值修正算法十分重要。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

用非交错网格求解曲线坐标系下的流动控制方程

本文采用非交错网格法,求解任意曲线坐标系下的流动控制方程。文中阐述了压力振荡的根本原因和相应的消除措施,导出了压力修正方程中控制体边界上曲线速度分量的修正表达式。在这些表达式中,出现了一个表征一个网格间距和两个网格间距压力梯度差分值的附加差值项。当出现压力振荡时,该差值很大,可以有效地消除振荡,而当压力场没有振荡时,这一差值又很小。作为计算方法可靠性的检验,本文分别计算了二维直通道和渐缩形混合管内的气流流动问题,结果是令人满意的。     

DSMC方法模拟微尺度流动时的特征参数选取

利用DSM C算法对微尺度平板附面层流动进行了数值模拟研究,分析了算法中3个主要特征参数:采样循环数、单位网格粒子数和子网格数对模拟结果的影响。分析结果表明,当采样循环数增加到一定程度后,计算结果趋于一致,考虑到计算机时,采样数取8左右为宜;单位网格粒子数过多或过少都会使计算机时增加,太少了还会影响计算精度,建议4~5个左右合适;子网格数对运算时间影响很大,分析认为每个网格内的子网格数为4时最为节约计算时间。最后通过微通道内流动计算进一步验证了本文参数选择方法的可行性。     

三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

二维非结构网格生成及其数值模拟

提出一种二维非结构网格生成算法以及一种新的绕翼型的数值计算格式。首先,在Delaunay生成算法的框架下结合了阵面推进算法的优点,从覆盖整个计算区域的三角形着手,给出判断准则,由这些三角形的边直接生成新点,让新点与符合要求的三角形的边的两点相连构成三角形,这样生成的三角形网格具有较好的几何性质;其次,针对二维N0012翼型跨声速绕流问题提出了一种新的计算格式。通过建立插值模板的概念,在每一个模板上运用最小二乘和MUSCL思想相结合的办法针对初始变量构造线性插值多项式,改进了原有的权系数的计算准则,将各插值多项式通过加权平均得到在整个模板集上通用的插值多项式,从而构造出高阶的有限体积格式;最后,计算结果表明本文工作具有一定实际应用价值。     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

GENERATION AND APPLICATION OF UNSTRUCTURED ADAPTIVE MESHES WITH MOVING BOUNDARIES

INTRODUCTIONMany engineering applications require thenumerical prediction of strongly unsteady flowsinteracting with moving bodies.A case in pointis that the detailed calculation of sphere motionsin an unsteady field induced by blast waves is akey issue to determine the drag force exerting up-on the sphere by the surrounding gas and,inturn,to reveal the mechanism of unsteady effectson the drag coefficient.To cope with such pro-blems,one should simulate unsteady flows withmoving boundary co…     

基于IPO的进气道扩压器电磁散射计算方法

进气道RCS目标特性在整个飞机中占有较大比例,但作为电大尺寸腔体,其RCS特性评估与外机身相比有其固有困难。本文研究了适合于任意腔体RCS评估的迭代物理光学算法,该方法基于磁场积分方程属于高频算法范畴,可以较精确地计算出进气道金属表面电流矢量分布并给出令人满意的RCS结果。通过软件开发实现了任意腔体RCS特性评估,并用标准算例考核了程序准确性。应用该软件对某型心线一定的金属进气道S弯扩压器开展计算研究,获得了不同宽高比扩压器在入射波长为0.03 m下的电磁散射特性变化规律。     

Ultrasonic Computed Tomography Imaging Method of Concrete Materials Based on Simulated Annealing Genetic Algorithm

Stability and accuracy of the imaging results are still unmet practical demands for ultrasonic computed tomography(CT)of concrete material.To address these issues,a CT technique based on simulated annealing genetic algorithm(SAGA)is presented in this work.Firstly,a natural weight matrix with clear physical meaning is introduced in the inverse algorithm and then a quadric broadening objective function is formed according to the propagation characteristics of ultrasound in concrete.After that,the simulated annealing(SA)searching is added to speed up the inverse process and to improve the convergence and stability of the algorithm.Finally,the optimal inverse imaging results have been achieved by variable ectopic adaptive genetic algorithm.The numerical simulation experiments have shown that the usage of the correct priori information and the excellent characteristic of SAGA in searching the global minimum value of the function have produced accurate and effective results with stable numerical values.The imaging resolution is improved and the imagining results reflecting the inner defections of the tested objects are more reliable and accurate.     

面向快速原型制造的海量数据模型重建(英文)

根据已有实物的测量数据进行模型重建 ,在机械产品逆向建模、计算机视觉、基于二维轮廓数据的生物外形重建等领域中具有重要应用价值。随着坐标测量设备的发展 ,获取包含被测物体更多细节的海量数据已非常方便 ,但大量的测量点却给模型重建带来了困难。本文首先提出了精度可控的海量数据自动简化算法。为了提高算法的效率 ,文中提出了一个数据集空间划分策略。根据简化后的数据集 ,应用步进立方体方法重建模型的三角网格曲面表示。由于种种原因 ,重建的三角网格模型常常含有不希望有的孔洞。为此 ,本文给出了一个算法产生形状优化的三角片以修补网格模型中的孔洞。经过孔洞修补 ,完全封闭的三角网格模型可以直接输出为快速原型制造中广泛应用的 STL文件。应用实例说明了本文的方法的可行性     

一种新的非结构动网格的生成方法

弹簧方法是一种常被采用的非结构动网格生成方法,文中采用一种新的非结构动网格生成方法即弹性体方法(Elastic solid method,ESM)。和弹簧方法相比,虽然弹性体方法的计算量大,但是可以处理大位移的情况,尤其善于处理网格扭转很厉害的情况。另外,采用弹性体方法可以减小网格重新生成的频率,从而可以减小映射误差,同时也可以减小网格重新生成带来的计算量。将两种方法的数值模拟结果进行对比,进一步表明该弹性体方法在处理动网格问题中的优势。