Cylindrical Wave Functions Research Articles

A fast algorithm for the two-dimensional Helmholtz transmission problem with large multiple scattering configurationsa).

We develop an efficient three-stage algorithm for simulating multiple acoustic scattering by two-dimensional configurations comprising large numbers of penetrable scatterers. Our approach is based on a boundary integral equation reformulation of the Helmholtz transmission partial differential equation, and a reduction of the boundary integral system for computationally efficient evaluation of wave interactions between scatterers. A key ingredient of our algorithm is to represent the interactions between scatterers using expansions of cylindrical wavefunctions. For large numbers of scatterers, this approach facilitates the application of the fast multipole method, leading to linear complexity of the algorithm with respect to the number of scatterers. Numerical results demonstrate the efficiency of our algorithm for configurations containing a few hundred to hundreds of thousands of individual scatterers.

Revisiting loop quantum gravity with selfdual variables: Hilbert space and first reality condition

Abstract We consider the quantization of gravity as an SL(2,C) gauge theory in terms of Ashtekar's selfdual variables and reality conditions for the spatial metric (RCI) and its evolution (RCII). We start from a holomorphic phase space formulation. It is then natural to push for a quantization in terms of holomorphic wave functions. Thus we consider holomorphic cylindrical wave functions over SL(2,C) connections. 

We use an overall phase ambiguity of the complex selfdual action to obtain Poisson brackets that mirror those of the real theory. We then show that there is a representation of the corresponding canonical commutation relations the space of holomorphic cylindrical functions. We describe a class of cylindrically consistent measures that implements RCI. We show that spin networks with SU(2) intertwiners form a basis for gauge invariant states. They are still mutually orthogonal, but the normalisation is different than for the Ashtekar-Lewandowski measure for SU(2). 

We do not consider RCII in the present article. Work on RCII is ongoing and will be presented elsewhere.

Scattering characteristics of a perfect electromagnetic conducting cylindrical object covered with a general Bi-isotropic layer

Abstract An analytic theory has been developed for the scattering of electromagnetic plane wave from a perfect electromagnetic conducting (PEMC) cylindrical object coated with a general bi-isotropic (BI) material. The proposed problem has been solved using cylindrical vector wave function expansion approach along with the application of tangential boundary conditions. Analytic expressions of the scattering coefficients have been derived in their simplest forms. It is seen that by proper selection of admittance of PEMC core, electromagnetic parameters of BI coating, and coating thickness, one can optimize the scattering characteristics for specific applications. It is shown that the specific types of BI and strong chiral-coated PEMC cylinders having certain coating thicknesses can be used to significantly enhance the co-polarized forward scattering while keeping the cross-polarized forward scattering very small. Such types of enhanced co-polarized forward scattering are preferred in point-to-point communication. Some interesting features have been discussed where co-polarized and cross-polarized backscattering may be suppressed, which find applications in radar engineering problems and stealth technology.

Vibration characteristics of eccentric annular plates

An effective theoretical method for investigating the vibration characteristics of eccentric annular plates is developed. The vibration equation of eccentric annular plates is solved by utilizing the Bessel function as the vibration mode function and combining it with the translational addition theorem for cylindrical vector wave functions. Substituting the boundary conditions into the vibration mode function yields the frequency equation. The natural frequencies and corresponding vibration modes of the plate are obtained by determining the natural frequency coefficients. Compared with the results obtained from the finite element method (FEM) and vibration experiment, the correctness of the mode shapes and natural frequencies of the plate obtained by the present method are validated. The influences of several parameters such as the eccentric distance, the structural parameters, and various boundary conditions on the vibration characteristics of the plates are discussed. An analysis on the mode composition of the eccentric annular plates is also presented. The research results are useful for the practical engineering applications of the eccentric annular plates.

Conical near-field acoustic holography based on cylindrical wave function of variable radius and nonconformal plane measurement

