Electromagnetic waves in dielectrics with memory

We analyse electromagnetic wave propagation in a dielectric with memory of the Maxwell-Hopkinson type. We show that the components of the electric and magnetic fields satisfy two different scalar wave equations and we first look for their harmonic plane wave solutions. Then we prove that dielectrics with memory can also support approximate Courant-Hilbert waves. We discuss the equations to be solved to get all the components of the electromagnetic field from a scalar solution from each wave equation and TE, TM harmonic plane waves are explicitly given.     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.     

Unified boundary integral equation for the scattering of elastic and acoustic waves: solution by the method of moments

A unified boundary integral equation (BIE) is developed for the scattering of elastic and acoustic waves. Traditionally, the elastic and acoustic wave problems are solved separately with different BIEs. The elastic wave case is represented in a vector BIE with the traction and displacement vectors as unknowns whereas the acoustic wave case is governed by a scalar BIE with velocity potential or pressure as unknowns. Although these two waves can be unified in the form of a partial differential equation, the unified form in its BIE counterpart has not been reported. In this work, we derive the unified BIE for these two waves and then show that the acoustic wave case can be derived from this BIE by introducing a shielding loss for small shear modulus approximation; hence only one code needs to be maintained for both elastic and acoustic wave scattering. We also derive the asymptotic Green's tensor for zero shear modulus and solve the corresponding vector equation. We employ the method of moments, which has been widely used in electromagnetics, as a numerical tool to solve the BIEs involved. Our numerical experiments show that it can also be used robustly in elastodynamics and acoustics.     

Computational simulation of wave propagation problems in infinite domains

This paper deals with the computational simulation of both scalar wave and vector wave propagation problems in infinite domains. Due to its advantages in simulating complicated geometry and complex material properties, the finite element method is used to simulate the near field of a wave propagation problem involving an infinite domain. To avoid wave reflection and refraction at the common boundary between the near field and the far field of an infinite domain, we have to use some special treatments to this boundary. For a wave radiation problem, a wave absorbing boundary can be applied to the common boundary between the near field and the far field of an infinite domain, while for a wave scattering problem, the dynamic infinite element can be used to propagate the incident wave from the near field to the far field of the infinite domain. For the sake of illustrating how these two different approaches are used to simulate the effect of the far field, a mathematical expression for a wave absorbing boundary of high-order accuracy is derived from a two-dimensional scalar wave radiation problem in an infinite domain, while the detailed mathematical formulation of the dynamic infinite element is derived from a two-dimensional vector wave scattering problem in an infinite domain. Finally, the coupled method of finite elements and dynamic infinite elements is used to investigate the effects of topographical conditions on the free field motion along the surface of a canyon.     

Low-frequency surface wave propagation along plane boundaries in fluid-saturated porous media

A method of the mechanics of a fluid-saturated porous medium is used to study the propagation of harmonic surface waves along the free boundary of such a medium, along the boundary between a porous medium and a fluid, and along the boundary between two porous half-spaces. It is shown that, at low frequencies (i.e., for waves with frequencies lower than the Biot characteristic frequency), the corresponding dispersion equations in zero-order approximation are reduced to the equations for an “equivalent” elastic medium. For the wave numbers of surface waves, corrections taking into account the generation of longitudinal waves of the second kind at the boundary are calculated. Examples of numerical solutions of dispersion equations for rock are presented.     

Lamb waves propagation in elastic plane layers with a joint strip

The Lamb waves are used for the ultrasonic characterization of welds because of their relative long-range propagation. In this paper, a simplified model of a weld-strip between two identical semi-infinite elastic layers is investigated. The reflected and transmitted ultrasonic fields are expressed by modal series whose coefficients are obtained by application of orthogonality relation. Comparisons with solutions obtained by finite elements wave propagation simulations are made. The energy balance between the incident and the scattered waves is also used to verify the accuracy of the obtained modal amplitudes.     

Radiation of scalar waves by the discontinuities in the antenna geometries

The electrical scalar potential satisfies the Helmholtz equation for time-harmonic waves and Aharonov-Bohm showed that the electromagnetic potentials have physical reality. For these reasons, the scalar potential represents scalar wave propagation in space. The gradient of the scalar wave leads to the electric field intensity, but it does not create the magnetic field. The criteria for the radiation of the scalar waves are studied on various linear antennas. The radiation integrals of the antennas are evaluated asymptotically by taking into account the edge effects. The radiation diagrams are plotted numerically and the phase velocities of the scalar waves are investigated in the near field.     

A Stiffness Reduction Method for efficient absorption of waves at boundaries for use in commercial Finite Element codes

