Search results

The time domain BEM/FEM coupling procedure is applied to 2‐D multi‐domain fluid–structure interaction problems. The fluid domain is acoustic and modeled by taking advantage of the BEM scheme that is suitable to either finite or infinite domains. The structure is modeled by elastodynamic finite elements that can be either linear or nonlinear. The input impact, which can be either plane waves or non‐plane waves, can either be forces acting directly on the fluid–structure system or be explosion sources in the fluid. The far field or near field explosion sources, which are difficult to be simulated for finite element analysis, are very easy to be simulated here by boundary element modeling as internal sources. The stability problem is solved by using the linear θ method, which makes the BEM scheme stable. The numerical results are compared with analytical solutions for two examples.

The purpose of this paper is to present a novel strategy used for acceleration of free-vibration analysis, in which the hierarchical matrices structure and Compute Unified Device Architecture (CUDA) platform is applied to improve the performance of the traditional dual reciprocity boundary element method (DRBEM).

The DRBEM is applied in forming integral equation to reduce complexity. In the procedure of optimization computation, ℋ-Matrices are introduced by applying adaptive cross-approximation method. At the same time, this paper proposes a high-efficiency parallel algorithm using CUDA and the counterpart of the serial effective algorithm in ℋ-Matrices for inverse arithmetic operation.

The analysis for free-vibration could achieve impressive time and space efficiency by introducing hierarchical matrices technique. Although the serial algorithm based on ℋ-Matrices could obtain fair performance for complex inversion operation, the CUDA parallel algorithm would further double the efficiency. Without much loss in accuracy according to the examination of the numerical example, the relative error appeared in approximation process can be fixed by increasing degrees of freedoms or introducing certain amount of internal points.

The paper proposes a novel effective strategy to improve computational efficiency and decrease memory consumption of free-vibration problems. ℋ-Matrices structure and parallel operation based on CUDA are introduced in traditional DRBEM.

An indirect symmetric Galerkin BEM (SGBEM) is applied to 2D potential problems in this paper. Based on the assumption that solutions from different methods should be the same, the hypersingular matrix appeared in SGBEM is approximately expressed by those matrices appeared in asymmetric Galerkin BEM (AGBEM). As only strong and weak singularities need to be solved, the problem becomes much simpler. The space derivatives of potential are expressed with a set of new meaning distributed flux, which will produce the same potential on the boundary position for Ω in the unbounded domain Ω+Ω′, so that hypersingularity will not appear for boundary points. Therefore, there is no need of C^1,α for the spatial interpolation function (no Galerkin integration can be used for this purpose). Formulations for both the steady‐state and time‐domain potential problems are given. Three numerical examples are analyzed to demonstrate the effectiveness and accuracy of the proposed indirect method.

The spectral boundary element method for solving two‐dimensional transient heat conduction problems is developed. This method is combined with the fast Fourier transform (FFT) to convert the solution between the time and frequency domains. The fundamental solutions in the frequency domain, required for the present method, are discussed. The resulting line integrations in the frequency domain are discretized using constant boundary elements and used in a Fortran boundary element program. Three examples are used to illustrate the accuracy and effectiveness of the method in both the frequency and time domains. First, the frequency domain solution procedure is verified using the steady‐state example of a semi‐infinite half space with a heat flux applied to a patch of the surface. This spectral boundary element method is then applied to the problem of a circular hole in an infinite solid subjected to a time‐varying heat flux, and solutions in both the frequency and time domains are presented. Finally, the method is used to solve the circular hole problem with a convection boundary condition. The accurary of these results leads to the conclusion that the spectral boundary element method in conjunction with the FFT is a viable option for transient problems. In addition, this spectral approach naturally produces frequence domain information which is itself of interest.

The Convolution Quadrature Method (CQM) has been widely applied to solve transient elastodynamic problems because of its stability and generality. However, the CQM suffers from the problems of huge memory requirement in case of direct implementation in time domain or CPU time in case of its reformulation in Laplace domain. The purpose of this paper is to combine the CQM with the pseudo-initial condition method (PICM) to achieve a good balance between memory requirement and CPU time.

The combined methods first subdivide the whole analysis into a few sub-analyses, which is dealt with the PICM, namely, the results obtained by previous sub-analysis are used as the initial conditions for the next sub-analysis. In each sub-analysis, the time interval is further discretized into a number of sub-steps and dealt with the CQM. For non-zero initial conditions, the pseudo-force method is used to transform them into equivalent body forces. The boundary face method is employed in the numerical implementation. Three examples are analyzed. Results are compared with analytical solutions or FEM results and the results of reformulated CQM.

