本文描述了如何在水槽中生成规则的床面形态,可视化流体在床面形态中的流动,并利用计算机模拟来仿真潜流带水流。计算机模拟结果与实验观测结果吻合良好。这种耦合模拟与实验的方法非常适合用于科研和教学目的。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
本文描述了如何在水槽中生成规则的床面形态,可视化流体在床面形态中的流动,并利用计算机模拟来仿真潜流带水流。计算机模拟结果与实验观测结果吻合良好。这种耦合模拟与实验的方法非常适合用于科研和教学目的。
沉积物孔隙空间与上覆水体之间的对流交换,在河流环境中称为潜流交换(hyporheic exchange),是驱动河流中溶质输运及许多重要生物地球化学过程的关键机制。为了通过可视化手段增进对这些过程的理解,我们在多智能体计算机建模平台 NetLogo 中构建了一个潜流流动模拟系统。该模拟展示了虚拟示踪剂流经覆盖有二维床面形态的河床的过程。模型将沉积物特性、水流条件和床面形态特征作为输入变量。我们展示了如何利用基于实测参数的实验室水槽实验观测结果来验证此类模拟的准确性。实验中将染料注入水槽的沉积物内,以可视化孔隙水的流动路径;在模拟中,虚拟示踪粒子被放置于相同位置以便对比。这种结合模拟与实验的方法已成功应用于本科及研究生教学实验中,使学生能够直观地观察河流与孔隙水之间的相互作用,并理解基于物理原理的流动模拟如何重现环境现象。学生通过透明水槽壁拍摄床面形态的照片,并将其与模拟中相同时刻的染料分布形态进行比较,结果显示出高度相似的变化趋势,从而帮助学生更深入地理解流动模式及其背后的数学模型。此外,该模拟还允许用户通过运行多个情景快速可视化各个输入参数的影响。这一方法也可用于科研工作中,用以阐释基本过程、关联界面通量与孔隙水输运关系,并支持基于过程的定量建模研究。
当地表水在溪流、河流或潮间带流动时,会产生水力梯度,驱动水流进入和渗出沉积物1在河流系统中,发生这种交换的河床沉积物部分被称为河床下渗带2,3该区域非常重要,因为许多营养物质和污染物在河床下渗带中被储存、沉积或转化。4-9示踪剂在沉积物中停留的时间称为滞留时间。滞留时间以及水流路径的位置均会影响转化过程。为了预测河流中溶质的迁移并解决由营养物质等物质扩散所引发的重大环境问题,需要更深入地理解影响沉积物流水过程的机制。例如, 近海缺氧10,11)。尽管潜流交换具有重要意义,但在本科阶段的水文学、流体力学和水力学课程中,往往未对其进行详细描述, 等等。 希望在课程中加入潜流交换内容的教育工作者,可能会发现使用能够清晰展示该过程的实验和数值可视化方法非常有益。
河道蜿蜒度、周围地下水位以及河床地形(即 河床沙坝、床面形态和生物成因土丘)均在不同程度上影响潜流带交换过程12-17。本研究聚焦于沙丘和波纹等床面形态,这类地形通常是影响潜流带水流的关键地貌特征14,15。我们通过数值模拟与实验室实验相结合的方法,对流经一系列规则床面形态的水流过程进行可视化研究。该模拟基于一系列前期研究成果,这些研究将潜流带水流路径与系统中易于观测的特征参数建立了关联15,18-21。由于这些研究构成了本模拟的科学基础,以下简要概述该理论的核心内容。床面形态地形 T(x) 由下式给出:
公式 1:

其中,H 为床面形态振幅的两倍,k 为波数,x 为平行于平均河床表面的纵向维度。图1 展示了此类床面形态地形的一个示例。

