Transient acoustic radiation from impulsively accelerated bodies by the finite element method

Bodies under impulsive motion, immersed in an infinite acoustic fluid, severely put to test any numerical method for the transient exterior acoustic problem. Such problems, in the context of the finite element method (FEM), are not well studied. FE modeling of such problems requires truncation of the infinite fluid domain at a certain distance from the structure. The volume of computation depends upon the extent of this domain as well as the mesh density. The modeling of the fluid truncation boundary is crucial to the economy and accuracy of solution and various boundary dampers have been proposed in the literature for this purpose. The second order damper leads to unsymmetric boundary matrices and this necessitates the use of an unsymmetric equation solver for large problems. The present paper demonstrates the use of FEM with zeroth, first and second order boundary dampers in conjunction with an unsymmetric, out of core, banded equation solver for impulsive motion problems of rigid bodies in an acoustic fluid. The results compare well with those obtained from analytical methods.     

边界元素法在集成电路CAD中的应用

本文采用边界元素法(BEM)计算任意多边形的二维电阻和复杂结构的二维电容。同一计算程序既能计算电阻又能计算电容,而只需作微小不同的后处理。对于边界元素法,由于仅需计算边界上的积分方程,其离散边界上的网格点数大大少于有限差分法及有限元素法所需的网格点数,使CPU执行时间显著减少,并简化了网格划分工作。此外,计算结果还指出边界元素法具有精度高和处理复杂边界能力强的优点。     

基于时域有限差分算法改进卷积完全匹配层的稳定性

为了克服时域有限差分算法中卷积完全匹配层对消逝波吸收效果差的缺点,提出一种在卷积完全匹配层后添加特殊吸收层的方法.在不增加物体与吸收层内层距离的情况下,通过调节特殊吸收层中两个衰减因子,使其为常数,并令吸收因子逐层从1增加到10,来增强吸收层对消逝波的吸收性能.平面波垂直入射到单层光子晶体的算例表明,添加了特殊吸收层的吸收边界在与散射体相距5个网格的情况下仍能够保持计算结果收敛,而传统的吸收边界则需要相距80个网格才能保证结果收敛,说明该方法提高了对消逝波的吸收性能.进一步在结构中采用此吸收边界来计算多层光子晶体的传输特性曲线,并将其与常规方法计算所得结果做比较,两种结果吻合较好.数值算例验证了该方法的有效性和正确性.     

二维多流体域耦合声场的光滑有限元解法

针对声学有限元分析中四节点等参单元计算精度低,对网格质量敏感的问题,将光滑有限元法引入到多流体域耦合声场的数值分析中,提出了二维多流体域耦合声场的光滑有限元解法。该方法在Helmholtz控制方程与多流体域耦合界面的声压/质点法向速度连续条件的基础上,得到二维多流体耦合声场的离散控制方程,并采用光滑有限元的分区光滑技术将声学梯度矩阵形函数导数的域内积分转换形函数的域边界积分,避免了雅克比矩阵的计算。以管道二维多流体域耦合内声场为数值分析算例,研究结果表明,与标准有限元相比,对单元尺寸较大或扭曲严重的四边形网格模型,光滑有限元的计算精度更高。因此光滑有限元能很好地应用于大尺寸单元或扭曲严重的网格模型下二维多流体域耦合声场的预测,具有良好的工程应用前景。      

平板大攻角平面绕流问题中涡生成的数值模拟

本文在涡团法的基础上结合有限元法及无限相似单元法在整个无界的流动区域上构造了Navier-Stokes方程的数值解,成功地克服了区域边界上凸角点邻域内解的奇性给数值计算带来的困难,同时克服了无穷远点给数值计算带来的困难。本文对平板大攻角平面绕流问题进行数值模拟,给出了初步的计算结果。     

Computational simulation of wave propagation problems in infinite domains