In practical acoustical measurement for large cones, maintaining coaxial and conformal between the holographic surface and reconstructed surface is usually tedious. To overcome this flaw, a conical near-field acoustic holography (NAH) based on cylindrical wave function of variable radius and nonconformal plane measurement is proposed in this paper. First, cylindrical wave function of variable radius (VR)–based modified statistically optimal cylindrical NAH (SOCNAH) method, entitled VR-SOCNAH, is proposed, which is more applicable to conical surfaces. Second, with the measured sound pressure data on the holographic plane, sound pressure on the conical conformal surface is first reconstructed using statistically optimal planar NAH (SOPNAH), and sound pressure on the reconstructed surface is then reconstructed using VR-SOCNAH. Third, fixed and variable radius measurement methods are comparatively studied, and the former is adopted. Furthermore, a measurement parameter optimization method based on orthogonal experiment is proposed to determine appropriate value ranges. Finally, the proposed NAH is applied to a test bed, and the robustness and accuracy of sound field reconstruction for large cones are significantly enhanced, which provides reliable scientific basis and engineering guidance for noise monitoring and control of mechanical equipment.

Electromagnetic Wave Scattering by a Multiple Core Model of Composite Cylindrical Wires at Oblique Incidence

A complex cylindrical structure consisting of a group of parallel stratified circular lossy dielectric cylinders, embedded in a dielectric circular cylindrical region and surrounded by unbounded dielectric space, is considered in this paper. The scattering of electromagnetic (EM) plane waves by the aforementioned configuration was studied; the EM waves impinged obliquely upon the structure and were arbitrarily polarized. The formulation used was based on the boundary-value approach coupled with the generalized separation of variables method. The EM field in each region of space was expanded in cylindrical wave-functions. Furthermore, the translational addition theorem of these functions was applied in order to match the EM field components on any cylindrical interface and enforce the boundary conditions. The end result of the analysis is an infinite set of linear algebraic equations with the wave amplitudes as unknowns. The system is solved by the truncation of series and unknowns and then matrix inversion; thus, we provide a semi-analytical solution for the scattered far-field and, as a consequence, for the scattering cross section of the complex cylindrical structure. The numerical results focus on calculations of the electric- and magnetic-field intensity of the far-field as well as of the total scattering cross section of several geometric configurations that fall within the aforementioned general structure. The effect of the geometrical and electrical characteristics of the structure on the scattered field was investigated. Specifically, the cylinders’ size and spacing, their conductivity and permittivity as well as the incidence direction were modified in order to probe how these variations are imprinted on scattering. Moreover, comparisons with previously published results, as well as convergence tests, were performed; all tests and comparisons proved to be successful.

An Asymptotic High-Frequency Solution for Scattering from an Electrically Wide Triangular Cavity

An asymptotic solution based on high-frequency approximations is proposed to determine the scattered waves from a wide empty isosceles triangular cavity. The modal method based on cylindrical wavefunction expansion with the physical optics technique is used to find analytical expressions for the unknown expansion coefficients and significantly improve the time efficiency of calculations. Some assumptions and simplifications are made to reduce the complexity of the problem while still being accurate for wide triangular cavities. Comparisons are achieved to illustrate the validity and time efficiency of the suggested solution.

Light scattering by an infinitely long circular cylinder above a plane interface

An exact analytical method is presented to calculate light scattering by an infinitely long circular cylinder above a plane interface. The transformation relations of rectangular and cylindrical vector wave functions are used to deal with the coupling scattering between the circular cylinder and plane interface. By applying the boundary conditions over the cylindrical surface, the unknown expansion coefficients of the scattered field by the circular cylinder in terms of cylindrical vector wave functions are determined. Mueller matrix elements are computed, and also compared with those measured for a plane wave with normal incidence.

Transmission characteristics of terahertz Bessel vortex beams through a multi-layered anisotropic magnetized plasma slab

