Fast High Order Algorithm for Three-Dimensional Helmholtz Equation Involving Impedance Boundary Condition with Large Wave Numbers

Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.

COMPACT FOUR-COMPONENT 2-D FDFD METHOD WITH EQUIVALENT SURFACE IMPEDANCE BOUNDARY CONDITION FOR MULTILAYER METAL-COATED WAVEGUIDE

A compact four-component two-dimensional （2-D） finite-difference frequency domain （FDFD） method with the equivalent surface impedance boundary condition is used to analyze the dispersion characteristics of multilayer metal-coated waveguides. According to the equivalent surface impedance boundary condition,the relationship between transverse field components on the boundary can be easily depicted. Once the eigen equation is solved,the propagation constant can be obtained as the eigen value for a given frequency. Results of the proposed method agaree well with those of high frequency structure simulator（HFSS）.

Transient Hydroelastic Response of VLFS by FEM with Impedance Boundary Conditions in Time Domain

A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastic responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.

Application of generalized impedance boundary conditions to sommerfeld half-space problem

Generalized impedance boundary conditions are employed to simplify the solution of the Sommerfeld half-space problem. An analytical expression is derived for the Hertz potential of a vertical electric dipole over the earth’s surface, in which the earth is assumed to be a layered media or homogeneous dissipative half-space. A Sommerfeld type integral in the potential function is expressed as the sum of two parts: a zeroth order Hankel function and an absolutely convergent series of Bessel functions. In addition, two expressions in closed form are obtained as the far-field and near-field approximation of the present result.

A High Frequency Boundary Element Method for Scattering by Convex Polygons with Impedance Boundary Conditions

We consider scattering of a time harmonic incident plane wave by a convex polygon with piecewise constant impedance boundary conditions.Standard finite or boundary element methods require the number of degrees of freedom to grow at least linearly with respect to the frequency of the incident wave in order to maintain accuracy.Extending earlier work by Chandler-Wilde and Langdon for the sound soft problem,we propose a novel Galerkin boundary element method,with the approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh with smaller elements closer to the corners of the polygon.Theoretical analysis and numerical results suggest that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency of the incident wave.

A Modified Helmholtz Equation with Impedance Boundary Conditions

Here considered is the problem of transient electromagnetic scattering from overfilled cavities embedded in an impedance ground plane.An artificial boundary condition is introduced on a semicircle enclosing the cavity that couples the fields from the infinite exterior domain to those fields inside.A Green’s function solution is obtained for the exterior domain,while the interior problem is solved using finite element method.Well-posedness of the associated variational formulation is achieved and convergence and stability of the numerical scheme confirmed.Numerical experiments show the accuracy and robustness of the method.

Efficient Computation of Scattering from Targets with Negative Impedance Surface

Current surface integral equations used for computing scattering from targets with negative impedance boundary condition（IBC）are not efficient.A modified surface dual integral equation（M-SDIE）for targets with negative IBC is presented.A pure imaginary number is used to balance the formulations.It is proved that the M-SDIE is accurate and efficient with three numerical examples.The first numerical example shows that the M-SDIE is accurate compared with Mie.The second example shows that the presented SIE is efficient.In the third example,a missile head is selected to present the computing power of the M-SDIE.All the examples show that the M-SDIE is an efficient algorithm for negative IBC.

MoM-based topology optimization method for planar metallic antenna design

The metallic antenna design problem can be treated as a problem to find the optimal distribution of conductive material in a certain domain. Although this problem is well suited for topology optimization method, the volumetric distribution of conductive material based on 3D finite element method (FEM) has been known to cause numerical bottlenecks such as the skin depth issue, meshed 'air regions' and other numerical problems. In this paper a topology optimization method based on the method of moments (MoM) for configuration design of planar metallic antenna was proposed. The candidate structure of the planar metallic antenna was approximately considered as a resistance sheet with position-dependent impedance. In this way, the electromagnetic property of the antenna can be analyzed easily by using the MoM to solve the radiation problem of the resistance sheet in a finite domain. The topology of the antenna was depicted with the distribution of the impedance related to the design parameters or relative densities. The conductive material (metal) was assumed to have zero impedance, whereas the non-conductive material was simulated as a material with a finite but large enough impedance. The interpolation function of the impedance between conductive material and non-conductive material was taken as a tangential function. The design of planar metallic antenna was optimized for maximizing the efficiency at the target frequency. The results illustrated the effectiveness of the method.

