网格和边界信息文件（fort.14）

网格和边界信息文件（fort.14）

该文件是运行ADCIRC模型所必需的。基本文件结构如下所示。输入数据的每一行由包含输入变量名称（粗体）的一行表示。空行只是为了提高可读性。循环指示输入的多行。条件输入由变量名后的if子句指示。通过热链接提供每个变量的定义。

AGRID

NE, NP

for k=1 to NP

JN, X(JN), Y(JN), DP(JN)

end k loop

for k=1 to NE

JE, NHY, NM(JE,1),NM(JE,2), NM(JE,3)

end k loop

NOPE

NETA

for k=1 to NOPE

NVDLL(k), IBTYPEE(k)

for j=1 to NVDLL(k)

NBDV(k,j)

end j loop

end k loop

NBOU

NVEL

for k=1 to NBOU

NVELL(k), IBTYPE(k)

for j=1,NVELL(k)

NBVV(k,j) ? include this line only if IBTYPE(k) = 0, 1, 2, 10, 11, 12, 20, 21, 22, 30

NBVV(k,j), BARLANHT(k,j), BARLANCFSP(k,j) include this line only if IBTYPE(k) = 3, 13, 23

NBVV(k,j), IBCONN(k,j), BARINHT(k,j), BARINCFSB(k,j), BARINCFSP(k,j) include this line only if IBTYPE(k) = 4, 24

NBVV(k,j), IBCONN(k,j), BARINHT(k,j), BARINCFSB(k,j), BARINCFSP(k,j), PIPEHT(k,j), PIPECOEF(k,j), PIPEDIAM(k,j), include this line only if IBTYPE(k) = 5, 25

end j loop

end k loop

法向磁通边界类型 NORMAL FLUX BOUNDARY TYPES

通过在连续性方程式中指定对法向边界通量积分的贡献（零或非零）并在动量方程式中指定（零或非零）法向速度，可以应用具有自由滑移的基本边界。此边界条件应满足全局意义上的通量平衡和每个边界节点处的法线通量。

通过在连续性方程中指定对法向边界通量积分的贡献（零或非零），并设置（零或非零）法向速度和零切向速度，而不是求解动量方程，可以应用无滑移的基本边界沿边界。此边界条件应正确地满足整体意义上的法向磁通平衡和每个边界节点处的切线速度为零。

通过在连续性方程中指定对法向边界通量积分的贡献为零或非零来应用自然边界。动量
方程式中的速度（法向或切向）没有限制。该边界条件应该在全局意义上正确满足通量平衡，但仅在无限分辨率的限制下满足每个边界节点的法向流量。

正常流边界条件的一般说明：

在列出任何内部（无法向流动或内部障碍的岛屿）边界段之前，应在整个域的外部以连续顺序列出所有外部（无法向外部流动，具有指定法向流动的外部边界和外部屏障）。与外部边界相交的内部障碍边界应单独指定，即使这会导致在两个边界中都指定了某些节点（请参见下文）。

应该通过将第一个节点重复为最后一个节点来关闭外部无法向流或完全围绕域（例如湖泊）的指定法向流边界。

应该通过将第一个节点重复为最后一个节点来关闭所有没有正常流动的内部边界（例如，孤岛）。

除非边界段是封闭的，否则始终开始列出两个边界连接的边界节点。

当外部指定法向流或外部屏障边界连接到外部非法向流边界时，外部指定法向流边界或外部屏障边界的初始分支用于确定
两个边界共同的节点处的法向和切线方向。

具有指定（非零）正常流边界条件和外部屏障边界的外部边界无法连接。它们必须由外部非正常流边界或高程指定边界分隔。

内部屏障边界可以与外部无法线边界相交。（例如，堤坝可能会从外部非正常流动边界伸出，在这种情况下，内部障碍边界上的2个节点（前后节点
）对于外部边界是公共的。）但是，必须处理公共外部节点在弱的意义上。ADCIRC将自动适应以下情况：

如果将滑移（IBTYPE（k） = 0）的外部无流动边界指定为必不可少，并且将滑移（IBTYPE（k） = 4）的内部屏障边界指定为必不可少，则公共外部边界节点将自动更改为自然无流边界节点（IBTYPE（k） = 20）。

–如果将外部无流动边界指定为无滑移的基本条件（IBTYPE（k） = 10）并且将内部屏障边界指定为有滑动滑移的基本条件（IBTYPE（k） = 4），则会自动更改公共外部边界节点到自然无流边界节点（IBTYPE（k） = 20）。

如果将外部无流动边界指定为带滑差的自然边界（IBTYPE（k） = 20），而将内部屏障边界指定为有滑动的插接自然边界（IBTYPE（k） = 4），则不会进行任何更改。

