Numerical simulation of oblique and multidirectional wave propagation and breaking on steep slope based on FEM model of Boussinesq equations

Creating a representative numerical simulation of the propagation and breaking of waves along slopes is an important problem in engineering design. Most studies on wave breaking have focused on the propagation of normal incident waves on gentle slopes. In practice, however, waves on steep slopes are obliquely incident or multidirectional irregular waves. In this paper, the eddy viscosity term is introduced to the momentum equation of the improved Boussinesq equations to model wave dissipation caused by breaking and friction, and a numerical model based on an unstructured finite element method (FEM) is established based on the governing equations. It is applied to simulate wave propagation on a steep slope of 1:5. Parallel physical experiments are conducted for comparative analysis that considered a large number of cases, including those featuring of normal and oblique incident regular and irregular waves, and multidirectional waves. The heights of the incident wave increase for different periods to represent different kinds of waves breaking. Based on examination, the effectiveness and accuracy of the numerical model is verified through a comprehensive comparison between the numerical and the experimental results, including in terms of variation in wave height, wave spectrum, and nonlinear parameters. Satisfactory agreement between the numerical and experimental values shows that the proposed model is effective in representing the breaking of oblique incident regular waves, irregular waves, and multidirectional incident irregular waves. However, the initial threshold of the breaking parameter η_t^(I) takes different values for oblique and multidirectional waves. This needs to be paid attention when the breaking of waves is simulated using the Boussinesq equations.     

Generating Boundary Conditions for the Calculation of Tsunami Propagation on Nested Grids

Some boundary conditions used to numerically simulate tsunami generation and propagation are studied. Special attention is given to generating boundary conditions thatmake it possible to simulate tsunami waves with desired characteristics (amplitude, time period and, in general, waveform). Since the water flow velocity in a propagating tsunami wave is uniquely defined by its height and ocean depth, one can simulate a wave propagating from the boundary into the simulation area. This can be done by specifying the wave height and water flow velocity on the boundary. This method is used to numerically simulate the propagation of a tsunami from the source to the coast on a sequence of refined grids. In this numerical experiment the wave parameters are transferred from the larger area to the subarea via boundary conditions. This method can also generate a wave that has certain characteristics on a specified line.     

An unstructured FEM model based on Boussinesq equations and its application to the calculation of multidirectional wave run-up in a cylinder group

In this paper, a numerical model based on the improved Boussinesq equations derived by Beji and Nadaoka [5] is first developed using unstructured finite element technique. A locally rotated coordinate system is introduced to improve the treatment for the fully reflective boundaries whose orientation does not coincide with the coordinate system. The Adams–Bashforth–Moulton predictor–corrector scheme is used for time integration. Typical examples are employed to validate the numerical model. Based on the developed model, multidirectional wave propagation through a cylinder group is numerically calculated and the effects of the wave directionality on the waves in the group and the wave run-up on the cylinders are investigated. Numerical results show that the wave directionality has considerable effect on the wave run-up in the cylinder group.     

多偏移距电磁波数值成像的高阶FDTD法研究

数值模拟是数据处理和反演解释的重要环节,为提高高频电磁波层析成像技术在工程勘察中的解译准确性,采用高阶时间域有限差分法(FDTD(2阶时间精度,4阶空间精度))模拟电磁波在溶洞地电模型下的数值成像,并采用多偏移距模式采集电磁波道集.分析了电磁波的波场快照图和波形剖面图,总结了电磁波在溶洞地质体中的传播规律,为实际工程勘察工作提供可靠依据.数值试验结果表明,FDTD(2,4)能高精度模拟电磁波在溶洞复杂介质中的数值成像和多偏移距采集模式具有高效性和灵活性.     

Wave generation in a computation domain

One of the possible methods to deal with the reflecting waves at the incident boundary in numerical modeling is to generate waves in the computation domain and absorb the outgoing waves at the incident boundary. A source function is introduced into the momentum equation of Boussinesq equations for generating wave in a computation domain in this paper. Typical numerical examples are given for the verification of the proposed method. Numerical examination for the wave diffraction through a breakwater gap shows that the proposed method is especially useful for multidirectional waves.     

一维高精度离散GDQ方法

