斜对角笛卡尔方法在流数值模拟中的应用.docxVIP

斜对角笛卡尔方法在流数值模拟中的应用 1 种新的不法边界处理方法—前言 在节水、环境和海洋工程的研究中，经常使用数学模型方法来模拟水流动力学。在计算水流数学模型时，由于不规则的复杂边界，如何处理完整的流量数学模型，直接影响整个模型计算的可靠性、精度和计算时间。这是值模拟的一个难题。选择复杂边界的处理方法与算法的选择有着同等重要性,目前在水利、海洋和环境工程中常用的不规则复杂边界处理方法主要有贴体坐标法、非结构性网格法和锯齿法。贴体坐标法由 Thompson(1974) 提出,该方法的主要优点是用曲线坐标逼近复杂边界,在边界模拟上具有较高的精度。但该方法在较复杂的边界上不易形成正交贴体网格;如果是非正交曲线坐标系,则会产生大量的交叉导数项,给求解带来很大的困难。非结构性网格法使用三角形和四边形计算单元,使得在复杂区域形成计算网格更加灵活和方便 (Baker,1983)。但非结构性网格法往往要建立一套标明网格节点之间关系的复杂系统,因此在编程计算处理时比较麻烦,特别是针对自由面流动和干湿边界变化时,这种困难更加突出。传统的笛卡尔坐标方法被应用的历史最长,由于其简单,易形成计算网格、正交等特点,在工程界应用得最为普遍,但传统的笛卡尔坐标方法在边界处形成锯齿状边界,导致边界的拟合精度不高,因此该方法在大尺度计算中应用较多(Chai, et al.1994)。本文介绍一种新的不规则边界处理方法—斜对角笛卡尔方法。该方法不仅保留了锯齿法正交、易形成计算网格、计算效率高等特点,同时又能有效地改善锯齿法拟合边界精度差的缺点,更能贴近实际,使边界流场更加合理。 2 斜对角基卡尔方法 一般而言,在笛卡尔坐标系统中,按常规方法对边界逼近,易形成锯齿状边界,从而导致各计算变量的误差,但利用斜对角笛卡尔方法后,则可显著改善计算的精度(Lin, et al, 1997; Fang, et al, 2000)。 2.1 拟合网格的形成 斜对角笛卡尔方法基本原则就是按就近逼近原则,把距实际水边线较近的网格节点连线以拟合实际水边线。生成网格时,根据边界是否固定,可采取不同的判定方法。当边界为固定边界时,可把固定边界离散为若干离散节点,确定各离散节点坐标,按距离最小原则选定固定边界上各离散节点的逼近网格节点,然后依序连接各网格节点即形成拟合边界;对于动边界,由于边界线变化不定,则可借助网格节点间水深的干湿变化和水深值,参照相应网格节点间的地形高差,即可判定距离实际水边线较近的网格节点,合并纵横向的拟合网格节点,即可形成拟合瞬时水边线。以某一正方形网格计算域为例(见图1),线L为实际水边线,点A、B、C、D、E、F、G、H为实际水边线与网格的交点,在斜对角笛卡尔法中,以距交点最近为原则,选取对应网格节点作为实际水边线的逼近点,根据这一原则,选取网格节点A′、B′、C′、D′、E′、F′、G′、H′分别逼近实际水边点A、B、C、D、E、F、G、H,依此类推,获得所有实际水边线与网格线交点的逼近网格节点,把各逼近网格节点依序连线便形成拟合水边线L′。 2.2 斜角栽培卡方法的优点 2.2.1 按表面活性剂法确定的水边线 由图 2a 可以看出,计算水边线与实际水边线拟合良好,其边界长度误差一般小于6%,仅为锯齿法长度误差的1/5 (Lin, et al,1997),且随网格间距的缩小而更为精细。图 2b 为按锯齿法 (也以距离最近为原则) 确定的计算水边线,对比图 2a、b 可知斜对角笛卡尔法所拟合的水边线比锯齿法拟合的水边线更为贴近实际水边线。事实上,斜对角法就是在锯齿法基础上通过斜对角连线添加一三角形网格单元,使边界形态更接近真实形态,且避免了锯齿状边界的产生,使流场计算更精确。 2.2.2 斜加角加卡尔法边界条件的应用 在矩形网格水域边界上运用斜对角笛卡尔法,计算时,只须对三角形边界网格单元进行处理,而三角形边界网格单元若把斜边中点视为一角点,则仍可视为四边形网格,可用统一方程求解,与内部计算网格无异,并不象贴体坐标拟合一样须对基本方程进行坐标变换而增添许多交叉导数项,因此,用斜对角笛卡尔法进行边界处理,不仅保留了离散简便、计算效率高的优点,同时又使边界拟合更加精细,避免了锯齿状边界对流场计算的影响。 3 边界条件及边界条件 水流数学模型基本控制方程如下 水流连续方程: ?Ζ?t+?(ΗU)?x+?(ΗV)?y=0?Z?t+?(HU)?x+?(HV)?y=0 水流运动方程: 式中: Z——水位 H——水深 f——哥氏(Coriolis)系数 U、V——分别为x,y向流速 C——谢才阻力系数 νt——紊动粘性系数 本文计算采用有限差分的迎风和中心差混合格式离散,根据涨退潮流向选取边界迎风方向,边界条件采用下节中所介绍的斜对角笛卡尔方法处理,动边界处理采用水面延伸法