如果将滑差（IBTYPE（k） = 0）指定为外部无流动边界，而滑差（IBTYPE（k） = 24）将内部障碍边界指定为自然边界，则不会进行任何更改。

如果将外部无流动边界指定为无滑动的基本条件（IBTYPE（k） = 10），将内部屏障边界指定为有滑动的自然边界（IBTYPE（k） = 24），则公共外部边界节点将自动更改滑移边界节点（IBTYPE（k） = 0）变为基本无流量。

如果将外部无流动边界指定为带滑差的自然边界（IBTYPE（k） = 20），而将内部屏障边界指定为有滑动滑触的自然边界（IBTYPE（k） = 24），则不会进行任何更改。

内部屏障边界不能与外部指定流边界段，外部屏障边界段或内部无法线边界相交。

对于所有正常流边界（即IBTYPE（k） = 0、1、2、3、4、10、11、12、13、20、21、22、23、24、30），连续性中的边界通量积分用适当的（零，指定或计算的）通量评估方程式。这是自然边界条件。对于自然法向流边界（IBTYPE（k） = 20,21,22,23,24），这是唯一使用的横向边界条件。

对于具有切向滑移（IBTYPE（k） = 0,1,2,3,4,10,11,12,13）的基本法向流边界，法向动量方程（通过重新定向x / y动量方程获得）消除法向角和切线方向），并通过将每单位宽度（零，指定或计算）的法向通量除以总水柱高度来设置法向速度。

对于没有切线滑移的基本法向流边界（IBTYPE（k） = 10,11，12,13），两个动量方程均被消除。切线速度设置为零，法向速度通过将单位宽度（零，指定或计算）的法向通量除以总水柱高度来设置。使用此边界条件需要格外小心，因为严格来说，这种类型的无滑动边界条件只有在模拟中使用了横向粘性项时才在数学上是合理的，而在横向边界层得到充分解析时则仅在物理上是合理的。

外部障碍边界注：（ IBTYPE（k） = 3，13，23）

垂直于外部屏障边界并通过外部屏障边界的每单位宽度的流出流量QN2（k，j）计算为：

情况1水位低于或等于屏障高度

QN2（k，j）= 0

案例2水位高于屏障高度

QN2(k,j) = -(2/3)*BARLANCFSP(k,j)*RBARWL*((2/3)*RBARWL*G)**0.5

where, RBARWL = ETA1(NBVV(k,j))-BARLANHT(k,j) = water height above the barrier

ETA1(NBVV(k,j)) = water level computed at the previous time step at node NBVV(k,j)

This formula is given by Leendertse (Aspects of SIMSYS2D ? A System for Two-Dimensional Flow Computation, Rand/R-3572-USGS, 1987) and is simply the formula for a broad crested weir (e.g., see Henderson, Open Channel Flow, section 6.6).

See also General Notes for Normal Flow Boundary Conditions

内部屏障边界注：（ IBTYPE（k） = 4，24）

内部屏障边界由一个细长的孤岛组成，该孤岛具有平行的正面和背面。成对的节点放置在边界的任一侧，以便在正面和背面的节点之间提供一对一的对应关系。假定流量从一个节点到另一侧的配对节点跨越边界。边界两侧的法向流量大小相等，符号相反（例如，正面的流出=背面的流入）。
内部屏障边界节点NBVV（k，j）及其配对节点IBCONN（k，j）的每单位宽度的正常流量QN2（k，j）计算为：

案例1水位低于或等于屏障两侧的屏障高度

QN2（k，j）= 0

案例2水位高于屏障高度，但屏障两侧相等

QN2（k，j）= 0

情况3的水位高于屏障的高度，但在亚临界水流越过屏障的情况下，正面比背面大。如果背面的屏障上方的水位高度大于正面的屏障上方的水位高度的2/3（即RBARWL2> 0.667 * RBARWL1），则会发生从前向后穿过屏障的亚临界流。

QN2（k，j）= -RAMP * BARINCFSB（k，j） * RBARWL2 *（2 * G *（RBARWL1-RBARWL2））** 0.5

案例4的水位高于障碍物高度，但正面越过障碍物，正面越多。如果背面的障碍物上方的水位高度小于或等于正面的障碍物上方的水位高度的2/3（即，RBARWL2 <0.667 * RBARWL1），则会发生从前向后穿过障碍物的超临界流。 ）。

QN2（k，j）=-（2/3）* RAMP * BARINCFSP（k，j） * RBARWL1 *（（2/3）* RBARWL1 * G）** 0.5

情况5的水位高于屏障的高度，但在亚临界水流穿过屏障时，背面的水位大于前面的水位。如果正面上的障碍物上方的水位高度大于背面上的障碍物上方的水位高度的2/3（即RBARWL1> 0.667 * RBARWL2），则会发生从后到前穿过障碍物的亚临界流。

