Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media

Size: px
Start display at page:

Download "Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media"

Transcription

1 INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys PII: S Anisotropic reference edia and the possible linearized approxiations for phase velocities of qs waves in weakly anisotropic edia Lin-Ping Song and Manfred Koch Departent of Geohydraulics and Engineering Hydrology School of Civil Engineering University of Kassel Kurt-Wolters-Strasse 4 D-409 Kassel Gerany Received 5 June 00 Published 4 Noveber 00 Online at stacks.iop.org/jphysd/5/007 Abstract Using perturbation theory we propose possible first-order approxiations of qs-wave phase velocities for a weakly anisotropic ediu under anisotropic reference edia. The chosen reference edia have the features of non-degenerate qs waves. It is found that the usually arbitrarily chosen zeroth-order eigenvector basis associated with the degenerate shear waves in an isotropic reference ediu are none other than the shear-wave polarization vectors in one such anisotropic reference ediu. It follows that the qs-wave phase velocities are approxiated as linear functions of the elastic constants of a weakly anisotropic ediu. The expressions illustrate that the qs-wave speeds have the ost probably sensitive coefficients C C C 66 [C +C 66 C + C ] [C +C C + C ] [C +C C + C ] and those quantifying the deviation fro orthorhobic syetry. Nuerical exaples in an orthorhobic ediu show that the present sipler approxiations are sufficiently applicable to propagation directions being close to syetry planes or around syetry axes and achieve in these cases the sae order accuracy of approxiation as the quadratic equations derived fro degenerate perturbation theory. In light of the used anisotropic reference edia the validity of the approxiations presented should be liited to anisotropic solids that possess at least one plane of syetry.. Introduction In anisotropic edia there are three bulk odes that propagate with distinct phase velocities and the particle displaceents for these wave odes are generally neither along nor perpendicular to the wave noral. These odes are usually called qp and qs waves and becoe pure P and S waves only in isotropic edia. The phase velocities and the polarization vectors of these three waves governed by the well-known Christoffel equation are generally very coplicated functions of the wave noral and the elastic constant tensor. On syetry planes of orthorhobic or higher syetry systes the coplexity of these expressions is reduced to a certain extent. Yet fro a theoretical point of view as well as for practical applications in fields as diverse as seisology and non-destructive evaluation it is quite useful to siplify such coplicated functions under soe reasonable assuptions or conditions. This ay not only offer better physical insight into the dependence of wave propagation but also substantially facilitate the solution of forward and inverse probles. Motivated by these actual needs a nuber of researchers have ade efforts to seek approxiations of the wave phase velocities in anisotropic edia with higher syetry e.g. transverse isotropy and arbitrary syetry by eploying linearizing ethods [ 9]. Aong these ethods the perturbation approach is a powerful tool. Recently a review on the various linearized approxiations of phase velocities 00-77/0/007+08$ IOP Publishing Ltd Printed in the UK 007

2 L-P Song and M Koch was given by Song et al [0]. Farra [] further discussed the use of the higher order perturbations on the approxiations of phase velocities and polarizations of qp and qs waves. So far the use of perturbation theory in the area of wave propagation is usually carried out by assuing an isotropic reference ediu which has the siplest analytical solution. In this approach the anisotropy is naturally treated as a perturbation fro isotropy. It follows that the first-order approxiation of the qp phase velocity along an arbitrary direction of propagation is explicitly expressed as a linear function of the elastic constants of weakly anisotropic edia. It has been deonstrated that the qp approxiation works very well even for fairly large degrees of anisotropy [ 7]. For the qp wave Mensch and Farra [] developed a Hailtonian ray tracing approach for a weakly orthorhobic ediu using an ellipsoidal reference ediu whose principal axes are parallel to the principal axes of orthorhobic syetry. Ettrich et al [] presented an averaging technique to copute a bestfitting reference ellipsoidal ediu to an arbitrary anisotropic ediu. For qs waves the degenerate case occurs for reference isotropic edia. For the general syetry case the derived first-order approxiations for qs-wave phase velocities along an arbitrary direction are coplicated non-linear functions of the elastic constants. Song et al [0] suggested qs-wave phase velocity expressions which are linear functions of the elastic constants assuing that wave propagation has a slight deviation fro a syetry plane of the orthorhobic ediu. Mensch and Rasolofosaon [6] presented the approxiations of qs-wave phase velocities by a coordinate-rotational technique which also take the sae for as those fro the firstorder degenerate perturbation theory. Notably they used a special anisotropic reference ediu for their derivation. In particular this reference ediu has characteristics that are isotropic with respect to a qp wave but ay have two different shear waves. Therefore such characteristics naturally spur us to apply in a different way fro Mensch and Rasolofosaon [6] nondegenerate perturbation theory to qs waves in such kinds of reference edia. As a result first-order approxiations of qs-wave phase velocities are obtained that are now expressed in the for of linear functions of elastic constants of an anisotropic ediu as given below. In the following sections we first describe the fundaental equations being used in the developent of the perturbation theory. Second the choices of special reference anisotropic edia are discussed in detail. Then the explicit linear expansions of the phase velocities are presented for general and orthorhobic edia. Finally the accuracy and liitation of the approxiations are illustrated by nuerical exaples.. Perturbation expansions of Christoffel equations: non-degenerate case Elastic wave propagation in hoogeneous anisotropic edia is described by the Christoffel equations [4]: Ɣa = σ a. In this equation Ɣ is the syetric Christoffel tensor whose eleents are Ɣ = C kl n k n l ; C kl is the fourth-rank tensor of elastic stiffnesses of the ediu with at ost independent constants and n j the coponent of wave noral n; a is the eigenvector or unit polarization vector of the equation; σ = ρv is the corresponding eigenvalue where ρ is the density of the ediu and v is the phase velocity. The characteristic equation of is Ɣ σi =0 where I is an identity atrix. This equation being cubic in σ generally has three distinct roots corresponding to a qp and two qs waves with utually orthogonal polarization vectors. Closed for expressions for the three roots can be obtained for an arbitrary anisotropic ediu and these siplify for certain special cases [5]. Now let us consider perturbation of the Christoffel equation. Assuing that we have a reference ediu with the elastic constants kl then the original elastic constant atrix C kl of the considered ediu can be decoposed into [ 5] C kl = kl + εc kl where εc kl is regarded as a perturbation and ε is a sall paraeter used to develop the perturbation series [6]. The perturbation of εc kl of the reference ediu leads to the perturbed Christoffel tensor Ɣ = +εɣ. So the perturbed Christoffel equation is + εɣ a = σ a 4 where the eigenvalue and eigenvector have expansions in powers of ε σ = a = a 0 + εσ + εa + 5a + 5b where =. Substituting equations 5a and 5b into equation 4 and equating the coefficients of powers of ε we obtain for the non-degenerate case the first-order eigenvalue corrections and the first-order eigenvector corrections [5 ]: a σ = q = a0t Ɣ a 0 a 0T Ɣ a 0 q σ 0 σ q 0 6a a 0 q 6b where q = and the superscript T eans vector transpose. One observes that the eigenvalue and eigenvector corrections are dependent on the unperturbed eigenvectors defined in a reference ediu [6]. Here the purpose of presenting equation 6b is to show that the first-order eigenvalue correction equation 6a is iplicitly validated under the condition σ 0 σ q 0 i.e. there is no singularity in a direction of propagation. Moreover we would like to ephasize the constraints that not only the perturbation of the eigenvalue should be sall equation 5a but the perturbation of the polarization as well equation 5b if an anisotropic reference ediu is used. This is why the 008

