各向异性超常材料中倒退波的传播研究

研究了完全各向异性超常材料中的倒退波传播现象,得到了在材料本征轴和传输轴成任意角度情形下倒退波形成的条件,分析了超常材料的介电张量和磁导率张量、电磁波的偏振方式对倒退波形成和传播的影响. 在此基础上,进一步分析了几种不同色散曲线关系的各向异性超常材料中倒退波的产生情况,获得了电磁波波矢和坡印亭矢量(能流)夹角的具体表达式和倒退波传播的一般性结论. 此外,还研究了近零介电常数超常材料中倒退波的传播特性,发现在此类超常材料中倒退波只能是完美倒退波. 关键词： 超常材料 负折射 倒退波 各向异性     

Stable reformulation of transfer matrix method for wave propagation in layered anisotropic media.

The numerical instability problem in the standard transfer matrix method has been resolved by introducing the layer stiffness matrix and using an efficient recursive algorithm to calculate the global stiffness matrix for an arbitrary anisotropic layered structure. For general anisotropy the computational algorithm is formulated in matrix form. In the plane of symmetry of an orthotropic layer the layer stiffness matrix is represented analytically. It is shown that the elements of the stiffness matrix are as simple as those of the transfer matrix and only six of them are independent. Reflection and transmission coefficients for layered media bounded by liquid or solid semi-spaces are formulated as functions of the total stiffness matrix elements. It has been demonstrated that this algorithm is unconditionally stable and more efficient than the standard transfer matrix method. The stiffness matrix formulation is convenient in satisfying boundary conditions for different layered media cases and in obtaining modal solutions. Based on this method characteristic equations for Lamb and surface waves in multilayered orthotropic media have been obtained. Due to the stability of the stiffness matrix method, the solutions of the characteristic equations are numerically stable and efficient. Numerical examples are given.     

Hybrid compliance-stiffness matrix method for stable analysis of elastic wave propagation in multilayered anisotropic media

This paper presents the hybrid compliance-stiffness matrix method for stable analysis of elastic wave propagation in multilayered anisotropic media. The method utilizes the hybrid matrix of each layer in a recursive algorithm to deduce the stack hybrid matrix for a multilayered structure. Like the stiffness matrix method, the hybrid matrix method is able to eliminate the numerical instability of transfer matrix method. By operating with total stresses and displacements, it also preserves the convenience for incorporating imperfect or perfect interfaces. However, unlike the stiffness matrix, the hybrid matrix remains to be well-conditioned and accurate even for zero or small thicknesses. The stability of hybrid matrix method has been demonstrated by the numerical results of reflection and transmission coefficients. These results have been determined efficiently based on the surface hybrid matrix method involving only a subset of hybrid submatrices. In conjunction with the recursive asymptotic method, the hybrid matrix method is self-sufficient without hybrid asymptotic method and may achieve low error level over a wide range of sublayer thickness or the number of recursive operations.     

Anisotropic effects in surface acoustic wave propagation from a point source in a crystal

Strong anisotropic effects in the propagation of surface acoustic waves (SAWs) from a point-like source are studied experimentally and theoretically. Nanosecond SAW pulses are generated by focused laser pulses and detected with a cw probe laser beam at a large distance from the source compared to the SAW wavelength, which allows us to resolve fine intricate features in SAW wavefronts. In our theoretical model, we represent the laser excitation by a localized impulsive force acting on the sample surface and calculate the far-field surface response of an elastically anisotropic solid to such a force. The model simulates the measured SAW waveforms very well and accounts for all experimentally observed features. Using the data obtained for the (111) and (001) surfaces of GaAs, we describe a variety of effects encountered in the SAW propagation from a point source in crystals. The most interesting phenomenon is the existence of cuspidal structures in SAW wavefronts resulting in multiple SAW arrivals for certain ranges of the observation angle. Cuspidal edges correspond to the phonon focusing directions yielding sharp peaks in the SAW amplitude. A finite SAW wavelength results in internal diffraction whereby the SAW wavefront spreads beyond the group velocity cusps. Degeneration of a SAW into a transverse bulk wave is another strong effect influencing the anisotropy of the SAW amplitude and making whole sections of the SAW wavefront including some phonon focusing directions unobservable in the experiment. The propagation of a leaky SAW mode (pseudo-SAW) is affected by a specific additional effect i.e. anisotropic attenuation. We also demonstrate that many of the discussed features are reproduced in powder patterns, a simple technique developed by us earlier for visualization of SAW amplitude anisotropy.Received: 17 June 2003, Published online: 15 October 2003PACS: 43.35.+d Ultrasonics, quantum acoustics, and physical effects of sound - 68.35.Gy Mechanical properties; surface strains - 62.65.+k Acoustical properties of solidsA.M. Lomonosov: On leave from the General Physics Institute, 117942 Moscow, Russia     

