多孔介质不混溶二相驱动问题后种新的差分解法及其收敛性分析

考虑不相混溶不可压缩二相驱动问题的数值解，提出了一种新的非均匀网格上的差分解法，用箱型差分格式解浓度方程，用块中心差分格式解压力方程，该法在时间和空间部分非均匀的情况下关于时间和空间具有二阶精度的误差估计。由于该法关于压力是五点差分格式，关于浓度是九点差分格式，求解工作量与传统的有限差分方法求解的一样，但比有限元方法的工作量要少得多，且该法计算稳定，易推广到三维问题中去。     

Poisson方程边值问题差分格式解的收敛性

给出poisson方程的边值问题在均匀网格剖分下的五点差分格式,并给出了差分格式解的收敛性及敛速估计.     

多孔介质中二维混溶驱动问题在时空局部加密复合网格上的有限差分格式

对多孔介质中不可压缩流体的混溶驱动问题,建立了其二维问题在时间和空间上进行局部加密的复合网格上的有限差分格式,压力方程采用5点差分格式近似,饱和度方程采用修正迎风格式,且在交界面上采用线性插值,并给出了误差估计．最后给出了数值算例．     

混凝土温湿度耦合方程组的数值解法

混凝土中热湿耦合方程组的数值解法是将混凝土中热湿耦合传导理论应用于工程实际的关键．算法在空间域采用有限元格式，在时间域采用两点差分格式．研究了在求解域内空间网格划分和时间域划分形式及时间域差分格式对数值解收敛性质和振荡性质的影响，选取了合适的数值计算格式，能使计算结果在收敛性、振荡性和结果精度方面都满足工程计算要求。     

NS方程的非结构化网格方法及其差分格式

采用同位非结构化网格上的有限体积法,对ＮａｖｉｅｒＳｔｏｋｅｓ 方程的ＳＩＭＰＬＥ 算法及差分格式进行了研究．提出了适用于非结构化网格的不用求解单元顶点变量值的二阶混合差分格式．该格式的优点在于：１） 减少了计算工作量;２） 避免了普通混合差分格式因使用简单的一阶迎风差分格式所引起的网格界面方向相关性的问题．最后,采用三角形网格,利用提出的方法及差分格式,对方腔内的驱动层流及绕翼型湍流进行了数值计算,计算结果与基准解或实验值的符合程度优于普通混合差分格式     

带有小参数的抛物型方程的多过渡点差分格式

讨论两类时间上带有小参数的抛物型方程,介绍了多过渡点的选取方法,依此法构造不等距差分格式并证明新的差分格式关于摄动参数是一阶一致收敛.多过渡点差分格式的收敛阶与Bakhvalov法相同,高于Shishkin网格法,但计算量比Bakhvalov法小得多,在实际应用中相当有效.     

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

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

多孔介质二相驱动问题压力方程Galerkin方法后处理的特征有限元方 …

对多孔介质中二相驱动问题提出了一种新的数值解法，即用常规有限元方法求解压力方程，经后处理后，再用特征有限元方法解浓度方程。该法不仅避免了用混合元法求解压力方程带来的困难，而且保持了特征有限元方法的优点，得到用标准有限元方法求解压力方程著不有得到的最优误差估计。     

对流扩散反应方程基于坐标变换的高阶紧致差分格式

基于非均匀网格,提出了一种求解一维定常对流扩散反应方程的高精度紧致差分格式。首先采用坐标变换方法将原方程由物理空间的非均匀网格转换为计算空间的均匀网格,然后给出一阶导数和二阶导数在均匀网格上的中心差分逼近式,并结合变换后的方程,得到了定常对流扩散反应方程具有四阶精度的紧致差分格式。最后,通过数值算例验证了该方法的精确性和高分辨率的特点。数值实验结果表明,对于所研究问题,该方法较不进行坐标变换而直接在物理域上建立的非均匀网格上的高阶紧致格式具有更高精度。     

两相渗流驱动问题有限元—特征线修正差分格式及理论分析

研究多孔介质中可压缩可混溶两相渗流驱动问题的计算方法：压力方程用有限元方法求解，饱合度方程用特征线修正有限差分方法求解，构造了全离散数值计算格式，证明了最佳收敛阶     

高次光滑边界单元的构建与性能

基于Lagrange插值多项式,实现了高次边界单元形函数系数的计算机自动生成.在已有闭合单元的基础上,充分利用粒子几何形状的特性构建了高次光滑边界单元.数值算例表明,与传统的二次单元和闭合单元相比,高次光滑边界单元能够较大地提高椭圆和椭球粒子数值模拟的计算精度与效率.     

直角坐标架下的精化不协调元法

在直角坐标架下直接插值，建立了满足收敛要求的不协调元法；其列式比等参不协调元法简单，直接、且不失一般性，由此建立的平面四边形单元精度高且可推出单元刚度阵显式。     

BCI－代数的换位元

在ＢＣＩ一代数中引入了换位元的概念．证明了换位元的若干重要性质，讨论了换位元与极大元、极小元以及半群之间的关系．     

厚板低阶广义协调矩形元

采用常内力状态下的广义协调条件，将剪切变量用结点挠度和转角表示，导出一个具有１２个自由度的厚板、薄板都通用的矩形弯曲单元。此单元的自由度少，精度高，能通过分片检验，不出现剪切闭锁现象，具有位移型单元简便实用的优点。     

钢筋混凝土平面单元研究

对目前已有钢筋混凝土有限元模型进行了研究分析,指出了这些模型存在的问题;同时,建立了较符合结构实际和材料的力学状态的钢筋混凝土有限元模型,导出了这些模型的单元刚度矩阵。这里,仅是有关钢筋混凝土平面单元的研究。     

精化杂交法及精化平面四边形单元

基于新的泛函、合理的变量假设及应变正交化，提出了称之为精化杂交 元的方法。精化杂交法可以使单元的应变能按假定的应变模式分解，由此得 到相应的分解的单元刚度矩阵，而且常常可以推出显式。精化杂交法有效地 提高了杂交应力元或广义杂交元的精度和计算效率。所建立的平面四边形精 化杂交元，可以作为对著名的Ｐｉａｎ单元的改进。算例表明，所建立的四边形 单元较已有的各类平面四边形单元具有更高的精度和计算效率。     

不规则波导本征值的边界元解

应用边界元法求解不规则波导的本征值,给出了混合元计算的Ｈ和Ｇ矩阵的计算公式,并通过算例验证了算法的正确性     

关于生物高分子元素活性中心分布规律的第四统计力学理论标度的研究(Ⅲ)——论各种化学元素成为生物高分子活性中心元素的先决条件

提出了各种化学元素成为有利于生命过程的生物高分子活性中心元素至少需要满足的八个条件：不能与环境中的OH^-和大量Cl^-发生沉淀；不能与各种有机酸和酰酮基团等过分络合；不能同氨基酸单体过分亲合；不能与生物蛋白分子中大量的酰胺基团过分络合；此外，有益于生命的元素离子还不能先同体内的—SH基团作用。本文还根据金属离子和不同基团的络合能力把各种元素分为有益元素群和有毒元素群，为从元素角度分析病因提供了原始科学依据。     

平面理性八节点曲边四边形有限元RCQ8

由于理性有限元用弹性力学的解插值并且在物理域内直接列式,因而等参元相比,具有物理意义明确和精度高的优点,本文采用最高至４次的平面弹性力学的多项式解对单元位移插值,推出了八节点曲边平面理性元ＲＣＱ８,若干算例表明,本单元可以避免剪切自锁和体积自锁,具有很高的精度和效率以及数值稳定性。