An accurately fast algorithm of calculating reflection/transmission coefficients

For the boundary between transversely isotropic media with a vertical axis of symmetry (VTI media), the interface between a liquid and a VTI medium, and the free-surface of an elastic half-space of a VTI medium, an accurately fast algorithm was presented for calculating reflection/transmission (R/T) coefficients. Specially, the case of post-critical angle incidence was considered. Although we only performed the numerical calculation for the models of the VTI media, the calculated results can be extended to the models of transversely isotropic media with a horizontal axis of rotation symmetry (HTI media). Compared to previous work, this algorithm can be used not only for the calculation of R/T coefficients of the boundary between ellipsoidally anisotropic media, but also for that between generally anisotropic media, and the speed and accuracy of this algorithm are faster and higher. According to the anisotropic parameters of some rocks given by the published literature, we performed the calculation of R/T coefficients by using this algorithm and analyzed the effect of the rock anisotropy on R/T coefficients. We used Snell’s law and the energy balance principle to perform verification for the calculated results. Supported by the Natural Science Foundation of Shaanxi Province, China (Grant No. 2007D15)     

TTI煤层群、相速度分析

以各向异性介质弹性矩阵为理论基础,结合TTI介质群、相速度相关公式,设计含有倾斜对称轴的横向各向同性煤层模型,通过模型试算,得到以下结论：随着各向异性参数ε、δ、γ变化,TTI型煤层相速度变化较为平缓,群速度变化较为剧烈,因此群速度决定着各向异性介质波场形状;在一定范围内,qP波的群、相速度随ε呈正相关变化;qSV波的群、相速度随δ呈负相关变化;qSH波的群、相速度只与各向异性参数γ有关,随γ呈正相关变化.     

Reflection and refraction of elastic wave at VTI-TTI media interface

We explore the physical phenomenon of acoustic waves induced at the interface between two different anisotropic rock media.Specifically,one medium is a transversely isotropic medium with a vertical axis of symmetry(VTI medium)and the other one is a transversely isotropic medium with a tilt axis of symmetry(TTI medium).By solving the Kelvin-Christoffel equation,an eighth-order polynomial is established for reflection and refraction angles,which is confirmed from SnelFs law.Three types of analytical expressions of the polarization coefficients of the induced waves are obtained corresponding to different incident angle regions.An effective algorithm has been developed for numerical analysis of the polarization coefficients.Applying characteristic anisotropic parameters reported in the literature,the influencing factors on reflection and refraction coefficients are analyzed,e.g.,the anisotropy,the tilt-angle of rock-layer,and the incident-angle.The calculated reflection and refraction coefficients have been rechecked for energy conservation.     

Method of projectors and the construction of Green's function of the wave equation

In the present article problems related to the propagation of waves in elastic anisotropic media with arbitrary types of symmetry are considered. Such problems are important for solid-body physics and for geophysics. An expansion of Green's function of the wave equation of the theory of elasticity is presented in the form of additive terms corresponding to the contributions of each of the three waves propagating in a solid body with designated anisotropic characteristics. An appropriate representation for the roots of the characteristic equation specifying the rate of wave propagation is presented. To illustrate the computation technique examples of certain types of media are considered. A representation is obtained for the static Green's function that does not require knowledge of the exact roots of the characteristic equation (assuming there is no degeneracy present).Moscow University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 4, pp. 101–111, April, 1995.     

Stiffness matrix invariants to validate the characterization of composite materials with ultrasonic methods

This paper presents a method of testing the ultrasonic measurements of the stiffness matrix, and the identification of the anisotropic behaviour, of composite materials. Some linear combinations of elastic constants are invariants for a rotation around an axis of symmetry. If the stacking sequence is the only parameter which changes in a set of long-fibre composites made of the superimposition of plies, the composites must own these invariants. So, PEEK-carbon fibre composite samples were constructed in this way to measure their elastic properties by immersion and contact ultrasonic methods, and to compare the results with predicted invariants. By changing the stacking sequence of plies three anisotropic models are tested: orthotropic, hexagonal and quadratic. Measurements of ultrasonic velocities in various planes of propagation permit the identification of the elastic constant and invariants. From the invariance of these linear combinations, the precision of the three-dimensional effective moduli can be estimated.     

Elastic moduli approximation of higher symmetry for the acoustical properties of an anisotropic material

The issue of how to define and determine an optimal acoustical fit to a set of anisotropic elastic constants is addressed. The optimal moduli are defined as those which minimize the mean-squared difference in the acoustical tensors between the given moduli and all possible moduli of a chosen higher material symmetry. The solution is shown to be identical to minimizing a Euclidean distance function, or equivalently, projecting the tensor of elastic stiffness onto the appropriate symmetry. This has implications for how to best select anisotropic constants to acoustically model complex materials.     