Two models of anisotropic propagation of a cardiac excitation wave

Propagation of the action potential in the real heart is direction-dependent (anisotropic). We propose two general physical models explaining this anisotropy on the cellular level. The first, “delay” model takes into account the frequency of the cell-cell transitions in different directions of propagation, assuming each transition requires some small time interval. The second model relies on the assumption that the action potential transmits to the next cell only from the area at the pole of the previous cell. We estimated parameters of both models by doing optical mapping and fluorescent staining of cardiac cell samples grown on polymer fiber substrate. Both models gave reasonable estimations, but predicted different behaviors of the anisotropy ratio (ratio of the highest and lowest wave velocities) after addition of the suppressor of sodium channels such as lidocaine. The results of the experiment on lidocaine effect on anisotropy ratio were in favor of the first, “delay” model. Estimated average cell-cell transition delay was 240 ± 80 μs, which is close to the characteristic values of synaptic delay.     

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.     

A finite difference method for a coupled model of wave propagation in poroelastic materials

A computational method for time-domain multi-physics simulation of wave propagation in a poroelastic medium is presented. The medium is composed of an elastic matrix saturated with a Newtonian fluid, and the method operates on a digital representation of the medium where a distinct material phase and properties are specified at each volume cell. The dynamic response to an acoustic excitation is modeled mathematically with a coupled system of equations: elastic wave equation in the solid matrix and linearized Navier-Stokes equation in the fluid. Implementation of the solution is simplified by introducing a common numerical form for both solid and fluid cells and using a rotated-staggered-grid which allows stable solutions without explicitly handling the fluid-solid boundary conditions. A stability analysis is presented which can be used to select gridding and time step size as a function of material properties. The numerical results are shown to agree with the analytical solution for an idealized porous medium of periodically alternating solid and fluid layers.     

Electromagnetic wave propagation from a point source in air through a medium with a negative refractive index

Refraction of an electromagnetic wave from a point source to a medium with a negative refractive index is considered. An electric Hertz dipole that is located in air and directed parallel to a plane interface is considered as the point source. It is rigorously shown that the electromagnetic radiation propagated from the dipole into the medium with the negative refractive index is focused within a certain region in this medium. The sizes of the focusing region are determined. As a result, it is pointed out that the diffraction limit cannot be overcome by using homogeneous materials with negative refractive indices.     

Forced wave propagation and energy distribution in anisotropic laminate composites

Elastodynamic response of anisotropic laminate composite structures subjected to a force loading is evaluated based on the integral representations in terms of Green's matrices. Explicit and asymptotic expressions for guided waves generated by a given source are then obtained from those integrals by means of series expansions and the residue technique. Unlike to conventional modal expansions, such representations keep information about the source, giving an opportunity for a quantitative near- and far-field analysis of generated waves. An effective computer implementation is achieved by the use of fast and stable algorithms for the Green matrix, pole, and residue calculations. The potential of the model is demonstrated by examples of anisotropy manifestation in the directivity of radiated waves. The effect of main energy outflow in the direction of either upper- or inner-ply orientation depending on the source size and frequency is discussed.     

Wave propagation in an anisotropic nickel-based superalloy

The effects of elastic anisotropy on ultrasound propagation in a nickel-based single crystal test component are studied using a 25 MHz focused probe in a water immersion system. Anisotropy gives rise to directionally dependent acoustic wavespeeds, beam steering, acoustic energy focusing and mode conversion for normal incidence. Transverse mode echoes are particularly strong in the vicinity of crystallographic directions in which the Gaussian curvature of the slowness surface is zero and divergence of the echo amplitude is predicted on the basis of the stationary phase approximation. There are other directions where the transverse mode echoes vanish for symmetry reasons. The longitudinal mode echo amplitude also shows significant variation with direction. Overall there is good agreement between the echo signal arrival times and amplitudes we measure and calculation. Progress in applying this technique to gas turbine blades is reported.     

