The discrete Moser Veselov algorithm for the free rigid body
|
|
- Ami Pitts
- 6 years ago
- Views:
Transcription
1 The discrete Moser Veselov algorithm for the free rigid body 1/100 Antonella Zanna and Robert McLachlan Report in Informatics nr. 255 (2003). To appear in Found. Comp. Math. Antonella Zanna University of Bergen, Norway anto
2 Report Documentation Page Form Approved OMB No Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 03 JAN REPORT TYPE N/A 3. DATES COVERED - 4. TITLE AND SUBTITLE The discrete MoserVeselov algorithm for the free rigid body 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) University of Bergen, Norway 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release, distribution unlimited 11. SPONSOR/MONITOR S REPORT NUMBER(S) 13. SUPPLEMENTARY NOTES See also ADM001749, Lie Group Methods And Control Theory Workshop Held on 28 June July 2004., The original document contains color images. 14. ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT UU a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified 18. NUMBER OF PAGES 31 19a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18
3 Overview The subject of this talk is the numerical solution of the free RB equations, M = [M, Ω], M = ΩJ + JΩ, 2/100 where M, Ω are skew-symmetric matrices and J is a diagonal matrix with positive entries. M is the matrix of body momenta Ω is the matrix of body angular velocity Often the above equations are associated with the equations that give the configuration of the body in the fixed frame, Q = QΩ, Q SO(N). The Discrete Moser Veselov description of the rigid body Integrability of the discrete algorithm ward error analysis of the the DMV algorithm Higher order integrable approximations Numerical experiments and comparisons with other methods Explicit methods for the 3 3 case
4 The Moser Veselov discrete version of the dynamics of a Rigid Body The Lagrangian of the continuous RB equations, is the kinetic energy, L = 1 2 tr(ω M) = 1 2 tr( Ω2 J ΩJΩ) = tr(ω JΩ), (1) where we take into account that Ω = Ω and that the trace is invariant under cyclic permutations. Following (Marsden, Pekarsky & Shkoller 1999), discretise Ω = X 1 Ẋ, where X SO(N) is the configuration of the body, using a finite difference approximation of the derivative, 3/100 Ω = X 1 Ẋ 1 h X k+1(x k+1 X k ), X k, X k+1 SO(N), which gives L 1 h 2 tr(j X k X k+1 J JX k+1x k X k X k+1 JX k+1x k ). Due to the orthogonality of the X k s and the cyclicity of the trace, the first and the last term cancel, and moreover, we can write L 2 h 2 tr(x kjx k+1). Up a scaling factor, this is precisely the discrete Lagrangian of Moser and Veselov.
5 Consider the fuctional S(X) determined by S = k tr(x k JX k+1) where X = {X k } with X k O(N) and J is a diagonal matrix. To obtain the stationary points of S, we consider tr(x k JXk+1) 1 tr(λ k (X k Xk I)), 2 (where Λ k = Λ k k is a Lagrange multiplier), and δs = 0 becomes k 4/100 X k+1 J + X k 1 J = Λ k X k, from which, multiplying by Xk on the right and taking into consideration the symmetry of Λ k, X k+1 JXk + X k 1 JXk = Λ k = Λ k = X k JXk+1 + X k JXk 1, (2) hence, the discrete analogue of the angular momentum in space, is conserved. m k = X k JX k 1 X k 1 JX k,
6 In the body variables, setting ω k = X k X k 1 O(N) and M k = X 1 k 1 m kx k 1 = ω k J Jω k so(n) (angular momentum w.r.t. the body), (2) becomes M k+1 = ω k M k ωk M k = ωk J Jω k. the discrete Euler Arnold equation. In the continuous limit: when t k = t 0 + kε, k = 0, 1, 2,..., X k = X(t k ) ω k = Xk X k 1 I εω(t k ), M k ε(jω + ΩJ) = εm(t k ), (3) 5/100 letting ε 0, one obtains the familiar Euler Arnold equations for the motions of the N-dimensional rigid body, M = [M, Ω] M = JΩ + ΩJ, Ω so(n).
7 To solve the discrete Euler Arnold equations (3): For k = 0, 1, 2,..., find ω k SO(N) such that M k = ω k J Jω k. Update M k+1 = ω k M k ω k. By construction, this algorithm preserves exactly momentum and energy (integrable map) is a second order approximation to the continuous rigid body preserves the standard Poisson structure of T so(n), 6/100 {f, g} = tr(m[f M, g M ]), f, g C (so(n)), where f M = ( f/ M i,j ). Note that Also the IMR is second order, preserves all the integrals of the continuous rigid body. Another much used method is a Lie Poisson integrator of McLachlan and Reich (LP2). For the 3 3 RB, it consists in splitting the Hamiltonian H = H 1 + H 2 + H 3 = m2 1 J 2 + J 3 + m2 2 J 1 + J 3 + m2 3 J 1 + J 2 and integrating explicitly (a la Strang) the vector fields of each split Hamiltonian. The method is second order, explicit, preserves the Poisson structure but does not preserve H.
8 Hamiltonian error Hamiltonian error 4 x IMR time 2 x McL splitting time 3 x MVlu Hamiltonian error time MVlu McL splitting IMR MVqr /100 Error in the Hamiltonian function H in the interval [0, 100] and for h = 1 2. The components of the vector m k for h = m = m 1 m 2 m 3 ˆm = M = 0 m 3 m 2 m 3 0 m 1 m 2 m 1 0
9 Integrability of the Moser Veselov equation Recall the Moser Veselov equation, 8/100 M k = ωk J Jω k, Mk = M k, ωk ω k = I. (4) in tandem with the update M k+1 = ω k M k ωk. The Moser Veselov equation (4) has not a unique solution; Lemma 1 (Moser and Veselov) Equation (4) is equivalent to the factorization (I λm k λ 2 J 2 ) = (ωk + λj)(ω k λj) Rewrite the above factorization as (I λm k λ 2 J 2 ) = A k (λ)a k (λ) A k (λ) = (ω k λj) and define the mapping with image point M k+1 such that By construction, (I λm k+1 λ 2 J 2 ) = A k (λ)a k (λ). (I λm k+1 λ 2 J 2 ) = A k (λ)(i λm k λ 2 J 2 )A 1 k (λ), hence the map is an integrable as it is an isospectral flow.
10 BEA for DMV Recall the DMV equations and the continuous RB equations 9/100 M k+1 = ω k M k ω k, M = [M, Ω] M k = ω k J Jω k, M = ΩJ + JΩ, where ω k I hω(t k ). We wish to write M k+1 = Φ h (M k ) = M k + h[m k, Ω k ] + h 2 d 2 + h 3 d 3 + h 4 d 4 +, and find the modified vector field M = [ M, Ω] + hf 2 ( M, Ω) + h 2 f 3 ( M, Ω) + h 3 f 4 ( M, Ω) + (5) such that Φ h (M k ) equals the solution M(t k+1 ) at time t k+1 = t 0 + (k + 1)h of the modified vector field (5). To find Φ k (h), we write ω k = exp( hω 0 h 2 Ω 1 h 3 Ω 2 h 4 Ω 3 h 5 Ω 4 + ), (6) where Ω 0, Ω 1, Ω 2,..., are skew-symmetric matrices computed so that ω k J Jω k = h(ω(t k )J + JΩ(t k )). (7)
11 we obtain h(ω(t k )J + JΩ(t k )) = h(ω 0 J + JΩ 0 ) + h 2 (Ω 1 J + JΩ (Ω2 0J JΩ 2 0)) + h 3 (Ω 2 J + JΩ [(Ω 0Ω 1 + Ω 1 Ω 0 ), J] (Ω3 0J + JΩ 3 0)) +. Comparing left and right-hand-sides, it is trivially observed that the order-h term disappears if Ω 0 = Ω (to simplify notation, we omit the dependence of Ω on t k ). In order to annihilate the h 2 -term, we require that Ω 1 J + JΩ (Ω2 0J JΩ 2 0) = 0. 10/100 Recall that M = ΩJ + JΩ and hence M = Ω J + JΩ. On the other hand, M = [M, Ω] = (Ω 2 J JΩ 2 ). Hence we can write O = Ω 1 J + JΩ M = Ω 1 J + JΩ (Ω J + JΩ ) and the identity is satisfied by if and only if Ω 1 = 1 2 Ω. (8) In general, the algorithm to derive Ω i, for i = 1, 2,..., is 1. Find the coefficient of h i+1 in (7) and set it equal to zero. This will give an equation of the type Ω i J +JΩ i = C i J +JC i +[D i, J]. Note that the terms C i J +JC i have an odd occurrence of the Ω j s, while the terms of the type [D i, J] have an even occurrence of the Ω j s. 2. Use the derivatives of M and Ω to express the term [D i, J] as C i J + J C i.
12 3. Deduce Ω i = C i + C i. Ω 0 = Ω Ω 1 = 1 2 Ω Ω 2 = 1 4 Ω 1 6 Ω3 Ω 3 = 1 8 Ω 1 24 (5Ω2 Ω + 2ΩΩ Ω + 5Ω Ω 2 ) 11/100 The functions Ω i Once the Ω i s are known, substituting back in M k+1 = ω k M kω k and using the well known identity exp(x)y exp( X) = exp adx Y = k=0 1 k! adk X(Y ), where ad X (Y ) = [X, Y ] and, recursively, ad k X(Y ) = [X, ad k 1 X (Y )], we find the expressions for the functions d i in terms of the Ω i 1, Ω i 2,..., Ω 0, d i = i ( 1) j j! j=1 k 1 +k 2 + +k j =i j ad Ωk1 ad Ωk2 ad Ωkj M, k 1,... k j {0, 1,..., i 1}. (9) d 2 = 1([M, 2 Ω ] + [[M, Ω], Ω]), d 3 = 1[M, 4 Ω ] + 1[[M, 4 Ω ], Ω] + 1[[M, Ω], 4 Ω ] + 1[[[M, Ω], Ω], Ω] 1[M, 6 6 Ω3 ], (10) d 4 =...,
13 Taylor expansion of the solution of the modified equation Consider d dtỹ = f(ỹ) + hf 2(ỹ) + h 2 f 3 (ỹ) +, where f(m) = [M, Ω] = [M, J 1 M] is the original vector field of the RB equations, where J is a linear operator, defined such that J Ω = ΩJ + JΩ = M. Putting ỹ(t) = M(t), we expand the solution of the above equation in a Taylor series and collect corresponding powers of h, ỹ(t + h) = M(t) + hf(m) + h (f 2 2 (M) + 1 ) 2! f f(m) + h (f 3 3 (M) + 1 2! (f f 2 (M) + f 2f(M)) + 1 ) 3! (f (f, f)(m) + f f f(m)) +, 12/100 where f is considered as a linear operator, f as a bilinear operator and so on and so forth. In our case, f (z)(m) = [z, J 1 M] + [M, J 1 z] = [z, Ω] + [M, J 1 z] f (z 1, z 2 )(M) = 2[z 1, J 1 z 2 ], and, since f is quadratic, f and all the other higher derivatives equal zero. At this point it is important to stress an important difference between the expressions for the modified vector field of (Hairer, Lubich & Wanner 2002) and ours. While the vector field discussed in (Hairer et al. 2002) is in R n, hence the f is a symmetric quadratic operator, this is not the case for our vector field which is on matrices, thus f (f f, f) f (f, f f).
14 This non-commutative case is discussed with more generality in (Munthe-Kaas & Krogstad 2002). However, we observe that all the terms containing combinations of f, f and f correspond simply to higher derivatives of f. The mixed terms are treated instead specifically. After some algebra, we have f 2 = d 2 1 2! f f(m) = O, 13/100 f 3 = d 3 1 3! (f (f, f)(m) + f f f(m)) = 1 12 [M, Ω [Ω, Ω ] 2Ω 3 ], f 4 = d 4 1 4! M (iv) 1 2! (f f 3 + f 3f) = O, (11) f 5 = d 5 1 M (v) 1 (f f 5! 2! 4 + f 4 f + 1 d (f 2! dt 3f + f f 3 )) = 1 [M, 80 Ω(iv) ] 1 [M, [Ω, 80 Ω ]] + 3 [M, 40 Ω5 Ω ΩΩ ] + 1 [M, 80 [Ω, Ω ]] 1 [M, 40 ΩΩ Ω] 1 [M, 20 Ω2 Ω + Ω Ω 2 ] + 1 [M, 20 [Ω3, Ω ]] 1 [M, 40 Ω 2 Ω + ΩΩ 2 + Ω[Ω, Ω ]Ω].
15 1 m / m m time The DMV solution of the RB equations (dotted line), the exact solution (solid line) and the trajectories corresponding to the modified vector fields f + h 2 f 3 (dashed line) and f + h 2 f 3 + h 4 f 5 (dash-dotted line) in the interval [0, 50] with h = 8 10.
16 Some important results about DMV Theorem 2 The DMV is time-reversible, hence f 2i = 0, i = 1, 2,.... Theorem 3 (Moser Veselov) In the 3 3 case, the DMV is a time-reparamtetrisation of the flow of the original vector field of the rigid body. Since the mapping preserves the underlying Poisson structure and all the integrals F i = c i of the system, it commutes with all commuting Hamiltonian flows generated by the F i s, M = {M, F i }. The nonsingular compact level sets T c = i (F i = c i ) consists of a finite union of 1-dimensional tori and on each torus the DMV mapping is a shift along the trajectory depending on the integral quantity H 2. Hence, the DMV solves the modified equation 15/100 M = (1 + h 2 τ 3 + h 4 τ h 2i τ 2i+1 + )[M, Ω], where h is the stepsize of integration and the τ 2i+1, for i = 1, 2,..., are constants.
17 Introduce the constants C J,i,j = J i 1J j 2 + J i 1J j 3 + J i 2J j 3, C J,i = C J,i,i C J = C J,1 = (J 1 + J 2 )(J 1 + J 3 )(J 2 + J 3 ) H 2 = (J 1 + J 2 )(J 1 + J 3 )(J 2 + J 3 )H C J m 2. 16/100 Theorem 4 One has and τ 5 = τ 3 = ((3 det(j)tr(j) + C J,2) m (3C J + tr(j 2 ))H 2 ), 1 ( (3tr(J 4 ) + 27C 40 4 J,2 + 15tr(J 2 )C J + 45 det(j)tr(j))h2 2 + (10C J, det(j)tr(j)c J + 10 det(j)tr(j)tr(j 2 ) + 2C J,2 tr(j 2 ) 28 det(j 2 )) m 2 2H 2 ) + (60 det(j 2 )C J + 3C J, det(j 2 )tr(j 2 ) + 15 det(j)(c J,2,3 + C J,3,2 )) m 4 2. Proof. By direct computation of f 3 and f 5.
18 Higher-order integrable methods For the original RB equations, scaling the initial condition is equivalent to scaling time. In our case, we know that DMV is a time-rescaling of the original RB equation. Therefore we wish to rescale the initial condition to obtain a better approximation of the unscaled original RB. 17/100 I.C. DMV h(ω(t k )J + JΩ(t k )) New I.C. DMV h(ω(t k )J+JΩ(t k )) 1+ τ 3 h 2 + τ 5 h 4 + We perform again the backward error analysis. We set now ω = exp( h Ω 0 h 2 Ω1 + ) and solve for the Ω i s as the skew-symmetric matrices that solve h(1 τ 3 h 2 + ( τ 2 3 τ 5 )h 4 + )(ΩJ + JΩ) = ω J J ω. (12) f 3 = d 3 1 3! M = τ 3 [M, Ω] + d 3 1 3! M = τ 3 [M, Ω] + f 3 = ( τ 3 + τ 3 )[M, Ω], hence, in order to have an order-four scheme, we must set f 3 = 0 which corresponds to the choice τ 3 = τ 3.
19 After further computations, one has f 5 = 0 τ 5 = τ 5 2τ 2 3. This value of τ 5 gives indeed a method of order six. The new proposed algorithms of order four and six are described below. The DMV4 algorithm: 18/ Compute τ 3 and set M 0 := M(t 0 )h/(1 + h 2 τ 3 ). 2. For k = 0, 1,..., n 1, find the unique w k as above such that M k = ω k J Jω k set M k+1 = ω k M k ω k end 3. Reconstruct M n := M n (1 + h 2 τ 3 )/h M(t n ). The DMV6 algorithm: 1. Compute τ 3, τ 5 and set τ 5 = τ 5 2τ 2 3 and M 0 := M(t 0 )h/(1 + h 2 τ 3 + h 4 τ 5 ). 2. For k = 0, 1,..., n 1, find the unique w k as above such that M k = ω k J Jω k set M k+1 = ω k M k ω k end 3. Reconstruct M n := M n (1 + h 2 τ 3 + h 4 τ 5 )/h M(t n ).
20 Some numerical experiments We consider with initial condition and matrix J given as J = m 0 = and compare the DMV explicit scheme with the Hamiltonian-splitting method LP2 of (McLachlan 1993) and the Implicit Midpoint Rule (IMR), H = m2 1 J 2 + J 3 + m2 2 J 1 + J 3 + m2 3 J 1 + J 2 = H 1 + H 2 + H 3 19/100 m k+1 = m k + hf( m k + m k+1 ), 2 where f(m) = m ( J) 1 m, J = J 2 + J J 1 + J J 1 + J 2.
21 / Error at T= LP2 IMR DMV DMV4 DMV Step size h Error versus step size computed at T = 100 for the methods LP2, IMR, DMV, DMV4, DMV6.
22 LP2 IMR DMV DMV4 DMV /100 Error at T= Floating point operations Floating point operations versus accuracy (T = 100) for the methods LP2, IMR, DMV, DMV4, DMV6. The roots of P (λ) are recomputed at each step, use QR with pivoting, (DMV 22 LP2 per step).
23 10 2 LP2 IMR 10 0 DMV DMV4 DMV6 RATTLE / Floating point operations versus accuracy (T = 100) for the methods LP2, IMR, DMV, DMV4, DMV6 and RATTLE6. The roots of P (λ) are computed once, use LU instead of QR (DMV 19
24 LP2 per step). Method h = 1 16 h = 1 2 h = 1.2 h = 2.2 h = 2.5 h = 4 LP e e e e e e+00 IMR e e e e e e-01 DMV e e-01 NaN NaN NaN NaN DMV e e e e e e-01 DMV e e e-02 NaN e e-01 23/100 Error for the various methods and selected step sizes
25 Connections with matrix Ricatti equations Consider the matrix equation M = XJ JX. (13) Cardoso & Leite (2003) shown that every solution of (13) (not necessarily orthogonal) is of the form X = (M/2 + S)J 1, for some symmetric matrix S. Furthermore, X is a orthogonal solution of (13) if and only if S is a symmetric solution of the Riccati equation S 2 + S(M/2) + (M/2) S (M 2 /4 + J 2 ) = 0. (14) Riccati equations are associated to symplectic matrices. In our case, the symplectic matrix is [ ] M I H sympl = 2 M 2 + J 2 M. (15) 4 2 If M J 2 is positive definite, it has been shown in (Cardoso & Leite 2003) that (14) has a unique solution S which is symmetric, positive definite, and such that the eigenvalues of W = M/2 + S have positive real parts. This matrix W is precisely the same matrix in Moser & Veselov (1991), from which one obtains ω = W J 1. 24/100
26 Algorithm(Cardoso & Leite 2003): Compute X, the unique solution of (13) in the special orthogonal group SO(n). 1. Find a real Schur form of H sympl, Q H sympl Q = [ T11 T 12 O T 22 ], (16) 25/100 where T 11 and T 22 are block upper-triangular matrices such that the real parts of the spectrum of T 11 are positive and the real parts of the spectrum of T 22 are negative definite. 2. Partition Q accordingly, Then, compute 3. Compute [ Q11 Q Q = 12 Q 21 Q 22 X = S = Q 21 Q ]. ( ) M 2 + S J 1. Some computational details Compute real Schur forms by QR iterations for eigenvalues (Golub & van Loan 1989) Cost: O((2N) 3 ) operations (implicit methods for ODEs: O(N 3 )) N being the dimension of M.
27 The case N = 3 26/100 In this case, it is possible to find an explicit spectral decomposition of H sympl (without the QR eigenvalue method) construct the real Schur decomposition (16) and hence X from the eigenstructure of H sympl. This yields an explicit numerical method for the reduced RB equations.
28 The eigenvalues of the matrix H sympl, H sympl = [ M 2 I M J 2 M 2 are the solutions of the quadratic eigenvalue problem P (λ) = det(λ 2 I λm J 2 ) = 0. Without loss of generality, we assume that J is diagonal, with entries J 1, J 2, J 3. Then, P (λ) = λ 6 λ 4 (J J J 2 3 m 2 12 m 2 13 m 2 23) + λ 2 (J 2 1 J J 2 1 J J 2 2 J 2 3 m 2 12J 2 3 m 2 13J 2 2 m 2 23J 2 1 ) J 2 1 J 2 2 J 2 3 = λ 6 λ 4 (tr(j 2 ) m 2 ) + λ 2 (C J,2 H 2 ) det(j 2 ). ] (17) (18) 27/100 C J,i,j = J i 1J j 2 + J i 1J j 3 + J i 2J j 3, C J,i = C J,i,i C J = C J,1 H 2 = (J 1 + J 2 )(J 1 + J 3 )(J 2 + J 3 )H C J m 2. Reduce to a cubic equation (compute the roots explicitely)
29 Schematical procedure Compute eigenvalues/eigenvectors of H sympl : H sympl [ Y1 Y 2 ] = [ Y1 Y 2 ] [ Λ + Λ ], Re Λ + 0, (the eigenvectors need not be orthogonal and may be complex). Y 1, Y 2 R 6 3, Λ ± R /100 Orthogonalize the eigenvectors (by Grahm-Schmidt or QR), [Y 1, Y 2 ] = QR, so that H sympl Q = QRΛR 1 is the complex Schur form. Reduce to a real Schur form by considering real/imaginary part (complex Givens rotation). Compute S = Q 21 Q 1 11, X = (M/2 + S)J 1. We don t need all the eigenvectors, just Y 1. Don t need R. Avoid complex arithmetic alltogether.
30 The numerical DMV algorithm m 0 hm 0 = h ˆm 0. Compute the eigenvalues of H sympl = H(M k ) solving for P (λ) = 0 as in (18). For t k = t 0 + kh, k = 0, 1, /100 Compute the (real) eigenvectors corresponding to Λ + Compute 3 quadratic eigenvectors (3 matrix factorizations, LU/QR, with pivoting). No need to compute explicitely L or Q. Compute the dependent eigenvectors. Orthogonalize the eigenspace By (modified) Grahm Schmidt or QR. Only the Q factor is needed. Compute S = Q 21 Q 1 11, ω k = (M/2 + S)J 1 Update M k+1 = ω k M kω k Rescale m N M N /h. This algorithm produces an explicit method that is about times more expensive than LP2, the explicit method of McLachlan and Reich.
31 Concluding remarks Explicit algorithms to solve for the N = 3 free rigid body The methods are up to 6th order, completely integrable, possible to increase to arbitrary order 30/100 The cost of the method is about times more expensive than the explicit LP2. The cheaper versions seem to be less stable expecially for large step-size and long time computations Reconstruction equations: Find the configuration X k+1 = X k ω k. The complexive order is still 2 but the error is generally smaller. ward error analysis reveals that some error terms cancel, however the components are reparametrized with 3 different time scales. To obtain higher order reconstructions, one can use interpolation of the vector field XΩ and its higher derivatives using the computed values of Ω k. Optimal step-size?
32 References Cardoso, J. R. & Leite, F. S. (2003), The Moser Veselov equation, Lin. Alg. Applic. Golub, G. H. & van Loan, C. F. (1989), Matrix Computations, 2nd edn, John Hopkins, Baltimore. Hairer, E., Lubich, C. & Wanner, G. (2002), Geometric Numerical Integration, number 31 in Springer Series in Computational Mathematics, Springer, Berlin. Marsden, J. E., Pekarsky, S. & Shkoller, S. (1999), Discrete Euler Poincaré and Lie Poisson equations, Nonlinearity 12, McLachlan, R. I. (1993), Explicit Lie Poisson integration and the Euler equations, Physical Review Letters 71, Moser, J. & Veselov, A. P. (1991), Discrete Version of Some Classical Integrable Systems and Factorization of Matrix Polynomials, Commun. Math. Phys. 139, Munthe-Kaas, H. Z. & Krogstad, S. (2002), On enumeration problems in Lie Butcher theory, Report in informatics, Institutt for informatikk, University of Bergen, Norway. To appear in Future Generation Computer Systems, Special issue. 31/100
Diagonal Representation of Certain Matrices
Diagonal Representation of Certain Matrices Mark Tygert Research Report YALEU/DCS/RR-33 December 2, 2004 Abstract An explicit expression is provided for the characteristic polynomial of a matrix M of the
More informationReport Documentation Page
Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationUse of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications
Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Joseph R. Gabriel Naval Undersea Warfare Center Newport, Rl 02841 Steven M. Kay University of Rhode
More informationP. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE
A wavelet-based generalization of the multifractal formalism from scalar to vector valued d- dimensional random fields : from theoretical concepts to experimental applications P. Kestener and A. Arneodo
More informationAttribution Concepts for Sub-meter Resolution Ground Physics Models
Attribution Concepts for Sub-meter Resolution Ground Physics Models 76 th MORS Symposium US Coast Guard Academy Approved for public release distribution. 2 Report Documentation Page Form Approved OMB No.
More informationHigh-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems
Aerodynamic Issues of Unmanned Air Vehicles Fluid-Structure Interaction High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Raymond E. Gordnier Computational Sciences
More informationReport Documentation Page
Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationREGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)
AL/EQ-TP-1993-0307 REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) John T. Mourand, John C. Crittenden, David W. Hand, David L. Perram, Sawang Notthakun Department of Chemical Engineering
More informationBroadband matched-field source localization in the East China Sea*
Broadband matched-field source localization in the East China Sea* Renhe Zhang Zhenglin Li Jin Yan Zhaohui Peng Fenghua Li National Laboratory of Acoustics, Institute of Acoustics, Chinese Academy of Sciences,
More informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationSystem Reliability Simulation and Optimization by Component Reliability Allocation
System Reliability Simulation and Optimization by Component Reliability Allocation Zissimos P. Mourelatos Professor and Head Mechanical Engineering Department Oakland University Rochester MI 48309 Report
More informationClosed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects Chris Harrison Centre for Marine
More informationThermo-Kinetic Model of Burning for Polymeric Materials
Thermo-Kinetic Model of Burning for Polymeric Materials Stanislav I. Stoliarov a, Sean Crowley b, Richard Lyon b a University of Maryland, Fire Protection Engineering, College Park, MD 20742 b FAA W. J.
More informationBabylonian resistor networks
IOP PUBLISHING Eur. J. Phys. 33 (2012) 531 537 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/33/3/531 Babylonian resistor networks Carl E Mungan 1 and Trevor C Lipscombe 2 1 Physics Department, US
More informationCRS Report for Congress
CRS Report for Congress Received through the CRS Web Order Code RS21396 Updated May 26, 2006 Summary Iraq: Map Sources Hannah Fischer Information Research Specialist Knowledge Services Group This report
More informationWavelets and Affine Distributions A Time-Frequency Perspective
Wavelets and Affine Distributions A Time-Frequency Perspective Franz Hlawatsch Institute of Communications and Radio-Frequency Engineering Vienna University of Technology INSTITUT FÜR NACHRICHTENTECHNIK
More informationReport Documentation Page
Inhibition of blood cholinesterase activity is a poor predictor of acetylcholinesterase inhibition in brain regions of guinea pigs exposed to repeated doses of low levels of soman. Sally M. Anderson Report
More informationAnalysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.
Concept V- Hull Floor Configuration in Two Dimensions Dr. Bijan Khatib-Shahidi & Rob E. Smith 10 November 2010 : Dist A. Approved for public release Report Documentation Page Form Approved OMB No. 0704-0188
More informationOcean Acoustics Turbulence Study
Ocean Acoustics Turbulence Study PI John Oeschger Coastal Systems Station 6703 West Highway 98 Panama City, FL 32407 phone: (850) 230-7054 fax: (850) 234-4886 email: OeschgerJW@ncsc.navy.mil CO-PI Louis
More informationSums of the Thue Morse sequence over arithmetic progressions
Sums of the Thue Morse sequence over arithmetic progressions T.W. Cusick 1, P. Stănică 2 1 State University of New York, Department of Mathematics Buffalo, NY 14260; Email: cusick@buffalo.edu 2 Naval Postgraduate
More informationUSMC Enlisted Endstrength Model
USMC Enlisted Endstrength Model MORS Workshop Personnel and National Security: A Quantitative Approach Working Group #2: Models for Managing Retention Molly F. McIntosh January 27, 2009 Report Documentation
More information7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationSuper-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography
More informationSuper-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography
More informationuniform distribution theory
Uniform Distribution Theory 9 (2014), no.1, 21 25 uniform distribution theory ON THE FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION Florian Luca Pantelimon Stănică ABSTRACT. Here, we show
More informationWeak Turbulence in Ocean Waves
Weak Turbulence in Ocean Waves Yuri V. Lvov Rensselaer Polytechnic Institute Troy, New York 12180-3590 Phone: (518) 276-6893 fax: (518) 276-4824 email: lvov@rpi.edu Award Number N000140210528 http://www.rpi.edu/~lvovy
More informationOn Applying Point-Interval Logic to Criminal Forensics
On Applying Point-Interval Logic to Criminal Forensics (Student Paper) Mashhood Ishaque Abbas K. Zaidi Alexander H. Levis Command and Control Research and Technology Symposium 1 Report Documentation Page
More informationA report (dated September 20, 2011) on. scientific research carried out under Grant: FA
A report (dated September 20, 2011) on scientific research carried out under Grant: FA2386-10-1-4150 First-principles determination of thermal properties in nano-structured hexagonal solids with doping
More informationTHE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS
Annales Univ. Sci. Budapest., Sect. Comp. 40 (2013) nn nnn THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS Florian Luca (Morelia, Mexico) Pantelimon Stănică (Monterey, USA) Dedicated to
More informationSW06 Shallow Water Acoustics Experiment Data Analysis
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. SW06 Shallow Water Acoustics Experiment Data Analysis James F. Lynch MS #12, Woods Hole Oceanographic Institution Woods
More informationMolecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'
Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG' S.A. Burke 1, J.W. Wright 2, M.K. Robinson, B. Bronk 3, and R.L. Warren 1 1 Battelle Dugway Operations 2 U.S. Army Soldier
More informationDIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION
CETN-I-6 3185 DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION PURPOSE : To present a parameterized model of a directional spectrum of the sea surface using an energy spectrum and a value for the
More informationQuantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)
Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) D.N. Kurk, T.E. Beegle, S.C. Spence and R.J. Swatski Report Documentation
More informationA GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS
A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS Michael Thiel*, Truong Vu Bang Giang and Achim Dreher Institute of
More informationScattering of Internal Gravity Waves at Finite Topography
Scattering of Internal Gravity Waves at Finite Topography Peter Muller University of Hawaii Department of Oceanography 1000 Pope Road, MSB 429 Honolulu, HI 96822 phone: (808)956-8081 fax: (808)956-9164
More informationCrowd Behavior Modeling in COMBAT XXI
Crowd Behavior Modeling in COMBAT XXI Imre Balogh MOVES Research Associate Professor ilbalogh@nps.edu July 2010 831-656-7582 http://movesinstitute.org Report Documentation Page Form Approved OMB No. 0704-0188
More informationPredictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance
Predictive Model for Archaeological Resources Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the
More informationAn Invariance Property of the Generalized Likelihood Ratio Test
352 IEEE SIGNAL PROCESSING LETTERS, VOL. 10, NO. 12, DECEMBER 2003 An Invariance Property of the Generalized Likelihood Ratio Test Steven M. Kay, Fellow, IEEE, and Joseph R. Gabriel, Member, IEEE Abstract
More informationRange-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Range-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects Robert I. Odom Applied Physics Laboratory
More informationGeneration and Propagation of Internal Solitary Waves on the Continental Shelf and Slope
Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope Roger H.J. Grimshaw Department of Mathematical Sciences Loughborough University Loughborough, LE11 3TU, UK phone:
More informationSensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts
Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts George Halliwell, MPO/RSMAS, University of Miami Alexander Barth, University
More informationEstimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight
Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight Dale P. Winebrenner Polar Science Center/Applied Physics Laboratory University of Washington
More informationDynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University
Dynamics of Droplet-Droplet and Droplet-Film Collision C. K. Law Princeton University The physical phenomena of droplet-droplet and droplet-film collision in the head-on orientation were studied experimentally
More informationDirect Numerical Simulation of Aeolian Tones
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 27-30, 2003 Direct Numerical Simulation of Aeolian Tones Osamu Inoue 1 1. Institute of Fluid Science, Tohoku University,2-1-1 Katahira,
More informationExtension of the BLT Equation to Incorporate Electromagnetic Field Propagation
Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation Fredrick M. Tesche Chalmers M. Butler Holcombe Department of Electrical and Computer Engineering 336 Fluor Daniel EIB Clemson
More informationImprovements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO
Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO William L. Dimpfl Space Science Applications Laboratory The Aerospace Corporation
More informationStudy of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method
Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Manohar D. Deshpande, NASA Langley Research Center, Hmapton, VA., USA Abstract
More informationFRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES
1 FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES Robert L. Street Environmental Fluid Mechanics Laboratory Department of Civil Engineering Stanford University Stanford, CA 94305-4020 650-723-4969;
More informationReal-Time Environmental Information Network and Analysis System (REINAS)
Calhoun: The NPS Institutional Archive Faculty and Researcher Publications Faculty and Researcher Publications 1998-09 Real-Time Environmental Information Network and Analysis System (REINAS) Nuss, Wendell
More informationPIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:
PIPS 3.0 Ruth H. Preller Naval Research Laboratory, Code 7322 Stennis Space Center, MS 39529 phone: (228) 688-5444 fax: (228)688-4759 email: preller@nrlssc.navy.mil Pamela G. Posey NRL Code 7322 Stennis
More informationDevelopment and Application of Acoustic Metamaterials with Locally Resonant Microstructures
Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures AFOSR grant #FA9550-10-1-0061 Program manager: Dr. Les Lee PI: C.T. Sun Purdue University West Lafayette, Indiana
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationDesign for Lifecycle Cost using Time-Dependent Reliability
Design for Lifecycle Cost using Time-Dependent Reliability Amandeep Singh Zissimos P. Mourelatos Jing Li Mechanical Engineering Department Oakland University Rochester, MI 48309, USA : Distribution Statement
More informationImpact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds
Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds Ramsey R. Harcourt Applied Physics Laboratory, University of Washington,
More informationHigh Resolution Surface Characterization from Marine Radar Measurements
DISTRIBUTION STATEMENT A: Distribution approved for public release; distribution is unlimited High Resolution Surface Characterization from Marine Radar Measurements Hans C. Graber CSTARS - University
More informationINFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS
INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS AMOS 2005 TECHNICAL CONFERENCE WORKSHOP 5 September, 2005 Maui, Hawaii M. Venner AFRL, Edwards AFB, CA M. Braunstein, L. Bernstein Spectral Sciences,
More informationRogue Wave Statistics and Dynamics Using Large-Scale Direct Simulations
Rogue Wave Statistics and Dynamics Using Large-Scale Direct Simulations Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationMixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains. Barry Cobb. Alan Johnson
Mixture Distributions for Modeling Lead Time Demand in Coordinated Supply Chains Barry Cobb Virginia Military Institute Alan Johnson Air Force Institute of Technology AFCEA Acquisition Research Symposium
More informationDetermining the Stratification of Exchange Flows in Sea Straits
Determining the Stratification of Exchange Flows in Sea Straits Lawrence J. Pratt Physical Oceanography Department, MS #21 Woods Hole Oceanographic Institution, Woods Hole, MA 02543 phone: (508) 289-2540
More informationPlaying Abstract games with Hidden States (Spatial and Non-Spatial).
Playing Abstract games with Hidden States (Spatial and Non-Spatial). Gregory Calbert, Hing-Wah Kwok Peter Smet, Jason Scholz, Michael Webb VE Group, C2D, DSTO. Report Documentation Page Form Approved OMB
More informationAssimilation of Synthetic-Aperture Radar Data into Navy Wave Prediction Models
Assimilation of Synthetic-Aperture Radar Data into Navy Wave Prediction Models David T. Walker Earth Sciences Group, Veridian ERIM International P.O. Box 134008, Ann Arbor, MI 48113-4008 phone: (734)994-1200
More informationModeling the Impact of Extreme Events on Margin Sedimentation
Modeling the Impact of Extreme Events on Margin Sedimentation Jasim Imran Department of Civil and Environmental Engineering, University of South Carolina, 3 Main Street, Columbia, SC 2928. phone: (83)
More informationEstimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight
Estimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight Dale P. Winebrenner Applied Physics Laboratory, Box 355640 University of Washington Seattle, WA 98195
More informationLAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT
LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT Robert C. Beardsley, Kenneth H. Brink, Richard Limeburner Clark Laboratory, Mail Stop 21 Woods Hole Oceanographic Institution Woods
More informationModels of Marginal Seas Partially Enclosed by Islands
Models of Marginal Seas Partially Enclosed by Islands Roxana C. Wajsowicz Dept. of Meteorology University of Maryland 3433 Computer and Space Science Building College Park, MD 20852 Phone: (301) 405-5396
More informationOptimizing Robotic Team Performance with Probabilistic Model Checking
Optimizing Robotic Team Performance with Probabilistic Model Checking Software Engineering Institute Carnegie Mellon University Pittsburgh, PA 15213 Sagar Chaki, Joseph Giampapa, David Kyle (presenting),
More informationREPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationCatalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids
Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids William M. Nelson,, PhD IL Waste Management & Research Center 2003 Joint Service Scientific Conference on Chemical & Biological Defense Research
More informationGrant Number: N IP To compare obtained theoretical results with NPAL experimental data.
Coherence of Low-Frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD4, 325 Broadway,
More informationHIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING
HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING R. P. Abbott, C. D. Boley, S. N. Fochs, L. A. Nattrass, J. M. Parker, A. M. Rubenchik, J. A. Smith, and R. M. Yamamoto* University of California
More informationFleet Maintenance Simulation With Insufficient Data
Fleet Maintenance Simulation With Insufficient Data Zissimos P. Mourelatos Mechanical Engineering Department Oakland University mourelat@oakland.edu Ground Robotics Reliability Center (GRRC) Seminar 17
More informationCurating and Maintenance of the Sediment Library and Dredge Collection
Curating and Maintenance of the Sediment Library and Dredge Collection Gerard C. Bond Lamont-Doherty Earth Observatory Deep-Sea Sample Repository of Columbia University Palisades, New York 10964 Phone:
More informationSTUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR
STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR William R. Creasy EAI Corporation, 1308 Continental Drive, Suite J, Abingdon MD 21009 David J. McGarvey, Jeffrey S. Rice, Richard O'Connor,
More informationExact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling
Applied Mathematical Sciences, Vol., 8, no. 4, 77-86 Exact Solution of a Constrained Optimization Problem in Thermoelectric Cooling Hongyun Wang Department of Applied Mathematics and Statistics University
More informationSMA Bending. Cellular Shape Memory Structures: Experiments & Modeling N. Triantafyllidis (UM), J. Shaw (UM), D. Grummon (MSU)
SMA Bending Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing
More informationMTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS
MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS Mingfu Li, Hsin-Min Peng, and Chia-Shu Liao Telecommunication Labs., Chunghwa
More informationGrant Number: N IP
Coherence of Low-frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD3, 325 Broadway,
More informationComputer Simulation of Sand Ripple Growth and Migration.
Computer Simulation of Sand Ripple Growth and Migration. Douglas J. Wilson OGI School of Environmental Science and Engineering at the Oregon Health and Sciences University 20000 N.W. Walker Road, Beaverton,
More informationMetrology Experiment for Engineering Students: Platinum Resistance Temperature Detector
Session 1359 Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector Svetlana Avramov-Zamurovic, Carl Wick, Robert DeMoyer United States Naval Academy Abstract This paper
More informationCoastal Engineering Technical Note
Coastal Engineering Technical Note CETN I-48 (12/91) Evaluation and Application of the Wave Information Study for the Gulf of Mexico INTRODUCTION The Wave Information Study (WIS) for the Gulf of Mexico
More informationImproved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models
Improved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models Gerbrant Ph. van Vledder ALKYON Hydraulic Consultancy & Research P.O.Box 248, 8300 AD
More informationGeneration and Propagation of Internal Solitary Waves on the Continental Shelf and Slope
Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope Roger H.J. Grimshaw Department of Mathematical Sciences Loughborough University Loughborough, LE11 3TU, UK phone:
More informationNormic continued fractions in totally and tamely ramified extensions of local fields
Normic continued fractions in totally and tamely ramified extensions of local fields Pantelimon Stănică Naval Postgraduate School Applied Mathematics Department, Monterey, CA 93943 5216, USA; email: pstanica@nps.edu
More informationA Bound on Mean-Square Estimation Error Accounting for System Model Mismatch
A Bound on Mean-Square Estimation Error Accounting for System Model Mismatch Wen Xu RD Instruments phone: 858-689-8682 email: wxu@rdinstruments.com Christ D. Richmond MIT Lincoln Laboratory email: christ@ll.mit.edu
More informationEffects of Constrictions on Blast-Hazard Areas
Effects of Constrictions on Blast-Hazard Areas So-young Song, Jae Woon Ahn, Jong Won Park Sang Hun. Baek, Soo Won Kim, Jun Wung Lee Agency for Defence Development Taejon, Republic of Korea ABSTRACT A series
More informationFlocculation, Optics and Turbulence in the Community Sediment Transport Model System: Application of OASIS Results
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Flocculation, Optics and Turbulence in the Community Sediment Transport Model System: Application of OASIS Results Emmanuel
More informationHurricane Wave Topography and Directional Wave Spectra in Near Real-Time
Hurricane Wave Topography and Directional Wave Spectra in Near Real-Time Edward J. Walsh NASA/Goddard Space Flight Center, Code 972 Wallops Flight Facility, Wallops Island, VA 23337 phone: (303) 497-6357
More informationMarginal Sea - Open Ocean Exchange
Marginal Sea - Open Ocean Exchange Michael A. Spall Mail Stop 21 Department of Physical Oceanography Woods Hole Oceanographic Institution Woods Hole, MA 02543-1541 phone: (508) 289-3342 fax: (508) 457-2181
More informationVortex Rossby Waves and Hurricane Evolution in the Presence of Convection and Potential Vorticity and Hurricane Motion
LONG-TERM GOALS/OBJECTIVES Vortex Rossby Waves and Hurricane Evolution in the Presence of Convection and Potential Vorticity and Hurricane Motion Michael T. Montgomery Department of Atmospheric Science
More informationAn Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight
An Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight Kevin B. Smith Code PH/Sk, Department of Physics Naval Postgraduate School Monterey, CA 93943
More informationVolume 6 Water Surface Profiles
A United States Contribution to the International Hydrological Decade HEC-IHD-0600 Hydrologic Engineering Methods For Water Resources Development Volume 6 Water Surface Profiles July 1975 Approved for
More informationExperimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin
LONG TERM GOALS Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin D.K. Perovich J.A. Richter-Menge W.B. Tucker III M. Sturm U. S. Army Cold Regions Research and
More informationCoupled Ocean-Atmosphere Modeling of the Coastal Zone
Coupled Ocean-Atmosphere Modeling of the Coastal Zone Eric D. Skyllingstad College of Oceanic and Atmospheric Sciences, Oregon State University 14 Ocean Admin. Bldg., Corvallis, OR 97331 Phone: (541) 737-5697
More informationLevinson s Theorem and the Nonlocal Saito Potential
Journal of Mathematics Research February, 21 Levinson s Theorem and the Nonlocal Saito Potential S.B. Qadri Department of Physics, Ohio State University Columbus, Ohio 4321, USA Naval Research Laboratory,
More informationParameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction
Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction Ming Li Horn Point Laboratory University of Maryland Center for Environmental Science 2020 Horn Point Road, Cambridge, MD 21613
More informationPeriodic Magnetoresistance Oscillations in Side-Gated Quantum Dots
Institute of Physics Publishing Journal of Physics: Conference Series 3 () 11 119 doi:1.1/17-59/3/1/9 NPMS-7/SIMD-5 (Maui 5) Periodic Magnetoresistance Oscillations in Side-Gated Quantum Dots T. Suzuki,
More informationLISA: THE LASER INTERFEROMETER SPACE ANTENNA
LISA: THE LASER INTERFEROMETER SPACE ANTENNA Masimo Tinto Jet Propulsion Laboratory, California Institute of Technology Abstract The Laser Interferometer Space Antenna (LISA) is a deep-space mission, jointly
More informationTesting Turbulence Closure Models Against Oceanic Turbulence Measurements
Testing Turbulence Closure Models Against Oceanic Turbulence Measurements J. H. Trowbridge Woods Hole Oceanographic Institution Woods Hole, MA 02543 phone: 508-289-2296 fax: 508-457-2194 e-mail: jtrowbridge@whoi.edu
More informationVLBA IMAGING OF SOURCES AT 24 AND 43 GHZ
VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ D.A. BOBOLTZ 1, A.L. FEY 1, P. CHARLOT 2,3 & THE K-Q VLBI SURVEY COLLABORATION 1 U.S. Naval Observatory 3450 Massachusetts Ave., NW, Washington, DC, 20392-5420,
More informationComparative Analysis of Flood Routing Methods
US Army Corps of Engineers Hydrologic Engineering Center Comparative Analysis of Flood Routing Methods September 1980 Approved for Public Release. Distribution Unlimited. RD-24 REPORT DOCUMENTATION PAGE
More information