3 Approxiations for phase velocities of qs waves use of an ellipsoidal anisotropic reference ediu should be interesting for directions close to a syetry plane as described below. In an isotropic reference ediu where the degenerate S-wave eigenvectors are undeterined the associated zeroth-order polarization vectors having to be constructed by the linear cobination of the two chosen vectors in a plane perpendicular to a wave propagation direction depend on the perturbed elastic paraeters and therefore the constraint on the polarization can be satisfied autoatically [].. Anisotropic reference ediu When using perturbation theory one has the freedo to choose a reference ediu. As a general rule the physical properties of the chosen reference ediu should be as close to those of the studied ediu as possible and at the sae tie a siple analytical solution should exist. Obviously a straightforward choice is to use an isotropic reference ediu and assue weak anisotropy of the ediu under study. However since an isotropic ediu has the two degenerate shear waves degenerate perturbation theory has to be used in such a case in order to develop the first-order approxiations for qs waves. As a result the associated approxiations then have to be expressed intricately as the square roots of a quadratic function and are thus highly nonlinear in the elastic constants of the ediu. Next we consider the other possible reference ediu having non-degenerate shear waves. One for of a reference ediu that has such desired properties was proposed by Mensch and Rasolofosaon [6]. Nevertheless in their developent the wave propagation characteristics were ignored for such a special reference ediu. Now we start with a reference ediu that possesses orthorhobic syetry. Using the Voigt atrix notation C αβ α β =...6; C αβ = C βα for the elastic constants C kl the possible reference elastic atrix is αβ = This atrix has seven independent elastic constants. Following Mensch and Rasolofosaon [6] the independent nuber in this atrix could be further reduced to four by iposing the conditions: = C0 C0 8a = C0 C0 8b = C0 C c Given the wave noral vector n the coponents of the corresponding Christoffel atrix are = C0 n + C0 66 n + C0 n = C0 66 n + C0 n + C0 n = C0 n + C0 n + C0 n = Ɣ0 = C0 = Ɣ0 = C0 = Ɣ0 = C0 66 n n n n n n. Now the reference ediu defined by the elastic atrix 7 and the conditions 8 is isotropic for the qp wave whose square phase velocity is in any direction and whose polarization vector is along the wave noral n. Denoting the polarization vectors of the two qs waves as ξ and η let us consider possible analytical solutions of the qs waves in the reference ediu. For the oent it suffices to write the following relations for the above elastic atrix: n in j = Ɣ0 ξ iξ j = η + C0 η + C0 66 η η iη j = ξ + C0 ξ + C0 66 ξ n iξ j = n iη j = 0 ξ iη j = ξ η + ξ η + 66 ξ η. 9 0 Note that the above derivations have used the cross-product relationships aong the three orthogonal eigenvectors n ξ and η. Fro 0 it is interesting to observe that the vectors n and ξ and n and η are orthogonal with respect to the reference Christoffel tensor. We call this property conjugate orthogonality by analogy with the conjugate concept in nuerical atheatics. To ake the two shear polarization vectors satisfy the conjugate orthogonal condition too we need to ipose new constraints together with conditions 8 on the elastic atrix. First let us set ξ iη j = ξ η + ξ η + 66 ξ η = 0. Then in connection with the vector orthogonal relation ξ η + ξ η + ξ η = 0 we have C0 ξ η + C0 66 ξ η = 0 a C0 ξ η + C0 66 ξ η = 0 b C0 66 ξ η + C0 66 ξ η = 0. c The relations iply the following four possible choices or cases: = C0 C0 66 ; where ξ = 0 and ξ η + ξ η = 0. = C0 66 C0 ; where ξ = 0 and ξ η + ξ η = 0. C0 = C0 66 ; where ξ = 0 and ξ η + ξ η = 0. 4 = C0 = C0 66 ; where ξ and η are arbitrarily chosen in a plane perpendicular to n. 009