Results demonstrate that the computation time and the storage requirement can be reduced significantly as compared to the CQM, by using the combined approach.

The combined methods can be successfully applied to the problems of long-time dynamic response, which requires a large amount of computer memory when CQM is applied, while preserving the CQM stability. If the number of time steps is high, then the accuracy of the proposed approach can be deteriorated because of the pseudo-force method.

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.

The purpose of this paper is the derivation and efficient implementation of surface impedance boundary conditions (SIBCs) for nonlinear magnetic conductors.

An approach based on perturbation theory is proposed, which expands to nonlinear problems the methods already developed by the authors for linear problems. Differently from the linear case, for which the analytical solution of the diffusion equation in the semi-infinite space for the magnetic field is available, in the nonlinear case the corresponding nonlinear diffusion equation must be solved numerically. To this aim, a suitable smooth map is defined to reduce the semi-infinite computational domain to a finite one; then the diffusion equation is solved by a Galerkin method relying on basis functions constructed via the push-forward of a Lagrangian polynomial basis whose degrees of freedom are collocated at Gauss–Lobatto nodes. The use of such basis in connection with a suitable under-integration naturally leads to mass-lumping without impacting the order of the method. The solution of the diffusion equation is coupled with a boundary element method formulation for the case of parallel magnetic conductors in terms of E and B fields.

The results are validated by comparison with full nonlinear finite element method simulations showing very good accordance at a much lower computational cost.

Limitations of the method are those arising from perturbation theory: the introduced small parameter must be much less than one. This implies that the penetration depth of the magnetic field into the magnetic and conductive media must be much smaller than the characteristic size of the conductor.

The efficient implementation of a nonlinear SIBC based on a perturbation approach is proposed for an electric and magnetic field formulation of the two-dimensional problem of current driven parallel solid conductors.

The purpose of this paper is to describe an extension of the boundary element method (BEM) and the dual reciprocity boundary element method (DRBEM) formulations developed for one- and two-dimensional steady-state problems, to analyse transient convection–diffusion problems associated with first-order chemical reaction.

The mathematical modelling has used a dual reciprocity approximation to transform the domain integrals arising in the transient equation into equivalent boundary integrals. The integral representation formula for the corresponding problem is obtained from the Green’s second identity, using the fundamental solution of the corresponding steady-state equation with constant coefficients. The finite difference method is used to simulate the time evolution procedure for solving the resulting system of equations. Three different radial basis functions have been successfully implemented to increase the accuracy of the solution and improving the rate of convergence.

The numerical results obtained demonstrate the excellent agreement with the analytical solutions to establish the validity of the proposed approach and to confirm its efficiency.

Finally, the proposed BEM and DRBEM numerical solutions have not displayed any artificial diffusion, oscillatory behaviour or damping of the wave front, as appears in other different numerical methods.

The purpose of this paper is to implement the Anderson acceleration for different formulations of eletromagnetic nonlinear problems and analyze the method efficiency and strategies to obtain a fast convergence.

The paper is structured as follows: the general class of fixed point nonlinear problems is shown at first, highlighting the requirements for convergence. The acceleration method is then shown with the associated pseudo-code. Finally, the algorithm is tested on different formulations (finite element, finite element/boundary element) and material properties (nonlinear iron, hysteresis models for laminates). The results in terms of convergence and iterations required are compared to the non-accelerated case.

The Anderson acceleration provides accelerations up to 75 per cent in the test cases that have been analyzed. For the hysteresis test case, a restart technique is proven to be helpful in analogy to the restarted GMRES technique.

The acceleration that has been suggested in this paper is rarely adopted for the electromagnetic case (it is normally adopted in the electronic simulation case). The procedure is general and works with different magneto-quasi static formulations as shown in the paper. The obtained accelerations allow to reduce the number of iterations required up to 75 per cent in the benchmark cases. The method is also a good candidate in the hysteresis case, where normally the fixed point schemes are preferred to the Newton ones.

A boundary element method for solving problems dealing with the dynamic analysis of thin elastoplastic flexural plates of arbitrary geometry and conditions inside the domain is proposed here. All possible edge boundary conditions, with any interior support conditions, such as isolated points (columns), lines (walls) or regions (patches) can be treated. The formulation by using the static fundamental solution of the problem leads to a system of boundary integral equations involving values of the layers along the edge. The solution of the problem with interior support conditions is achieved by an elimination of the unknown boundary layers. Subsequently, a descritization leads to a system of simultaneous algebraic equations which is solved numerically. A step‐by‐step time integration algorithm is employed to evaluate the dynamic inelastic response of the plate. Several examples are presented to illustrate the efficiency of the method.