流动模拟二阶精度的简易保持研究

首先分析了各类低阶和高阶差分格式,指出二阶精度格式具有综合优势．接着以二维对流扩散方程为例,推导出延迟修正的CDS格式,并定义为DCDS格式．通过标量绕静点输运、自然对流及混合对流等实例,并与UDS,CDS及QUICK格式比较,论证了DCDS格式不仅能维持计算结果的二阶精度,而且极大地约束了CDS格式所固有的振荡和越界行为,边界处理也很简单．总之,在程序的效率、调试、边界处理、计算精度及稳定性等因素权衡下．DCDS格式是保持二阶数值精度的简易途径。     

Cosine软件包组件参数程序Coslatc燃耗功能的验证

Cosine软件包组件程序Coslatc中的燃耗计算是反应堆组件参数计算的重要功能之一,Coslatc的计算精度很大程度上影响了堆芯物理的计算结果。通过对单栅元到组件的模型进行计算,将Coslatc的计算结果与国际基准题以及基于ORIGEN-S的SCALE程序中的TRITON模块计算结果进行对比进行数值验算。通过计算结果表明,Cosine软件包组件程序Coslatc具有单栅元到组件再推广到堆芯的计算能力,其燃耗计算精度基本和SCALE程序中的TRITON模块计算结果相当。     

谱元方法求解对流扩散方程及其稳定性分析

探讨了一维对流扩散方程的一种高精度数值解法,该解法在空间上采用了Chebyshev谱元方法,在时间上结合了半隐式Adams方法。通过数值算例验证了解法的可行性,利用特征分析法得到了对流扩散方程谱元求解时不同离散形式的稳定性条件,并对数值求解的稳定性进行了预测。通过时间步长、网格划分、对流项和黏性项插值阶数的影响研究表明:耦合Chebyshev谱元方法和半隐式Adams方法在求解对流扩散方程时能够获得高精度的数值解;时间离散时Adams方法的黏性项采用一阶插值形式、对流项采用二阶插值形式,在未增加计算量的同时能够获得较大的稳定区域和较高的计算精度。     

一种同位网格上的迎风Galerkin 有限元方法

通过 Green定理将对流项变量从微分算子中分离出来 ,从插值函数入手引入迎风格式 ,是对强对流问题 Galerkin有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速—压同位网格公式 ,构成了一种对高 Reynolds数流体流动数值模拟比较有效的有限元方法。数值试验表明 ,采用该方法能较好地提高计算精度。     

一种求解反应对流扩散方程的稳定化谱元方法

针对谱元方法求解二维非稳态反应对流扩散方程中出现的稳定性问题,提出了一种稳定的高精度数值方法。该方法在空间上将Chebyshev谱元方法和一致逼近迎风方法相结合,时间上采用分步θ-格式。通过解析解算例验证了该方法的精度及数值稳定性,并对含有不同类型边界层的反应对流扩散问题进行了求解。研究表明:一致逼近迎风项的增加扩大了谱元方法求解反应对流扩散方程的稳定域,在对流项及反应项占优时保持了数值解的高精度;对于含有边界层的复杂反应对流扩散问题,数值解在整个计算区域内获得了一致收敛的结果。研究工作对谱元方法在反应对流扩散问题高精度数值求解中的应用提供参考。     

下落薄液膜中孤立子波强化传热的机理

发展了分析下落薄液膜中孤立子波热输运现象的边界元方法，散度理论被用于压力补偿函数边界元公式的非线性对流容积积分，结果使速度和温度梯度消失，全部公式用速度和温度来写，回避了对流离散的困难，大大节约了计算内程和计算量，并改进了计算精度。用法向应力平衡条件和打靶迭代法计算孤立子波形状。文中给出孤立子波内和波纹状膜流中的流场、温度场、沿着壁面和液体交界面处的Ｎｕｓｓｅｌｔ数变化。在孤立子波回流区中存在大的     

高等级公路施工放样数据现场快速计算技术

叙述高级编程计算器按模块化程序的编制的技巧与步骤.并以线元分析法原理计算线元任意点坐标为计算模型,快速计算路线放样数据的方法.     

三维层流自然对流换热边界层方程的数值解法

采用坐标变换及凯勒单元法求解强制对流换热边界层方程已比较成熟,但对自然对流换热的求解还未见报道。为了检验此方法用于求解自然对流换热的计算精度,采用该方法数值求解了竖窄条三维层流自然对流换热边界层方程,计算结果与前人的实验结果相符。结果表明,用此方法求解三维层流自然对流换热边界层方程是可行的。     

高阶精度格式在火灾模拟中的应用

对流作用在火灾诱导的烟气运动中占主导地位。为了更好地研究火灾诱导的烟气运动,该文构造了一种新的4阶精度的有限差分格式,用以离散火灾烟气运动的控制方程——低M ach数可压缩N av ier-S tokes方程的对流项。将此格式加入到美国国家标准技术局开发的火灾模拟软件F ireD ynam ics S im u lator(FDS)中,对典型的油池火灾和单室火灾的问题进行了数值模拟。将速度和温度分布的计算结果与实验结果和二阶精度的FDS程序计算结果相比较,结果表明该文提出的格式在较少的网格下得到了比二阶精度格式更接近于实验值的结果。这说明在火灾模拟中,相对于二阶精度的格式来说,高精度格式是一种更好的选择。     

沥青路面温度和荷载耦合分析

采用罚函数法处理本质边界条件,编制了无网格法程序。针对工程实例,建立了无网格法计算模型并按顺序耦合方式进行了模拟分析。计算结果、现场实测数据和有限元分析结果吻合较好,计算精度高于有限元方法。结果表明:无网格法应用于沥青路面结构温度场和应力场的计算分析是可行的。     