4 L-P Song and M Koch Note that for the last 4 and the second possible choices the elastic atrix 7 under conditions 8 corresponds to a well-known isotropic reference ediu that will not be considered further here and the reference ediu of Mensch and Rasolofosaon [6] see their equation 7 respectively. For the first three cases we can now easily write down the eigenvalues and eigenvectors of their corresponding Christoffel equation : Case = n = n n n a = 66 n + n + C 0 n ξ = n n 0 n n b = η = n n n n n c n n where n = n + n /. Case = n = n n n a = n + n + C 0 n ξ = n 0 n n n b = η = n n n n n c n n where n = n + n /. Case = n = n n n 4a = n + n + C 0 n ξ = 0 n n n n 4b = η = n n n n n 4c n n where n = n + n /. For all three cases we can specify the eigenvalues associated with the P wave and the two S waves as P = S = ax = in S respectively by definition the subscript S represents the fast shear wave and S the slow shear wave. We observe that the solutions are siple and analytical. They represent one pure longitudinal wave and two pure nondegenerate shear waves. The polarization vector of one of the pure shear waves always lies in one of the syetry planes for exaple in the x y plane for case. In particular we can give a physical interpretation of the arbitrarily chosen unit vectors for the degenerate transverse waves coonly used in an isotropic reference ediu [7 0 ]. Indeed they are exactly the shear polarization vectors for case. Siilarly one ay choose one of the other two sets of shear eigenvectors for the cases and for the degenerate perturbation case. The solution characteristics of the three anisotropic reference edia suggest that a reasonable choice of the depend on that syetry plane which is the closest to the region of interest. Chadwick and Norris [7] studied the conditions under which the slowness surfaces of an anisotropic elastic aterial consist of aligned ellipsoids. Interestingly by adding an elastic constant constraint C = C +C to the conditions C discussed by Chadwick and Norris [7] their results for the eigenvalues S and eigenvectors P of the Christoffel equation could be siplified as those of case. The set of eigenvalues for each of the cases here describes the slowness surfaces with two spherical sheets and one ellipsoidal sheet. Following Chadwick and Norris [7] a ediu defined by the elastic atrix 7 and subjected to both the conditions 8 and one of the first three further possible choices on the reference shear oduli ay in a generalized sense be regarded as one class of ellipsoidal reference ediu. 4. Specification of the reference elastic constants In this section we would like to specify the reference elastic constants in ters of a weakly anisotropic ediu. Like Mensch and Rasolofosaon [6] we ay siply set the qp and qs wave oduli along a specially interested direction of wave propagation e.g. the vertical or the horizontal direction as reference constants. They ay be set as the averaging quantities for instance = C + C + C = C0 = C + C 66 = C 66 for case. More rigorously we can use Fedorov s ethod [8] to derive the reference elastic constants by iniizing a function F that is the su of the squared difference between the Christoffel atrices Ɣ and of the anisotropic and reference edia respectively: F = Ɣ = ij= ij= Ɣ + ij= ij= Ɣ 5 where the sybol denotes the averaging operation over all directions of n in order to obtain a relation that is solely dependent on the aterial properties. Through this iniization procedure and soe algebraic anipulations we obtain the reference elastic constants of the elastic atrix under the conditions 8 5 = C + C + C +C +C 66 +C +C +C +C 6a 70 = 4C +4C +4C 66 C +8C +8C +C +C 8C 6b 00

5 Approxiations for phase velocities of qs waves 70 = 4C +4C +4C 66 +8C C +8C +C 8C +C 6c = 4C +4C +4C 66 +8C +8C C 8C +C +C. 6d Equation 6a is just the longitudinal odulus already given in an isotropic reference ediu [8]. Then for case averaging any two e.g. and C0 of three equations 6b 6d would produce the other reference shear coponent required: 70 = 70 = 9C +9C +4C 66 +C +C +8C +C 8C 8C 7 where the sybol denotes the arithetic average. Of course equation 7 can also be obtained by setting = a priori in the iniization of function F 5. The derivation for the cases = 66 and = 66 is siilar. In addition the averaging of C0 and C0 66 in equations 6b 6d yields the shear odulus in an isotropic reference ediu [8]: 5 = C + C + C 66 + C + C + C C + C + C. 8 The reference ediu specified by such a iniization procedure is called an isotropic/anisotropic replaceent ediu IRM/ARM [6 8 8]. To quantify the intrinsic deviation of the considered ediu fro a reference ediu and taking at the sae tie the effects of wave propagation into account we ay use Ɣ 9 Ɣ where denotes the nor and Ɣ and are the Christoffel atrices of the considered ediu and a reference ediu respectively. Mensch and Rasolofosaon [6] directly used the elastic constant atrices to do this. 5. Expansions for the phase velocities of qs waves Using the eigenvalues in one of the anisotropic reference edia and their non-degeneracy as discussed in section we can directly write down the phase velocities to first-order: ρvqp = P + Ɣ n i n j = Ɣ n i n j 0a ρv = axɣ ξ i ξ j Ɣ η i η j 0b ρv qs = inɣ ξ i ξ j Ɣ η i η j 0c where the corrected shear-wave eigenvalues follow fro the relationships + Ɣ ξ i ξ j = Ɣ ξ i ξ j + Ɣ η i η j = Ɣ η i η j. For coparisons the first-order expressions for the qs-wave phase velocities under a degenerate isotropic reference ediu are also listed as [0 ] [ ρvqs = Ɣ ξ i ξ j + Ɣ η i η j ] ± Ɣ ξ i ξ j Ɣ η i η j +4Ɣ ξ i η j where + and before the square root ter correspond to and qs waves respectively. Substituting the expanded expressions for the coponent Ɣ and the corresponding unperturbed polarization vectors in one of the anisotropic reference edia e.g. in equation for case into equation 0 we get the following explicit firstorder approxiate forulae for the phase velocities in general anisotropy: ρvqp = { C n 4 + C n 4 + C n 4 +C +C 66 n n +C +C n n +C +C n } n + { 4C 6 n n +4C 5 n n +4C 6 n n +4C 4 n n +4C 5 n n +4C 4 n n +4C 4 +C 56 n n n +4C 5 +C 46 n n n +4C 6 +C 45 n n } n a ρv = axɣ ξ i ξ j Ɣ η i η j b ρv qs = inɣ ξ i ξ j Ɣ η i η j c where Ɣ ξ i ξ j = n { C66 n 4 + [ C + C C +C 66 ] n n + C n n + C n n +[ n n C 6 n n C 6 n n + C 46 n n + C 56 n n + n n n C 4 n + C 5 n C 45 n ]} Ɣ η i η j = C n + C n n + [ C +C 66 C + C ] n n n n + { [C + C C +C ]n +[C + C C +C ]n } n +4n n n [ n C6 n + C 6n C 45 + C 6 n ] +C45 n n n +n n n n n n [C4 +C 56 n + C 5 +C 46 n ]+n n n C5 n + C 4 n +n n n C 5 n + C 4n. It can be seen that the qp wave in equation a is ost sensitive to the 5 elastic constants or cobinations C C C C +C C +C 66 C +C and C 5 C 6 C 4 C 6 C 4 C 5 C 4 +C 56 C 5 +C 46 and C 6 + C 45 [ 6 8] the latter 9 paraeters of which involve the elastic constants quantifying the deviation fro orthorhobic syetry. The expressions in equations b c illustrate that the qs-wave speeds are ost probably sensitive to the elastic constants or cobinations C C C 66 [C +C 66 C + C ] [C +C C + C ] [C +C C + C ] and those elastic constants appearing in the latter 9 qp-wave sensitive paraeters. 0

