A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
|
|
- Job Weaver
- 5 years ago
- Views:
Transcription
1 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 J.A Mackenzie and A. Nicola Department of Mathematics, University of Strathclyde, 6 Richmond St, Glasgow. jam@maths.strath.ac.uk In this paper we consider the adaptive numerical solution of Hamilton-Jacobi HJ equations φ t + Hφ,...,φ d =, φ, = φ, where =,..., d IR d, t >. HJ equations arise in many practical areas such as differential games, mathematical finance, image enhancement and front propagation. It is well known that solutions of are Lipschitz continuous but derivatives can become discontinuous even if the initial data is smooth. There is a close relation between HJ equations and hyperbolic conservation laws. With this in mind, it not surprising to find that many of the numerical methods used to solve HJ equations are motivated by conservative finite difference or finite volume methods for conservation laws. An increasingly popular approach to solve hyperbolic conservation laws is the discontinuous Galerkin DG finite element method. Recently, Hu and Shu [] proposed a DG method to solve HJ equations by first rewriting as a system of conservation laws w i t + Hw i =, i =,...,d, w, = φ, where w = φ. The usual DG formulation would be obtained if w belonged to a space of piecewise polynomials. However, we note that w i, i =,...,d are not independent due to the restriction that w = φ. In [] a least squares procedure was used to enforce this condition whereas in [4] a suitable basis was used such that w h = φ h The aim of this paper is to consider the use of the DG method of Hu and Shu [] to solve HJ equation using a moving mesh method based on the solution of moving mesh PDEs []. The governing equation are transformed to include the effect of the movement of the mesh and this is done in such a way that the conservation properties of the original equation are not lost. The adaptive mesh is driven by a monitor function which is shown to be non-singular in the presence of solution discontinuities.. DG discretisation of Hamilton-Jacobi equations... HJ equation in D Let us consider the Hamilton-Jacobi equation in one dimension { φt + Hφ =,,t a,b,t] φ, = φ, with appropriate boundary conditions. If u = φ and we differentiate 3 with respect to then u satisfies the hyperbolic conservation law { ut + Hu =,,t a,b,t] 4 u, = u. 3
2 We assume that the physical domain Ω p = [a,b] is the image of a computational domain Ω c = [,] which is obtained via the time-dependent mapping = ξ,t. We will consider a discretisation of the conservative form in computational coordinates { ξ u + Hu ẋu ξ =, ξ,t,,t], 5 uξ, = u ξ. The aim is to discretise 5 in space using a DG method.... The DG discretisation { } N For each partition ξ j+ of Ω c, we denote j= I j = ξ j,ξ j+, j = ξ j+ ξ j, j =,...,N. We use a uniform mesh to cover the computational domain such that h = ξ j+ ξ j = N, ξ j = ξ j+ + ξ j, j =,...,N. 6 If we multiply 5 by an arbitrary smooth function v, then integrating by parts over the interval I j we obtain ξ uv dξ Hu,ẋvξ dξ + Hu,ẋv =, 7 I j I Ij j where and Hu, ẋv = Huξ j+ Ij Hu,ẋ = Hu ẋu, 8,t,ẋξ j+ vξ j+ Huξ j,t,ẋξ j vξ j. 9 We will approimate the eact solution u of 7 and the smooth function v locally by polynomials of degree at most k. The numerical approimations u h and v h belong to the finite dimensional space } Vh {v k = : v Ij P k I j, j =,...,N, where P k I j denotes the space of polynomials on I j of degree at most k. In the discontinuous Galerkin numerical method u h and v h are discontinuous at the points ξ j+, j =,...,N, and the question is how to evaluate their values in 9. We set v h j+ = lim v h ξ, and v + ξ ξ h = lim v j h ξ. j+ ξ ξ + j Furthermore, we replace the nonlinear flu function Hu h,ẋ defined in 8 by a numerical flu function that depends on the values of u h from the left and right at the point,t, that is Hu h,ẋ ξj+ t = H j+ In this paper we have used the local La-Friedrichs flu Hp,q,ẋ = ξ j+ u h ξ,t,u j+ h ξ +,t,ẋ j+ ξj+ t. Hp,ẋ + Hq,ẋ cq p, 3
3 where c = ma minp,q u map,q H u. The discontinuous Galerkin space approimation is given as the solution of the weak formulation ξ u h v h dξ Huh,ẋv h ξ dξ + H j+ v h ξj+ H j v + h ξj =, 4 I j I j and the initial condition is given by the projection u h ξ,v h ξ dξ = u ξv h ξ dξ, j =,...,N, 5 I j I j for all v h V h.... Temporal integration The discontinuous Galerkin discretisation gives rise to a system of ODEs which can be written as { d dt ξu h = L h u h, in,t] 6 u h = u h, where u h is the initial approimation for the vector of the degrees of freedom u h, and ξ depends on the meshes at time level t n and t n+ as defined earlier. To integrate 6 numerically we use the second-order TVD Runge-Kutta method n+ ξ u h = n ξ un h + t Lun h, n+ ξ u n+ h = n ξ un h + n+ ξ u h + 7 t Lu h.... Recovery of φ h If φ h = u h, then φ h will be determined on each interval I j up to a constant. Following [], the constant can be retrieved in two ways:. require that I j φ h t + Hu h d, j =,...,N. 8. use 8 to update only the leftmost element I and then since φ h = u h we have..3. HJ equations in D In two-dimensions we have φ h j,t = φ h,t + j u h,t d. 9 φ t + Hφ,φ y =,,y Ω p, t >, where Ω p lr is the physical domain. To generate an adaptive mesh we have seen earlier that it is useful to regard the physical domain Ω p as the image of a computational domain Ω c under the invertible maps = ξ,η,t, y = yξ,η,t, and ξ = ξ,y,t, η = η,y,t, 3
4 where =, y and ξ = ξ, η are the physical and computational coordinates, respectively. A mesh covering Ω p is obtained by applying the mapping given by to a partitioning of Ω c. Given that the mesh can move then it is necessary to epress the Eulerian -fied temporal derivative in in terms of the Lagrangian derivative along the trajectory of the moving mesh. If a dot denotes differentiation with respect to time with ξ,η fied then becomes φ ẋφ ẏφ y + Hφ,φ y =, where ẋ,ẏ represents the mesh velocity. If we let u = φ ξ and v = φ η, then by differentiating with respect to ξ and η we arrive at the system u t + H ξ u,v,ẋ,ẏ =, 3 where v t + H η u,v,ẋ,ẏ =, 4 Hu,v,ẋ,ẏ = Hξ u + η v,ξ y u + η y v ẋξ u + η v ẏξ y u + η y v. 5 Using the vectorial notation u = u v H, F =, F = H we can rewrite 3 and 4 as the system of conservation laws u t + [F u] ξ + [F u] η =. 6 We apply a discontinuous Galerkin method to solve this coupled system 6 for u,v = ξ φ, where the discretisation takes place in the computational domain Ω c. A recovery procedure will then be used to obtain an approimation of φ. To obtain a weak formulation we take the inner product of 6 with a test function v h V h, integrate over K T h and replace u by its approimation u h V h. That is d dt K u h ξ,t v h ξ dξ + and using integration by parts we obtain d u h ξ,t v h ξ dξ dt K + e K e K, [F ξ v h + F η v h ] dξ =, 7 { } n e,k F v h ξ intk + n e,k F v h ξ intk ds [F v h ξ + F v h η ] dξ =, K where n e,k = [n e,k,n e,k ] is the unit normal to the edge e of the element K, intk denotes the value taken from the within the element K, and etk denotes the value taken from the eterior of the element K. Net we define the normal flu F e,k u,ẋ,ẏ = n e,k F + n e,k F = [ ne,k H n e,k H 8 ]. 9 Note that F e,k is not well defined on the element edge as u h is discontinuous there. Therefore, we replace it by the numerical flu function F e,k u h ξ intk,u h ξ etk,ẋ,ẏ. 3 4
5 Again we will use the La-Friedrichs flu F e,k a,b,ẋ,ẏ = [F e,ka + F e,k b + α e,k a b], 3 where F α e,k = γρ n e,k u + n F e,k, 3 u and γ and ρ is the spectral radius of the normal Jacobian matri. We now discuss the basis used when u h is a discontinuous bilinear function. Following Hu and Shu [] we assume that the approimate solution of the Hamilton-Jacobi equation φ h is piecewise discontinuous and quadratic. Therefore, in the cell K ij we assume where φ h t = φt + φ ξ tµ i + φ η tν j + φ ξη tµ i ν j +φ ξξ t µ i + φ ηη t νj 3, 33 3 µ i ξ = ξ ξ i ξ, and ν j η = η η j. η The gradient in computational space ξ φ h therefore takes the form [ ξ ξ φ h = φ ξt + ξ φ ξηtν j + 4 ξ φ ] ξξtµ i η φ ηt + η φ ξηtµ i + 4 η φ, 34 ηηtν j and hence there are five unknowns that determine ξ φ h. A suitable basis for V = {v,v,v 3,v 4,v 5 : v K = ξ φ, φ P K, K T h } is given by b = [ ] [, b = ] [ µi, b 3 = ] [, b 4 = ν j ] [ νj, b 5 = ξ µ i η ]. 35 The local basis functions 35 can easily be shown to be orthogonal and hence the mass matri is diagonal. 3. Moving mesh equations Let us consider =, T to be a point in the physical domain Ω p and ξ = ξ,ξ T to be a point in the computational domain Ω c. The mapping from the computational domain to the physical domain is assumed to satisfy the moving mesh PDE [] τ t = a i,i + b i i= a i G a j ξ i ξ j i,j= a i G a j, 36 ξ j ξ i i,j= where a i = ξ i and a i,j = a i G a j, b i = i,j= a i G a j. 37 ξ j 5
6 Dirichlet boundary conditions for the above system are obtained by solving a one-dimensional MMPDE. If p Ω p and c Ω c denote the physical and computational boundary segments with arc-lengths l and l c respectively, then the mesh on p is the solution of τ s t = M + M ζ ζ M s ζ, ζ,l c, 38 with s = and sl c = l. Here M is the one-dimensional projection of the two-dimensional monitor function along the boundary. That is, if t is a unit tangent vector along the boundary then Ms,t = t T Gt Choice of monitor function The choice of a suitable adaptivity criterion for problems which develop shocks is highly nontrivial. It is important that the mesh points move smoothly towards regions where discontinuities eist so that the O error at the shock is localised within one or two mesh elements. For robustness, it is also important that the rate of convergence of the method is not severely impaired by the non-uniformity of the moving mesh. To obtain a bounded monitor function we consider the choice [3] M,t = + u, α = α β ma u, 39 and β is a user-chosen parameter. It is clear that this monitor function satisfies the bounds M,t + β. If the mesh is found by equidistribution of the monitor function then if M ecept in the cell containing then + β min b a + + β N min. 4 Therefore, the minimum mesh spacing is min b a N + β 4 which is approimately a factor of + β smaller than that using a uniform mesh. For two dimensional problems the monitor function used is of Winslow type where [ ] w G = w and w = + u α, where α = β ma Ωp u. This monitor function behaves in a similar fashion to the one-dimensional monitor in the presence of solution discontinuities. 4. Numerical eperiments 4... Eample The first eample considered is Burgers equation φ t + φ + =, < <, t > φ, = cosπ.85, with periodic boundary conditions. 4 6
7 u h u h φ h.4 φ h Figure : DG numerical solution and eact solution + obtained on a uniform stationary mesh left, and on a moving mesh right at time t = 3/π Eample We consider the two-dimensional Burgers equation with periodic boundary conditions. 5. Conclusions φ t + + φ + φ y =,,y, 43 φ,y, = cosπ + y/, 44 We have presented a DG method for the solution of nonlinear Hamilton-Jacobi equations using an adaptive moving mesh refinement strategy. A suitable mesh refinement criterion was used which does not become singular when the solution become non-smooth. Numerical eperiments in one and two dimensions show that method works well and efficiently delivers high resolution solutions using relatively coarse meshes. References [] C. Hu and C.W. Shu. A Discontinuous Galerkin Finite Element Method for Hamilton- Jacobi Equations. SIAM Journal on Scientific Computing, :666 69,
8 .5.5 φ.5 φ.5.5 y.5 y Figure : Approimations of φ top and φ bottom at t =.5/π using uniform mesh left and using moving mesh right, with N = y.5 y Figure 3: Adapted meshes at t =.5/π left, and t =.5/π right obtained with N = 4 for Eample. [] W. Huang. Practical aspects of formulation and solution of moving mesh partial differential equations. Journal of Computational Physics, 7: ,. [3] J.M. Stockie, J.A. Mackenzie, and R.D. Russell. A Moving Mesh Method for One- Dimensional Hyperbolic Conservation Laws. SIAM Journal on Scientific Computing, 5:79 83,. [4] F. Li and C.W. Shu. Reinterpretation and simplified implementation of a discontinuous Galerkin method for Hamilton-Jacobi equations. Applied Mathematics Letters, 8:4 9, 5. 8
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations Feng Zheng, Chi-Wang Shu 3 and Jianian Qiu 4 Abstract In this paper, a new type of finite difference Hermite weighted essentially
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
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 informationAn adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions
An adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions H.-Z. Tang, Tao Tang and Pingwen Zhang School of Mathematical Sciences, Peking University, Beijing
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
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 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 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 informationSemi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations
Semi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations Steve Bryson, Aleander Kurganov, Doron Levy and Guergana Petrova Abstract We introduce a new family of Godunov-type
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationIn particular, if the initial condition is positive, then the entropy solution should satisfy the following positivity-preserving property
1 3 4 5 6 7 8 9 1 11 1 13 14 15 16 17 18 19 1 3 4 IMPLICIT POSITIVITY-PRESERVING HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR CONSERVATION LAWS TONG QIN AND CHI-WANG SHU Abstract. Positivity-preserving
More informationarxiv: v1 [math.na] 2 Feb 2018
A Kernel Based High Order Eplicit Unconditionally Stable Scheme for Time Dependent Hamilton-Jacobi Equations Andrew Christlieb, Wei Guo, and Yan Jiang 3 Abstract. arxiv:8.536v [math.na] Feb 8 In this paper,
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
More informationA Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations
DOI.7/s95-4-9856-7 A Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations Liang Wu Yong-Tao Zhang Received: 7 September 23 / Revised: 23 February 24 / Accepted: 5
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
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 informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationA discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations
Journal of Computational Physics 3 (7) 39 415 www.elsevier.com/locate/jcp A discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations Yingda Cheng, Chi-Wang Shu *
More informationShooting methods for numerical solutions of control problems constrained. by linear and nonlinear hyperbolic partial differential equations
Shooting methods for numerical solutions of control problems constrained by linear and nonlinear hyperbolic partial differential equations by Sung-Dae Yang A dissertation submitted to the graduate faculty
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationNumerical discretization of tangent vectors of hyperbolic conservation laws.
Numerical discretization of tangent vectors of hyperbolic conservation laws. Michael Herty IGPM, RWTH Aachen www.sites.google.com/michaelherty joint work with Benedetto Piccoli MNCFF 2014, Bejing, 22.5.2014
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 informationA STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS *
Journal of Computational Mathematics Vol.8, No.4, 010, 55 568. http://www.global-sci.org/jcm doi:10.408/jcm.1003-m0016 A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
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 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 informationFunctions with orthogonal Hessian
Functions with orthogonal Hessian B. Dacorogna P. Marcellini E. Paolini Abstract A Dirichlet problem for orthogonal Hessians in two dimensions is eplicitly solved, by characterizing all piecewise C 2 functions
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationA Moving-Mesh Finite Element Method and its Application to the Numerical Solution of Phase-Change Problems
A Moving-Mesh Finite Element Method and its Application to the Numerical Solution of Phase-Change Problems M.J. Baines Department of Mathematics, The University of Reading, UK M.E. Hubbard P.K. Jimack
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationA central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
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 informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationA Nonlinear PDE in Mathematical Finance
A Nonlinear PDE in Mathematical Finance Sergio Polidoro Dipartimento di Matematica, Università di Bologna, Piazza di Porta S. Donato 5, 40127 Bologna (Italy) polidoro@dm.unibo.it Summary. We study a non
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationEfficient Runge-Kutta Based Local Time-Stepping Methods
Efficient Runge-utta Based Local Time-Stepping Methods by Alex Ashbourne A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Mathematics
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More information3.3.1 Linear functions yet again and dot product In 2D, a homogenous linear scalar function takes the general form:
3.3 Gradient Vector and Jacobian Matri 3 3.3 Gradient Vector and Jacobian Matri Overview: Differentiable functions have a local linear approimation. Near a given point, local changes are determined by
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 5, No., pp. 544 573 c Society for Industrial and Applied Mathematics ARBITRARILY HIGH-ORDER ACCURATE ENTROPY STABLE ESSENTIALLY NONOSCILLATORY SCHEMES FOR SYSTEMS OF CONSERVATION
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 informationThe Fast Sweeping Method. Hongkai Zhao. Department of Mathematics University of California, Irvine
The Fast Sweeping Method Hongkai Zhao Department of Mathematics University of California, Irvine Highlights of fast sweeping mehtod Simple (A Gauss-Seidel type of iterative method). Optimal complexity.
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationAccepted Manuscript. A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
Accepted Manuscript A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang, Hongkai Zhao PII: S-()- DOI:./j.jcp... Reference: YJCPH To
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationLagrangian submanifolds and generating functions
Chapter 4 Lagrangian submanifolds and generating functions Motivated by theorem 3.9 we will now study properties of the manifold Λ φ X (R n \{0}) for a clean phase function φ. As shown in section 3.3 Λ
More informationAn homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
More informationNUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES
J. Appl. Math. & Computing Vol. 21(2006), No. 1-2, pp. 603-622 NUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES MINYOUNG YUN Abstract. The level-set method is quite efficient
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationJournal of Computational Physics
Journal of Computational Physics 3 () 3 44 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A local discontinuous Galerkin method
More informationSurface Integrals (6A)
Surface Integrals (6A) Surface Integral Stokes' Theorem Copright (c) 2012 Young W. Lim. Permission is granted to cop, distribute and/or modif this document under the terms of the GNU Free Documentation
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationFifth order multi-moment WENO schemes for hyperbolic conservation laws
J. Sci. Comput. manuscript No. (will be inserted by the editor) Fifth order multi-moment WENO schemes for hyperbolic conservation laws Chieh-Sen Huang Feng Xiao Todd Arbogast Received: May 24 / Accepted:
More informationDOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE
DOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE MARTIN J. GANDER AND RONALD D. HAYNES Abstract. Moving mesh methods based on the equidistribution principle are powerful
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
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 informationcomputations. Furthermore, they ehibit strong superconvergence of DG solutions and flues for hyperbolic [,, 3] and elliptic [] problems. Adjerid et al
Superconvergence of Discontinuous Finite Element Solutions for Transient Convection-diffusion Problems Slimane Adjerid Department of Mathematics Virginia Polytechnic Institute and State University Blacksburg,
More informationWave breaking in the short-pulse equation
Dynamics of PDE Vol.6 No.4 29-3 29 Wave breaking in the short-pulse equation Yue Liu Dmitry Pelinovsky and Anton Sakovich Communicated by Y. Charles Li received May 27 29. Abstract. Sufficient conditions
More informationDuality and dynamics in Hamilton-Jacobi theory for fully convex problems of control
Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control RTyrrell Rockafellar and Peter R Wolenski Abstract This paper describes some recent results in Hamilton- Jacobi theory
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More information