三维最小二乘无网格方法精度与误差分析

为了得到最小二乘无网格方法在求解三维Euler方程时的计算精度及其应用于三维流场求解的准确性,通过数学方法对该无网格方法在正交等间距离散点分布情况下的计算精度进行了论证,给出了基于正交等间距和非正交均匀布点的4个三维验证算例.并将该文计算结果与其他网格方法得到的计算结果和理论结果进行了比较,两种方法得到的计算结果吻合较好,与理论结果的误差在合理范围以内.结果表明将最小二乘无网格方法应用于三维流场的求解是完全可行的.     

非结构四面体网格上三维非线性扩散方程的有限体积法

为了有效地数值模拟科学和工程中有广泛应用的非线性扩散方程,在三维线性扩散方程非结构四面体网格的有限体积法的基础上,提出了一个计算非结构四面体网格上非线性扩散方程的有限体积法。方法采用网格单元中心作为计算节点,相对于网格点的方法,计算量减少了一半。用L agrange因子法得到网格点上的值,考虑了网格中心点和网格点的相对位置,更适应大变形的网格。利用算子分裂,使计算更加简单。用N ew ton-B iCG STAB法来求解得到非线性方程组。数值结果表明:该方法具有二阶精度、保持通量守恒、对大变形的网格适应性强。     

非均匀裂缝性油藏蒸汽驱的自适应网格计算

为提高非均匀裂缝性油藏流动数值计算的效率,自适应网格算法对裂缝和基质采用了不同网格结构,避免了基质渗透率非均匀情况下基质饱和度空间不连续给网格的粗化、细化带来的严重数学问题,并采用重正化方法计算裂缝粗网格等效渗透率。二维蒸汽驱算例的计算结果表明,相比于全精细网格计算,自适应网格算法不仅能够很好地提高计算效率,其计算精度也得到了较好的保证。     

网格计算系统中服务部分的研究与实现

网格计算是近年来研究的一个热点,它旨在使动态多机构虚拟组织中的资源得到共享和协同问题得到解决。在介绍网格及网格计算的概念,分析网格计算系统的基础上,提出了网格软件系统中服务部分的实现构想,并给出了应用实现实例。     

自动生成四边形网格的方法及其在数值模拟中的应用

采用间接方法生成四边形网格,首先利用改进的两点前沿推进法把计算区域剖分成三角形网格,然后采用插点和细分的技术生成单元全部是四边形的网格,通过边互换、删点和局部插点技术进一步光滑平顺,得到适用于数值计算的网格。剖分结果表明,该方法能够在任意二维平面区域内自动生成全四边形网格,并能生成光滑过渡的局部加密网格和贴体性较好的边界层网格。该方法具有算法简单,计算量少的特点。利用所生成的网格对计算传热学中的典型算例-方腔自然对流进行求解,计算结果与基准解吻合,网格质量能够满足数值分析计算的要求。     

一种三角形表面网格质量综合度量方法

在计算机图形学中,有限元等科学计算与分析方法得到了越来越广泛的应用．而有限元单元网格质量的好坏将直接影响有限元计算过程的收敛性以及有限元解的误差．通过对计算机图形学中常用的三个表面网格模型进行分析,利用模糊数学的方法,给出了一种三角形网格质量综合度量方法．利用该方法得到的评价指标可以综合、快速地评价一个有限元网格模型的质量．     

基于大涡模拟方法的闭口箱梁气动三分力系数识别

以苏通长江大桥主桥闭口箱梁断面为例,采用大涡模拟方法的动态Smagorinsky模式,计算了雷诺数与风洞试验相当的断面气动三分力系数,并与实验值及雷诺时均方法的结果进行了对比.同时也进行了二维计算域的计算,并对不同三维计算域尺度及网格分辨率进行了计算对比,计算得到了不同计算域及网格对计算结果的影响情况.     

基于子网格材质平均的高阶FDTD(2,4)方法

文章提出了一种基于子网格材质平均的高阶FDTD(2,4)方法,通过对原网格内外环路所围区域进行二次剖分,分别得到内外环路的平均电磁参数;再通过加权平均得到等效电参数,避免了积分运算,简化了建模过程.FDTD(2,4)运算仍在原网格上进行,所以没有增加计算量.该方法用于求解三维介质球散射,数值结果表明,在相同计算量的条件下,该方法具有更高的精度.     

多重网格在黏性流动最小二乘等几何模拟中的应用

针对最小二乘等几何方法模拟黏性流动时条件数大、迭代法收敛速度慢的问题,提出了基于多重网格技术的加速方法。计算中自动生成一系列疏密不同的网格,在最密网格上用最小二乘等几何方法将Navier-Stokes方程离散为代数方程组,用多重网格方法作为独立求解器或共轭梯度法的预处理器迭代求解所得到的代数方程组。对雷诺数为100、400、1 000和2 500的顶盖驱动流进行了数值模拟,计算中进行23次迭代可使方程组的余量降低10个数量级,流动特征量的计算误差在1%以内。计算结果表明,通过多重网格技术加速迭代,提高了最小二乘等几何方法模拟黏性流动的计算效率。     