热连轧层流冷却的数学模型

本文以武钢１７００热连轧厂层流冷却现场实测数据为基础，通过传热学的理论分析，确定了对流换热系数是解决问题的关键，并对７种模型结构形式进行了比较，研制了精度较高的对流换热系数数学模型，在此基础上编写了模拟武钢１７００连轧层流冷却过程的程序，取得了较高的模拟精度。     

面向对象的三维对流扩散问题有限元程序设计

根据三维对流扩散问题的有限元分析,采用面向对象的程序设计语言C ,建立向量类、矩阵类、单元类、结点类、材料类和有限元方法类等,实现三维对流扩散问题的有限元分析程序开发,对相关的对流扩散问题进行数值模拟.结果表明与用结构化语言开发的传统有限元程序相比,程序更易重用、维护、扩充,并且可以融入用C 编制的大型通用有限元科学和工程计算软件.     

无相间物质传递化学驱浓度方程算子分裂隐式解法

为了改进化学驱数学模型 U TCHEM显式求解浓度方程计算速度慢、计算结果精度低的缺点 ,研究了隐式求解组份浓度方程的方法。根据具有无相间物质传递关系的化学驱油藏流体渗流过程满足的相行为 ,推导出了化学驱数学模型 UTCHEM物质守恒方程的等价形式 :饱和度方程和组份浓度方程。利用算子分裂技术将组份浓度方程分裂为扩散方程和对流方程 ,隐式交替求解对流方程和扩散方程得到组份浓度方程的隐式解。扩散方程采用隐式局部一维格式差分离散 ,利用追赶法求解 ;对流方程选用了隐式迎风格式差分 ,并且结合油藏模拟问题的流场是有势场的特点 ,实现了对流方程隐式差分显式求解。所建立的隐式求解浓度方程的方法提高了计算精度 ,可以加大计算时间步长 ,加快计算速度。     

一类四阶发展方程的混合有限体积格式

针对一类带对流项的四阶发展方程,运用特征线方法,直接从积分守恒形式出发导出了混合有限体积格式,给出了格式的截断误差,证明了该格式按离散L^2模对时间和空间具有一阶精度．数值算例表明,该格式计算效果良好．     

光电测距导线测量的精度探讨

光电测距导线测量的精度是当前测绘界普遍关心的课题之一。本文利用PC——1500机编制的固定程序,按坐标变量平差法对单一导线进行严密平差,计算时根据不同的测角与测距精度,拟定了六种典型方案,分别计算了待定点位置的误差元素,如纵、横坐标误差,误误差椭圆元素以及方向差和边长误差等等。然后根据计算数据对平差后导线点的点位误差分布,导线边方向误差的特点进行探讨。同时还可按不同方案的计算结果,分析导线测量所能达到的精度,从而为导线的设计或制定相应规范提供参考。     

按垂直温差计算折光系数的方法之二

本文根据Chaughey和Palmer(1979)的研究成果“对流边界层(即混合层)热通量H近似地随高度而线性减小”,提出了一种新的按垂直温差计算折光系数的方法。这种方法应用于山区或丘陵地区的三角高程测量中,能够提高折光系数和三角高程测量的精度。     

三维无界元和有限元的耦合方法

本文提出了三维无界元和有限元的耦合方法,即用无界元来模拟远场效应,用有限元模拟近场效应.同时,给出了三维四点无界模型及形函数、位移函数,无界元的单刚矩阵等整套计算公式,编制了三维无界元和等参有限元的耦合程序.该程序具有提高计算精度和节省机时、节省存贮量等优点,是解决岩土工程问题的有效工具.     

构造对流扩散方程高精度非均匀网格格式的指数变换方法

提出了一种基于非均匀差分网格,构造求解对流扩散方程的高精度格式的指数变换方法.引入指数函数,将对流扩散方程变换为扩散反应方程,消除了数值求解中较难处理的对流项.采用优化差分方法推导出扩散反应方程基于非均匀网格的高精度差分格式,进而通过逆变换得到对流扩散方程的高精度格式.理论分析表明,该方法具有3至4阶精度,当计算区域为均匀网格时取得4阶精度.数值实例表明,在相同的非均匀网格系统中,此方法的计算精度明显优于传统的隐式差分方法.在水环境的实际模拟计算中,根据物理量的变化规律灵活地调整非均匀网格的间距,不仅能增强高精度差分方法的实用性,而且可以取得比均匀网格方法更为精确的计算结果.     

基于Additive Runge-Kutta方法的激波聚焦起爆高精度数值模拟

基于详细氢氧化学动力学模型,建立了描述氢氧爆轰的多组分反应欧拉方程组. 针对建立的反应欧拉方程组,数值方法上采用3阶Additive Runge-Kutta方法对时间项进行积分,采用5阶精度的加权本质无振荡(WENO)格式对空间对流项进行离散,自主研发了大规模高精度计算程序. 该程序能够处理化学反应源项引起的刚性问题,且能节省计算时间和计算内存. 对半球型、半椭球型、圆锥型3种结构形式凹面腔内的激波聚焦起爆过程进行了数值模拟,数值模拟研究得到了不同结构形式凹面腔内的激波聚焦起爆过程.      

分区粘弹性介质样条边界元反分析数值解法

通过改造非均质样条边界元列式方法,建立了分区粘弹性介质边界元反分析的数值求解格式,该数值解法与均匀弹性体边界元法相类似,不需作分区组合计算,简化了计算程序,减少了计算工作量、提高了计算精度和解题效率,为解决非均质围岩反分析中遇到的困难奠定了基础。