اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA regularized shallow-water waves system with slip-wall boundary conditions in a basin: Theory and numerical analysis
63
0
0.0
(
0
)
نشر من قبل
Dimitrios Mitsotakis
تاريخ النشر
2020
مجال البحث
الهندسة المعلوماتية
والبحث باللغة
English
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
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.
قيم البحث
اقرأ أيضاً
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 wi
th 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 cla
ss 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 esti
mate 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.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