Commercially available Finite Element packages are being used increasingly for modelling elastic wave propagation problems. Demand for improved capability has resulted in a drive to maximise the efficiency of the solver whilst maintaining a reliable solution. Modelling waves in unbound elastic media to high levels of accuracy presents a challenge for commercial packages, requiring the removal of unwanted reflections from model boundaries. For time domain explicit solvers, Absorbing Layers by Increasing Damping (ALID) have proven successful because they offer flexible application to modellers and, unlike the Perfectly Matched Layers (PMLs) approach, they are readily implemented in most commercial Finite Element software without requiring access to the source code. However, despite good overall performance, this technique requires the spatial model to extend significantly outside the domain of interest. Here, a Stiffness Reduction Method (SRM) has been developed that operates within a significantly reduced spatial domain. The technique is applied by altering the damping and stiffness matrices of the system, inducing decay of any incident wave. Absorbing region variables are expressed as a function of known model constants, helping to apply the technique to generic elastodynamic problems. The SRM has been shown to perform significantly better than ALID, with results confirmed by both numerical and analytical means.     

Analysis of the space-time dynamics of acoustic fields on the surface of a solid

The space-time dynamics of an acoustic field produced by a piezoelectric transducer in a pulsed mode is studied. The detection of acoustic fields is achieved using a Doppler laser interferometer. It is shown that, for a pictorial representation of the dynamics of a pulsed process, it is convenient to use the patterns of instantaneous spatial field distributions within the scanning area, the observation of which at successive instants makes it possible to trace the acoustic field variations on a time scale considerably smaller than the period of the ultrasonic wave. Experimental data demonstrating the process of phase propagation along the sample boundary as a function of time are presented. They are in good agreement with theoretical results obtained by using various methods of acoustic field calculation and different scalar potential distributions over the transducer surface. It is shown that the velocity of phase propagation along the sample boundary, which is mainly determined by the wave front curvature of the elastic wave incident on the sample surface, can considerably exceed the wave velocity in the unbounded medium.     

两侧有固体层负载时板中Lamb波的传播

本文研究了薄板二面有固体导负载时板中Ｌａｍｂ波的传播，从弹性波理论出发并结合应的边界条件，导出板中Ｌａｍｂ波的色散方程，数值计算表示，不管作为自由状态时板中Ｌａｍｂ波相速（板厚取定时）是大于或小于外层固体的声表面波波速，板中对称及反对称模式的Ｌａｍｂ波相速都随着外层固体层厚度增加而变化并且渐近于外层固体的声表面波波速，数值计算又表明，对很薄的板，板中对称及反对称模式的相速均随负载板的厚度呈线性变化     

Potential surface waves at the boundary of a metal with a magnetoactive plasma at finite pressure

The propagation of surface-type potential waves along the interfacial boundary of a plasma with an ideally conducting metal in an external magnetic field perpendicular to the boundary is examined. It is shown that a necessary condition for the existence of these waves in the system is a finite gas kinetic pressure. Dispersion relations for these waves and expressions for the penetration depth of the wave fields into the plasma are obtained, and they are studied numerically for various plasma parameters. The frequency region for propagation of these waves is found. It is also shown that in a nonzero external magnetic field a system of this kind has a range of frequencies in which the wave is a generalized surface wave. Zh. Tekh. Fiz. 69, 30–33 (November 1999)     

Solutions of scalar and electromagnetic wave equations in the metric of gravitational and electromagnetic waves

The wave equation for a scalar field ? and vector potentialA* are solved in the background metric of a gravitational wave. The corresponding solutions when the metric is generated by a plane electromagnetic wave, is obtained from these solutions. The solution for the scalar wave is discussed in detail. It is found that because of the interaction, two new waves are generated in the lower order approximations. One of them has the same phase dependence as the original wave while the other shows a transient character. There is no interaction when the waves are along the same direction.     

Excitation of an elastic half-space by a buried line source of conical waves

A formal solution is obtained to the problem of a buried line source of conical waves propagating at a constant phase velocity c in an isotropic elastic half space. By applying the boundary conditions at the free surface, it is determined that the reflected field, in addition to the incident field, requires addition of a scalar potential and two components of the vector potential. The latter is in contrast to the case of cylindrical waves where only one component of the vector potential is needed. The formal solution for the conical wave source goes over to that for the two dimensional cylindrical wave case in the limit of infinite phase velocity c.     

Elastic nonlinearity parameters of tetragonal crystals

The equations of motion for the propagation of finite amplitude elastic waves in crystals of tetragonal symmetry have been derived starting from the expression for the elastic strain energy. The equations have been solved for a finite amplitude sinusoidal wave propagating along the pure mode directions which are [100], [110] and [001] for the tetragonal group TI. The solutions corresponding to longitudinal wave propagation yield expressions for the amplitudes of the fundamental and generated second harmonic for these directions in terms of certain combinations of second and third order elastic constants of the medium. The results will aid the experimenter to determine these constants using ultrasonic harmonic generation technique.     

On Virtual Phonons,Photons, and Electrons

