Elements of Dynamic Ray tracing
|
|
- Jasper Jacobs
- 5 years ago
- Views:
Transcription
1 Elements of Dynamic Ray tracing ErSE360 Contents List of Figures 1 1 Introduction Coordinates systems 1 Cartesian coordinates 3 Ray coordinates γ 1 γ s 3 3 Ray-centered coordinate system q 1 q s 3 3 Coordinate transformations 6 31 Matrices 7 3 Jacobian 9 33 Hamiltonian in orthogonal curvilinear coordinates 10 4 Dynamic Ray Tracing 1 41 DRT in ray-centered Coordinates 1 4 DRT in Cartesian coordinates Paraxial traveltimes Algorithm for paraxial traveltime 15 5 Amplitude continuation Ray tubes 16 5 Geometrical spreading Coustics 18 6 Remarks 19 7 Exercise 19 References 19 List of Figures 1 Position vector in ray-centered coordinates 4 Difference between Frenet frame and ray-centered coordinates 7 3 Perpendicular planes and ray tube 16 4 Coustic points 18 1
2 ERSE360 1 Introduction Ray Tracing RT is an indispensable tool in applications which are concerned with the propagation of waves whether it is light or sound In seismic metho for geophysics we deal with sound waves propagating through the earth Apart from modeling these waves through wave equations elastic or acoustic there are many occasions when it is feasible to deal with traveltimes of the waves rather than the full wavefield This lea us to the evaluation of the rays which are lines connecting a source and any observation point such that they are always perpendicular to the wavefront Rays are a high frequency approximation to the wave equation In the last course equations for such an analysis which are eikonal and transport equations have been covered so we shall not into those derivations here Appropriate references are provided however Solution to the RT equations gives us kinematic part of the wave that is it tells us the path traced by a particular point on the wavefront and the time it took to reach from say point S to point R It gives us no information about the amplitude of the wave at the observation point For the amplitude at a particular point on the ray we solve the transport equation Solving each of these equations numerically for each point in the velocity model is expensive In these notes we shall consider a method of ray tracing called Dynamic Ray Tracing DRT which not only reduces the computation time for the traveltimes but also allows for the calculation of amplitudes along rays Albeit this is an approximate method The basic idea is to trace rays using conventional RT but doing so sparsely that is not covering every point in the model These are the so called central rays Then we calculate the traveltimes away in the so called paraxial vicinity from these central rays and this is where the approximation comes in The farther we are away from the central rays the less accurate are the traveltimes These are called the paraxial traveltimes Their accuracy depen of course on the complexity of the velocity model In this course we shall assume the medium to be smoothly varying and with no viscosity In that case the decay of the amplitudes of the waves is only due to geometrical spreading As we shall see an analytical solution to the transport equation involves a quantity J called the ray jacobian and the geometrical spreading is given by J Ray jacobian can also be expressed in terms of a matrix Q which can be calculated by DRT therefore DRT can give us the decay in amplitude due to geometrical spreading Paraxial traveltimes can be expressed in terms of matrices P and Q which are explicitly calculated in DRT Therefore with DRT they can also be evaluated Paraxial traveltimes can be evaluated explicitly using the Paraxial Ray Tracing PRT but we don t need to do it since DRT can take care of both traveltimes and amplitudes We shall also see that both PRT and DRT are identical as far as equations of motion are concerned The only difference is in the physical interpretation of the quantities being calculated In case of PRT the quantities calculated are the position and slowness components along rays in the paraxial vicinity of the central ray called paraxial rays In the case of DRT the calculated quantities are P and Q along the paraxial rays Coordinates systems In everyday life we face problems and our approach towar solving them depen on our frame of mind No matter what frame of mind we are in the nature of the problem is the same but our approach has a profound impact on the solution Same is true in physics We can choose any reference frame coordinates for a theory without altering the physics behind it but different coordinates will facilitate different features Therefore in our endeavor to study the DRT it would
3 ELEMENTS OF DYNAMIC RAY TRACING 3 be of prime importance to understand different coordinate systems which we can use to formulate the theory 1 Cartesian coordinates Let us start with the obvious Cartesian coordinates The RT equations with arclength s as a parameter are given by dx = V p dp 1 1 = V = 1 V where p = τ τ is the traveltime and V is the velocity Rays are the curves defined by vector x s which is function of the arc length Then from the first of equations 1 that the tangent to these curves is parallel to p therefore p gives the direction of travel of the ray First and second of equations 1 are a set of coupled first order ODEs To solve them we would need initial conditions that is three components of x 0 and three components of p 0 From eikonal equation we know that p = 1 V therefore only two of the the three components of p 0 are independent Once x and p have been solved for third one of equations 1 can be directly integrated to give the traveltime τ along the ray For spherical polar coordinates see [Červený 001] Ray coordinates γ 1 γ s Since we are tracing individual rays and in the end we want to be able to calculate the traveltimes in the neighborhoo of each of these rays it is important that we have a way of distinguishing one ray from the other Therefore we need some parameter or parameters which would stay constant along the ray Also they should be both necessary and sufficient to define a ray uniquely From the Cartesian coordinates above we can see that staying constant requirement is met by the initial conditions But since two different rays can have the same x 0 for a point source every ray will have the same starting point the three initial position component are not required We are left with the three components of initial direction p 0 but we saw above that only two of these are independent Therefore we conclude that in 3D only two independent parameters are required to uniquely identify a ray These two can be the angles specifying the initial direction of the ray also called the emergence angles It is obvious that indeed the initial direction is unique to every ray These two parameters for D only one is required can be any functions of the initial conditions but they have to be independent Let us call these two independent parameters γ 1 and γ Then consider a coordinate system with first two axes as γ 1 and γ and the third as s that is the arclength This is called the ray coordinate system In the ray coordinate system every point on s = constant plane would represent a single ray therefore is would indispensable if we have to calculate for example a rate of change of some quantity with respect to the rays themselves As we shall see this is precisely what we have to do in order to calculate geometrical spreading 3 Ray-centered coordinate system q 1 q s The two coordinate systems discussed above are such that they define every ray that it a point on every ray is a point in either one of those coordinates systems Now we shall consider a coordinate system which is attached to one particular ray Since orthogonal coordinate systems are easier to deal with it is natural that we look for coordinates which form and orthogonal system For a start we take the unit vector along one of the axes to be t s dx / dx = V p tangent to the ray itself and the other two unit vectors n and b to be perpendicular to t Coordinate along t is represented by s along n by q 1 and along b by q therefore a point is represented by q 1 q s One of the non-tangent unit vectors say n is defined as n s = dt / dt = 1 dts K K is the curvature of the ray therefore it is perpendicular to t s Then we define b as b s = t s n s The coordinate system
4 4 ERSE360 defined by t s n s and b s is called the Frenet frame We would also need their derivatives with respect to s which are a b c dt s = K s n s dn s = T s b s K s t s db s = T s n s These are the well known Frenet formulae [Thomas and Finney 1996] Now we want to calculate the distance from a point R on the ray to R in the Frenet frame According to Figure 1 the e R q q 1 R e 1 Plane perpendicular to the ray r q q s 1 Ray r00 s O 1 Figure 1 Position vector in ray-centered coordinates radius vector to the origin of the new coordinate system which is at R is r 0 0 s therefore the radius vector to R is 1 3 r q 1 q s = r 0 0 s + q 1 n + q b Length squared of a differential element of a vector in Cartesian coordinates is given as l = dx dx = dx + dy + dz In any other coordinate system say q 1 q q 3 square of the length dx becomes dx dx = dq i dx dx dq j dq i dq j g ij dq i dq j dq i is a unitary vector in the direction of q i and g ij is the so called metric tensor This means that if the new coordinate system is orthogonal we would have i = j in the new formula for square of the length Therefore in general when the new coordinate system is orthogonal square of the length becomes l = g ii dqi g 11dq1+g dq+g 33 dq3 that is the metric tensor does not have off-diagonal terms Now the differential of the vector in 1 For now consider the unit vectors e1 and e in Figure 1 to be n and b respectively The difference between these two is illustrated in Figure
5 ELEMENTS OF DYNAMIC RAY TRACING 5 equation 3 is 4 dr = dr dq 1 dq 1 {}}{ ndq 1 + dr dq dq {}}{ bdq + dr dq 3 dq 3 {[ }}{ t + dn q 1 + db ] q Length of dr can be calculated by taking the dot product of equation 4 by itself but let us first check whether our unitary vectors are orthogonal or not For equation 4 the cross terms dr dq 1 dr dq are zero because the vectors n and b are orthogonal But the off-diagonal terms involving q 3 are not zero for example dr dq 1 dr = T q which would in general be non-zero To see what went wrong we note that the third unitary vector is not in the direction of s This can be seen from the fact that it is a sum of t and other vectors which are not in the direction of t as is evident from equations Therefore when calculate the off-diagonal element of g ij with either i or j equal to three the dot product is not zero The problem is not t because it is orthogonal to dr n and dr dq b it is the other two vectors which are given by equations Therefore the Frenet frame is not a suitable choice for our purpose we require some other coordinates The solution is simple you might have guessed already We just need to make the derivatives of n and b parallel to t while keeping all three of them mutually perpendicular To distinguish from the Frenet frame notation we call our new vectors e 1 e and t Since we do not have to change t we keep the same notation We also set the new vectors according to 5 e 1 e = t e t = e 1 t e 1 = e This will ensure that we have a right-handed system of orthogonal coordinates Now to make de1 and de parallel to t we write de I s 6 = a I s t s I = 1 Where a s is some scalar function of s Throughout these notes indices in upper case for example I here run from 1 to To find a I we take the derivative of e t = 0 with respect to s This gives us de I t = e I dt By substituting this into the dot product of equation 6 with t we get a I = e I dt Then we use this value of a I in equation 6 again and get de I = e I dt 7 t We know that t = V p therefore dt = V dp + p dv Then since both e I are perpendicular to t hence also to p equation 7 becomes de I = V e I dp 8 pv Finally by using dp = 1 V = V V in equation 8 we get de I 9 = e I V p ei V = V A numerical solution of this differential equation will determine e I We need only to determine only of the e I the other one can be calculated by equation 5 p can be calculated by tracing a t dq 1
6 6 ERSE360 ray using conventional method this is the central ray The initial value for say e 1 at s = s 0 should be such that e 1 s 0 p s 0 = 0 and e 1 s 0 e 1 s 0 = 1 Now we shall write equation 4 in terms of our new coordinates that is 10 dr dq 3 dq 3 dr dr dq dq 1 1 dq dq {{}}{{}}{ [ }}{ dr = e 1 dq 1 + e dq + t + de 1 q 1 + de ] q [ e1 v e v = e 1 dq 1 + e dq q 1 + v v [ = e 1 dq 1 + e dq q 1 v + q v v q 1 v q q 1=0q =0 q ] t q 1=0q =0 The dot product of equation 10 does not have any off-diagonal term for g ij because all three vectors involved are orthogonal to each other Thus we have 11 where 1 dr dr = g 11 dq 1 + g dq + g 33 dq 3 = dq 1 + dq + h g ii h i and h 1 = 1 h = 1 h 3 = 1 + q I v v q I ] t I = 1 q 1=0q =0 h i are called the scale factors From this analysis we see that our new coordinate system is indeed orthogonal This new coordinate system with unit vectors e 1 e t and coordinates q 1 q s is called the ray-centered coordinate system There were two problems with the Frenet frame First one was apparent that it they were not orthogonal Second more subtle problem was that the unit vectors t n and b depended only on the local properties of the curve ray in our case This means that we would have had no way of distinguishing one ray from another if the tangent t happened to be the same for different rays at any point In the ray-centered coordinate system we overcome this problem by using the differential equation 9 to calculate the e I As we have to use the initial value p s 0 to determine e I and we know that the initial direction determines the ray uniquely the ray-centered coordinate system will depend on a particular ray rather than only the value of the tangent at a particular point The difference between Frenet frame and ray-centered coordinates as they move along the ray is illustrated in Figure As we saw above the unit vector t is common to both It is possible that at some point S on the ray e 1 e and n b are the same But as the ray moves they get rotated with respect to each other Note that e 1 e and n b would always lie in the same plane Of course we can apply simple rotation to get the ray centered coordinates from Frenet frame but we shall not go into that procedure here For details see [Červený 001] 3 Coordinate transformations Having studied the coordinate systems we would also need to study how to move back and forth between them As we shall see coordinate transformations play a very important role in not only simplifying the equations but also determining some important properties of the rays and their amplitudes We shall first look at the transformation matrices then the transformation of volume elements Jacobians of transformation and then the Hamiltonian systems in generalized coordinates
7 ELEMENTS OF DYNAMIC RAY TRACING 7 n e 1 e = n 1 t e b R t Ray S e = b Plenes perpendicular to the ray Figure Difference between Frenet frame and ray-centered coordinates 31 Matrices Throughout these notes the matrices with a hat example ˆQ would represent 3 3 matrices and without a hat example Q matrices matrices will represent the first two coordinates example q 1 and q of any system ˆQ Ray to Ray-centered Coordinates This transformation can be written as qi γ j By i and j can be regarded as rows and columns this transformation can be written as a matrix 13 ˆQ = q i γ j = q 1 γ 1 q 1 γ q 1 γ 3 q γ 1 q γ q γ 3 q 3 γ 1 q 3 γ q 3 γ 3 If we regard vectors as column matrices that is dˆγ = γ 1 γ γ 3 T and dˆq = q 1 q q 3 T equation 13 can be written as dˆq = ˆQdˆγ qi 14 or dq i = dˆγ j = Q ij dˆγ j γ j If we calculate ˆQ along a ray say R q 1 and q are zero This is because e I are perpendicular to t = vp which is tangent to the ray Moreover q 1 and q stay unchanged along the ray that is they are always zero Therefore we have q I γ 3 = 0 Then equation 13 becomes ˆQ R = Q 11 Q Q 1 Q 0 Q 31 Q 3 Q 33 Therefore we have R 16 det ˆQ R = Q 33 det Q R = q 3 γ 3 det Q R
8 8 ERSE360 where Q is the matrix which consists of first two rows and columns of rhs of equation 15 If we take the third coordinate to be the same in both ray and ray-centered coordinates we have det ˆQ R = det Q R Ĥ Ray-centered to Cartesian Coordinates This transformation can be written as dx i = H ij dq j = x i q j dq j or dˆx = Ĥdˆq For this matrix also we would need to determine it only along the ray R orthonormal along the ray that is The matrix Ĥ is 19 Ĥ 1 R = ĤT R det Ĥ R = 1 ˆQx Ray to Cartesian coordinates Again we shall consider the matrix only along the ray Like the other ones matrix ˆQ x is given by 0 dˆx = ˆQ x dˆγ or dx i = Q x ij dγ j = x i γ j dγ j By inserting the value of dˆq from equation 14 into equation 18 we get dˆx = Ĥ R ˆQ R dˆγ Then by comparing this with equation 0 we get 1 ˆQ x R = Ĥ R ˆQ R ˆP Ray coordinates to ray-centered components of slowness vector This transformation matrix transforms the ray coordinates to the components of slowness vector in the ray-centered coordinates The slowness vector p in ray-centered coordinates is given by p q = τ + τ + τ q 1 q q 3 where τ is the traveltime and the superscript q denotes the fact that the vector is in ray-centered coordinates The transformation along the ray R is written as dp q p q i τ 3 i = P ij dγ j = dγ j = dγ j or dˆp q = γ j q i γ ˆPdˆγ j R as 4 R ˆP x Ray to Cartesian components of slowness vector This matrix is defined along the ray R dp i = P x ij dγ j = pi γ j dγ j = R τ x i γ j dγ j or dˆp = ˆP x dˆγ R No superscript on the components p i means that they are in Cartesian coordinates The relation between ˆP x and ˆP is ˆP x R = Ĥ R ˆP R 5 M Second derivative of traveltime M is a matrix which can be represented in terms of P and Q as 6 M = τ x I x J = PQ 1 I J = 1 p q is not the physical slowness vector The components of the physical slowness vector p are related to the components of p q by p q h 1 i = p i
9 ELEMENTS OF DYNAMIC RAY TRACING 9 3 Jacobian After determining the transformation matrices we shall now turn our attention to the Jacobians of transformation Jacobians come across when we transform volume elements from one coordinate system to another Although details can be found in any text on calculus or classical differential geometry for example [Arfken and Weber 005] we shall briefly discuss the main features which will be of use to us A Jacobian is the determinant of the transformation For our purpose we look at the transformation from ray coordinates to Cartesian coordinates We already know the matrix for this transformation that is ˆQ x Then the Jacobian of transformation is det ˆQ x We set γ 3 u where u can be any monotonically increasing parameter along the ray like arclength s or traveltime τ Then we have x 1 x 1 x 1 J u = det ˆQ x γ 1 γ u x = x x 7 γ 1 γ u x 3 γ 1 γ u x 3 x 3 Since along the ray dx du = dx du = g 1 u dx dx and we know that = t equation 7 becomes x 1 x 1 γ 1 γ u t 1 u t 1 t t 3 J u = gu 1 x x γ 1 γ u t = g u 1 x 1 x x 3 8 u γ γ u γ u u x 3 x 3 x γ 1 γ u t 3 1 x x 3 γ 1 γ u u 1 γ u 1 u The last determinant in equation 8 is a vector triple product have 9 where Ω u = x γ 1 J u = 1 g u Ω u t = J g u x γ 1 x γ and J is called the ray Jacobian Since u x γ I x γ t Therefore we u are vectors in the directions of γ 1 and γ the direction of vector area Ω u is perpendicular to the plane u = constant It s magnitude Ω u is the area in u = constant plane Moreover the differential element of an area is given by 30 dω u = x x dγ 1 dγ = Ω u dγ 1 dγ γ 1 γ u Because wavefront is always perpendicular to the ray tangent we can get the projection of dω u on to the wavefront by 31 dω u t dω Then by substituting equation 30 in 9 we get 3 J = dω dγ 1 dγ Transformation of a volume element dxdydx is then given as [ x dxdydz = gu 1 x ] 33 t dγ 1 dγ du = gu 1 J dγ 1 dγ du γ 1 γ Before we leave this section note that 34 u J = g u J u
10 10 ERSE360 and especially J s = J 33 Hamiltonian in orthogonal curvilinear coordinates The Hamiltonian H x is a function of coordinates 3 say x i and momenta p i For our purpose it is the eikonal Since the eikonal equation remains constant along the rays H x x i p i is constant Rays can be parametrized by any monotonically increasing parameter u for example arclength s or traveltime τ H x x i p i will then be a function of u implicitly that is through x i and p i In the last course we derived RT systems for such parameters by calculating dh du Now we shall generalize that technique for curvilinear orthogonal coordinates say ξ 1 ξ ξ 3 These can be any coordinates as long as they are orthogonal to each other First we shall need to represent the eikonal equation in terms of ξ 1 ξ ξ 3 Since the eikonal equation is τ τ = V ξ 1 ξ ξ 3 we need to represent in terms of the curvilinear coordinates Let us write the expressions for length squared and in curvilinear coordinates which are 35a 35b = h 1dξ1 + h dξ + h 3dξ3 1 1 = e 1 + e h 1 ξ 1 h 1 + e 3 ξ h 3 ξ 3 where h i are the so called scale factors As we saw before h i = g ii remember that our curvilinear coordinates are orthogonal so g ij = 0 for i j Here e i represent the unitary vectors in the direction of ξ i Using equation 35b the eikonal in orthogonal curvilinear coordinates becomes 36 τ τ = 1 h τ h τ + 1 h τ 1 3 = 3 V ξ 1 ξ ξ 3 where τ i = τ ξ i Remember that the physical components of slowness vector in ξ i direction are 1 τ h i ξ i τ i is just an abbreviation Now we shall generalize the Hamiltonian As we know the eikonal can be expressed in different ways to be considered as a Hamiltonian and from these expression we can obtain RT equations for different parameters u [Slawinski 003] Here we shall write a general form for the Hamiltonian which is [ H ξ i τ i = 1 1 n h τ h τ + 1 ] n/ 37 h τ3 1 3 V n = 0 For practice you can prove that for n = 0 equation 37 becomes lim H ξ i τ i = 1 1 n=0 ln h τ h τ h τ3 + ln V 3 Before going any further with orthogonal curvilinear coordinates let us recall that for a Hamiltonian H x x i p i the ray RT system with a motonically increasing parameter u is given as [Slawinski 003 Červený 001] 39 dx i du = Hx p i dp Hx i du = x i du = p H x i p i 3 We will use the superscript x to denote Cartesian coordinates
11 ELEMENTS OF DYNAMIC RAY TRACING 11 If the Hamiltonian is written as H x x i p i = n 1 [ p k p k n/ 1/V n ] by inserting this into equations 39 the corresponding RT system is given by dx i du = p kp k n/ 1 dp i p i du = n x i V n For n = 0 it is 41 dx i = p kp k n/ 1 p i dp i = ln V x i du = p kp k n/ = V n = 1 As you can easily verify from equations 40 the case n = 1 correspon to u s and n = to u σ = τv For n = 0 equation 41 gives the RT system with u τ To derive a RT system for orthogonal curvilinear coordinates we just need to replace x i with ξ i p i with τ i note that the new components which replace p i are not the physical components of slowness vector and use equation 37 for H in equation 39 Then we get 4 dξ i du = τ An/ 1 i h i i du = 1 n ξ i n A n/ 1 τk h k V h k=1 k ξ i du = An/ = 1 V n If we are considering the Hamiltonian in orthogonal curvilinear coordinates for the case where n = 0 we have to use equation 38 in 39 then we get 43 dξ i = τ A 1 i h i i = ln V ξ i + A 1 3 τ k h k=1 k h k ξ i = 1 note that u has been changed to τ for n = 0 case So far we have dealt with RT systems in coordinates which were functions or some parameter u now we shall consider a situation where u is the third coordinates This is the case in ray and ray-centered coordinate systems where γ 3 u and q 3 u respectively An advantage of using such coordinate systems is that the conventional RT equations can be reduced from six to four To how this comes about consider equation 39 in Cartesian coordinates Set x 3 u so that the fist of these equations becomes dxi dx 3 = Hx p i When i = 3 in this equations we have dx3 dx 3 hence we conclude that H x must be of the form 44 H x x i p i = p 3 + H xr x i p I = 0 i = 1 3 and I = 1 = Hx p 3 = 1 and where H xr x i p I is the so called reduced Hamiltonian Since we can solve the eikonal equation for p 3 the reduced Hamiltonian can be evaluated by 45 H xr x i p I = p 3 To get the RT equations we note that the first of equations 39 has to be evaluated only for i = 1 because for i = 3 we have dx3 dx 3 = 1 which would be just 1 all along the ray The second one of equations 39 for i = 3 becomes dp3 dx 3 = Hx x 3 = 1 this is because of equation 44 Here also the derivative is just 1 all along the ray Therefore this one too has to be evaluated only for i = 1 Note however that p 3 is calculated by equation 45 Finally the third of equations 39 becomes H dx 3 = p 3 +p xr I p I this is because of Hx x 3 = 1 and equation 44 Thus we get our RT equations for a coordinate system where the third coordinate is a monotonously increasing parameter along
12 1 ERSE360 the ray These are 46 dx I dx 3 = HR p I dp I dx 3 = HxR x I H xr = p 3 + p I dx 3 p I As before to get the equation 46 in orthogonal curvilinear coordinates we first replace x i with ξ i and p i with τ i Equation for the reduced Hamiltonian corresponding to equation 44 then becomes 47 H ξ i τ i = τ 3 + H R ξ i τ I = 0 i = 1 3 and I = 1 and H R ξ i τ I is then given by 48 H R ξ i τ I = τ 3 Then the RT equations corresponding to equations 46 become 49 dξ I dξ 3 = HR τ I I dξ 3 = HR ξ I H R = τ 3 + τ I dξ 3 τ I Before leaving this section let us consider the general Hamiltonian from equation 37 for n = 1 in that case u ξ 3 s The reduced Hamiltonian becomes [ 1 H R ξ i τ I = τ 3 = h 3 V τ 1 h τ ] 1/ 50 1 h By inserting equation 50 in 49 we get the RT equations in orthogonal curvilinear coordinates that is 51 dξ I dξ 3 = 1 τ 3 h3 h I τ I I = τ 3 h 3 + h 3 dξ 3 h 3 ξ I τ 3 [ V h 1 τ1 h τ ξ I ξ I ξ I dξ 3 = h 3 V [ 1 V h 1 τ 1 + h τ ] ] 1/ = h 3 V τ Dynamic Ray Tracing Now we have gathered the tools necessary to understand DRT We shall look at the DRT first in ray-centered and then Cartesian coordinates 41 DRT in ray-centered Coordinates Calculation of the paraxial traveltimes are easier to do in the ray-centered coordinates Ray-centered coordinates are such that the third axis is the parameter which increases monotonously in the direction of the ray As we saw in the last section to derive a RT system for that case we need the reduced Hamiltonian Before going any further let us introduced some notation 5 v s = [V q 1 q s] q1=0q =0 v i s = [ V q 1 q s / q i ] q1=0q =0 v ij s = [ V q 1 q s / q i q j ]q 1=0q =0
13 ELEMENTS OF DYNAMIC RAY TRACING 13 Equations 5 are the velocity and its first and second partial derivatives along the ray that is q 1 q = 0 We already know the scale factors for ray-centered coordinates from equation 1 In our new notation h becomes 53 h = 1 + v 1 v I q I By using the scale factors for ray-centered coordinate system equation 1 and setting ξ 1 ξ ξ 3 to q 1 q s we get from equation 50 we get [ ] 1/ H R 1 54 q i τ I = τ 3 = h V q 1 q s τ 1 τ Similarly from equation 51 we get the RT equations in ray-centered coordinates that is 55a 55b dξ I = h τ I dξ 3 τ 3 [ I = h 1 dξ 3 τ 3 q I 1 V + τ 3 ] h h 3 q I Equations 55 are exact There is a problem however as far as numerical computation is concerned The rhs of equations 55b are zero along the ray that is when q 1 q = 0 but this is as a consequence of it being a difference of two terms which individually are not zero even when q 1 q = 0 Due to limited accuracy in numerical computations this difference might not always be zero as required theoretically This would introduce instability Since we are interested in the paraxial vicinity of the central ray we shall look for distances away from the central ray only up to the second order in q I that is in the so called quadratic vicinity Remember it was noted in the introduction that DRT is an approximate method this is where the approximation comes in First we do a Taylor series expansion of the velocity V about the central ray that is V = v + v I q I + 1 v IJq I q J With this expansion of V equation 53 for h and keeping only quadratic terms in q I our approximate reduced Hamiltonian becomes [ H R q i τ I = τ 3 = v v 1 v IJ q I q J 1 τ 1 + τ ] 56 Then from equations 49 the RT system becomes dq I = vτ I I = v v IJ q J and because τ i = h i p q i we get the final form for RT system in ray-centered coordinates that is 57 dq I = vpq I dp q I = v v IJ q J I J = 1 where the superscript q denotes the fact that the slowness components are in ray-centered coordinate system This is the Paraxial Ray Tracing PRT system in ray-centered coordinates We are a short way away from DRT now We need to calculate the rate of change of q I and p I given by the PRT equations 57 with respect to the rays themselves Since the rays can be identified by the first two axes of the ray coordinate system γ 1 γ We also note that since the third axis of ray coordinate system is monotonically increasing parameter which can be the arclength s This γ I and d means that to γ I this gives is d qi 58 γ K commute with each other Then we differentiate equations 57 with respect = v pq I γ K d p q I = v q J v IJ I J K = 1 γ K γ K
14 14 ERSE360 From equations 14 and 3 we know that q I γ K Q and pq I γ K P Then if we define another matrix V such that V IJ = v IJ = V q 1 q s / q I q J q 1=0q =0 equation 58 becomes 59 dq = vp dp = v VQ This the Dynamic Ray Tracing DRT system in ray-centered coordinates As we shall see the matrices P and Q calculated by DRT are essential in calculating paraxial traveltimes and geometrical spreading Geometrical spreading depen on the ray Jacobian which can be expressed in terms of Q according to equation 7 By differentiating the first of equations 59 with respect to s we can get a second order ODE for matrix Q only which is 60 v d Q v dq + VQ = 0 s Calculation of paraxial traveltimes in ray-centered coordinates involves a matrix M introduced in equation 6 A first order nonlinear ODE for M can be derived which is 61 dm + vm + v V = 0 4 DRT in Cartesian coordinates To get the DRT system in ray-centered coordinates we differentiated the PRT system also in ray-centered coordinates with respect to γ I along the ray Similarly to get DRT in Cartesian coordinates we differentiate equation 39 with respect to γ I H x depen on γ I implicitly that is through x i and p i For simplicity we shall take only one of the γ I say for I = 1 and write it as γ Then the derivative of the first of equations 39 x i γ for example is d x j du p x i x j γ + p x i p x j that slowness components are in Cartesian coordinates If we define 6 = H x xi Q x i γ u=const H x and we get the DRT in Cartesian coordinates which is p x j γ where the superscript x signifies the fact p x i γ u=cons P x i 63 where dq x i du dp x i du = Ax ij Qx j = Cx ij Qx j + B x ij P x j D x ij P x j 64 A x ij = H x p x i x j B x C x ij = H x x i x j ij = H x B x ij p x i p x j = H x x i p x j
15 Therefore for u σ we use H x = 1 coordinates that is 65 ELEMENTS OF DYNAMIC RAY TRACING 15 [ d dσ Qx i = P x i p x k px ] k 1/V and get the DRT in equations in Cartesian d dσ P x i = 1 x i x j 1 V Q x j 43 Paraxial traveltimes Now we shall look at traveltimes in the paraxial vicinity of the central ray These are called paraxial traveltimes As we can imagine ray-centered coordinates would give the simplest expression for paraxial traveltimes τ q 1 q s where q 1 and q are assumed to be small Of course on the central ray q 1 and q are zero Therefore the traveltime on the central ray is τ 0 0 s An equation for τ q 1 q s is obtained by its Taylor series expansion about a point say 0 0s on the ray only up to the second order in q I that is 66 τ q 1 q s = τ 0 0 s + qmq T q q 1 q T The matrix M can be written in terms of P and Q according to equation 6 It would be of practical importance to expand the traveltime in terms of Cartesian coordinates which is 67 τ x 1 x x 3 = τ x 10 x 0 x 30 + dxp x + dxm x dx T where dx x 1 x 10 x x 0 x 3 x 30 T and p x = τ The point x 10 x 0 x 30 is on the central ray and x 1 x x 3 in the paraxial vicinity Matrix can be evaluated in terms of P x and Q x with an equation analogous to 6 that is 68 M x = P x Q x 1 where matrices Q x and P x are given by equations 0 and 4 respectively P x and Q x can be calculated using equations Algorithm for paraxial traveltime Finally we are in a position to outline the steps required to calculate the paraxial traveltimes First in the ray-centered coordinates 1 Determine the central ray of the ray-centered coordinate system by solving RT equations 39 in Cartesian coordinates Determine the traveltimes and p along those rays 3 Solve for the vectors e I of the ray-centered coordinate system using equation 9 along the central-ray 4 Determine the propagation of matrices P and Q using equation 59 5 Calculate the paraxial traveltimes using equation 66 In Cartesian coordinates the steps are as follows: 1 Determine the central ray of the ray-centered coordinate system by solving RT equations 39 in Cartesian coordinates Determine the traveltimes 3 Determine the propagation of matrices P x and Q x using equation 65 4 Calculate the paraxial traveltimes using equation 67
16 16 ERSE360 D D C C Central Ray Paraxial Ray A=A N a Vectors N and t are the normals to the planes ABCD and A B C D respectively t B B dω τ τ+ b Ray tube dω Central Ray Figure 3 Perpendicular planes and ray tube 5 Amplitude continuation It was mentioned in the introduction that DRT can be used to calculated the geometrical spreading in this section we shall see how Geometrical spreading comes as a consequence of initial energy being distributed over the wavefront Since energy is conserved but the area of wavefront gets larger as the wave propagates same amount of energy gets spread over larger and larger area Energy at a point of the wavefront is related to the amplitude at that point therefore as energy per unit area is decreased so does the amplitude of the wave This suggests that geometrical spreading can be defined in terms of the area of wavefront To do that we should first understand the concept of ray tubes 51 Ray tubes Consider a central ray calculated by the standard RT equations in Cartesian coordinates We also take other rays in the paraxial vicinity of the central ray to make a bundle of rays as shown in Figure 3a This is called the ray tube The face of the ray tube represents area perpendicular to the direction of tangent to the central ray Therefore the face of the ray tube represents the plane of constant u where u is the monotonically increasing parameter along the central ray If u τ the face represents a plane tangent to the wavefront This is because wavefront is the surface of constant τ For u s or u σ the plane of constant u will in general be different from the wavefront and therefore will not be perpendicular to the ray-tangent direction p Figure 3b However we can take the dot product of the vector area element with the tangent to the central ray that is t = V p see equation 31 Since we are interested in the area of the wavefront only we shall take the u τ in the analysis below which would mean that the face of our ray tube will always be tangent to the wavefront As we shall see in order to calculate the amplitudes along a ray we need to evaluate the divergence of the slowness vector p that is p = τ We recall from elementary calculus that divergence is nothing but the net outward flux per unit volume Therefore we select a differential volume element from our ray tube Figure 3a and write the divergence of p as 69 1 p = lim D 0 D A surf p da surf where D = δxδyδz is the volume of the ray tube and A surf is the boundary of D which is the surface A surf is the vector surface area of the ray tube whose direction is always outward normal
17 ELEMENTS OF DYNAMIC RAY TRACING 17 to the surface and magnitude is the surface area Surface area of the ray tube can be divided into three parts areas of the two faces and area of the side of the ray tube Since the sides of the ray tube are the paraxial rays themselves the net outward flux of p from the sides is zero that is A sides p da sides = 0 Faces of our ray tube are tangent to the wavefront therefore the vector area elements of the face da face is nothing but dω t where dω is given by equation 3 Therefore equation 69 becomes 70 p = 1 gτ 1 Jdγ 1 dγ We know that du = g u therefore gτ 1 71 Therefore we have 7 p = 1 V J τ = p = 1 V J [ dω p t τ+ dω p t τ ] = V Also p t = V 1 Then equation 70 becomes [ JV 1 τ+ JV 1 ] τ d J = 1 V J d J V 5 Geometrical spreading Now let us see how the amplitude of the wave changes along a ray Obviously we will need the transport equation to do that The transport equation for an acoustic isotropic medium with variable density ρ is [Červený 001 ] 73 τ Ψ ρ + Ψ ρ τ = 0 where Ψ is the amplitude along the ray We know that p = τ = V 1 d Ψ t Also is nothing but the component of gradient in the direction tangent to the ray Therefore we can write d Ψ Ψ ρ = t ρ With these relations equation 73 can be written as d Ψ + V Ψ τ = 0 ρ ρ 74 By multiplying equation 73 throughout by Ψ ρ we get d Ψ Ψ 75 + V τ = 0 ρ ρ Ψ Let ρ Φ then the solution of equation 75 is 76 Φ s = Φ s 0 e R s s 0 V τs where τ depen on s implicitly that is through x y z equation 7 to get 77 Φ s = Φ s 0 e R s s 0 V J d J V ρ Then we substitute for τ from or by evaluating the integral in the exponent [ ] J s0 /V s 0 J s0 V s 78 Φ s = Φ s 0 exp ln = Φ s 0 J s /V s J s V s 0
18 18 ERSE360 Ψ Now we restore ρ Φ and get 79 Ψ s = ρ s Ψ s 0 ρ s0 Ψs Therefore the amplitude at a point s on a ray is ρs 80 J s0 V s J s V s 0 ρ s V s J s 0 Ψ s = Ψ s 0 ρ s 0 V s 0 J s Because the quantities which are functions of the initial point s 0 remain constant along the ray we can gather them into one constant that is C s 0 Ψ s 0 Js0 ρs 0V s 0 Then equation 80 becomes ρ s V s 81 Ψ s = C s 0 J s This is the equation of the amplitude variation along a ray and it is inversely proportional to the square root of ray Jacobian Therefore in analogy to the conventional geometrical spreading which is inversely proportional to the radial distance here we define it as J Finally we would like to express equation 81 in terms of quantities which can be directly evaluated by DRT To do that we note that from equations 7 and 9 J = J s = det Q ˆ x From equation 1 we have det ˆQ x s = det Ĥ s det ˆQ s along the ray But due to equations 17 and 19 det ˆQ x s = det Q s Therefore we can way that J = det Q s Thus our final expression for amplitude continuation along a ray is 8 ρ s V s Ψ s = C s 0 det Q s The central ray can be obtained by conventional RT Therefore ρ s and V s along the central ray are known Matrix Q s can be calculated along the ray by DRT equations 59 a Coustic point of the first order b Coustic point of the second order Figure 4 Coustic points 53 Coustics We do not intend to cover coustics in these notes but a few remarks would not be out of place In analyzing the ray tubes we assumed that the rays in the tube will never
19 ELEMENTS OF DYNAMIC RAY TRACING 19 cross each other This means that we assumed that the area of the face of our ray tube will never become zero But in reality the rays do cross each other When the rays intersect along a line the point on the rays which intersect that line is called coustic point of the first order Figure 4a and when they intersect each other at a point it is called the coustic point of second order Figure 4b It is clear from Figure 4 that as coustics are formed the area of our ray tube goes to zero that is Ω = Jdγ 1 dγ = 0 This means that at coustics J = det Q = 0 From equations 81 and 8 we can see that at coustics amplitude goes to infinity Therefore coustics are characterized by vanishing of the Jacobian For further reading see [Červený 001] 6 Remarks DRT is a much richer subject than what could be presented here Some of the tricks like rotation of the Frenet frame to fit with the ray-centered coordinate system were not discussed in these notes The purpose here was to give the most basic elements of DRT Most of the derivations in these notes have been taken from [Červený 001] 7 Exercise 1 Prove equation 19 Hint use equation 10 Prove equation 5 3 Prove equation 6 4 Derive equation 51 in ray-centered coordinates for n = 0 1 in equation 37 Is there any difference or not explain 5 Show that rhs of equations 55b are zero along the ray q 1 q = 0 6 Derive equation 61 7 Derive the DRT equations in Cartesian coordinates for u s and u τ 8 Why does equation 66 not contain the linear term in q References [Arfken and Weber 005] Arfken G B and Weber H J 005 Mathematical metho for physicists Sixth Edition [Slawinski 003] Slawinski M A 003 Seismic waves and rays in elastic media 34 [Thomas and Finney 1996] Thomas G B and Finney R L 1996 Calculus and analytic geometry Sixth Edition [Červený 001] Červený V 001 Seismic ray theory
Gaussian beams in inhomogeneous anisotropic layered structures
Gaussian beams in inhomogeneous anisotropic layered structures Vlastislav Červený 1 ) and Ivan Pšenčík 2 ) 1) Charles University, Faculty of Mathematics and Physics, Department of Geophysics, Ke Karlovu
More informationCovariant Formulation of Electrodynamics
Chapter 7. Covariant Formulation of Electrodynamics Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 11, and Rybicki and Lightman, Chap. 4. Starting with this chapter,
More informationSummary: Curvilinear Coordinates
Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 10 1 Summary: Curvilinear Coordinates 1. Summary of Integral Theorems 2. Generalized Coordinates 3. Cartesian Coordinates: Surfaces of Constant
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationChapter 11. Special Relativity
Chapter 11 Special Relativity Note: Please also consult the fifth) problem list associated with this chapter In this chapter, Latin indices are used for space coordinates only eg, i = 1,2,3, etc), while
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationSeismic ray theory. Summary. Introduction
Seismic ray theory Vlastislav Červený 1 ) and Ivan Pšenčík 2 ) 1 ) Charles University, Faculty of Mathematics and Physics, Department of Geophysics, Prague, Czech Republic. E-mail: vcerveny@seis.karlov.mff.cuni.cz.
More informationIn a uniform 3D medium, we have seen that the acoustic Green s function (propagator) is
Chapter Geometrical optics The material in this chapter is not needed for SAR or CT, but it is foundational for seismic imaging. For simplicity, in this chapter we study the variable-wave speed wave equation
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationIII. TRANSFORMATION RELATIONS
III. TRANSFORMATION RELATIONS The transformation relations from cartesian coordinates to a general curvilinear system are developed here using certain concepts from differential geometry and tensor analysis,
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationFrequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix
Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Ludě Klimeš Department of Geophysics, Faculty of Mathematics Physics, Charles University, Ke Karlovu 3, 121 16 Praha
More informationPHYS Handout 6
PHYS 060 Handout 6 Handout Contents Golden Equations for Lectures 8 to Answers to examples on Handout 5 Tricks of the Quantum trade (revision hints) Golden Equations (Lectures 8 to ) ψ Â φ ψ (x)âφ(x)dxn
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationLecture: Lorentz Invariant Dynamics
Chapter 5 Lecture: Lorentz Invariant Dynamics In the preceding chapter we introduced the Minkowski metric and covariance with respect to Lorentz transformations between inertial systems. This was shown
More informationChapter 0 of Calculus ++, Differential calculus with several variables
Chapter of Calculus ++, Differential calculus with several variables Background material by Eric A Carlen Professor of Mathematics Georgia Tech Spring 6 c 6 by the author, all rights reserved - Table of
More information4 Relativistic kinematics
4 Relativistic kinematics In astrophysics, we are often dealing with relativistic particles that are being accelerated by electric or magnetic forces. This produces radiation, typically in the form of
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationReview of Vector Analysis in Cartesian Coordinates
Review of Vector Analysis in Cartesian Coordinates 1 Scalar: A quantity that has magnitude, but no direction. Examples are mass, temperature, pressure, time, distance, and real numbers. Scalars are usually
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationHOMEWORK 2 SOLUTIONS
HOMEWORK SOLUTIONS MA11: ADVANCED CALCULUS, HILARY 17 (1) Find parametric equations for the tangent line of the graph of r(t) = (t, t + 1, /t) when t = 1. Solution: A point on this line is r(1) = (1,,
More informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationChapter 1. Vector Algebra and Vector Space
1. Vector Algebra 1.1. Scalars and vectors Chapter 1. Vector Algebra and Vector Space The simplest kind of physical quantity is one that can be completely specified by its magnitude, a single number, together
More informatione 2 = e 1 = e 3 = v 1 (v 2 v 3 ) = det(v 1, v 2, v 3 ).
3. Frames In 3D space, a sequence of 3 linearly independent vectors v 1, v 2, v 3 is called a frame, since it gives a coordinate system (a frame of reference). Any vector v can be written as a linear combination
More informationDistances, volumes, and integration
Distances, volumes, and integration Scribe: Aric Bartle 1 Local Shape of a Surface A question that we may ask ourselves is what significance does the second fundamental form play in the geometric characteristics
More informationLaplace equation in polar coordinates
Laplace equation in polar coordinates The Laplace equation is given by 2 F 2 + 2 F 2 = 0 We have x = r cos θ, y = r sin θ, and also r 2 = x 2 + y 2, tan θ = y/x We have for the partials with respect to
More informationMechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Mechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras Lecture 08 Vectors in a Plane, Scalars & Pseudoscalers Let us continue today with
More informationEngineering Mechanics Prof. U. S. Dixit Department of Mechanical Engineering Indian Institute of Technology, Guwahati Kinematics
Engineering Mechanics Prof. U. S. Dixit Department of Mechanical Engineering Indian Institute of Technology, Guwahati Kinematics Module 10 - Lecture 24 Kinematics of a particle moving on a curve Today,
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationAN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES
AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES Overview Throughout the first few weeks of the semester, we have studied vector calculus using almost exclusively the familiar Cartesian x,y,z coordinate
More informationPhysics 209 Fall 2002 Notes 5 Thomas Precession
Physics 209 Fall 2002 Notes 5 Thomas Precession Jackson s discussion of Thomas precession is based on Thomas s original treatment, and on the later paper by Bargmann, Michel, and Telegdi. The alternative
More informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationMATH 423/ Note that the algebraic operations on the right hand side are vector subtraction and scalar multiplication.
MATH 423/673 1 Curves Definition: The velocity vector of a curve α : I R 3 at time t is the tangent vector to R 3 at α(t), defined by α (t) T α(t) R 3 α α(t + h) α(t) (t) := lim h 0 h Note that the algebraic
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More informationCovariant Formulation of Electrodynamics
Chapter 7. Covariant Formulation of Electrodynamics Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 11, and Rybicki and Lightman, Chap. 4. Starting with this chapter,
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationInverse Kinematics. Mike Bailey. Oregon State University. Inverse Kinematics
Inverse Kinematics Mike Bailey mjb@cs.oregonstate.edu inversekinematics.pptx Inverse Kinematics Forward Kinematics solves the problem if I know the link transformation parameters, where are the links?.
More informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (SPHERICAL POLAR COORDINATES) Question 1 a) Determine with the aid of a diagram an expression for the volume element in r, θ, ϕ. spherical polar coordinates, ( ) [You may not
More informationSuperpositions of Gaussian beams and column Gaussian packets in heterogeneous anisotropic media
Superpositions of Gaussian beams and column Gaussian packets in heterogeneous anisotropic media Luděk Klimeš Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu,
More information5.3 Surface Theory with Differential Forms
5.3 Surface Theory with Differential Forms 1 Differential forms on R n, Click here to see more details Differential forms provide an approach to multivariable calculus (Click here to see more detailes)
More information(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3
Math 127 Introduction and Review (1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3 MATH 127 Introduction to Calculus III
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationMatrices and Deformation
ES 111 Mathematical Methods in the Earth Sciences Matrices and Deformation Lecture Outline 13 - Thurs 9th Nov 2017 Strain Ellipse and Eigenvectors One way of thinking about a matrix is that it operates
More informationAmplitude calculations for 3-D Gaussian beam migration using complex-valued traveltimes
CWP-648 Amplitude calculations for 3-D Gaussian beam migration using complex-valued traveltimes Norman Bleistein 1 & S. H. Gray 2 1 Center for Wave Phenomena, Dept. of Geophysics, Colo. Sch. of Mines,
More informationKinematics (2) - Motion in Three Dimensions
Kinematics (2) - Motion in Three Dimensions 1. Introduction Kinematics is a branch of mechanics which describes the motion of objects without consideration of the circumstances leading to the motion. 2.
More informationAstr 2320 Tues. May 2, 2017 Today s Topics Chapter 23: Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s
Astr 0 Tues. May, 07 Today s Topics Chapter : Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s Field Equations The Primeval Fireball Standard Big Bang Model Chapter
More informationSOME PROBLEMS YOU SHOULD BE ABLE TO DO
OME PROBLEM YOU HOULD BE ABLE TO DO I ve attempted to make a list of the main calculations you should be ready for on the exam, and included a handful of the more important formulas. There are no examples
More information2.5 Sound waves in an inhomogeneous, timedependent
.5. SOUND WAVES IN AN INHOMOGENEOUS, TIME-DEPENDENT MEDIUM49.5 Sound waves in an inhomogeneous, timedependent medium So far, we have only dealt with cases where c was constant. This, however, is usually
More informationLecture 10 - Moment of Inertia
Lecture 10 - oment of Inertia A Puzzle... Question For any object, there are typically many ways to calculate the moment of inertia I = r 2 dm, usually by doing the integration by considering different
More informationexample consider flow of water in a pipe. At each point in the pipe, the water molecule has a velocity
Module 1: A Crash Course in Vectors Lecture 1: Scalar and Vector Fields Objectives In this lecture you will learn the following Learn about the concept of field Know the difference between a scalar field
More informationClassical differential geometry of two-dimensional surfaces
Classical differential geometry of two-dimensional surfaces 1 Basic definitions This section gives an overview of the basic notions of differential geometry for twodimensional surfaces. It follows mainly
More informationMo, 12/03: Review Tu, 12/04: 9:40-10:30, AEB 340, study session
Math 2210-1 Notes of 12/03/2018 Math 2210-1 Fall 2018 Review Remaining Events Fr, 11/30: Starting Review. No study session today. Mo, 12/03: Review Tu, 12/04: 9:40-10:30, AEB 340, study session We, 12/05:
More informationLecture D4 - Intrinsic Coordinates
J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 Lecture D4 - Intrinsic Coordinates In lecture D2 we introduced the position, velocity and acceleration vectors and referred them to a fixed cartesian coordinate
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationManifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds
MA 755 Fall 05. Notes #1. I. Kogan. Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds Definition 1 An n-dimensional C k -differentiable manifold
More informationREVIEW - Vectors. Vectors. Vector Algebra. Multiplication by a scalar
J. Peraire Dynamics 16.07 Fall 2004 Version 1.1 REVIEW - Vectors By using vectors and defining appropriate operations between them, physical laws can often be written in a simple form. Since we will making
More informationVectors and Fields. Vectors versus scalars
C H A P T E R 1 Vectors and Fields Electromagnetics deals with the study of electric and magnetic fields. It is at once apparent that we need to familiarize ourselves with the concept of a field, and in
More informationPhysics 6303 Lecture 3 August 27, 2018
Physics 6303 Lecture 3 August 27, 208 LAST TIME: Vector operators, divergence, curl, examples of line integrals and surface integrals, divergence theorem, Stokes theorem, index notation, Kronecker delta,
More informationMathematical Notation Math Calculus & Analytic Geometry III
Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and should be emailed to the instructor at james@richland.edu.
More informationLecture 13 The Fundamental Forms of a Surface
Lecture 13 The Fundamental Forms of a Surface In the following we denote by F : O R 3 a parametric surface in R 3, F(u, v) = (x(u, v), y(u, v), z(u, v)). We denote partial derivatives with respect to the
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationSOLUTIONS TO EXERCISES FOR MATHEMATICS 133 Part 1. I. Topics from linear algebra
SOLUTIONS TO EXERCISES FOR MATHEMATICS 133 Part 1 Winter 2009 I. Topics from linear algebra I.0 : Background 1. Suppose that {x, y} is linearly dependent. Then there are scalars a, b which are not both
More information2t t dt.. So the distance is (t2 +6) 3/2
Math 8, Solutions to Review for the Final Exam Question : The distance is 5 t t + dt To work that out, integrate by parts with u t +, so that t dt du The integral is t t + dt u du u 3/ (t +) 3/ So the
More information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. V9. Surface Integrals Surface
More informationII. Unit Speed Curves
The Geometry of Curves, Part I Rob Donnelly From Murray State University s Calculus III, Fall 2001 note: This material supplements Sections 13.3 and 13.4 of the text Calculus with Early Transcendentals,
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More information2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.
Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and
More informationCourse Outline. 2. Vectors in V 3.
1. Vectors in V 2. Course Outline a. Vectors and scalars. The magnitude and direction of a vector. The zero vector. b. Graphical vector algebra. c. Vectors in component form. Vector algebra with components.
More information(You may need to make a sin / cos-type trigonometric substitution.) Solution.
MTHE 7 Problem Set Solutions. As a reminder, a torus with radii a and b is the surface of revolution of the circle (x b) + z = a in the xz-plane about the z-axis (a and b are positive real numbers, with
More informationThe integral test and estimates of sums
The integral test Suppose f is a continuous, positive, decreasing function on [, ) and let a n = f (n). Then the series n= a n is convergent if and only if the improper integral f (x)dx is convergent.
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More informationContravariant and Covariant as Transforms
Contravariant and Covariant as Transforms There is a lot more behind the concepts of contravariant and covariant tensors (of any rank) than the fact that their basis vectors are mutually orthogonal to
More informationn=0 ( 1)n /(n + 1) converges, but not
Math 07H Topics for the third exam (and beyond) (Technically, everything covered on the first two exams plus...) Absolute convergence and alternating series A series a n converges absolutely if a n converges.
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationAn Analytic Basis for Heat Conduction in a Fuel Rods
An Analytic Basis for Heat Conduction in a Fuel Rods Dan Shields December 15, 211 Introduction It is important in nuclear science and engineering to understand heat flow in and around fuel in reactors.
More informationKinematics of fluid motion
Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx
More informationLecture 6.1 Work and Energy During previous lectures we have considered many examples, which can be solved using Newtonian approach, in particular,
Lecture 6. Work and Energy During previous lectures we have considered many examples, which can be solved using Newtonian approach, in particular, Newton's second law. However, this is not always the most
More informationPhysics 106a/196a Problem Set 7 Due Dec 2, 2005
Physics 06a/96a Problem Set 7 Due Dec, 005 Version 3, Nov 7, 005 In this set we finish up the SHO and study coupled oscillations/normal modes and waves. Problems,, and 3 are for 06a students only, 4, 5,
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationWilliam P. Thurston. The Geometry and Topology of Three-Manifolds
William P. Thurston The Geometry and Topology of Three-Manifolds Electronic version 1.1 - March 00 http://www.msri.org/publications/books/gt3m/ This is an electronic edition of the 1980 notes distributed
More informationContents. Motivation. 1 di 7 23/03/ :41
1 di 7 23/03/2015 09:41 From Wikipedia, the free encyclopedia In mathematics, orthogonal coordinates are defined as a set of d coordinates q = (q 1, q 2,..., q d ) in which the coordinate surfaces all
More informationPhysics 6303 Lecture 2 August 22, 2018
Physics 6303 Lecture 2 August 22, 2018 LAST TIME: Coordinate system construction, covariant and contravariant vector components, basics vector review, gradient, divergence, curl, and Laplacian operators
More informationTwo And Three Dimensional Regions From Homothetic Motions
Applied Mathematics E-Notes, 1(1), 86-93 c ISSN 167-51 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Two And Three Dimensional Regions From Homothetic Motions Mustafa Düldül Received
More informationPhysics/Astronomy 226, Problem set 4, Due 2/10 Solutions. Solution: Our vectors consist of components and basis vectors:
Physics/Astronomy 226, Problem set 4, Due 2/10 Solutions Reading: Carroll, Ch. 3 1. Derive the explicit expression for the components of the commutator (a.k.a. Lie bracket): [X, Y ] u = X λ λ Y µ Y λ λ
More informationPhys460.nb Back to our example. on the same quantum state. i.e., if we have initial condition (5.241) ψ(t = 0) = χ n (t = 0)
Phys46.nb 89 on the same quantum state. i.e., if we have initial condition ψ(t ) χ n (t ) (5.41) then at later time ψ(t) e i ϕ(t) χ n (t) (5.4) This phase ϕ contains two parts ϕ(t) - E n(t) t + ϕ B (t)
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationChapter 7. Kinematics. 7.1 Tensor fields
Chapter 7 Kinematics 7.1 Tensor fields In fluid mechanics, the fluid flow is described in terms of vector fields or tensor fields such as velocity, stress, pressure, etc. It is important, at the outset,
More informationMATH 260 Homework assignment 8 March 21, px 2yq dy. (c) Where C is the part of the line y x, parametrized any way you like.
MATH 26 Homework assignment 8 March 2, 23. Evaluate the line integral x 2 y dx px 2yq dy (a) Where is the part of the parabola y x 2 from p, q to p, q, parametrized as x t, y t 2 for t. (b) Where is the
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationMath 32B Discussion Session Week 10 Notes March 14 and March 16, 2017
Math 3B iscussion ession Week 1 Notes March 14 and March 16, 17 We ll use this week to review for the final exam. For the most part this will be driven by your questions, and I ve included a practice final
More informationMath 234. What you should know on day one. August 28, You should be able to use general principles like. x = cos t, y = sin t, 0 t π.
Math 234 What you should know on day one August 28, 2001 1 You should be able to use general principles like Length = ds, Area = da, Volume = dv For example the length of the semi circle x = cos t, y =
More informationENGI Gradient, Divergence, Curl Page 5.01
ENGI 94 5. - Gradient, Divergence, Curl Page 5. 5. The Gradient Operator A brief review is provided here for the gradient operator in both Cartesian and orthogonal non-cartesian coordinate systems. Sections
More informationGRE Math Subject Test #5 Solutions.
GRE Math Subject Test #5 Solutions. 1. E (Calculus) Apply L Hôpital s Rule two times: cos(3x) 1 3 sin(3x) 9 cos(3x) lim x 0 = lim x 2 x 0 = lim 2x x 0 = 9. 2 2 2. C (Geometry) Note that a line segment
More information2 Feynman rules, decay widths and cross sections
2 Feynman rules, decay widths and cross sections 2.1 Feynman rules Normalization In non-relativistic quantum mechanics, wave functions of free particles are normalized so that there is one particle in
More information