GDQ method is a kind of high order accurate numerical methods developed several years ago, which have been successfully used to simulate the solution of smooth engineering problems such as structure mechanics and incompressible fluid dynamics. In this paper, extending the traditional GDQ method, we develop a new kind of discontinuous GDQ methods to solve compressible flow problems of which solutions may be discontinuous. In order to capture the local features of fluid flows, firstly, the computational domain is divided into many small pieces of subdomains. Then, in each small subdomain, the GDQ method is implementedand some kinds of numerical flux limitation conditions will be required to keep the correct flow direction. At the boundary interface between subdomains, we also use some kind of flux conditions according to the flow direction. The numerical method obtained by the above steps has the advantages of high order accuracy and easy to treat boundary conditions. It can simulate perfectly nonlinear waves such as shock, rarefaction wave and contact discontinuity. Finally, the numerical experiments on one dimensional Burgers equation and Euler equations are given.The numerical results verify the validation of the method.     

用CT-VOF法在数值波试验槽中生成线性和非线性波

为渡水槽中波的模拟和传播提出了二维的数值模型．假设流动的流体为粘性、不可压缩的,并将Navier-Stokes方程和连续性方程作为控制方程．用标准的k-ε模型来模拟紊流流动;用交错网格的有限差分法,离散化Navier-Stokes方程;并用简化的标识和单元(SMAC)方法进行求解．使用活塞型波发生器生成并传播波;数值渡水槽的端部采用敞开式的边界条件．为了证明模型的有效性,进行了一些标准的试验,如顶盖驱动的方腔测试试验、单向的常速度场试验以及干燥河床上的溃坝试验．为了论证方法的性能及其精度,将所生成波的结果与已有波理论的结果进行比较．最后,采用群集技术(CT)生成网格,并提出最佳的网格生成条件．     

A Kind of Discrete Non-Reflecting Boundary Conditions for Varieties of Wave Equations

Abstract In this paper, a new kind of discrete non-reflecting boundary conditions is developed.It can be usedfor a variety of wave equations such as the acoustic wave equation, the isotropic and anisotropic elastic waveequations and the equations for wave propagation in multi-phase media and so on.In this kind of boundaryconditions,the composition of all artifical reflected waves,but not the individual reflected ones,is consideredand eliminated.Thus, it has a uniform formula for different wave equations.The velocity C_A of the composedreflected wave is determined in the way to make the reflection coefficients minimal,the value of which depends onequations.In this psper,the construction of the boundary conditions illustrated and C_A is found,numericalresults are presented to illustrate the effectiveness of the boundary conditions.     

Long distance wave computation using nonlinear solitary waves

A recently developed method is described to propagate short wave equation pulses over indefinite distances and through regions of varying indices of refraction, including multiple reflections. The method, “Wave Confinement”, utilizes a newly developed nonlinear partial differential equation (pde) that propagates basis functions according to the wave equation. These basis functions are generated as stable solitary waves where the discretized equation can be solved without any numerical dissipation. The method can also be used to solve for harmonic waves in the high frequency (Eikonal) limit, including multiple arrivals. The solution involves discretizing the wave equation on a uniform Eulerian grid and adding a simple nonlinear “Confinement” term. This term does not change the amplitude (integrated through each point on the pulse surface) or the propagation velocity, or arrival time, and yet results in capturing the waves as thin surfaces that propagate as thin nonlinear solitary waves and remain ∼2-3 grid cells in thickness indefinitely with no numerical spreading. A new feature described in this paper involves computing scattering of short pulses from complex objects such as complete aircraft. A simple “immersed surface” approach is used, that utilizes the same uniform grid as the propagation and avoids complex, body fitted or adaptive grid schemes.The new method should be useful in areas of wave propagation, from radar scattering and long distance communications to cell phone transmission.     

Finite element modelling of surface waves in a channel

A variational formulation of the vertically-integrated differential equations for free surface wave motion is presented. A finite element model is derived for solving this nonlinear system of hydrodynamic equations. The time integration scheme employed is discussed and the results obtained demonstrate its good stability and accuracy.Several applications of the model are considered: the first problem is an open channel of uniform depth and the second an open channel of linearly varying depth. The ‘inflow’ boundary condition is prescribed in terms of the velocity which represents a wavemaker and/or a flow source, while the ‘outflow’ boundary condition is specified in terms of the water elevation. The outflow condition is adjusted for two cases, a reflecting boundary (finite channel) and a non-reflecting boundary (open-ended channel). The latter boundary condition is examined in some detail and the results obtained show that the numerical model can produce the non-reflecting boundary that is similar to the analytical radiation condition for waves. Computational results for a third problem, involving wave reflection from a submerged cylinder, are also presented and compared with both experimental data and analytical predictions.The simplicity and the performance of the computational model suggest that free surface waves can be simulated without excessively complicated numerical schemes. The ability of the model to simulate outflow boundary conditions properly is of special importance since these conditions present serious problems for many numerical algorithms.     

