Unstructured finite volume discretization of two‐dimensional depth‐averaged shallow water equations with porosity

This paper deals with the numerical discretization of two‐dimensional depth‐averaged models with porosity. The equations solved by these models are similar to the classic shallow water equations, but include additional terms to account for the effect of small‐scale impervious obstructions which are not resolved by the numerical mesh because their size is smaller or similar to the average mesh size. These small‐scale obstructions diminish the available storage volume on a given region, reduce the effective cross section for the water to flow, and increase the head losses due to additional drag forces and turbulence. In shallow water models with porosity these effects are modelled introducing an effective porosity parameter in the mass and momentum conservation equations, and including an additional drag source term in the momentum equations. This paper presents and compares two different numerical discretizations for the two‐dimensional shallow water equations with porosity, both of them are high‐order schemes. The numerical schemes proposed are well‐balanced, in the sense that they preserve naturally the exact hydrostatic solution without the need of high‐order corrections in the source terms. At the same time they are able to deal accurately with regions of zero porosity, where the water cannot flow. Several numerical test cases are used in order to verify the properties of the discretization schemes proposed. Copyright © 2009 John Wiley & Sons, Ltd.     

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

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

Analysis of efficient preconditioned defect correction methods for nonlinear water waves

Robust computational procedures for the solution of non‐hydrostatic, free surface, irrotational and inviscid free‐surface water waves in three space dimensions can be based on iterative preconditioned defect correction (PDC) methods. Such methods can be made efficient and scalable to enable prediction of free‐surface wave transformation and accurate wave kinematics in both deep and shallow waters in large marine areas or for predicting the outcome of experiments in large numerical wave tanks. We revisit the classical governing equations are fully nonlinear and dispersive potential flow equations. We present new detailed fundamental analysis using finite‐amplitude wave solutions for iterative solvers. We demonstrate that the PDC method in combination with a high‐order discretization method enables efficient and scalable solution of the linear system of equations arising in potential flow models. Our study is particularly relevant for fast and efficient simulation of non‐breaking fully nonlinear water waves over varying bottom topography that may be limited by computational resources or requirements. To gain insight into algorithmic properties and proper choices of discretization parameters for different PDC strategies, we study systematically limits of accuracy, convergence rate, algorithmic and numerical efficiency and scalability of the most efficient known PDC methods. These strategies are of interest, because they enable generalization of geometric multigrid methods to high‐order accurate discretizations and enable significant improvement in numerical efficiency while incuring minimal storage requirements. We demonstrate robustness using such PDC methods for practical ranges of interest for coastal and maritime engineering, that is, from shallow to deep water, and report details of numerical experiments that can be used for benchmarking purposes. Copyright © 2014 John Wiley & Sons, Ltd.     

M-lump,rogue waves,breather waves,and interaction solutions among four nonlinear waves to new (3+1)-dimensional Hirota bilinear equation

In this research article, we study a new (3+1)-dimensional Hirota bilinear equation which can describe the dynamics of ion-acoustic wave and Alvin wave of small but finite amplitude in plasma physics and describe the propagation process of nonlinear waves in shallow water. First, we apply two methods to study the equation, namely the Hirota bilinear method and long-wave limit method M-lump solution, and line rogue waves are reported. Furthermore, we investigate the velocity, propagation trajectory, and interaction phenomenon of M-lump solution(M=2,3). Then, based on the multi-solitons, two cases of high-order breather solution are constructed by selecting some special parameters. Finally, four types interaction solutions are successfully obtained by employing long-wave limit method and selecting some special parameters. More importantly, we explore physical collision phenomenon of the interaction between nonlinear waves. In order to better illustrate the characteristics of the interaction solutions, the results are shown in three-dimensional plots and numerical simulation. To our knowledge, all of the obtained solutions in this article are novel. The results of this article may be provide an important theoretical basis for explaining some nonlinear phenomena in the field of fluid mechanics and shallow water.

     

Traveling Wave Solutions for a Class of One-Dimensional Nonlinear Shallow Water Wave Models

In this paper we consider a class of one-dimensional nonlinear shallow water wave models that support weak solutions. We construct new traveling wave solutions for these models. Moreover, we show that these new traveling wave solutions are stable.     

Numerical boundary conditions for globally mass conservative methods to solve the shallow‐water equations and applied to river flow

A revision of some well‐known discretization techniques for the numerical boundary conditions in 1D shallow‐water flow models is presented. More recent options are also considered in the search for a fully conservative technique that is able to preserve the good properties of a conservative scheme used for the interior points. Two conservative numerical schemes are used as representatives of the families of explicit and implicit numerical methods. The implementation of the different boundary options to these schemes is compared by means of the simulation of several test cases with exact solution. The schemes with the global conservation boundary discretization are applied to the simulation of a real river flood wave leading to very satisfactory results. Copyright © 2005 John Wiley & Sons, Ltd.     

