Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA regularized shallow-water waves system with slip-wall boundary conditions in a basin: Theory and numerical analysis
63
0
0.0
(
0
)
Added by
Dimitrios Mitsotakis
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The simulation of long, nonlinear dispersive waves in bounded domains usually requires the use of slip-wall boundary conditions. Boussinesq systems appearing in the literature are generally not well-posed when such boundary conditions are imposed, or if they are well-posed it is very cumbersome to implement the boundary conditions in numerical approximations. In the present paper a new Boussinesq system is proposed for the study of long waves of small amplitude in a basin when slip-wall boundary conditions are required. The new system is derived using asymptotic techniques under the assumption of small bathymetric variations, and a mathematical proof of well-posedness for the new system is developed. The new system is also solved numerically using a Galerkin finite-element method, where the boundary conditions are imposed with the help of Nitsches method. Convergence of the numerical method is analyzed, and precise error estimates are provided. The method is then implemented, and the convergence is verified using numerical experiments. Numerical simulations for solitary waves shoaling on a plane slope are also presented. The results are compared to experimental data, and excellent agreement is found.
rate research
Read More
We propose a new iterative scheme to compute the numerical solution to an over-determined boundary value problem for a general quasilinear elliptic PDE. The main idea is to repeatedly solve its linearization by using the quasi-reversibility method with a suitable Carleman weight function. The presence of the Carleman weight function allows us to employ a Carleman estimate to prove the convergence of the sequence generated by the iterative scheme above to the desired solution. The convergence of the iteration is fast at an exponential rate without the need of an initial good guess. We apply this method to compute solutions to some general quasilinear elliptic equations and a large class of first-order Hamilton-Jacobi equations. Numerical results are presented.
We consider approximating the solution of the Helmholtz exterior Dirichlet problem for a nontrapping obstacle, with boundary data coming from plane-wave incidence, by the solution of the corresponding boundary value problem where the exterior domain is truncated and a local absorbing boundary condition coming from a Pade approximation (of arbitrary order) of the Dirichlet-to-Neumann map is imposed on the artificial boundary (recall that the simplest such boundary condition is the impedance boundary condition). We prove upper- and lower-bounds on the relative error incurred by this approximation, both in the whole domain and in a fixed neighbourhood of the obstacle (i.e. away from the artificial boundary). Our bounds are valid for arbitrarily-high frequency, with the artificial boundary fixed, and show that the relative error is bounded away from zero, independent of the frequency, and regardless of the geometry of the artificial boundary.
The paper proposes a new, conservative fully-discrete scheme for the numerical solution of the regularised shallow water Boussinesq system of equations in the cases of periodic and reflective boundary conditions. The particular system is one of a class of equations derived recently and can be used in practical simulations to describe the propagation of weakly nonlinear and weakly dispersive long water waves, such as tsunamis. Studies of small-amplitude long waves usually require long-time simulations in order to investigate scenarios such as the overtaking collision of two solitary waves or the propagation of transoceanic tsunamis. For long-time simulations of non-dissipative waves such as solitary waves, the preservation of the total energy by the numerical method can be crucial in the quality of the approximation. The new conservative fully-discrete method consists of a Galerkin finite element method for spatial semidiscretisation and an explicit relaxation Runge--Kutta scheme for integration in time. The Galerkin method is expressed and implemented in the framework of mixed finite element methods. The paper provides an extended experimental study of the accuracy and convergence properties of the new numerical method. The experiments reveal a new convergence pattern compared to standard Galerkin methods.
Energy estimates of the shallow water equations (SWEs) with a transmission boundary condition are studied theoretically and numerically. In the theoretical part, using a suitable energy, we begin with deriving an equality which implies an energy estimate of the SWEs with the Dirichlet and the slip boundary conditions. For the SWEs with a transmission boundary condition, an inequality for the energy estimate is proved under some assumptions to be satisfied in practical computation. Hence, it is recognized that the transmission boundary condition is reasonable in the sense that the inequality holds true. In the numerical part, based on the theoretical results, the energy estimate of the SWEs with a transmission boundary condition is confirmed numerically by a finite difference method (FDM). The choice of a positive constant c0 used in the transmission boundary condition is investigated additionally. Furthermore, we present numerical results by a Lagrange-Galerkin scheme, which are similar to those by the FDM. From the numerical results, it is found that the transmission boundary condition works well numerically.
We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poissons and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poissons problem and the linear peridynamic solid model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
