A field expansions method for scattering by periodic multilayered media

The interaction of acoustic and electromagnetic waves with periodic structures plays an important role in a wide range of problems of scientific and technological interest. This contribution focuses upon the robust and high-order numerical simulation of a model for the interaction of pressure waves generated within the earth incident upon layers of sediment near the surface. Herein described is a boundary perturbation method for the numerical simulation of scattering returns from irregularly shaped periodic layered media. The method requires only the discretization of the layer interfaces (so that the number of unknowns is an order of magnitude smaller than finite difference and finite element simulations), while it avoids not only the need for specialized quadrature rules but also the dense linear systems characteristic of boundary integral/element methods. The approach is a generalization to multiple layers of Bruno and Reitich's "Method of Field Expansions" for dielectric structures with two layers. By simply considering the entire structure simultaneously, rather than solving in individual layers separately, the full field can be recovered in time proportional to the number of interfaces. As with the original field expansions method, this approach is extremely efficient and spectrally accurate.     

Solution of the scalar wave equation over very long distances using nonlinear solitary waves: Relation to finite difference methods

The linear wave equation represents the basis of many linear electromagnetic and acoustic propagation problems. Features that a computational model must have, to capture large scale realistic effects (for over the horizon or “OTH” radar communication, for example), include propagation of short waves with scattering and partial absorption by complex topography. For these reasons, it is not feasible to use Green’s Function or any simple integral method, which neglects these intermediate effects and requires a known propagation function between source and observer. In this paper, we describe a new method for propagating such short waves over long distances, including intersecting scattered waves. The new method appears to be much simpler than conventional high frequency schemes: Lagrangian “particle” based approaches, such as “ray tracing” become very complex in 3-D, especially for waves that may be expanding, or even intersecting. The other high frequency scheme in common use, the Eikonal, also has difficulty with intersecting waves.Our approach, based on nonlinear solitary waves concentrated about centroid surfaces of physical wave features, is related to that of Whitham [1], which involves solving wave fronts propagating on characteristics. Then, the evolving electromagnetic (or acoustic) field can be approximated as a collection of propagating co-dimension one surfaces (for example, 2-D surfaces in three dimensions). This approach involves solving propagation equations discretely on an Eulerian grid to approximate the linear wave equation. However, to propagate short waves over long distances, conventional Eulerian numerical methods, which attempt to resolve the structure of each wave, require far too many grid cells and are not feasible on current or foreseeable computers. Instead, we employ an “extended” wave equation that captures the important features of the propagating waves. This method is first formulated at the partial differential equation (PDE) level, as a wave equation with an added “confining” term that involves both a positive and a negative dissipation. Once we have the stable PDE, the discrete formulation is simply a multidimensional PDE with (stable) perturbations caused by the discretization. The resulting discrete solution can then be low order and very simple and yet remain stable over arbitrarily long times. When discretized and solved on an Eulerian grid, this new method allows far coarser grids than required by conventional resolution considerations, while still accounting for the effects of varying atmospheric and topographic features. An important point is that the new method is in the same form as conventional discrete wave equation methods. However, the conventional solution eventually decays, and only the “intermediate asymptotic” solution can be used. Simply by adding an extra term, we show that a nontrivial true asymptotic solution can be obtained. A similar solitary wave based approach has been used successfully in a different problem (involving “Vorticity Confinement”), for a number of years.     

一维周期结构散射近-远场外推的两种方法

首先讨论了自由空间中一维周期结构近-远场外推的Floquet模方法,即对FDTD计算所得散射近场中输出边界上散射场进行级数展开,求出各阶Floquet模的复数幅值,再求得远区场。接着介绍了求解一维周期结构远区散射场的周期Green函数方法,即根据FDTD计算所得散射近场中输出边界上散射场,求得等效面电磁流后,再借助周期Green函数进行外推得到远区场。两种方法均仅用一个周期单元内的散射近场进行外推。计算结果验证了上述两种方法外推的有效性。     

