Highly oscillatory quadrature and its applications
|
|
- Abner Banks
- 5 years ago
- Views:
Transcription
1 UNIVERSITY OF CAMBRIDGE CENTRE FOR MATHEMATICAL SCIENCES DEPARTMENT OF APPLIED MATHEMATICS & THEORETICAL PHYSICS Highly oscillatory quadrature and its applications Arieh Iserles & Syvert Nørsett
2 Report Documentation Page Form Approved OMB No Public reporting burden for the collection of information is estimated to average 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, 25 Jefferson Davis Highway, Suite 24, 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.. REPORT DATE 3 JAN REPORT TYPE N/A 3. DATES COVERED - 4. TITLE AND SUBTITLE Highly oscillatory quadrature and its applications 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 Cambridge Centre For Mathematical Sciences Department Of Applied Mathematics & Theoretical Physics 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES). SPONSOR/MONITOR S ACRONYM(S) 2. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release, distribution unlimited. SPONSOR/MONITOR S REPORT NUMBER(S) 3. SUPPLEMENTARY NOTES See also ADM749, Lie Group Methods And Control Theory Workshop Held on 28 June 24 - July 24., The original document contains color images. 4. ABSTRACT 5. SUBJECT TERMS 6. SECURITY CLASSIFICATION OF: 7. LIMITATION OF ABSTRACT UU a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified 8. NUMBER OF PAGES 33 9a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-8
3 Part I: Quadrature of HiOsc integrals Suppose that the kernel K(x, ω) L[[,] R + ] oscillates rapidly for ω, e.g. The Fourier oscillator e iωx ; The irregular oscillator e iωg(x), where g is real; The singular oscillator e iω x y / x y α, where y [,] and α < ; The Bessel oscillator J ν (ωx), ν R; The Airy oscillator Ai( ωx). We are interested in approximating I[f] = f(x)k(x, ω)dx. 2
4 What s wrong with Gaussian quadrature? Gauss Christoffel quadrature: I[f] Q GC [f] = φ(x)dx = ν l= b l f(c l )K(c l, ω), where φ is the (ν )-degree polynomial that interpolates f(x)k(x, ω) at the quadrature nodes c < c 2 < < c ν in [,]. In particular: c, c 2,..., c ν zeros of P ν (2x ) the quadrature is of maximal order 2ν: this is the Gauss Legendre quadrature. Suppose for simplicity that K(x, ω) = e iωx. For ω and fixed x the value e iωx is, to all intents and purposes, a random number on the complex unit circle z =. Therefore, for fixed ν Q GC [f] = ν l= b l f(c l )e iωc l O(), ω. On the other hand, Riemann Lebesgue implies that lim ω I[f] =, f L [,]. 3
5 e(+iω)x dx = e+iω + iω 3.5 ν=2 ν=3 ν= ω ω ω 3.5 ν= ν=5 ν= ω ω ω 4
6 Generalising the Filon method Instead of interpolating the integrand f(x)k(x, ω) at the quadrature nodes, we interpolate the values of f(x) there by the polynomial φ: I[f] Q F [f]= φ(x)k(x, ω)dx = ν l= b l (ω)f(c l ). Note that the weights depend on the frequency ω. Filon Legendre: ν = 2, c = [ 2 6 3, ] ω Q F [e x ] I[e x ] ω Q F [e x ] I[e x ] ω Just two quadrature points... 5
7 ... but we can do even better! Choose instead Lobatto points: ν = 2, c = [,]: ω Q F [e x ] I[e x ] ω 2 Q F [e x ] I[e x ] ω The wonder-method is f(x)eiωx dx b (ω)f() + b 2 (ω)f(), where b (ω)= iω + eiω ( iω) 2, b 2 (ω)= eiω iω eiω ( iω) 2. But why does it work so well? 6
8 Asymptotic expansion Let K(x, ω) = e iωg(x), where g is real and smooth. In addition, we require that g in [,]. Integrating by parts, I[f]= iω = iω f(x) g (x) de iωg(x) dx [ iωg() f() g () iω I[d(f/g )/dx]. We continue by induction. Let m= e dx eiωg() f() g () σ (x)=f(x), σ m+ (x)= d σ m (x) dx g (x), m Z +. Then, in the limit, [ σm () I[f] ( iω) m+ g () eiωg() σ m() g () eiωg() ] ]. 7
9 We have and so on. σ =f, σ = g g 2 f + g f, σ 2 = 3g 2 gg g 4 Asymptotic quadrature Let s Q A s [f]= m= ( iω) m+ f 3 g g 3 f + g 2 f [ σm () g () eiωg() σ m() g () eiωg() The method uses s derivatives of f and Q A s [f] I[f] O( ω s ), ω. For g(x) = x we have s Q A s [f] = m= ( iω) m+[f(m) () e iω f (m) ()]. ]. 8
10 Filon-type methods Given nodes c, < < c ν and n,..., n ν N, we choose the unique polynomial φ of degree l n l such that for all l =,..., ν φ (j) (c l ) = f (j) (c l ), j =,..., n l. A Filon-type method is Q F [f] = I[ φ] ω ω 3 Q F [e x ] I[e x ] ω BLUE: Q A 2 ; GREEN: Q F, ν = 2, n = n 2 = 2, c =, c 2 = ; MAGENTA: Q F, ν = 3, n = n 3 = 2, n 2 =, c =, c 2 = 2, c 3 =. 9
11 THEOREM If c =, c ν = and n = n ν = s then Q F [f] I[f] O ( ω s ) when ω. Proof Since Q F [f] I[f] = I[ φ f] and φ (j) () = f (j) (), φ (j) () = f (j) (), for j =,,..., s, the proof follows from the asymptotic expansion of I[ φ f]. Thus, Filon has the same asymptotic order as the asymptotic method. Typically it has a smaller error constant, which can be further decreased, by adding extra nodes in (,). All this is true as long as there are no stationary points of the oscillator in [,], i.e. g in the interval.
12 Stationary points Suppose first that g (y) =, g (y), for some y [,] and that g elsewhere. Naive integration by parts breaks down, since division by g introduces a polar singularity. An alternative is the method of stationary phase (Cauchy, Stokes, Kelvin), except that, while requiring nasty contour integration, it does not deliver all the information we need. Instead, let and µ (ω) = eiωg(x) dx de iωg(x) I[f]=f(y)µ (ω) + f(x) f(y) iω g (x) dx =f(y)µ (ω) + [ e iωg()f() f(y) iω g () ] e iωg()f() f(y) g () [ ] d f(x) f(y) iω dx g e iωg(x) dx. (x) dx
13 We continue by induction. Letting we have σ (x)=f(x), σ m+ (x)= d σ m (x) σ m (y) dx g, m N, (x) I[f] µ (ω) + m= m= ( iω) m+ σ m (y) ( iω) m [ e iωg()σ m() σ m (y) g () e iωg()σ m() σ m (y) g () ]. The van der Corput lemma µ (ω) = O ( ω /2). Therefore, using the first s derivatives at, y and gives an asymptotic method with an asymptotic error of O ( ω s 3 2 ). Easy generalisation to several stationary points and to g (y) = = g (r) (y) =, g (r+) (y), r. 2
14 Filon again... The Filon method can be generalised to cater for stationary points. Again, the idea is to interpolate to f and its derivatives at {, y, y 2,..., y n,}, where y,..., y n are the stationary points g(x) = x( x) y= 2, s=.2 ω 3/2 Q A [e x ] I[e x ] ω ω 3/2 Q F [e x ] I[e x ] c =, c 2 = 4, c 3 = 2,.5.6 c 4 = 3 4, c 5 =..4 ω 3/2 Q F [e x ] I[e x ].5 c =, c 2 = 2, c 3 = ω ω 3
15 Exotic oscillators Most results can be extended to more exotic oscillators, for example J ν (ωx) and Ai( ωx) but, clearly, much remains to be done. Computation of special functions In the stage of tentative ideas: using Filon quadrature for fast computation of special functions (e.g. hypergeometric and Bessel functions) for large arguments, (hopefully) more precise than using standard asymptotic formulæ. Singular integrals (Hermann Brunner, AI & SPN) Similar techniques have been applied to the kernels e iω x y x y α and e iω x y log x y, where α < and y [,]. The asymptotic expansion is more difficult, mainly since σ m need not be smooth at y, but an important observation is that singularities play similar role to stationary points. 4
16 Multivariate integrals This is perhaps the most fascinating chapter of our work! RESULT Let Ω R d be a compact domain with piecewise-linear boundary, Ω Then, as long as g(x) in clω, Ω f(x)eiωg(x) dv m= ( iω) m+d k a m,k [f], where each functional a m,k, a periodic function in ω, depends just on f and its first m derivatives at the kth vertex. 5
17 Consequently, we again have two options, both resulting in an O ( ω s d) quadrature: either truncate the asymptotic expansion at m = s or replace f by Hermite interpolation to the function and its first s directional derivatives at the vertices. RESULT 2 Let Ω R d be a compact domain with piecewise-smooth boundary and without cusps. Suppose again that there are no critical points in clω, i.e. that g there. Then = iω Ω f(x)eiωg(x) dv iω Ω f(x) g(x) 2n(x) g(x)e iωg(x) ds Ω f(x) g(x) 2 g(x)eiωg(x) dv, where n is the outward unit normal. This can be converted into an asymptotic expansion, a Stokes-type theorem, pushing the integral from Ω to the boundary. All this can be extended to cater for nondegenerate critical points x Ω, where g(x ) =, det g(x ). 6
18 Part II: Fredholm equations of the 2nd kind (Hermann Brunner, AI & SPN) Consider the problem K[φ](y) = λφ(y) g(y), y [,], where g is given, and λ σ(k). K[φ](y) = φ(x)eiω x y dx A naive approach Cover [,] with the grid = y < y < < y N < y N =. Let φ k φ(y k ) and replace integrals with Filon. We obtain a linear system of the form N l= b k,l (ω)φ l = λφ k g k, k =,,..., N. This will not work, since the solution φ also oscillates with frequency ω, and this means that our asymptotics break down. 7
19 An alternative We seek complex numbers λ m and complex-valued functions φ m s.t. Since K[φ](y)= we have But K[φ m ] = λ m φ m. y φ(x)eiω(y x) dx+ y φ(x)eiω(x y) dx, dk[φ](y) =iω dy d 2 K[φ](y) dy 2 dk[φ](y) dy [ y φ(x)eiω(y x) dx ] y φ(x)eiω(x y) dx =(iω) 2 K[φ](y) + 2iωφ(y). = λφ (y), d 2 K[φ](y) dy 2, = λφ (y). 8
20 Therefore λφ = (iω) 2 λφ + 2iωφ. Let then θ(ω) = ω 2 2iω λ, Moreover, φ + θ 2 φ =. λφ ()= dk[φ]() dy λφ ()= dk[φ]() dy = iωλφ(), = iωλφ(). The condition at y = results (up to normalization) in φ(x) = (θ ω)e iθx + (θ + ω)e iθx. Note however that θ depends on the unknown eigenvalue λ. 9
21 Using the boundary condition at y = we obtain (θ ω) 2 e iθ = (θ + ω) 2 e iθ. Therefore (θ ω)e 2 iθ = ±(θ + ω)e 2 iθ, Taking the plus sign we obtain the transcendental equation iθtan θ 2 = ω, while the minus sign yields iθcot θ 2 = ω. The solutions of these equations interlace: the first has a solution for Re θ (2mπ,(2m + )π) and the second in Re θ ((2m + )π,(2m + 2)π). We observe that The real part of θ behaves like O(m), The imaginary part of θ is O() and small. 2
22 The values of θ m /(2π) in the complex plane for m 2. How is this going to help? Let (with greater generality) K[f](y) = f(x)k(x, y)dx, y [,]. Then the Hilbert Schmidt theory tells us that K has a countable number of distinct eigenvalues and eigenfunctions {λ m, φ m }. Let f, g = f(x)g(x)dx be the standard real L 2 inner product. 2
23 For m n λ m φ m (y) = λ m φ m, φ n = By symmetry, also φ m(y)k(x, y)dx φ m(x)φ n (x)k(x, y)dxdy. λ n φ m, φ n = φ m(x)φ n (x)k(x, y)dxdy and we just deduce L 2 orthogonality of the eigenfunctions, φ m, φ n =, m n. Note that, is not a positive definite inner product: it is complex-valued and it is entirely possible that f, f =, f. However, φ m, φ m =, and that s all we need. 22
24 A spectral method We expand Therefore, where f(y) = f m = m= f m φ m (y). g m λ λ m, m, g m = g, φ m φ m, φ m = g, φ m 2(λ 2 m 2iω ω2 ). We thus need to compute g, φ m for a large number of ms. However, if then g, φ m =(θ m ω) θ m = α m iβ m + (θ m + ω) g(x)e(β m+iα m )x dx g(x)e (β m+iα m )x dx. 23
25 Recall: while α m 2πm is large, β m is small. Moreover, g is nonoscillatory. Therefore All the integrals can be computed very fast and accurately by either the asymptotic method or a Filon-type method. An ongoing challenge is to generalize all this to other Fredholm kernels, e.g. K[f](y) = for γ > and f(x)xγ e iω x y dx for γ (,). K[f](y) = f(x)eiω x y x y γdx 24
26 Part III: Solving HiOsc differential equations We commence from the linear ODE y = A(t)y, t, y() = y. Suppose that its solution oscillates fast, e.g. that all the eigenvalues of A live in cl C and there are large eigenvalues on ir. Standard numerical methods perform very poorly, the reason being that the principal error term of a p th-order classical method is of the form h p+ D p+ (y(t N )), where D p+ is a linear combination of elementary differentials of order p +. y(t) oscillates with frequency ω y (p+) (t) ω p+ y(t), hence D p+ is very large! 25
27 An alternative: Change of variables To time-step from t N to t N+ = t N + h, set y(t) = e (t tn)ã x(t t N ), t t N, where à = A(t N+ ). Then 2 where x = B(t)x, t, x() = y N, B(t) = e tã [A(t) Ã]e tã. Since e ±tã oscillates rapidly, so does B(t). We have already seen that high oscillation can be turned to our advantage. The main idea is to invert the reason for the failure of classical methods: Integrate, don t differentiate! Specifically, for an s-fold integral and B(x) a product of s terms from {B(x ),..., B(x s )}, B(x)dx s dx O ( ω s). 26
28 The Magnus method Letting we have Ω(t)= t B(x)dx t 2 t + 4 t x x(t) = e Ω(t) x, [B(x 2), B(x )]dx 2 dx x x2 x x Thus, repeated integration.... [[B(x 3), B(x 2 )], B(x )]dx 3 dx 2 dx [B(x 3),[B(x 2 ), B(x )]]dx 3 dx 2 dx An advantage of Magnus: If A(t) lives in a Lie algebra g then y(t) evolves on a homogeneous space M, acted upon by the corresponding Lie group G. Using Magnus (with or without change of variables) ensures y N M, N. 27
29 An example: The Airy equation Let y + ty =. h=/ Runge Kutta 4 y N y(t N ) O(h 4 t 3/4 ) h=/ h=/ h=/ RK Gauss Legendre 4 y N y(t N ) O(h 4 t 3/4 ).5 h=/ h=/
30 x x x x x x Modified Magnus, h = 4, y N y(t N ) O(h 3 t /4 ) 29
31 Calculating integrals: We use Filon: all (multivariate) integrals can be calculated to high precision using just B() and B(h). As before, high oscillation helps computation! A disadvantage of Magnus: We need to calculate two exponentials per step: e ±hã and e Ω(h). This can be problematic when the dimension is large. However, while Ω is typically unstructured, this is not the case with Ã. Suppose that the ODE originates in a semidiscretized PDE. Then often à is block Toeplitz and e ±Ã can be calculated very fast by FFT. The challenge is thus to do away with the need for the calculation of e Ω(h). 3
32 The Neumann method To avoid the calculation of the second exponential, we abandon Magnus in favour of the Neumann expansion x(t) = where N (t) I and t x m= N k (t)y N, xm N m (t) = B(x ) B(x m )dx m dx. Because of high oscillation, N m (h) O((h/ω) m ), hence very rapid convergence. Multivariate integrals can be computed in a very small number of function evaluations, similarly to Magnus integrals. Again, high oscillation of B means that Filon methods are very precise. Numerical results for Airy are virtually identical to Magnus, but the method comes into its own for HiOsc PDEs, e.g. the Schrödinger equation. 3
33 Nonlinear equations Suppose that y = A(t)y + g(y) is highly oscillatory. Transforming as before, but letting à = A(t N ), we have Let x = B(t)x + e tã g(t N + t,e tã x). Φ = B(t)Φ, t, Φ() = I. Note that we can evaluate Φ by either Magnus or Neumann. Then x(t)=φ(t)y N t + Φ(t ξ)eξã g(t N + ξ,e ξã x(ξ))dξ. This motivates the waveform relaxation approach, x [] (t) y N, x [m+] (t)=φ(t)y N + t Φ(t ξ)eξã g(t N + ξ,e ξã x [m] (ξ))dξ. 32
34 Next steps... Filon without derivatives: Work in progress. Letting nodes depend on ω, it is possible to obtain arbitrary degree of error attenuation without using derivatives; Exotic oscillators: For starters, how to compute f(x)sin(ω sin πx)dx? Multivariate HiOsc integrals: What are the implications of the Stokes-type theorem? Volterra HiOsc equations: The current approach doesn t scale up e.g. to singular kernels; HiOsc PDEs: Much further work required for specific PDEs, e.g. Schrödinger and Hamilton Jacobi; Stochastic DEs: Perhaps... 33
Report 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 informationDiagonal 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 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 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 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 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 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 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 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 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 informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
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 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 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 informationMaximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers
Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers Curtis R. Menyuk based on the PhD dissertation of: Dr. Jonathan Hu now at Princeton University 1 Report Documentation
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 informationCalculation of Highly Oscillatory Integrals by Quadrature Method
Calculation of Highly Oscillatory Integrals by Quadrature Methods Krishna Thapa Texas A&M University, College Station, TX May 18, 211 Outline Purpose Why Filon Type Quadrature Methods Our Expectations
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationA Numerical Approach To The Steepest Descent Method
A Numerical Approach To The Steepest Descent Method K.U. Leuven, Division of Numerical and Applied Mathematics Isaac Newton Institute, Cambridge, Feb 15, 27 Outline 1 One-dimensional oscillatory integrals
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 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 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 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 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 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 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 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 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 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 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 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 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 informationNUMERICAL SOLUTIONS FOR OPTIMAL CONTROL PROBLEMS UNDER SPDE CONSTRAINTS
NUMERICAL SOLUTIONS FOR OPTIMAL CONTROL PROBLEMS UNDER SPDE CONSTRAINTS AFOSR grant number: FA9550-06-1-0234 Yanzhao Cao Department of Mathematics Florida A & M University Abstract The primary source of
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 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 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 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 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 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 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 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 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 informationLong-Range Underwater Sound Propagation: Environmental Variability, Signal Stability and Signal Coherence
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Long-Range Underwater Sound Propagation: Environmental Variability, Signal Stability and Signal Coherence Michael G. Brown
More informationReduced Basis and Stochastic Modeling of Liquid Propellant Rocket Engine as a Complex System W. A. Sirignano, A. Sideris, S. Menon, R. Munipalli, D.
Reduced Basis and Stochastic Modeling of Liquid Propellant Rocket Engine as a Complex System W. A. Sirignano, A. Sideris, S. Menon, R. Munipalli, D. Ota, D. R. Kassoy The treatment of combustion and flow
More informationThe Mechanics of Bubble Growth and Rise in Sediments
The Mechanics of Bubble Growth and Rise in Sediments PI: Bernard P. Boudreau Department of Oceanography Dalhousie University Halifax, Nova Scotia B3H 4J1, Canada phone: (902) 494-8895 fax: (902) 494-3877
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 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 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 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 informationContract No. N C0123
Particle Size Distribution and Optical Volume Scattering Function in the Mid and Upper Water Column of Optically Deep Coastal Regions: Transport from the Bottom Boundary Layer Y. C. Agrawal Sequoia Scientific,
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 informationHYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA
HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling By BIROL KARA, ALAN WALLCRAFT AND JOE METZGER Naval Research Laboratory, Stennis Space Center, USA MURAT GUNDUZ Institute
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 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 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 informationANALYSIS OF LOCALIZED HIGH MAGNETIC SUSCEPTIBILITY ZONES AT JEFFERSON PROVING GROUND, INDIANA
ANALYSIS OF LOCALIZED HIGH MAGNETIC SUSCEPTIBILITY ZONES AT JEFFERSON PROVING GROUND, INDIANA Ryan E. North*, Eric W. Smith, and Jose L. Llopis U. S. Army Engineer Research and Development Center Geotechnical
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 informationMass Transport by Second Mode Internal Solitary Waves
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Mass Transport by Second Mode Internal Solitary Waves Alan Brandt, PI Johns Hopkins Univ. Applied Physics Laboratory Laurel,
More informationNonlinear Fourier Analysis of Shallow Water Waves
Nonlinear Fourier Analysis of Shallow Water Waves A.R. Osborne Dipartimento di Fisica Generale, Università di Torino Via Pietro Giuria 1 10125 Torino, Italy phone: (390) 11-670-7455 or (390) 11-822-4739
More informationDISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Understanding Impacts of Outflow on Tropical Cyclone Formation and Rapid Intensity and Structure Changes with Data Assimilation
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 informationEFFECTS OF LOCAL METEOROLOGICAL VARIABILITY ON SURFACE AND SUBSURFACE SEISMIC-ACOUSTIC SIGNALS
EFFECTS OF LOCAL METEOROLOGICAL VARIABILITY ON SURFACE AND SUBSURFACE SEISMIC-ACOUSTIC SIGNALS Jason R. McKenna and Mihan H. McKenna U.S. Army Engineer Research & Development Center, Vicksburg, MS 39180
More informationHIGH PERFORMANCE CONTROLLERS BASED ON REAL PARAMETERS TO ACCOUNT FOR PARAMETER VARIATIONS DUE TO IRON SATURATION
HIGH PERFORMANCE CONTROLLERS BASED ON REAL PARAMETERS TO ACCOUNT FOR PARAMETER VARIATIONS DUE TO IRON SATURATION Jorge G. Cintron-Rivera, Shanelle N. Foster, Wesley G. Zanardelli and Elias G. Strangas
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 informationINFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS)
AFRL-MN-EG-TP-2006-7403 INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS) C. Michael Lindsay, National Research Council, Post Doctoral Research Associate Mario E. Fajardo
More informationMine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)
Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST) Marcelo H. Garcia Department of Civil and Environmental Engineering University of Illinois at Urbana-Champaign 205 North Mathews
More informationExpressions for the Total Yaw Angle
ARL-TN-0775 SEP 2016 US Army Research Laboratory Expressions for the Total Yaw Angle by Benjamin A Breech NOTICES Disclaimers The findings in this report are not to be construed as an official Department
More informationPredicting Tropical Cyclone Formation and Structure Change
Predicting Tropical Cyclone Formation and Structure Change Patrick A. Harr Department of Meteorology Naval Postgraduate School Monterey, CA 93943-5114 Telephone: (831)656-3787 FAX:(831)656-3061 email:
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 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 information