非均质复合材料力学性能的确定性多尺度计算方法

针对具有多尺度特性的非均质复合材料进行结构强度分析,提出了一种确定性多尺度计算方法——有限细胞法(FCM)。该方法通过矩阵凝聚和插值近似得到粗网格的刚度矩阵,在计算得到粗网格解后,通过回代技术可获得全域细网格解。为验证FCM的计算精度和计算效率,构造了多个数值算例,并与理论和有限元(FEM)结果进行了对比,结果发现:相对于FEM,FCM可显著扩大计算规模,提高计算速度,且易于并行;相对于其他多尺度计算方法,FCM在构造粗网格刚度矩阵过程中无需引入微观边界条件,计算精度高,易于降尺度计算。因此,FCM是分析复合材料结构强度的一种有效的多尺度算法,同时该方法易于推广至多尺度的非线性结构分析、热分析、地下水渗流分析等领域。     

埕东油田铬冻胶堵剂的改进

本文分析了埕东油田铬冻胶堵剂的优缺点,研究了铬东胶成冻时间和冻胶强度的影响因素,提出了铬冻胶堵剂的改进配方,并讨论了铬冻胶由半液法堵剂改为双液法堵剂的有关问题。     

配合物的合成及其吸氧规律性研究

本文合成了一系列配合物,分别进行了吸氧,脱氧性能研究,探讨了配合物吸氧规律性     

DNA sequence and expression of the B95-8 Epstein-Barr virus genome

The complete (172,282 base pairs) nucleotide sequence of the B95-8 strain of Epstein-Barr virus has been established using the dideoxynucleotide/M13 sequencing procedure. Many RNA polymerase II promoters have been mapped and the mRNAs from these promoters have been assigned to the latent or early/late productive virus cycles. Likely protein-coding regions have been identified and three of these have been shown to encode a ribonucleotide reductase, a DNA polymerase and two surface glycoproteins.     

在氢键系统中激发的质子孤立子的量子力学特征

根据我们所提出的在氢键系统中的新哈密顿函数，并且使用完整的量子力学方法，本文得到了该系统中激发的质子孤立子的动力学方程组，依据这组运动方程得到了质子孤立子的量子力学特征，并得出这类质子孤立子的运动服从经典运动规律，其它有趣的性质已同时得到。     

均布载荷下变宽度悬臂薄板的大挠度问题

应用拖带坐标（Ｃｏ—ｍｏｖｉｎｇＣｏｏｒｄｉｎａｔｅ）描述法［１］建立了一类变宽度悬臂薄板大挠度问题的简化计算方法，编制了相应的通用电子计算机程序，并对数值算例进行了实验验证，理论计算与实验结果相符。     

电除尘器应用于CFB锅炉仿真设计的研究

在分析循环流化床锅炉尾部烟尘特性及其对电除尘器性能影响的基础上,以其配套的电除尘器为研究对象,建立了一系列仿真数学模型.并以原始工况参数为依据,对电除尘器进行仿真设计,提出一系列改进措施.     

基于犹豫度下区间直觉模糊信息的决策方法及其应用

对区间直觉模糊集的度量方法进行了研究．首先介绍了区间直觉模糊集的相关概念,在充分利用犹豫度信息的情况下,构造了一种新的得分函数,并证明了其相关性质;同时,将其与几类精确函数作了比较,并结合加权算术平均算子和加权几何平均算子,提出了一种基于区间直觉模糊信息的决策方法;最后,通过实例分析验证了决策方法的有效性和合理性．     

钉螺在动水中运动特性研究

通过试验，研究了钉螺在动水中爬行运动和钉螺在动水中的起动流速，推导了钉螺在开厣状态下起动流速公式；还根据试验结果描述了不同流速区段钉螺运动特性，分析了钉螺在渠系中扩散方式，从而为防止钉螺扩散提供了理论依据。     

单振子辐射器在圆柱状分层生物组织中的场分布

以牙根分层生物组织为例，采用近似模型，导出了单振子辐射器在圆柱状分层生物介质中的近场公式。通过数值计算，作出了电场分布的曲线，进一步分析比较了各层组织中场强的相对大小，从理论上探讨了微波治疗牙根尖周病的原理，理论分析结果与临床实践相符。     

二阶递归数字滤波器的极限环分析

The problems of limit cycles in second-order recursive digital filters with a rounding off quantizer have been studied more systematically in this paper. The properties of limit cycles have been summarized,the distribution regions of various limit cycles on the coefficient planes(β1,β2)have been pointed out, and the linear stability conditions of second-order recursive digital filters have been derived, All conclusions in this paper have been verified by computer simulations.