Anisotropic photonic crystals: Generalized plane wave method and dispersion symmetry properties

This article presents a generalized vector plane wave expansion method, applicable to isotropic and anisotropic periodic dielectric media of arbitrary geometry and dimension. The influence of anisotropic material orientation on the symmetry properties of photonic crystal dispersion surface is discussed. It is shown that the overall Brillouin zone symmetry is formed by the intersection of the photonic crystal lattice symmetry and the symmetry determined by the anisotropic material orientation. This work explains how to define the irreducible Brillouin zone of a two-dimensional anisotropic photonic crystal and demonstrates that doing it correctly allows one to avoid erroneous results, when calculating band gap diagrams of anisotropic photonic crystals. With the help of the methods presented, the possibility of controlling the band gaps of anisotropic photonic crystals by means of external electric field is shown.     

Deviatoric stress: a nuisance or a gold mine?

Both synchrotron radiation and deviatoric stress were once considered to be nuisances. Now synchrotron radiation is one of the most important tools available to scientists of all disciplines and deviatoric stress is one of the most useful aspects of x-ray diffraction at extreme conditions. Samples in high-pressure devices are under true hydrostatic pressure only when surrounded by a fluid, thus limiting true hydrostatic pressure studies at ambient temperatures to pressures below about 11?GPa. Elevated temperature is able to extend this limit but has rarely been used for this purpose. Instead, noble gases have been used as pressure media as their solids are especially soft. Deviatoric stress and resultant anisotropic elastic strain in solid samples and solid media have led to many subtle errors in determinations of elastic properties and crystal structures, especially in the days before it was realized that they could be measured and were potentially a valuable source of information. In recent years, measuring anisotropic elastic strain by x-ray diffraction has provided new insights into materials strength, elastic properties, crystal structures, mechanisms of phase transitions, slip systems, lattice preferred orientation, and, of course, ways to make corrections when deviatoric stress is indeed a nuisance.     

Elastic properties of spherically anisotropic piezoelectric composites

Effective elastic properties of spherically anisotropic piezoelectric composites, whose spherically anisotropic piezo-electric inclusions are embedded in an infinite non-piezoelectric matrix, are theoretically investigated. Analytical solutions for the elastic displacements and the electric potentials under a uniform external strain are derived exactly. Taking into account of the coupling effects of elasticity, permittivity and piezoelectricity, the formula is derived for estimating the effective elastic properties based on the average field theory in the dilute limit. An elastic response mechanism is revealed, in which the effective elastic properties increase as inclusion piezoelectric properties increase and inclusion dielectric properties decrease. Moreover, a piezoelectric response mechanism, of which the effective piezoelectric response vanishes due to the symmetry of spherically anisotropic composite, is also disclosed.     

Some elastic properties of solid nematics

The general approach to study the properties of the mechanical deformations of solid nematics, which are the macroscopic homogeneous elastic media having the rotational symmetry of the nematic liquid crystals is proposed. The stress tensor, the Young modulus and the Poisson ratios for the parallel and perpendicular homogeneous orientations of nematic molecules relative to the axis of external forces influence are obtained by the varying of the free energy of mechanical deformation. It is shown that these constants have the anisotropic character and the experiments for the direct measurement of five elasticity coefficients entering the free energy expression are proposed.     

Thin-film heteroepitaxy by the formation of the dilatation dipole ensemble

It is shown that in substantially anisotropic media, such as a crystal with cubic lattice symmetry, the analogous dilatation centers can strongly attract each other considerably decreasing the total elastic energy. Such attractive centers form stable objects of a new type, the elastic dilatation dipoles. By the example of heteroepitaxy of the film of silicon carbide SiC on a silicon substrate, the calculations of elastic energy of the system are performed. It is shown that the elastic energy can relax completely only due to the ensemble of dilatation dipoles. A new growth method of SiC on Si is suggested and implemented experimentally. In this method, the elastic energy relaxes due to the ensemble of elastic dilatation dipoles. It is shown experimentally that the heteroepitaxial SiC films grown do not contain cracks and lattice-misfit dislocations despite the tremendous difference in lattice parameters. Thick low-defect GaN and AlN layers are grown on the SiC/Si templates obtained. The laboratory model of the operating light-emitting diode is originally obtained based on these structures.     

Phase and group velocity considerations for dynamic modulus measurement in anisotropic media

