流体示踪质点轨迹线的计算方法与动态显示技术研究.pdfVIP

．1392． 第六届全国泥沙基本理论研究学术讨论会论文集 一_———————————————————————————————————————————————一 流体示踪质点轨迹线的计算方法与动态显示技术 汤立群1 申锦瑜2 (1．中国水利水电科学研究院北京 100044； 2．Dalhouse 4J1) University，Halifax，NS，Canada，B3H 摘要：本义简要介绍了三重双向嵌套三维环流数值模型，简单分析了模型计算结果及模型计算区域水 流运动的基本特征。在此基础上，着重介绍了流体示踪质点轨迹线的计算方法。以模型计算结果为水流 动力条件，采用具有4阶精度的Runge—Kutta方法计算出给定质点的轨迹线，并以此模拟悬浮颗粒在不 同水深处的运动轨迹。提出颗粒保留率概念，并计算分析了小模型计算域(Inner)内4个给定子区域颗粒 保留率的年内变化过程，对工程泥沙的研究具有实际意义。利用Gri等专业绘图软件生成并播放悬浮颗 粒在流体中的运动过程。 关键词：示踪质点 轨迹线 悬浮颗粒保留率 动态显示 1前言 流体质点的轨迹可以直观地反映悬浮颗粒或无游泳能力的生物幼体在流体中的运动 状况，有利于对流体和固体物质运动机理的认识。通常用数学模型计算得到的模拟结果 描述的是Euler物理量场，它可以清楚地描述各点处矢量在不同时刻的变化，但无法描 述固定质点随时间的变化，即质点的轨迹。在许多研究场合下，人们希望能将Euler场 转化为Lagrangian场以便直观地描述流体质点或泥沙颗粒的输移规律。本文利用三维数 值模型结果，采用具有4阶精度的Runge—Kutta方法计算出给定质点的轨迹线，并以此 模拟悬浮颗粒或无游泳能力阶段的生物在不同水深处的运动轨迹，可为河流泥沙工程规 划设计或渔业管理提供直观的数字图形，并利用Gri等专业绘图软件生成并播放泥沙颗 粒或生物幼体在流体中的运动过程。 2模型及模拟结果 2．1模型简介 在考虑风应力、热通量、海表面温盐场及开边界海流输送作用的基础上，建立了三 重双向嵌套海洋环流三维数值计算模型【61，其基本特点是：大一级模型的变量被插值到 小一级模型的网格上以提供小一级模型的边界条件；小一级模型变量同样被插值到公共 区域大一级模型的网格上作为模型的信息反馈；用OSPM的方法，通过大一级模型和小 作者简介：汤立群，男，1962年出生，教授级高级工程师。河口海岸研究室主任，主要从事河口海岸与河流泥沙 模拟、土壤侵蚀与流域产沙及水土保持研究。 十、河流模拟技术 ·1393· 一级模型的动量方程来交换信息(温盐场)。模型被成功地应用于加勒比海西海岸地区浅 水珊瑚礁生态系统水流动力条件的模拟计算。模型的计算结果如体积输送流函数、多年 平均海表面流场、垂直断面流速及海温等值线、计算域各点不同水深处的海温随时间变 化过程等都得到了满意的验证效果【1】。 2．2模拟结果 模型被应用到加勒比海西海岸浅水地区，用来模拟计算珊瑚礁生态系统的水流动 力连接条件。模型连续计算6年，因最初几个月模型的计算稳定性较差，通常第一年 计算结果不予采纳，后5年的模拟结果被用来分析沿岸珊瑚礁地区的水流特征和温盐 场分布。 图1表示的是加勒比海表层(5m)处的月平均流场，2月、5月、8月和11月分别表 示冬、春、夏和秋4个季节，反映了该地区表层流场的季节性变化特征。在每一幅图中， Caribbean Sea，网格尺寸为19 Outer、Middle和Inner分别表示大(Western km)、中 BarrierReefs (Meso·．american SystemMBRS，网格尺寸为6 MBRS， km)