二维浅水水流的一种新型三角形网格FVM计算格式

将三角形网格良好的边界拟合性能与有限体积法守恒性质好、求解速度快的优点相结合，构造了适用于三角形网格的有限体积计算模式，模式以交替方向隐式(ADI)求解．采用长江三峡工程河段水流资料对模型进行验证计算，将实测值与计算值作了对比，结果令人满意．     

非结构四面体网格上三维非线性扩散方程的有限体积法

为了有效地数值模拟科学和工程中有广泛应用的非线性扩散方程,在三维线性扩散方程非结构四面体网格的有限体积法的基础上,提出了一个计算非结构四面体网格上非线性扩散方程的有限体积法。方法采用网格单元中心作为计算节点,相对于网格点的方法,计算量减少了一半。用L agrange因子法得到网格点上的值,考虑了网格中心点和网格点的相对位置,更适应大变形的网格。利用算子分裂,使计算更加简单。用N ew ton-B iCG STAB法来求解得到非线性方程组。数值结果表明:该方法具有二阶精度、保持通量守恒、对大变形的网格适应性强。     

多重网格法与有限体积法在流体力学中的应用

多重网格技术是一种非常有效的计算机数值计算方法，本文采用多重网格中ＦＡＳ格式进行数值实验，计算加速效果十分明显；同时，结合矢量及有限体积法，大大提高了计算机捕捉激波的质量。     

非正交曲线同位网格三维水动力数值模型

建立一个自由表面三维水动力学计算模型，平面采用贴体坐标变换，竖直方向σ坐标变换，基于有限体积法，同位网格布置即流速，压力布置于同一控制点上，自由表面由二维水深积分方程中的水位方程求解，运用SIMPLEC法求解速度－压力耦合问题，动水压力校正方程校正压力和流速场，标准k-ε紊流模计算紊动粘性系统以封闭Reynolds应力项，运用该数值模糊模拟S形弯道水流三维流动。     

弹性力学中三种有限体积法方案的对比

为了研究弹性力学中三种典型的有限体积法方案(FVM1,FVM2和FVM3)的精度差异问题,采用理论推导和算例验证的方法,建立了基于有限体积法原理的力平衡积分方程,利用三角形网格进行离散化,在计算出内部单元和边界单元对力平衡贡献的基础上,对每一种有限体积法方案确立了与有限单元法方案(FEM)刚度方程组相类似的线性方程组,对比所建立的各个方程组的待定系数.结果表明:方案FVM2与FEM精度相同,即为最优方案;在不考虑集中力的情况下,方案FVM3的精度也与FEM相同.最后以一个数值算例,对这三种有限体积法方案进行了对比分析,验证了上述分析.该结果对工程计算数值模拟具有一定的参考价值和指导意义.     

有限体积法实现变系数椭圆型方程数值解

研究了变系数椭圆型偏微分方程的有限体积法,该方法将研究区域划分为一系列不重复的分割区域,并且每个网格点都包含在一个分割区域,再用待求的偏微分方程对每个分割区域进行积分,便可得到一组离散方程。基于这些离散方程,采用matlab编程达到数值实现的目的。最后,通过数值实例展示了有限体积法的计算精度,并得出了一些普遍且有益的结论。     

流体动力学中不连续源项大Reynolds数问题的有限体积法

讨论不连续源项大Reynolds数问题,引入2种非均匀网格技术,应用有限体积法构造了数值方法:Sh-ishkin型有限体积法和多过渡点型有限体积法.数值实验显示,2种方法都较好地模拟了边界层和内部层的性质,具有一致收敛的计算效果;多过渡点型有限体积法明显优于Farrell所提出的Shishkin有限差分法.     

一种Lagrange格式及重映算法

从积分形式的二维Lagrange流体力学方程组出发，用有限体积格式进行计算，考虑压力梯度分布对速度和能量改变的影响，构造了在两个控制体上的动量方程的计算格式。在重映算法上，采用积分重映的方法，针对不规则的四边形网格，根据非结构网格的ENO插值的思想，构造线性插值多项式，由于利用其插值点自适应选取的特性，在物理量变化剧烈区域选取最光滑区域的点来插值，虽然不能保持单调性，但只允许出现非常小的振荡。数值结果表明了该方法的可行性。     

TVD有限体积法在喷流计算中的应用

该文根据ＷａｎｇＪ．Ｃ．Ｔ．的二阶ＴＶＤ有限体积法求解守恒型Ｅｕｌｅｒ方程，对原文中的数值格式进行了部分修改。计算网格单元是非均匀的任意四边形。文中计算了楔形风洞中的斜激波反射问题和无粘喷流近场流动结构问题。数值计算结果表明，该文方法具有激波分辨率高、计算稳定、伪振荡小等优点。同时该文方法可以处理任意复杂边界的二维及轴对称流动问题的计算，也可以将该文的有限体积法思想推广到三维具有复杂边界的流动的计     