This paper deals with the computational simulation of both scalar wave and vector wave propagation problems in infinite domains. Due to its advantages in simulating complicated geometry and complex material properties, the finite element method is used to simulate the near field of a wave propagation problem involving an infinite domain. To avoid wave reflection and refraction at the common boundary between the near field and the far field of an infinite domain, we have to use some special treatments to this boundary. For a wave radiation problem, a wave absorbing boundary can be applied to the common boundary between the near field and the far field of an infinite domain, while for a wave scattering problem, the dynamic infinite element can be used to propagate the incident wave from the near field to the far field of the infinite domain. For the sake of illustrating how these two different approaches are used to simulate the effect of the far field, a mathematical expression for a wave absorbing boundary of high-order accuracy is derived from a two-dimensional scalar wave radiation problem in an infinite domain, while the detailed mathematical formulation of the dynamic infinite element is derived from a two-dimensional vector wave scattering problem in an infinite domain. Finally, the coupled method of finite elements and dynamic infinite elements is used to investigate the effects of topographical conditions on the free field motion along the surface of a canyon.     

Reducing the calculation workload of the Green function for electromagnetic scattering in a Schwarzschild gravitational field

When the finite difference time domain(FDTD) method is used to solve electromagnetic scattering problems in Schwarzschild space-time, the Green functions linking source/observer to every surface element on connection/output boundary must be calculated.When the scatterer is electrically extended, a huge amount of calculation is required due to a large number of surface elements on the connection/output boundary.In this paper, a method for reducing the calculation workload of Green function is proposed.The Taylor approximation is applied for the calculation of Green function.New transport equations are deduced.The numerical results verify the effectiveness of this method.     

A multiscale hp-FEM for 2D photonic crystal bands

A multiscale generalised hp-finite element method (MSFEM) for time harmonic wave propagation in bands of locally periodic media of large, but finite extent, e.g., photonic crystal (PhC) bands, is presented. The method distinguishes itself by its size robustness, i.e., to achieve a prescribed error its computational effort does not depend on the number of periods. The proposed method shows this property for general incident fields, including plane waves incident at a certain angle to the infinite crystal surface, and at frequencies in and outside of the bandgap of the PhC. The proposed MSFEM is based on a precomputed problem adapted multiscale basis. This basis incorporates a set of complex Bloch modes, the eigenfunctions of the infinite PhC, which are modulated by macroscopic piecewise polynomials on a macroscopic FE mesh. The multiscale basis is shown to be efficient for finite PhC bands of any size, provided that boundary effects are resolved with a simple macroscopic boundary layer mesh. The MSFEM, constructed by combing the multiscale basis inside the crystal with some exterior discretisation, is a special case of the generalised finite element method (g-FEM). For the rapid evaluation of the matrix entries we introduce a size robust algorithm for integrals of quasi-periodic micro functions and polynomial macro functions. Size robustness of the present MSFEM in both, the number of basis functions and the computation time, is verified in extensive numerical experiments.     

An improved immersed boundary-lattice Boltzmann method for simulating three-dimensional incompressible flows

The recently proposed boundary condition-enforced immersed boundary-lattice Boltzmann method (IB-LBM) [14] is improved in this work to simulate three-dimensional incompressible viscous flows. In the conventional IB-LBM, the restoring force is pre-calculated, and the non-slip boundary condition is not enforced as compared to body-fitted solvers. As a result, there is a flow penetration to the solid boundary. This drawback was removed by the new version of IB-LBM [14], in which the restoring force is considered as unknown and is determined in such a way that the non-slip boundary condition is enforced. Since Eulerian points are also defined inside the solid boundary, the computational domain is usually regular and the Cartesian mesh is used. On the other hand, to well capture the boundary layer and in the meantime, to save the computational effort, we often use non-uniform mesh in IB-LBM applications. In our previous two-dimensional simulations [14], the Taylor series expansion and least squares-based lattice Boltzmann method (TLLBM) was used on the non-uniform Cartesian mesh to get the flow field. The final expression of TLLBM is an algebraic formulation with some weighting coefficients. These coefficients could be computed in advance and stored for the following computations. However, this way may become impractical for 3D cases as the memory requirement often exceeds the machine capacity. The other way is to calculate the coefficients at every time step. As a result, extra time is consumed significantly. To overcome this drawback, in this study, we propose a more efficient approach to solve lattice Boltzmann equation on the non-uniform Cartesian mesh. As compared to TLLBM, the proposed approach needs much less computational time and virtual storage. Its good accuracy and efficiency are well demonstrated by its application to simulate the 3D lid-driven cubic cavity flow. To valid the combination of proposed approach with the new version of IBM [14] for 3D flows with curved boundaries, the flows over a sphere and torus are simulated. The obtained numerical results compare very well with available data in the literature.     