一维周期结构散射近-远场外推的两种方法

 首先讨论了自由空间中一维周期结构近-远场外推的Floquet模方法,即对FDTD计算所得散射近场中输出边界上散射场进行级数展开,求出各阶Floquet模的复数幅值,再求得远区场。接着介绍了求解一维周期结构远区散射场的周期Green函数方法,即根据FDTD计算所得散射近场中输出边界上散射场,求得等效面电磁流后,再借助周期Green函数进行外推得到远区场。两种方法均仅用一个周期单元内的散射近场进行外推。计算结果验证了上述两种方法外推的有效性。     

Resonant field amplification near the RWM stability boundary in a tokamak

Evolution of a resistive wall mode (RWM)—a magnetic field perturbation produced by a plasma and partially stabilized by a conducting wall—is considered. It is assumed that there is a small resonant harmonic in the spectrum of the static error field. It is shown that the effect of this harmonic on the dynamics of stable RWMs increases as the plasma approaches the RWM stability boundary. The error field is “amplified” during the transition through this boundary. The smaller the rotation velocity of the perturbation and the longer the time during which the plasma stays near the stability boundary, the stronger this amplification is.     

Electromagnetic waves in uniaxial crystals with metallized boundaries: Mode conversion,pure reflections,and bulk polaritons

A theory is constructed for the reflection of plane electromagnetic waves in uniaxial crystals with a positive definite permittivity tensor and an arbitrarily oriented metallized boundary. The problem is solved both for general-position orientations corresponding to three-partial reflection and for special conditions allowing two-partial reflections: mode conversions when the incident and reflected waves belong to different sheets of the refraction surface and “pure” reflections when both waves belong to the same sheet. The space of pure reflections is shown to be formed by two types of optical-axis orientations: arbitrary directions in the plane of the crystal surface and in the plane of incidence. The configurations of the conversion surface for optically positive and negative crystals are investigated. A subspace of pure reflections that transform into one-partial bulk polaritons with the energy flux parallel to the surface at grazing incidence has been found. The domain of existence of such bulk eigenmodes is bounded by two “lines” of solutions. These are any directions along the boundary containing the optical axis for ordinary polaritons and the direction along the projection of the optical axis onto the surface at an arbitrary orientation of the axis with respect to the boundary for extraordinary polaritons.     

Electromagnetic field properties in the vicinity of a massive wormhole

It is proved that not only massless but also traversable massive wormholes can have electromagnetic “hair.” An analysis is also presented of the passage from a traversable wormhole to the limit of a Reissner-Nordström black hole, with the corresponding disappearance of “hair.” A general method is developed for solving stationary axisymmetric Maxwell’s equations in the field of a massive, spherically symmetric wormhole. As a particular example of application of the method, a solution is found to the axisymmetric magnetostatic problem for a current loop in the field of the Bronnikov-Ellis-Morris-Thorne wormhole.     

A high accuracy algorithm for 3D periodic electromagnetic scattering

We develop a highly accurate numerical method for scattering of 3D electromagnetic waves by doubly periodic structures. We approximate scattered fields using the Müller boundary integral formulation of Maxwell’s equations. The accuracy is achieved as singularities are isolated through the use of partitions of unity, leaving smooth, periodic integrands that can be evaluated with high accuracy using trapezoid sums. The removed singularities are resolved through a transformation to polar coordinates. The method relies on the ideas used in the free space scattering algorithm of Bruno and Kunyansky.     

Nonlinear wave propagation through a ferromagnet with damping in (2+1) dimensions

We investigate how dissipation and nonlinearity can affect the electromagnetic wave propagating through a saturated ferromagnet in the presence of an external magnetic field in (2+1) dimensions. The propagation of electromagnetic waves through a ferromagnet under an external magnetic field in the presence of dissipative effect has been studied using reductive perturbation method. It is found that to the lowest order of perturbation the system of equations for the electromagnetic waves in a ferromagnet can be reduced to an integro-differential equation.     

[Russian Text Ignored]

The stability of the electromagnetic plasma confinement by powerful external s-polarized pump waves is considered. The parametric excitation of standing electromagnetic waves along the plasma boundary with frequencies close to the frequency of the pump wave leads to a periodic density modulation of the plasma boundary. The density disturbances along the direction of the external wave field are connected to the excitation of transverse p-polarized surface waves while the modulation in the direction perpendicular to the pump field are created by the parametric interaction between the external wave and s-polarized trapped leaking oscillations. Only when the leaking waves are excited the scale length of the modulation is larger than half the free space wave length of the incident radiation.     

