A Moving Pseudo-Boundary MFS for Three-Dimensional Void Detection

We propose a new moving pseudo-boundary method of fundamental solutions (MFS) for the determination of the boundary of a three-dimensional void (rigid inclusion or cavity) within a conducting homogeneous host medium from overdetermined Cauchy data on the accessible exterior boundary. The algorithm for imaging the interior of the medium also makes use of radial spherical parametrization of the unknown star-shaped void and its centre in three dimensions. We also include the contraction and dilation factors in selecting the fictitious surfaces where the MFS sources are to be positioned in the set of unknowns in the resulting regularized nonlinear least-squares minimization. The feasibility of this new method is illustrated in several numerical examples.     

The Method of Fundamental Solutions for Solving Exterior Axisymmetric Helmholtz Problems with High Wave-Number

In this paper, we investigate the method of fundamental solutions (MFS) for solving exterior Helmholtz problems with high wave-number in axisymmetric domains. Since the coefficient matrix in the linear system resulting from the MFS approximation has a block circulant structure, it can be solved by the matrix decomposition algorithm and fast Fourier transform for the fast computation of large-scale problems and meanwhile saving computer memory space. Several numerical examples are provided to demonstrate its applicability and efficacy in two and three dimensional domains.     

A multilayer method of fundamental solutions for Stokes flow problems

The method of fundamental solutions (MFS) is a meshless method for the solution of boundary value problems and has recently been proposed as a simple and efficient method for the solution of Stokes flow problems. The MFS approximates the solution by an expansion of fundamental solutions whose singularities are located outside the flow domain. Typically, the source points (i.e. the singularities of the fundamental solutions) are confined to a smooth source layer embracing the flow domain. This monolayer implementation of the MFS (monolayer MFS) depends strongly on the location of the user-defined source points: On the one hand, increasing the distance of the source points from the boundary tends to increase the convergence rate. On the other hand, this may limit the achievable accuracy. This often results in an unfavorable compromise between the convergence rate and the achievable accuracy of the MFS. The idea behind the present work is that a multilayer implementation of the MFS (multilayer MFS) can improve the robustness of the MFS by efficiently resolving different scales of the solution by source layers at different distances from the boundary. We propose a block greedy-QR algorithm (BGQRa) which exploits this property in a multilevel fashion. The proposed multilayer MFS is much more robust than the monolayer MFS and can compute Stokes flows on general two- and three-dimensional domains. It converges rapidly and yields high levels of accuracy by combining the properties of distant and close source points. The block algorithm alleviates the overhead of multiple source layers and allows the multilayer MFS to outperform the monolayer MFS.     

Adaptive multilayer method of fundamental solutions using a weighted greedy QR decomposition for the Laplace equation

The mixed boundary value problem of the Laplace equation is considered. The method of fundamental solutions (MFS) approximates the exact solution to the Laplace equation by a linear combination of independent fundamental solutions with different source points. The accuracy of the numerical solution depends on the distribution of source points. In this paper, a weighted greedy QR decomposition (GQRD) is proposed to choose significant source points by introducing a weighting parameter. An index called an average degree of approximation is defined to show the efficiency of the proposed method. From numerical experiments, it is concluded that the numerical solution tends to be more accurate when the average degree of approximation is larger, and that the proposed method can yield more accurate solutions with a less number of source points than the conventional GQRD.     

Complex energy plane and topological invariant in non-Hermitian systems

Non-Hermitian systems as theoretical models of open or dissipative systems exhibit rich novel physical properties and fundamental issues in condensed matter physics. We propose a generalized local–global correspondence between the pseudo-boundary states in the complex energy plane and topological invariants of quantum states. We find that the patterns of the pseudo-boundary states in the complex energy plane mapped to the Brillouin zone are topological invariants against the parameter deformation. We demonstrate this approach by the non-Hermitian Chern insulator model. We give the consistent topological phases obtained from the Chern number and vorticity. We also find some novel topological invariants embedded in the topological phases of the Chern insulator model, which enrich the phase diagram of the non-Hermitian Chern insulators model beyond that predicted by the Chern number and vorticity. We also propose a generalized vorticity and its flipping index to understand physics behind this novel local–global correspondence and discuss the relationships between the local–global correspondence and the Chern number as well as the transformation between the Brillouin zone and the complex energy plane. These novel approaches provide insights to how topological invariants may be obtained from local information as well as the global property of quantum states, which is expected to be applicable in more generic non-Hermitian systems.     

Singular meshless method using double layer potentials for exterior acoustics

Time-harmonic exterior acoustic problems are solved by using a singular meshless method in this paper. It is well known that the source points cannot be located on the real boundary, when the method of fundamental solutions (MFS) is used due to the singularity of the adopted kernel functions. Hence, if the source points are right on the boundary the diagonal terms of the influence matrices cannot be derived. Herein we present an approach to obtain the diagonal terms of the influence matrices of the MFS for the numerical treatment of exterior acoustics. By using the regularization technique to regularize the singularity and hypersingularity of the proposed kernel functions, the source points can be located on the real boundary and therefore the diagonal terms of influence matrices are determined. We also maintain the prominent features of the MFS, that it is free from mesh, singularity, and numerical integration. The normal derivative of the fundamental solution of the Helmholtz equation is composed of a two-point function, which is one of the radial basis functions. The solution of the problem is expressed in terms of a double-layer potential representation on the physical boundary based on the potential theory. The solutions of three selected examples are used to compare with the results of the exact solution, conventional MFS, boundary element method, and Dirichlet-to-Neumann finite element method. Good numerical performance is demonstrated by close agreement with other solutions.     

