Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMathematical analysis of wave radiation by a step-like surface
118
0
0.0
(
0
)
Added by
Wangtao Lu
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Authors
Wangtao Lu
Ask ChatGPT about the research
No Arabic abstract
This paper proposes, for wave propagating in a globally perturbed half plane with a perfectly conducting step-like surface, a sharp Sommerfeld radiation condition (SRC) for the first time, an analytic formula of the far-field pattern, and a high-accuracy numerical solver. We adopt the Wiener-Hopf method to compute the Green function for a cracked half plane, a background for the perturbed half plane. We rigorously show that the Green function asymptotically satisfies a universal-direction SRC (uSRC) and radiates purely outgoing at infinity. This helps to propose an implicit transparent boundary condition for the scattered wave, by either a cylindrical incident wave due to a line source or a plane incident wave. Then, a well-posedness theory is established via an associated variational formulation. The theory reveals that the scattered wave, post-subtracting a known wave field, satisfies the same uSRC so that its far-field pattern is accessible theoretically. For a plane-wave incidence, asymptotic analysis shows that merely subtracting reflected plane waves, due to non-uniform heights of the step-like surface at infinity, from the scattered wave in respective regions produces a discontinuous wave satisfying the uSRC as well. Numerically, we adopt a previously developed perfectly-matched-layer (PML) boundary-integral-equation method to solve the problem. Numerical results demonstrate that the PML truncation error decays exponentially fast as thickness or absorbing power of the PML increases, of which the convergence relies heavily on the Green function exponentially decaying in the PML.
rate research
Read More
By improving the trace finite element method, we developed another higher-order trace finite element method by integrating on the surface with exact geometry description. This method restricts the finite element space on the volume mesh to the surface accurately, and approximates Laplace-Beltrami operator on the surface by calculating the high-order numerical integration on the exact surface directly. We employ this method to calculate the Laplace-Beltrami equation and the Laplace-Beltrami eigenvalue problem. Numerical error analysis shows that this method has an optimal convergence order in both problems. Numerical experiments verify the correctness of the theoretical analysis. The algorithm is more accurate and easier to implement than the existing high-order trace finite element method.
This paper studies the PML method for wave scattering in a half space of homogeneous medium bounded by a two-dimensional, perfectly conducting, and locally defected periodic surface, and develops a high-accuracy boundary-integral-equation (BIE) solver. Along the vertical direction, we place a PML to truncate the unbounded domain onto a strip and prove that the PML solution converges linearly to the true solution in the physical subregion of the strip with the PML thickness. Laterally, we divide the unbounded strip into three regions: a region containing the defect and two semi-waveguide regions, separated by two vertical line segments. In both semi-waveguides, we prove the well-posedness of an associated scattering problem so as to well define a Neumann-to-Dirichlet (NtD) operator on the associated vertical segment. The two NtD operators, serving as exact lateral boundary conditions, reformulate the unbounded strip problem as a boundary value problem onto the defected region. Due to the periodicity of the semi-waveguides, both NtD operators turn out to be closely related to a Neumann-marching operator, governed by a nonlinear Riccati equation. It is proved that the Neumann-marching operators are contracting, so that the PML solution decays exponentially fast along both lateral directions. The consequences culminate in two opposite aspects. Negatively, the PML solution cannot exponentially converge to the true solution in the whole physical region of the strip. Positively, from a numerical perspective, the Riccati equations can now be efficiently solved by a recursive doubling procedure and a high-accuracy PML-based BIE method so that the boundary value problem on the defected region can be solved efficiently and accurately. Numerical experiments demonstrate that the PML solution converges exponentially fast to the true solution in any compact subdomain of the strip.
In this work, by introducing the seismic impedance tensor we propose a new Rayleigh wave dispersion function in a homogeneous and layered medium of the Earth, which provides an efficient way to compute the dispersion curve -- a relation between the frequencies and the phase velocities. With this newly established forward model, based on the Mixture Density Networks (MDN) we develop a machine learning based inversion approach, named as FW-MDN, for the problem of estimating the S-wave velocity from the dispersion curves. The method FW-MDN deals with the non-uniqueness issue encountered in studies that invert dispersion curves for crust and upper mantle models and attains a satisfactory performance on the dataset with various noise structure. Numerical simulations are performed to show that the FW-MDN possesses the characteristics of easy calculation, efficient computation, and high precision for the model characterization.
We consider the numerical algorithm for the two-dimensional time-harmonic elastic wave scattering by unbounded rough surfaces with Dirichlet boundary condition. A Nystr{o}m method is proposed for the scattering problem based on the integral equation method. Convergence of the Nystr{o}m method is established with convergence rate depending on the smoothness of the rough surfaces. In doing so, a crucial role is played by analyzing the singularities of the kernels of the relevant boundary integral operators. Numerical experiments are presented to demonstrate the effectiveness of the method.
Numerical approximation of a general class of nonlinear unidirectional wave equations with a convolution-type nonlocality in space is considered. A semi-discrete numerical method based on both a uniform space discretization and the discrete convolution operator is introduced to solve the Cauchy problem. The method is proved to be uniformly convergent as the mesh size goes to zero. The order of convergence for the discretization error is linear or quadratic depending on the smoothness of the convolution kernel. The discrete problem defined on the whole spatial domain is then truncated to a finite domain. Restricting the problem to a finite domain introduces a localization error and it is proved that this localization error stays below a given threshold if the finite domain is large enough. For two particular kernel functions, the numerical examples concerning solitary wave solutions illustrate the expected accuracy of the method. Our class of nonlocal wave equations includes the Benjamin-Bona-Mahony equation as a special case and the present work is inspired by the previous work of Bona, Pritchard and Scott on numerical solution of the Benjamin-Bona-Mahony equation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
