Source author record

Jinyou Xiao

Jinyou Xiao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2016arXiv

Resolvent sampling based Rayleigh-Ritz method for large-scale nonlinear eigenvalue problems

A new algorithm, denoted by RSRR, is presented for solving large-scale nonlinear eigenvalue problems (NEPs) with a focus on improving the robustness and reliability of the solution, which is a challenging task in computational science and engineering. The proposed algorithm utilizes the Rayleigh-Ritz procedure to compute all eigenvalues and the corresponding eigenvectors lying within a given contour in the complex plane. The main novelties are the following. First and foremost, the approximate eigenspace is constructed by using the values of the resolvent at a series of sampling points on the contour, which effectively circumvented the unreliability of previous schemes that using high-order contour moments of the resolvent. Secondly, an improved Sakurai-Sugiura algorithm is proposed to solve the projected NEPs with enhancements on reliability and accuracy. The user-defined probing matrix in the original algorithm is avoided and the number of eigenvalues is determined automatically by provided strategies. Finally, by approximating the projected matrices with the Chebyshev interpolation technique, RSRR is further extended to solve NEPs in the boundary element method, which is typically difficult due to the densely populated matrices and high computational costs. The good performance of RSRR is demonstrated by a variety of benchmark examples and large-scale practical applications, with the degrees of freedom ranging from several hundred up to around one million. The algorithm is suitable for parallelization and easy to implement in conjunction with other programs and software.

preprint2016arXiv

Solving large-scale nonlinear eigenvalue problems by rational interpolation approach and resolvent sampling based Rayleigh-Ritz method

Numerical solution of nonlinear eigenvalue problems (NEPs) is frequently encountered in computational science and engineering. The applicability of most existing methods is limited by matrix structures, property of eigen-solutions, size of the problem, etc. This paper aims to break those limitations and to develop robust and universal NEP solvers for large-scale engineering applications. The novelty lies in two aspects. First, a rational interpolation approach (RIA) is proposed based on the Keldysh theorem for holomorphic matrix functions. Comparing with the existing contour integral approach (CIA), the RIA provides the possibility to select sampling points in more general regions and has advantages in improving accuracy and reducing computational cost. Second, a resolvent sampling scheme using the RIA is proposed for constructing reliable search spaces for the Rayleigh-Ritz procedure, based on which a robust eigen-solver, denoted by RSRR, is developed for solving general NEPs. RSRR can be easily implemented and parallelized. The advantages of the RIA and the performance of RSRR are demonstrated by a variety of benchmark and practical problems.

preprint2015arXiv

An efficient adaptive frequency sampling scheme for large-scale transient boundary element analysis

The frequency-domain approach (FDA) to transient analysis of the boundary element method, although is appealing for engineering applications, is computationally expensive. This paper proposes a novel adaptive frequency sampling (AFS) algorithm to reduce the computational time of the FDA by effectively reducing the number Nc of sampling frequencies. The AFS starts with a few initial frequencies and automatically determines the subsequent sampling frequencies. It can reduce Nc by more than 2 times while still preserving good accuracy. In a porous solid model with around 0.3 million unknowns, 4 times reduction of Nc and the total computational time is successfully achieved.

preprint2014arXiv

A fast directional BEM for large-scale acoustic problems based on the Burton-Miller formulation

In this paper, a highly efficient fast boundary element method (BEM) for solving large-scale engineering acoustic problems in a broad frequency range is developed and implemented. The acoustic problems are modeled by the Burton-Miller boundary integral equation (BIE), thus the fictitious frequency issue is completely avoided. The BIE is discretized by using the Nyström method based on the curved quadratic elements, leading to simple numerical implementation (no edge or corner problems) and high accuracy in the BEM analysis. The linear systems are solved iteratively and accelerated by using a newly developed kernel-independent wideband fast directional algorithm (FDA) for fast summation of oscillatory kernels. In addition, the computational efficiency of the FDA is further promoted by exploiting the low-rank features of the translation matrices, resulting in two- to three-fold reduction in the computational time of the multipole-to-local translations. The high accuracy and nearly linear computational complexity of the present method are clearly demonstrated by typical examples. An acoustic scattering problem with dimensionless wave number $kD$ (where $k$ is the wave number and $D$ is the typical length of the obstacle) up to 1000 and the degrees of freedom up to 4 million is successfully solved within 10 hours on a computer with one core and the memory usage is 24 GB.

preprint2014arXiv

A fast directional boundary element method for high frequency acoustic problems in three dimensions

A highly efficient fast boundary element method (BEM) for solving large-scale engineering acoustic problems in a broad frequency range is developed and implemented. The acoustic problems are modeled by the Burton-Miller boundary integral equation (BIE), thus the fictitious frequency issue is completely avoided. The BIE is discretized by using the collocation method with piecewise constant elements. The linear systems are solved iteratively and accelerated by using a newly developed kernel-independent wideband fast directional algorithm (FDA) for fast summation of oscillatory kernels. In addition, the computational efficiency of the FDA is further promoted by exploiting the low-rank features of the translation matrices. The high accuracy and nearly linear computational complexity of the present method are clearly demonstrated by typical examples. An acoustic scattering problem with dimensionless wave number $kD$ (where $k$ is the wave number and $D$ is the typical length of the obstacle) up to 1000 and the degrees of freedom up to 4 million is successfully solved within 4 hours on a computer with one core and the memory usage is 24.7 GB.

preprint2013arXiv

An efficient method for evaluating BEM singular integrals on curved elements with application in acoustic analysis

The polar coordinate transformation (PCT) method has been extensively used to treat various singular integrals in the boundary element method (BEM). However, the resultant integrands of the PCT tend to become nearly singular when (1) the aspect ratio of the element is large or (2) the field point is closed to the element boundary; thus a large number of quadrature points are needed to achieve a relatively high accuracy. In this paper, the first problem is circumvented by using a conformal transformation so that the geometry of the curved physical element is preserved in the transformed domain. The second problem is alleviated by using a sigmoidal transformation, which makes the quadrature points more concentrated around the near singularity. By combining the proposed two transformations with the Guiggiani's method in [M. Guiggiani, \emph{et al}. A general algorithm for the numerical solution of hypersingular boundary integral equations. \emph{ASME Journal of Applied Mechanics}, 59(1992), 604-614], one obtains an efficient and robust numerical method for computing the weakly-, strongly- and hyper-singular integrals in high-order BEM with curved elements. Numerical integration results show that, compared with the original PCT, the present method can reduce the number of quadrature points considerably, for given accuracy. For further verification, the method is incorporated into a 2-order Nyström BEM code for solving acoustic Burton-Miller boundary integral equation. It is shown that the method can retain the convergence rate of the BEM with much less quadrature points than the existing PCT. The method is implemented in C language and freely available.

preprint2013arXiv

Efficiency improvement of the frequency-domain BEM for rapid transient elastodynamic analysis

The frequency-domain fast boundary element method (BEM) combined with the exponential window technique leads to an efficient yet simple method for elastodynamic analysis. In this paper, the efficiency of this method is further enhanced by three strategies. Firstly, we propose to use exponential window with large damping parameter to improve the conditioning of the BEM matrices. Secondly, the frequency domain windowing technique is introduced to alleviate the severe Gibbs oscillations in time-domain responses caused by large damping parameters. Thirdly, a solution extrapolation scheme is applied to obtain better initial guesses for solving the sequential linear systems in the frequency domain. Numerical results of three typical examples with the problem size up to 0.7 million unknowns clearly show that the first and third strategies can significantly reduce the computational time. The second strategy can effectively eliminate the Gibbs oscillations and result in accurate time-domain responses.