Radiation impedance matrices for rectangular interfaces within rigid baffles: Calculation methodology and applications

The generation of acoustic disturbances in a fluid of semi-infinite extent by the motion of a circular piston surrounded by a plane rigid baffle has been studied quite extensively (see (1), (2), (3), (4) and further references given in these papers). Attention has been devoted mainly to the case in which the piston executes a harmonic oscillation of small amplitude, and only comparatively recently has Oberhettinger (2) demonstrated how the time-harmonic solution can be used to solve the more general problem in which the normal velocity of the piston is an arbitrary function of time. The purpose of the present paper is two-fold. Firstly, we point out that for arbitrary normal motion of the piston the " baffled piston problem " can be solved directly, and in a particularly simple manner, by means of a technique involving integral transforms which has been applied by Mitra (5) and Eason (6) to the study of shear wave propagation in an elastic half-space. Secondly, we give a more detailed account than appears to be available in the literature of the structure of the sound pulse generated by the arbitrary normal motion of a baffled piston.

In developing accurate and stable discretization schemes for integral equations, important steps are the choices of the basis and testing functions. In the electromagnetic surface integral equations, the unknown quantities, the equivalent electric and magnetic surface current densities, are by default divergence conforming and therefore should be expanded with the divergence‐conforming basis functions. The choice of the testing functions however is less obvious and depends on the choice of the basis functions and the properties of the integral operators. We develop stable and well‐defined discretization procedures where the operators are tested either their range or on the dual of their range. This requires the use of non‐conventional finite‐elements spaces, called ‘dual spaces’, and non‐traditional testing procedures. We show with numerical examples that in most cases, testing the operators on the dual of their range leads to improved solution accuracy compared with testing the operators on their range. For the integral equations of the first kind, this agrees with the conventional Galerkin's scheme, but for the integral equations of the second kind, this requires the use of Petrov–Galerkin type discretization schemes with appropriate curl‐conforming (dual) testing functions. Copyright © 2012 John Wiley & Sons, Ltd.

Reconstruction of unknown objects requires efficient inversion for measured electromagnetic scattering data. In frequency-domain integral equation method for reconstructing dielectric objects, the volume integral equations (VIEs) are involved because the imaging domain with both true unknown objects and part of background is inhomogeneous. When solving the forward scattering integral equation (FSIE), the Nyström method is used since the traditional method of moments (MoM) with the Schaubert-Wilton-Glisson (SWG) basis function may not be convenient due to the inhomogeneity of the imaging domain. The benefits of the Nyström method include the simple implementation without using any basis and testing functions and lower requirement on geometrical discretization, so it is very suitable for inhomogeneous problems. When solving the inverse scattering integral equation (ISIE), the Gauss-Newton minimization approach (GNMA) with a multiplicative regularization method (MRM) is employed. Numerical examples for reconstructing three-dimensional dielectric objects are presented to illustrate the inversion approach.

Diagonalizability of operators occurring in linear transport theory is discussed from a general point of view. In particular, the operator A−1T, where T is the multiplication operator in L2 (−1, 1) and A is given by a formula of the type A f = f − Σ o n aj pj , is investigated. Diagonalization of this operator which is connected with one-group neutron transport is carried out in the general case that the coefficients aj are arbitrary complex numbers. Also, a peculiarity of multi-group theory, where the operator involved has a multiple continuous spectrum, is pointed out. A correct interpretation of the main result of that theory is provided.