Mathematical and numerical studies on meshless methods for exterior unbounded domain problems

The method of fundamental solutions (MFS) is an efficient meshless method for solving boundary value problems in an exterior unbounded domain. The numerical solution obtained by the MFS is accurate, while the corresponding matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper basis functions is proposed by the introduction of the modified Trefftz method (MTM). The concrete expressions of the corresponding condition numbers are given in mathematical forms and the solvability by these methods is mathematically proven. Thereby, the optimal parameter minimizing the condition number is also mathematically given. Numerical experiments show that the condition numbers of the matrices corresponding to the MTM and the MMFS are reduced and that the numerical solution by the MMFS is more accurate than the one by the conventional method.     

The Method of Fundamental Solutions for Solving Convection-Diffusion Equations with Variable Coefficients

A meshless method based on the method of fundamental solutions (MFS) is proposed to solve the time-dependent partial differential equations with variable coefficients. The proposed method combines the time discretization and the one-stage MFS for spatial discretization. In contrast to the traditional two-stage process, the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations. The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easier to find than the traditional approach. The numerical results show that the one-stage approach is robust and stable.     

二维声场预测的快速多极基本解法

传统基本解法在二维大规模模型的声场求解过程中,系统方程形成和求解的计算量正比于自由度N的二次方O(N2)和三次方O(N3),求解效率低;为此,引入快速多极子算法并采用广义极小残差法迭代求解,提出一种用于二维声场预测的快速多极基本解法。对无限长圆柱体及二维类车体辐射模型的仿真结果表明,当N为3000时,分别采用快速多极基本解法与传统基本解法求解所需的时间比值约为百分之四,且N越大比值越小;最终实现系统方程的形成和求解的计算量降低到正比于自由度O(N),提高了对二维大规模模型声场预测计算效率。      

Method of fundamental solutions with optimal regularization techniques for the Cauchy problem of the Laplace equation with singular points

The purpose of this study is to propose a high-accuracy and fast numerical method for the Cauchy problem of the Laplace equation. Our problem is directly discretized by the method of fundamental solutions (MFS). The Tikhonov regularization method stabilizes a numerical solution of the problem for given Cauchy data with high noises. The accuracy of the numerical solution depends on a regularization parameter of the Tikhonov regularization technique and some parameters of the MFS. The L-curve determines a suitable regularization parameter for obtaining an accurate solution. Numerical experiments show that such a suitable regularization parameter coincides with the optimal one. Moreover, a better choice of the parameters of the MFS is numerically observed. It is noteworthy that a problem whose solution has singular points can successfully be solved. It is concluded that the numerical method proposed in this paper is effective for a problem with an irregular domain, singular points, and the Cauchy data with high noises.     

An adaptive greedy technique for inverse boundary determination problem

In this paper, the method of fundamental solutions (MFS) is employed for determining an unknown portion of the boundary from the Cauchy data specified on parts of the boundary. We propose a new numerical method with adaptive placement of source points in the MFS to solve the inverse boundary determination problem. Since the MFS source points placement here is not trivial due to the unknown boundary, we employ an adaptive technique to choose a sub-optimal arrangement of source points on various fictitious boundaries. Afterwards, the standard Tikhonov regularization method is used to solve ill-conditional matrix equation, while the regularization parameter is chosen by the L-curve criterion. The numerical studies of both open and closed fictitious boundaries are considered. It is shown that the proposed method is effective and stable even for data with relatively high noise levels.     

Function projective synchronization of identical and non-identical modified finance and Shimizu–Morioka systems

In this paper, function projective synchronizations (FPS) of identical and non-identical modified finance systems (MFS) and Shimizu?CMorioka system (S-MS) are studied via active control technique. The technique is applied to construct a response system which synchronizes with a given drive system for a desired function relation between identical MFS, identical S-MS and between MFS and S-MS. The results are validated via numerical simulations.     

The method of fundamental solutions for 2D and 3D Stokes problems

A numerical scheme based on the method of fundamental solutions (MFS) is proposed for the solution of 2D and 3D Stokes equations. The fundamental solutions of the Stokes equations, Stokeslets, are adopted as the sources to obtain flow field solutions. The present method is validated through other numerical schemes for lid-driven flows in a square cavity and a cubic cavity. Test results obtained for a rectangular cavity with wave-shaped bottom indicate that the MFS is computationally efficient than the finite element method (FEM) in dealing with irregular shaped domain. The paper also discusses the effects of number of source points and their locations on the numerical accuracy.     

Three Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear FGMs with Kirchhoff and Laplace Transformation

