解非线性抛物型方程的多层网格扰动迭代法

提出一种数值求解非线性抛物型方程初边值问题的多层网格扰动迭代法；该方法有效地结合了多层网格方法和扰动迭代方法，在固定的时间网格层上该方法有二阶敛速，渐近最优；整体计算量为Ｏ（ＭＮｌ），其中Ｍ是时间计算层数目，Ｎｌ是空间分划细网层节点变量个数；计算误差不传播，且解决了迭代初值的选择问题。     

解非线性抛物型方程的多层网格扰动迭代法

提出一种数值求解非线性抛物型方程初边值问题的多层网格扰动迭代法；该方法有效地结合了多层网格方法和扰动迭代方法，在固定的时间网格层上该方法有二阶敛速，渐近最优；整体计算量为Ｏ（ＭＮ＿ｔ），其中Ｍ是时间计算层数目，Ｎ＿ｔ是空间分划细网层节点变量个数；计算误差不传播，且解决了迭代初值的选择问题。     

非线性椭圆问题的紧致差分格式及瀑布两网格法

为了讨论来源于科学工程问题的二维非线性椭圆问题的离散格式及其数值解法。首先,将泊松方程的四阶紧致差分格式推广至二维非线性椭圆问题,提出了紧致差分(CFD)格式,基于CFD格式,选取合适的步长,形成粗网格层和细网格层。在粗网格层上,使用牛顿法求得对应的非线性方程的高精度数值解;在细网格层上,运用插值算子将粗网格上的数值解进行插值,得到细层上较好的初始值,并再次使用牛顿法进行求解,提出了CFD格式下的瀑布两网格(CTG)法。数值实验表明:提出的CFD格式具有四阶计算精度,CTG法迭代步数少、计算时间短。     

薄壳动力分析的三维半显式迭代算法

利用三维变分差分方法研究薄壳的动力分析。针对显式迭代格式最大稳定时间步长过小，而隐式迭代格式计算量大且精度不足这一问题，构造了一种半显式迭代格式（即关于厚度方向隐式、而关于其余两个方向显式），它的最大稳定时间步长较显式迭代格式有很大的提高，而计算量并未显著增加。算例的数值结果表明，这种半显式迭代格式具有较高的精度，它的时间步长满足计算薄壳波动问题的要求。     

HLLC方法在二维非结构动网格上的应用

以自适应非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的Godunov一阶精度方法求解Euler方程，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，空间和时间都是一阶精度，对含有运动边界的多体干扰流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步进行多体分离数值模拟打下了基础。     

校正Cahn-Hilliard方程的混合有限元两层网格方法

针对校正Cahn-Hilliard方程的非线性、四阶导数以及小参数等特点,提出将混合有限元法与两层网格法相结合的混合有限元两层网格方法;该数值方法由2步完成,第1步在粗网格上用隐式混合有限元方法求解一个四阶非线性系统,第2步在细网格上求解2个线性系统,然后给出所提方法的稳定性分析与收敛性证明,并通过数值实验对理论分析进行验证。结果表明,理论与实际算例结果相一致,并在计算过程中达到了降阶与缩短计算时间的目的,说明了所提方法的有效性与可行性。     

四阶Adams-Bashforth组合公式的预估-校正方法的对比试验

本文研究了由4阶显式的Adams-Bashforth公式与同阶隐式的Adams-Moulton,Hamming和Gear公式组合构造了预估-校正方法,对它们进行了数值对比试验,获得了Adams-Bashforth-Hamming预估-校正方法比其它两种方法的计算结果稳定.     

基于二维泊松方程六阶紧致格式的多重网格方法

利用六阶紧致差分格式、结合多重网格V循环算法求解了二维泊松方程的Dirichlet边值问题，并用不同的松驰算子与四阶精度格式的多重网格方法进行了比较，计算结果表明，该方法在不明显增加计算量的前提下较四阶精度格式的多重网格方法具有更好的精确度和收敛阶，且ZLGS迭代不论对四阶精度还是对六阶精度格式的多重网格算法，都是一种较其他松弛算子更加有效的“光滑剂”。     

基于自由液面预测的非线性液体晃动问题的数值模拟

针对一类矩形贮液箱的非线性晃动问题,提出一种基于有限差分算法的数值模拟方法.该方法首先采用σ变换将流体区域转化为矩形域;然后在时间迭代的过程中,运用自由液面运动学条件的显格式预测出下一时层的自由液面形状,并对边界条件中的非线性项作近似处理,由此得到描述下一时层各状态量的离散线性方程;最后根据数值解来对下一时层的自由液面进行校正.使用该方法计算了自由晃动、受迫晃动和行波3种情况,数值计算结果表明该方法具有良好的收敛性和准确性,能够有效地实现非线性晃动过程.此外,这种预测校正格式能够在保证算法稳定性的前提下,灵活地处理包含运动边界的数学问题.     

Cahn-Hilliard方程的时间双层网格混合有限元方法