A conjugated infinite element method for half-space acoustic problems

Many acoustic problems (especially in environmental acoustics) involve half-space domains bounded by a plane subjected to normal admittance boundary conditions. In the "low" frequency domain, the numerical treatment of such problems usually relies on boundary element methods based on a particular Green's function suited for the half-(admittance) plane. In the present paper, an alternative hybrid finite/infinite element scheme is proposed. The method relies on a direct treatment of nonhomogeneous boundary conditions along infinite element edges (or faces). The procedure is validated through comparisons with an available reference solution.     

移动边界问题的时空域正则区域迭代配点法

考虑热传导方程的移动边界问题,其定解区域随着时间而变化。构造一种时空域上的高精度数值算法求解1+1维移动边界问题。在时空域上假设一个初始移动边界位置,构成移动边界问题的不规则计算区域,选择一个适当的正则区域（矩形区域）完全覆盖所计算的不规则区域,在正则区域上利用移动边界约束条件和固定边界条件,采用时空域重心插值配点法求解1+1维扩散方程,得到正则区域上扩散方程数据。采用二维重心插值计算假设移动边界上函数关于时间偏导数的数值,进而利用一维重心插值配点法求解移动界面控制常微分方程,得到新的假设移动界面位置。重复上述流程,最终得到问题的数值解和移动界面的最终位置。通过典型数值算例验证所建立的数值方法的有效性和数值计算精度。     

Element-free Galerkin method for a kind of KdV equation

The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Galerkin (EFG) method which is based on the moving least-squares approximation. A variational method is used to obtain discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem. It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.     

An h-adaptive finite element solver for the calculations of the electronic structures

In this paper, a framework of using h-adaptive finite element method for the Kohn–Sham equation on the tetrahedron mesh is presented. The Kohn–Sham equation is discretized by the finite element method, and the h-adaptive technique is adopted to optimize the accuracy and the efficiency of the algorithm. The locally optimal block preconditioned conjugate gradient method is employed for solving the generalized eigenvalue problem, and an algebraic multigrid preconditioner is used to accelerate the solver. A variety of numerical experiments demonstrate the effectiveness of our algorithm for both the all-electron and the pseudo-potential calculations.     

Numerical solution of nonlocal hydrodynamic Drude model for arbitrary shaped nano-plasmonic structures using Nédélec finite elements

Nonlocal material response distinctively changes the optical properties of nano-plasmonic scatterers and waveguides. It is described by the nonlocal hydrodynamic Drude model, which – in frequency domain – is given by a coupled system of equations for the electric field and an additional polarization current of the electron gas modeled analogous to a hydrodynamic flow. Recent attempt to simulate such nonlocal model using the finite difference time domain method encountered difficulties in dealing with the grad–div operator appearing in the governing equation of the hydrodynamic current. Therefore, in these studies the model has been simplified with the curl-free hydrodynamic current approximation; but this causes spurious resonances. In this paper we present a rigorous weak formulation in the Sobolev spaces H(curl) for the electric field and H(div) for the hydrodynamic current, which directly leads to a consistent discretization based on Nédélec’s finite element spaces. Comparisons with the Mie theory results agree well. We also demonstrate the capability of the method to handle any arbitrary shaped scatterer.     

Vibro-acoustic design sensitivity analysis using the wave-based method