高层建筑物风场和风压力的数值计算

针对高层建筑风场和风压力问题，选用了标准Ｋ－ε紊流模型和交错网格的有限体积法，对两项高层建筑物风场和压力场进行了数值模拟，其结果与经验计算结果进行比较，数模结果比较满意。     

四阶抛物型积分-微分方程三角网格上的混合有限体积元方法

采用混合体积元方法在三角网格上求解一类四阶抛物型积分-微分方程的初边值问题,构造了问题的半离散混合体积元格式,得到了误差估计结果．     

椭球外区域问题的自然边界元法

主要讨论了椭球面外部区域Laplace方程的自然边界元法.首先引入椭球坐标,通过分离变量法导出了Poisson积分公式和自然积分算子的无穷级数显示表达式.这样原无界区域问题就归化为椭球面上边界积分方程,然后再数值求解该积分方程.给出了该积分方程的变分问题的适定性和逼近解的误差估计,且该误差估计不仅依赖于网格参数而且依赖与级数截断后的项数.数值例子说明了该方法的有效性和理论分析的正确性.     

偏微分方程高效高精度数值方法研究

针对偏微分方程类型的最优控制问题、多孔介质渗流驱动问题、地下水流的非线性反应扩散方程、对流占有的对流扩散方程、Volterra积分微分方程等阐述混合有限元方法高精度后处理技术、具有超收敛性质的计算格式和高效自适应网格局部加密算法;扩张混合有限元快速收敛的两层网格算法;迎风差分格式的高效自适应移动网格算法;具有高精度的谱...     

矩形网格抛物型问题的质量集中有限元方法

就一类典型的抛物型问题——热传导方程,研究矩形网格上质量集中有限元方法的有关性质.首先给出了矩形单元上双线性有限元基函数的积分公式,在此基础上讨论质量集中有限元方法的误差估计.研究表明,矩形网格上的质量集中有限元方法具有与普通的有限元方法同等的逼近精度,但却具有更少的计算量,并且在一定条件下可以保持极值性质.最后给出了在矩形网格上质量集中有限元方法保持极值性质的剖分条件.     

一种模拟岩体裂纹扩展的三角单元网格开裂技术

基于三角网格的几何特征,提出一种利用有限元方法模拟岩体裂纹扩展的三角单元网格开裂技术。该方法选取三角网格进行单元离散,采用远场围线积分计算裂尖应力强度因子,由最大周向应力准则确定裂纹扩展方向,最后通过开裂单元的网格分裂或节点移动,实现裂纹扩展的数值模拟。以有限宽中心裂纹板、曲线翼型裂纹扩展和含孔洞多裂隙岩体的裂纹扩展为例进行模拟验证。结果表明:在该方法中,裂纹可以直接劈开一个单元,或沿单元边界扩展,因此裂纹能够不受初始网格的限制沿任意路径扩展;与现有的网格重构算法相比,该方法只须对裂尖局部单元进行网格开裂或节点移动,更加简便、高效,该方法还具有较好的适用性,能够准确模拟拉伸、压剪等复杂应力状态下的裂纹萌生和扩展。     

笛卡尔切割单元法在湖泊二维流场模拟中的应用

针对湖泊几何边界复杂的特点,为了高效并准确地计算流场,采用了基于有限体积法的笛卡尔切割单元法.该方法在流场主要区域采用笛卡尔网格,而对与固体边界相交处的单元采用特殊的处理.这样既保留了笛卡尔网格求解速度快的优点,同时在边界处又保证了精度,避免了传统差分法处理该类问题需采用阶梯形网格逼近而导致的失真问题.通过与一个标准算例的贴体网格计算结果对比,显示计算结果一致.最后运用笛卡尔切割单元法,求解了杭州西湖的二维流场,结果表明该方法适用于具有复杂几何边界的湖泊流场求解,且精度较高.     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

激波诱导圆柱运动数值模拟

该文采用基于含有运动边界自适应非结构网格点的显式有限体积方法求解Euler方程，提出了一种运动边界的数值处理方法，对激波诱导的圆柱运动进行了数值模拟。由于物体的运动，物体边界附近的网格将发生变形。对运动边界附近严重变形的网格采用局部网格重新生成的方法进行处理。采用HLLC方法来处理运动边界问题，计算并比较了气流流过静止圆柱和圆柱在静止流场中运动2种情形。结果表明，应用HLLC方法处理运动边界是可行的。     

柴油机涡流室内流场计算的有限容积法

阐述柴油机涡流室内空气流场数值计算的有限容积法的主要内容，并运用有限容积法对柴油机压缩冲程中涡流室内流场进行数值计算．且将计算值与激光多普勒（ＬＤＡ）测量实验值进行比较，结果吻合