Modeling the scalar wave equation with Nyström methods
|
|
- Edgar Francis
- 5 years ago
- Views:
Transcription
1 GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 200 ; P. T151 T158, 7 FIGS / Modeling the calar wave equation with Nytröm method Jing-Bo Chen 1 ABSTRACT High-accuracy numerical cheme for modeling of the calar wave equation baed on Nytröm method are developed in thi paper. Space i dicretized by uing the peudopectral algorithm. For the time dicretization, Nytröm method are ued. A fourth-order ymplectic Nytröm method with peudopectral patial dicretization i preented. Thi cheme i compared with a commonly ued econdorder cheme and a fourth-order nonymplectic Nytröm method. For a typical time-tep ize, the econd-order cheme exhibit patial diperion error for long-time imulation, while both fourth-order cheme do not uffer from thee error. Numerical comparion how that the fourth-order ymplectic algorithm i more accurate than the fourth-order nonymplectic one. The capability of the ymplectic Nytröm method in approximately preerving the dicrete energy for long-time imulation i alo demontrated. INTRODUCTION Seimic modeling i an important foundation of exploration eimology. High-accuracy eimic modeling cheme are increaing in demand a computing capacity increae. In thi context, high-accuracy cheme are developed uing Nytröm time-tepping method and peudopectral patial dicretization. Numerical cheme for modeling of the calar wave equation involve dicretization of both pace and time. Thi dicretizing proce can be accomplihed in two way. The firt way i to dicretize both pace and time imultaneouly. The finite-difference method Claerbout, 1985 belong to thi category. The econd way i firt to dicretize pace to obtain a ytem of ordinary differential equation ODE with time a a variable and then to contruct numerical cheme by dicretizing time for the ytem of ODE. The peudopectral method Gazdag, 1981 and the finite-element method Ciarlet and Lion, 1991 fall into thi econd category. To form a unified viewpoint and facilitate our undertanding of numerical imulation cheme, we will accept that all numerical modeling cheme can be obtained by the econd approach to dicretization. Therefore, to contruct modeling cheme for the calar wave equation, we can firt dicretize pace by finite-difference method, peudopectral method, or finite-element method. Second, we numerically integrate the reulting ytem of ordinary differential equation by dicretizing time. Note that the obtained ytem of ODE i of econd order. When numerically olving thi ytem directly, we need to approximately evaluate additional tarting value in addition to initial value. The higher the order of the cheme, the more complex the evaluation of the additional tarting value. Nytröm method are numerical method deigned for econdorder differential equation. Nytröm 1925 firt conidered thi a a implification of Runge-Kutta method. Hairer et al ytematically developed Nytröm method. There are two benefit in uing Nytröm method. Firt, we do not need to evaluate the additional tarting value. Second, by introducing an intermediate variable, the calar wave equation can be cat into a Hamiltonian ytem. Uing Nytröm method, we can develop the correponding ymplectic method. The ymplectic method have remarkable capability for long-time computation. Thi i becaue the ymplectic propertie guarantee that the numerical olution evolve in the ame ytem a the olution of the original continuou differential equation Feng, 1993; Sanz-Serna and Calvo, In thi paper, I preent the peudopectral method and explore high-order time dicretization. I follow thi with a dicuion of Nytröm method. I then develop a fourth-order ymplectic Nytröm method with peudopectral patial dicretization. Finally, I illutrate the performance of thi cheme with numerical experiment that compare it with a econd-order cheme and a fourth-order nonymplectic Nytröm method. PSEUDOSPECTRAL METHODS A mentioned above, we can ue finite-difference method, peudopectral method, or finite-element method for the pace dicretization. I ue peudopectral method, for which the error reult Manucript received by the Editor July 18, 2005; revied manucript received February, 200; publihed online September 5, ChineeAcademy of Science, Intitute of Geology and Geophyic, P.O. Box 9825, Beijing , China. chenjb@mail.igca.ac.cn. 200 Society of Exploration Geophyicit. All right reerved. T151
2 T152 Chen mainly from the time dicretization. For higher accuracy, we need higher-order time dicretization. The iue of computing additional tarting value become more evident in thi ituation. Peudopectral method for modeling of the calar wave equation were preented by Gazdag The main point of the peudopectral method i to expre the wavefield under conideration in term of a complete et of orthogonal bai function whoe derivative are known exactly. In practice, Fourier peudopectral method are uually employed. Conider the calar wave equation 2 u t 2 = c2 2 u x u y u z 2, where u x,y,z,t i the wavefield and c x,y,z i the velocity. Let u = u 1,1,1,..., u Nx,N y,n z T, where T repreent tranpoe and u i,j,l are the wavefield value at dicrete location, i.e., u i,j,l u i x, j y, l z,t ; i = 1,..., N x ; j = 1,..., N y ; l = 1,..., N z. The value x, y, and z are grid increment in the x-, y-, and z-direction; and N x, N y, and N z are the number of grid line in the x-, y-, and z-direction, repectively. The emi-dicrete ytem reulting from the peudopectral method for equation 1 i 1 d 2 u dt 2 = c2 F 1 w F u, 2 where F and F 1 repreent 3D forward and invere finite Fourier tranform, repectively, and w = w 1,1,1,..., w Nx,N y,n z T with w i,j,l = k 2 xi + k 2 yj + k 2 zl, where k xi, k yj, and k zl are dicrete wavenumber in the x-, y-, and z-direction, repectively. The aterik denote array multiplication between vector. For example, uppoe that p = p 1,p 2,..., p m and q = q 1,q 2,..., q m, then p q = p 1 q 1,p 2 q 2,..., p m q m. To obtain the final modeling cheme, we further dicretize equation 2 in time. A tandard econd-order time difference i often ued and give u t 0 +2 t =2u t 0 + t u t 0 + c 2 F 1 w F u t 0 + t t 2, 3 where t i the time-tep ize and t 0 i the initial time. A uual, we take t 0 = 0. To compute u 2 t uing equation 3, we need to know u 0 and u t, which can be computed from the initial condition of equation 1: u 0 and u 0 / t. We can evaluate u 0 directly from u 0. To obtain u t, a third-order cheme baed on a Taylor erie expanion wa given by Gazdag 1981 a 3 u t = i u 0 t i i=0 t i. 4 i! In equation 4, u 0 and u 0 / t are known, 2 u 0 / t 2 i defined by equation 1, and 3 u 0 / t 3 i obtained by 3 u 2 3 = c2 t x y z 2 u t, which i obtained by ubtituting u/ t for u in equation 1. Finally, we can compute u t directly from u t. The main error in cheme equation 3 i aociated with the time differencing. Thu, to improve accuracy, we need to ue a higherorder time difference. A natural choice i the fourth-order time difference u 4 t =1u 3 t 30u 2 t +1u t u 0 c 2 F 1 w F u 2 t t 2. However, we can how that cheme 5 i unconditionally untable by uing the tandard pectral analyi. Uing the Taylor erie and the wave equation, another approach for high-order time differencing wa preented by Etgen 198. But thi approach till require the computation of additional tarting value. For example, for a fourthorder method, we need a fifth-order cheme to compute the tarting value intead of the third-order cheme 4. We are faced with a cumberome procedure for computing tarting value. Therefore, an alternative approach i needed. NYSTRÖM METHODS The baic idea of Nytröm 1925 method i a follow: Firt introduce an intermediate variable to reduce the econd-order wave equation to an equivalent firt-order ytem; then apply Runge-Kutta method to the firt-order ytem and implify by taking advantage of the pecial form of the firt-order ytem. The obtained Nytröm method do not need to evaluate tarting value and have coniderably le computational cot than Runge-Kutta method applied directly to the firt-order ytem for econd-order equation uch a equation 2. For detail, refer to Hairer et al For our purpoe, conider a econd-order ytem of ordinary differential equation written in the form d 2 y dt 2 = f y. ANytröm method for ytem read Z i = f y 0 + c i tz 0 + t 2 a ij Z j, i = 1, 2,...,, j=1 y 1 = y 0 + tz 0 + t 2 b iz i, i=1 z 1 = z 0 + t b i Z i, i=1 where z = dy/dt, y 0 = y 0, z 0 = z 0, y 1 y t, z 1 z t and where c i, a ij, b i, and b i are contant that determine the order of the method. The numerical olution y 1 i obtained through auxiliary variable Z i in cheme 7. The tarting value y 0 and z 0 = dy 0 /dt in cheme 7 are both known, and no additional tarting value are needed. Scheme 7 require computation of z 1, but thi i eay to accomplih becaue the main evaluation of Z i have been done previouly in the computation of y 1. The order condition via the algebraic equation atified by c i, a ij, b i, and b i have been obtained by uing tree theory in Hairer et al By olving thee algebraic equation, we can obtain the correponding method. Conider a table fourth-order explicit method Qin and Zhu, 1991 : 5 7
3 Modeling with Nytröm method T153 c 1 = 3+ 3 b 1 = b 1 = 3 2 3, c 2 = 3 3, b 2 = 3+ 3 a 21 = 2 3 3, a 32 =,, c 3 = 3+ 3,, b 3 = 1+ 3,, b 2 = 1 2, b 3 = 3+2 3, a 11 = a = a 13 = a 22 = a 23 = a 31 = a 33 =0. With thee coefficient, cheme 7 become Z 1 = f y tz 0, method. Therefore, Nytröm method can be applied directly. We introduce a variable v = du/dt and apply Nytröm method 8 to equation 2 to obtain V 1 = c 2 F 1 w F u tv 0, V 2 = c 2 F 1 w F u tv t 2 V 1, V 3 = c 2 F 1 w F u tv 0 + t2 V 2, u 1 = u 0 + tv 0 + t V 3, V V 2 Z 2 = f y tz t 2 Z 1, v 1 = v 0 + t V V V 3. 9 Z 3 = f y tz 0 + t2 Z 2, y 1 = y 0 + tz 0 + t Z 3, z 1 = z 0 + t Z Z 2 Z Z Z 3. Another benefit of Nytröm method i that we can develop ymplectic algorithm if equation poee a ymplectic tructure. In the pat, we only focued on the accuracy of the numerical method when olving differential equation numerically. However, with the tremendou progre in computer technology and numerical analyi theory in recent decade, computation ha become a third category of method in cientific reearch in addition to theory and experiment. Not only i the accuracy of numerical method expected, but tructure-preerving propertie alo are required. Solution to differential equation uually preerve variou tructure uch a ymplectic tructure, multiymplectic tructure, and variou geometric tructure. When olving differential equation numerically, ome numerical method alo preerve thee tructure they are uually called tructure-preerving algorithm, while other violate them. The tructure-preerving method have remarkable capability for long-time computation. The calar wave equation 1 ha a claical Hamiltonian tructure Chen, Therefore, we can develop the correponding tructure-preerving method eeappendix A. A PSEUDOSPECTRAL METHOD WITH FOURTH-ORDER NYSTRÖM TIME DIFFERENCE Now we return to the iue of higher-order time difference for equation 2. Equation 2 i a econd-order ytem of ordinary differential equation obtained from equation 1 by uing peudopectral 8 Scheme 9 i a table fourth-order explicit Nytröm method with accuracy of O t 4 Qin and Zhu, 1991; Hairer et al., Becaue we ue peudopectral patial dicretization, the patial accuracy i of exponential order O exp x Fornberg, 199. Therefore, the total accuracy of cheme 9 i O t 4 + exp x.a pointed out by Hairer et al. 1993, the Nytröm method cheme 9 i coniderably more efficient than Runge-Kutta method. In eimic modeling baed on cheme 9, we compute u 1 and v 1 from known u 0 and v 0 = du 0 /dt. Then we repeat thi proce with u 0 and v 0 replaced by u 1 and v 1, repectively. Higher-order Nytröm time difference are eaily available, and they have the ame format a cheme 9 but with more V i. NUMERICAL EXPERIMENTS In thi ection, we perform numerical experiment to tet the feaibility and performance of the preented cheme. For clarity and implicity, we only conider the 2D cae. A in Gazdag 1981, we ue the initial condition and u x,z,t =0 = exp x 2 + z z 0 2 u 0 =0. t The grid increment are x = z = 50 m. Here, z 0 i a contant that indicate the poition of the ource. In the following numerical example, we ue periodic boundary condition. In the firt example, we ue cheme 9 to imulate wave propagation in a heterogeneou medium. The velocity model Figure 1 conit of three layer with velocitie c = 1000 m/, 2000 m/, and 3000 m/; a traight interface; and a curved interface. The ource i et at x =0m,z = 200 m. The time-tep ize t i Figure 2 how the wavefield at t = 0, 0.438, 0.93, 1.8, 1.5, and The imulated reflection and tranmiion phenomena are qualitatively correct.
4 T154 Chen We now make numerical comparion between cheme 9 and the commonly ued cheme 3. In thi example, we ue the velocity c = 3000 m/ and the time tep ize t = The tarting value for equation 3 i computed through equation 4, which require evaluating the econd-order patial derivative of the initial condition. Therefore, when uing cheme like equation 3 in wave imulation, we uually require that the initial condition have a certain degree of moothne, which ometime i not atified in practice for example, dicrete initial condition. In contrat, when uing cheme 9, we Figure 1. The velocity model ued for the numerical experiment. It conit of three layer with velocitie c = 1000 m/, 2000 m/, and 3000 m/, repectively; a traight interface; and a curved interface. have no pecial requirement on the initial condition, which make cheme like 9 more practical. In Figure 3, we how the impule repone for cheme 3 and 9 after 200, 500, and 5000 time tep, repectively. For time tep equal to or le than 200, impule repone for both cheme 3 and 9 do not exhibit patial diperion error. After 500 time tep, the patial diperion become evident in the wavefront computed with cheme 3, which contrat with the clean wavefront computed with cheme 9 Figure 3c and d. After 5000 time tep t =30, the wavefront computed with cheme 3 ha blurred very badly, but the wavefront computed with cheme 9 i till harp Figure 3e and f. Thee obervation indicate that for a typical time-tep ize of 0.00, the commonly ued cheme 3 i not uitable for lengthy imulation. Now we perform imulation by cheme 3 with maller time-tep ize of 0.004, 0.003, 0.002, and To compute the wavefront at 30, the number of time tep i 7500, 10,000, 15,000, and 30,000, repectively. Figure 4 how the imulation reult. We ee that diperion error decreae a time-tep ize decreae. The reult with time-tep ize of i imilar to that obtained with cheme 9 with a time-tep ize of 0.00 Figure 3d. Thi indicate that to obtain the impule repone at 30, cheme 9 with a time-tep ize of 0.00 i more efficient than cheme 3 with a time-tep ize of We can explain thi a follow: Although for each time tep cheme 9 need to compute three pair of FFT forward and invere, which i three time a many a that of the cheme 3, the timetep ize ued in cheme 9 i ix time a large a that ued in cheme 3. Becaue the main computational cot in both cheme 9 and 3 i the computation of FFT, the overall computational cot of cheme 9 i only half that of cheme 3. To demontrate the advantage of the ymplectic method, we now make numerical comparion between the fourth-order ymplectic Nytröm cheme 9 with the following fourth-order nonymplectic Nytröm cheme: V 1 = c 2 F 1 w F u 0, V 2 = c 2 F 1 w F u tv t2 V 1, V 3 = c 2 F 1 w F u 0 + tv t2 V 2, u 1 = u 0 + tv 0 + t 2 1 V V 2, v 1 = v 0 + t 1 V V V Figure 2. Wavefield evolution with time computed with cheme 9. Ditance and time are meaured in meter and econd, repectively. Figure 5 how the impule repone computed with cheme 10 after 5000 time tep. The time-tep ize i We ee that the reult i imilar to that computed with the ymplectic cheme 9 Figure 3f. To examine more cloely the reult obtained by thee two fourth-order cheme, amplitude curve at a fixed point x = 100 m, z = 100 m over the time near time tep are hown in Figure a. We ee dicrepancy in the amplitude curve computed by the two fourth-order cheme. In Figure b, another amplitude curve i added, which i computed by the nonymplectic cheme 10 with a maller time tep ize of It can be een that the amplitude curve computed by the ymplectic cheme 9 with a time-tep ize of 0.00 agree well with the amplitude curve com-
5 Modeling with Nytröm method T155 puted by the nonymplectic cheme 10 with a maller time-tep ize of Thi indicate that the fourth-order ymplectic algorithm cheme 9 i more accurate than the fourth-order nonymplectic algorithm cheme 10 for the ame time-tep ize. To how how the amplitude curve computed by the nonymplectic cheme 10 with maller time-tep ize approach the amplitude curve computed by the ymplectic cheme 9 with a time-tep ize of 0.00, we plot the following amplitude-difference curve Figure c : AD t t = S 0.00 t NS t t, where AD t t denote the amplitude-difference curve; S 0.00 t, the amplitude curve computed by the ymplectic cheme 9 with a timetep ize of 0.00 ; and NS t t, the amplitude curve computed by the nonymplectic cheme 10 with a time-tep ize of t. Here we take t = 0.00, 0.004, 0.002, and We can ee that the magnitude of the curve of the amplitude difference diminihe with decreaing t. From the above numerical comparion, we can conclude that to achieve approximately the ame accuracy, the ymplectic algorithm cheme 9 i much more efficient than the nonymplectic algorithm cheme 10, conidering the time-tep ize ued. Thi i a great advantage. Another important advantage of a ymplectic algorithm i it ability to approximately preerve energy for lengthy imulation Feng, 1993; Sanz-Serna and Calvo, For the wave equation u tt = c 2 u xx + u zz with a periodic boundary condition, the true olution preerve the energy v M 2 + c 2 u 2 x + u 2 z dxdz, where v = u t and M i the integral domain. Now in our numerical experiment, we monitor the dicrete energy 2 G n = v n i,l 2 + c u n i+1,l i,l 2 + c u n i,l+1 z n u i,l n u i,l x 2 x z, 2 Figure 4. Impule repone at 30. The reult are computed with cheme 3 with time-tep ize of 0.004, 0.003, 0.002, and 0.001, repectively. Figure 3. Impule repone computed with cheme 3 and 9. a and b 200 time tep; c and d 500 time tep; e and f 5000 time tep. View a, c, and e are computed with cheme 3; b, d, and f are computed with cheme 9. Figure 5. Impule repone after 5000 time tep computed with cheme 10 with a tep-time ize of 0.00.
6 T15 Chen Figure 7. Dicrete energy curve G n veru n for the ymplectic cheme 9 and the nonymplectic cheme 10. a N = 5000; b N = 100,000. half of the energy i lot in the energy curve for the nonymplectic cheme 10. Thi demontrate the capability of the ymplectic algorithm for approximately preerving the dicrete energy for lengthy computation. CONCLUSIONS Figure. Amplitude and amplitude-difference curve at a fixed point x = 100 m, z = 100 m over the time near 30. a Amplitude curve computed by the ymplectic cheme 9 with t = 0.00 blue line and the nonymplectic cheme 10 with t = 0.00 red line. b Amplitude curve computed by the ymplectic cheme 9 with t = 0.00 blue line ; the nonymplectic cheme 10 with t = 0.00 red line ; and the nonymplectic cheme 10 with t = green line. c Amplitude-difference curve: t = 0.00 red line ; t = blue line ; t = yellow line ; and t = green line. where v n i,l v i x,l z,n t, u n i,l u i x,l z,n t, and n = 0,1,2,..., N. Here, N i the number of time tep. Figure 7a how the dicrete energy curve G n a a function n for both cheme 9 and 10, with N = For the ymplectic cheme 9, the energy curve fluctuate about a contant energy. For the nonymplectic cheme 10, a numerical lo in energy i oberved. The correponding reult for N = 100,000 i hown in Figure 7b. For thi very long computation, the energy curve for the ymplectic cheme 9 till fluctuate about a contant energy while nearly Nytröm time differencing with peudopectral patial dicretizaiton i preented for modeling of the calar wave equation. Two fourth-order Nytröm cheme with peudopectral patial dicretizaiton are developed; one i ymplectic and the other i nonymplectic. Numerical experiment demontrate three reult. Firt, thee fourth-order cheme can be ued for long-time imulation with a time-tep ize of 0.00, wherea the commonly ued econd-order cheme exhibit evere diperion error for thi computation time and time-tep ize. Second, the fourth-order cheme with a timetep ize of 0.00 are more efficient than the econd-order cheme with a mall time-tep ize of And third, the ymplectic algorithm i more accurate than the nonymplectic algorithm and ha a better capability of approximately preerving the dicrete energy for lengthy imulation. Preerving propagation energy i only one feature of preerving ymplectic tructure. The ymplectic algorithm alo have other feature uch a better accuracy, which can be explained by the theory of backward error analyi Hairer et al., According to thi theory, the olution of both ymplectic and nonymplectic algorithm for a Hamiltonian ytem formally atify a perturbed ytem.
7 Modeling with Nytröm method T157 However, the perturbed ytem atified by the ymplectic algorithm retain the Hamiltonian tructure; the ytem atified by the nonymplectic algorithm doe not. ACKNOWLEDGMENTS I would like to thank profeor Arcone and anonymou reviewer for valuable uggetion and correction, that led to a great improvement of thi paper. Thi work i upported by the National Natural Science Foundation of China under grant and alo partially by the Chinee Academy of Science with Key Project of Knowledge Innovation KZCX1-SW-18. du = v, dt dv dt = c2 F 1 w F u. Further, we can reformulate equation A-3 a d dt v u = 0 I 0 H I u, H v A-3 A-4 APPENDIX A HAMILTONIAN STRUCTURE In thi appendix, I preent the Hamiltonian tructure and the correponding ymplectic algorithm for the calar wave equation. Introducing an intermediate variable v = du/dt, equation 1 i equivalent to u t = v, c2 v t = 2 Equation A-1 can be reformulated a where H = 1 2 v 2 + c u x u y u z 2. t v u = u 2 x u 0 H H, v A-1 A-2 + y u 2 + z u 2 dxdy i the Hamiltonian ytem, and denote the variational derivative. Equation A-2 i an infinite-dimenional Hamiltonian ytem. The evolution of the wavefield with time i characterized by a ymplectic tructure du Ù dvdxdy, which can be viewed a an antiymmetric quadratic form. For the exact definition of ymplectic tructure, a well a the infinite-dimenional Hamiltonian ytem and variational derivative, ee Olver In olving equation A-2 numerically, we can ue finite-difference, finite-element, or peudopectral method to dicretize the pace and obtain a finite-dimenional ytem. Equation 2, which i obtained with a peudopectral method, can be cat into a finite-dimenional Hamiltonian ytem. Firt, introducing v = du/dt, we can rewrite equation A-2 a where I i an identity matrix and H = 1 2 vt v + c 2 u T Du. The matrix D i a econd-order pectral differential matrix which atifie F 1 w F u = Du. For detail about pectral differential matrix, ee Chen and Qin Equation A-4 i a tandard finite-dimenional Hamiltonian ytem. A finite-dimenional Hamiltonian ytem i jut a ytem of ordinary differential equation that ha the form of equation A-4. The true olution of equation A-4 preerve the ymplectic tructure du Ù dv. In other word, if we uppoe that the olution of equation A-4 i P = F P 0, where P = u,v T and P 0 = u t =0,v t =0 T, then the fact that F P 0 preerve the ymplectic tructure du Ù dv i equivalent to F P 0 atifying F P 0 J F T P 0 = J, where J = 0 I I 0. Here, F/ P 0 i the Jacobian of the vector-valued function F P 0. A function atifying the above equality i alo called a ymplectic mapping. Thu, it i concluded that the time evolution of the eimic wavefield i a ymplectic mapping. In eimic modeling, we need to olve equation A-4 numerically. Suppoe that F P 0 i the numerical olution of equation A-4 obtained by ome numerical algorithm. Of coure, the numerical olution F P 0 i the approximation of the true olution F P 0. Doe the numerical olution alo have other propertie? In fact, ome numerical olution are ymplectic mapping while other are not. Anumerical method i called a ymplectic algorithm if the reulting numerical olution i a ymplectic mapping. In the pat, we mainly focued on the accuracy of the numerical olution to a Hamiltonian ytem; it tructure-preerving propertie were not conidered eriouly. With the introduction of ymplectic algorithm and it great ucce in many phyical field, tructurepreerving algorithm have developed into a very active and promiing reearch area. The Hamiltonian framework developed here provide a bai for applying ymplectic algorithm in eimic modeling. Baed on thi framework, we hould develop numerical olution of equation A-4 that are alo ymplectic mapping, i.e., ymplectic algorithm. Variou ymplectic algorithm for Hamiltonian ytem have been developed. We can directly apply the ymplectic algorithm to equation A-4. Baed on a finite-difference patial dicretization, Hamiltonian ytem alo can be obtained Luo
8 T158 Chen et al., However, thi kind of patial dicretization uually caue patial diperion error. Now conider a Nytröm formulation for equation A-4: V i = c 2 F 1 w F u 0 + c i tz 0 + t 2 a ij V j, u 1 = u 0 + tv 0 + t 2 b iv i, i=1 If the coefficient in cheme A-5 atify b i = b i 1 c i, i = 1,...,, j=1 b i b j a ij = b j b i a ji, i, j = 1,...,, v 1 = v 0 + t b i V i. i=1 A-5 then cheme A-5 i a ymplectic algorithm Sanz-Serna and Calvo, It i eay to check that the coefficient ued in cheme 9 atify the above condition; therefore, cheme 9 i a ymplectic algorithm. REFERENCES Chen, J. B., 2004, Multiymplectic geometry for the eimic wave equation: Communication in Theoretical Phyic, 41, Chen, J. B., and M. Z. Qin, 2001, Multiymplectic Fourier peudopectral method for the nonlinear Schrödinger equation: Electronic Tranaction on NumericalAnalyi,, Ciarlet, P. G., and J. L. Lion, 1991, Handbook of numerical analyi: North- Holland. Claerbout, J. F., 1985, Imaging the earth interior: Blackwell Scientific Publication, Inc. Etgen, J. T., 198, High-order finite-difference revere migration with the two-way non-reflecting wave equation: SEP-48: Stanford Univerity. Feng, K., 1993, The collected work of Feng Kang: Defence Indutrial Pre. Fornberg, B., 199, A practical guide to peudopectral method: Cambridge Univ. Pre. Gazdag, J., 1981, Modeling of the acoutic wave equation with tranform method: Geophyic, 4, Hairer, E., C. Lubich, and G. Warnner, 2002, Geometric numerical integration: Springer-Verlag, Berlin. Hairer, E., S. P. Nøett, and G. Warnner, 1993, Solving ordinary differential equation I: Springer-Verlag, Berlin. Luo, M. Q., H. Liu, and Y. M. Li, 2001, Hamiltonian decription and ymplectic method of eimic wave propagation: Chinee Journal of Geophyic, 44, 0 8. Nytröm, E. J., 1925, Über die numeriche integration von differentialgleichungen: Acta Societe Scientiarum Fennicae, 50, Olver, P. J., 1993, Application of Lie group to differential equation: Springer-Verlag, New York. Qin, M. Z., and W. J. Zhu, 1991, Canonical Runge-Kutta-Nytröm method for econd order ODE : Computational Mathematic Application, 22, Sanz-Serna, J. M., and M. Calvo, 1994, Numerical Hamiltonian problem: Chapman and Hall.
High-order time discretizations in seismic modeling
GEOPHYSICS, VOL. 72, NO. 5 SEPTEMBER-OCTOBER 27; P. SM115 SM2, 5 FIGS. 1.119/1.275424 High-order time dicretization in eimic modeling Jing-Bo Chen 1 ABSTRACT Seimic modeling play an important role in exploration
More information696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of
Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationReformulation of Block Implicit Linear Multistep Method into Runge Kutta Type Method for Initial Value Problem
International Journal of Science and Technology Volume 4 No. 4, April, 05 Reformulation of Block Implicit Linear Multitep Method into Runge Kutta Type Method for Initial Value Problem Muhammad R., Y. A
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationMulti-dimensional Fuzzy Euler Approximation
Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com
More informationDiagonally Implicit Symplectic Runge-Kutta Methods with Special Properties
Appl. Math. Inf. Sci. 9, No. 1L, 11-17 (2015) 11 Applied Mathematic & Information Science An International Journal http://dx.doi.org/10.12785/ami/091l02 Diagonally Implicit Symplectic Runge-Kutta Method
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationEE Control Systems LECTURE 6
Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More information722 Chen Xiang-wei et al. Vol. 9 r i and _r i are repectively the poition vector and the velocity vector of the i-th particle and R i = dm i dt u i; (
Volume 9, Number 10 October, 2000 1009-1963/2000/09(10)/0721-05 CHINESE PHYSICS cfl 2000 Chin. Phy. Soc. PERTURBATION TO THE SYMMETRIES AND ADIABATIC INVARIANTS OF HOLONOMIC VARIABLE MASS SYSTEMS * Chen
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationHybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification
www.ccenet.org/ma Modern Applied Science Vol. 6, No. ; February Hybrid Projective Dilocated Synchronization of Liu Chaotic Sytem Baed on Parameter Identification Yanfei Chen College of Science, Guilin
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
More informationEstimating floor acceleration in nonlinear multi-story moment-resisting frames
Etimating floor acceleration in nonlinear multi-tory moment-reiting frame R. Karami Mohammadi Aitant Profeor, Civil Engineering Department, K.N.Tooi Univerity M. Mohammadi M.Sc. Student, Civil Engineering
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationAccuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / 11-23 BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal,
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationUSPAS Course on Recirculated and Energy Recovered Linear Accelerators
USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular
More informationTRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL
GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National
More informationEXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal
EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS Otto J. Roech, Hubert Roth, Aif Iqbal Intitute of Automatic Control Engineering Univerity Siegen, Germany {otto.roech,
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationLinear Motion, Speed & Velocity
Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding
More informationarxiv: v2 [nucl-th] 3 May 2018
DAMTP-207-44 An Alpha Particle Model for Carbon-2 J. I. Rawlinon arxiv:72.05658v2 [nucl-th] 3 May 208 Department of Applied Mathematic and Theoretical Phyic, Univerity of Cambridge, Wilberforce Road, Cambridge
More informationComputers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order
Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More informationAlternate Dispersion Measures in Replicated Factorial Experiments
Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More informationDetermination of the local contrast of interference fringe patterns using continuous wavelet transform
Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,
More informationUnbounded solutions of second order discrete BVPs on infinite intervals
Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b
More informationChapter 7. Root Locus Analysis
Chapter 7 Root Locu Analyi jw + KGH ( ) GH ( ) - K 0 z O 4 p 2 p 3 p Root Locu Analyi The root of the cloed-loop characteritic equation define the ytem characteritic repone. Their location in the complex
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationFiltered Leapfrog Time Integration with Enhanced Stability Properties
Journal of Applied Mathematic and Phyic, 6, 4, 4-7 Publihed Online July 6 in SciRe. http://www.cirp.org/ournal/amp http://dx.doi.org/.46/amp.6.474 Filtered Leapfrog Time Integration with Enhanced Stability
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationLocal Fractional Laplace s Transform Based Local Fractional Calculus
From the SelectedWork of Xiao-Jun Yang 2 Local Fractional Laplace Tranform Baed Local Fractional Calculu Yang Xiaojun Available at: http://workbeprecom/yang_iaojun/8/ Local Fractional Laplace Tranform
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More informationSOLVING THE KONDO PROBLEM FOR COMPLEX MESOSCOPIC SYSTEMS
SOLVING THE KONDO POBLEM FO COMPLEX MESOSCOPIC SYSTEMS V. DINU and M. ÞOLEA National Intitute of Material Phyic, Bucharet-Magurele P.O. Box MG-7, omania eceived February 21, 2005 Firt we preent the calculation
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationGeneral Field Equation for Electromagnetism and Gravitation
International Journal of Modern Phyic and Application 07; 4(5: 44-48 http://www.aacit.org/journal/ijmpa ISSN: 375-3870 General Field Equation for Electromagnetim and Gravitation Sadegh Mouavi Department
More informationGNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase
GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,
More informationStochastic Neoclassical Growth Model
Stochatic Neoclaical Growth Model Michael Bar May 22, 28 Content Introduction 2 2 Stochatic NGM 2 3 Productivity Proce 4 3. Mean........................................ 5 3.2 Variance......................................
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationEffects of vector attenuation on AVO of offshore reflections
GEOPHYSICS, VOL. 64, NO. 3 MAY-JUNE 1999); P. 815 819, 9 FIGS., 1 TABLE. Effect of vector attenuation on AVO of offhore reflection J. M. Carcione ABSTRACT Wave tranmitted at the ocean bottom have the characteritic
More informationRELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS
www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin
More informationThe Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor
The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura
More informationGreen-Kubo formulas with symmetrized correlation functions for quantum systems in steady states: the shear viscosity of a fluid in a steady shear flow
Green-Kubo formula with ymmetrized correlation function for quantum ytem in teady tate: the hear vicoity of a fluid in a teady hear flow Hirohi Matuoa Department of Phyic, Illinoi State Univerity, Normal,
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More information84 ZHANG Jing-Shang Vol. 39 of which would emit 5 He rather than 3 He. 5 He i untable and eparated into n + pontaneouly, which can alo be treated a if
Commun. Theor. Phy. (Beijing, China) 39 (003) pp. 83{88 c International Academic Publiher Vol. 39, No. 1, January 15, 003 Theoretical Analyi of Neutron Double-Dierential Cro Section of n+ 11 B at 14. MeV
More informationRepresentation Formulas of Curves in a Two- and Three-Dimensional Lightlike Cone
Reult. Math. 59 (011), 437 451 c 011 Springer Bael AG 14-6383/11/030437-15 publihed online April, 011 DOI 10.1007/0005-011-0108-y Reult in Mathematic Repreentation Formula of Curve in a Two- and Three-Dimenional
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationTHE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER
Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationDYNAMIC MODELS FOR CONTROLLER DESIGN
DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that
More informationA BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH
Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan
More informationτ ), respectively. Then
Downloaded /0/ to 9.7.47.34. Reditribution ubject to SEG licene or copyright; ee Term of Ue at http://library.eg.org/ Quality factor affect channel wave propagation in 3D iotropic vicoelatic medium Hui
More informationPredicting the shock arrival time using 1D-HD solar wind model
Article Geophyic April 1 Vol.55 No.11: 153158 doi: 1.17/11434-9-61-8 SPECIAL TOPICS: Predicting the hock arrival time uing 1D-HD olar wind model ZHANG Ying 1,*, CHEN JingYi & FENG XueShang 1 Intitute of
More informationSTOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND
OPERATIONS RESEARCH AND DECISIONS No. 4 203 DOI: 0.5277/ord30402 Marcin ANHOLCER STOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND The generalized tranportation problem
More informationNoether symmetry and non-noether conserved quantity of the relativistic holonomic nonconservative systems in general Lie transformations
Vol 16 No 11, November 2007 c 2007 Chin. Phy. Soc. 1009-1963/2007/1611/3182-05 Chinee Phyic and IOP Publihing Ltd Noether ymmetry and non-noether conerved quantity of the relativitic holonomic nonconervative
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationarxiv: v2 [math.nt] 30 Apr 2015
A THEOREM FOR DISTINCT ZEROS OF L-FUNCTIONS École Normale Supérieure arxiv:54.6556v [math.nt] 3 Apr 5 943 Cachan November 9, 7 Abtract In thi paper, we etablih a imple criterion for two L-function L and
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationLTV System Modelling
Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction
More informationResearch Article Reliability of Foundation Pile Based on Settlement and a Parameter Sensitivity Analysis
Mathematical Problem in Engineering Volume 2016, Article ID 1659549, 7 page http://dxdoiorg/101155/2016/1659549 Reearch Article Reliability of Foundation Pile Baed on Settlement and a Parameter Senitivity
More informationSimulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System
Proceeding of the 7th IAPRI World Conference on Packaging Simulation Study on the Shock Propertie of the Double-Degree-of-Freedom Cuhioning Packaging Sytem Xia Zhu, Qiaoqiao Yan, Xiaoling Yao, Junbin Chen,
More informationwhere F (x) (called the Similarity Factor (SF)) denotes the function
italian journal of pure and applied mathematic n. 33 014 15 34) 15 GENERALIZED EXPONENTIAL OPERATORS AND DIFFERENCE EQUATIONS Mohammad Aif 1 Anju Gupta Department of Mathematic Kalindi College Univerity
More informationTarzan s Dilemma for Elliptic and Cycloidal Motion
Tarzan Dilemma or Elliptic and Cycloidal Motion Yuji Kajiyama National Intitute o Technology, Yuge College, Shimo-Yuge 000, Yuge, Kamijima, Ehime, 794-593, Japan kajiyama@gen.yuge.ac.jp btract-in thi paper,
More informationFast Convolutional Sparse Coding (FCSC)
Fat Convolutional Spare Coding (FCSC) Bailey ong Department of Computer Science Univerity of California, Irvine bhkong@ic.uci.edu Charle C. Fowlke Department of Computer Science Univerity of California,
More informationInvariance of a Partial Differential Equation of Fractional Order under the Lie Group of Scaling Transformations
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 8197 1998 ARTICLE NO AY986078 Invariance of a Partial Differential Equation of Fractional Order under the Lie Group of Scaling Tranformation Evelyn
More informationSome Sets of GCF ϵ Expansions Whose Parameter ϵ Fetch the Marginal Value
Journal of Mathematical Reearch with Application May, 205, Vol 35, No 3, pp 256 262 DOI:03770/jin:2095-26520503002 Http://jmredluteducn Some Set of GCF ϵ Expanion Whoe Parameter ϵ Fetch the Marginal Value
More informationThe type 3 nonuniform FFT and its applications
Journal of Computational Phyic 206 (2005) 1 5 Short Note The type 3 nonuniform FFT and it application June-Yub Lee a, *, Lelie Greengard b a Department of Mathematic, Ewha Woman Univerity, 11-1 Daehyundong,
More informationFrequency dependent attenuation and dispersion in patchysaturated
Frequency dependent attenuation and diperion in patchyaturated porou rock Huixing Zhang, ritopher A. Innanen. Key Lab of Submarine Geocience and Propecting Technique, MOE, Ocean Univerity of China;. Department
More informationTheoretical study of the dual harmonic system and its application on the CSNS/RCS
Theoretical tudy of the dual harmonic ytem and it application on the CSNS/RCS Yao-Shuo Yuan, Na Wang, Shou-Yan Xu, Yue Yuan, and Sheng Wang Dongguan branch, Intitute of High Energy Phyic, CAS, Guangdong
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More information