QN2（k，j）= RAMP * BARINCFSB（k，j） * RBARWL1 *（2 * G *（RBARWL2-RBARWL1））** 0.5

案例6的水位高于屏障高度，但背面超前，而超临界流则越过屏障。如果正面上的障碍物上方的水位高度小于或等于背面上的障碍物上方的水位高度的2/3，则会发生从屏障到背面的超临界流动（即RBARWL1 <0.667 * RBARWL2 ）。

QN2(k,j) = (2/3)*RAMP*BARINCFSP(k,j)*RBARWL2*((2/3)*RBARWL2*G)**0.5

where

RBARWL1 = ETA1(NBVV(k,j))-BARINHT(k,j) = water height above the barrier on the front side of the barrier

RBARWL2 = ETA1(IBCONN(k,j))- BARINHT(k,j) = water height above the barrier on the back side of the barrier

ETA1(NBVV(k,j)) = water level computed at the previous time step on the front side of the barrier

ETA1(IBCONN(k,j)) = water level computed at the previous time step on the back side of the barrier

These formulae are given by Leendertse (Aspects of SIMSYS2D ? A System for Two-Dimensional Flow Computation, Rand/R-3572-USGS, 1987) and are simply the formulae for a broad crested weir (e.g., see Henderson, Open Channel Flow, section 6.6).

See also General Notes for Normal Flow Boundary Conditions

带交叉障碍管的内部障碍边界注：（IBTYPE（k） = 5,25）

此类型与IBTYPE（k） = 5或25不同。此类型与IBTYPE（k）= 4和24不同之处在于，屏障包含具有指定高度，管道系数和管道直径的交叉屏障管道。在技​​术出版物“ ADCIRC沿海水动力法典中的内部泄漏屏障正常流边界 ”中详细描述了这种边界类型的公式。

因为此边界类型与IBTYPE（k） = 4或24 有很多共同点，所以这里仅描述差异。在具有交叉屏障管节点NBVV（k，j）及其成对节点IBCONN（k，j）的内部屏障边界处，每单位宽度的正常流量QN2（k，j）计算为：

情况1：内部屏障两侧的水位低于横向屏障管顶部的高度

QN2（k，j）= 0

情况2：内部屏障两侧的水位相等

QN2（k，j）= 0

情况3：内部屏障正面的水位高度大于背面屏障的水位高度；前侧的水位高度大于交叉障碍管的顶部高度；背面的水位低于管冠高度

QN2(k,j) = -RAMP*(0.25*pi*D^2)*(2*G*RBARWL1/(1+PIPECOEFR(k,j)))^0.5

where

RBARWL1 = ETA2(NBVV(k,j))-PIPEHTR(k,j)

情况4：内部屏障正面的水位高于背面的水位；前侧的水位高度大于交叉障碍管的顶部高度；以及水位高于管道顶部高度的背面。

QN2(k,j) = -RAMP*0.25*pi*PIPEDIAMR^2*(2*G*(RBARWL1-RBARWL2)/PIPECOEF(k,j))^0.5

where

RBARWL1 = ETA2(NBVV(k,j))-PIPEHTR(k,j)

RBARWL2 = ETA2(IBCONN(k,j))-PIPEHTR(k,j)

情况5：内部屏障背面的水位高于正面的水位；背面的水位高度大于交叉障碍管的顶部高度；正面的水位低于管冠高度。

QN2(I)= RAMP*0.25*pi*(PIPEDIAMR)^2 * (2*G*RBARWL2 / (1 +/PIPECOEF(k,j)))^0.5

where

RBARWL2 = ETA2(IBCONN(k,j))-PIPEHTR(k,j)

情况6：内部屏障背面的水位高于正面的水位；背面的水位高度大于交叉障碍管的顶部高度；且水位高于管道顶部高度。

QN2(I)= RAMP**0.25*pi*(PIPEDIAMR(k,j))^2*(2*G*(RBARWL2-RBARWL1)/PIPECOEFR(k,j))^0.5

where

RBARWL1 = ETA2(NBVV(k,j))-PIPEHTR(k,j)

RBARWL2 = ETA2(IBCONN(k,j))-PIPEHTR(k,j)

可以在以下参考资料中找到此边界类型的原理和实现的详细信息：

Westerink，JJ，RA Luettich和A.Militello，2001，ADCIRC沿海水动力法中的内部屏障正常流边界泄漏，沿海和水利工程技术说明ERDC / CHL CHETN-IV-32，美国陆军工程师研究与发展中心，密西西比州维克斯堡，2001年2月，第28页。单击此处进行发布。