针对数值求解Cahn-Hilliard方程时非线性项引起的时间耗时问题,提出了时间双层网格混合有限元方法.首先,在时间粗网格上,通过非线性牛顿迭代方法求解非线性混合有限元系统,其中空间离散采用混合有限元方法,时间离散采用隐式欧拉格式;其次,基于初始迭代数值解和拉格朗日插值公式,在时间细网格上求解线性混合有限元系统;最后,分析了该方法的稳定性和误差估计,并通过数值算例进行验证.结果表明,与传统的混合有限元方法相比,该方法可以节省计算时间.     

一维动边界问题的变时间步长方法

本文介绍了求解一维常物性动边界问题的MVTS方法,对其迭代方法作了较为彻底的改进,大大减少了计算时间。文中还将VTS方法推广到一维变物性动边界问题,提出了直接计算法,与目前工程上较为通行的变换法作了比较,结果表明直接法比变换法有更大的优越性,可以大大减少计算工作量。     

电机通风系统风量问题的新型迭代解法

研究了含风压源的复杂风网中风量的求解方法，提出了统一的风网矩阵模型和迭代公式，并给出了两个算例。结果表明，在迭代初值任意选取的情况下，经过几次迭代就可使各支路的风量收敛于稳定值，模型适用于复杂通风结构电机的风量计算。     

一种迭代未知输入观测器设计

针对具有未知输入的线性系统讨论了基于迭代控制思想的未知输入观测器设计方法.针对每一步迭代,由于当前未知输入已知,因而可以通过设计常规的Luenberger观测器得到当前步的状态估计和输出估计.在此基础上采用D-型迭代学习控制的思想,通过当前输出步误差和未知输入估计值提出下一步的未知输入的迭代估计方法,并利用此次未知输入估计值再次设计Luenberger观测器估计下一步的状态.之后,对迭代算法的收敛性进行了分析,并给出了迭代算法收敛的充分条件.最后,对一个实际模型进行仿真,验证了算法的有效性.     

二维抛物型方程初边值问题拟多重网格预处理迭代法

将求解二维椭圆方程边值问题的拟多重网格预处理迭代法推广到二维抛物型方程中去,采用Crank—Nicolson格式来离散二维抛物型方程．由于网格节点顺序对迭代格式的构造至关重要,因此对每一时间层上的Z层网格节点按照旋转红一黑序进行排序．数值试验表明,此方法迭代次数较SOR法有明显减少,迭代解与精确解的误差值相对较低,收敛速度较快．因此,在求解二维抛物型方程初边值问题中拟多重网格预处理迭代法是一种很有效的方法．     

中子输运问题源项反演的反幂法

对于包括裂变反应在内的中子输运源项反演问题，研究关于源项有效倍增因子的惫一特征值问题的求解．基于球谐函数展开和有限差分离散，给出了中子输运方程的源项反演逼近的反幂算法，该方法的优势是在适当的初值条件下可以显著提高计算速度．计算结果表明，在对有效倍增因子有较好的预先估计值的情况下，反幂法迭代3步，误差就为0．04545％，而乘幂法迭代20步，误差为0．109％，由此可以看出反幂法计算速度更快，计算结果更精确。     

带子矩阵约束的二次逆特征值问题的最小二乘埃尔米特广义斜哈密顿矩阵迭代解

针对带子矩阵约束的二次逆特征值问题的最小二乘埃尔米特广义斜哈密顿结构矩阵解问题,给出了一种共枙梯度迭代算法。首先提出了带子矩阵约束的二次逆特征值问题的最小二乘问题及其最佳逼近问题;然后分别给出了基于共轭梯度的迭代算法,证明了算法的收敛性。对于任意初始约束矩阵,在不存在舍入误差的情况下,用该迭代算法可以在有限步迭代中得到迭代解。最后,给出了一个数值实例,数值实例证明了所提算法的有效性。     

一种电力网络潮流计算的综合算法

针对牛顿-拉夫逊法对初值要求严格,迭代速度快的特点,利用电力网的结构特点,使用高斯-塞得尔迭代法的第一次迭代结果作为牛顿-拉夫逊法的计算初值。这样既解决了牛顿-拉夫逊法对初值要求高的问题,又提高了收敛速度。计算结果表明,综合算法在迭代次数和收敛速度上有优势。     

带有固定步长的非单调信赖域方法

该文提出一种新的非单调信赖域方法．当试探步不能被接受时，算法沿着试探步的方向求得下一个迭代点，其中步长利用固定公式计算．这种方法既避免了重复求解信赖域子问题，又减少了线搜索方法计算函数值的次数．该文采用的非单调策略是基于张洪超和Hanger（2004）出的非单调线搜索技术．在适当的条件下，证明了算法的全局收敛性及超线性收敛性．最后给出了初步的数值实验结果．     

一类半线性椭圆方程的二重网格差分算法

应用二重网格差分算法处理了一类半线性椭圆问题。无需求细网格上的非线性解,对粗网格(可以很粗)上的数值解在细网格上进行几次线性修正即可,且重复算法的最后一步可以按粗网格步长任意阶地逼近细网格上的非线性解。算法提高了计算效率但不降低精度,有数值算例加以验证。