Search results

The boundary integral equation (BIE) numerical technique is applied to several practical pressure vessels and piping problems. Axisymmetric and three‐dimensional formulations of the BIE method for linear elastic stress analysis are reviewed. Isoparametric quadratic elements which exhibit excellent modelling capabilities are used to discretize the surfaces. Several three‐dimensional and axisymmetric structures are analysed.

This paper presents a modification of the classical boundary integral equation method (BIEM) for two‐dimensional potential boundary‐value problem. The proposed modification consists in describing the boundary geometry by means of Hermite curves. As a result of this analytical modification of the boundary integral equation (BIE), a new parametric integral equation system (PIES) is obtained. The kernels of these equations include the geometry of the boundary. This new PIES is no longer defined on the boundary, as in the case of the BIE, but on the straight line for any given domain. The solution of the new PIES does not require boundary discretization as it can be reduced merely to an approximation of boundary functions. To solve this PIES a pseudospectral method has been proposed and the results obtained compared with exact solutions.

The purpose of this paper is to apply rectangular Bézier surface patches directly into the mathematical formula used to solve boundary value problems modeled by Laplace's equation. The mathematical formula, called the parametric integral equation systems (PIES), will be obtained through the analytical modification of the conventional boundary integral equations (BIE), with the boundary mathematically described by rectangular Bézier patches.

The paper presents the methodology of the analytic connection of the rectangular patches with BIE. This methodology is a generalization of the one previously used for 2D problems.

In PIES the paper separates the necessity of performing simultaneous approximation of both boundary shape and the boundary functions, as the boundary geometry has been included in its mathematical formalism. The separation of the boundary geometry from the boundary functions enables to achieve an independent and more effective improvement of the accuracy of both approximations. Boundary functions are approximated by the Chebyshev series, whereas the boundary is approximated by Bézier patches.

The originality of the proposed approach lies in its ability to automatic adapt the PIES formula for modified shape of the boundary modeled by the Bézier patches. This modification does not require any dividing the patch into elements and creates the possibility for effective declaration of boundary geometry in continuous way directly in PIES.

The purpose of this paper is to develop an easy-to-implement and accurate fast boundary element method (BEM) for solving large-scale elastodynamic problems in frequency and time domains.

A newly developed kernel-independent fast multipole method (KIFMM) is applied to accelerating the evaluation of displacements, strains and stresses in frequency domain elastodynamic BEM analysis, in which the far-field interactions are evaluated efficiently utilizing equivalent densities and check potentials. Although there are six boundary integrals with unique kernel functions, by using the elastic theory, the authors managed to accelerate these six boundary integrals by KIFMM with the same kind of equivalent densities and check potentials. The boundary integral equations are discretized by Nyström method with curved quadratic elements. The method is further used to conduct the time-domain analysis by using the frequency-domain approach.

Numerical results show that by the fast BEM, high accuracy can be achieved and the computational complexity is brought down to linear. The performance of the present method is further demonstrated by large-scale simulations with more than two millions of unknowns in the frequency domain and one million of unknowns in the time domain. Besides, the method is applied to the topological derivatives for solving elastodynamic inverse problems.

An efficient KIFMM is implemented in the acceleration of the elastodynamic BEM. Combining with the Nyström discretization based on quadratic elements and the frequency-domain approach, an accurate and highly efficient fast BEM is achieved for large-scale elastodynamic frequency domain analysis and time-domain analysis.

The purpose of this paper is to solve generic magnetostatic problems by BEM, by studying how to use a boundary integral equation (BIE) with the double layer charge as unknown derived from the scalar potential.

Since the double layer charge produces only the potential gap without disturbing the normal magnetic flux density, the field is accurately formulated even by one BIE with one unknown. Once the double layer charge is determined, Biot‐Savart's law gives easily the magnetic flux density.

The BIE using double layer charge is capable of treating robustly geometrical singularities at edges and corners. It is also capable of solving the problems with extremely high magnetic permeability.

The proposed BIE contains only the double layer charge while the conventional equations derived from the scalar potential contain the single and double layer charges as unknowns. In the multiply connected problems, the excitation potential in the material is derived from the magnetomotive force to represent the circulating fields due to multiply connected exciting currents.

The paper gives the description of boundary element method (BEM) with subdomains for the solution of convection—diffusion equations with variable coefficients and Burgers’ equations. At first, the whole domain is discretized into K subdomains, in which linearization of equations by representing convective velocity by the sum of constant and variable parts is carried out. Then using fundamental solutions for convection—diffusion linear equations for each subdomain the boundary integral equation (in which the part of the convective term with the constant convective velocity is not included into the pseudo‐body force) is formulated. Only part of the convective term with the variable velocity, which is, as a rule, more than one order less than convective velocity constant part contribution, is left as the pseudo‐source. On the one hand, this does not disturb the numerical BEM—algorithm stability and, on the other hand, this leads to significant improvement in the accuracy of solution. The global matrix, similar to the case of finite element method, has block band structure whereas its width depends only on the numeration order of nodes and subdomains. It is noted, that in comparison with the direct boundary element method the number of global matrix non‐zero elements is not proportional to the square of the number of nodes, but only to the total number of nodal points. This allows us to use the BEM for the solution of problems with very fine space discretization. The proposed BEM with subdomains technique has been used for the numerical solution of one‐dimensional linear steady‐state convective—diffusion problem with variable coefficients and one‐dimensional non‐linear Burgers’ equation for which exact analytical solutions are available. It made it possible to find out the BEM correctness according to both time and space. High precision of the numerical method is noted. The good point of the BEM is the high iteration convergence, which is disturbed neither by high Reynolds numbers nor by the presence of negative velocity zones.