In anisotropic media, the direction of energy propagation does not necessarily coincide with the wave normal, i.e. the energy flux vector does not coincide with the wave normal. Since, experimentally, one measures group velocity not phase velocity, one must therefore be careful in interpreting ultrasonic wave speed measurements in anisotropic media. This is of particular importance in elastic property reconstruction where acoustic velocity measurements are used as the basis for determining anisotropic material properties. In this work, the consequences of energy flux deviation from the wave normal are considered for typical experimental geometries. Particular attention is devoted to developing appropriate relationships between the phase velocity and ultrasonic transit time measurements, as these relations are most useful for elastic property reconstruction. In all the cases considered, it is shown that the phase velocity can be directly calculated from appropriate time delay measurements.     

Angular auto-correlation functions in molecular crystals and liquids: application to incoherent neutron scattering law

The effect of the local hindered molecular motion in the incoherent neutron scattering spectra in crystals and liquids is studied theoretically on the basis of the angular auto-correlation functions symmetrized on the dynamical point groups. The extended angular jump model simulates the motion. The microscopic properties of the matter, the time scale and the point symmetry of the molecule motion as well as the site symmetry of the molecule, are taken into account with the help of the dynamical variables of the model. The incoherent neutron scattering function is anisotropic in a monocrystalline sample. The scattering function consists of elastic plus quasi-elastic components. The shape of the quasi-elastic scattering spectrum is expressed by the sum of the weighted Lorentzian curves symmetrized with respect to the non-identical irreducible representations of the molecule motion point symmetry group. The elastic part of the intensity is increased by the contribution arising from the molecule motion of the identity representation symmetry. The scattering expression related to the molecule motion symmetry of a perfect cubic group is similar to the well-known expression derived for the scattering from the spherical top molecules exhibiting rotation diffusion.     

Spatial structure of an individual Mn acceptor in GaAs

The wave function of a hole bound to an individual Mn acceptor in GaAs is spatially mapped by scanning tunneling microscopy at room temperature and an anisotropic, crosslike shape is observed. The spatial structure is compared with that from an envelope-function, effective mass model and from a tight-binding model. This demonstrates that anisotropy arising from the cubic symmetry of the GaAs crystal produces the crosslike shape for the hole wave function. Thus the coupling between Mn dopants in GaMnAs mediated by such holes will be highly anisotropic.     

Specific propagation directions of acoustic waves for piezoelectric and nonpiezoelectric media

I reanalyze the problem of the existence of longitudinal normals inside the symmetry planes of piezoelectric crystals belonging to the symmetry class mm2. The equations determining components of longitudinal normals situated outside the symmetry planes for media of this symmetry are discussed. It is proven that nonpiezoelectric media of rhombic symmetry could have 4 or 8 distinct acoustic axes. Examples of nonpiezoelectric elastic media of monoclinic symmetry without acoustic axes are given. The method of determination of the components of acoustic axes for piezoelectric media of arbitrary symmetry is presented. With the help of this method, I discuss the problem of acoustic axes for piezoelectric media of the symmetry class mm2. The text was submitted by the author in English.     

Some Nonlinear Rock Behavior Effects

The nonlinear rock behavior effects observed in loading diagrams are analyzed which are usually ignored in conventional models of elastoplastic media. The initial deformation stage and unloading of rock samples are considered. The nonlinear behavior on these loading stages is interpreted from the viewpoint of partial closure of cracks initiated during deformation beyond the elastic limit or in earlier loading history. Phenomenological relations are derived to account for the discussed nonlinear effects in numerical modeling. The postcritical deformation stage corresponding to the stage of strain localization and main crack formation is studied. Corrections are made to provide a more accurate determination of model parameters.     

单向复合纤维材料弹性系数的超声测量

本文将单向复合纤维材料视作六角对称晶体，认为其弹性特征由五个独立的弹性系数来描述。利用Christoffel方程，我们得到这一各向异性介质中不同方向的声波速度与弹性系数的关系．我们还建立了一套计算机控制的实验测量系统，通过它对纤维板样品中的声速进行测量，反演得到了五个独立的弹性系数C11、C33、C13、C14和C12。理论计算得到的声速和实验测量值吻合较好，表明我们的方法是有实用意义的。     

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.     

Acoustic Features of the Elastic Anisotropic Properties of Rock Samples Along the Section of the Kola Superdeep Borehole (SD-3)

Acoustical Physics - The article compares the elastic anisotropic properties of samples recovered from the Kola Superdeep Borehole (SD-3 ). Samples were recovered from the Proterozoic and Archean...