粘弹性波动方程的H1-Galerkin时空混合有限元分裂格式

构造一维粘弹性波动方程的H¹-Galerkin时空有限元分裂格式.这种新的分裂格式在时空两个方向同时利用有限元离散,具有H¹-Galerkin混合有限元方法和时空有限元方法的优点,如在不受LBB相容性条件限制的同时能够高精度逼近流体的压力和达西速度,有限元空间可以利用不同次数的多项式空间,能同时得到时间和空间两个变量的形式高阶精度等.通过构造时空投影算子并讨论其相关逼近性质,证明了解的存在唯一性和稳定性,给出混合时空有限元解的误差估计,给出数值算例验证了理论推导结果的合理性和算法的有效性,并和传统H¹-Galerkin方法做比较,得到了更小的误差和超收敛阶.     

一个分数阶生长-抑制系统的参数反问题

本文研究一个分数阶生长-抑制线性系统模型及其参数反问题.首先利用Laplace逆变换得到正问题解的唯一存在性.其次,考虑一个利用内点观测数据确定微分阶数与衰减率的反问题,应用极值原理在Laplace像空间中证明反演的唯一性.最后,基于正问题的有限差分解,应用同伦正则化算法进行数值反演.计算结果表明算法的收敛性及反问题的数值稳定性.     

二维振荡叶栅非定常粘性流动数值模拟

采用显式四步Runge-Kutta格式，结合Baldwon-Lomax紊流模型求解Navier-Stokes方程，借助运动网格技术，完成了对二维振荡叶栅非定常粘性流动的数值模拟。为了加速求解过程，引入了变系数隐式残差光顺方法，取得了较好效果。数值结果与已公布的数据有很好的一致性。     

多级轴流压气机叶栅内三维紊流流场的数值模拟

对多级轴流压气机叶栅内三维紊流流场进行数值模拟,采用高精度高分辨率的三阶ENN格式以保证对激波的捕捉和对紊流特征的正确模拟,利用LU-SGS隐式解法提高了计算速度,从而构成了一种既准确又高效的多级跨声速轴流压气机紊流流动数值求解系统。重点研究了动静网格的交接方法及相应的动静交接面处理模型。计算了某带进口导叶的三级轴流压气机,计算结果与实验数据吻合较好。     

对“二阶迎风TVD数值格式在非平衡高超声速钝体绕流中的运用”一文的讨论

对“二阶迎风ＴＶＤ数值格式在非平衡高超声速钝体绕流中的运用”一文的讨论胡光初（中国空气动力研究与发展中心，四川绵阳６２１０００）１９９６１１１１收到．为了方便讨论，将我们于１９９４年３月在《力学学报》第２期上发表的“二阶迎风ＴＶＤ数值格式在非平衡...     

混合层强化混合的数值研究

受 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验的启发，采用高阶精度的差分格式，通过数值模拟的方法，研究了二维混合层及限于两平板间的二维混合层（二维受限混合层）入口处加振动对提高混合层混合效率的作用．计算结果表明：对二维混合层，振动的频率越低，在混合层中产生的大尺度涡结构的尺度越大，在频率很低时，涡具有相似性；对限于两平板间的二维混合层，在一定的振动频率下，混合层中产生的涡较大而且破碎得也较好，这将有利于混合．这一结论与 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验观测到的结果是一致的．     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

三维N-S方程计算隐式有限体积法

利用NND有限差分格式，发展了一种新的完全隐式的有限体积数值方法，以求解与时间相关的N－S方程．对通过单元体界面的无粘流和粘性流通量均作隐式处理．对绕流钝锥体和不同攻角的气动辅助实验飞行器的高超声速粘性流和化学反应流获得了定常数值解．对流加热率和流场电子密度的计算值与实验数据进行了比较，符合较好，证实了本方法的精确性．     

EFK方程一个新的低阶非协调混合有限元方法的高精度分析

对Extended Fisher-Kolmogorov(EFK)方程,利用EQ_1~(rot)元和零阶RaviartThomas(R-T)元建立了一个新的非协调混合元逼近格式.首先,证明了半离散格式逼近解的一个先验估计并证明了逼近解的存在唯一性.在半离散格式下,利用上述两种元的高精度分析结果以及这个先验估计,在不需要有限元解u_h属于L~∞的条件下,得到了原始变量u和中间变量v=-?u的H~1-模以及流量p=u的(L~2)~2-模意义下O(h~2)阶的超逼近性质.同时,借助插值后处理技术,证明了上述变量的具有O(h~2)阶的整体超收敛结果.其次,建立了一个新的线性化向后Euler全离散格式并证明了其逼近解的存在唯一性.另一方面,通过对相容误差和非线性项采取与传统误差分析不同的新的分裂技巧,分别导出了以往文献中尚未涉及的关于u和v在H~1-模以及p在(L~2)~2-模意义下具有O(h~2+τ)阶的超逼近性质,进一步地,借助插值后处理技术,得到了上述变量的整体超收敛结果.这里h和τ分别表示空间剖分参数和时间步长.最后,给出了一个数值算例,计算结果验证了理论分析的正确性.     

一类新的(2n-1)点二重动态逼近细分

利用正弦函数构造了一类新的带有形状参数ω的(2n-1)点二重动态逼近细分格式.从理论上分析了随n值变化时这类细分格式的C~k连续性和支集长度;算法的一个特色是随着细分格式中参数ω的取值不同,相应生成的极限曲线的表现张力也有所不同,而且这一类算法所对应的静态算法涵盖了Chaikin,Hormann,Dyn,Daniel和Hassan的算法.文末附出大量数值实例,在给定相同的初始控制顶点,且极限曲线达到同一连续性的前提下和现有几种算法做了比较,数值实例表明这类算法生成的极限曲线更加饱满,表现力更强.