状态空间精细法计算弹性地基板的动力响应

采用状态空间精细法计算弹性地基板的动力响应量。对空间域采用有限元计算，对时间域采用状态空间精细法。在计算过程中，不解线性方程组，直接进行精细迭代计算，程序简单，效率高。给出２个数值算例，对数值结果与直接积分法进行了比较，其计算结果比较接近。     

二维抛物型方程精细积分法与差分法比较

可用单内点子域精细积分法，求解二维抛物型方程初值问题。当单内点精细积分中的传递函数即指数函数用Taylor展开式的一阶近似来替代时，精细积分转化为差分方程。研究这一对应关系，使各种常见差分格式均找到对应的单点精细积分格式，并在单点精细积分的一般公式中获得统一表达式。     

解波动方程的精细积分法及其数值稳定性分析

将精细积分法用于求解波动方程。详细论述了精细积分法的数值方法，并给出了相应的计算公式。数值算例表明，用精细积分法得到的解与精确解十分吻合，比有限差分法具有更高的精度。同时，推导了解波动方程精细积分法的稳定性条件。与有限差分法相比，精细积分法有更好的数值稳定性。精细积分法的计算公式适用于求解实际工程问题的波动方程，并易于推广应用到二维和三维波动方程的数值求解。     

时滞抛物型方程的拟小波精细积分法

对时滞抛物型方程初值问题提出了采用拟小波精细积分法进行计算,采取拟Shannon尺度函数为权函数,利用小波配点法对空间域离散,将时滞抛物型方程转化为常微分方程组,然后用高效的精细积分法求解时滞常微分方程组。这种方法的优点是精确度高、稳定性好。数值算例表明,本文提出的拟小波精细积分法具有很高的精度,因而是一种有效的数值方法。     

分层半空间中Love表面波的精细方法

精细积分法是求解线性常微分方程两端边值问题和初值问题的精细算法.应用精细积分法(PIM)和扩展Wittrick-Williams(W-W)算法求解了横观各向同性、分层半空间中的Love表面波问题.岩层是由分层介质置于半无限空间上组成.Love表面波对应于波数-频率域线性常微分方程的本征值问题.利用本征值计数技术,扩展W-W算法可以不遗漏地找到所有本征值,得到计算机精度意义下的精确解.     

解波动方程的精细积分法及其数值稳定性分析

将精细积分法用于求解波动方程。详细论述了精细积分法的数值方法,并给出了相应的计算公式。数值算例表明,用精细积分法得到的解与精确解十分吻合,比有限差分法具有更高的精度。同时,推导了解波动方程精细积分法的稳定性条件。与有限差分法相比,精细积分法有更好的数值稳定性。精细积分法的计算公式适用于求解实际工程问题的波动方程,并易于推广应用到二维和三维波动方程的数值求解。     

关于结构动力学中精细积分法的研究

介绍了精细积分法的基本原理及精细积分法的稳定性,通过给出的一个算例,比较了用精细积分法与中心差分法、Newmark法的计算结果。结果表明用精细积分法得到的位移响应最接近理论值,而且当取相同位数有效数字时,与理论值完全一致。     

对流扩散方程的精细积分并行算法

讨论了一维扩散方程的全域精细积分和子域精细积分的并行算法，给出了对流扩散方法的子域精细积分并行算法。子域精细积分考虑了细积分法高精度的特点，又避免了全域积分的大矩阵运算；春精度优于单点精织积分法，同时子域精细积分很容易实行并行计算。算例表明了精细积分并行算法有良好的并行加速比和效率。     

波动方程的无网格-精细积分方法

将无网格法和精细积分用于波动方程的计算.在空间上用无网格法进行离散,用修正变分原理处理本征边界条件;在时间域上用精细积分法求解动力学方程,然后给出两个波动方程的算例.数值结果表明此方法是稳定、精确的.     

拱坝随机振动分析中的精细直接积分法

为了解决传统方法计算拱坝随机振动效率较低的问题 ,将虚拟激励法和精细直接积分方法结合 ,针对虚拟激励法中荷载为振荡函数的特点提出了振荡函数的精细直接积分法并应用于拱坝地基动力相互作用的有限元边界元无限边界元 (FE- BE- IBE)耦合模型中。以某拱坝为例对Wilson- θ法、Cotes积分法和振荡函数的精细直接积分法进行了比较计算。结果表明 ,精细直接积分方法可以在保证计算精度的前提下增大计算时步的步长 ,提高计算效率 30～4 0倍 ;精细直接积分法中振荡函数的精细直接积分法效率更优于 Cotes积分方法     