6 L-P Song and M Koch For higher order syetry edia naely orthorhobic edia the equations b c for qs waves can be further reduced to ρvqs = axɣ ξ i ξ j Ɣ η i η j a ρvqs = inɣ ξ i ξ j Ɣ η i η j b where Ɣ ξ i ξ j = n [ C66 n 4 + [ C + C C +C 66 ] n n + C n n + C n ] n Ɣ η i η j = C n + C n n + [ C +C 66 C + C ] n n n n + {[ n C + C C +C ] n + [C + C C +C ] n }. Equations a b have been suggested by Song et al [0] who considered that the ter Ɣ ξ i η j in equation ay be sufficiently sall to be negligible when the wave vector deviates only slightly fro the x y syetry plane because its expansion has the leading factor n n n as shown below: Ɣ ξ i η j = n n n n { C n C n + C C + C +C 66 n n +[C +C C +C ]n }. 4 Note that the above expansions are based on the eigenvectors ξ and η in case. Under the reference ediu in case the approxiate expressions of qs-wave phase velocities in a weakly orthorhobic ediu will take the for ρvqs = axɣ ξ i ξ j Ɣ η i η j 5a ρvqs = inɣ ξ i ξ j Ɣ η i η j 5b where Ɣ ξ i ξ j = n [ C n 4 + [C + C C +C ] n n + C n n + C 66n ] n Ɣ η i η j = C 66 n + C n n + [ C +C C + C ] n n n n + { n [C + C C +C 66 ]n +[C + C C +C ]n }. As stated previously we ay choose the eigenvectors ξ and η in case as a zeroth-order eigenvector basis associated with the degenerate shear waves in an isotropic reference ediu. With this choice the expansion of the ter Ɣ ξ i η j in equation now is Ɣ ξ i η j = n n n n { C n C n + C 66 C + C +C n [ n + C +C C +C 66 ] n } 6 and also has the leading factor n n n. This again illustrates that the ter Ɣ ξ i η j in equation ay be sufficiently sall to be negligible when the wave vector is near the x z syetry plane. These two sets of approxiations and 5 for weakly orthorhobic edia assuing that the polarization of one of the qs waves is not too uch perturbed fro that of a pure ode for the reference edia of cases and show that the qs-wave phase velocities propagating near the x y plane or the x z plane will be ostly influenced by the different elastic constants or cobinations. 6. Nuerical exaples We test the validity of the linearized approxiations for the qs waves using the orthorhobic ediu of Farra [] which has the density-noralized elastic constants in units of k s : C = 0.8 C =. C = 8.5 C =.6C =.9 C 66 = 4.C =.C =.9 and C =.7. The corresponding IRM ediu to be used for coparison in the coputation has elastic constants = 0.04 C0 = 4.0 as calculated by equations 6a and 8. The corresponding ARM ediu coputed by equations 6a 7 and 6d has elastic constants = 0.04 = = = 4.4. The intrinsic deviations of the orthorhobic ediu fro the IRM and ARM edia are 9.5% and 9.% respectively. Here the intrinsic IRM deviation of 9.5% is saller that of 0.6% in Farra [] who used the elastic tensor for the calculation. First let us inspect the plot at the polar angle easured fro the vertical z-axis of 60. In figure a both the linear equation and the quadratic approxiations equation for fast quasi-shear wave are good with average relative errors of only 0.% and 0.05% respectively. One also notices that the quadratic equation is a better approxiation in the iddle aziuthal range than the linear one. Copared to the fast quasi-shear wave the approxiations for the qs wave slow quasi-shear wave have higher errors with.09% for the linear and 0.80% for the quadratic equation. The approxiate phase velocities of the qs wave are larger than the exact ones in the aziuthal range This iplies the need for soe aount of corrections fro the higher order perturbations see details in Farra []. Here we only focus on the issues of the first-order approxiation of qs-wave phase velocities. When the polar angle is gradually increased towards the x y plane the above described differences between the approxiate and the exact phase velocities are reduced significantly and the approxiate curves get closer and closer to the exact ones as shown in figures b d. In particular the sipler linear equations could yield satisfactory approxiations like the quadratic equations. For the linear approxiations the errors at a polar angle of 70 are now decreased to 0.% and 0.66% for and qs waves respectively; and for a polar angle of 80 the errors are further reduced to 0.07% and 0.0% for the two qs waves respectively. Even within the polar angle range the errors in the linear approxiations of waves are saller than those in the quadratic ones as depicted in table. It is observed that the direction factor n contained in the linear and quadratic equations has an iportant effect on the accuracy of their approxiations. When the wave propagates near the z-axis it ay cause nuerical instabilities. For exaple at a polar angle of 0 the linear equation equation cannot properly approxiate qs waves possibly due to the relatively saller values of n and the large deviation of the polarization vector fro that of a pure ode see equation. However using the linear equations 5a 5b which now contain the direction factor n and correspond to one reference shear polarization in the x z plane see equation the qs waves can be well approxiated like using the quadratic equations as illustrated in figure a. At this polar angle the average errors in 0

