用格子Boltzmann方程模拟浅水波问题

提出了用格子Ｂｏｌｔｚｍａｎｎ方程（ＬＢＥ）模拟浅水波问题的方法．通过无粘气体运动方程与浅水波方程的比较，确定了ＬＢＥ方法中平衡态的形式，使宏观方程与浅水波方程一致．计算了二维浅水波的一个问题，数值结果与精确解做了比较．     

求解浅水波方程的并行物理信息神经网络算法

双曲守恒律方程是一类比较特殊的偏微分方程,其数值求解方法的研究一直是一个热点问题,一个显著特性是即使初始条件是光滑的,其解也可能会发展成间断。浅水波方程作为非线性双曲守恒律方程,由于间断解的存在,其精确求解存在很大困难。针对浅水波方程数值求解问题,本文基于PINN(Physics informed neural networks)反问题网络结构构造新的网络,构造的网络结构包括两个并行的神经网络,其中一个网络与已知状态数据(熵稳定格式加密求出)相关,另一个网络与方程本身相关。利用已知速度数据结合浅水波方程本身求解未知水深,最终通过一些数值算例验证网络的可行性。结果表明,新的网络结构可用于浅水波方程求解,利用速度数据可以较为精确地推算出水深。     

复杂地形条件下WENO-Roe方法在浅水模拟中的改进

在分析浅水方程与二维空气运动方程差异的基础上,在不改变原有浅水方程形式的前提下,提出了局部水位法离散连续方程,并针对动量方程中底坡项提出了更具普遍意义的离散方法,保证了浅水方程离散后的平衡性。通过不规则地形下潮流以及混合流流动的模拟,得到的计算结果符合物理实际,与精确解符合良好,最大相对误差不超过4%,验证了此方法在复杂地形上的平衡性,同时本方法又具有良好的间断捕捉能力与稳定性。     

关于采用高阶Boussinesq方程求解波浪散射问题的注记

本文讨论了采用高阶Boussinesq方程模拟波浪散射时对基本速度变量位置的局部光滑处理方法。通过光滑局部基本速度变量的取值深度，减小其高阶导数项的量值、加快级数收敛速度进而改善模型方程求解深水波浪散射问题的能力。对于底部边界具有一阶导数不连续的情况，通过局部光滑．可以将基本速度变量取值深度尖角转化为圆角过渡，从而改善速度分布。对于其它任意变化的底部边界，为了减少高阶底坡导数项的影响，在曲率和高阶底坡导数项与斜率具有相同量级的情况下亦需要对基本速度变量的取值深度局部光滑。数值计算结果表明本文提出的光滑技术可以很好地改善Boussinesq方程模拟浅水波和深水波在斜坡地形上散射问题的能力。     

非均匀介质散射问题的体积分方程数值解法

将非均匀介质视为某一均匀背景介质中的扰动，可建立用均匀背景介质格林函数作基本解的体积分方程．给出了配置法求解体积分方程的数值方法，首先解得扰动域内各点以速度扰动为权的波场函数，然后回代计算得到观测面上各接收点的散射波场．与边界元法和Ｂｏｒｎ近似法计算结果比较表明该方法具有很高的精度，可得到穿过非     

Boussinesq方程的发展形式与Airy波理论差异性研究

应用势流理论，采用递推函数方法推导出一个新形式的Ｂｏｕｓｉｎｅｓｑ方程。通过对新方程的参数设置，可以讨论出Ｂｏｕｓｓｉｎｅｓｑ方程发展趋势和不同的发展形式。对浅水波动的描述方程，Ｂｏｕｓｓｉｎｅｓｑ方程的发展趋势为适用水深范围的拓展。拓展应用范围的大小则由其方程频散特征向Ａｉｒｙ波频散解逼近程度来决定。而Ｂｏｕｓｉｎｅｑ方程又不同于Ａｉｒｙ波，主要原因是Ｂｏｕｓｓｉｎｅｓｑ方程中含有线性频散项，Ａｉｒｙ波则只是长波首项近似，无线性频散项。其频散特征为精确的线性频散解。对实际水波传播而言，Ａｉｒｙ波理论的局限性是不言而喻的。     

两层流体中二维非线性界面波的演化方程

给出两水平固壁间两层不可压缩理想流体中二维非线性界面波的演化方程，首先建立出这个演化方程，并由此方程在一定条件下得到二维非线性界面长波满足的近似方程，然后从理论上证明这个长波近似方程包含了以下两个描述一阶界面升高的著名的浅水孤立波方程；Korteweg-de Vries(KdV)方程和Kadomtsev-Petvishvili(KP)方程，所得特殊结果与前人的一致，表明所建立的二维非线性界面波演化方程正确且具有一般性。     

