Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA periodic FM-BEM for solving the acoustic transmission problems in periodic media
56
0
0.0
(
0
)
Added by
Wenhui Meng
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
This paper presents a new fast multipole boundary element method (FM-BEM) for solving the acoustic transmission problems in 2D periodic media. We divide the periodic media into many fundamental blocks, and then construct the boundary integral equations in the fundamental block. The fast multipole algorithm is proposed for the square and hexagon periodic systems, the convergence of the algorithm is analyzed. We then apply the proposed method to the acoustic transmission problems for liquid phononic crystals and derive the acoustic band gaps of the phononic crystals. By comparing the results with those from plane wave expansion method, we conclude that our method is efficient and accurate.
rate research
Read More
This paper is devoted to studying a type of contact problems modeled by hemivariational inequalities with small periodic coefficients appearing in PDEs, and the PDEs we considered are linear, second order and uniformly elliptic. Under the assumptions, it is proved that the original problem can be homogenized, and the solution weakly converges. We derive an $O(epsilon^{1/2})$ estimation which is pivotal in building the computational framework. We also show that Robin problems--- a special case of contact problems, it leads to an $O(epsilon)$ estimation in $L^2$ norm. Our computational framework is based on finite element methods, and the numerical analysis is given, together with experiments to convince the estimation.
Consider using the right-preconditioned generalized minimal residual (AB-GMRES) method, which is an efficient method for solving underdetermined least squares problems. Morikuni (Ph.D. thesis, 2013) showed that for some inconsistent and ill-conditioned problems, the iterates of the AB-GMRES method may diverge. This is mainly because the Hessenberg matrix in the GMRES method becomes very ill-conditioned so that the backward substitution of the resulting triangular system becomes numerically unstable. We propose a stabilized GMRES based on solving the normal equations corresponding to the above triangular system using the standard Cholesky decomposition. This has the effect of shifting upwards the tiny singular values of the Hessenberg matrix which lead to an inaccurate solution. Thus, the process becomes numerically stable and the system becomes consistent, rendering better convergence and a more accurate solution. Numerical experiments show that the proposed method is robust and efficient for solving inconsistent and ill-conditioned underdetermined least squares problems. The method can be considered as a way of making the GMRES stable for highly ill-conditioned inconsistent problems.
Accurate numerical simulations of interaction between fluid and solid play an important role in applications. The task is challenging in practical scenarios as the media are usually highly heterogeneous with very large contrast. To overcome this computational challenge, various multiscale methods are developed. In this paper, we consider a class of linear poroelasticity problems in high contrast heterogeneous porous media, and develop a mixed generalized multiscale finite element method (GMsFEM) to obtain a fast computational method. Our aim is to develop a multiscale method that is robust with respect to the heterogeneities and contrast of the media, and gives a mass conservative fluid velocity field. We will construct decoupled multiscale basis functions for the elastic displacement as well as fluid velocity. Our multiscale basis functions are local. The construction is based on some suitable choices of local snapshot spaces and local spectral decomposition, with the goal of extracting dominant modes of the solutions. For the pressure, we will use piecewise constant approximation. We will present several numerical examples to illustrate the performance of our method. Our results indicate that the proposed method is able to give accurate numerical solutions with a small degree of freedoms.
In this paper, an efficient iterative method is proposed for solving multiple scattering problem in locally inhomogeneous media. The key idea is to enclose the inhomogeneity of the media by well separated artificial boundaries and then apply purely outgoing wave decomposition for the scattering field outside the enclosed region. As a result, the original multiple scattering problem can be decomposed into a finite number of single scattering problems, where each of them communicates with the other scattering problems only through its surrounding artificial boundary. Accordingly, they can be solved in a parallel manner at each iteration. This framework enjoys a great flexibility in using different combinations of iterative algorithms and single scattering problem solvers. The spectral element method seamlessly integrated with the non-reflecting boundary condition and the GMRES iteration is advocated and implemented in this work. The convergence of the proposed method is proved by using the compactness of involved integral operators. Ample numerical examples are presented to show its high accuracy and efficiency.
We adapt a symmetric interior penalty discontinuous Galerkin method using a patch reconstructed approximation space to solve elliptic eigenvalue problems, including both second and fourth order problems in 2D and 3D. It is a direct extension of the method recently proposed to solve corresponding boundary value problems, and the optimal error estimates of the approximation to eigenfunctions and eigenvalues are instant consequences from existing results. The method enjoys the advantage that it uses only one degree of freedom on each element to achieve very high order accuracy, which is highly preferred for eigenvalue problems as implied by Zhangs recent study [J. Sci. Comput. 65(2), 2015]. By numerical results, we illustrate that higher order methods can provide much more reliable eigenvalues. To justify that our method is the right one for eigenvalue problems, we show that the patch reconstructed approximation space attains the same accuracy with fewer degrees of freedom than classical discontinuous Galerkin methods. With the increasing of the polynomial order, our method can even achieve a better performance than conforming finite element methods, such methods are traditionally the methods of choice to solve problems with high regularities.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
