国家自然科学基金(10172044) 作品数:6 被引量:13 H指数:3 相关作者: 周春华 伍贻兆 更多>> 相关机构: 南京航空航天大学 更多>> 发文基金: 国家自然科学基金 中国航空科学基金 更多>> 相关领域: 理学 航空宇航科学技术 更多>>
非结构网格上Euler方程有限体积解法的改进 被引量:4 2003年 对二维非结构三角形网格上 Euler方程有限体积解法的格点格式进行了一些改进 ,重点在于提高数值解的精度 ,细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件 ;发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步 Runge-Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果 。 周春华 伍贻兆关键词:非结构网格 EULER方程 翼型 不可压Navier-Stokes方程基于误差估算的网格自适应解法 被引量:1 2005年 首先导出了广义S tokes方程Petrov-G a lerk in有限元数值解的当地事后误差估算公式;以非连续二阶鼓包(bum p)函数空间为速度、压强误差的近似空间,该估算基于求解当地单元上的广义S tokes问题。然后,证明了误差估算值与精确误差之间的等价性。最后,将误差估算方法应用于N av ier-S tokes环境,以进行不可压粘流计算中的网格自适应处理。数值实验中成功地捕获了多强度物理现象,验证了本文所发展的方法。 周春华关键词:有限元 一种新的并行自适应网格有限元算法 被引量:4 2003年 给出了一种新的适用于流体力学问题的并行自适应有限元算法。首先 ,基于初始稀网格上获得的事后误差估算值 ,应用反复谱对剖分方法对初网格进行划分 ,使各子域上总体误差近似相等 ,从而解决并行自适应计算中的负载平衡问题。然后在各处理器上独立地求解整体问题 ,并进行指定子域上的网格自适应处理。最后将各子域上的自适应网格组合成一个整体网格 ,应用基于粘接元技术的区域分裂法在该网格上获得最终解。 周春华 Randolph E Bank关键词:流体力学 网格划分 不可压N-S方程的基于粘接元的区域分裂解法 2004年 首先,简单介绍了基于粘接元的无重叠区域分裂方法.这种方法利用变分原理,非常适合有限元近似.然后,着重讨论了这种区域分裂方法在求解不可压Navier Stokes方程中的应用,具体包括等价变分公式的建立、通过算子分裂的时间离散、区域分裂情形下广义Stokes问题的共轭梯度迭代求解方法、空间的有限元离散.最后,以数值实验结果验证了这种区域分裂方法应用于不可压Navier Stokes方程求解时的可靠性. 周春华关键词:NAVIER-STOKES方程 STOKES问题 等价 共轭梯度 流动数值模拟中一种并行自适应有限元算法 被引量:4 2006年 给出了一种流动数值模拟中的基于误差估算的并行网格自适应有限元算法.首先,以初网格上获得的当地事后误差估算值为权,应用递归谱对剖分方法划分初网格,使各子域上总体误差近似相等,以解决负载平衡问题.然后以误差值为判据对各子域内网格进行独立的自适应处理.最后应用基于粘接元的区域分裂法在非匹配的网格上求解N-S方程.区域分裂情形下N-S方程有限元解的误差估算则是广义Stokes问题误差估算方法的推广.为验证方法的可靠性,给出了不可压流经典算例的数值结果. 周春华关键词:有限元 不可压流数值模拟中基于事后误差估算的网格自适应方法 首先导出广义Stokes方程Petrov-Galerkin有限元数值解的当地事后误差估算公式。以非连续二阶鼓包(bump)函数空间为速度、压强误差的近似空间,该估算基于求解当地单元上广义Stokes问题。然后,将误差估算... 周春华关键词:不可压NAVIER-STOKES方程 有限元 文献传递 虚拟区域法及其在流体力学中的应用 被引量:1 2005年 讨论了一类椭圆型算子Dirichlet问题的一种基于Lagrange乘子的虚拟区域方法;由此导出的鞍点问题用共轭梯度法迭代求解.为加速迭代收敛,构建了合适的预处理器.着重考虑了这种方法在不可压粘性流动数值模拟中的应用.通过基于算子分裂的劋laMarchukYanenko时间离散格式,将虚拟区域情形下的不可压NavierStokes方程分裂成非线性对流扩散方程、准Stokes方程和虚拟区域情形下的线性椭圆型方程三个子问题.给出了绕固定和运动圆二维流动的数值实验结果. 周春华关键词:LAGRANGE乘子 算子分裂 有限元