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.

Axial acoustic radiation force on an elastic spherical shell near an impedance boundary for zero-order quasi-Bessel-Gauss beam

Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery,which requires the precise prediction of the acoustic radiation force under various circumstances to improve the system efficiency.The acoustic radiation force exerted by a zero-order quasi-Bessel-Gauss beam on an elastic spherical shell near an impedance boundary is theoretically and numerically studied in this study.By means of the finite series method and the image theory,a zero-order quasi-Bessel-Gauss beam is expanded in terms of spherical harmonic functions,and the exact solution of the acoustic radiation force is derived based on the acoustic scattering theory.The acoustic radiation force function,which represents the radiation force per unit energy density and per unit cross-sectional surface,is especially investigated.Some simulated results for a polymethyl methacrylate shell and an aluminum shell are provided to illustrate the behavior of acoustic radiation force in this case.The simulated results show the oscillatory property and the negative radiation force caused by the impedance boundary.An appropriate relative thickness of the shell can generate sharp peaks for a polymethyl methacrylate shell.Strong radiation force can be obtained at small half-cone angles and the beam waist only affects the results at high frequencies.Considering that the quasi-Bessel-Gauss beam possesses both the energy focusing property and the non-diffracting advantage,this study is expected to be useful in the development of acoustic tweezers,contrast agent micro-shells,and drug delivery applications.

Tunable acoustic radiation pattern assisted by effective impedance boundary

The acoustic wave propagation from a two-dimensional subwavelength slit surrounded by metal plates decorated with Helmholtz resonators （HRs） is investigated both numerically and experimentally in this work. Owing to the presence of HRs, the effective impedance of metal surface boundary can be manipulated. By optimizing the distribution of HRs, the asymmetric effective impedance boundary will be obtained, which contributes to generating tunable acoustic radiation pattern such as directional acoustic beaming. These dipole-like radiation patterns have high radiation efficiency, no finger- print of sidelobes, and a wide tunable range of the radiation pattern directivity angle which can be steered by the spatial displacements of HRs.

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.

Acoustic radiation force on a rigid cylinder between two impedance boundaries in a viscous fluid

Acoustofluidic technology combines acoustic and microfluidic technologies to realize particle manipulation in microchannels driven by acoustic waves,and the acoustic radiation force(ARF)with boundaries is important for particle manipulation in an acoustofluidic device.In the work reported here,the ARF on a free cylinder immersed in a viscous fluid with an incident plane wave between two impedance boundaries is derived analytically and calculated numerically.The influence of multiple scattering between the particle and the impedance boundaries is described by means of image theory,the finite-series method,and the translational addition theorem,and multiple scattering is included partly in image theory.The ARF on a free rigid cylinder in a viscous fluid is analyzed by numerical calculation,with consideration given to the effects of the distances from cylinder edge to boundaries,fluid viscosity,cylinder size,and boundary reflectivity.The results show that the interaction between the two boundaries and the cylinder makes the ARF change more violently with different frequencies,while increasing the viscosity can reduce the amplitude of the ARF in boundary space.This study provides a theoretical basis for particle manipulation by the ARF in acoustofluidics.

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.

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.

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.

Axisymmetric fundamental solutions for a finite layer with impeded boundaries

Axisymmetric fundamental solutions that are applied in the consolidation calculations of a finite clay layer with impeded boundaries were derived. Laplace and Hankel integral transforms were utilized with respect to time and radial coordinates, respectively in the analysis. The derivation of fundamental solutions considers two boundary value problems involving unit point loading and ring loading in the vertical. The solutions are extended to circular distributed and strip distributed normal load. The computation and analysis of settlements, vertical total stress and excess pore pressure in the consolidation layer subject to circular loading are presented.

Reflection of plane waves from a micropolar thermoelastic solid half-space with impedance boundary conditions

In this paper,the governing equations of linear,isotropic,homogeneous and generalized micropolar thermoelasticity are specialized in a plane.The governing equations are solved for plane harmonic wave solutions.Two separate velocity equations are obtained which indicate the existence of four plane waves with distinct speeds.A problem on reflection of plane waves from a thermally insulated/isothermal surface is considered with impedance boundary conditions.Appropriate potentials for incident and reflected waves are formulated which satisfy the boundary conditions at a plane surface.Relations between reflection coefficients as well as the expressions of energy ratios for various reflected waves are obtained.For illustration,the reflection coefficients and energy ratios of reflected waves are computed for relevant material parameters of an aluminium-epoxy composite.Effect of impedance parameters on all reflected waves is shown graphically at each angle of incidence.

ONE-DIMENSIONAL CONSOLIDATION OF LAYERED SOILS WITH IMPEDED BOUNDARIES UNDER TIME-DEPENDENT LOADINGS

On the basis of Terzaghi's one-dimensional consolidation theory, the variation of effective stress ratio in layered saturated soils with impeded boundaries under time-dependent loading was studied. By the method of Laplace transform, the solution was presented. Influences of different kinds of cyclic loadings and impeded boundaries conditions were discussed. Through numerical inversion of Laplace transform, useful illustrations were given considering several common time-dependent loadings. Pervious or impervious boundary condition is just the special case of the problem considered here. Compared with average index method,the results from the method illustrated are more accurate.

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.

An Unpreconditioned Boundary-Integral for Iterative Solution of Scattering Problems with Non-Constant Leontovitch Impedance Boundary Conditions

This paper concerns the electromagnetic scattering by arbitrary shaped three dimensional imperfectly conducting objects modeled with non-constant Leontovitch impedance boundary condition.It has two objectives.Firstly,the intrinsically wellconditioned integral equation(noted GCSIE)proposed in[30]is described focusing on its discretization.Secondly,we highlight the potential of this method by comparison with two other methods,the first being a two currents formulation in which the impedance condition is implicitly imposed and whose the convergence is quasioptimal for Lipschitz polyhedron,the second being a CFIE-like formulation[14].In particular,we prove that the new approach is less costly in term of CPU time and gives a more accurate solution than that obtained from the CFIE formulation.Finally,as expected,It is demonstrated that no preconditioner is needed for this formulation.

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.

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.

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 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.