The transmission of terahertz (THz) Bessel vortex beams through a multi-layered anisotropic magnetized plasma slab is investigated by using a hybrid method of cylindrical vector wave functions (CVWFs) and Fourier transform. On the basis of the electromagnetic boundary conditions on each interface, a cascade form of expansion coefficients of the reflected and transmitted fields is obtained. Taking a double Gaussian distribution of the plasma density as an example, the influences of the applied magnetic field, the incident angle and polarization mode of the incident beams on the magnitude, OAM mode and polarization of the transmitted beams are analyzed in detail. The results indicate that the applied magnetic field has a major effect upon the polarization state of the transmitted fields but not upon the transmitted OAM spectrum. The incident angle has a powerful influence upon both the amplitude profile and the OAM spectrum of the transmitted beam. Furthermore, for multiple coaxial vortex beams, an increase of the maximum value of the plasma density causes more remarkable distortion of both the profile and OAM spectrum of the transmitted beam. This research makes a stable foundation for the THz OAM multiplexing/demultiplexing technology in a magnetized plasma environment.

Mitigation of railway-induced vibrations by using periodic wave impeding barriers

Railway-induced vibrations can cause significant environmental issues. This paper proposes an efficient analytical method to investigate the mitigation of railway-induced vibrations by using periodic barriers in a layered half-space. The general solutions for layered ground and multiple inclusions are derived by using the potential decomposition and multiple scattering theory. The conversion equation between cylindrical and exponential functions and the addition theorem are introduced to achieve the transformation between plane and cylindrical wave functions and the translation between cylindrical wave functions. Combined with the transfer matrix method, the fundamental solution for the soil-inclusion dynamic interaction in a layered half-space is derived. The railway train and track are subsequently coupled to the ground-inclusion system. Numerical results demonstrate that the phononic crystal effect induced by the periodic distribution of barriers improves the mitigation efficiency at high frequencies. The increase in the number, size, and stiffness of barriers can give a higher mitigation efficiency in a wider frequency range. The mitigation efficiency of periodic barriers can be guaranteed when their depth is shorter than half the Rayleigh wavelength in the considered frequency range. Owing to the scattering of waves at layer interfaces, the periodic barriers beneath the track have a higher efficiency than those located next to the track, which does not appear in the homogeneous half-space. The performance of periodic barriers is significantly affected by the soil stiffness of the upper shallow layer, while it is less affected by the soil stiffness of the bottom stiffer layer.

Effect of a perfectly conducting corner space on the PAFs for an absorptive dielectric circular cylinder

The photophoretic asymmetry factors (PAFs) for an absorptive dielectric circular cylinder, located near a perfectly conducting and totally reflecting corner space are derived and computed. The method used in this analysis relies on the modal expansion method in cylindrical coordinates, the classical method of images, and the translational addition theorem of cylindrical wave functions. Initially, the components of the internal electric field vector are obtained stemming from an analysis of the scattering. Subsequently, the solution is used to integrate to the normalized intensity function over the cylinder’s volume to obtain the longitudinal (L) and transverse (T) PAFs. Both TM- and TE-polarized plane progressive waves with arbitrary incidence (in the polar plane) are considered. Attention is given to varying the dimensionless size parameter of the cylinder, the angle of incidence of the incoming waves, and the dimensionless distance parameters from the corner space. Numerical examples illustrate the analysis and demonstrate the net effect of the totally reflecting corner space on the L- and T-PAFs, where negative, positive, and neutral values have been predicted. The results are relevant in applications related to the emergence of the photophoretic force and torque on an absorptive particle located near surfaces and topics in electromagnetic/optical scattering, particle manipulation and assembly, optically bound matter, light–matter interactions, and photopheresis.

Electromagnetic wave beam propagation through a chiral slab

A general scheme for investigating the electromagnetic wave beam reflection by and transmission through a chiral slab is proposed. The electromagnetic wave beams within different regions are described by appropriate cylindrical vector wave function expansions, and their expansion coefficients are determined by using the boundary conditions and projection method. Under illumination by a fundamental Gaussian beam, a Hermite–Gaussian beam and a doughnut mode beam, numerical results of the normalized field intensity distributions are provided and discussed briefly.

Dyadic Green’s function for the graphene–dielectric stack with arbitrary field and source points