A macroscopic realization of the peculiar virtual particles is presented. The classical Helmholtz and the Schrödinger equations are differential equations of the same mathematical structure. The solutions with an imaginary wave number are called evanescent modes in the case of elastic and electromagnetic fields. In the case of non-relativistic quantum mechanical fields they are called tunneling solutions. The imaginary wave numbers point to strange consequences: The waves are non-local, they are not observable, and they are described as virtual particles. During the last two decades QED calculations of the solutions with an imaginary wave number have been experimentally confirmed for phonons, photons, and electrons. The experimental proofs of the predictions of non-relativistic quantum mechanics and the Wigner phase time approach for the elastic, electromagnetic and Schrödinger fields will be presented in this article. The results are zero time in the barrier and an interaction time (i.e. a phase shift) at the barrier interfaces. The measured tunneling time scales approximately inversely with the particle energy. Actually, the tunneling time is given only by the barrier boundary interaction time, as zero time is spent inside a barrier.     

Probability of boundary conditions in quantum cosmology

One of the main interest in quantum cosmology is to determine boundary conditions for the wave function of the universe which can predict observational data of our universe. For this purpose, we solve the Wheeler–DeWitt equation for a closed universe with a scalar field numerically and evaluate probabilities for boundary conditions of the wave function of the universe. To impose boundary conditions of the wave function, we use exact solutions of the Wheeler–DeWitt equation with a constant scalar field potential. These exact solutions include wave functions with well known boundary condition proposals, the no-boundary proposal and the tunneling proposal. We specify the exact solutions by introducing two real parameters to discriminate boundary conditions, and obtain the probability for these parameters under the requirement of sufficient e-foldings of the inflation. The probability distribution of boundary conditions prefers the tunneling boundary condition to the no-boundary boundary condition. Furthermore, for large values of a model parameter related to the inflaton mass and the cosmological constant, the probability of boundary conditions selects an unique boundary condition different from the tunneling type.     

Effect of the surface free energy on the behaviour of surface and guided waves

The aim of this paper is to evaluate the influence of the surface free energy upon the propagation of the eigenmodes of structures, by studying successively (a) the Rayleigh wave for an elastic half-space, (b) the Lamb waves for an elastic layer, and (c) the guided modes for a tri-layer structure (e.g., metal/adhesive/metal). The surface free energy is a parameter which appears in the jump conditions of stresses and displacements at each interface, and which consequently modifies the eigenmodes, solutions of the boundary conditions system. As expected, the Rayleigh wave is dispersive and its velocity increases when the surface free energy increases. In the same way, the velocity of Lamb waves also increases except at normal angle of propagation where the surface free energy does not arise. Moreover, near the Rayleigh angle, the behaviour of the A₀ and S₀ Lamb modes varies strongly according to the surface free energy. Similar results are observed for the tri-layer structure.     

Acoustic properties of a dielectric with cavities in a polarizing electric field

The effect of an external polarizing electric field on the shear wave propagation in a centrosymmetric crystal with electrostriction, whose body is penetrated with parallel cylindrical cavities (pores), is considered. The cavities are distributed throughout the crystal at random and with a low density. The waves are assumed to be polarized along the cavity generatrices, and the wave propagation occurs in the elastic isotropy plane, which is orthogonal to the cavity axes. The external field is assumed to be axial. Possibility of controlling the propagation of shear waves by the polarizing field is demonstrated for the case of metallized cavity surfaces.     

充水有限长圆柱薄壳声散射:Ⅰ.理论

研究水中充水有限长圆柱薄壳的声散射。采用弹性薄壳理论结合边界条件导出了有限长圆柱薄壳倾斜入射时散射声场的解析解,内部充水的作用反映在系统总阻抗中增加了内部流体阻抗。数值计算表明,在反向散射声场中,刚性散射部分仅在正横附近有比较大的贡献,而斜入射时起主要作用的是弹性散射部分。散射波相速度频散曲线的计算表明,与内部真空的情况相比内部充水后除正横附近的纵波和剪切波以外,增加了一组内部流体引起的附加波,其对散射声场贡献非常重要。在计算的频率-角度谱图中,内部流体附加波呈现"碗"形共振曲线。      

The general integral expressions for on-axis nonparaxial vectorial spherical waves diffracted at a circular aperture

Based on the vectorial Rayleigh-Sommerfeld integrals, the general propagation integral expressions for on-axis nonparaxial vectorial spherical wave diffracted at a circular aperture are derived. The results are strict integral formulae for the light field on the axis, valid for either strong or weak focusing, both far and near zones, and for the systems in which the size of the aperture is comparable to or smaller than the wavelength. Thus, it has the advantage for general application. For convergent spherical waves, the numerical calculation results are compared with those obtained by using the integral formulae of scalar paraxial approximation and scalar nonparaxial approximation, confirming the consistence in the situation of scalar approximation.