Ultrasonic wave propagation across a thin nonlinear anisotropic layer between two half-spaces

Boundary conditions and perturbation theory are combined to create a set of equations which, when solved, yield the reflected and transmitted wave forms in the case of a thin layer of material that is perfectly bonded between two isotropic half-spaces. The set of perturbed boundary conditions is created by first using the fully bonded boundary conditions at each of the two interfaces between the thin layer and the half-spaces. Then, by restricting the layer's thickness to be much smaller than an acoustic wavelength, perturbation theory can be used on these two sets of boundary equations, producing a set of equations which effectively treat the thin layer as a single interface via a perturbation term. With this set of equations, the full range of incident and polar angles can be considered, with results general enough to use with a layer that is anisotropic, nonlinear, or both anisotropic and nonlinear. Finally the validity of these equations is discussed, comparing the computer simulation results of this theory to results from standard methods, and looking at cases where the results (or various properties of the results) are known or can be predicted.     

Acoustic propagation from a spiral wave front source in an ocean environment

A spiral wave front source generates a pressure field that has a phase that depends linearly on the azimuthal angle at which it is measured. This differs from a point source that has a phase that is constant with direction. The spiral wave front source has been developed for use in navigation; however, very little work has been done to model this source in an ocean environment. To this end, the spiral wave front analogue of the acoustic point source is developed and is shown to be related to the point source through a simple transformation. This makes it possible to transform the point source solution in a particular ocean environment into the solution for a spiral source in the same environment. Applications of this transformation are presented for a spiral source near the ocean surface and seafloor as well as for the more general case of propagation in a horizontally stratified waveguide.     

Theory of electromagnetic wave propagation in superlattices with optically anisotropic layers

A formalism is developed to find the photon dispersion relations in superlattice systems having layers with low optical symmetry or having magnetic layers. This formalism is an exact solution of the first order Maxwell's equations including all the information for the anisotropic optical response tensors and including the coupling of the TE and TM modes. Based on a 4×4 matrix approach for solving complicated reflection and transmission problems in stratified anisotropic media and employing a plane wave expansion of the field components to take into account the periodicity of the superlattices, the photon dispersion relation can be obtained numerically with a simple algorithm. This result is useful in predicting the absence of certain electromagnetic modes along the superlattice axis, and in identifying observed resonances with a particular excitations of the system.     

Moment-screen method for wave propagation in a refractive medium with random scattering

Direct numerical solution of a parabolic equation (PE) for the second moment of the sound field in a refracting medium with random scattering is described. The method determines the mean-square sound pressure without requiring generation of random realizations of the propagation medium. The second-moment matrix is factored into components that are independently propagated with a conventional PE algorithm. A moment screen is periodically applied to attenuate the coherence of the wavefield, much as phase screens are often applied in the method of random realizations. An example involving upwind and downwind propagation in the near-ground atmosphere shows that the new direct method converges to an accurate solution faster than the method of random realizations and is particularly well suited to calculations at low frequencies.     

Love wave propagation in functionally graded piezoelectric material layer

An exact approach is used to investigate Love waves in functionally graded piezoelectric material (FGPM) layer bonded to a semi-infinite homogeneous solid. The piezoelectric material is polarized in z-axis direction and the material properties change gradually with the thickness of the layer. We here assume that all material properties of the piezoelectric layer have the same exponential function distribution along the x-axis direction. The analytical solutions of dispersion relations are obtained for electrically open or short circuit conditions. The effects of the gradient variation of material constants on the phase velocity, the group velocity, and the coupled electromechanical factor are discussed in detail. The displacement, electric potential, and stress distributions along thickness of the graded layer are calculated and plotted. Numerical examples indicate that appropriate gradient distributing of the material properties make Love waves to propagate along the surface of the piezoelectric layer, or a bigger electromechanical coupling factor can be obtained, which is in favor of acquiring a better performance in surface acoustic wave (SAW) devices.     

Non-iterative bi-directional wave propagation method for treating reflections

We describe a bi-directional wave propagation method based on the split-step non-paraxial collocation method. The method is non-iterative and does not require any matrix inversion or diagonalization. The method is applied to reflections from waveguide discontinuities, terminations as well as multilayer reflecting structures. The results obtained compare very well with those obtained using computationally intensive methods such as the finite element method (FEM) and the finite-difference time-domain (FDTD) method.