In this paper, the dyadic Green’s function for a graphene–dielectric stack is formulated based on the scattering superposition method. To this end, the scattering Green’s function in each layer is expanded in terms of cylindrical vector wave functions with unknown coefficients. Using the Kronecker delta function in the field expansion, it is considered that the field and source points lie in the arbitrary layers. Afterward, recurrence relations to calculate the unknown expansion coefficients are derived by applying the impedance boundary condition at the interface of a graphene sheet surrounded by two adjacent dielectric layers. The verification of the calculated coefficients is conducted by using them in the analysis of graphene-based structures with different numbers of layers, including (1) free-standing frequency-selective surfaces and (2) parallel plates with graphene walls. A potential application of our proposed structure is investigating the interaction of donor–acceptor pairs residing in the arbitrary layers of the graphene–dielectric stack with a desired number of layers.

Quadratic cross-sections in the multiple scattering by a pair of liquid cylinders insonified by arbitrary-shaped acoustical sheets

Stemming from the law of energy conservation applied to scattering, generalized exact partial-wave series expressions are derived for the extrinsic and intrinsic acoustic scattering, extinction and absorption cross-sections for a pair of fluid/liquid-like (viscous) cylinders of arbitrary radii. The incident insonifying field is of arbitrary shape such that any structured wavefront in two-dimensions is accomodated by this formalism, in contrast with the case of plane waves. The modal expansion method in combination with the translational addition theorem for cylindrical wave functions are used to derive the exact analytical expressions for the quadratic (nonlinear) cross-sections, and numerical computations for a non-paraxial focused Gaussian acoustical sheet chosen as an example illustrate the analysis. The results clearly show the difference between the behavior of the extrinsic and intrinsic cross-sections. The formalism presented here is generalized for any acoustical sheet such that Airy, Hermite-Gaussian and other wavefronts (in 2D) can be considered provided the appropriate beam-shape coefficients are used. The analogy with the optical counterpart is also noted.

Cutoff Wavenumbers of Multilayered Gyrotropic Circular Waveguides

A coupled-field volume integral equation (CFVIE) method is developed for the calculation of the normalized cutoff wavenumbers of circular metallic walled waveguides having concentric continuously varying highly inhomogeneous gyrotropic (i.e., gyroelectric and gyromagnetic) infill. The normalized cutoff wavenumbers are obtained as the roots of a determinantal equation formed by solving the CFVIEs using the cylindrical Dini series expansion (CDSE) method where the unknown fields inside the waveguide are expanded by entire domain orthogonal Dini-type vectorial basis functions. To account for the electric boundary condition (BC) on waveguide's circular perfect electric conducting (PEC) surface, two modified 2-D tensorial Green's functions (GFs), expanded in cylindrical vector wave functions (CVWFs), are employed in the kernels of the CFVIEs. These modified 2-D tensorial GFs are constructed by enforcing, on their dyadic form, the satisfaction of the electric BC. The CDSE, along with the modified 2-D tensorial GFs, allow for the analytical integration of the volumetric-type integrals and the reduction of the CFVIEs to a set of algebraic equations. We exhaustively demonstrate the validity of the CFVIE-CDSE by a series of comparisons on the normalized cutoff wavenumbers: we first construct the solutions for obtaining the normalized cutoff wavenumbers in homogeneous gyrotropic waveguides by the separation of variables method (SVM), and second, we employ HFSS commercial software for two-layered isotropic and three-layered gyroelectric loaded waveguides. We characterize the type of modes, i.e., TE/TM or hybrid HE/EH, for each configuration presented, and discuss the efficiency of the CFVIE-CDSE method.

Complex WGM frequencies of gyroelectric cylindrical resonators