Nonlinear long-term behaviour of spherical shallow thin-walled concrete shells of revolution

The nonlinear long-term buckling behaviour (creep buckling) of spherical shallow, thin-walled concrete shells of revolution (including domes) subjected to sustained loads is investigated herein. A thorough understanding of their nonlinear time-dependent behaviour, as well as the development of comprehensive analytical models for their analysis, has hitherto not been fully established and further studies are required. A nonlinear axisymmetric theoretical model, which accounts for the effects of creep and shrinkage, and which considers the ageing of the concrete material and the variation of the internal stresses and geometry in time, is developed for this purpose. The governing field equations are derived using variational principles, equilibrium requirements, and integral-type constitutive relations. A systematic step-by-step procedure is used for the solution of the integral-type governing equations. First, the nonlinear short-term behaviour is studied to provide a benchmark for the long-term analysis. Different theories for the analysis of the shell structure are examined for this purpose and compared with results obtained by the finite element method. A numerical study, which highlights the capabilities of the nonlinear theoretical model and which provides insight into the nonlinear long-term behaviour of shallow concrete domes, is presented. The results show that long-term effects are critical for the design and structural safety of shallow, thin-walled concrete domes, and so these effects need to be fully understood and quantifiable.     

Comparing analytical and numerical solution of nonlinear two and three‐dimensional hydrostatic flows

New test cases for frictionless, three‐dimensional hydrostatic flows have been derived from some known analytical solutions of the two‐dimensional shallow water equations. The flow domain is a paraboloid of revolution and the flow is determined by the initial conditions, the nonlinear advective terms, the Coriolis acceleration and by the hydrostatic pressure. Wetting and drying is also included. Some specific properties of the exact solutions are discussed under different hypothesis and relative importance of the forcing terms. These solutions are proposed for testing the stability, the accuracy and the efficiency of numerical models to be used for simulating environmental hydrostatic flows. The computed solutions obtained with a semi‐implicit finite difference—finite volume algorithm on unstructured grid are compared with the corresponding analytical solutions in both two and three space dimension. Excellent agreement are obtained for the velocity and for the resulting water surface elevation. Comparison of the computed inundation area also shows a good agreement with the analytical solution with degrading accuracy observed when the inundation area becomes relatively large and for long simulation time. Copyright © 2006 John Wiley & Sons, Ltd.     

浅水区海底管道周围海床渗透压力计算研究

我国海上油田开采起步较晚,大部分油田处于浅水区,因此,在设计管道时,应充分考虑由浅水区波浪引起的管道周围海床渗流力。根据浅水波相关假设,考虑自由水面非线性影响,推导出椭圆余弦波的波面方程,在此基础上进一步得到一个关于速度势的表达式,并根据该表达式得出作用于海床表面的波压公式。考虑海床土的压缩性,推导出一阶近似椭圆余弦波作用下浅水区埋置管道周围海床的渗流压力解析解,最后将计算结果与大型水槽试验及以往研究成果作对比。结果表明,在椭圆余弦波的作用下,由一阶椭圆余弦波理论得到的计算结果与试验结果规律基本一致,与相似工况下的现有理论成果数值基本相同,具有一定的可行性和工程价值。     

A modified Galerkin/finite element method for the numerical solution of the Serre‐Green‐Naghdi system

A new modified Galerkin/finite element method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low‐order Lagrange finite element spaces, despite the fact that the system contains third order spatial partial derivatives for the depth averaged velocity of the fluid. After studying the efficacy and the conservation properties of the new numerical method, we proceed with the validation of the new numerical model and boundary conditions by comparing the numerical solutions with laboratory experiments and with available theoretical asymptotic results. Copyright © 2016 John Wiley & Sons, Ltd.     

A meshless numerical wave tank for simulation of nonlinear irregular waves in shallow water

Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd.     

Linear and geometrically nonlinear analysis of non-uniform shallow arches under a central concentrated force

In this paper an integral equation solution to the linear and geometrically nonlinear problem of non-uniform in-plane shallow arches under a central concentrated force is presented. Arches exhibit advantageous behavior over straight beams due to their curvature which increases the overall stiffness of the structure. They can span large areas by resolving forces into mainly compressive stresses and, in turn confining tensile stresses to acceptable limits. Most arches are designed to operate linearly under service loads. However, their slenderness nature makes them susceptible to large deformations especially when the external loads increase beyond the service point. Loss of stability may occur, known also as snap-through buckling, with catastrophic consequences for the structure. Linear analysis cannot predict this type of instability and a geometrically nonlinear analysis is needed to describe efficiently the response of the arch. The aim of this work is to cope with the linear and geometrically nonlinear problem of non-uniform shallow arches under a central concentrated force. The governing equations of the problem are comprised of two nonlinear coupled partial differential equations in terms of the axial (tangential) and transverse (normal) displacements. Moreover, as the cross-sectional properties of the arch vary along its axis, the resulting coupled differential equations have variable coefficients and are solved using a robust integral equation numerical method in conjunction with the arc-length method. The latter method allows following the nonlinear equilibrium path and overcoming bifurcation and limit (turning) points, which usually appear in the nonlinear response of curved structures like shallow arches and shells. Several arches are analyzed not only to validate our proposed model, but also to investigate the nonlinear response of in-plane thin shallow arches.     