This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials (FGMs). The three methods are, respectively, the method of fundamental solution (MFS), the boundary knot method (BKM), and the collocation Trefftz method (CTM) in conjunction with Kirchhoff transformation and various variable transformations. In the analysis, Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions. The proposed MFS, BKM and CTM are mathematically simple, easy-to-programming, meshless, highly accurate and integration-free. Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered, and the results are compared with those from meshless local boundary integral equation method (LBIEM) and analytical solutions to demonstrate the efficiency of the present schemes.     

The fractal energy measurement and the singularity energy spectrum analysis

The singularity exponent (SE) is the characteristic parameter of fractal and multifractal signals. Based on SE, the fractal dimension reflecting the global self-similar character, the instantaneous SE reflecting the local self-similar character, the multifractal spectrum (MFS) reflecting the distribution of SE, and the time-varying MFS reflecting pointwise multifractal spectrum were proposed. However, all the studies were based on the depiction of spatial or differentiability characters of fractal signals. Taking the SE as the independent dimension, this paper investigates the fractal energy measurement (FEM) and the singularity energy spectrum (SES) theory. Firstly, we study the energy measurement and the energy spectrum of a fractal signal in the singularity domain, propose the conception of FEM and SES of multifractal signals, and investigate the Hausdorff measure and the local direction angle of the fractal energy element. Then, we prove the compatibility between FEM and traditional energy, and point out that SES can be measured in the fractal space. Finally, we study the algorithm of SES under the condition of a continuous signal and a discrete signal, and give the approximation algorithm of the latter, and the estimations of FEM and SES of the Gaussian white noise, Fractal Brownian motion and the multifractal Brownian motion show the theoretical significance and application value of FEM and SES.     

Improvement of metal–ferroelectric–silicon structures without buffer layers between Si and ferroelectric films

Si-based metal–ferroelectric–semiconductor (MFS) structures without buffer layers between Si and ferroelectric films have been developed by depositing SrBi₂Ta₂O₉ (SBT) directly on n-type (100)-oriented Si. Some effective processes are adopted to improve the electrical properties of these MFS structures. Contrary to the conventional MFS structures with top electrodes directly on ferroelectrics, our MFS structures have been developed with thin dense SiO₂ films deposited between ferroelectric films and top electrodes. Due to the SiO₂ films, the leakage current densities of MFS structures are reduced to 2×10^-8 A/cm² under the bias of 5 V. The C-V electrical properties of the MFS structures are greatly improved after annealing at 400 °C in N₂ ambient for 1 h. The C-V memory windows are increased to 3 V, which probably results from the decrease of the interface trap density at the Si/SBT interface. Received: 7 September 1999 / Accepted: 24 November 1999 / Published online: 2 August 2000     

An Improved Formulation of Singular Boundary Method

This study proposes a new formulation of singular boundary method (SBM) to solve the 2D potential problems, while retaining its original merits being free of integration and mesh, easy-to-program, accurate and mathematically simple without the requirement of a fictitious boundary as in the method of fundamental solutions (MFS). The key idea of the SBM is to introduce the concept of the origin intensity factor to isolate the singularity of fundamental solution so that the source points can be placed directly on the physical boundary. This paper presents a new approach to derive the analytical solution of the origin intensity factor based on the proposed subtracting and adding-back techniques. And the troublesome sample nodes in the ordinary SBM are avoided and the sample solution is also not necessary for the Neumann boundary condition. Three benchmark problems are tested to demonstrate the feasibility and accuracy of the new formulation through detailed comparisons with the boundary element method (BEM), MFS, regularized meshless method (RMM) and boundary distributed source (BDS) method.     

Galerkin Formulations of the Method of Fundamental Solutions

In this paper, we introduce two Galerkin formulations of the Method of Fundamental Solutions (MFS). In contrast to the collocation formulation of the MFS, the proposed Galerkin formulations involve the evaluation of integrals over the boundary of the domain under consideration. On the other hand, these formulations lead to some desirable properties of the stiffness matrix such as symmetry in certain cases. Several numerical examples are considered by these methods and their various features compared.     

Importance of the (nabla) D term in frequency-resolved optical diffusion imaging

The effects of the approximation DD=0 that is often used in frequency-resolved optical diffusion imaging are examined. It is shown that this approximation can affect the performance of integral-equation-based approaches to optical diffusion imaging, such as the Born iterative method and the distorted Born iterative method. The approximation introduces errors into the calculation of data used in simulations, which can lead to misleading evaluations of reconstruction algorithms. Numerical calculations show the magnitude of these effects and the appearance of artifacts in reconstructed images when conventional inversion algorithms are applied to more accurately calculated data.     

Solving unsteady Schr?dinger equation using the improved element-free Galerkin method

By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrödinger equation. In the IEFG method, the two-dimensional (2D) trial function is approximated by the IMLS approximation, the variation method is used to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. Because the number of coefficients in the IMLS approximation is less than in the moving least-square (MLS) approximation, fewer nodes are needed in the entire domain when the IMLS approximation is used than when the MLS approximation is adopted. Then the IEFG method has high computational efficiency and accuracy. Several numerical examples are given to verify the accuracy and efficiency of the IEFG method in this paper.