Featured Cover Zhang, Chao; Fu, Zhuojia; Zhang, Yaoming
International journal for numerical methods in engineering,
02/2023, Letnik:
124, Številka:
3
Journal Article
Recenzirano
The cover image is based on the Research Article A novel high‐order collocation indirect boundary element method based on the Leis formulation for three‐dimensional high frequency exterior acoustic ...problems by Chao Zhang et al., https://doi.org/10.1002/nme.7138.
Using the standard approach of the Finite Element Method, numerical values of the spatial derivatives of the primal variable present low accuracy. It is due to the order reduction in the ...interpolation functions used for the approximation of the derivatives. Aiming to overcome this problem, this work uses the integral equation of the Boundary Element Method to recalculate new values of internal variables, using the nodal boundary values obtained by the Finite Element solution. The mathematical foundation for this procedure is based on the idea that the boundary integral equation is equivalent to a weighted residual sentence, and its reuse implies a new minimization of numerical errors. This strategy was previously used within the Boundary Element context to successfully recalculate nodal variables in scalar problems governed by Laplace and Poisson Equations, also used successfully for solving linear elastic problems expressed by Navier's Equation. Here, to confirm the consistency of the proposed model, computational tests are performed, in which the standard FEM results are compared with those obtained by the proposed procedure. Internal potential derivatives and also internal potential values were recalculated. To assess the quality of the solution, the results are compared to a benchmark.
A new hybrid finite element method-boundary element method (FEM-BEM) scheme is proposed for the solution of the nonlinear 2-D Laplace problem. The novelty is an original approach of the BEM where the ...domain integrals are eliminated at the discrete level by using the FEM approximation of the fundamental solution at every node of the related mesh in the linear regions. The implementation of this FEM-Green approach requires less computational burden than the standard BEM. The coupling with FEM is straightforward and appears to be more natural. The validity of the method is examined through numerical examples.
A novel method of combination of wavelet‐based boundary element method (WBEM) with frequency‐independent fundamental solutions is proposed to determine the band structures of fluid–solid phononic ...crystals (PCs) with square and triangular lattices. Integral equations established are based on the frequency‐independent fundamental solutions, which can avoid nonlinear eigenvalue problems and reduce computing time. Domain integral terms arising from the use of frequency‐independent fundamental solutions are handled with the radial integration method (RIM) and dual reciprocity method (DRM), respectively. The results show the lower precision in high frequency domain of using RIM to handle domain integral terms than that of using DRM, which can be solved by increasing Gauss point. The B‐spline wavelet on the interval and wavelet coefficients are applied to approximate the physical boundary conditions. It is proved that coupling conditions between matrix and scatterers and Bloch theorem are also applicable to wavelet coefficients. Some small matrix entries generated by wavelet vanishing moment characteristics are truncated by the provided matrix compression technique, and the influence of compressed matrices on the results is studied. Furthermore, the final Eigen equations constructed are modified to avoid numerical instability. Some examples are provided to demonstrate the accuracy and efficiency of the proposed method.
A boundary-domain integral formulation inevitably arises when the boundary element method (BEM) is applied for solving the differential equation that governs linear elastostatic problems with body ...forces. Although the domain integrals introduced by the body forces can be evaluated using internal cells this destroys the boundary-only meshing feature of BEM and makes the integration processes inefficient. This paper shows that these problems can be solved more efficiently without using internal cells by the Radial Basis Integration Method (RBIM) which employs a meshless quadrature obtained by performing boundary-only offline computations. Using RBIM, weakly singular domain integrals can be computed via a simple selective quadrature procedure, whereas, strong singular domain integrals may be computed using two schemes: The first, is based on the same selective quadrature procedure, and the second is based on the singularity separation scheme. The present method, unlike classical radial integration method (RIM), does not have problems in calculating integrals in concave or multiple-connected domains. The results obtained in some 2D linear elastostatic problems with arbitrary body forces show that this method can be as accurate as RIM but less time-consuming than the latter. This method could be applied to other engineering problems involving source terms.
This paper presents an acoustic topology optimization approach using isogeometric boundary element methods based on subdivision surfaces to optimize the distribution of sound adsorption materials ...adhering to structural surfaces. The geometries are constructed from triangular control meshes through Loop subdivision scheme, and the associated Box-spline functions that generate limit smooth subdivision surfaces are employed to discretize the acoustic boundary integral equations. The effect of sound-absorbing materials on the acoustic response is characterized by acoustic impedance boundary conditions. The optimization problem is formulated in the framework of Solid Isotropic Material with Penalization methods and the sound absorption coefficients on elements are selected as design variables. The potential of the proposed topology optimization approach for engineering prototyping is illustrated by numerical examples.
•Isogeometric boundary element method is applied to exterior acoustic analysis.•Subdivision surfaces are used in both CAD and numerical analysis.•The sound absorbing material distribution is optimized.•Burton-Miller method is adopted to prevent instabilities.•Adjoint variable method is employed for sensitivity analysis.
The paper presents a novel approach for multi-frequency acoustic topology optimization of sound-absorption materials. In this work, the isogeometric boundary element method based on subdivision ...surfaces is used to solve Helmholtz equations. To avoid time-consuming frequency sweep, we adopt a series expansion method to decouple the frequency-dependent terms from the integrand in the boundary element method, including the terms associated with the impedance boundary conditions that were introduced to model the absorption materials. Moreover, the second-order Arnoldi (SOAR) approach is employed to reduce the order of the systems. Three dimensional numerical examples were given to demonstrate the effectiveness of the proposed algorithm.
•Series expansions are utilized for frequency decoupling of acoustic BEM.•The model order is reduced by the second-order Arnoldi (SOAR) algorithm.•The isogeometric BEM is employed to integrate CAD and numerical analysis.•Multi-frequency topology optimization is performed with the isogeometric BEM.
A nonconformal twofold domain decomposition method (TDDM) based on the hybrid finite element method-boundary element method (FEM-BEM) is proposed for analyzing 3-D multiscale composite structures. ...The proposed TDDM starts by partitioning the composite object into a closed exterior boundary domain and an interior volume domain. The interior and exterior boundary value problems are coupled to each other through the Robin transmission conditions (RTCs). Both domains are then independently decomposed into subregions to facilitate computation. Specifically, FEM-DDM with the second order transmission conditions (SOTCs) is employed for the interior domain, and BEM-discontinuous Galerkin (BEM-DG) based on the combined field integral equation (CFIE) is applied for the exterior boundary domain. The proposed TDDM allows for nonconformal discretization between any touching subdomains. Without the introduction of a stabilization term that relies on a line integral over intersection of nonmatching meshes and relevant terms involving surface-line integrals, the proposed TDDM provides an effective domain decomposition (DD) preconditioner for the global system. Numerical examples are presented, and the comparisons of the simulation results with FEM-BEM confirm the validity and accuracy of TDDM. Moreover, its ability to model practical large-scale and multiscale targets is also demonstrated.
•A BEM method for multi-physics electrical contact on rough surface is presented.•Both film resistance and constriction resistance are conducted in the model.•Linear relationship between contact ...conductance and contact stiffness is verified.•Probability density function of contact pressure and temperature rise is studied.
In this study, a numerical strategy based on the boundary element method (BEM) is developed to simulate the electrical contact behavior on rough surfaces. By using the Fourier transform, the solutions of the electric potential, temperature rise, and displacement fields in the Fourier space are derived with finite boundaries, in which the electro-thermo-mechanical coupling effect is considered. With the consideration of the contact condition, an iterative algorithm is proposed to seek the solutions of rough surface electrical contact under a given contact force and total electric potential drop. The contact pressure, real contact area, electrical contact resistance, electric current, and temperature rise distribution are predicted. For the validation, the electrical contact results of a spherical punch-half space contact and a rough surface-elastic substrate contact are calculated based on the finite element method. The good agreement between the results from the BEM-based method and the finite element method indicates the accuracy of the presented method. Then, the effect of the total electric potential and the thermal expansion on the electrical contact behavior, the relationship between the contact conductance and contact stiffness, and the probability density function of the contact pressure and temperature rise distributions are evaluated using the proposed method.
•Simultaneous equations can overcome the non-uniqueness of boundary element method.•The non-uniqueness of eigenfrequency can be overcome by adding complex damping.•The complex damping method is more ...accurate than Burton-Miller method.•Singular integrals can be calculated indirectly by matrix substitution.•The accuracy of boundary element method can be improved by using indirect matrix.
In this paper, a high-precision combined coupled double boundary element method (CCD-BEM) with a full wavenumber unique solution is proposed by combining the boundary element method (BEM) equation with the virtual indirect BEM (VIBEM) equation and using the equivalent relationship between their coefficient matrices. This method only requires overcoming the non-uniqueness of the VIBEM equation to obtain the unique solution of full-wave number, thus avoiding the failure of interior-point configuration in the combined Helmholtz integral equation formulation (CHIEF) method and the calculation of hyper-singular integral in Burton-Miller method. In addition, the singular matrix in BEM can be calculated indirectly using the equivalent relationship between the BEM and VIBEM coefficient matrix, which completely avoids the direct calculation of weak singular or singular integral. The numerical results of acoustic radiation and scattering show that CCD-BEM has a lower condition number of the coefficient matrix and higher calculation accuracy than BEM.