图 1. 用户定义和控制的参数设置。 在界面中,示踪粒子以通量加权的方式在水体/沉积物界面处释放,并在沉积物中被追踪。如果显示路径?设为“开启”,则水体中的示踪粒子会标记其经过的位置,从而显示其运动路径。当示踪粒子返回上覆水体时,系统的示踪粒子总数会发生变化,此时若重新投放?设为“关闭”,则该变化将被记录。累积滞留时间分布图通过绘制沉积层中剩余示踪粒子数量与初始数量之比随时间的变化来反映这一过程。如果重新投放?设为“开启”,则离开系统的示踪粒子将按照与初始粒子相同的通量加权方式被重新补充,此时累积分布图将被禁用。请点击此处查看该图的高清版本。
| 参数名称 | 单位 | 定义 | 界面 | 鼠标投放 | ||
| Lambda (λ) | cm | 床面形态的波长(见图1) | ![]() | ![]() | ||
| BedformHeight (H) | cm | 床面形态振幅的两倍(见图1) | ![]() | ![]() | ||
| BedDepth (D) | cm | 沉积物的深度(见图1) | ![]() | ![]() | ||
| HydrCond (K) | cm/s | 水力传导率 | ![]() | ![]() | ||
| Porosity (θ) | 孔隙度 | ![]() | ![]() | |||
| ChannelVelocity (U) | cm/s | 地表水或河道中的平均流速 | ![]() | ![]() | ||
| Depth (d) | cm | 水深(见图1) | ![]() | ![]() | ||
| Slope (S) | 床面形态和水面的坡度 | ![]() | ![]() | |||
| NumParticles | 释放到系统中的粒子数量 | ![]() | ||||
| TimeX (Time1, Time2..) | min | 每次颜色变化发生的时间 | ![]() | |||
| 模拟按钮 | 定义 | 界面 | 鼠标投放 | |||
| Setup | 使用所示参数设置模拟 | ![]() | ![]() | |||
| go/stop | 启动和停止模拟 | ![]() | ![]() | |||
| Step | 点击“Step”将推进一个时间步长。这允许用户放慢程序运行速度,精确观察100秒内发生的情况。 | ![]() | ||||
| clear paths | 清除屏幕上所有蓝色粒子轨迹 | ![]() | ![]() | |||
| Advance to next time | 使程序运行至下一个颜色变化时间(TimeX) | ![]() | ||||
| mouse-drop | 必须先点击此按钮,才能通过点击地下位置在地下放置粒子。 | ![]() | ||||
| show-paths? | 如果“show-paths?”为“开启”状态,水粒子将留下蓝色轨迹以显示其经过的路径(见图1)。 | ![]() | ![]() | |||
| re-drop? | 如果“re-drop?”为“开启”状态,每当粒子离开系统时,将以通量加权的方式重新投放粒子,此时累积图将无法正常工作。当粒子离开潜流带时,若“re-drop?”为“关闭”状态,则系统中的粒子数量会减少(见图1)。 | ![]() | ||||
表1. 潜流带参数与模拟控制. 本表列出了用户可调节的每个参数、按钮和滑块,并附有其定义。
在本模拟中,有两种过程在砂层中引发流体速度。第一种是水流与床面形态相互作用所致。由床面形态在水-沉积物界面处引发的速度水头也近似呈正弦变化,且相对于床面形态本身相位偏移四分之一个波长22。界面处速度水头函数的振幅已通过测量估算为16:
公式 2:

其中 U 为水面平均流速,g 为重力常数,d 为水深(见图1)。流速水头函数由此给出:
公式 3:

然后可利用该水头函数,通过求解具有恒定砂层深度的拉普拉斯方程,计算基于床面形态的渗流速度函数分量20。孔隙水速度的第二个分量由系统坡度 S 决定,该坡度对应于重力水头梯度,导致顺流方向的流动与
成正比。最终的孔隙水速度函数为:
公式 4:

公式 5:

其中 u 是纵向速度分量, v 是垂直速度分量, K 是沉积物的平均水力传导系数,是沉积物的平均孔隙度 y 是纵坐标,而 D 是沉积物的深度。
创建了粒子追踪模拟,该模拟使用 NetLogo 建模语言和模拟平台23。这两种实现方式(Mousedrop.nlogo 和 Interface.nlogo)采用相同的模拟核心,并利用这些方程对潜流带水流进行建模。主要区别在于示踪粒子的初始位置。Mousedrop 允许用户在地下介质中的任意位置放置模拟示踪剂。通过地下流速方程 4 和 5 驱动示踪剂运动,以模拟染料注入实验。而在 Interface 中,示踪剂始终以通量加权的方式放置于地表/地下界面处。这种设置模拟了溶解态和悬浮物质从地表水向孔隙水输送的过程,这对于理解潜流带交换至关重要。随后,示踪剂在地下介质中运移,直至再次返回地表水体。在水槽中追踪染料路径,并利用 NetLogo 模拟其路径,只要观测期间流动条件和床面形态保持稳定,即可获得流场的流线图。Interface.nlogo 可生成累积滞留时间分布,该分布以时间为函数,显示在任意时刻仍保留在沉积物中的示踪粒子数量与初始时刻放置的示踪粒子总数之比。
如近期一项文献综述所述24,教育研究领域对于动手实验与模拟实验及计算机模型的相对优势仍存在相当大的争议。一方面,一些人认为“动手实践是学习的核心”25,并警示称,以节约成本为由的论点可能会推动用计算机模拟取代动手实验活动,从而不利于学生对知识的理解26另一方面,一些科学与工程教育领域的研究者认为,模拟实验至少与传统的动手实验同样有效27,或探讨计算机模拟在促进以学生为中心的“发现式学习”中的优势28尽管尚未达成共识,但许多研究人员认为,理想情况下,计算机模拟应作为动手实验的补充,而非取代动手实验29,30科学与工程教育领域也已开展了一些举措,旨在将物理实验和真实世界的感知与现象的计算机模拟同步结合;参见, 例如“双焦点建模”31.
通过与物理系统及其对应的计算机模拟系统进行互动,学生可以更深入地理解科学概念以及科学研究过程。本实验要求学生开展溶质运移实验,演示由重力和河床形态引起的潜流交换流现象,并将他们自身的实验装置与结果与同一现象的计算机模拟进行比对。这种对比有助于实现重要的学习目标,促进对科学方法的深入讨论,以及对模型/理论构建与通过数据收集进行实证验证之间相互作用的理解。在完成该对比后,学生还可利用计算机模拟的优势,通过调整模型参数,快速探索多种不同的替代情景。
访问受限。请登录或开始试用以查看此内容。
1. 模拟软件
2. 水槽演示
3. 模拟