对溶质运移模型的并行化处理

在了解和广泛收集前人研究资料的基础上,将并行算法引入到求解首采区卤水动态二维模型中关于溶质运移的问题中.把溶质运移方程按时间分裂法分成5个子方程,针对这5个子方程讨论了全域精细积分和子域精细积分的并行算法,给出了对流和扩散方程的子域精细积分并行算法.子域精细积分考虑了精细积分法高精度的特点,又避免了全域积分的大矩阵运算,其精度优于单点精细积分法.     

用改进的一点配置法计算多孔催化剂有效因子*(I)方法的建立

用改进的一点配置法求解多孔催化剂颗粒扩散－反应边值问题，从而获得催化剂颗粒有效因子的解。改进的一点配置法将整个一维区间划分为几个子区间，然后在每个子区间内用一点配置法求解非线性扩散－反应边值问题。与一级不可逆反应精确解比较，用本文方法计算结果的精度十分令人满意。     

扩展有限元中非连续区域的一种积分方案

为促进和推广扩展有限元在实际复杂问题中的应用,提出一种适用于扩展有限元非连续区域积分的简易积分方案。利用基于该积分方案的扩展有限元程序,结合粘聚裂纹模型实现了对三点弯梁开裂过程的模拟。结果表明,由此积分方案得出的结果与常规分片积分法方案得到的结果以及已有研究结果十分吻合。与已有积分方案相比,此简易积分方案既能保证精度和收敛性,同时无需复杂的分片积分过程,也无需考虑零能模式的控制措施,大大地降低了建模的复杂性和成本。     

化学驱软件中化学平衡的并行计算

为了探索化学驱油藏数值模拟软件高效的整体并行化方案 ,对 DQCHEM2 .0软件中化学平衡计算部分进行并行化改造 ,设计了两种区域分解并行化的方式。一种方式是将整个区域按 CPU数进行分解 ,每个 CPU计算一个子域 ;另一种方式是按照额定的子域中网格单元的数量进行分解 ,每个 CPU将计算多个子域。比较而言 ,前者数据传输量少 ,而后者有利于负载平衡。在并行程序中 ,采用 MPI消息传递库实现数据的传输。测试结果显示了较好的局部并行效率。 8个 CPU的加速比达到 6.4。对测试数据的分析表明 ,两种区域分解方法适用于不同的情况     

基于非重叠区域分解算法的大地电磁法二维正演模拟

区域分解算法采用分而治之的思想,将大规模问题转化为若干个小问题进行求解,已成为大规模数值计算领域的常用算法之一。将非重叠区域分解算法引入到大地电磁法二维正演模拟中。首先将整个求解区域分解为多个互不重叠的子域,子域之间共享边界元素;然后对每个子域采用有限差分进行离散,采用Schur补偿算法解耦得到共享边界节点上的未知数,并作为子域问题的边界条件得到关于子域内部节点的线性方程组;最后,利用直接求解算法对上述方程组进行求解,实现了大地电磁法二维正演。该算法的准确性和可行性通过多个地电模型的对比试算得到了验证。此外,还统计分析了采用不同子域分区方式和分解个数时的计算耗时,结果表明子域分区的方式对计算效率影响不大,但子域分解个数的影响则较大,进行区域分解时需要选择合适的子域个数。     

Helmholtz方程的径向基配点不重叠Schwarz交替法

将径向基函数配点法和不重叠型Schwarz交替法结合用于求解Helmholtz方程.该方法把求解大规模问题转化为求解多个小的子区域问题,克服了在求解大规模问题时用一般的全域径向基配点法所带来的配置矩阵为非对称满阵,且高度病态的问题.首先给出具体算法,然后给出算法的收敛性,最后通过数值算例得出相应结论.     

RTI在基于MDA仿真中的应用

为满足HLA(High Level Architecture)自身发展的需要,将模型驱动架构MDA(Model Driven Architecture)和高层体系结构HLA相结合,提出了基于MDA的HLA仿真框架。分析了这种框架下对于主题域的划分,为进一步研究RTI（Run-Time Infrastructure）的应用问题奠定基础。与一般的RTI划分不同,这种框架下,RTI不是作为模型转换的平台,而是作为服务域的子域,这样更有利于模型的重用。针对这种划分,讨论了RTI与应用域模型的集成应用问题。该设计是MDA在HLA中应用的一种新思路,它使仿真业务模型与HLA/RTI的耦合性降低,可提高其重用能力。