We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
|
|
- Anastasia White
- 5 years ago
- Views:
Transcription
1 We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll Research) SUMMARY We have evelope a novel approach for solving multi-frequency frequency omain wave equation. The approach is base on efficient Krylov subspace approximants an projection-base moel reuction techniques. We have consiere polynomial Krylov an extene Krylov subspaces for approximating the solution given by stability-correcte resolvent. Our numerical examples inicate that polynomial Krylov subspace allows to obtain solution for the whole a priori given frequency range at the cost of solution for minimal frequency (for that frequency range) obtaine using unpreconitione BiCGStab solver. Extene Krylov subspace has been shown to improve the convergence by proviing more uniform rate for the whole frequency range.
2 Introuction Developing fast an robust forwar-moeling methos for frequency-omain wave problems is not only a significant topic by itself, it is also of great importance for full wave inversion. When solutions for multiple frequencies are require, conventional workflow consists of solving iscretize frequencyomain problems for each frequency separately. Typically, however, iscretization gris for 3D problems consist of up to a billion noes an even with state-of-the-art preconitioners such a workflow results in rather computationally intensive tasks. Moel reuction is a well-establishe tool, allowing us to efficiently obtain solutions in the time- or frequency-omains by projecting the large-scale ynamical system on a small Krylov or rational Krylov subspace. For lossy iffusion ominate problems, for example, moel reuction has been shown to provie significant spee ups (see Zaslavsky et al. (211)). For seismic exploration, we note that problems involving lossless meia in unboune omains effectively behave as lossy ones since infinity can be viewe as an absorber of outgoing waves. In fact, two of the authors showe that these problems, polynomial Krylov subspace (PKS) moel reuction outperforms the finite-ifference time-omain metho on large time intervals (see Druskin an Remis (213)). In this paper we exten the PKS approach an use extene Krylov subspaces (EKS) for reuce-orer moel construction (Druskin an Knizhnerman (1998)). An EKS is generate by the system matrix an its inverse an reuce-orer moels taken from such a space usually converge much faster than PKS reuce-orer moels especially if a wie frequency range containing small frequencies is of interest. To compute the moels in an efficient manner, we use a moifie version of the EKS algorithm propose by Jagels an Reichel (29). Specifically, we generate a complex-orthogonal basis of the EKS via short-term recurrences by exploiting the complex-symmetric structure of the system matrix. Problem formulation Consier the multiimensional Helmholtz equation Au + ω 2 u = b. (1) In this equation, A is a self-ajoint nonnegative partial ifferential equation operator on an unboune omain that has an absolutely continuous spectrum. We note that via a proper change of variables, all frequency-omain acoustic an elastic fiel equations can be written in a form as given by Eq. (1). We now iscretize (1) using a secon-orer gri in the interior of the computational omain an use a PML for omain truncation. As a result, we obtain a matrix à N A, where à N C N N is complex symmetric. In practice, the orer of this matrix can be up to billion or even more. When the computational omain is truncate using a conventional time-omain perfectly matche layers (PML) formulation (Berenger (1994)), operator à N becomes frequency-epenent. Inee, that is sufficient for traitional preconitione solvers that treat each frequency one by one. In our approach, however, we inten to reuse the same operator for multi-frequency computations. We therefore follow Druskin an Remis (213) an Druskin et al (213) an apply a fixe-frequency PML with optimal iscrete stretching allowing low cost error control for a prescribe frequency interval. It is tempting to substitute straightforwarly the approximate operator à N in (1) an take the solution u N (ω) =(à N + ω 2 I) 1 b N as an approximation to u. Here, we note that it is out goal to solve Eq. (1) for multiple frequencies employing a fixe-frequency PML such that it is still possible to transform the frequency-omain solution back to the time-omain. However, using the non-hermitian matrix à N
3 instea of its exact Hermitian counterpart A, qualitatively changes the behavior of the solution on the complex plane. In particular, the symmetry relation u(ω)=u( ω) breaks for the approximate solution an the time-omain transform of u N (ω) is unstable. Fortunately, it is shown in Druskin an Remis (213) that we can correct for these efects by using the stabilize approximation to Eq. (1) in the form ũ N (ω)= 1 2 [ B 1 N (B N + iωi) 1 + B 1 ( N BN + iωi ) ] 1 b N, B N = ( ) 1/2 à N (2) Moel reuction an Krylov subspace methos With the stability-correcte fiel approximations available, we can now construct reuce-orer moels base on EKS in the usual way. Specifically, the moels are rawn from the Krylov subspace K m1,m 2 = span{ã m 1+1 N b N,...,à 1 N b N,b N,à N b N...,à m 2 1 N b N }. We note that the PKS K m = span{b N,à N b N,...,à m 1 N b N } correspons to K 1,m2 an PKS reuce-orer moels with application to stability-correcte solutions were investigate in Druskin an Remis (213). The avantage of such a PKS moel-orer reuction approach is that its computational costs are essentially the same as the costs of m 2 steps of the explicit finite-ifference time-omain metho or m 2 steps of the unpreconitione bi-cg metho for the single frequency Helmholtz equation. However, the reuce-orer moels taken from the PKS may not provie us with the fastest convergence if wie frequency ranges with small enough frequencies are of interest. For such problems, we therefore resort to an EKS reuce-orer moeling approach. An EKS can be seen as a special case of a rational Krylov subspace with one expansion point at zero an one at infinity. The action of à 1 N on a vector is require to generate a basis for such a space. This essentially amounts to solving a Poisson-type equation for which efficient solution techniques are available. Computing matrix-vector proucts with the inverse of the system matrix is generally still more expensive than computing matrix-vector proucts with the system matrix itself, however, an from a computational point of view we therefore prefer to eal with EKS K m1,m 2 with m 1 < m 2. In aition, the basis vectors shoul be constructe via short-term recurrence relations, since storage of all basis vectors is generally not practical for large-scale applications. In Jagels an Reichel (29), the authors evelope such an EKS algorithm in which the orthogonal bases V k(i+1) =[v,v 1,...,v i,v 1,...,v k+1,...,v ik ] for the sequence of subspaces K 1,i+1 K 2,2i+1... K k,ki+1 are inee generate via short term recurrences. Here, i is an integer that allows us to optimize the accuracy an computational costs. In this paper, we moify this algorithm an generate complex-orthogonal basis vectors of the EKS by exploiting the complex-symmetric structure of matrix à N. Denote = k(i + 1) an let V C N be matrix with columns being the generate basis vectors. Then all iterations can be summarize into the equation à N V = V H + z e T, where z = h +1, v k + h +2, v ik+1 an h ij is the (i, j) entry of matrix H. Furthermore, matrix H is a pentaiagonal matrix that satisfies D H = V T ÃNV, where D is iagonal matrix with entries δ,δ 1,...,δ i,δ 1,...,δ k+1,...,δ ik. The entries of matrix H can easily be obtaine in explicit form from the pentaiagonal matrix given in Jagels an Reichel (29). The frequency-omain EKS reuceorer moel is given by ũ (ω)= 1 [ 2 δ V B 1 (B + iωi ) 1 +V B 1 ( ) ] 1 B + iωi e 1, (3)
4 45 4 BiCGStab solver 8 7 BiCGStab solver Real part of solution Imaginary part of solution Receiver Receiver Figure 1 SEG/EAGE Salt moel. Real an imaginary parts of the solution for a frequency of 7.5 Hz compute using the BiCGStab an PKS solvers with 1,6 matrix-vector multiplications an 1,8 iterations, respectively. The solutions are almost inistinguishable. ( ) 1/2, where B = D 1/2 H D 1/2 I R is the ientity matrix, an e 1 is its first column. With the help of our moifie EKS algorithm, we now have reuce the computation of functions of a very large finite-ifference matrix à N to the action of functions of a much smaller five-iagonal matrix D 1/2 H D 1/2 times a skinny matrix V. Once both of these matrices have been compute, simulation of the frequency omain curve can be one rather cheaply. Numerical experiments First, we have consiere the 3D SEG/EAGE Salt moel an benchmarke our solver against an inepenently evelope unpreconitione BiCGStab algorithm that is a competitive conventional iterative Helmholtz solver (see Pan et al (212) for etails). Here, we took the case m 1 = 1 which correspons to PKS. Fig. 1 shows excellent agreement between two approaches. Then we note that the PKS approach requires one matrix-vector multiplication per iteration, while BiCGStab nees two. Since this part constitutes the most computationally intensive part of the iteration process, it makes sense to compare the performance of these two methos in terms of matrix-vector multiplications rather than in terms of iterations. We have plotte the convergence rates of BiCGStab against the PKS moel-orer reuction metho on Fig. 2 (left). Clearly, for this single frequency problem, both rates are rather close. We note that both approaches converge faster for higher frequencies an convergence slows own for lower frequencies. Inee, Fig. 2 (right) shows how the reuce-orer moeling metho converges for ifferent frequencies in the range from 2.5 Hz to 7.5 Hz. However, the principal ifference between the PKS metho an BiCGStab is that the former approach obtains solutions for the whole given frequency range at the convergence cost of BiCGStab for the lowest (from that range) frequency. Next, we consier what aing negative powers gives us in terms of performance. In Fig. 3 (left), we have plotte a number of convergence curves (with respect to increasing k) for approximants obtaine using EKS K k,ki+1 for ifferent values of fixe i (i = correspons to PKS) an on a frequency range running from 2.5 Hz to 7.5 Hz. As is clear from this figure, aing negative powers (i < ) visibly improves the convergence of the Krylov subspace metho. Inee, while PKS converges faster for higher frequencies, EKS improves convergence for smaller frequencies an, consequently, convergence is more uniform on the entire frequency range. To confirm that point, in Fig. 3 (right), we have plotte convergence curves for the EKS metho on a frequency range with expane lower part. As one can observe, PKS significantly slows own while EKS performs as efficient as for a narrower frequency range. Conclusions We have evelope a powerful moel-orer reuction tool for solving large-scale multi-frequency wave problems. The PKS metho allows obtaining solutions for a whole range of frequencies at the cost of
5 Convergence for 7.5Hz frequency BiCGStab for fixe frequency Iterations 2 Iterations 3 Iterations Exact solution Error 1 3 Response Number of matrix vector multiplications Frequency, Hz Figure 2 SEG/EAGE Salt moel. Convergence rates of BiCGStab an PKS at a fixe frequency of 7.5 Hz (left) an convergence behavior of PKS for a frequency range running from 2.5 Hz to 7.5 Hz (right) Frequency range 2.5Hz to 7.5Hz i=3 i=5 i=7 i=inf 1 2 Frequency range.1hz to 7.5Hz 1 L2 Error L2 Error Iteration number i=3 1 6 i=5 i=7 i=inf Iteration number Figure 3 SEG/EAGE Salt moel. Convergence of EKS for ifferent values of i an for frequencies ranging from 2.5 Hz to 7.5 Hz (left) an from.1 Hz to 7.5 Hz (right). solving a single-frequency problem with the BiCGStab solver. Moreover, EKS moel-orer reuction significantly outperform polynomial reuce-orer moeling when solutions for small frequencies are require. We also note that the EKS approach can be applie to wave problems in the time-omain. Furthermore, our approach is not limite to secon-orer schemes in the interior. In fact, since the optimal iscrete PML has spectral accuracy (see Druskin an Remis (213); Druskin et al (213)), it woul be preferable to use optimal gris (Asvaurov et al (2)) or spectral methos for interior part. Acknowlegements We thank Dr. Guangong Pan for proviing us with results for his BiCGStab finite-ifference solver. References Asvaurov, S, Druskin, V., Knizhnerman, L. [2] Application of the ifference Gaussian rules to the solution of hyperbolic problems. J. Comp. Phys., 158, Berenger, J. P. [1994] A perfectly matche layer for the absorption of electromagnetic waves. J. Comput. Phys., 114, Druskin, V., Guettel, S., Knizhnerman, L. [213] Near-optimal perfectly matche layers for inefinite Helmholtz problems. MIMS perprint , University of Manchester. Druskin, V., Knizhnerman, L. [1998] Extene Krylov subspaces: approximation of the matrix square root an relate functions. SIAM J. Matrix Anal. Appl., 19(3), Druskin, V., Remis, R. [213] A Krylov stability-correcte coorinate stretching metho to simulate wave propagation in unboune omains. SIAM J. Sci. Comput., 35(2), B376 B4. Jagels, C., Reichel, L. [29] The extene Krylov subspace metho an orthogonal Laurent polynomials. Linear Algebra Appl., 431, Pan, G., Abubakar, A., Habashy, T. [212] An effective perfectly matche layer esign for acoustic fourth-orer frequency-omain finite-ifference scheme. Geophys. J. Int., 188, Zaslavsky, M., Druskin, V., Knizhnerman, L. [211] Solution of 3D time-omain electromagnetic problems using optimal subspace projection. Geophysics, 76, F339 F351.
Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques
1 Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques From Resonators to Open Domains Rob Remis Delft University of Technology November 3, 2017 ICERM Brown University 2 Thanks
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff
More informationNONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES
NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES Vassili N. Kolokoltsov November 21, 2006 To appear in Journal of Statistical Physics. Abstract Semigroups of positivity preserving linear
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationRegularized extremal bounds analysis (REBA): An approach to quantifying uncertainty in nonlinear geophysical inverse problems
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L03304, oi:10.1029/2008gl036407, 2009 Regularize extremal bouns analysis (REBA): An approach to quantifying uncertainty in nonlinear geophysical inverse problems
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
More informationA Spectral Method for the Biharmonic Equation
A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationA simplified macroscopic urban traffic network model for model-based predictive control
Delft University of Technology Delft Center for Systems an Control Technical report 9-28 A simplifie macroscopic urban traffic network moel for moel-base preictive control S. Lin, B. De Schutter, Y. Xi,
More informationSpace-time Linear Dispersion Using Coordinate Interleaving
Space-time Linear Dispersion Using Coorinate Interleaving Jinsong Wu an Steven D Blostein Department of Electrical an Computer Engineering Queen s University, Kingston, Ontario, Canaa, K7L3N6 Email: wujs@ieeeorg
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationKNN Particle Filters for Dynamic Hybrid Bayesian Networks
KNN Particle Filters for Dynamic Hybri Bayesian Networs H. D. Chen an K. C. Chang Dept. of Systems Engineering an Operations Research George Mason University MS 4A6, 4400 University Dr. Fairfax, VA 22030
More informationConvergence of Random Walks
Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationA Constructive Inversion Framework for Twisted Convolution
A Constructive Inversion Framework for Twiste Convolution Yonina C. Elar, Ewa Matusiak, Tobias Werther June 30, 2006 Subject Classification: 44A35, 15A30, 42C15 Key Wors: Twiste convolution, Wiener s Lemma,
More information2Algebraic ONLINE PAGE PROOFS. foundations
Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationDot trajectories in the superposition of random screens: analysis and synthesis
1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationSYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS
SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationP. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,
More informationUltra-thin Acoustic Metasurface-Based Schroeder Diffuser
Ultra-thin Acoustic Metasurface-Base Schroeer Diffuser Yifan Zhu, Xuong Fan, Bin Liang *, Jianchun Cheng *, an Yun Jing * Key Laboratory of Moern Acoustics, MOE, Institute of Acoustics, Department of Physics,
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationMonte Carlo Methods with Reduced Error
Monte Carlo Methos with Reuce Error As has been shown, the probable error in Monte Carlo algorithms when no information about the smoothness of the function is use is Dξ r N = c N. It is important for
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationNon-Linear Bayesian CBRN Source Term Estimation
Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationMultiple Rank-1 Lattices as Sampling Schemes for Multivariate Trigonometric Polynomials
Multiple Rank-1 Lattices as Sampling Schemes for Multivariate Trigonometric Polynomials Lutz Kämmerer June 23, 2015 We present a new sampling metho that allows the unique reconstruction of (sparse) multivariate
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
More informationCOMPONENT mode synthesis (CMS) methods have been
AIAA JOURNAL Vol. 56, No. 1, January 018 Consiering the Higher-Orer Effect of Resiual Moes in the Craig Bampton Metho Jaemin Kim, Seung-Hwan Boo, an Phill-Seung Lee Korea Avance Institute of Science an
More informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
More informationEVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES
MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION
More informationNonlinear acoustic imaging via reduced order model backprojection
Nonlinear acoustic imaging via reduced order model backprojection Alexander V. Mamonov 1, Vladimir Druskin 2, Andrew Thaler 3 and Mikhail Zaslavsky 2 1 University of Houston, 2 Schlumberger-Doll Research
More informationInvestigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model
ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationSimple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor
DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor
More informationDeriving ARX Models for Synchronous Generators
Deriving AR Moels for Synchronous Generators Yangkun u, Stuent Member, IEEE, Zhixin Miao, Senior Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract Parameter ientification of a synchronous generator
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationTransactions on Engineering Sciences vol 5, 1994 WIT Press, ISSN
Generalize heat conuction transfer functions in two imensional linear systems of assigne geometry F. Marcotullio & A. Ponticiello Dipartimento i Energetica- Universit i L 'Aquila Monteluco i Roio - 67100
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More informationCalculus and optimization
Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function
More informationu!i = a T u = 0. Then S satisfies
Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace
More informationVIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK
AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationOptimal Control of Spatially Distributed Systems
Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist
More information6 Domain decomposition and fictitious domain methods with distributed Lagrange multipliers
M _ Thirteenth International Conference on Domain Decomposition Methos Eitors:. Debit M.Garbey R. Hoppe J. Périaux D. eyes Y. uznetsov c 2 DDM.org Domain ecomposition an fictitious omain methos with istribute
More informationData-to-Born transform for inversion and imaging with waves
Data-to-Born transform for inversion and imaging with waves Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann Arbor,
More informationExponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics
Exponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics Murthy N. Guddati North Carolina State University November 9, 017 Outline Part 1: Impedance Preserving Discretization
More informationLeft-invariant extended Kalman filter and attitude estimation
Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More information