基于斜对角笛卡尔方法的分蓄洪区水流数学模型的研究.pdfVIP

基于斜对角笛卡尔方法的分蓄洪区 水流数学模型研究 穆锦斌张小峰张艳霞屈文谦 (武汉大学水赍澡§林电击睁畔学国家重点实嘘蟊j蚓4XXTt2) 籀要分蓄洪区运用前，采用数学模型对区内洪水演进作出预报是一项必不可 少的工作．对于复杂边界问题，本文引进了斜对角笛卡尔方法，较传统的笛卡尔法， 该方法能更好的贴近实际水边线．在此基础上，对原有的平面二维分蓄洪区水流数学 模型进行了改进．经检验，模型能更好模拟复杂边界处水流运动规律，拟合精度较高． 改进后的模型运用到湖北龙骨湖，沉底湖分蓄洪区实际计算中，取得了令人满意的结 果． 关键词分蓄洪区，复杂边界，斜对角笛卡尔方法，数学模型 1前言 随着经济的迅速发展和城市的急刷扩张，越来越多的分洪区己被交通道路、工厂、企 业、乡村甚至城镇所占据，分洪区内的人口、财产不断增长。分洪区一经启用，将造成财 产的重大损失。如何保证人员安全疏散安置，尽可能减小财产损失，是分洪区安全建设研 究重要的课题之一，其中，对分洪区内洪水演进作出准确预报是一项必不可少的工作。平 面¨二维数学模型能较好地模拟水流运动规律¨_：1，再现分洪区内洪水演进情况．它已成为防 御洪水灾害中不可或缺的重要工具。 在二维水流运动的数值计算中，求解区域边界的复杂性给数值模拟带来了很大的困难， 如何处理不规则的复杂边界直接影响整个模型计算的可靠性及精度p州。传统的笛卡尔坐标 方法由于其简单，计算网格正交等特点，已被广泛应用，但此方法在边界处易形成锯齿状 边界，导致边界拟合精度不高。为了增强该方法对非规则区域的适应性，本文引进了一种 新的不规则边界处理方法一斜对角笛卡尔方法16】。该方法不仅保留了传统的笛卡尔坐标方 法的简单、正交、易形成计算网格、计算效率高等特点，同时又改善了锯齿边界拟合精度 差的缺点，更能贴近实际，使边界流场更加合理，对于复杂边界的分洪区洪水演进计算有 较高的精度。 2斜对角笛卡尔方法概述 斜对角笛卡尔方法基本原则嘲是按就近逼近原则，把距实际水边线较近的网格节点连 成线以拟合实际水边线。生成网格时，根据边界是否固定，采取不同的判定方法。当边界为 固定边界时，可把固定边界离散为若干离散节点，确定各离散节点坐标，按距离最小原则 选定固定边界上各离散节点的逼近网格节点，然后依序 连接各网格节点形成拟合边界；对于动边界，由于边界 ／ B’fC’1 ，夕 R≮ 线变化不定，可借助网格节点间水深的干湿变化和水深 B、 C ≮? 值，参照相应网格节点间的地形高差，判定距离实际水 D1 瓤 边线较近的网格节点，合并纵横向的拟合网格节点，形 F1 宴-；边界jEL／ 成拟合瞬时水边线。 如图l所示，线￡为实际水边线，点彳、B、C、 ＼ D、E、，、G、日、，为实际水边线与网格的交点． I、 网格节点∥、口’、cI、D’、E’、F、G’、日’、，’ ∥ 分别逼近实际水边点4、嚣、C、D、E、F、G、