Domain Decomposition and Model Reduction of Systems with Local Nonlinearities
|
|
- Cory Warren
- 5 years ago
- Views:
Transcription
1 Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation moel reuction an omain ecomposition to erive reuce orer moels with guarantee error bouns for systems of iscretize partial ifferential equations (PDEs) with a spatially localize nonlinearities. Domain ecomposition techniques are use to ivie the problem into linear subproblems an small nonlinear subproblems. Balance truncation is applie to the linear subproblems with inputs an outputs etermine by the original in- an outputs as well as the interface conitions between the subproblems. The potential of this approach is emonstrate for a moel problem. Introuction Moel reuction sees to replace a large-scale system of ifferential equations by a system of substantially lower imension that has nearly the same response characteristics. This paper is concerne with moel reuction of systems of iscretize partial ifferential equations (PDEs) with spatially localize nonlinearities. In particular, we are intereste in constructing reuce orer moels for which the error between the input-to-output map of the original system an the input-to-output map of the reuce orer moel can be controlle. Balance truncation is a particular moel reuction technique ue to [6], which for linear time invariant systems leas to reuce orer moels which approximate the original input-to-output map with a user controlle error [, 6]. Although extensions of balance truncation to nonlinear systems have been propose, see, e.g., Kai Sun, Matthias Heinenschloss, Danny C. Sorensen Department of Computational an Applie Mathematics, MS-34, Rice University, 6 Main Street, Houston, Texas 77-89, USA. leinsun/heinen/sorensen@rice.eu Rolan Glowinsi University of Houston, Department of Mathematics, 6 P. G. Hoffman Hall, Houston, Texas , USA. rolan@math.uh.eu
2 Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen [9, 4], there are no bouns available for the error between the input-to-output map of the original system an that of the reuce orer moel. Proper Orthogonal Decomposition (POD) is often use for moel reuction of nonlinear systems. Error bouns are available for the error between the so-calle snapshots an the reuce orer moel, see, e.g., [, 3], but no bouns for the error between the input-tooutput map of original system an that of the reuce orer moel, unless the socalle snapshot set reflects all possible inputs. Our approach uses omain ecomposition techniques to ivie the problem into linear subproblems an small nonlinear subproblems. Balance truncation is applie only to the linear subproblems with inputs an outputs etermine by the original inan outputs as well as the interface conitions between the subproblems. We expect that this combination of omain ecomposition an balance truncation leas to a substantial reuction of the original problem if the nonlinearities are localize, i.e., the nonlinear subproblems are small relative to the other subomains, an if the interfaces between the subproblems are relatively small. To eep our presentation brief, we consier a moel problem which couples the D Burgers equation to two heat equations. This is motivate by problems in which one is primarily intereste in a nonlinear PDE which is pose on a subomain an which is couple to linear PDEs on surrouning, larger subomains. The linear PDE solution on the surrouning subomains nees to be compute accurately enough to provie acceptable bounary conitions for the nonlinear problem on the inner subomain. Such situations arise, e.g., in regional air quality moels. Our wor is also relate to [4], which is an example paper which iscusses the coupling of linear an nonlinear PDEs, but no imension reuction is applie. Domain ecomposition an POD moel reuction for flow problems with moving shocs are iscusse in []. POD moel reuction is applie on the subomains away from the shoc. The paper [8] iscusses a ifferent moel reuction technique for secon orer ynamical systems with localize nonlinearities. The papers [, ] an [] iscuss ifferent moel reuction an substructuring techniques for secon orer ynamical systems an moel reuction of interconnect systems respectively. The Moel Problem Let Ω = 3 = Ω, where Ω = (, ), Ω = (,) an Ω 3 = (,) an let T > be given. Our moel problem is given by y ρ t (x,t) µ y (x,t) = S (x,t), (x,t) Ω (,T ), (a) y (x,) = y (x), x Ω, =,3, (b) y (,t) =, y 3 (,t) = t (,T ), (c)
3 Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities 3 y ρ t (x,t) µ y + y y (x,t) =, (x,t) Ω (,T ), () with the following interface conitions y (x,) = y (x), x Ω, (e) y (,t) = y (,t), y (,t) = y 3 (,t), t (,T ), (a) y µ (,t) = µ y (,t), µ y (,t) = µ y 3 3 (,t), t (,T ). (b) We assume that the forcing functions S, S 3 are given by S = n s i= b i (x)u i (t), =,3. (3) To obtain the wea form of () an (), we multiply the ifferential equations (a, ) by test functions v i H (Ω i ), i =,,3, respectively, integrate over Ω i, an apply integration by parts. Using the bounary conitions (c, h) this leas to ρ y v x + µ t Ω ρ y v x + µ t Ω y Ω Ω y v x µ y v x + v = S v x, =,3, (4a) Ω Ω y Ω y y v x µ v =. (4b) If v H (Ω ), =,3, satisfy v ( ) =, v 3 () =, then (c), (4a) imply y ( ) µ y 3 () µ 3 = y S v x + ρ y v x + µ Ω t Ω Ω y 3 y 3 v 3 x µ 3 Ω 3 Ω 3 = S 3 v 3 x ρ 3 Ω 3 t v x, v 3 x. (a) (b) If v H (Ω ) satisfies v ( ) = an v () =, then (4b) implies µ y ( ) y v = ρ y v x µ t Ω Ω x y Ω y v x. Finally, if v H (Ω ) satisfies v ( ) = an v () =, then (4b) implies µ y () y v = ρ y v x + µ t Ω Ω x + y Ω y v x. (c) () The ientities () are use to enforce the interface conitions (). We iscretize the ifferential equations in space using piecewise linear functions. We subivie Ω j, j =,,3, into subintervals. Let x i enote the subinterval enpoints an let v i be the piecewise linear basis function with v i (x i ) = an v i (x j ) = or all j i. We efine the following inex sets
4 4 Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen I I = {i : x i [, )}, I I = {i : x i (,)}, I I 3 = {i : x i (,]}, I Γ = {i : x i = }, I Γ 3 = {i : x i = }. Given y i for i I Γ IΓ 3, we compute functions y (t,x) = y i (t)v i (x) + y i (t)v i (x), =,3, (6a) i I I i I Γ y (t,x) = i I I y i (t)v i (x) + y i (t)v i (x) + y i (t)v i (x), i I Γ i I3 Γ where in (6a) we use I Γ = I Γ if = an IΓ = I3 Γ if = 3, as solutions of ρ y v i x + µ t Ω Ω y x v ix = S v i x,i I I, =,3, Ω ρ y v i x + µ t Ω y x v ix + y y v i x =,i I. I Ω Ω (6b) If we set y I = (y i) i I I, =,,3, y Γ j = (y i) i I Γ j, j {,3}, y Γ = (y Γ,yΓ 3 )T, an u = (u i ) i=,...,ns, =,3 (cf. (3)), the previous ientities can be written as M II M II t yi + A II y I + M IΓ t yi + A II y I + M IΓ M II 3 t yi 3 + A II 3 y I 3 + M IΓ 3 t yγ + A IΓ y Γ = B I u, (7a) t yγ + A IΓ y Γ + N I (y I,y Γ ) =, (7b) t yγ 3 + A IΓ 3 y Γ 3 = B I 3u 3. (7c) By construction, the functions y j, j =,,3, in (6) satisfy (a). To enforce (b) we insert the ientities (), (6) into (b). The resulting conitions can be written as M Γ I M Γ I 3 t yi + A Γ I y I + (M Γ Γ + M Γ Γ ) t yγ + (A Γ Γ + A Γ Γ )y Γ (8a) +M Γ I t yi + A Γ I y I + N Γ (y I,y Γ ) = B Γ u, (8b) t yi 3 + A Γ 3 I y I 3 + (M Γ 3 Γ + M Γ 3 Γ ) t yγ 3 + (A Γ 3 Γ + A Γ 3 Γ )y Γ 3 (8c) +M Γ 3 I t yi + A Γ 3 I y I + N Γ 3(y I,y Γ 3) = B Γ 3 u 3. (8) To summarize, our iscretization of () an () is given by (7) an (8). As outputs we are intereste in the solution of the PDE at the spatial locations ξ =,ξ =,ξ 3 =. Thus the output equations are y (t,ξ ) = i I I y i (t)v i (ξ ), =,,3, which can be written as z I (t) = CI jy I (t), where CI R II, =,,3.
5 Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities 3 Balance Truncation Moel Reuction Given E R n n symmetric positive efinite, A R n n, B R n m, C R q n, an D R q m, we consier linear time invariant systems in state space form E t y(t) = A y(t) + Bu(t), t (,T ), y() = y, (9a) z(t) = C y(t) + Du(t), t (,T ). (9b) Projection methos for moel reuction generally prouce n r matrices V,W with r n an with W T E V = I r. One obtains a reuce form of equations (9) by setting y = V ŷ an projecting (imposing a Galerin conition) so that W T [E V ŷ(t) AV ŷ(t) Bu(t)] =, t (,T ). t This leas to a reuce system of orer r with matrices Ê = W T E V = I, A = W T AV, B = W T B, C = C V, an D = D. Balance reuction is a particular techniqe for constructing the projecting matrices V an W, see, e.g., [, 6]. One first solves the controllability an the observability Lyapunov equation A PE + E PA T + BB T = an A T QE + E QA + C T C T =, respectively. Uner the assumptions of stability, controllability an observability, the matrices P,Q are both symmetric an positive efinite. There exist methos to compute (approximations of) P = UU T an Q = LL T in factore form. In the large scale setting the factorization is typically a low ran approximation. See, e.g., [8, 7]. The balancing transformation is constructe by computing the singular value ecomposition U T E L = ZSY T an then setting W = UZ r, V = LY r, where S r = iag(σ,σ,...,σ r ) is the r r submatrix of S = S n. The singular values σ j are in ecreasing orer an r is selecte to be the smallest positive integer such that σ r+ < τσ where τ > is a prespecifie constant. The matrices Z r,y r consist of the corresponing leaing columns of Z,Y. It is well nown [6] that A must be stable an that for any given input u we have z ẑ L u L (σ r σ n ), () t yγ, AIΓ where ẑ is the output (response) of the reuce moel. Moel reuction techniques for infinite imensional systems are reviewe in, e.g., [3]. We want to apply balance truncation moel to the linear subsystems an 3 in (7) an (8). We nee to ientify the input-output relations for these subsystems in the context of the couple system to ensure that balancing techniques applie to these subsystems leas to a reuce moel for the couple system with error bouns. To ientify the appropriate input-output relations, we focus on subsystem. Examination of (7a,b) an (8a) shows that M IΓ yγ an BI u are the inputs into system an C Γ I are the outputs. Hence, if yi, MΓ I t yi + AΓ IyI
6 6 Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen then we nee to apply moel reuction to M II M IΓ = an M Γ I =, () t yi = A II y I A IΓ y Γ + B I u (a) z I = C I y I, z Γ = A Γ I y I. (b) The system () is exactly of the form (9) an we can apply balance truncation moel reuction to obtain M II t ŷi = ÂII ẑ I = ĈI ŷi, ŷi ÂIΓ y Γ + B I u (3a) ẑ Γ = ÂΓ I ŷi. (3b) Subsystem 3 can be reuce analogously. The reuce moel for the couple nonlinear system (7) an (8) is now obtaine by replacing the subsystem matrices for subsystems an 3 by their reuce matrices. Whether the balance truncation error boun () can be use to erive an error boun between the original couple problem (7) an (8) an its reuce moel is uner investigation. In our finite element iscretization we use mass lumping to obtain (). However other iscretizations, such as spectral elements or iscontinuous Galerin methos satisfy () irectly, see [, ]. 4 Numerical Results We subivie Ω j into equiistant subintervals of length h = /N, =,,3, an we use piecewise linear basis functions.the size of the system (7), (8) is 9(N + N 3 )+N +. The parameters in the PDE are ρ =, =,,3, an µ =., µ =., µ =.. For subsystem an 3 we compute low-ran approximate solutions of the controllability an observability Lyapunov equations using the metho escribe in [8]. We truncate such that σ r+ < τσ, where τ = 4. The sizes of the full an of the reuce orer moels for various iscretization parameters are shown in Table. The subsystems an 3 reuce substantially an the size of the subsystem limits the amount of reuction achieve overall. For example, for N = N 3 = the subsystems an 3 are each reuce in size from 8 to. The size of the couple system is reuce from 36 + N to 3 + N. Next, we compare the system output given forcing functions S = u (t), S 3 = u 3 (t) (cf., (3)) with u (t) = sin(3t)(.8t/t ), u 3(t) = sin(t)(.3 +.7t/T ) on (,T ) = (,). The full orer moel (7), (8) an the corresponing reuce orer moel are solve using the moifie θ-scheme [7, 9] with (macro) time step t = T /. Figure shows the outputs, i.e., the approximate solution of the PDE at ξ =,ξ =,ξ 3 =. The left plot in Figure shows the solution of the reuce orer iscretize PDE. The solution of the iscretize PDE is visually inistinguishable
7 Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities 7 Table Dimension of the full an of the reuce orer moels for various iscretization parameters N, N, N3 an τ = 4. N = N3 N size of full orer moel size of reuce orer moel Fig. Outputs,, 3 of the full orer system corresponing to the iscretization N = N = N3 = are given by, an, respecitively. Outputs,, 3 of the reuce orer system are given by otte, ashe an soli lines, respectively. 4 t $!!$ t!! x t!! x Fig. Solution of the reuce orer iscretize PDE (left) an error between the solution of the iscretize PDE an the reuce orer system (right) for iscretization N = N = N3 =. from the solution of the reuce orer iscretize PDE, as inicate by the size of the error shown in the right plot in Figure. The error is larger in the right subomain because the PDE solution is positive an the avection term in () avects the solution to the right. Our numerical results inicate that the coupling of balance truncation reuction for linear time variant subsystems with spatially localize nonlinear moels leas to a couple reuce orer moel with an error in the input-to-output map that is comparable to the error ue to balance truncation moel reuction applie to the linear subsystems alone. The efficiency of the approach epens on the size of the interface an on the size of the localize nonlinearity. Investigations for higher imensional problems are unerway to explore the overall gains in efficiency.
8 8 Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Acnowlegements This research was supporte in part by NSF grant ACI-36 an AFOSR grant FA References. Antoulas, A.C.: Approximation of Large-Scale Systems. SIAM, Philaelphia (). Bai, Z., Liao, B.S.: Towars an optimal substructuring metho for moel reuction. In: J. Dongarra, K. Masen, J. Wasniewsi (es.) Applie Parallel Computing, pp Springer Lecture Notes in Computer Science, Vol. 373 (6) 3. Curtain, R.F.: Moel reuction for control esign for istribute parameter systems. In: Research irections in istribute parameter systems (Raleigh, NC, ), Frontiers Appl. Math., vol. 7, pp. 9. SIAM, Philaelphia, PA (3) 4. Discacciati, M., Quarteroni, A.: Convergence analysis of a subomain iterative metho for the finite element approximation of the coupling of Stoes an Darcy equations. Comput. Vis. Sci. 6(-3), 93 3 (4). Givoli, D., Barbone, P.E., Patlasheno, I.: Which are the important moes of a subsystem? Internat. J. Numer. Methos Engrg. 9(), (4) 6. Glover, K.: All optimal Hanel-norm approximations of linear multivariable systems an their L -error bouns. Internat. J. Control 39(6), 93 (984) 7. Glowinsi, R.: Finite element methos for incompressible viscous flow. Hanboo of Numerical Analysis, Vol. IX. North-Hollan, Amsteram (3) 8. Gugercin, S., Sorensen, D.C., Antoulas, A.C.: A moifie low-ran Smith metho for largescale Lyapunov equations. Numer. Algorithms 3(), 7 (3) 9. Hahn, J., Egar, T.F.: An improve metho for nonlinear moel reuction using balancing of empirical Gramians. Computers an Chemical Engineering 6, (). Hesthaven, J.S., Warburton, T.: Noal Discontinuous Galerin Methos: Analysis, Algorithms, an Applications. Springer-Verlag, Berlin (8). Hinze, M., Volwein, S.: Proper orthogonal ecomposition surrogate moels for nonlinear ynamical systems: Error estimates an suboptimal control. In: P. Benner, V. Mehrmann, D.C. Sorensen (es.) Dimension Reuction of Large-Scale Systems, Lecture Notes in Computational Science an Engineering, Vol. 4, pp Springer-Verlag, Heielberg (). Karniaais, G.E., Sherwin, S.J.: Spectral/hp Element Methos for Computational Flui Dynamics, secon en. Oxfor University Press () 3. Kunisch, K., Volwein, S.: Galerin proper orthogonal ecomposition methos for a general equation in flui ynamics. SIAM J. Numer. Anal. 4(), 49 (electronic) () 4. Lall, S., Marsen, J.E., Glavaši, S.: A subspace approach to balance truncation for moel reuction of nonlinear control systems. Int. J. Robust Nonlinear Control (6), 9 3 (). Lucia, D.J., King, P.I., Beran, P.S.: Reuce orer moeling of a two-imensional flow with moving shocs. Computer an Fluis 3, (3) 6. Moore, B.C.: Principal component analysis in linear systems: controllability, observability, an moel reuction. IEEE Trans. Automat. Control 6(), 7 3 (98) 7. Penzl, T.: Algorithms for moel reuction of large ynamical systems. Linear Algebra Appl. 4(-3), (6) 8. Segalman, D.J.: Moel reuction of systems with localize nonlinearities. Journal of Computational an Nonlinear Dynamics (3), (7) 9. Ture, S., Rivin, L., Hron, J., Glowinsi, R.: Numerical stuy of a moifie time-stepping θ-scheme for incompressible flow simulations. J. Sci. Comput. 8(-3), (6). Vanenorpe, A., Dooren, P.V.: On moel reuction of interconnecte systems. In: W. Schilers, H. van er Vorst (es.) Moel Orer Reuction: Theory, Research Aspects an Applications, Mathematics in Inustry, pp Springer Verlag, Heielberg (7)
Math 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
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 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 informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
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 informationApproximate Reduction of Dynamical Systems
Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,
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 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 informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
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 informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
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 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 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 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 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 informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
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
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 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 informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationAdaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements
Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan
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 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 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 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 informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
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 informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
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 informationOn Decentralized Optimal Control and Information Structures
2008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 2008 FrC053 On Decentralize Optimal Control an Information Structures Naer Motee 1, Ali Jababaie 1 an Bassam
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 informationWUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationDesign of an Industrial Distributed Controller Near Spatial Domain Boundaries
Design of an Inustrial Distribute Controller Near Spatial Domain Bounaries Stevo Mijanovic, Greg E. Stewart, Guy A. Dumont, an Michael S. Davies ABSTRACT This paper proposes moifications to an inustrial
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 informationBalancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling
Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
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 informationMinimum-time constrained velocity planning
7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università
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 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 informationGeneralized Nonhomogeneous Abstract Degenerate Cauchy Problem
Applie Mathematical Sciences, Vol. 7, 213, no. 49, 2441-2453 HIKARI Lt, www.m-hikari.com Generalize Nonhomogeneous Abstract Degenerate Cauchy Problem Susilo Hariyanto Department of Mathematics Gajah Maa
More informationARBITRARY NUMBER OF LIMIT CYCLES FOR PLANAR DISCONTINUOUS PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO ZONES
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 228, pp. 1 12. ISSN: 1072-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu ARBITRARY NUMBER OF
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationExponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity
Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationAnalyzing Tensor Power Method Dynamics in Overcomplete Regime
Journal of Machine Learning Research 18 (2017) 1-40 Submitte 9/15; Revise 11/16; Publishe 4/17 Analyzing Tensor Power Metho Dynamics in Overcomplete Regime Animashree Ananumar Department of Electrical
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 information1. Introduction. In this paper we consider the nonlinear Hamilton Jacobi equations. u (x, 0) = ϕ (x) inω, ϕ L (Ω),
SIAM J. UMER. AAL. Vol. 38, o. 5, pp. 1439 1453 c 000 Society for Inustrial an Applie Mathematics SPECTRAL VISCOSITY APPROXIMATIOS TO HAMILTO JACOBI SOLUTIOS OLGA LEPSKY Abstract. The spectral viscosity
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationNATURAL BOUNDARY ELEMENT METHOD FOR THREE DIMENSIONAL EXTERIOR HARMONIC PROBLEM WITH AN INNER PROLATE SPHEROID BOUNDARY
NATURAL BOUNDARY ELEMENT METHOD FOR THREE DIMENSIONAL EXTERIOR HARMONIC PROBLEM WITH AN INNER PROLATE SPHEROID BOUNDARY Hong-ying Huang De-hao Yu State Key Laboratory of Scientific Engineering Computing,
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 informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationTMA4195 Mathematical modelling Autumn 2012
Norwegian University of Science an Technology Department of Mathematical Sciences TMA495 Mathematical moelling Autumn 202 Solutions to exam December, 202 Dimensional matrix A: τ µ u r m - - s 0-2 - - 0
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
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 informationFrom Local to Global Control
Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department,
More informationOn Characterizing the Delay-Performance of Wireless Scheduling Algorithms
On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,
More informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
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 informationInterconnected Systems of Fliess Operators
Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear
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 informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
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 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 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 informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationCriteria for Global Stability of Coupled Systems with Application to Robust Output Feedback Design for Active Surge Control
Criteria for Global Stability of Couple Systems with Application to Robust Output Feeback Design for Active Surge Control Shiriaev, Anton; Johansson, Rolf; Robertsson, Aners; Freiovich, Leoni 9 Link to
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 informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationGlobal Solutions to the Coupled Chemotaxis-Fluid Equations
Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria
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 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 informationA transmission problem for the Timoshenko system
Volume 6, N., pp. 5 34, 7 Copyright 7 SBMAC ISSN -85 www.scielo.br/cam A transmission problem for the Timoshenko system C.A. RAPOSO, W.D. BASTOS an M.L. SANTOS 3 Department of Mathematics, UFSJ, Praça
More informationPerturbation Analysis and Optimization of Stochastic Flow Networks
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
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 informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationBalance laws on domains with moving interfaces. The enthalpy method for the ice melting problem.
Balance laws on omains with moving interfaces. The enthalpy metho for the ice melting problem. Gowin Kakuba 12th March 2008 Outline 1 Balance equations on omains with moving bounaries Introuction Rankine-Hugoniot
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
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 informationA new identification method of the supply hole discharge coefficient of gas bearings
Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico
More information