Printed circuits in bounded media encompass a wide range of practical structures such as discontinuities in waveguides, planar circuits embedded in shielded multilayered media or even two-dimensional printed periodic structures. The Electromagnetic (EM) modeling of printed circuits in layered bounded media is performed via an Integral Equation (IE) technique. Green's functions (GFs) are specially defined to satisfy both the Boundary Conditions (BCs) imposed by the layered media and by the transverse boundary enclosing the entire structure. Finally, a system of IEs on the equivalent sources can be solved numerically by means of the Method of Moments (MoM). Each of the problems enumerated above has already been solved by other authors using IE-MoM techniques. Nevertheless, our formulation introduces a unified approach applicable to all the aforementioned problems. Due to the symmetry presented by a bounded layered media, the GF problem can be reduced into a two-dimensional transverse boundary problem and a one-dimensional transmission line problem in the normal direction. Both problems can be treated independently. This thesis proposes and fully develops an efficient technique that encompasses different laterally bounded multilayered problems with a seamless transition between them. The method is based on a modal representation of the transverse boundary problem and on the expansion of the equivalent surface currents by zero-curl & constant-charge Basis Functions (BFs). It offers a unified and versatile approach that, on one hand eliminates redundancy in the formulation and on the other hand simplifies each particular problem to the evaluation of constant coefficients or basic line integrals. Analytical solutions can be found for the combination of linear subsectional basis functions in rectangular and circular Perfect Electric Conductor (PEC) boundaries as well as for periodic lattices. This thesis then solves the problem of transmission line model in the longitudinal direction by proposing an efficient algorithm that guarantees numerical stability under a variety of known critical conditions where other already known formulations fail. In addition, it introduces alternate equivalent expressions of this formulation that allow new interpretations of the problem. Due to its practical interest, the method is applied for the EM modeling of multilayered boxed printed circuits. This motivated the implementation of a dedicated software tool for the efficient analysis of these topologies including losses. Extensive numerical experiments have been carried out to assess the validity of the aforementioned theory and some properties of test-structures (losses, mesh, etc).

A novel equivalence principle algorithm (EPA) domain decomposition method previously for the surface integral equation is explored for volume integral equation in this work. The whole object is decomposed into several subdomains, the equivalence spheres are automatically produced to enclose each subdomain. The Rao–Wilton–Glisson (RWG) and body of revolution (BoR) basis functions are defined on the equivalence sphere surface. The equivalence currents defined on the RWG basis functions are evaluated with multilevel fast multipole algorithm, and then the currents defined on the RWGs can be projected onto the BoR basis functions. Compared with surface integral equation EPA-BoR algorithm, the ratio of the number of average original unknowns in each subdomain with respect to the equivalence unknowns is relatively larger, which leads to significant reduction of computations. Numerical simulations validate the accuracy and efficiency of the proposed method.

RCS data is the most important source for SAR processing. For electrically large targets, high accuracy and efficiency calculation is a very difficult topic. Higher order moment method is a Galerkin type version of MOM based on integral equation. It uses entire domain basis function in relatively large patches. Compared with sub-domain basis function, such as RWG, this method can reduce unknowns distinctly. It provides possibilities for electrically large models simulation. This paper adopts bilinear quadrilateral elements as the modelling method and higher order hierarchical polynomial expansion as the basis function. This kind of basis function can use different orders on different elements in order to get the highest efficiency. Based on this theory, codes were made to test its accuracy and efficiency. NASA almond and sphere, two classical metallic RCS test targets, are solved with the codes. Simulation results are identical with the measured and the analytical data.

This paper presents an overview of several integral equation (IE) formulations where the use of Legendre basis functions is beneficial. Specifically, we introduce the basis functions and deal with formulations for surface IE, volume IE, volume-surface IE, and layered media, respectively. In some cases, the high complexity of a problem arises due its large size in terms of wavelengths which can be handled with the MoM-PO approach described. Finally, we present a numerical example solved with two different formulations.

Hydrodynamic parameters for a two-body axisymmetric system

The Galerkin moment method based on triangular sinusoidal currents is proposed. Its basis function is a modified version of Rao-Wilton-Glisson (RWG) vector basis function when the meshed triangles are modeled by bilinear interpolated parametric surfaces under particular oblique linear local coordinates. Because of the adoption of sinusoidal factor in the current basis function, the fourfold impedance integrals reduce into double integrals, while the double potential integrals reduce into only one dimensional integrals. The integrand of the remaining double impedance integral involves complex exponential integrals, which does not bring any computation benefit to the total quadrature time for the moment. However, the remaining one-dimensional potential integrals make it rather effective to do scattering modal analysis or near-zone energy flow analysis after solving the matrix equation. In either scene, there is a requirement for frequent calculation of the near- or far-zone fields of complex surface objects. The proposed method is applied to scattering problems of conducting square plate and sphere illuminated by a plane wave. The computed surface currents agree very well with results via RWG Galerkin method, and the method also has a little superior convergence rate.