Quelques modèles diffusifs capillaires de type KortewegSome diffusive capillary models of Korteweg type.

The purpose of this Note is to propose new diffusive capillary models of Korteweg type and discuss their mathematical properties. More precisely, we introduce viscous models which provide some additional information on the behavior of the density close to vacuum. We actually prove that if some compatibility conditions between diffusion and capillarity are satisfied, some extra regularity information on a quantity involving the density is available. We obtain a non-trivial equality deduced from the special structure of the momentum equation. This Note generalizes to some extent the authors' previous works on the Korteweg model (with constant capillary coefficient) and on the shallow water equation. To cite this article: D. Bresch, B. Desjardins, C. R. Mecanique 332 (2004).     

Nonlinear dynamic behaviors of viscoelastic shallow arches

The nonlinear dynamic behaviors of nonlinear viscoelastic shallow arches sub- jected to external excitation are investigated. Based on the d＇Alembert principle and the Euler-Bernoulli assumption, the governing equation of a shallow arch is obtained, where the Leaderman constitutive relation is applied. The Galerkin method and numerical in- tegration are used to study the nonlinear dynamic properties of the viscoelastic shallow arches. Moreover, the effects of the rise, the material parameter and excitation on the nonlinear dynamic behaviors of the shallow arch viscoelastic shallow arches may appear to have a are investigated. The results show that chaotic motion for certain conditions.     

Fluid surface oscillations in a shallow two-dimensional basin

The two-dimensional oscillation of a free liquid surface in a shallow basin of arbitrary cross-section is considered It is assumed that the fluid motion is not influenced by the earth rotation and that the fluid is ideal, Since the linearized equations are not valid in the region of the shoreline, they are used only in the center region of the basin. In the proximity of the shoreline the method of characteristics is used as applied to the nonlinear shallow water equations. The bottom of the basin in this region is approximated by a straight line. Since it is implicitly assumed that the center region drives the shallow wedge at the shoreline, but conversely the shore region does not affect the center region, the usefulness of the solution is confined to the first few cycles only. The analytical solution was checked by applying a numerical method of solution for the entire basin. The agreement was satisfactory over the first two cycles.     

A high‐order PIC method for advection‐dominated flow with application to shallow water waves

A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.     

Nonlinear static and dynamic analysis for laminated,annular, spherical caps of moderate thickness

Based on Timoshenko-Mindlin kinematic hypothesis, the shallow shell theory is extended to include the transverse shear deformation for the nonlinear axisymmetric dynamic analysis of the symmetric cross-ply shallow spherical shell. Using the orthogonal point collocation method and the Newmark scheme, an iterative solution is formulated. The numerical results for the nonlinear static and dynamic responses and dynamic buckling of these shallow spherical shells with circular holes under uniformly distributed static or dynamic normal impact loads are presented and compared with available data.     

Nonlinear Multigrid Methods for Numerical Solution of the Unsaturated Flow Equation in Two Space Dimensions

Picard and Newton iterations are widely used to solve numerically the nonlinear Richards’ equation (RE) governing water flow in unsaturated porous media. When solving RE in two space dimensions, direct methods applied to the linearized problem in the Newton/Picard iterations are inefficient. The numerical solving of RE in 2D with a nonlinear multigrid (MG) method that avoids Picard/Newton iterations is the focus of this work. The numerical approach is based on an implicit, second-order accurate time discretization combined with a second-order accurate finite difference spatial discretization. The test problems simulate infiltration of water in 2D unsaturated soils with hydraulic properties described by Broadbridge–White and van Genuchten–Mualem models. The numerical results show that nonlinear MG deserves to be taken into consideration for numerical solving of RE.     

破碎带波浪的数值模拟

基于一组色散关系得到改进的完全非线性Boussinesq方程建立了一个波浪模型可以模拟近岸水域的波浪变浅、破碎以及在海滩上的爬高等多种变形。波浪破碎引起的能量衰减是在动量方程中引入一个在空间和时间上都只作用于波前的涡粘项来模拟。动海岸线边界用窄缝法处理。波浪爬高用非线性浅水方程推导的非破碎波浪在斜坡上爬高的解析解来验证。本模型还模拟了波浪在斜坡上不同类型的破碎变形过程,并将其波高和平均水位的沿程变化和物理模型实验的结果比较,两者符合良好。