Numerical discretization of tangent vectors of hyperbolic conservation laws.
|
|
- Gervase Hensley
- 5 years ago
- Views:
Transcription
1 Numerical discretization of tangent vectors of hyperbolic conservation laws. Michael Herty IGPM, RWTH Aachen joint work with Benedetto Piccoli MNCFF 2014, Bejing,
2 Contents Originally on multi-phase fluid flow on networks (extension of coupling conditions - NumHyp 2013 ) Recent results on tangent vectors in collaboration with B. Piccoli (Rutgers) Tangent vectors are used for optimization problems where the dynamics is governed by 1 d hyperbolic systems Tangent vectors are a first order sensitivity calculus for solutions u 0 u(t, x) Examples are control of compressor stations in gas networks, control of gates in open canals, parameter identification problems subject to transport problems,...
3 This presentation References and theoretical problem Example of tangent vectors for Burgers Numerical challenges Approximation of evolution of tangent vectors Numerical examples
4 References Concept of tangent vectors present one possibility to compute and analyze sensitivity of systems on 1-d hyperbolic balance laws Existing calculus for spatially one dimensional systems of conservation laws Introduced by Bressan/Marson (1995), extended by Bressan and co workers (1997,2007), Bianchini (2000), Piccoli and co workers (2000 ) Lipschitz continuous dependence for 1 d system Crasta/Bressan/Piccoli (2001) Related for scalar convex case: Ulbrich (2001), Giles (1996), Ulbrich/Giles (2010), Zuazua/Castro et al (2008 ). Application of sensitivity equations for general inverse problems, optimal control problems, and control / controllability questions,...
5 Theoretical challenge t u + x f (u) = 0, u(t = 0, x) = u 0 (x) Solution operator for nonlinear conservation laws u(t, ) = S t u 0 generically non differentiable on L 1 No classical calculus for first order variations of S t u 0 with respect to u 0 No characterizing conditions for parameter identification problems or optimal control problems min J(u(T, x))dx subject to t u+ x f (u) = 0, u(t = 0, x) = u 0 (x) u 0
6 Example for Burger s equation t u + x u 2 2 = 0, u(t = 0, x) = uɛ 0 = ɛ x χ [0.1] variations of initial data and interest in the behavior of the solution S t u0 ɛ with respect to h at e.g. ɛ u ɛ 0 S t u ɛ+h 0 S t u ɛ 0 + hv + O(h) Exact solution u ɛ (t, x) = (1+ɛ)x wrt ɛ Shock position is 1 + ɛt 1+ɛt χ [0, 1+ɛt] is Lipschitz in L1 t = 0: a first order approximation in L 1 exists and is v(t, x) = lim ɛ+h u0 (x) u0 ɛ(x) h 0 ɛ = ɛ x χ [0,1] But for any t > 0 the jump depends on ɛ + h and the previous limit does not define any function in L 1
7 Idea of tangent vectors for Burger s equation lim h 0 u ɛ+h (t, x) u ɛ (t, x) dx = O(1)+ 1 1+(ɛ+h)t h h = O( 1 h ) 1+ɛt Use weak formulation to compute the limit h 0 for arbitrary value of ɛ lim h 0 lim h 0 lim h 0 u ɛ+h (t, x) u ɛ (t, x) φ(x)dx = h x (1 + (ɛ + h)t)(1 + ɛt) χ [0, 1+ɛt] (x)φ(x)dx+ (ɛ + h)x 1 + (ɛ + h)t dx t (ɛ + h)t φ(y) (ɛ + h)y 1 + (ɛ + h)t, y [ 1 + ɛt, 1 + (ɛ + h)t]
8 Idea of tangent vectors (cont d) lim h 0 lim h 0 x (1 + (ɛ + h)t)(1 + ɛt) χ [0, 1+ɛt] (x)φ(x)dx+ t (ɛ + h)t φ(y) (ɛ + h)y 1 + (ɛ + h)t, y [ 1 + ɛt, 1 + (ɛ + h)t] A suitable differential of u ɛ (t, x) at ɛ may therefore consist of two components: an L 1 part and a measure located at the jump of u ɛ S t u ɛ+h 0 = u ɛ+h (t, ) = S t u ɛ 0 + h ( v(t, ) + δ s ɛ (t)( )[u ɛ ]ξ(t) )
9 Idea of tangent vectors for Burger s equation (2/2) lim h 0 lim h 0 x (1 + (ɛ + h)t)(1 + ɛt) χ [0, 1+ɛt] (x)φ(x)dx+ t (ɛ + h)t φ(y) (ɛ + h)y 1 + (ɛ + h)t, y [ 1 + ɛt, 1 + (ɛ + h)t] v(t, x) is the L 1 part consists of the a.e. pointwise limit u ɛ+h (t, x) u ɛ (t, x) v(t, x) = lim = h 0 h x (1 + ɛt) 2 ξ(t) is a real number and is the variation in the shift of the shock position s ɛ s ɛ+h (t) s ɛ (t) t ξ(t) = lim = h 0 h ɛt, u ɛ+h (t, ) = u ɛ (t, ) + h ( v(t, ) + δ s ɛ (t)( )[u ɛ ]ξ(t) ) (v, ξ) is called generalized tangent vector
10 Generalized tangent vectors (v, ξ) v(t, x) is the L 1 part consists of the a.e. pointwise limit u v(t, x) = lim ɛ+h (t,x) u ɛ (t,x) h 0 h ξ(t) is a real number and is the variation in the shift of the shock position s ɛ s ξ(t) = lim ɛ+h (t) s ɛ (t) h 0 h u ɛ+h (t, ) = u ɛ (t, ) + h ( v(t, ) + δ s ɛ (t)( )[u ɛ ]ξ(t) ) Question Which variations ɛ u0 ɛ lead to well defined tangent vectors? Different settings possible. Here, we follow B/M with u0 ɛ piecewise Lipschitz with a finite number of isolated discontinuities A tangent vector a t = 0 may lead to a tangent vector at time t by solving suitable evolution equations (first order variations of hyperbolic system). For hyperbolic balance laws and suitable variations this is proven up to shock interaction (B/M) and extended in case of conservation laws (C/B/P).
11 Suitable variations of initial data leading to well defined tangent vectors (v, ξ) A suitable variation u h of u 0 (x) is within an L 1 equivalence class to u h (x) = u 0 (x)+hv(x)+ ξ i <0 [u 0 ]χ [xi +hξ,x i ](x) ξ i >0[u 0 ]χ [xi,x i +hξ](x) for some v L 1 and ξ R N where N is the number of isolated discontinuities located at x i in u 0 (x).
12 Theoretical results (Theorem 2.2 (B/M) For previous variations the evolution equation for the tangent vector is well defined up to any point in time when two shocks coincide Result holds for system of balance laws u t + f (u) x = h(t, x, u), A(u) = Df (u).
13 Numerical implementation of tangent vectors Requires knowledge on exact shock positions for evaluation of the evolution of ξ(t) Requires solution of compatibility condition between waves of different families Evolution of the L 1 part v influences evolution of shift ξ(t) Non conservative system for v
14 Approximation to the problem for numerical implementation Presentation on the simplest case of a 1 d scalar conservation law y (1) t + f (y (1) ) x = 0, y (1) (t = 0, x) = u 0 (x) Replace equation by Jin Xin relaxation approximation for 0 < µ << 1 ( ) 0 1 y t + a 2 y 0 x = 1 ( ) 0 µ f (y (1) ) y (2) y (1) (0, x) = u 0 (x), y (2) (0, x) = f (u 0 (x)) Relaxation is a viscous approximation (( y (1) t + f (y (1) ) x = µ a 2 f (y (1) ) 2) ) y x (1) System is linear hyperbolic with stiff source term x.
15 Tangent vectors for relaxation system y t + ( ) 0 1 a 2 y 0 x = 1 ( ) 0 µ f (y (1) ) y (2) L 1 part of tangent vector does not depend on y (1) and is linear system ( ) 0 1 v(0, ) = v( ), v t + a 2 v 0 x = 1 ( ) 0 µ f (y (1) )v (1) v (2), v = (v (1), v (2) ) outside the discontinuities of y N discontinuities and l j jth eigenvector and i v = v(x i (t)+, t) v(x i (t), t) jump across ithe discontinuities ξ i (t) = ξ i and l j ( i v + i y x ξ i ) = 0 j k.
16 Variation of general cost functional wrt initial data Given y d L 1 (R), T > 0 given and y = y (1) solution to relaxation system J(y(, ), y d ( )) = χ I (x) (y(t, x) y d (x)) 2 dx Variation of J with respect to variations in u 0 generating the tangent vectors (v, ξ) and a fixed number N of discontinuities u0 J(y, y d ) = 2 ( ) χ I (x) y (1) (T, x) y d (x) v (1) (T, x)dx+ ( ( ) y (1) (T, x i +) y d (x i +) + N(u 0 ) i=1 ( )) y (1) (T, x i ) y d (x i ) i y (1) (T, )ξ i (T )
17 Variation of cost and tangent vectors used for maximal descent ( u0 J(y, y d ) = 2 χ I (x) y (1) ) (T, x) y d (x) v (1) (T, x)dx+ N(u 0 ) ( ( y (1) ) ( (T, x i +) y d (x i +) + y (1) )) (T, x i ) y d (x i ) i y (1) (T, )ξ i (T ) i=1 Update of control u 0 using gradient based information with stepsize ρ > 0 ũ 0 (x) = u 0 (x) ρv(0, x) N(u 0 ) i=1 χ [xi +ρξ i (0),x i+1 +ρξ i+1 (0)](x)u 0 (x i +) Choice of v for maximal descent in J ( ) v (1) (T, x) = y (1) (T, x) y d (x), v (2) (T, x) = 0 Choice of ξ for maximal descent in J ξ i (T ) = (( y (1) ) (T, x i +) y d (x i +) + ( y (1) )) (T, x i ) y d (x i ) i y (1) (T )
18 Numerical implementation of tangent vectors i=1 ( ) 0 1 v t + a 2 v 0 x = 1 ( µ ( u0 J(y, y d ) = 2 χ I (x) y (1) ) (T, x) y d (x) v (1) (T, x)dx+ N(u 0 ) ( ( y (1) ) ( (T, x i +) y d (x i +) + y (1) )) (T, x i ) y d (x i ) i y (1) (T, )ξ i (T ) 0 f (y (1) )v (1) v (2) ), ξ i (t) = ξ i and l j ( i v + i y x ξ i ) = 0 j k. Introduce a spatial grid (x i ) n i=1 and intervals I i = [x i 1, x 2 i+ 1 ] 2 Choose a control u 0 a piecewise constant on I i leading to ξ R 2n possible shock variations ( ) 0 1 Diagonalise matrix a 2 and derive tangent vectors for 0 system in diagonal form Splitting of transport and source term integration
19 Advantage of relaxation system over nonlinear system v t + ( ) 0 1 a 2 v 0 x = 1 ( µ 0 f (y (1) )v (1) v (2) ), ξ i (t) = ξ i and l j ( i v + i y x ξ i ) = 0 j k. Due to the linearity of system the equations for L 1 and real part ξ decouple Nonlinear source term does not influence evolution of t ξ(t) Variations of initial data on each grid cell I i leading to shock variations ξ i at each grid cell and an L 1 variation piecewise constant v i on each grid cell Evaluation of source term may lead to additional discontinuities and therefore additional shock variations Apply a splitting technique and use piecewise constant control in each component (v (1) 0, v (2) 0 ) spanning two consecutive intervals
20 Diagonalise the system in variables η with source term S and discretization according to the splitting above using a transformation T Choose control η 0 in first component constant on even and in the second on odd gridpoints Choose x and t according to the CFL condition with constant one Exact evaluation of first order Upwind scheme for the transport part Admissible controls and discretization of tangent vectors (v, ξ) v t + ( ) 0 1 a 2 v 0 x = 1 ( µ 0 f (y (1) )v (1) v (2) ), t η (1) + a x η (1) = 0, t η (2) = 0, t η (1) = S(η (1), η (2) ), t η (2) = S(η (1), η (2) ), t η (2) a x η (2) = 0, t η (1) = 0, (η (1) 0 ) 2i = (η (1) 0 ) 2i+1, (η (2) 0 ) 2i 1 = (η (2) 0 ) 2i, i = 0,..., n 2
21 Admissible controls and discretization of tangent vectors (v, ξ) (cont d) t η (1) + a x η (1) = 0, t η (2) = 0, t η (1) = S(η (1), η (2) ), t η (2) = S(η (1), η (2) ), t η (2) a x η (2) = 0, t η (1) = 0, x 2i 1 (t n ) = x 2i 1 (0) + a t n, x 2i (t n ) = x 2i (0) a t n, i = 0,..., n, Splitting does not need introduce additional shocks since exact Upwind leads to piecewise constant functions on even and odd grid cells Shock position at time t n is given x i (t n ) = xi 0 ± at n depending if discontinuity at even or odd cell boundary Condition l j ( i v + i y x ξ i ) = 0 j k, automatically satisfied since there is only one discontinuity at each cell boundary at even grid points in the first and odd in the second component
22 Tests and algorithm for control problems with cost J min y 0 =(y (1) 0,y (2) 0 ) J sbj to relaxation system Unregularized cost J(y(, ), y d ( )) = χ I (x) (y(t, x) y d (x)) 2 dx Update of control u 0 using gradient based information with stepsize ρ > 0 ũ 0 (x) = u 0 (x) ρv(0, x) N(u 0) i=1 Choice of v for maximal descent in J v (1) (T, x) = Choice of ξ for maximal descent in J ξ i (T ) = χ [xi +ρξ i (0),x i+1 +ρξ i+1 (0)](x)u 0 (x i +) ( ) y (1) (T, x) y d (x), v (2) (T, x) = 0 (( ) ( )) y (1) (T, x i +) y d (x i +) + y (1) (T, x i ) y d (x i ) i y (1) (T
23 Algorithm for computing a gradient 1. Set terminal time T > 0, a 2 max(f (y (1) ) 2 and choose an y (1) equidistant spatial discretization( with N x gridpoints. ) Choose t = x a and k = 0. Let η0,i k = (η (1) 0,i )k, (η (2) 0,i )k for i = 0,..., N x, be an arbitrary initial control such that η0 k is admissible Let (y d ) i be a discretization of the given function y d ( ). 2. Solve forward relaxation equations with (y 0 ) i := T η0,i k to obtain η Nt i = T 1 y Nt i. 3. Set initial data ϕ 0 i and shock variations ξ i where x i (t Nt ) are given by equation shock positions 4. Solve sensitivity equations backwards in time for initial data ṽi 0 := T ϕ 0 i to obtain ϕ Nt i = T 1 ṽ Nt i. 5. Update the current iterate η0,i k by setting η k+1 0,i := η i 0 by ηi 0 := η0,i k. 6. Provided that J(T η, y d ) is sufficiently small we terminate. Otherwise set k k + 1 and continue with step (2).
24 Numerical results: Grid convergence for splitting
25 Numerical results: Convergence for a linear system without source term
26 Numerical results: Convergence table using (v, ξ) components
27 Numerical results: Convergence table using only v components
28 Numerical results: Relaxation for Burger s equation
29 Numerical results: Convergence
30 Thank you for your attention.
NUMERICAL METHOD FOR THE COMPUTATION OF TANGENT VECTORS TO 2 2 HYPERBOLIC SYSTEMS OF CONSERVATION LAWS
NUMERICAL METHOD FOR THE COMPUTATION OF TANGENT VECTORS TO HYPERBOLIC SYSTEMS OF CONSERVATION LAWS MICHAEL HERTY AND BENEDETTO PICCOLI Abstract. We are interested in the development of a numerical method
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationA model for a network of conveyor belts with discontinuous speed and capacity
A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica - 6.03.2018 work in collaboration with M. Pfirsching, S. Goettlich
More informationAn optimal control problem for a parabolic PDE with control constraints
An optimal control problem for a parabolic PDE with control constraints PhD Summer School on Reduced Basis Methods, Ulm Martin Gubisch University of Konstanz October 7 Martin Gubisch (University of Konstanz)
More informationM. HERTY, CH. JÖRRES, AND B. PICCOLI
EXISTENCE OF SOLUTION TO SUPPLY CHAIN MODELS BASED ON PARTIAL DIFFERENTIAL EQUATION WITH DISCONTINUOUS FLUX FUNCTION M. HERTY, CH. JÖRRES, AND B. PICCOLI Abstract. We consider a recently [2] proposed model
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationNumerical Methods for the Optimal Control of Scalar Conservation Laws
Numerical Methods for the Optimal Control of Scalar Conservation Laws Sonja Steffensen, Michael Herty, and Lorenzo Pareschi RWTH Aachen University, Templergraben 55, D-52065 Aachen, GERMANY {herty,steffensen}@mathc.rwth-aachen.de
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationOn a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws
On a Lyapunov Functional Relating Shortening Curves and Viscous Conservation Laws Stefano Bianchini and Alberto Bressan S.I.S.S.A., Via Beirut 4, Trieste 34014 Italy. E-mail addresses: bianchin@mis.mpg.de,
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 (Many figures from C. M. Bishop, "Pattern Recognition and ") 1of 254 Part V
More informationAPM526 CLASSROOM NOTES PART IB. APM526, Spring 2018 Last update: MAR 12
APM526 CLASSROOM NOTES PART IB APM526, Spring 2018 Last update: MAR 12 1 PARTICLE METHODS 1. Particle based methods: Large systems of ODEs for individual particles. 2. Partial differential equations: Equations
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More informationOn finite time BV blow-up for the p-system
On finite time BV blow-up for the p-system Alberto Bressan ( ), Geng Chen ( ), and Qingtian Zhang ( ) (*) Department of Mathematics, Penn State University, (**) Department of Mathematics, University of
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationConvex Optimization. Problem set 2. Due Monday April 26th
Convex Optimization Problem set 2 Due Monday April 26th 1 Gradient Decent without Line-search In this problem we will consider gradient descent with predetermined step sizes. That is, instead of determining
More informationRecitation 1. Gradients and Directional Derivatives. Brett Bernstein. CDS at NYU. January 21, 2018
Gradients and Directional Derivatives Brett Bernstein CDS at NYU January 21, 2018 Brett Bernstein (CDS at NYU) Recitation 1 January 21, 2018 1 / 23 Initial Question Intro Question Question We are given
More informationA maximum principle for optimally controlled systems of conservation laws
RENDICONTI del SEMINARIO MATEMATICO della UNIVERSITÀ DI PADOVA ALBERTO BRESSAN ANDREA MARSON A maximum principle for optimally controlled systems of conservation laws Rendiconti del Seminario Matematico
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationAN ALTERNATING DESCENT METHOD FOR THE OPTIMAL CONTROL OF THE INVISCID BURGERS EQUATION IN THE PRESENCE OF SHOCKS
December 7, 2007 10:10 WSPC/INSTRUCTION FILE Castro- Palacios- Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company AN ALTERNATING DESCENT METHOD FOR THE OPTIMAL CONTROL
More informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
More informationDownloaded 04/15/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 51, No. 1, pp. 126 136 c 2013 Society for Industrial and Applied Mathematics A NOTE ON ADJOINT ERROR ESTIMATION FOR ONE-DIMENSIONAL STATIONARY BALANCE LAWS WITH SHOCKS JOCHEN
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationA re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws
A re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c,3 a Department of Applied Mathematics and National
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
More informationInstability of Finite Difference Schemes for Hyperbolic Conservation Laws
Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Alberto Bressan ( ), Paolo Baiti ( ) and Helge Kristian Jenssen ( ) ( ) Department of Mathematics, Penn State University, University
More informationConservation laws and some applications to traffic flows
Conservation laws and some applications to traffic flows Khai T. Nguyen Department of Mathematics, Penn State University ktn2@psu.edu 46th Annual John H. Barrett Memorial Lectures May 16 18, 2016 Khai
More informationSOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS
SOME PRELIMINARY NOTES ON HYPERBOLIC CONSERVATION LAWS KAYYUNNAPARA T. JOSEPH CENTRE FOR APPLICABLE MATHEMATICS TATA INSTITUTE OF FUNDAMENTAL RESEARCH BANGALORE 5665 1. Introduction Many of the physical
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationAn alternating descent method for the optimal control of the inviscid Burgers equation in the presence of shocks
An alternating descent method for the optimal control of the inviscid Burgers equation in the presence of shocks Carlos CASTRO, Francisco PALACIOS and Enrique ZUAZUA June 20, 2007 Abstract We introduce
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationRegularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws
1 Regularity and Lyapunov stabilization of weak entropy solutions to scalar conservation laws Sébastien Blandin, Xavier Litrico, Maria Laura Delle Monache, Benedetto Piccoli and Alexandre Bayen Abstract
More informationMeasure-valued - strong uniqueness for hyperbolic systems
Measure-valued - strong uniqueness for hyperbolic systems joint work with Tomasz Debiec, Eduard Feireisl, Ondřej Kreml, Agnieszka Świerczewska-Gwiazda and Emil Wiedemann Institute of Mathematics Polish
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationIntersection Models and Nash Equilibria for Traffic Flow on Networks
Intersection Models and Nash Equilibria for Traffic Flow on Networks Alberto Bressan Department of Mathematics, Penn State University bressan@math.psu.edu (Los Angeles, November 2015) Alberto Bressan (Penn
More informationNetwork Newton. Aryan Mokhtari, Qing Ling and Alejandro Ribeiro. University of Pennsylvania, University of Science and Technology (China)
Network Newton Aryan Mokhtari, Qing Ling and Alejandro Ribeiro University of Pennsylvania, University of Science and Technology (China) aryanm@seas.upenn.edu, qingling@mail.ustc.edu.cn, aribeiro@seas.upenn.edu
More informationVariational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University
Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic
More informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationVariational formulation of entropy solutions for nonlinear conservation laws
Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for
More informationPrimal/Dual Decomposition Methods
Primal/Dual Decomposition Methods Daniel P. Palomar Hong Kong University of Science and Technology (HKUST) ELEC5470 - Convex Optimization Fall 2018-19, HKUST, Hong Kong Outline of Lecture Subgradients
More informationSecond Order ODEs. Second Order ODEs. In general second order ODEs contain terms involving y, dy But here only consider equations of the form
Second Order ODEs Second Order ODEs In general second order ODEs contain terms involving y, dy But here only consider equations of the form A d2 y dx 2 + B dy dx + Cy = 0 dx, d2 y dx 2 and F(x). where
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationPartial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations
Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationOptimal Control and Vanishing Viscosity for the Burgers Equation
7 Optimal Control and Vanishing Viscosity for the Burgers Equation C. Castro, F. Palacios, and E. Zuazua 3 Universidad Politécnica de Madrid, Spain; carlos.castro@upm.es Universidad Politécnica de Madrid,
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationOn Asymptotic Variational Wave Equations
On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationarxiv: v3 [math.ds] 22 Feb 2012
Stability of interconnected impulsive systems with and without time-delays using Lyapunov methods arxiv:1011.2865v3 [math.ds] 22 Feb 2012 Sergey Dashkovskiy a, Michael Kosmykov b, Andrii Mironchenko b,
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 1: Introduction to Multigrid 2000 Eric de Sturler 1 12/02/09 MG01.prz Basic Iterative Methods (1) Nonlinear equation: f(x) = 0 Rewrite as x = F(x), and iterate x i+1
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationarxiv: v2 [math.na] 24 Mar 2016
arxiv:1504.04107v2 [math.na] 24 Mar 2016 Strong stability preserving explicit linear multistep methods with variable step size Yiannis Hadjimichael David I. Ketcheson Lajos Lóczi Adrián Németh March 21,
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationProximal-like contraction methods for monotone variational inequalities in a unified framework
Proximal-like contraction methods for monotone variational inequalities in a unified framework Bingsheng He 1 Li-Zhi Liao 2 Xiang Wang Department of Mathematics, Nanjing University, Nanjing, 210093, China
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More information, b = 0. (2) 1 2 The eigenvectors of A corresponding to the eigenvalues λ 1 = 1, λ 2 = 3 are
Quadratic forms We consider the quadratic function f : R 2 R defined by f(x) = 2 xt Ax b T x with x = (x, x 2 ) T, () where A R 2 2 is symmetric and b R 2. We will see that, depending on the eigenvalues
More informationApproximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle
Malaya J. Mat. 4(1)(216) 8-18 Approximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle B. C. Dhage a,, S. B. Dhage a and S. K. Ntouyas b,c, a Kasubai,
More informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
More informationNumerisches Rechnen. (für Informatiker) M. Grepl P. Esser & G. Welper & L. Zhang. Institut für Geometrie und Praktische Mathematik RWTH Aachen
Numerisches Rechnen (für Informatiker) M. Grepl P. Esser & G. Welper & L. Zhang Institut für Geometrie und Praktische Mathematik RWTH Aachen Wintersemester 2011/12 IGPM, RWTH Aachen Numerisches Rechnen
More information