Conventional element-based methods, such as the finite element method (FEM) and boundary element method (BEM), require mesh refinements at higher frequencies in order to converge. Therefore, their applications are limited to low frequencies. Compared to element-based methods, the wave-based method (WBM) adopts exact solutions of the governing differential equation instead of simple polynomials to describe the dynamic response variables within the subdomains. As such, the WBM does not require a finer division of subdomains as the frequency increases in order to exhibit high computational efficiency. In this paper, the design sensitivity formulation of a semi-coupled structural-acoustic problem is implemented using the WBM. Here, the results of structural harmonic analyses are imported as the boundary conditions for the acoustic domain, which consists of multiple wave-based subdomains. The cross-sectional area of each beam element is considered as a sizing design variable. Then, the adjoint variable method (AVM) is used to efficiently compute the sensitivity. The adjoint variable is obtained from structural reanalysis using an adjoint load composed of the system matrix and an evaluation of the wave functions of each boundary. The proposed sensitivity formulation is subsequently applied to a two-dimensional (2D) vehicle model. Finally, the sensitivity results are compared to the finite difference sensitivity results, which show good agreement.     

Numerical treatment of acoustic problems with the smoothed finite element method

We incorporated a cell-wise acoustic pressure gradient smoothing operation into the standard compatible finite element method and extended the smoothed finite element method (SFEM) for 2D acoustic problems. This enhancement was especially useful for dealing with the problem of an arbitrary shape with violent distortion elements. In this method, the domain integrals that involve shape function gradients can be converted into boundary integrals that involve only shape functions. Restrictions on the shape elements can be removed, and the problem domain can be discretized in more flexible ways. Numerical results showed that the proposed method achieved more accurate results and higher convergence rates than the corresponding finite element methods, even for violently distorted meshes. The most promising feature of SFEM is its insensitivity to mesh distortion. The superiority of the method is remarkable, especially when solving problems that have high wave numbers. Hence, SFEM can be beneficially applied in solving two-dimensional acoustic problems with severely distorted elements, which, in practice, have more foreground than regularity mesh.     

A new method for true and spurious eigensolutions of arbitrary cavities using the combined Helmholtz exterior integral equation formulation method

Integral equation methods have been widely used to solve interior eigenproblems and exterior acoustic problems (radiation and scattering). It was recently found that the real-part boundary element method (BEM) for the interior problem results in spurious eigensolutions if the singular (UT) or the hypersingular (LM) equation is used alone. The real-part BEM results in spurious solutions for interior problems in a similar way that the singular integral equation (UT method) results in fictitious solutions for the exterior problem. To solve this problem, a Combined Helmholtz Exterior integral Equation Formulation method (CHEEF) is proposed. Based on the CHEEF method, the spurious solutions can be filtered out if additional constraints from the exterior points are chosen carefully. Finally, two examples for the eigensolutions of circular and rectangular cavities are considered. The optimum numbers and proper positions for selecting the points in the exterior domain are analytically studied. Also, numerical experiments were designed to verify the analytical results. It is worth pointing out that the nodal line of radiation mode of a circle can be rotated due to symmetry, while the nodal line of the rectangular is on a fixed position.     

Element-free Galerkin （EFG） method for a kind of two-dimensional linear hyperbolic equation

The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin (EFG) method which is based on the moving least-square approximation for the test and trial functions. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for hyperbolic problems needs only the scattered nodes instead of meshing the domain of the problem. It neither requires any element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for two-dimensional hyperbolic problems is investigated by two numerical examples in this paper.     

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extraction algorithm is developed to map the image on the geometric domain. Identification algorithm for the location of nodes in polygon area is proposed to determine the state of the node. To promote the average quality of the mesh and the efficiency of mesh generation, a novel force-based mesh smoothing algorithm is proposed. One test case and a typical electromagnetic calculation are used to testify the effectiveness and efficiency of the proposed algorithm. The results demonstrate that the proposed algorithm can produce a high-quality mesh with less iteration.