On the conductivity of composites with two-dimensional periodic structure

A method is suggested of successive solution of the problem on the conductivity of two-dimensional periodic systems with inclusions of arbitrary shape. The complex potential outside of the inclusions is expressed in terms of the Weierstrass zeta function and its derivatives. The field induced on a separate inclusion is described using the matrix of multipole polarizabilities. The “joining” of potentials is performed at a distance such that R < ρ < a, where R is the characteristic dimension (maximum “radius”) of the inclusion and a is the half-period of the lattice. The approach suggested enables one to find exact virial expansions for the conductivity of other effective characteristics of similar systems as well.     

Dynamical System Approach to a Coupled Dispersionless System: Localized and Periodic Traveling Waves

We investigate the dynamical behavior of a coupled dispersionless system describing a current-conducting string with infinite length within a magnetic field. Thus, following a dynamical system approach, we unwrap typical miscellaneous traveling waves including localized and periodic ones. Studying the relative stabilities of such structures through their energy densities, we find that under some boundary conditions, localized waves moving in positive directions are more stable than periodic waves which in contrast stand for the most stable traveling waves in another boundary condition situation.     

On the electronic structure of liquid metals

We consider the electronic structure of liquid metals as a problem of multiple scattering. The total scattering operator is developed in terms of the reaction operators of a single scattering event, having a ‘standing wave’ boundary condition. This method has several advantages with respect to the approach where the total scattering operator is developed in terms of the scattering operators of a single scattering. Since the reaction operators are hermitian, the occurring perturbation expansions have a real perturbation parameter, which is more easily handled. Furthermore, if in the expansion of the total scattering operator, only terms up to the second order are maintained, this approximate expression satisfies already the ‘optical theorem’, just as the exact total scattering operator. Finally, spurious damping terms are eliminated at the outset. For very dense systems we obtain practically non-attenuated waves and this can be understood as if the damping due to a single scattering event is compensated by the reaction of the surrounding scatterers.     

Enhancement of the Local Electromagnetic Field over Planar “Particles” Formed on the Surface of a Polar Crystal

Small “particles” of an open surface were formed on a SiC polar crystal with openings in the opaque metal mask covering the sample. The dimensions of the holes were about surface phonon polariton wavelength. Such a sample was irradiated with an electromagnetic wave ( λ = 10.68μm) at a frequency close to the lattice resonance of SiC. A significant enhancement in the field amplitude of surface phonon polariton waves was detected over such “particles” compared to the amplitude over an infinite open surface of SiC. Such a phenomenon, observed by us in the IR band, is similar to plasmon resonance on small metal particles in the visible band, but the lateral resolution of the ASNOM used (no worse than 30 nm at 10 μm) makes the obtained field distribution more detailed. The maps of the local field amplitude and phase obtained on SiC surface with ASNOM are in a good quantitative agreement with simulations using the Green’s function.     

Penetration of electromagnetic waves into magnetoactive plasma

Landau studied the longitudinal field penetration into a semi-bounded plasma considering specularly reflecting electrons at the boundary. Kurilko and Popov and then Gorman examined the same problem considering a fraction “p” of all incident electrons reflecting specularly and the remaining portion difusely. Lastly, Mason treated the same problem at length and gave solutions for specular and diffuse reflections. Here the penetration of transverse electromagnetic waves into magnetoactive semi-bounded plasma has been studied and the electric field for specular reflection of electrons is obtained.     

Band structure of collective modes and effective properties of binary magnonic crystals

