径向Hartree-Fock-Slater方程的新的数值解法

本文提出径向Hartree-Fock-Slater(HFS)方程的一种新的数值解法。使用二阶或四阶有限差分法,把HFS方程离散化,构成一组拟线性代数本征值问题。其系数矩阵与解有关,需用迭代方法求近似解。本工作对数值方法做了改进,使得数值处理比较完整和统一。已编成计算机程序,对Cu⁺和Mn²⁺进行了计算。所得结果与文献[1,2]的结果一致,表明此法的稳定性和精确性良好。还对一系列原子和离子问题进行过计算,证明此法能成功地用于中性原子和正离子的计算。     

肺毛细血管与组织间物质交换的双重介质模型

本文在文献[4]的基础上利用双重孔隙介质模型研究了肺毛细血管与组织间的物质交换。利用小参数展开的分析方法和数值计算方法求解了非线性的控制方程。通过比较说明在一般的生理条件下小参数展开法所得的解已能充分描述肺毛细血管和组织间的物质交换状况。最后文中还利用所得解分析了肺水肿形成的机理。     

基于Excel的化学计算方法

对化学计算中遇到的求解方程或方程组的复杂计算问题,利用Excel提供的数值计算和绘制函数图像的功能,通过单变量求解和规划求解能使问题得到圆满解决。与其它软件比较,该方法对初值的选择和计算过程中迭代方法的选择、迭代次数、精度、误差等都可以进行自由设置,更能体现出门槛低、简单、快捷的特点。     

一种计算Feshbach共振态寿命的新方法

在自然碰撞坐标下构建偏分势能面, 利用数值传播方法求解沿反应坐标的核运动方程, 然后用过渡态波函数的相移因子构造反应体系共振态寿命矩阵. 这是一种直接计算化学反应散射共振寿命的量子散射方法. 用此方法计算了I+HI(υ)→IH(υ’)+I体系的第一散射共振态寿命, 所得数值与Neumark 的高分辨阈值光分离光谱实验的结果相一致.     

分子阻障内旋转的配分函数

本文用WIGNER-KIRKWOOD微扰展开法来推导分子阻障内旋转配分函数的近似表示式,准确到h[2]项,所得公式只包含两个参量(即分子自由旋转配分函数Qf和势垒高度UO与kT之比值)并可用BESSEL函数来表示,数值计算结果表明,在内旋转转动惯量较大,温度较高这一实用上最为重要的范围内,根据本文中所得公式,对任意的Qf值和U0/kT值,可直接算得阻障势能所产生的熵的修正值(Sir-Srir)/k,比用数值求解法较为简单而准确,如计及微扰展开的高次项,还可以在整个温度范围内取得良好的数值结果。不仅如此,本文所得出的解析公式容易推广到多维运动的情况中去,因此较PITZER的数值求解法具有更普遍的意义。     

用Matlab求解化工原理计算问题的教学初探

以管路计算中的设计型问题和精馏过程计算中的操作型问题为例,分析了求解问题的思路,列出了详尽的Matlab计算源程序。计算过程及结果表明,利用Matlab求解化工原理计算问题具有数值稳定性好、使用方便、高效、快速等特点。     

大型线性方程组的迭代求解

将求解大型标准特征值问题的Davidson方法推广到求解大型线性方程组. 推广后的方法既可以用于求解对称问题, 又可以用于求解非对称问题. 数值计算结果显示, 该方法对这两类问题的求解分别优于共轭梯度方法和双共轭梯度方法.     

可变误差多面体法同时测定镧及其它轻稀土含量

报导了以化学性质极其相似的La、Ce、Pr、Nd、Sm为分析对象, 采用五个已知组成的标样, 通过全选主元的Gauss消去法求得单组分的吸光常数。据此再用可变误差多面体法求解样品中各组分的浓度, 实践证明, 此法不仅减小了实验误差和手工计算时间,而且考虑了组分间的相互作用, 能获得比较准确的结果, 可变误差多面体法也是解决这类问题更为合适的方法。     

计算机求解氰乙基醋酸纤维素酸性水解动力学参数

从氰乙基醋酸纤维素在 pH=0.43的酸性条件下水解后的表观乙酰基含量数据,根据非线性最小二乘法原理,用计算机求解了各个反应速率常数,所得结果与实验值符合较好。     

可变误差多面体法同时测定镧及其它轻稀土含量

报导了以化学性质极其相似的La、Ce、Pr、Nd、Sm为分析对象, 采用五个已知组成的标样, 通过全选主元的Gauss消去法求得单组分的吸光常数。据此再用可变误差多面体法求解样品中各组分的浓度, 实践证明, 此法不仅减小了实验误差和手工计算时间,而且考虑了组分间的相互作用, 能获得比较准确的结果, 可变误差多面体法也是解决这类问题更为合适的方法。     

结晶型聚合物注塑过程冷却模型的数值模拟