7 Approxiations for phase velocities of qs waves a Velocity k/s...0 Polar angle at 60 qs b Velocity k/s...0 Polar angle at 70 qs c Velocity k/s Aziuth deg...0 Polar angle at 80 d Velocity k/s Aziuth deg Polar angle at 90 qs qs Aziuth deg Aziuth deg Figure. The exact and approxiate linear equations a and b and quadratic phase velocity curves of qs waves in the orthorhoboic aterial for the aziuths of 0 90 at the fixed polar angles a 60 b 70 c 80 and d 90. Table. Averaged relative errors as a percentage between the exact and approxiate linear equations a and b and quadratic phase velocities of qs waves. For each polar angle the calculations are carried out at the aziuths of Polar angle - - -qs -qs the linear approxiations are 0.0% and 0.05% for and qs waves respectively being exactly the sae as those of the quadratic approxiations. For the polar angles of 0 and 40 the results in figures b c fro both the linear 5a 5b and the quadratic equations have siilar orders of approxiation accuracy; however they are less satisfactory due to the singularity between this polar angle interval [] than those at the polar angle of 0 shown in figure a. 7. Conclusions This paper discussing the choices of a set of ellipsoidal reference edia explored the possible first-order approxiations of qs-wave phase velocities for a weakly anisotropic ediu by non-degenerate perturbation theory. The derived approxiations appear as linear functions of the elastic constants of an anisotropic ediu. The corresponding expressions show that the qs-wave speeds are ost probably sensitive to the elastic coefficients C C C 66 [C +C 66 C + C ] [C +C C + C ] [C +C C + C ] and those constants quantifying the deviation fro orthorhobic syetry. In addition we could give a physical interpretation that the usually arbitrarily chosen zeroth-order eigenvector basis associated with the degenerate shear waves in an isotropic reference ediu are none other than the shear-wave polarization vectors in one of anisotropic reference edia discussed here. In for the present expressions are different fro the usual first-order degenerate approxiation of qs-wave phase velocities which are the square roots of a quadratic function of the ediu elastic constants. We see that the perturbed approxiations ay be siple or coplex depending on the reference ediu used. Of course the validity of the associated approxiations ight be liited to a certain range. The nuerical tests in the orthorhobic ediu deonstrate that the sipler linear approxiation is sufficiently applicable to propagation directions close to syetry planes or around syetry axes and does in these cases achieve the sae order of accuracy of approxiation as 0

8 L-P Song and M Koch a Velocity k/s Polar angle at 0 b Velocity k/s Polar angle at 0 qs.95 qs Aziuth deg Aziuth deg c Polar angle at 40 Velocity k/s.96.9 qs Aziuth deg Figure. The exact and approxiate linear equations 5a and 5b and quadratic phase velocity curves of qs waves in the orthorhoboic aterial for the aziuths of 0 90 at the fixed polar angles a 0 b 0 and c 40. the quadratic equations. However in view of the anisotropic reference edia used the valid use of these approxiations should be liited to aterials having at least one plane of syetry. Note that the quadratic approxiations can be used in a relatively large range of angles as they have the property of rotational invariance in a plane perpendicular to the wave noral [0]. In suary this study provides a theoretical and practical basis to possibly siplify inverse probles involving aterial characterization in non-destructive testing or seisology under qs-wave arrival easureent slightly deviating fro syetry planes/axes. Acknowledgents We are very grateful to V Farra for any valuable discussions and coents in the course of this work. We used her progra to calculate the exact phase velocities of qs waves in an orthorhobic ediu. We express our sincere thanks to A G Every I Psencik and the two anonyous reviewers whose coents and suggestions helped to iprove the anuscript and to clarify certain issues. LPS is thankful to T B Guttea for his help with soe issues regarding wording. References [] Backus G E 965 J. Geophys. Res [] Crapin S 98Wave Motion 4 [] Every A G and Sachse W 99 Ultrasonics 0 4 [4] Gassann F 964 PAGEOPH 58 6 [5] Jech J and Psencik I 989 Geophys. J. Int [6] Mensch T and Rasolofosaon P 997 Geophys. J. Int. 8 4 [7] Psencik I and Gajewski D 998 Geophys [8] Sayers C M 994 Geophys. J. Int [9] Thosen L 986 Geophysics [0] Song L-P Every A G and Wright C J. Phys. D: Appl. Phys [] Farra V 00 Geophys. J. Int [] Mensch T and Farra V 999 Geophys. J. Int. 8 [] Ettrich N Gajewski D and Kashtan B 00 Geophys. Prospecting 49 [4] Auld B A 990 Acoustic Fields and Waves in Solids vol Malabar FL: Krieger [5] Every A G 980 Phys. Rev. B 746 [6] Butkov E 968 Matheatical Phys. Reading MA: Addison-Wesley chapter 5 [7] Chadwick P and Norris A N 990 Q. J. Mech. Appl. Math [8] Fedorov F I 968 Theory of Elastic Waves in Crystals New York: Plenu 04

Data-Driven Imaging in Anisotropic Media

Data-Driven Imaging in Anisotropic Media 18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands

More information

Chapter 6 1-D Continuous Groups

Chapter 6 1-D Continuous Groups Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:

More information

Physics 139B Solutions to Homework Set 3 Fall 2009

Physics 139B Solutions to Homework Set 3 Fall 2009 Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about

More information

Feature Extraction Techniques

Feature Extraction Techniques Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

The Weierstrass Approximation Theorem

The Weierstrass Approximation Theorem 36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined

More information

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co

More information

A Simple Regression Problem

A Simple Regression Problem A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where

More information

Physics 215 Winter The Density Matrix

Physics 215 Winter The Density Matrix Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it

More information

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN

More information

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction

More information

3.3 Variational Characterization of Singular Values