Well-Posedness of a Problem of Propagation of Nonlinear Acoustic-Gravity Waves in the Atmosphere Generated by Surface Pressure Variations

Currently there are many international microbarograph networks for high-resolution recording of wave pressure variations on the Earth’s surface. This arouses interest in wave propagation in the atmosphere generated by atmospheric pressure variations. A full system of nonlinear hydrodynamic equations for atmospheric gases with lower boundary conditions in the form of wavelike pressure variations on the Earth’s surface is considered. Since the wave amplitudes near the Earth’s surface are small, linearized equations are used in the analysis of well-posedness of the problem. With the help of a wave energy functional method, it is shown that in the non-dissipative case the solution to the boundary value problem is uniquely determined by the variable pressure field on the Earth’s surface. The corresponding dissipative problem is well-posed if, in addition to the pressure field, appropriate conditions on the velocity and temperature on the Earth’s surface are given. In the case of an isothermal atmosphere, the problem admits analytical solutions that are harmonic in the variables x and t. A good agreement between the numerical and analytical solutions is obtained. The study shows that the temperature and density can rapidly vary at the lower boundary of the boundary value problem. An example of solving the three-dimensional problem with variable pressure on the Earth’s surface taken from experimental observations is given. The developed algorithms and computer programs can be used to simulate atmospheric waves generated by pressure variations on the Earth’s surface.     

Nonlinear Wave–Current Interactions in Shallow Water

We study here the propagation of long waves in the presence of vorticity. In the irrotational framework, the Green–Naghdi equations (also called Serre or fully nonlinear Boussinesq equations) are the standard model for the propagation of such waves. These equations couple the surface elevation to the vertically averaged horizontal velocity and are therefore independent of the vertical variable. In the presence of vorticity, the dependence on the vertical variable cannot be removed from the vorticity equation but it was however shown in 1 that the motion of the waves could be described using an extended Green–Naghdi system. In this paper, we propose an analysis of these equations, and show that they can be used to get some new insight into wave–current interactions. We show in particular that solitary waves may have a drastically different behavior in the presence of vorticity and show the existence of solitary waves of maximal amplitude with a peak at their crest, whose angle depends on the vorticity. We also show some simple numerical validations. Finally, we give some examples of wave–current interactions with a nontrivial vorticity field and topography effects.     

Analysis of Guided Lamb Wave Propagation (GW) in Honeycomb Sandwich Panels

The paper aims to introduce the guided lamb wave propagation (GW) in a honeycomb sandwich panels to be used in the health monitoring applications. Honeycomb sandwich panels are well-known as lightweight structures with a good stiffness behavior and a wide range of applications in different industries. Due to the complex geometry and complicated boundary conditions in such a structure, the development of analytical solutions for describing the wave propagation and the interaction of waves with damages is hardly possible. Therefore dimensional finite element simulations have been used to model GW for different frequency ranges and different sandwich panels with different geometrical properties. The waves, which are highly dispersive, have been excited by thin piezoelectric patches attached to the surface of the structure. In the first step, the honeycomb panel has been simplified as an orthotropic layered continuum medium. The required material data have been calculated by applying a numerical homogenization method for the honeycomb core layer. The wave propagation has been compared in the homogenized model with the real geometry of a honeycomb sandwich panel. Such calculations of high frequency ultrasonic waves are costly, both in creating a proper finite element model as well as in the required calculation time. In this paper the influence of changes in the geometry of the sandwich panel on the wave propagation is presented. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A simple and robust boundary treatment for the forced Korteweg–de Vries equation

In this paper, we propose a simple and robust numerical method for the forced Korteweg–de Vries (fKdV) equation which models free surface waves of an incompressible and inviscid fluid flow over a bump. The fKdV equation is defined in an infinite domain. However, to solve the equation numerically we must truncate the infinite domain to a bounded domain by introducing an artificial boundary and imposing boundary conditions there. Due to unsuitable artificial boundary conditions, most wave propagation problems have numerical difficulties (e.g., the truncated computational domain must be large enough or the numerical simulation must be terminated before the wave approaches the artificial boundary for the quality of the numerical solution). To solve this boundary problem, we develop an absorbing non-reflecting boundary treatment which uses outward wave velocity. The basic idea of the proposing algorithm is that we first calculate an outward wave velocity from the solutions at the previous and present time steps and then we obtain a solution at the next time step on the artificial boundary by moving the solution at the present time step with the velocity. And then we update solutions at the next time step inside the domain using the calculated solution on the artificial boundary. Numerical experiments with various initial conditions for the KdV and fKdV equations are presented to illustrate the accuracy and efficiency of our method.     