Asymptotic closed-form expressions for calculating complex transverse electric (TE)/transverse magnetic (TM) whispering gallery mode (WGM) frequencies in homogeneous gyroelectric circular cylindrical resonators of infinite length are derived. In addition, a volume integral equation (VIE)-cylindrical Dini series expansion method is extended to support the prediction of complex WGMs for continuously varying highly inhomogeneous gyroelectric circular cylindrical resonators. To this end, the entire domain orthogonal vectorial basis of VIE is extended to support very large indices of the involved Dini-type cylindrical vector wave functions via asymptotic closed-form expressions. This way, the eigenbasis required to solve the VIE becomes free of numerical instabilities arising when very large orders of the involved Bessel functions are employed. The complex frequencies obtained by the asymptotic closed-form expressions for the case of the homogeneous gyroelectric resonator, as well as those obtained by the VIE when the multilayered gyroelectric resonator is reduced to one layer, are validated by comparisons with the complex roots extracted by numerically solving the TE/TM characteristic equations obtained from the separation of variables method, using complex root finding techniques. We demonstrate the calculation of very high-order WGM frequencies for cylindrical resonators composed of homogeneous and highly inhomogeneous permittivity profiles. This asymptotic theory constitutes a rigorous tool that may serve for verifying method of analytical regularisation–based numerical solutions for other non-circular inhomogeneous cylinders, and for interpreting experimental data for applications such as WGM lasing, refractometric sensing, and magneto-optic coupling.

Three-dimensional analytical model for vibrations from a tunnel embedded in an unsaturated half-space

An analytical model for predicting the vibrations from an underground railway tunnel embedded in an unsaturated half-space is proposed. The tunnel lining is modeled as an infinite Flugge cylindrical shell, and the unsaturated soil is modeled as a three-phase medium comprising solid grains and pores containing water and air. By using the transformation properties between the plane wave functions and the cylindrical wave functions, the model is coupled based on the boundary conditions. The developed model is validated by comparison with existing tunnel models, and the effect of saturation on the dynamic response of the tunnel–soil system is demonstrated through a case study.

Electromagnetic beam propagation through a gyrotropic anisotropic circular cylinder

A general scheme is proposed for investigating the propagation of an incident electromagnetic beam through a gyrotropic anisotropic circular cylinder. The electromagnetic fields in different regions are expanded in terms of appropriate cylindrical vector wave functions, and the unknown expansion coefficients are determined by using the boundary conditions and projection scheme. For incidence of a fundamental Gaussian beam, Hermite–Gaussian beam and zero-order Bessel beam, normalized field intensity distributions are evaluated and discussed briefly.

Рассеяние наклонно падающей плоской звуковой волны упругим цилиндром с неоднородным покрытием, находящимся вблизи плоской поверхности

В статье рассматривается задача о рассеянии плоской монохроматической звуковой волны, падающей произвольным образом на упругий круговой цилиндр с радиальнонеоднородным покрытием в присутствии плоской поверхности (абсолютно жесткой и акустически мягкой). Методом мнимых источников с применением теорем сложения для цилиндрических волновых функций получено аналитическое решение задачи. Волновые поля в содержащей среде и однородном упругом цилиндре находятся в виде разложений по волновым цилиндрическим функциям, а для нахождения полей смещений в неоднородных покрытиях построена краевая задача для системы обыкновенных дифференциальных уравнений второго порядка.Проведены численные расчеты частотных и угловых характеристик рассеянного поля для упругих однородных цилиндров с покрытием и без него, находящихся вблизи подстилающей плоскости. Выявлено существенное влияние непрерывно-неоднородных упругих покрытий на звукоотражающие свойства упругих цилиндрических тел.

Propagation of a terahertz Bessel vortex beam through a homogeneous magnetized plasma slab

This paper provides an analytic method to study propagation characteristics of a linearly polarized Bessel vortex beam through a homogeneous magnetized plasma slab. The incident Bessel vortex beam, as well as the reflected, transmitted and internal fields are expanded in terms of cylindrical vector wave functions (CVWFs). The effects of plasma thickness, electron density and magnetic induction strength on the contour profiles of the reflected and transmitted beams and orbital angular momentum (OAM) spectra are analyzed and discussed in detail. In particular, the magnetic induction strength has a significant impact on the polarization of the transmitted beam, but not on OAM state distribution. The channel capacity of THz OAM multiplexing decreases with an increase of plasma thickness and electron density.