澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法．应用于澳门水域潮流的模拟,计算结果与实测吻合较好．     

求解具有复杂地形二维浅水方程的修正HLL格式

采用非结构网格的有限体积方法,对具有复杂地形的二维浅水方程进行数值模拟.采用HLL近似Riemann解计算界面数值通量,基于三角形网格,底坡源项采用简单的斜底模型离散,保证了地形的离散精度,摩阻源项采用全隐方式求解以保证格式的稳定性.采用多维重构及多维限制器的方法获得空间二阶精度的格式,时间离散采用三阶Runge-Kutta法以获得高阶时间精度.为保证格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.数值计算的结果验证此格式具有良好的高精度捕捉间断的能力,可以应用到地形复杂的二维浅水问题计算中去.     

一维和二维溃坝波的混合有限分析解

针对水深平均二维浅水方程高分辨率模型,采用非交错网格下的二维对流扩散方程的混合有限分析法,运用SIMPLE算法结合Rhie and Chow的动量插值技术计算二维溃坝水流．一维理论解和其他格式的二维数值结果符合良好,表明非交错网格下的混合有限分析法能正确有效地求解溃坝问题．进而对部分溃决情况下的溃坝波传播和绕流问题进行了数值模拟,给出了精细的数值结果,揭示了复杂的运动特性,进一步表明混合有限分析法具有较强的溃坝波的捕捉能力,是求解溃坝流动的有效方法之一．     

黏塑性流体运动的二维浅水方程模拟

利用垂向平均浅水方程模型描述黏塑性流体的运动,并基于剪切流动区切应力分布与明渠均匀层流切应力分布相等的假定,推导了浅水方程中切应力项的摩阻系数表达式.建立了黏塑性流体运动平面二维非结构化网格数值模型,空间离散采用有限体积法显式求解,时间积分采用具有二阶精度的Runge-Kutta格式求解.数值模型模拟结果与已有黏塑性流体流动实验结果进行了比较,二者有较好的一致性.     

控制体有限元方法中质量加权辽风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散，为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式，通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法地回流问题数值模拟的     

控制体有限元方法中质量加权迎风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散．为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式．通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法对回流问题数值模拟的假扩散．     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

控制体有限元方法与流向型插值格式

给出在同位网格上对二维定常不可压层流数值模拟的控制体有限元方法及其流向型插值格式的计算过程,并对几个层流流动问题进行数值计算。     

高精度WENO和MWENO格式溃坝波模拟

采用高精度WENO和MWENO格式对二维浅水方程变量进行数值重构,研究分析高精度格式捕捉溃坝波的能力和分辨率.对浅水方程,采用Roe-Riemann格式和波传算法进行数值离散,并与以上2种重构方法两两结合成4种计算方法.采用Fortran语言分别对其进行编程计算.计算结果显示,高精度格式在捕捉溃坝瞬时基波延续推进的过程较好.在4种计算方法计算结果中,MWENO型波传格式模拟结果最好,捕捉溃坝波能力强,分辨率很高.     

高精度格式溃坝波模拟

采用高精度WENO和MWNEO格式对二维浅水方程变量进行数值重构,研究分析高精度格式捕捉溃坝波的能力和分辨率.对浅水方程,采用Roe-Riemann格式和波传算法进行数值离散,并与以上2种重构方法两两结合成4种计算方法.采用Fortran语言分别对其进行编程计算.计算结果显示,高精度格式在捕捉溃坝瞬时基波延续推进的过程较好.在4种计算方法计算结果中,MWENO型波传格式模拟结果最好,捕捉溃坝波能力强,分辨率很高.     

二维无结构网格浅水波方程的高阶非振荡有限体积法

基于二维浅水波方程,对无结构网格给出了一种三阶精度非振荡有限体积方法,方法的主要思想是先对每一个三角形单元采用最小二乘的思想构造一个二次插值多项式,而在计算交界面的流通量时采用了两点高斯公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对圆柱形溃坝问题和倾斜水跃问题进行了数值试验,得到了满意的结果。     

求解带源项双曲守恒律方程的IGVC格式

研究带源项双曲守恒律方程的IGVC(改进的群速度控制)格式。为了控制数值方法带来的非物理振荡,引入IGVC格式,考察其TV性质,对于源项我们采用分裂方法处理,该方法用于求解一维浅水方程,数值结果是令人满意的。     