Numerical Study of Nonlinear Scattering Characteristics of SH0 Waves Encountering Cracks in Prestressed Plates北大核心CSCD

超声导波因具有传播距离远、能量衰减小等优点在结构健康监测领域中被广泛关注.厘清结构中导波与损伤作用后的散射规律,对于传感器阵列的设计和信号分析均具有重要意义.通过发展的数值方法,研究了受载结构中零阶水平剪切波(SH 0波)与微裂纹作用的接触声非线性作用规律.在双势谱方法的基础上,进一步通过mortar方法将谱单元和有限单元进行了耦合,以充分利用谱元法计算导波传播效率高的优点和有限元在离散复杂结构中的优势.利用该方法计算了板壳结构在自由状态和受载状态下SH 0波与不同角度微裂纹作用的非线性散射场.结果表明,SH 0波与裂纹作用后的二次谐波散射场关于裂纹面近似对称分布,并且单轴预应力不会改变二次谐波散射场的对称性,仍可以通过散射场的分布来确定微裂纹的取向.     

Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves. In this paper, we will study two implicit finite difference schemes for the simulation of waves. They are the weighted alternating direction implicit (ADI) scheme and the locally one-dimensional (LOD) scheme. The approximation errors, stability conditions, and dispersion relations for both schemes are investigated. Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme. Moreover, the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time. In order to improve computational efficiency, numerical algorithms based on message passing interface (MPI) are implemented. Numerical examples of wave propagation in a three-layer model and a standard complex model are presented. Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.     

Singular boundary method based on time‐dependent fundamental solutions for active noise control

Active noise control is an efficient strategy of noise control. A numerical wave shielding model to inhibit wave propagation, which can be considered as an extension of traditional active noise control, is established using the singular boundary method using time‐dependent fundamental solutions in this study. Two empirical formulas to evaluate the origin intensity factors with Dirichlet and Neumann boundary conditions are derived respectively. In comparison with other similar numerical methods, the method can obtain highly accurate results using very few boundary nodes and small CPU time. These meet the major technical requirements of simulation of active noise control. The subsequent numerical experiments show that the proposed model can shield efficiently from the wave propagation for both inner and exterior problems. By applying the newly derived empirical formulas, the CPU time of the singular boundary method is further reduced significantly, which makes the method a competitive new and efficient meshless method. In addition, the singular boundary method makes active noise control in an online manner via time‐dependent fundamental solutions as its basis functions.     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.     

Boundary conditions for wave propagation problems

Wave propagation simulation requires a correct implementation of boundary conditions to avoid numerical instabilities. Similar problems are posed by domain decomposition methods where the aim is to find the correct modeling of physical phenomena across the interfaces separating the subdomains. The technique described here is based on physical grounds since it relies on the fact that the wave equation can be decomposed into incoming and outgoing wave modes at the boundary. The result is a modified wave equation for the boundaries which automatically includes the boundary condition. The boundary treatment is applied to a realistic problem of ultrasonic wave propagation through a vertical interface separating an anelastic solid at the surface. The results show that the method correctly describes the anelastic properties of the Rayleigh wave in the presence of a strong contrast in the material properties.     

Analysis of the Nonlinear Shallow Water Equations Over Nonplanar Topography

The role played by the beach bottom profile on coastal inundation phenomena is analyzed here by means of approximate analytical solutions of the nonlinear shallow water equations (NSWEs) over uneven bottoms. These are obtained by only using the assumptions of small waves at the seaward boundary and small topographic forcing. Our work, built on the Carrier and Greenspan [ 1 ] hodographic transformation and on the solution of the boundary value problem (BVP) for the NSWEs proposed by Antuono and Brocchini [ 2 ], focuses on the propagation of nonlinear non-breaking waves over quasi-planar beaches. Since the terms associated with the perturbed bottom only appear in the second-order perturbed solutions, the breaking conditions for the planar-beach bathymetry also predict well the breaking occurring on the nonplanar beaches analyzed here. The most important results, concerning the shoreline position and the near-shoreline velocity, are given for both pulse-like and periodic input waves propagating over two types of nonplanar bathymetries. The solution proposed here is a fundamental benchmark for any numerical and theoretical analyzes concerned with estimates of wave run-up on beaches of complex shape.