In this paper a theoretical study of the band structure of collective modes of binary ferromagnetic systems formed by a submicrometric periodic array of cylindrical cobalt nanodots partially or completely embedded into a permalloy ferromagnetic film is performed. The binary ferromagnetic systems studied are two-dimensional periodic, but they can be regarded as three-dimensional, since the magnetization is non uniform also along the z direction due to the contrast between the saturation magnetizations of the two ferromagnetic materials along the thickness. The dynamical matrix method, a finite-difference micromagnetic approach, formulated for studying the dynamics in one-component periodic ferromagnetic systems is generalized to ferromagnetic systems composed by F ferromagnetic materials. It is then applied to investigate the spin dynamics in four periodic binary ferromagnetic systems differing each other for the volume of cobalt dots and for the relative position of cobalt dots within the primitive cell. The dispersion curves of the most representative frequency modes are calculated for each system for an in-plane applied magnetic field perpendicular to the Bloch wave vector. The dependence of the dispersion curves on the cobalt quantity and position is discussed in terms of distribution of effective “surface magnetic charges” at the interface between the two ferromagnetic materials. The metamaterial properties in the propagative regime are also studied (1) by introducing an effective magnetization and effective “surface magnetic charges” (2) by describing the metamaterial wave dispersion of the most representative mode in each system within an effective medium approximation and in the dipole-exchange regime. It is also shown that the interchange between cobalt and permalloy does not necessarily lead to an interchange of the corresponding mode dispersion. Analogously to the case of electromagnetic waves in two-dimensional photonic crystals, the degree of localization of the localized collective modes is expressed in terms of an energy concentration factor.     

合成冷/热射流控制超声速边界层流动稳定性

为了探究超声速边界层流动稳定性及其转捩控制机理,提出基于合成冷/热射流的边界层速度-温度耦合控制方法,并通过数值模拟研究了Ma=4.5超声速平板边界层不稳定波的传播,采用线性稳定性理论中的时间模式分析了壁面吹吸、射流温度、扰动频率、扰动振幅等对不稳定波控制效果的影响.结果表明:无射流控制时,边界层内同时存在不稳定的第一模态扰动波和第二模态扰动波,且二维波形式的第二模态占主导地位;壁面吹吸作用下,仅出现更加不稳定的第二模态,第一模态被抑制;速度-温度耦合控制下,射流温度对扰动模态的不稳定区域大小及扰动增长率影响显著,射流温度与来流温度不同时,温度的脉动使得流动转捩为湍流的速度加快,边界层速度型更加饱满,抗干扰能力增强,流动稳定性提高;高频的吹吸扰动对流场的控制效果优于低频扰动,扰动频率超过400 Hz时,第二模态扰动波时间增长率降低,扰动分量对边界层速度剖面和温度剖面的修正加快,第二模态更加稳定;扰动振幅减小为主流速度的1%时,仅出现时间增长率较小的第二模态,控制效果较好,进一步减小时,第一模态重新出现,并且波数范围与第二模态先重合后分离,对应的时间增长率先增加后减小.研究结果为边界层转捩控制技术提供了新的思路.     

Study of wave propagation in various kinds of plasmas using adapted simulation methods,with illustrations on possible future applications

The understanding of wave propagation in turbulent magnetized plasmas can be rather complex, particularly if they are inhomogeneous and time-dependent. Simulation can be a useful tool for wave propagation studies, provided that the “model” equations take into account the characteristics of the medium relevant for the studied problem and that the numerical scheme including boundary conditions is stable and accurate enough. The choices for the model equations and the corresponding schemes are analyzed and discussed as a function of various parameters, such as the order of the numerical scheme and the number of grid points per wavelength. A quick review of the up-to-date numerical developments is given on the sheath boundary conditions and on the perfect matching layer in anisotropic media. Possible developments of plasma diagnostics conclude this state-of-the-art of simulations of electromagnetic waves in plasmas.     

Graphene under perpendicular incidence of electromagnetic waves: Gaps and band structure

We study the energy spectrum for the problem of graphene's carriers under perpendicular incidence of electromagnetic waves. To obtain the spectrum, we solve the Dirac equation under such a field. Then a generalized Mathieu equation is obtained. The limiting cases of strong and weak fields, as well as long and short wavelengths, are analyzed. The energy spectrum is obtained numerically by using the Whittaker method, and a perturbation method is used to obtain the analytical separatrix curve between allowed and forbidden states. The energy spectrum presents bands, and a gap at the Fermi energy. The gap is linear on the electromagnetic wave intensity.