Fast computation of scattering from 3D complex structures by MLFMA

This paper introduces the research work on the extension of multilevel fast multipole algorithm （MLFMA） to 3D complex structures including coating object, thin dielectric sheet, composite dielectric and conductor, cavity. The impedance boundary condition is used for scattering from the object coated by thin lossy material. Instead of volume integral equation, surface integral equation is applied in case of thin dielectric sheet through resistive sheet boundary condition. To realize the fast computation of scattering from composite homogeneous dielectric and conductor, the surface integral equation based on equivalence principle is used. Compared with the traditional volume integral equation, the surface integral equation reduces greatly the number of unknowns. To computc conducting cavity with electrically large aperture, an electric field integral equation is applied. Some numerical results are given to demonstrate the validity and accuracy of the present methods.

4-component 2-D CFDFD method in analysis of lossy circular waveguide with fractal rough surface

In this paper, equivalent surface impedance boundary condition （ESIBC）, which takes fractal parameters （D, G） into SIBC, is implemented in the 4-component 2-D compact finite difference frequency domain （2-D CFDFD） method to an- alyze the propagation characteristics of lossy circular waveguide with fractal rough surface based on Weierstrass-Mandelbrot （W-M） function. Fractal parameters’ effects on attenuation constant are presented in the 3 mm lossy circular waveguide, and the attenuation constants of the first three modes vary monotonically with scaling constant （G） and decrease as the fractal dimension （D） increasing.

Uniqueness in Determining a Ball with a Single Incoming Wave

We prove that a ball with the impedance boundary condition is uniquely determined by the far-field pattern corresponding to an incident plane wave at one given wavenumber and one given incident direction. In the uniqueness proof, the impedance parameter in the impedance boundary condition is unknown.

A METHOD FOR SOLVING THE INVERSE SCATTERING PROBLEM FOR SHAPE AND IMPEDANCE

The inverse problem considered in this paper is to determine the shape and the impedance of an obstacle from a knowledge of the time-harmonic incident field and the phase and amplitude of the far field pattern of the scattered wave in two-dimension. Single-layer potential is used to approach the scattered waves. An approximation method is presented and the convergence of the proposed method is established. Numerical examples are given to show that this method is both accurate and easy to use.

A Method for Solving the Inverse Scattering Problem for Shape and Impedance of Crack

The inverse problem considered in this paper is to determine the shape and the impedance of crack from a knowledge of the time-harmonic incident field and the corresponding far field pattern of the scattered waves in two-dimension.The combined single-and double-layer potential is used to approach the scattered waves.As an important feature,this method does not require the solution of u and δu/δv at each iteration.An approximate method is presented and the convergence of this method is proven.Numerical examples are given to show that this method is both accurate and simple to use.

IBC/FEM Analysis of Electromagnetic Scatter of Cavities Coated with Layered Medium

The Leontovich impedance boundary condition(IBC)is combined with the edge-based finite element method(FEM)in this paper to analyze the electromagnetic(EM)scattering of cavities coated with a multilayered dielectric.The IBC on the surface of the medium and the boundary integral equation on the aperture of the cavity are transformed into the third boundary condition,and then the functional of the boundary value problem is obtained.The surface impedance of the layered dielectric is calculated by the generalized reflection coefficient;hence,the multireflection of the EM wave in the dielectric is involved.As a result,the IBC is improved.Numerical results are presented,which demonstrate that the presented IBC/FEM approach is accurate and convenient for the analysis of EM scattering of open-ended cavities coated with the dielectric.