基于简单WENO-RKDG的溃坝问题数值模拟

对溃坝问题水流间断面的高精度、高分辨率数值模拟是水动力学的重要内容。简单加权本质无振荡(WENO)限制器由"问题单元"及其相邻单元的解重构"问题单元"的解,从而抑制数值解的非物理振荡,能够很好地模拟间断问题。本文详细介绍了简单WENO限制器的基本原理和过程。将简单WENO限制器-Runge-Kutta间断Galerkin方法应用于二维浅水控制方程的求解中,对二维矩形明渠中大坝瞬间全溃、局部溃塌所致的水流运动进行了数值模拟,并将数值计算结果与理论分析进行了比较。计算结果表明,方法能够清晰地捕捉到溃坝全过程中的间断,没有非物理的振荡现象发生,简单WENO限制器-RKDG方法能够很好地模拟溃坝波的演进过程。     

变深度浅水域中非定常船波

以Green—Naghdi(G—N)方程为基础，采用波动方程／有限元法计算船舶经过变深度浅水域时非定常波浪特性．把运动船舶对水面的扰动作为移动压强直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用．以Series60 CB＝0．6船为算例，给出自由面坡高，波浪阻力在船舶经过一个水下凸包时变化规律，并与浅水方程的结果进行了比较．计算结果表明，当船舶经过凸包时，波浪阻力先增加，后减少，并逐渐趋于正常．同时发现，当船速小于临界速度时(Fr＝√gh＜1．0)，G—N方程给出的船后尾波比浅水方程的结果明显，波浪阻力也比浅水方程的结果有所提高，频率散射必须考虑．当船速大于临界速度时(Fr＝√gh＞1．0)，G—N方程的计算结果与浅水方程差别不大，频率散射的影响可以忽略．     

数值模拟在钱塘江涌潮分析中的应用──Ⅰ.数值计算方法

利用数值模拟的方法对钱塘江涌潮从杭州湾口开始形成、发展直至消失的全过程进行了深入全面的描写．从杭州湾口到钱塘江出口,采用二维圣维南浅水波方程描述水波的运动,而在钱塘江河内采用一维圣维南浅水波方程描写涌潮的发展过程．详细描述用于一维和二维圣维南方程计算钱塘江涌潮的数值计算方法,首次把无结构网格上的NND格式应用于求解二维圣维南方程,并给出了详细的推导过程．对上下游水边界分别采用无反射边界条件和特征线方法,而对于动边界问题本文也给出了相应的处理方法．     

Boundary element calculation of shallow water waves

A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled lattice Boltzmann model for the shallow water‐contamination system

A coupled numerical method for the direct simulation of shallow water dynamics and pollutant transport is formulated and implemented. The conservation equations of shallow water dynamics equations and the convection–diffusion equations are solved using the lattice Boltzmann (LB) method. The local equilibrium distribution of the pollutant has no terms of second order in flow velocity. And the relaxation time of the pollutant deviates from a constant for the flows with variable free surface water depth. The numerical tests show that this scheme strictly obeys the conservation law of mass and momentum. Excellent agreement is obtained between numerical predictions and analytical solutions in the pure diffusion problem and convection–diffusion problem. Furthermore, the influences on the accuracy of the lattice size and the diffusivity are also studied. The results indicate that the variation in the free surface water depth cannot affect the conservation of the model, and the model has the ability to simulate the complex topography problem. The comparison shows that the LB scheme has the capacity to solve the complex convection–diffusion problem in shallow water. Copyright © 2008 John Wiley & Sons, Ltd.     

Modified shallow water equations which admit the propagation of discontinuous waves over a dry bed

A method for modeling the propagation of discontinuous waves over a dry bed using the first approximation of shallow water theory is proposed. The method is based on a modified conservation law of total momentum that takes into account the concentrated momentum losses due to the formation of local turbulent vortex structures in the fluid surface layer at a discontinuous-wave front. A quantitative estimate of these losses is obtained by deriving the shallow water equations from the Navier-Stokes equations with allowance for viscosity, which has a rapidly increasing effect in the turbulent flow regions described by discontinuous waves. The stability of the discontinuous waves admitted by the modified system of conservation laws of shallow water theory is examined. As an example, a comparative analysis is performed of the solutions of the dam-break problem obtained for the classical and modified shallow water models. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 22–43, November–December, 2007     