The purpose of this paper is to detect the degenerate scale of a 2D bending plate analytically and numerically.

To avoid the time-consuming scheme, the influence matrix of the boundary element method (BEM) is reformulated to an eigenproblem of the 4 by 4 matrix by using the scaling transform instead of the direct-searching scheme to find degenerate scales. Analytical degenerate scales are derived from the boundary integral equation (BIE) by using the degenerate kernel only for the circular case. Numerical results of the direct-searching scheme and the eigen system for the arbitrary shape are also considered.

Results using three methods, namely, analytical derivation, the direct-searching scheme and the 4 by 4 eigen system, are also given for the circular case and arbitrary shapes. Finally, addition of a constant for the kernel function makes original eigenvalues (2 real roots and 2 complex roots) of the 4 by 4 matrix to be all real. This indicates that a degenerate scale depends on the kernel function.

The analytical derivation for the degenerate scale of a 2D bending plate in the BIE is first studied by using the degenerate kernel. Through the reformed eigenproblem of a 4 by 4 matrix, the numerical solution for the plate of an arbitrary shape can be used in the plate analysis using the BEM.

To provide a time domain formulation for reconstruction of transient currents flowing in massive parallel conductors from magnetic data collected in the dielectric space surrounding the conductors.

A boundary integral equation (BIE) formulation involving Mitzner's and Rytov's high order surface impedance boundary conditions (SIBCs) is used. Input data of the inverse problem are the magnetic fields at given locations near the conductors. In order to validate the inversion algorithm, the magnetic field data are computed solving the direct problem with FEM for given current waveforms.

The improvement in reconstruction accuracy of the new time domain BIE formulation employing high order SIBCs has been demonstrated numerically in a simple test case. The range of validity of the technique has been extended to current pulses of longer duration and the computational burden has shown to increase only by a factor of 4.

The proposed formulation can be compared with other possible formulations, both in the time and in the frequency domain.

Based on this formulation a new current sensing technique is proposed for realization of low cost current sensors based on magnetic sensor arrays.

The inverse problem addressed in the paper has been solved for the first time.

The purpose of this paper is to supply a numerical analysis tool to solve eddy currents induced in nonlinear materials such as steel by boundary element method (BEM), and then apply it to design and analysis of power devices.

Utilizing integral formulas derived on the basis of rapid attenuation of the electromagnetic fields, the paper formulates eddy currents in steel. In the formulation, nonlinear terms are regarded as virtual sources, which are improved iteratively with the electromagnetic fields on the surface. The periodic electromagnetic fields are expanded in Fourier series and each harmonic is analyzed by BEM. The surface and internal electromagnetic fields are obtained numerically one after the other until convergence by the Newton‐Raphson method.

It is confirmed that this approach gives accurate solutions with meshes much larger than the skin depth and therefore is adequate to apply to a large‐scale application.

The eddy current is formulated by utilizing the impedance boundary condition in order to meet a large‐scale application, and so solutions near the edge are poor. In the case of better solutions being required, some modifications are necessary.

To lessen computer memory consumption, the parallel component of the currents to the steel surface is analyzed as a 2D problem and the normal component is obtained from the parallel component. One 2D equation for one analyzing region is discretized by dividing the region into layers adaptively and then solved. Next, another is solved sequentially. This method gives a compatible numerical analysis tool with finite element method.

This paper is concerned with the application of radial basis function networks (RBFNs) as interpolation functions for all boundary values in the boundary element method (BEM) for the numerical solution of heat transfer problems. The quality of the estimate of boundary integrals is greatly affected by the type of functions used to interpolate the temperature, its normal derivative and the geometry along the boundary from the nodal values. In this paper, instead of conventional Lagrange polynomials, interpolation functions representing these variables are based on the “universal approximator” RBFNs, resulting in much better estimates. The proposed method is verified on problems with different variations of temperature on the boundary from linear level to higher orders. Numerical results obtained show that the BEM with indirect RBFN (IRBFN) interpolation performs much better than the one with linear or quadratic elements in terms of accuracy and convergence rate. For example, for the solution of Laplace's equation in 2D, the BEM can achieve the norm of error of the boundary solution of O(10⁻⁵) by using IRBFN interpolation while quadratic BEM can achieve a norm only of O(10⁻²) with the same boundary points employed. The IRBFN‐BEM also appears to have achieved a higher efficiency. Furthermore, the convergence rates are of O(h^1.38) and O(h^4.78) for the quadratic BEM and the IRBFN‐based BEM, respectively, where h is the nodal spacing.