Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation
|
|
- Bryan Lindsey
- 6 years ago
- Views:
Transcription
1 BULLETIN of the Malaysan Mathematcal Scences Socety Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave Propagaton L.H. Wryanto Department of Mathematcs, Insttut Teknolog Bandung Jalan Ganesha 10 Bandung, Indonesa leo@dns.math.tb.ac.d Abstract. A model of nterfacal-wave propagaton s solved numercally based on one-dmensonal tme doman Boussnesq equatons, usng a predctor-corrector method. The numercal procedure s able to show the effect of non-lnearty of the model by observng the wave speed of soltary waves. The procedure s then used to smulate the wave propagaton generated, on the left boundary, by a snusodal functon Mathematcs Subject Classfcaton: 76B07 (76B55) Key words and phrases: Boussnesq equatons, nterfacal wave, predctor-corrector method. 1. Introducton The exstence of a class of nonlnear waves n a densty-stratfed medum has been demonstrated by several nvestgators, such as Benney [1], Benjamn [2], Ono [6], and Kubota, Ko & Dobbs [5]. They derved a sngle equaton dfferently as the model of ther nternal waves. For two flud system Segur & Hammack [7] and Cho & Camassa [3] derved the model of the nterfacal waves n a KdV-type equaton. Another model descrbng the nterfacal waves s Boussnesq type. The dervaton of the model can be seen n Grmshaw & Pudjaprasetya [4] who use Hamltonan formulaton. Most of the equatons were then solved analytcally to descrbe soltary waves. Numercal solutons of the equatons are an alternatve way to observe the wave propagaton. Wryanto [8] solved the KdV equaton for the nterfacal wave by fnte dfference method. The procedure was able to smulate the wave propagaton n one drecton and ts deformaton at the nterface that can be explaned as a class of one-drecton waves from the model of Boussnesq type. In ths paper the propagaton of the nterfacal waves s observed numercally from the model of Boussnesq equatons. Ths numercal approach s able to show the propagaton for more general type of waves than just soltary, such as snusodal. The procedure s constructed by predctor-corrector method, and wll be used to Receved: October 11, 2004; Revsed: March 16, 2005.
2 164 L.H. Wryanto show the characterstc such as n analytcal works of soltary waves,.e. the wave breaks up nto two, and each wave propagates n dfferent drecton. Meanwhle, to preserve the soltary wave propagates n one drecton, t s requred a relaton for the ntal elevaton and velocty. In lnear case, the coeffcent q c s determned quanttatvely, and t s larger than the result from the analytcal formulaton as the nonlnear effect of Boussnesq equatons. The wave speed correspondng to the one-drecton wave s another quantty that can be determned followng q c. In presentng ths paper, we frst descrbe the Boussnesq equatons as the model of nterfacal wave n secton 2. To provde the condton of one-drecton wave, we derve the relaton between the elevaton and velocty based on depth average n secton 3. In the next secton, detals of the numercal procedure are gven, and followed by presentng the numercal computatons n secton 5. All descrbed n ths paper s fnally concluded n secton Boussnesq equatons The problem s the moton of gravty wave at the nterface of two-flud system havng dfferent densty ρ 1 for the upper flud and ρ 2 for the lower flud as shown n Fgure 1. Note that we then use varables wth subscrbe 1 for the upper flud and subscrbe 2 for the lower flud. In the undsturbed fluds, the depth of each layer s h 1 and h 2 measured from the nterface that s chosen as the horzontal x-axs. Therefore the nterfacal wave s presented as the elevaton of the nterfacal y(x, t) from the undsturbed level, where t s as the tme varable. Meanwhle, the flud system s bounded above and bellow by flat plane. Fgure 1. Sketch of the coordnates The fluds are assumed to be ncompressble and nvscd, and the flow s rrotatonal. Ths allows us to formulate the problem n potental functons φ 1 and φ 2 satsfyng Laplace equaton (2.1) 2 φ = 0 for = 1, 2 n each flud doman, the knematcs condtons (2.2) (2.3) (2.4) φ 1x = 0 on y = h 1 φ 2x = 0 on y = 2 y t + φ x y x = φ x on y = y(x, t) from upper and lower fluds, and dynamc condton on the nterface (2.5) ρ 1 {φ 1t + 1 } { 2 (φ2 1x + φ 2 1y) + gy = ρ 2 φ 2t + 1 } 2 (φ2 2x + φ 2 2y) + gy
3 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave Propagaton 165 Equaton (2.1) followed (2.2)-(2.5) s approxmated nto the Boussnesq equatons for the nterfacal wave as derved n [4], u t + g(ρ 1 ρ 2 )y x + 2νuu x = 0 (2.6) y t + h1h2 ρ 1h 2+ρ 2h 1 u x + 2ν(yu) x + 2βu xxx = 0 where g s acceleraton due to gravty, (2.7) ν = 1 ρ 2 h 2 1 ρ 1h (ρ 1 h 2 + ρ 2 h 1 ) 2, β = h2 1 h2 2 ρ 1 h 1 + ρ 2 h 2 6 (ρ 1 h 2 + ρ 2 h 1 ) 2 and u s defned as (2.8) u = (ρ 2 φ 2 ρ 1 φ 1 ) x on y = y(x, t). For convenence, Equaton (2.6) s non-dmensonalzed by takng h 1 as the unty of length and gh 1 as the unty of speed. Therefore two parameters nvolved n the equatons are (2.9) h = h 2 h 1, ρ = ρ 2 ρ 1. The equatons are then wrtten n normalzed varables whch wll show the domnatng terms n the equatons, by defnng (2.10) η = y a, ξ = x a, τ = t a, ϑ = u a where a s a reference value of y such as the order of ampltude. These transform (2.6) to the form (2.11) ϑ τ + ( (ρ 1)η + a(ρ ) h2 ) 2(h + ρ) 2 ϑ = 0 ξ ( h (2.12) η τ + ρ + h ϑ + a(ρ h2 ) (h + ρ) 2 ηϑ + h2 a(ρh + 1) 3(ρ + h) 2 ϑ ξξ For small a, (2.11) and (2.12) can be reduced to a smple wave equaton of ϑ or η by elmnatng one of both unknowns, and can be solved analytcally by separatng the ndependent varable. On the other hand, the effect of the nonlnear and thrd dervatve terms wll be observed by solvng (2.11) and (2.12) numercally. 3. One-drecton wave Equatons (2.11) and (2.12) are able to explan wave propagaton n two drectons. In specal relaton of the ntal condtons ϑ and η, the equatons gve a wave travellng only n one drecton. We derve that relaton by consderng the depth-average horzontal velocty defned as (3.1) φ1x = 1 1 aη (3.2) φ2 x = 1 h + aη 1 aη φ 1x dy φ 2x dy n non-dmensonal varables correspondng to (2.10). ) for the upper flud for the lower flud ξ = 0
4 166 L.H. Wryanto The ntegraton n (3.1) and (3.2) can be expressed n η derved as followed. For the lower flud, we start from Laplace equaton (2.1) ntegrated on the thckness of the lower flud (3.3) aφ 2xx + φ 2yy dy = 0. From the knematcs condtons (2.3) and (2.4), and relaton (3.4) (3.3) becomes (3.5) x φ 2x dy = a φ 2x aη η x + x φ 2x dy + η t = 0. φ 2xx dy The next step, we suppose that the wave travels to the rght. Ths s expressed n η as (3.6) η t + cη x = 0 wth wave speed (3.7) c = h(ρ 1) h + ρ obtaned from the leadng order of (2.11) and (2.12). We then substtute (3.6) to (3.5), and ntegrate wth respect to x, gvng (3.8) φ 2x dy = cη The constant of ntegraton s zero as η, φ 2x 0 for x. Relaton (3.8) s then used to elmnate the ntegral n (3.2), so that we have (3.9) φ2x = cη h + aη. When the above procedure s followed for the upper flud, (13) becomes (3.10) φ1x = cη aη 1. Results (3.9) and (3.10) are used to approxmate u n Boussnesq equatons to obtan one-drecton travelng wave,.e. (3.11) u ( ρ h + aη 1 ) cη. aη 1
5 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave Propagaton Numercal procedure In nvestgatng the propagaton of the nterfacal waves, the model n form of Boussnesq equatons s solved numercally. A numercal procedure s developed for (2.11) and (2.12). A predctor-corrector method s chosen as t can be appled to the equatons explctly. To smplfy n explanng the method, (2.11) and (2.12) are expressed n form ϑ τ = F(η, ϑ) where η τ = G(η, ϑ) (4.1) ( F(η, ϑ) = (ρ 1)η + a(ρ ) h2 ) 2(h + ρ) 2 ϑ2 ξ (4.2) ( h G(η, ϑ) = ρ + h ϑ + a(ρ h2 ) (h + ρ 2 ) ηϑ + h2 a(ρh + 1) 3(ρ + h) 2 ϑ ξξ The equatons are fnte-dfferenced on a grd n ξ and τ. The dscrete ndependent varables are expressed as ξ = ξ, τ = n τ wth level n referrng to nformaton at the present, known tme level. The predctor step s the thrd-order explct Adam-Bashforth scheme gven by (4.3) ϑ n+1 (4.4) η n+1 = ϑ n + τ 12 = η n + τ 12 [ 23F n 16F n 1 [ 23G n 16G n 1 + 5F n 2 ] + 5G n 2 ] The predcted values (4.3) and (4.4) are then used to evaluate F n+1 and G n+1 from (4.1) and (4.2) whch wll be used to correctng ϑ n+1 and η n+1 above. The scheme for t s the fourth-order Adams-Moulton method, gven by (4.5) ϑ n+1 (4.6) η n+1 = ϑ n + τ 24 = η n + τ 24 [ 9F n+1 [ 9G n F n 5F n G n 5G n 1 + F n 2 ] ) + G n 2 ]. Ths scheme s also calculated explctly. In runnng the numercal procedure, the predctor step requests two prevous tme levels at n 1 and n 2. Ths must be treated specally at n = 0 and 1 as we need values before the ntal ones. Ths problem can be solved by defnng the same value at the ntal condton. Another problem s boundary condtons and some values of ϑ n and η n outsde the observaton doman. We suppose to observe the wave propagaton n ξ [0, L] and we dscretze as ξ = ξ, ξ = L, = 0, 1, 2,, M. M. ξ
6 168 L.H. Wryanto Table 1. Data of the left-gong wave (a) and the rght-gong wave (b). t x a m (a) t x a m (b) The fnte dfference for the dervatve (up to thrd) of ϑ and η n the observaton doman needs values from = 2 to = M +2. As the values outsde the doman are provded every tme level n, we can overcome ths problem by lnearzng those values of both ϑ and η to the values n the doman. Ths can be done as the assumpton that the doman s relatvely long so that the boundary gves very small effect and can be neglected to the wave whch s propagatng far from the boundary. 5. Numercal results The numercal procedure descrbed n the prevous secton s used to observe the propagaton of nterfacal waves for varous ntal values ϑ(ξ, 0) and η(ξ, 0). Most results presented here are obtaned from calculaton usng ξ = 0.1, τ = 0.02, a = 0.1; and waves are plotted n non-dmensonal varables y versus x. The propagaton and deformaton of waves are shown by plottng the result of calculaton at some tme levels t n the same coordnates, but we shft the values y upper for hgher tme levels. Fgure 2 s typcal nterfacal waves when the ntal condtons are soltary n the form (5.1) y(x, 0) = A sech 2 [b(x xo)] for constants A, b, xo, and the horzontal velocty u = 0. The plot s the numercal result for quanttes h = 0.8, ρ = 1.5, and A = 0.15, b = 1.0, xo = From the calculaton we found that the soltary wave breaks up nto two, and then each wave travels n dfferent drecton. The ampltude and the wave speed are calculated from Table 1a representng tme t, poston of wave peak x and the ampltude a m at that tme for the left-gong wave, and Table 1b for the other wave. We present the data for the last 8 curves as shown n Fgure 2 where the ntal wave has broken up nto two waves. The thrd column of both tables ndcates that the both waves have same ampltude, but t decreases by ncreasng tme. Our calculaton shows that up to t = the maxmum heght reaches Meanwhle, plot of the frst column versus second one produces a lne x = 0.4t for Table 1a, and x = 0.4t for Table 1b. These are the poston of the wave peak at tme t.
7 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave Propagaton 169 Fgure 2. Propagaton of nterfacal wave y(x, 0) = 0.15 sech 2 (x 67.08), u(x,0) = 0 for h = 0.8, ρ = 1.5 (a) (b) Fgure 3. Propagaton of wave for h = 0.8, ρ = 1.5 wth (a) u(x,0) = 0.7 y(x,0), (b) u(x, 0) = 1.5 y(x,0) The coeffcent of the lne represents the wave speed,.e. both waves have the same wave speed 0.4. The negatve sgn ndcates left drecton. Now when the above ntal soltary wave s followed by u(x, 0) = q y(x, 0), wth constant q, the wave wll break up nto two havng dfferent ampltude. Fgure 3a shows the wave deformaton for q = 0.7. The rght-gong wave has wave speed c R = and ampltude Meanwhle the other wave has wave speed c L = and ampltude For larger q the left-gong wave has smaller ampltude, and q c = s the crtcal value where no left-gong wave. Ths crtcal value s dfferent wth the coeffcent of the lnearzed of (3.11),.e , snce we solve the Boussnesq equatons nvolvng the nonlnear terms. Therefore, the dfference s the corrector factor of the lnearzed (3.11). For q > q c the left-gong wave appears wth negatve ampltude. We show ths n Fgure 3b for q = 1.5. In comparson wth the KdV-type equaton, we found that the numercal soluton of Boussnesq equatons can explan the character n KdV equaton by settng u(x, 0) = q c y(x, 0). For soltary waves deformng to some soltons, our procedure s able to show them usng the ntal elevaton as gven n (5.1) for a sutable value b wth respect to the ampltude. Fgure 4a shows a soltary wave deformng to two soltons. We obtan that result as the soluton of Bossnesq equatons usng ntal condton y(x, 0) = 0.4 sech 2 [0.5(x 10.0)] u(x, 0) = y(x, 0)
8 170 L.H. Wryanto (a) (b) (c) Fgure 4. Wave deformaton from soltary to two soltons (a), three soltons (b), and appearng a tran of waves (c). Table 2. The crtcal value q c and wave speed c for varous h and ρ. ρ = 1.2 ρ = 1.5 h q c c (a) h q c c (b) and parameters h = 0.8, ρ = 1.5. For smaller values b the ntal soltary wave can deforms to more soltons, and larger values b a tran of small waves appears behnd the man soltary. We show these n Fgure 4b and 4c, correspondng to b = 0.2 and b = 1.7 respectvely. For other values h and ρ, the crtcal value q c was determned as descrbed above. The correspondng wave propagates to the rght wth wave speed c as gven n Table 2. The ampltude of the wave was observed durng propagatng and we found that t s relatvely constant around a m = For larger values h, our numercal procedure s not able to gve smooth waves, snce ths ncreases the coeffcent of ϑ ξξξ n (12), and ths term domnates the equaton.
9 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave Propagaton 171 Fgure 5. Snusodal waves generated on the left as the boundary and propagatng along the nterface. Fgure 5 shows the propagaton of a typcal snusodal wave comng n the observaton doman havng ρ = 1.5 and h = 0.8. At the begnnng the nterface s undsturbed, and the left boundary s assumed to be a wave generator wth ( u(0, t) = y(0, t) = 0.15 sn(ct) ρ y(0, t) + h 1 y(0, t) 1 ) cy(0, t) where c s the ncomng wave speed defned n (3.7). From our calculaton, we obtaned that the elevaton of the nterface forms waves propagatng to the rght, and deformaton occurs as shown n Fgure 5. The crest of the snusodal deforms to soltary-lke wave and breaks to two waves. Ths s also obtaned from the KdV model, see [8]. 6. Conclusons Numercal solutons of Boussnesq equatons for nterfacal wave have been presented. The model was solved by a predctor-corrector method. We demonstrated the solutons as the wave propagaton for typcal soltary and snusodal. The ntal condton of the horzontal velocty plays an mportant role n observng the wave deformaton, manly break up to two waves travellng n dfferent drecton. In case where the ntal velocty u s lnear to the ntal elevaton, we found that the wave travels only n one drecton for a certan value of the coeffcent q c of that lnear relaton. The class of these one-drecton waves can explan the character of solutons n KdV model. The non-lnear effect of Boussnesq equatons s shown n values q c and wave speed c whch are dfferent from the lnearzed (3.11) and (3.7) for varous h and ρ. Acknowledgment. The research of ths report was supported by QUE (qualty undergraduate educaton) project for Mathematcs Insttut Teknolog Bandung. The author thanks Dr. Pudjaprasetya for stmulatng dscussons. References [1] D. J. Benney, Long nonlnear waves n flud flows, J. Math. Physcs. 45 (1966), [2] T. B. Benjamn, Internal waves of permanent form n fluds of great depth, J. Flud Mech. 29 (1967),
10 172 L.H. Wryanto [3] W. Cho and R. Camassa, Weakly nonlnear nternal waves n a two-flud system, J. Flud Mech. 313 (1996), [4] R. Grmshaw and S.R. Pudjaprasetya, Hamltonan formulaton for the descrpton of nterfacal soltary waves, Nonlnear processes n geophyscs 5 (1998), [5] T. Kubota, D. R. S. Ko, and L. D. Dobbs, Weakly-nonlnear, long nternal gravty waves n stratfed fluds of fnte depth, J. Hydronoutcs (1978), [6] H. Ono, Algebrac soltary wave n stratfed fluds, J. Phys. Soc. Japan 39 (1975), [7] H. Segur and J. L. Hammack, Solton models of long nternal waves, J. Flud Mech. 118 (1982), [8] L. H. Wryanto, A fnte dfference method on KdV equaton of nternal waves, (n preparaton).
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationThe effect of a background shear current on large amplitude internal solitary waves
PHYSICS OF FLUIDS 8, 03660 006 The effect of a background shear current on large ampltude nternal soltary waves Wooyoung Cho a Department of Mathematcal Scences, Center for Appled Mathematcs and Statstcs,
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationOn the internal soliton propagation over slopeshelf
Journal of Physcs: Conference Seres PAPER OPEN ACCESS On the nternal solton propagaton over slopeshelf topography To cte ths artcle: Albert Sulaman 7 J. Phys.: Conf. Ser. 856 3 Vew the artcle onlne for
More informationPROPERTIES OF SURFACE AND INTERNAL SOLITARY WAVES. Kei Yamashita 1 and Taro Kakinuma 2
PROPERTIES OF SURFACE AND INTERNAL SOLITARY WAVES Ke Yamashta and Taro Kaknuma Numercal solutons o surace and nternal soltary waves are obtaned through a new method where advecton equatons on physcal quanttes
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationAn Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors
An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationarxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationJournal of Computational Physics
Journal of Computatonal Physcs 30 (0) 0 030 Contents lsts avalable at ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp An teratve method to solve a regularzed model
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More information( ) + + REFLECTION FROM A METALLIC SURFACE
REFLECTION FROM A METALLIC SURFACE For a metallc medum the delectrc functon and the ndex of refracton are complex valued functons. Ths s also the case for semconductors and nsulators n certan frequency
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationAP Physics 1 & 2 Summer Assignment
AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationModeling acoustic transducer surface waves by Transmission Line Matrix method
Modelng acoustc transducer surface waves by Transmsson Lne Matrx method Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Peter-Chrstan Eccardt Semens AG, CT MS, München Wllam
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics
ECOOMICS 35*-A Md-Term Exam -- Fall Term 000 Page of 3 pages QUEE'S UIVERSITY AT KIGSTO Department of Economcs ECOOMICS 35* - Secton A Introductory Econometrcs Fall Term 000 MID-TERM EAM ASWERS MG Abbott
More informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More information), it produces a response (output function g (x)
Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More informationNONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM
Advanced Steel Constructon Vol. 5, No., pp. 59-7 (9) 59 NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM M. Abdel-Jaber, A.A. Al-Qasa,* and M.S. Abdel-Jaber Department of Cvl Engneerng, Faculty
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More information