Shallow Water Flow Analysis with Moving Boundary Technique Using Least-squares Bubble Function

This paper presents a computational simulation method for a river problem. For the actual flow problem, it is necessary to compute flow velocity, water elevation and water region at the same time. For the basic formulation, the unsteady shallow water equations are used. As the numerical approach, implicit FEM is proposed by bubble function. To control numerical stability and accuracy, LSBF (Least-Squares Bubble Function) is used to solve the finite element equations. Also, the fixed boundary technique is combined to deal with wet and dry areas in the moving finite element mesh. Some numerical tests are shown to check this method.     

Finite volume solution to integrated shallow surface–saturated groundwater flow

This paper describes development of an integrated shallow surface and saturated groundwater model (GSHAW5). The surface flow motion is described by the 2‐D shallow water equations and groundwater movement is described by the 2‐D groundwater equations. The numerical solution of these equations is based on the finite volume method where the surface water fluxes are estimated using the Roe shock‐capturing scheme, and the groundwater fluxes are computed by application of Darcy's law. Use of a shock‐capturing scheme ensures ability to simulate steady and unsteady, continuous and discontinuous, subcritical and supercritical surface water flow conditions. Ground and surface water interaction is achieved by the introduction of source‐sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The numerical solutions and coupling algorithms are explained. The model has been applied to 1‐D and 2‐D test scenarios. The results have shown that the model can produce very accurate results and can be used for simulation of situations involving interaction between shallow surface and saturated groundwater flows. Copyright © 2005 John Wiley & Sons, Ltd.     

A mode split,Godunov‐type model for nonhydrostatic,free surface flow

A previously developed model for nonhydrostatic, free surface flow is redesigned to improve computational efficiency without sacrificing accuracy. Both models solve the Reynolds averaged Navier–Stokes equations in a fractional step manner with the pressure split into hydrostatic and nonhydrostatic components. The hydrostatic equations are first solved with an approximate Riemann solver. The hydrostatic solution is then corrected by including the nonhydrostatic pressure and requiring the velocity field to obey the incompressibility constraint. The original model requires the solution of a Riemann problem at every cell face for each vertical layer of cells, which is computationally expensive. The redesigned model instead solves the shallow water (long wave) equations for the hydrostatic solution. Vertical shear is computed by subtracting the shallow water equations from the full three dimensional equations, which removes the hydrostatic thrust terms. Therefore, the required fluxes may be more efficiently computed with velocity based upwind differencing rather than solving a Riemann problem in each vertical layer of cells. This approach is termed mode splitting and has been used in hydrostatic coastal and ocean circulation models, but not surf zone models. Numerical predictions are compared with analytical solutions and experimental data to show that the mode split model is as accurate as the original model, but requires significantly less computational effort especially for large numbers of cell layers. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.     

Lattice Boltzmann modeling of shallow water flows over discontinuous beds

Numerical modeling of shallow water flows over discontinuous beds is presented. The flows are described with the shallow water equations and the equations are solved using the lattice Boltzmann method (LBM) with single relaxation time (Bhatnagar–Gross–Krook‐LBM (BGK‐LBM)) and the multiple relaxation time (MRT‐LBM). The weighted centered scheme for force term together with the bed height for a bed slope is described to improve simulation of flows over discontinuous bed. Furthermore, the resistance stress is added to include the local head loss caused by flow over a step. Four test cases, one‐dimensional tidal over regular bed and steps, dam‐break flows, and two‐dimensional shallow water flow over a square block, are considered to verify the present method. Agreements between predictions and analytical solutions are satisfactory. Furthermore, the performance and CPU cost time of BGK‐LBM and MRT‐LBM are compared and studied. The results have shown that the lattice Boltzmann method is simple and accurate for simulating shallow water flows over discontinuous beds. This demonstrates the capability and applicability of the lattice Boltzmann method in modeling shallow water flows on bed topography with a discontinuity in practical hydraulic engineering. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparison between the theory and the numerical experiment in the problem of dam break on a jump of the cross-sectional area of a rectangular channel

Within the framework of the first approximation of the spatially one-dimensional shallow water theory the problem of flow generated by the dam break on a jump of the cross-sectional area of a rectangular channel is solved in the case in which the upper pool is wider than the lower one. The self-similar solutions constructed contain an euristic parameter related with the amount of the total flow energy lost on the cross-sectional area jump. The parameter is determined by means of comparing the one-dimensional solutions with the results of the numerical modeling of the problem on the basis of spatially two-dimensional equations of the shallow water theory.