2维Ginzburg-Landau方程的分裂LOD高阶紧致格式

采用分裂技巧研究了2维的Ginzburg-Landau方程构造高效的数值格式.把2维Ginzburg-Landau方程变成线性和非线性问题以避免求解耦合的非线性方程组.为减少存储量和计算量,对线性问题进一步运用局部1维方法,把它分解为2个1维问题求解.所得到的数值格式具有高效、高精度等数值特征.最后,用数值算例模拟了2维Ginzburg-Landau方程所描述的物理现象,新方法具有较大的优越性.     

基于高阶谱差分和嵌套网格方法的三维无粘流动数值模拟

数值方法采用高阶谱差分格式,单元边界上的通量采用基于Roe格式的近似黎曼解,在嵌套网格上求解三维Euler方程.嵌套边界上使用了一种高阶紧致的插值格式.对三维亚音速定常无粘流动进行数值模拟,使用嵌套网格与单块网格得到的计算结果吻合很好.     

二维浅水流动数学模型的多解性评价

二维浅水流动的多解性现象经常呈现于实际工程中．数值模拟业已成为二维浅水流动分析的主要手段．所采用的控制方程也必须能正确反映水流运动在该方面的本质．与零方程模型相应的含紊动粘性项二维浅水方程组具有多解性;无对流项和紊动粘性项的二维浅水方程组解具有唯一性;有对流项而无紊动粘性项的二维浅水方程组只有离散时引入人工粘性项解才是有限的;与双方程模型相应的含紊动粘性项二维浅水方程组具有多解性也是必然的．     

风生环流的准三维数值模拟

本文在水平流速垂向分布表达式采用二次多项式的基础上,对二维浅水方程组中的动量方程进行修正,建立了风生环流的准三维数学模型.应用破开算子法求解考虑动量修正的深度平均方程,对破开后的两个分步分别采用隐式格式和预估-校正格式离散,得到了准三维数值解.验证计算表明该数值模式能用于模拟河口,水库和港湾内的风生环流,算法的稳定性好,收敛快.     

基于GLS／SUPG的三维注射成形充模过程数值模拟

为稳定地求解速度、压力、温度并准确跟踪流动前沿,实现三维塑料注射成形充模过程数值模拟,采用迦辽金最小平方法(GLS)消除因速度和压力的不当插值组合导致的数值震荡问题,建立了速度与压力同次插值的对称、稳定的集成有限元计算格式;用流线迎风Petrov伽辽金法(SUPG)消除因能量方程中的对流项而导致的数值震荡问题,建立了能量场的稳定有限元计算格式.比较了不同网格密度对计算结果的影响,并通过与文献中实验结果及现有的商业分析软件Moldflow的模拟结果对比,表明提出的数值模拟方法具有良好的稳定性和精度.     

一类二维非线性发展方程的交替方向有限体积元方法

采用交替方向有限元方法的离散思想,利用差分方法中的Douglas格式,给出了求解一类二维非线性发展方程的基于双线性插值的交替方向有限体积元法.文中对该方法进行了误差估计,并用数值例子验证了方法的有效性.     

QUICK和乘方格式在顶盖驱动方腔流动数值计算中的比较

推导了非均分网络系统中QUICK格式附加源项的表达式，并用QUICK和乘方两种差分格式，分别对不同网格数下的二维、三维顶盖驱动方腔流动进行了数值计算。通过与文献中基准解的比较，考察了两种对流差分格式的数值预测性能。     

新型高分辨率格式及其在CFD的应用

 利用小波奇异分析构造了一种新型的高分辨率高效计算格式,借助小波奇异分析所找到的流场中奇异点区域和光滑区域,合理选择高分辨率、高精度WENO格式和高精度中心差分格式进行流场的离散求解;利用该格式对超声速二维前台阶流动问题、跨声速RAE 2822典型翼型的二维外流场、跨声速VKI-LS59涡轮叶栅的二维流场以及NASA Rotor 37跨声速压气机转子和NASA Rotor 67跨声速风扇转子的三维流场进行数值模拟,计算结果表明,采用这种新型格式可以用于内外流复杂流场计算。数值计算还发现,在采用相同网格的情况下,本文格式不仅具有通常WENO格式高分辨率的特点,而且计算效率较WENO格式高出3～5倍。