对HDPE平板状注塑件的冷却过程进行了数值模拟 .基于小分子结晶性物质的冷却模型 ,建立了固相、液相及固液相界面的控制方程 ,方程中引入热焓E和参数T 后 ,将带耦合条件的热传导方程转变为非线性抛物型方程 ,按控制容积法将其离散后 ,用数学软件Matlab获得了PE的三维温度场 ,这为探讨带相转变的热传导问题提供了一种有效的数值方法 .通过对制件芯层温度变化的模拟值与实验值的比较 ,比较了不同相变温度下理论结果的准确性 ,结果表明 ,在引入一定过冷度后 ,将相变温度定义在低于平衡熔点以下时 ,该方法的模拟值具有更好的准确性 .同时 ,模型分析了不同工艺条件对制件芯层冷却过程的影响 ,通过相变平台出现的先后以及终止时间的不同 ,说明充模温度主要影响熔体阶段的冷却以及相变出现的先后 ,而模具温度主要影响固体阶段的冷却以及相变平台终止的先后 .     

SO2基态势能面和振转能级的理论研究

在键长-键角坐标下,精确求解SO₂的核振动方程,并通过与实验数据比较来优化势能参数,由所得势能面计算得到38个振动能级,与实验值相比,均方根误差为0.93cm^-1.计算了³⁴SO₂的部分振动能级以及³²SO₂的J=6以下的部分振动能级,所得结果均与实验值较为吻合     

聚乙烯熔体的离散松弛时间谱与熔融指数的关系

该文通过动态线性小振幅剪切震荡实验所得的聚乙烯熔体储能模量Ｇ′（ω）和耗能模量Ｇ″（ω）数据,采用最小二乘法线性回归、正则法和非线性回归法分别计算得到离散松弛时间谱,比较三种方法得到结果的差别,讨论计算参数、温度和材料对离散松驰时间谱影响．     

自动调节光谱响应的高精度多组分分析研究: 人工神经网络分光光度法

本文论证了人工神经网络(ANN)对多组分光度分析的应用。探讨了人工神经网络算法及其影响因素。用此法做五组分体系(维生素B1, B2, B6, C和菸酰胺)紫外光度分析, 并与CPA方法作比较。该方法的特点在于可由调节不同光谱范围的线性及非线性响应的计算过程来自动校正。研究结果证实了该方法的优越性。     

金属化合物价电子结构的多体理论研究

本文发展了一种在量子化学计算得到的分子轨道基础上求解金属化合物价电子体系多电子Scbrdinger方程的方法,利用这种方法对一系列金属化合物的电子光谱进行了计算,所得结果与实验大体符合。     

凝并和成核机理下颗粒尺度分布的Monte Carlo求解

颗粒的凝并和成核现象影响其尺度分布,现有的MonteCarlo方法描述颗粒尺度分布的时间演变过程存在若干困难.提出了一种新的多重MonteCarlo(MMC)算法,基于时间驱动,利用加权的虚拟颗粒的思想,在模拟过程中保持虚拟颗粒总数不变和计算区域体积不变.利用该算法对“常凝并核,一阶成核”的情况下颗粒尺度分布的时间演变过程进行了数值求解,所得结果与数值解相符,表明MMC算法具有高且稳定的计算精度.另外,MMC算法由于跟踪比实际颗粒数目少得多的虚拟颗粒而具有较低的计算代价.     

OCS电子基态势能面与振动光谱的理论研究

本文采用键长-键角内标系下的自洽场-组态相互作用方法精确计算了OCS分子的振动高激发态能级,并结合实验观测到的振动能级利用非线性最小二乘法优化电子基态势能函数中的势能参数。由优化所得的势能面计算出的振动激发态能级与50个实验观测到的振动能级比较,标准偏差为0.08cm^-^1。此外,还用该势能面计算了OCS同位素分子的振动能级,计算结果与实验值也十分吻合。     

泛函迭代法研究等同平行平板双电层间的相互作用

利用泛函分析理论中的迭代法, 分别计算了等同平行平板双电层在128、77、和25 mV电位下的相互作用能. 并以数值法所得结果为参照, 在各电位下分别与Debye-Huckel(DH)线性近似法、Langmuir 近似法所得的结果进行比较. 结果表明, DH线性近似法和Langmuir近似法均只能分别局限于极低或极高电位, 而泛函迭代法不但有简单的解析表达式, 而且在各种电位下都能得到较满意的结果.     

用分光光度法测定络合物真正组成及其稳定常数的一种新的图解法

测定络合物真正组成的目的、意义以及用分光光度法来判别单核与聚核络合物的方法等,文献中已有报导。这些方法中最主要的一种类型是基于稀释对络合物离解度的影响而建立的数学模型。本文提出的方法虽然也是属于这一类型的,然而与其不同的是,此法是把不稳定常数K当作诸多自变量的函数而建立的一种新的数学模型,由此可推出一个不定方程式,求解便可得到络合物的真正组成及其稳定常数。本法推导严格,应用简便,结果准确。     

两种C~6~0双炔衍生物的结构、光谱和二阶非线性光学性质的 理论研究

用INDO/2和INDO/SCI方法,计算了C~6~0C(C=CH)~2基态电子结构和电子光谱,所得结果与实验结果基本一致.在此基础上,设计了C60(C=CH)~2,并用同样方法计算了电子结构和光谱.在正确的电子光谱基础上,用ZINDO-SOS方法计算了2个分子的二阶非线性光学系数β~i~j~k和β~μ,并对其结果进行了分析和讨论.结论是,乙炔基与C~6~0(C=CH)~2有更大的非线性光学系数,是有希望的非线性光学备选材料分子。