3.3 Variational Characterization of Singular Values 3.3. Variational Characterization of Singular Values 61 3.3 Variational Characterization of Singular Values Since the singular values are square roots of the eigenvalues of the Heritian atrices A A and

More information

III.H Zeroth Order Hydrodynamics

III.H Zeroth Order Hydrodynamics III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t

More information

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization

Support Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering

More information

Pattern Recognition and Machine Learning. Artificial Neural networks

Pattern Recognition and Machine Learning. Artificial Neural networks Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lessons 7 20 Dec 2017 Outline Artificial Neural networks Notation...2 Introduction...3 Key Equations... 3 Artificial

More information

Electromagnetic fields modeling of power line communication (PLC)

Electromagnetic fields modeling of power line communication (PLC) Electroagnetic fields odeling of power line counication (PLC) Wei Weiqi UROP 3 School of Electrical and Electronic Engineering Nanyang echnological University E-ail: 4794486@ntu.edu.sg Keyword: power line

More information

In this chapter, we consider several graph-theoretic and probabilistic models

In this chapter, we consider several graph-theoretic and probabilistic models THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions

More information

Explicit Analytic Solution for an. Axisymmetric Stagnation Flow and. Heat Transfer on a Moving Plate

Explicit Analytic Solution for an. Axisymmetric Stagnation Flow and. Heat Transfer on a Moving Plate Int. J. Contep. Math. Sciences, Vol. 5,, no. 5, 699-7 Explicit Analytic Solution for an Axisyetric Stagnation Flow and Heat Transfer on a Moving Plate Haed Shahohaadi Mechanical Engineering Departent,

More information

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University

More information

2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL

2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL (Del I, teori; 1 p.) 1. In fracture echanics, the concept of energy release rate is iportant. Fro the fundaental energy balance of a case with possible crack growth, one usually derives the equation where

More information

Analytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order

Analytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order Geophysical Prospecting, 2006, 54, 651 661 Analytical solution for the electric potential in arbitrary anisotropic layered edia applying the set of Hankel transfors of integer order E. Pervago 1,A.Mousatov

More information

13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization

13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization 3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

Ph 20.3 Numerical Solution of Ordinary Differential Equations

Ph 20.3 Numerical Solution of Ordinary Differential Equations Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing

More information

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2

More information

Curious Bounds for Floor Function Sums

Curious Bounds for Floor Function Sums 1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International

More information

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics

More information

A GENERAL FORM FOR THE ELECTRIC FIELD LINES EQUATION CONCERNING AN AXIALLY SYMMETRIC CONTINUOUS CHARGE DISTRIBUTION

A GENERAL FORM FOR THE ELECTRIC FIELD LINES EQUATION CONCERNING AN AXIALLY SYMMETRIC CONTINUOUS CHARGE DISTRIBUTION A GENEAL FOM FO THE ELECTIC FIELD LINES EQUATION CONCENING AN AXIALLY SYMMETIC CONTINUOUS CHAGE DISTIBUTION BY MUGU B. ăuţ Abstract..By using an unexpected approach it results a general for for the electric

More information

(a) As a reminder, the classical definition of angular momentum is: l = r p

(a) As a reminder, the classical definition of angular momentum is: l = r p PHYSICS T8: Standard Model Midter Exa Solution Key (216) 1. [2 points] Short Answer ( points each) (a) As a reinder, the classical definition of angular oentu is: l r p Based on this, what are the units

More information

INNER CONSTRAINTS FOR A 3-D SURVEY NETWORK

INNER CONSTRAINTS FOR A 3-D SURVEY NETWORK eospatial Science INNER CONSRAINS FOR A 3-D SURVEY NEWORK hese notes follow closely the developent of inner constraint equations by Dr Willie an, Departent of Building, School of Design and Environent,

More information

ma x = -bv x + F rod.

ma x = -bv x + F rod. Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous

More information

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas

More information

MA304 Differential Geometry

MA304 Differential Geometry MA304 Differential Geoetry Hoework 4 solutions Spring 018 6% of the final ark 1. The paraeterised curve αt = t cosh t for t R is called the catenary. Find the curvature of αt. Solution. Fro hoework question

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion

More information

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Chapter 1: Basics of Vibrations for Simple Mechanical Systems Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,

More information

Sharp Time Data Tradeoffs for Linear Inverse Problems

Sharp Time Data Tradeoffs for Linear Inverse Problems Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used

More information

Ch 12: Variations on Backpropagation

Ch 12: Variations on Backpropagation Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith

More information

P032 3D Seismic Diffraction Modeling in Multilayered Media in Terms of Surface Integrals

