差分模型网络嵌套技术在工程潮流计算中的应用.pdfVIP

差分模型网格嵌套技术在工程潮流 计算中的应用 辛文杰 (南京水利科学研究院，210024) 擒要 本丈对差分模型的嗣格嵌套技术作了介绍，提日{流速梯度项疣f／西是不司忽视的边界控制 条件，并给出了采用线性插值获取如／尻’值的处理方法，通过实际算例证明嗣格嵌套技术在实际工 程潮流计算中的作用。 美键词差分模型嵌套边界 1前 言 在对河口或海湾的平面潮流场进行数值计算时，人们常采用有限差分格式。差分模型 具有格式简单、计算效率高和稳定性好等优点，其缺点主要是难以适应复杂多变的地形边 界，网格剖分的灵活性较差，尤其是对1：程边界的模拟，只能以缩小空间步长为代价获得 对局部边界的近似概化，结果往往使计算网格大大增加，时间步长也随之缩小，使模型的 计算效率明显降低。 为了克服这一缺陷，有人把差分格式改为变步长，也有人通过坐标变换实现边界的拟 合等等。这些处理方法都是可行的，但同时也增加了计算的复杂性和数据处理的麻烦。本 文提出一科比较简单易行的网格嵌套技术，可以在保持原有差分方程面貌的情况下，通过 政套边界的参数传递，实现大范围疏网格整体模型与小区域细网格局部模型的有效连接， 从而既满足了边界控制条件，又达到了1[程概化要求，并且操作起来比较灵活。该方法在 ～些工程潮流的模拟计算中已取得丁令人满意的结果。 2格嵌套域的确定 对于需要加密细化的局部模型来说．如何确定合适的计算域是首先要考虑的问题。一 般而言，其开边界的位置应在基本不受工程边界反馈影响的地方为宜，过远则降低计算效 率，过近则需要给定1二程后的边界控制条件，而这类条件通常是未知的，唯一的途径只有 从嵌套域外的整体模型中获取。因此，在确定局部模型的计算域时，应事先在整体模型中 根据简化的工程边界计算出其影响范围及幅度，再根据研究区域的地形地貌特点和其它技 术要求，划定模型范围。如果确定的计算域开边界仍在工程影响范围之内，划需要针对不 同的工程布置，从整体模型中取得相应的边界条件。 一，，{一 3嵌套网格的划分 嵌套在整体模型中的局部模型，剖分网格的：k小只跟所需概化的工程尺度有关a密集 的网格可以更恰当地模拟出工程建筑物的几何形态，反之则差一些，与整体模型网格的大 小没有直接关系。但在实践中，粗、细网格间距比不宜过大，比值Jo△J／As‘10较为台 适(As’代表粗网格，△J代表细网格L 嵌套网格的坐标方向应与整体模型保持～致，否则在进行边界条件的处理和传递时会 遇到许多困难。图1是嵌套模型网格划分的简单示意。 ，一， —’7、j ．， ＼ |、 ＼ f 擀#＼ p’ ＼ V＼一 图l崴套网格划分不意图 4嵌套边界的控制 4．，边界控I时羹件 嵌套在整体模型内的局部模型，在开边界上仅考虑潮位或流速控制条件往往是不够的， 应增加边界上流速梯度的控制条件。原因是： (1) 由于局部模型区域较小，由流速形成的惯性水头比潮波比降形成的重力水头 对边界控制的作用更为重要： (2) 局部模型开边界的流速并不能f{i}证平行于边界法线方向，而且纵向流速也随 力 C吖 r为边界切向)； 因此，嵌套边界的控制条件t应包括潮位f，流速“和流速梯度af／加才是完备的。 图2是在某海域采用嵌套模型计算潮流场时有无巧k／两边界条件的对比情况。(a)是边界无 砌／两条件，局部模型与整体模型在对应点的流速过程偏差较大；(b)是增加了巩／斑控 制条件，计算和实测的流速过程吻合较好。该算