Accurate electromagnetic (EM) modeling for packaging structures with lossy interconnects requires to consider the finite conductivity of the conducting interconnects. In the integral equation approach, the conducting loss is usually accounted for through an approximate surface impedance when the skin depth of current is small, but this approximation may not be valid for large skin depth caused by low frequencies or small conductivities. In this paper, we treat the lossy interconnects as dielectric-like media and use the two-region electric field integral equations (EFIEs) to characterize the EM feature of the interconnects. The EFIEs are solved with the method of moments in which the Rao-Wilton-Glisson (RWG) basis function and dual basis function (DBF) are used to represent the electric and magnetic current densities, respectively. The joint use of RWG and DBF can remove the need of wisely formulating integral equations and selecting testing schemes and also alleviate the low-frequency effect in a wideband modeling for multiscale interconnect structures. Numerical examples are presented to demonstrate the approach and its robustness has been verified.

Acoustic wave scattering by scatterers in layered media is studied using a new method that combines the hybrid ray–mode and boundary integral equation methods. The latter is used to formulate the scattering process and the former is used to provide the Green’s function of the layered environment. Following standard procedures of the boundary integral equation method, integral equations are formulated first for field distributions along the boundaries of the obstacles. By expressing these unknown distributions in terms of appropriate basis functions, the integral equations are then reduced to algebraic equations that are solved numerically. In the integral equations, it is necessary to compute the Green’s function of the layered environment for various arrangements of locations of source and receiver. None of these conventional approaches (rays, modes, and spectral integration) to evaluate the Green’s function are satisfactory for all possible arrangements. Only the newly developed hybrid ray–mode method is best suited for this purpose because it combines rays and modes self-consistently within a single framework and optimizes the advantages of each. Numerical implementation illustrates these aspects.

The improved hybrid discretisation of volume integral equation (VIE) combined with the characteristic basis function method is proposed to analyse the scattering from complex dielectric objects in this study. To mitigate the excessive number of discrete elements which arise from the conformal discretisation on the complex dielectric objects, the hybrid discretisation is developed. For the conventional hybrid discretisation, the conformal Schaubert–Wilton–Glisson (SWG) basis functions are used to discretise the homogeneous regions, meanwhile, the non‐conformal half‐SWG basis functions are employed to discretise the boundary regions which separate different materials or multiscale structures. The hybrid discretisation takes advantage of the merits of both conformal and non‐conformal basis functions and improves the computational efficiency for solving the VIE. In our study, a kind of edge basis function defined within boundary tetrahedral elements substitutes for the SWG basis functions to model the boundary regions, leading to the further reduction on the number of basis functions. Moreover, the characteristic basis functions extracted from the low‐level basis functions are employed to construct the reduced matrix of which the size is considerably smaller than that of the system matrix derived from the VIE. Several numerical results are presented to demonstrate the accuracy and efficiency of the proposed method.

In this article several numerical methods are used for solving the first kind of integral equations. Most of the works are based on projection methods. One of the most attractive developments in recent years has been the use of wavelets as basis functions in projection methods. Wavelet technology makes it possible to create very efficient algorithms in comparison with known regularizing algorithms. We obtained a numerical solution using the proposed method.

In this paper, a novel realization of the Integral Equation in combination with the fast Fourier transform for the CFIE is established by Fitting both the Green's function and its Gradient onto the nodes of a uniform Cartesian grid. The new method has been compared with several existing popular FFT-based methods, including the AIM, the IE-FFT, and the p-FFT. The accuracy of the proposed method is significantly higher than other FFT-based methods, and the method is not sensitive to both the grid spacing and the expansion order. The outstanding merit of the proposed method is that the fitting procedure is independent of the basis functions. Therefore, when the higher order basis functions would be adopted in the method of moments, only one fitting procedure for the Green's function and its gradient on a basis function support is needed to meet all of basis functions defined on this support. Some numerical examples are provided in this paper to demonstrate the accuracy and efficiency of the proposed method.