Lecture 16. Wednesday, May 25, φ t + V φ = 0 (1) In upwinding, we choose the discretization of the spatial derivative based on the sign of u.
|
|
- Anissa Lamb
- 6 years ago
- Views:
Transcription
1 Lecture 16 Wednesday, May 25, 2005 Consider the linear advection equation φ t + V φ = 0 (1) or in 1D φ t + uφ x = 0. In upwinding, we choose the discretization of the spatial derivative based on the sign of u. if u > 0, we use φ x (e.g., D φ = φi φi 1 x ) if u < 0, we use φ + x (e.g., D + φ = φi+1 φi x ) if u = 0, do nothing Using ENO we construct higher order approximations to φ x φ + x. Here we describe the WENO (weighted essentially non-oscillatory) method, which gives a better approximation of φ x than ENO. The following is Osher Fedkiw, Hamilton-Jacobi WENO When calculating (φ x ) i, the third order accurate HJ ENO scheme uses a subset of {φ i 3, φ i 2, φ i 1, φ i, φ i+1, φ i+2 } depending on how the stencil is chosen. In fact, there are exactly three possible HJ ENO approximations to (φ x ) i. Defining v 1 = D φ i 2, v 2 = D φ i 1, v 3 = D φ i, v 4 = D φ i+1 v 5 = D φ i+2 allows us to write φ 1 x = v 1 3 7v v 3 6 (2) φ 2 x = v v v 4 3 (3) 1
2 φ 3 x = v v 4 6 v 5 6 (4) as the three potential HJ ENO approximations to φ x. The goal of HJ ENO is the choose the single approximation with the least error by choosing the smoothest possible polynomial interpolation of φ. In [3], Liu et. al. pointed out that the ENO philosophy of picking exactly one of three cidate stencils is overkill in smooth regions where the data is well behaved. They proposed a Weighted ENO (WENO) method that takes a convex combination of the three ENO approximations. Of course, if any of the three approximations interpolate across a discontinuity, it is given minimal weight in the convex combination in order to minimize its contribution the resulting errors. Otherwise, in smooth regions of the flow, all three approximations are allowed to make a significant contribution in a way that improves the local accuracy from third order to fourth order. Later, Jiang Shu [2] improved the WENO method by choosing the convex combination weights in order to obtain the optimal fifth order accuracy in smooth regions of the flow. In [1], following the work on HJ ENO in [5], Jiang Peng extended WENO to the Hamilton-Jacobi framework. This Hamilton-Jacobi WENO or HJ WENO scheme turns out to be very useful for solving equation 1 as it reduces the errors by more than an order of magnitude over the third order accurate HJ ENO scheme for typical applications. The HJ WENO approximation of (φ x ) i is a convex combination of the approximations in equations 2, 3 4 given by φ x = ω 1 φ 1 x + ω 2 φ 2 x + ω 3 φ 3 x (5) where the 0 ω k 1 are the weights with ω 1 +ω 2 +ω 3 = 1. The key observation for obtaining high order accuracy in smooth regions is that weights of ω 1 =.1, ω 2 =.6 ω 3 =.3 give the optimal fifth order accurate approximation to φ x. While this is the optimal approximation, it is only valid in smooth regions. In nonsmooth regions, this optimal weighting can be very inaccurate we are better off with digital (ω k = 0 or ω k = 1) weights that choose a single approximation to φ x, i.e. the HJ ENO approximation. Reference [2] pointed out that setting ω 1 =.1 + O(( x) 2 ), ω 2 =.6 + O(( x) 2 ) ω 3 =.3 + O(( x) 2 ) still gives the optimal fifth order accuracy in smooth regions. In order to see this, we rewrite these as ω 1 =.1 + C 1 ( x) 2, ω 2 =.6 + C 2 ( x) 2 ω 3 =.3 + C 3 ( x) 2 plug them into equation 5 to obtain.1φ 1 x +.6φ 2 x +.3φ 3 x (6) C 1 ( x) 2 φ 1 x + C 2 ( x) 2 φ 2 x + C 3 ( x) 2 φ 3 x (7) 2
3 as the two terms that are added up to give the HJ WENO approximation to φ x. The term given by equation 6 is the optimal approximation which gives the exact value of φ x plus an O(( x) 5 ) error term. Thus, if the term given by equation 7 is O(( x) 5 ), then the entire HJ WENO approximation is O(( x) 5 ) in smooth regions. To see that this is the case, first note that each of the HJ ENO φ k x approximations gives the exact value of φ x, denoted φ E x, plus an O(( x) 3 ) error term (in smooth regions). Thus, the term in equation 7 is C 1 ( x) 2 φ E x + C 2 ( x) 2 φ E x + C 3 ( x) 2 φ E x (8) plus an O(( x) 2 )O(( x) 3 ) = O(( x) 5 ) term. Since, each of the C k are O(1), as is φ E x, this appears to be an O(( x) 2 ) term at first glance. However, since ω 1 +ω 2 +ω 3 = 1, we have C 1 +C 2 +C 3 = 0 implying that the term in equation 8 is identically zero. Thus, the HJ WENO approximation is O(( x) 5 ) in smooth regions. Note that [3] obtained only fourth order accuracy, since they chose ω 1 =.1 + O( x), ω 2 =.6 + O( x) ω 3 =.3 + O( x). In order to define the weights, ω k, we follow [1] estimate the smoothness of the stencils in equations 2, 3 4 as S 1 = (v 1 2v 2 + v 3 ) (v 1 4v 2 + 3v 3 ) 2 (9) S 2 = (v 2 2v 3 + v 4 ) (v 2 v 4 ) 2 (10) S 3 = (v 3 2v 4 + v 5 ) (3v 3 4v 4 + v 5 ) 2 (11) respectively. Using these smoothness estimates, we define α 1 = α 2 =.1 (S 1 + ɛ) 2 (12).6 (S 2 + ɛ) 2 (13) α 3 =.3 (S 3 + ɛ) 2 (14) with ɛ = 10 6 max { v1, 2 v2, 2 v3, 2 v4, 2 v5 2 } (15) where the term is set to avoid division by zero in the definition of the α k. This value for epsilon was first proposed by Fedkiw et al. [4] where the 3
4 first term is a scaling term that aids in the balance between the optimal fifth order accurate stencil the digital HJ ENO weights. In the case that φ is an approximate signed distance function, the v k which approximate φ x are approximately equal to one so that the first term in equation 15 can be set to This first term can then absorb the second term yielding ɛ = 10 6 in place of equation 15. Since the first term in equation 15 is only a scaling term, it is valid to make this v k 1 estimate in multidimensions as well. A smooth solution has small variation leading to small S k. If the S k are small enough compared to ɛ then equations 12, become α 1.1ɛ 2, α 2.6ɛ 2 α 3.3ɛ 2 exhibiting the proper ratios for the optimal fifth order accuracy. That is, normalizing the α k to obtain the weights ω 1 = α 1 (16) ω 2 = ω 3 = α 2 (17) α 3 (18) gives (approximately) the optimal weights of ω 1 =.1, ω 2 =.6 ω 3 =.3 when the S k are small enough to be dominated by ɛ. Nearly optimal weights are also obtained when the S k are larger than ɛ as long as all the S k are approximately the same size as is the case for sufficiently smooth data. On the other h, if the data is not smooth as indicated by large S k, then the corresponding α k will be small compared to the other α k s giving that particular stencil limited influence. If two of the S k are relatively large, then their corresponding α k s will both be small the scheme will rely most heavily on a single stencil similar to the digital behavior of HJ ENO. In the unfortunate instance that all three of the S k are large, the data is poorly conditioned none of the stencils are particularly useful. This case is problematic for the HJ ENO method as well, but fortunately it usually occurs only locally in space time allowing the methods to repair themselves after the situation subsides. (φ + x ) i is constructed with a subset of {φ i 2, φ i 1, φ i, φ i+1, φ i+2, φ i+3 }. Defining v 1 = D + φ i+2, v 2 = D + φ i+1, v 3 = D + φ i, v 4 = D + φ i 1 v 5 = D + φ i 2 allows us to use equations 2, 3 4 as the three HJ ENO approximations to (φ + x ) i. Then the HJ WENO convex combination is given by equation 5 with the weights given by equations 16, References [1] Jiang, G.-S. Peng, D., Weighted ENO Schemes for Hamilton Jacobi Equations, SIAM J. Sci. Comput. 21, (2000). 4
5 [2] Jiang, G.-S. Shu, C.-W., Efficient Implementation of Weighted ENO Schemes, J. Comput. Phys. 126, (1996). [3] Liu, X.-D., Osher, S. Chan, T., Weighted Essentially Non-Oscillatory Schemes, J. Comput. Phys. 126, (1996). [4] Fedkiw, R., Merriman, B., Osher, S., Simplified Upwind Discretization of Systems of Hyperbolic Conservation Laws Containing Advection Equations, J. Comput. Phys. 157, (2000). [5] Osher, S. Shu, C.-W., High Order Essentially Non-Oscillatory Schemes for Hamilton-Jacobi Equations, SIAM J. Numer. Anal. 28, (1991). 5
A 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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationFinite 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 informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
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 information2 Introduction Standard conservative schemes will admit non-physical oscillations near material interfaces, which can be avoided by the application of
A Conservative Approach to the Multiphase Euler Equations without Spurious Pressure Oscillations Xu-Dong Liu Ronald P. Fedkiw Stanley Osher June 8, 998 Abstract In [4], the pressure evolution equation
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationSemi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations
IMA Journal of Numerical Analysis 005 5, 113 138 doi: 10.1093/imanum/drh015 Semi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations STEVE BRYSON Programme in Scientific
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 new type of Hermite WENO schemes for Hamilton-Jacobi equations
A new type of Hermite WENO schemes for Hamilton-Jacobi equations Jun Zhu and Jianxian Qiu Abstract In this paper, a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes is designed
More informationSOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD
ITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N. 38 2017 (184 193 184 SOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD A. Tongxia Li Department of Mathematics Hulunbeier
More informationEssentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws
NASA/CR-97-0653 ICASE Report No. 97-65 Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws Chi-Wang Shu Brown University Institute for Computer
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 numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
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 informationA Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem
A Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem Frédéric Gibou Ronald Fedkiw April 27, 2004 Abstract In this paper,
More informationStrong Stability Preserving Time Discretizations
AJ80 Strong Stability Preserving Time Discretizations Sigal Gottlieb University of Massachusetts Dartmouth Center for Scientific Computing and Visualization Research November 20, 2014 November 20, 2014
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
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 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 informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
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 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 informationNUMERICAL SOLUTION TO NONLINEAR TRICOMI EQUATION USING WENO SCHEMES
Eighth Mississippi State - UAB Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conf. 9 ), pp. 35 44. ISSN: 7-669. URL: http://ejde.math.txstate.edu
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationA Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws
J Sci Comput (011) 46: 359 378 DOI 10.1007/s10915-010-9407-9 A Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws Fei Teng Li Yuan Tao Tang Received: 3 February 010 / Revised:
More informationA Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations
J Sci Comput (00) 45: 404 48 DOI 0.007/s095-009-9340-y A Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations Fengyan Li Sergey Yakovlev Received: 6 August 009 / Revised: November 009 /
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 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 informationDevice Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models
Device Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models Carlo Cercignani, Irene M. Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract This paper describes benchmark comparisons for
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 informationApplicability of the High Field Model: A Preliminary Numerical Study
Applicability of the High Field Model: A Preliminary Numerical Study Carlo Cercignani,IreneM.Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract In a companion presentation, we have discussed the theory
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 informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
More informationSynchronization of Weighted Essentially Non-Oscillatory Methods
45th AIAA Aerospace Sciences Meeting and Exhibit, 8 11 January 2007, Reno NV Synchronization of Weighted Essentially Non-Oscillatory Methods Ellen M. Taylor and M. Pino Martín Princeton University, Princeton
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationBuilding Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations
Building Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations Adam Oberman University of Texas, Austin http://www.math.utexas.edu/~oberman Fields Institute Colloquium January 21,
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 informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
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 informationA LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS
A LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS SHI JIN AND STANLEY J. OSHER Abstract. We develop a level set method for the
More informationccsd , version 1-9 Oct 2005
THE EXPLICIT SIMPLIFIED INTERFACE METHOD FOR COMPRESSIBLE MULTICOMPONENT FLOWS. BRUNO LOMBARD AND ROSA DONAT Abstract. This paper concerns the numerical approximation of the Euler equations for multicomponent
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationICES REPORT Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws
ICES REPORT 17-29 October 2017 Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws by Todd Arbogast, Chieh-Sen Huang, Xikai Zhao The Institute for Computational Engineering
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
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 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 informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
More informationEfficient Implementation of Weighted ENO Schemes*
JOURNAL OF COMPUTATIONAL PHYSICS 126, 202 228 (1996) ARTICLE NO. 0130 Efficient Implementation of Weighted ENO Schemes* GUANG-SHAN JIANG AND CHI-WANG SHU Division of Applied Mathematics, Brown University,
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 informationConstruction of Shapes Arising from the Minkowski Problem Using a Level Set Approach
Construction of Shapes Arising from the Minkowski Problem Using a Level Set Approach Li-Tien Cheng June 15, 2006 Department of Mathematics University of California, San Diego 9500 Gilman Dr. La Jolla,
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
More informationarxiv: v3 [math.na] 31 Jan 2017
Embedded WENO: a design strategy to improve existing WENO schemes arxiv:1607.05019v3 [math.na] 31 Jan 017 Bart S. van Lith a,1, Jan H.M. ten Thije Boonamp a, Wilbert L. IJzerman a,b a Department of Mathematics
More informationCompact Central WENO Schemes for Multidimensional Conservation Laws
arxiv:math/9989v [math.na] 3 Nov 999 Compact Central WENO Schemes for Multidimensional Conservation Laws Doron Levy Gabriella Puppo Giovanni Russo Abstract We present a new third-order central scheme for
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
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 informationThree-dimensional high order large-scale numerical investigation on the air explosion
Three-dimensional high order large-scale numerical investigation on the air explosion Cheng Wang a, JianXu Ding a, Chi-Wang Shu b,, Tao Li a a State Key Laboratory of Explosion Science and Technology,
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationA higher order frozen Jacobian iterative method for solving Hamilton-Jacobi equations
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (26), 62 6227 Research Article A higher order frozen Jacobian iterative method for solving Hamilton-Jacobi equations Ebraheem O. Alzahrani a,
More informationA PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS
A PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS ROBERTO CROCE, MICHAEL GRIEBEL, AND MARC ALEXANDER SCHWEITZER Abstract. In this paper we present
More informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationTHIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS
Palestine Journal of Mathematics Vol. 5(Special Issue: ) (6), 8 37 Palestine Polytechnic University-PPU 6 THIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS A.
More informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
More informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationDownloaded 12/03/12 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No., pp. 93 938 c Society for Industrial and Applied Mathematics Downloaded /3/ to 35..9.73. Redistribution subject to SIAM license or copyright; see http://www.siam.org/journals/ojsa.php
More informationA semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws
A semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c a Department of Applied Mathematics, National Sun Yat-sen University,
More informationA New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods
A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods Raymond J. Spiteri Steven J. Ruuth Technical Report CS-- May 6, Faculty of Computer Science 65 University Ave.,
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 informationDepartment of Mathematics University of California Santa Barbara, Santa Barbara, California,
The Ghost Fluid Method for Viscous Flows 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 9435 Email:fedkiw@cs.stanford.edu Xu-Dong Liu Department of Mathematics
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
More informationThe Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1
The Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 94305
More informationModeling and simulation of urban air pollution from the dispersion of vehicle exhaust: A continuum modeling approach. Abstract
Modeling and simulation of urban air pollution from the dispersion of vehicle exhaust: A continuum modeling approach Liangze Yang 1, Tingting Li 2, S.C.Wong 3, Chi-Wang Shu 4 and Mengping Zhang 5 Abstract
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 informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationMaximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments
Maximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments By Xiangxiong Zhang 1 and Chi-Wang Shu 1 Department of Mathematics, Brown University,
More informationModeling full-mach-range cavitating flow with sharp interface model
Center for Turbulence Research Proceedings of the Summer Program 28 83 Modeling full-mach-range cavitating flow with sharp interface model By X. Y. Hu, N. A. Adams, E. Johnsen AND G. Iaccarino In this
More informationTHE CLOSEST POINT METHOD FOR TIME-DEPENDENT PROCESSES ON SURFACES
THE CLOSEST POINT METHOD FOR TIME-DEPENDENT PROCESSES ON SURFACES by Colin B. Macdonald B.Sc., Acadia University, 200 M.Sc., Simon Fraser University, 2003 a thesis submitted in partial fulfillment of the
More informationIMA Preprint Series # 1993
AN ADAPTIVE METHOD WITH RIGOROUS ERROR CONTROL FOR THE HAMILTON-JACOBI EQUATIONS. PART I: THE ONE-DIMENSIONAL STEADY STATE CASE By Bernardo Cockburn and Bayram Yenikaya IMA Preprint Series # 993 ( September
More informationSimplied Discretization of Systems of Hyperbolic. Conservation Laws Containing Advection Equations. Ronald P. Fedkiw. Barry Merriman.
Simplied Discretization of Systems of Hyperbolic onservation Laws ontaining dvection Equations Ronald P. Fedkiw arry Merriman Stanley Osher ugust 4, 999 bstract The high speed ow of complex materials can
More information