访问受限。请登录或开始试用以查看此内容。
结合实验使用模拟,可使学生观察理想化的数学模型与更复杂的实际系统之间的相似性和差异性。 图4 展示了染料注射照片与之比较的示例 小鼠坠落 模拟。初始照片用于确定时间零点时模拟染料示踪剂的位置,随后运行34.2分钟的模拟,并与该时刻拍摄的照片进行比较。总体而言,模型在该时间段内对染色水运动的再现效果非常出色。位于床面形态背流侧的第一个染料团块,在模拟和实验系统中均从沉积物中逸出。第二个染料团块发生拉长并向下移动,在扩散过程中形成新月形,使得部分示踪剂从原始位置下游逸出,另一部分则向上游逸出。最后一个染料团块向上游传播,且部分示踪剂进一步深入沉积物内部。这表明在床面形态下确实发生了潜流交换,且潜流交换流的模式与床面形态的几何特征相关。模拟与实验结果的高度一致性在一级近似水平上验证了模型方程的有效性。该过程还明确表明,潜流交换是一个显著的过程,其规模与床面形态的大小成比例,且几乎有一半的孔隙水在床面形态下方向上游流动。然而,仔细观察后仍可发现观测到的染料输运与模拟结果之间存在细微差异。模拟结果比实际染料分布更平滑,且渗入沉积物的深度较浅。这些差异源于测量误差以及由不规则床面形态几何结构、沉积...
访问受限。请登录或开始试用以查看此内容。
水槽演示与颗粒示踪模拟相结合,为不同背景的受众全面介绍了潜流带水流过程。各级参与者均可通过直观方式观察到由河床形态引起的潜流带水交换现象,以及河床形态下方地下水流路径的显著变异性。这些方法既可作为本科生或中小学学生孔隙水流动的简单演示,也可在研究生课程中结合河流水力学、泥沙输运及潜流带交换机制的深入讲解而使用。无论教学层次如何,利用这一简单的可视化模型作为互动技术,均能使学生对这些复杂而重要的现象形成比单纯依靠抽象理论和讨论更为深入的理解。
在使用这些方法时,物理系统与模拟之间的差异不应被视为“错误”,而应视为一种“可教学的契机”。 即,这是引发深入讨论并最终促进学习的起点。应引导学生思考一系列问题,包括:误差的来源有哪些(来自模型、测量以及实验操作过程)?其中哪些因素可能导致模拟结果与实际观察之间的差异?在构建模型时做出了哪些简化假设?微小差异的重要性如何?它们是否意味着模型是“错误的”?正如统计学家乔治·博克斯(George Box)的名言所说:“本质上,所有模型都是错误的,但其中一些是有用的。”34 一个优秀的科学模型能够捕捉到系统中的...
访问受限。请登录或开始试用以查看此内容。
作者无任何利益冲突需要披露。
本材料基于美国国家科学基金会资助的项目,资助号为 EAR-0810270、EAR-1215898 和 EAR-1344280,以及一项 NSF 研究生科研奖学金。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 水槽 | 工程实验室设计 | 定制 | 带有透明侧壁的实验室水槽,用于持续运行24-48小时。或者,可根据我们补充材料中提供的指南,以低于300美元的成本建造小型教学用水槽。 |
| 流量计 | Rosemount | 8800涡街 | 安装在水槽的循环回路内部 |
| 沙子 | 美国硅砂 | F30 | 实验级沙子,用于在水槽内铺设10-20厘米厚的沙层 |
| 染料 | 食品公司样品 | 水溶性食品级染料,配制成水溶液。深色如红色、蓝色和绿色效果最佳。(避免使用丙二醇基的食用色素。) | |
| 注射器 | HSW | 4100.000V0 | 5-10 ml,例如 HSW Norm-Ject 两件式一次性注射器 |
| 移液针头 | Cadence Science | 7942 | 14号,6英寸钝头, 用于将染料注入沙层深处 |
| 数码相机 | 任意品牌 | 配备稳固三脚架的数码相机。(可使用延时摄影相机以快速、等间隔采集数据。)我们使用的是Nikon D7000。 | |
| 直尺 | 任意品牌 | 透明材质为佳。 | |
| 卷尺 | 任意品牌 | ||
| Netlogo软件 | CCL | http://ccl.northwestern.edu/netlogo/ | |
| Mousedrop.nlogo | Netlogo Commons | 4259 | http://modelingcommons.org/browse/one_model/4259 |
| Interface.nlogo | Netlogo Commons | 4258 | http://modelingcommons.org/browse/one_model/4258 |
访问受限。请登录或开始试用以查看此内容。