P032 3D Seismic Diffraction Modeling in Multilayered Media in Terms of Surface Integrals P032 3D Seisic Diffraction Modeling in Multilayered Media in Ters of Surface Integrals A.M. Aizenberg (Institute of Geophysics SB RAS, M. Ayzenberg* (Norwegian University of Science & Technology, H.B.

More information

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths

More information

Homotopy Analysis Method for Nonlinear Jaulent-Miodek Equation

Homotopy Analysis Method for Nonlinear Jaulent-Miodek Equation ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol. 5, No.,, pp. 8-88 Hootopy Analysis Method for Nonlinear Jaulent-Miodek Equation J. Biazar, M. Eslai Departent of Matheatics, Faculty

More information

On Poset Merging. 1 Introduction. Peter Chen Guoli Ding Steve Seiden. Keywords: Merging, Partial Order, Lower Bounds. AMS Classification: 68W40

On Poset Merging. 1 Introduction. Peter Chen Guoli Ding Steve Seiden. Keywords: Merging, Partial Order, Lower Bounds. AMS Classification: 68W40 On Poset Merging Peter Chen Guoli Ding Steve Seiden Abstract We consider the follow poset erging proble: Let X and Y be two subsets of a partially ordered set S. Given coplete inforation about the ordering

More information

Pattern Recognition and Machine Learning. Artificial Neural networks

Pattern Recognition and Machine Learning. Artificial Neural networks Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2016 Lessons 7 14 Dec 2016 Outline Artificial Neural networks Notation...2 1. Introduction...3... 3 The Artificial

More information

Xiaoming Mao. Department of Physics and Astronomy, University of Pennsylvania. Collaborators: Tom Lubensky, Ning Xu, Anton Souslov, Andrea Liu

Xiaoming Mao. Department of Physics and Astronomy, University of Pennsylvania. Collaborators: Tom Lubensky, Ning Xu, Anton Souslov, Andrea Liu Xiaoing Mao Departent of Physics and Astronoy, University of Pennsylvania Collaborators: To Lubensky, Ning Xu, Anton Souslov, Andrea Liu Feb., 009 What is isostaticity? Isostatic systes are at the onset

More information

Intelligent Systems: Reasoning and Recognition. Artificial Neural Networks

Intelligent Systems: Reasoning and Recognition. Artificial Neural Networks Intelligent Systes: Reasoning and Recognition Jaes L. Crowley MOSIG M1 Winter Seester 2018 Lesson 7 1 March 2018 Outline Artificial Neural Networks Notation...2 Introduction...3 Key Equations... 3 Artificial

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of

More information

General Properties of Radiation Detectors Supplements

General Properties of Radiation Detectors Supplements Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent

More information

Interactive Markov Models of Evolutionary Algorithms

Interactive Markov Models of Evolutionary Algorithms Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary

More information

lecture 36: Linear Multistep Mehods: Zero Stability

lecture 36: Linear Multistep Mehods: Zero Stability 95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,

More information

RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS

RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di

More information

Supporting Information for Supression of Auger Processes in Confined Structures

Supporting Information for Supression of Auger Processes in Confined Structures Supporting Inforation for Supression of Auger Processes in Confined Structures George E. Cragg and Alexander. Efros Naval Research aboratory, Washington, DC 20375, USA 1 Solution of the Coupled, Two-band

More information

Principal Components Analysis

Principal Components Analysis Principal Coponents Analysis Cheng Li, Bingyu Wang Noveber 3, 204 What s PCA Principal coponent analysis (PCA) is a statistical procedure that uses an orthogonal transforation to convert a set of observations

More information

3.8 Three Types of Convergence

3.8 Three Types of Convergence 3.8 Three Types of Convergence 3.8 Three Types of Convergence 93 Suppose that we are given a sequence functions {f k } k N on a set X and another function f on X. What does it ean for f k to converge to

More information

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in

More information

Optimal nonlinear Bayesian experimental design: an application to amplitude versus offset experiments

Optimal nonlinear Bayesian experimental design: an application to amplitude versus offset experiments Geophys. J. Int. (23) 155, 411 421 Optial nonlinear Bayesian experiental design: an application to aplitude versus offset experients Jojanneke van den Berg, 1, Andrew Curtis 2,3 and Jeannot Trapert 1 1

More information

Mutual capacitor and its applications

Mutual capacitor and its applications Mutual capacitor and its applications Chun Li, Jason Li, Jieing Li CALSON Technologies, Toronto, Canada E-ail: calandli@yahoo.ca Published in The Journal of Engineering; Received on 27th October 2013;

More information

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents

More information

Analytical solution for nonlinear Gas Dynamic equation by Homotopy Analysis Method

Analytical solution for nonlinear Gas Dynamic equation by Homotopy Analysis Method Available at http://pvau.edu/aa Appl. Appl. Math. ISSN: 932-9466 Vol. 4, Issue (June 29) pp. 49 54 (Previously, Vol. 4, No. ) Applications and Applied Matheatics: An International Journal (AAM) Analytical

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 5 Lecture Oscillations (Chapter 6) What We Did Last Tie Analyzed the otion of a heavy top Reduced into -diensional proble of θ Qualitative behavior Precession + nutation Initial condition

More information

Moveout approximation for P waves in a homogeneous VTI medium

Moveout approximation for P waves in a homogeneous VTI medium Moveout approximation for P waves in a homogeneous VTI medium Véronique Farra 1 and Ivan Pšenčík 2 1 Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR 7154 CNRS,

More information

The Transactional Nature of Quantum Information

The Transactional Nature of Quantum Information The Transactional Nature of Quantu Inforation Subhash Kak Departent of Coputer Science Oklahoa State University Stillwater, OK 7478 ABSTRACT Inforation, in its counications sense, is a transactional property.

More information

IN A SENSE, every material is a composite, even if the

IN A SENSE, every material is a composite, even if the IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 47, NO. 11, NOVEMBER 1999 2075 Magnetis fro Conductors and Enhanced Nonlinear Phenoena J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart,

More information

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NIKOLAOS PAPATHANASIOU AND PANAYIOTIS PSARRAKOS Abstract. In this paper, we introduce the notions of weakly noral and noral atrix polynoials,

More information

Some Perspective. Forces and Newton s Laws

Some Perspective. Forces and Newton s Laws Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of

More information

The linear sampling method and the MUSIC algorithm

The linear sampling method and the MUSIC algorithm INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probles 17 (2001) 591 595 www.iop.org/journals/ip PII: S0266-5611(01)16989-3 The linear sapling ethod and the MUSIC algorith Margaret Cheney Departent

More information

EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS

EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS Matheatical and Coputational Applications, Vol. 5, No., pp. 96-07, 00. Association for Scientific Research EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI

More information

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION International Journal of Aerospace and Lightweight Structures Vol. 3, No. 1 (2013) 109 133 c Research Publishing Services DOI: 10.3850/S201042862013000550 A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS

More information

Numerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda

Numerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji

More information

Generalized r-modes of the Maclaurin spheroids

Generalized r-modes of the Maclaurin spheroids PHYSICAL REVIEW D, VOLUME 59, 044009 Generalized r-odes of the Maclaurin spheroids Lee Lindblo Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena, California 9115 Jaes R. Ipser

More information

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay

More information

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles An Inverse Interpolation Method Utilizing In-Flight Strain Measureents for Deterining Loads and Structural Response of Aerospace Vehicles S. Shkarayev and R. Krashantisa University of Arizona, Tucson,

More information

ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION

ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION The 4 th World Conference on Earthquake Engineering October -7, 8, Beijing, China ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION S. Li C.H. Zhai L.L. Xie Ph. D. Student, School of

More information

In this chapter we will study sound waves and concentrate on the following topics:

In this chapter we will study sound waves and concentrate on the following topics: Chapter 17 Waves II In this chapter we will study sound waves and concentrate on the following topics: Speed of sound waves Relation between displaceent and pressure aplitude Interference of sound waves

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

Golden ratio in a coupled-oscillator problem

Golden ratio in a coupled-oscillator problem IOP PUBLISHING Eur. J. Phys. 28 (2007) 897 902 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/28/5/013 Golden ratio in a coupled-oscillator proble Crystal M Mooran and John Eric Goff School of Sciences,

More information

A new type of lower bound for the largest eigenvalue of a symmetric matrix

A new type of lower bound for the largest eigenvalue of a symmetric matrix Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box

More information

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control An Extension to the Tactical Planning Model for a Job Shop: Continuous-Tie Control Chee Chong. Teo, Rohit Bhatnagar, and Stephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachusetts

More information

Department of Physics Preliminary Exam January 3 6, 2006

Department of Physics Preliminary Exam January 3 6, 2006 Departent of Physics Preliinary Exa January 3 6, 2006 Day 1: Classical Mechanics Tuesday, January 3, 2006 9:00 a.. 12:00 p.. Instructions: 1. Write the answer to each question on a separate sheet of paper.

More information

Time Evolution of Matter States

Time Evolution of Matter States Tie Evolution of Matter States W. M. Hetherington February 15, 1 The Tie-Evolution Operat The tie-evolution of a wavefunction is deterined by the effect of a tie evolution operat through the relation Ψ

More information

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

P016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion

P016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion P016 Toward Gauss-Newton and Exact Newton Optiization for Full Wavefor Inversion L. Métivier* ISTerre, R. Brossier ISTerre, J. Virieux ISTerre & S. Operto Géoazur SUMMARY Full Wavefor Inversion FWI applications

More information

Quantum Chemistry Exam 2 Take-home Solutions

Quantum Chemistry Exam 2 Take-home Solutions Cheistry 60 Fall 07 Dr Jean M Standard Nae KEY Quantu Cheistry Exa Take-hoe Solutions 5) (0 points) In this proble, the nonlinear variation ethod will be used to deterine an approxiate solution for the

More information

Least Squares Fitting of Data

Least Squares Fitting of Data Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a

More information

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS * IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS

More information

Comparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations

Comparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah

More information

Explicit solution of the polynomial least-squares approximation problem on Chebyshev extrema nodes

Explicit solution of the polynomial least-squares approximation problem on Chebyshev extrema nodes Explicit solution of the polynoial least-squares approxiation proble on Chebyshev extrea nodes Alfredo Eisinberg, Giuseppe Fedele Dipartiento di Elettronica Inforatica e Sisteistica, Università degli Studi

More information

Supplemental Material for Correlation between Length and Tilt of Lipid Tails

Supplemental Material for Correlation between Length and Tilt of Lipid Tails Suppleental Material for Correlation between Length and Tilt of Lipid Tails Ditry I. Kopelevich and John F. Nagle I. RESULTS FOR ALTERNATIVE DIRECTOR DEFINITIONS A. Alternative Director Definitions The

More information

Synchronization in large directed networks of coupled phase oscillators

Synchronization in large directed networks of coupled phase oscillators CHAOS 16, 015107 2005 Synchronization in large directed networks of coupled phase oscillators Juan G. Restrepo a Institute for Research in Electronics and Applied Physics, University of Maryland, College

More information

Bernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind

Bernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol., No., 6, pp.-9 Bernoulli Wavelet Based Nuerical Method for Solving Fredhol Integral Equations of the Second Kind S. C. Shiralashetti*,

More information

Biostatistics Department Technical Report

Biostatistics Department Technical Report Biostatistics Departent Technical Report BST006-00 Estiation of Prevalence by Pool Screening With Equal Sized Pools and a egative Binoial Sapling Model Charles R. Katholi, Ph.D. Eeritus Professor Departent

More information

Solving initial value problems by residual power series method

Solving initial value problems by residual power series method Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi

More information

Four-vector, Dirac spinor representation and Lorentz Transformations

Four-vector, Dirac spinor representation and Lorentz Transformations Available online at www.pelagiaresearchlibrary.co Advances in Applied Science Research, 2012, 3 (2):749-756 Four-vector, Dirac spinor representation and Lorentz Transforations S. B. Khasare 1, J. N. Rateke

More information

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique

More information

Dynamic analysis of frames with viscoelastic dampers: a comparison of damper models

Dynamic analysis of frames with viscoelastic dampers: a comparison of damper models Structural Engineering and Mechanics, Vol. 41, No. 1 (2012) 113-137 113 Dynaic analysis of fraes with viscoelastic dapers: a coparison of daper odels R. Lewandowski*, A. Bartkowiak a and H. Maciejewski

More information

Hysteresis model for magnetic materials using the Jiles-Atherton model

Hysteresis model for magnetic materials using the Jiles-Atherton model Hysteresis odel for agnetic aterials using the Jiles-Atherton odel Predrag Petrovic Technical faculty Svetog Save 65 32 Cacak, pegi@ei.yu Nebojsa itrovic Technical faculty Svetog Save 65 32 Cacak, itar@tfc.tfc.kg.ac.yu

More information

Phononic band gaps of elastic periodic structures: A homogenization theory study

Phononic band gaps of elastic periodic structures: A homogenization theory study PHYICA REVIEW B 75, 5414 27 Phononic band gaps of elastic periodic structures: A hoogenization theory study Ying-Hong iu, 1,3, * Chien C. Chang, 1,2,3, Ruey-in Chern, 2,3 and C. Chung Chang 1,3 1 Division

More information

Symbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm

Symbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat

More information

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes

More information