近岸区波浪破碎引起的沿岸流计算.pdfVIP

近岸区波浪破碎引起的沿岸流计算 李绍武 (天津大学建筑工程学院水港系300072) 柴山知也 (日本国横滨国立大学工学部土木工学科240) 摘要 本文诣在探讨一种具有广泛适用性的、基于Bau锚i∞sq方程的计算近岸区破波带内沿岸流计 算的数值方法．文中讨论了数值计算中的边界处理问题和底部摩阻对沿岸流的影响，提出了用“数值水 泵”加缓冲区来处理衡边界的方法。诣在减小侧边界对沿岸流分布的影响，获得均匀沿岸流分布．并减 小破波带以外回流。文中运用Viss丘(1991)的实验室资料对沿岸流计算结果进行了对证． 关键词Boassinesq方程波浪破碎沿岸流数值水泵 l前言 (1972)从理论上解决了，其公式在世界范围得到应用。在工程实践中，经常遇到弯曲岸线、 底部地形比较复杂．特别是近岸区存在水工结构物的情况，给沿岸流的计算带来困难。 在有些情况下甚至无法运用该公式。因此，工程上需要一种适用性更广泛的方法。这里 应力的空间分布，即动量流在一个波周期内平均后的余量，也是二阶项。因此，从理论 et a1． 时．计算模式中边界条件的设定也需要进一步的研究。由于理论上的困难，数值计算中 常将侧边界设为全反射边界，这势必会对计算域内的波浪场和流场产生影响。本文提出 了一种减小侧边界影响的方法，同时讨论摩阻系数对沿岸流的影响。 2 理论模式 2．1基本方程 Boussinesq方程为适用于浅水区波浪变形、且具有二阶精度的波动方程。经典的 一349— 守恒型二维Boussinesq方程可表示为(李绍武等，1999) 挈+孥+孥：o， (1) 护r 护z 矛y ∞ 箬+旦ex㈤kO)+岳㈢+驴瓦@r／芬Q厄虿=一争磊≥， 篝+去(铡D+云吲D+妒筹Y略g厕一13炉黯。毋t 毋x＼ j ay＼ )。 a 谢砷、… 号妒r ㈣ 式中，奴和g分别为x和Y方向的单宽流量分董，r／为瞬时水位，^为静止水深， D=^+q为实际水深，g为重力加速度，L为阻力系数n 方程(2)和(3)右边项为考虑水质点垂向加速度而导出的色散项，区别于Airy波理 论。近年来为拓展方程在深水区的应用，还出现了不同的方程形式。其本质区别是对 方程色教特性的处理。本文重点在破波带附近，属浅水区，故以上述经典形式为基础 进行研究，其结论可延伸到其他形式的Boussinesq方程。 为了模拟波浪在近岸区破碎后的衰减，需要在基本方程中引入动量衰减项。有关 兹就破波模型的主要结论简述如下。方程(2)和(3)中引入动量衰减项后成为 等+丢[矧+南[警卜筹专以孵 ㈤ 。土蚓+土型一土。：．血 P 0l p 0Y 3 8zat2 辔+去(引+寿[譬)憎篝+争岛厕㈤ ：三判+三剑一土^：宝1 分量．可表示为 吾吨【筹+矧一扛， ㈣ 弘c，譬， (7) 式中，i，为垂向平均紊动粘滞系数，f为垂向平均紊动动能，巳为系数，巧。为Kxon∞ker 算符，U，和U』为垂向平均流速分量。 紊动粘滞系数i，和紊动能量f通过求解如下々方程确定。 (s，、。 jt