High-Order Accurate, Low Numerical Diffusion Methods for Aerodynamics

Total Page:16

File Type:pdf, Size:1020Kb

High-Order Accurate, Low Numerical Diffusion Methods for Aerodynamics ARTICLE IN PRESS Progress in Aerospace Sciences 41 (2005) 192–300 www.elsevier.com/locate/paerosci High-order accurate, lownumerical diffusion methods for aerodynamics John A. Ekaterinaris Foundation for Research and Technology FORTH, Institute of Applied and Computational Mathematics IACM, FORTH/IACM, 71110 Heraklion, Crete, Greece Abstract In recent years numerical methods have been widely used to effectively resolve complex flow features of aerodynamics flows with meshes that are reasonable for today’s computers. High-order numerical methods were used mainly in direct numerical simulations and aeroacoustics. For many aeronautical applications, accurate computation of vortex- dominated flows is important because the vorticity in the flow field and the wake of swept wings at an incidence and rotor blades largely determines the distribution of loading. The main deficiency of widely available, second-order accurate methods for the accurate computation of these flows is the numerical diffusion of vorticity to unacceptable levels. Application of high-order accurate, low-diffusion numerical methods can significantly alleviate this deficiency of traditional second order methods. Furthermore, higher-order space discretizations have the potential to improve detached eddy simulation predictions of separated flows with significant unsteadiness. Recently developed high-order accurate finite-difference, finite-volume, and finite-element methods are reviewed. These methods can be used as an attractive alternative of traditional low-order central and upwind computational fluid dynamics methods for improved predictions of vortical and other complex, separated, unsteady flows. The main features of these methods are summarized, from a practical user’s point of view, their applicability and relative strength is indicated, and examples from recent applications are presented to illustrate their performance on selected problems. r 2005 Elsevier Ltd. All rights reserved. Contents 1. Introduction . 194 1.1. Prior work on high-order numerical methods . 196 2. Governing equations . 198 2.1. Conservative form of the NS equations . 198 2.2. Entropy splitting . 199 2.3. Time integration of the Euler and NS equations . 200 2.3.1. Explicit temporal schemes. 201 2.3.2. Implicit schemes. 202 2.3.3. Time accurate solutions with multigrid. 202 E-mail address: [email protected]. 0376-0421/$ - see front matter r 2005 Elsevier Ltd. All rights reserved. doi:10.1016/j.paerosci.2005.03.003 ARTICLE IN PRESS J.A. Ekaterinaris / Progress in Aerospace Sciences 41 (2005) 192–300 193 3. High-order finite-difference schemes . 203 3.1. Centered compact schemes . 203 3.2. Boundary closures for high-order schemes . 204 3.3. Simultaneous evaluation of the first and second derivative with compact schemes . 205 3.4. Modified high-order finite-difference schemes . 206 3.4.1. Upwind high-order schemes . 206 3.4.2. Dispersion-relation-preserving (DPR) scheme . 207 3.5. Spectral-type filters. 207 3.6. Characteristic-based (ACM) filters . 208 3.6.1. ENO and WENO ACM filters . 210 3.7. Finite-volume compact difference-based schemes. 211 3.8. Results with ACM filters . 212 3.9. Results with compact schemes and filters . 218 4. ENO and WENO schemes . 227 4.1. High-order reconstruction . 228 4.2. Approximation in one dimension. 228 4.3. One-dimensional conservative approximation of the derivative . 229 4.4. ENO reconstruction . 230 4.5. 1D finite-volume ENO scheme . 231 4.6. 1D finite-difference ENO scheme . 232 4.7. WENO approximation . 232 4.8. Multidimensional ENO and WENO reconstruction. 234 4.8.1. Finite-volume reconstruction for Cartesian mesh. 234 4.8.2. Two-dimensional reconstruction for triangles . 235 4.8.3. Multidimensional finite-difference ENO . 236 4.9. Optimization of WENO schemes . 236 4.10. Compact WENO approximation . 237 4.11. Hybrid compact-WENO scheme for the Euler equations . 239 4.12. Applications of ENO and WENO . 240 5. The discontinuous Galerkin (DG) method . 256 5.1. DG space discretization . 257 5.2. Element bases . 258 5.3. Arbitrary high-order DG schemes . 260 5.4. Analysis of the DG method for wave propagation . 260 5.5. Dissipative and dispersive behavior of high-order DG method . 261 5.6. Limiting of DG expansions . 261 5.6.1. Rectangular elements . 262 5.6.2. Triangular elements . 262 5.6.3. Component-wise limiters. 263 5.7. DG stabilization operator . 264 5.8. DG space discretization of the NS equations . 264 5.9. The DG variational multiscale (VMS) method . 265 5.10. Implicit time marching of DG discretizations . 269 5.11. p-type multigrid for DG . 270 5.12. Results with the DG method. 270 6. The spectral volume (SV) method . 276 6.1. Fundamentals of the SV method . 277 6.2. SV method in one dimension. 279 6.2.1. Reconstruction for SV in one dimension . 279 6.3. Polynomial reconstruction in triangular SVs . 281 6.4. Linear, quadratic, and cubic SV reconstructions . 283 6.5. Spectral volume partitions. 284 6.6. Multidimensional limiters . 284 6.7. The SV method for two-dimensional systems . 285 6.8. Multidimensional TVD and TVB limiters.
Recommended publications
  • Conservative High Order Semi-Lagrangian Finite Difference WENO Methods for Advection in Incompressible Flow Jing-Mei Qiu1 and Ch
    Conservative high order semi-Lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu1 and Chi-Wang Shu 2 Abstract In this paper, we propose a semi-Lagrangian finite difference formulation for approximat- ing conservative form of advection equations with general variable coefficients. Compared with the traditional semi-Lagrangian finite difference schemes [4, 21], which approximate the advective form of the equation via direct characteristics tracing, the scheme proposed in this paper approximates the conservative form of the equation. This essential difference makes the proposed scheme naturally conservative for equations with general variable coef- ficients. The proposed conservative semi-Lagrangian finite difference framework is coupled with high order essentially non-oscillatory (ENO) or weighted ENO (WENO) reconstructions to achieve high order accuracy in smooth parts of the solution and to capture sharp inter- faces without introducing spurious oscillations. The scheme is extended to high dimensional problems by Strang splitting. The performance of the proposed schemes is demonstrated by linear advection, rigid body rotation and swirling deformation. As the information is propa- gating along characteristics, the proposed scheme does not have CFL time step restriction of the Eulerian method, allowing for a more efficient numerical realization for many application problems. Keywords: advection in incompressible flow; conservative scheme; semi-Lagrangian meth- ods; WENO reconstruction. 1Department of Mathematical and Computer Science, Colorado School of Mines, Golden, 80401. E-mail: [email protected]. Research supported by NSF grant number 0914852 and Air Force Office of Scientific Research grant number FA9550-09-1-0344. 2Division of Applied Mathematics, Brown University, Providence, 02912.
    [Show full text]
  • PROGRAM and ABSTRACTS
    LASERENANUM E´ R I C A II Octavo Encuentro de Analisis´ Numerico´ de Ecuaciones Diferenciales Parciales Departamento de Matem´aticas,Universidad de La Serena La Serena, Chile, Enero 14 - 16, 2015 PROGRAM and ABSTRACTS Contents 1 Introduction 2 2 Wednesday, January 14 3 3 Thursday, January 15 4 4 Friday, January 16 5 5 Abstracts 6 1 1 Introduction The Octavo Encuentro de An´alisisNum´ericode Ecuaciones Diferenciales Parciales has been organized in sequential talks of 45 and 30 minutes length (40 and 25 minutes of presentation, respectively, and 5 minutes for questions and comments). All the talks will be given at Salon´ Multiuso, CETECFI (3th floor), Facultad de Ingenier´ıa, Universi- dad de la Serena. In the following pages we describe the corresponding program. In case of a multi-authored contribution, the speaker is underlined. The organizers acknowledge financial support by: • CONICYT-Chile through Project Anillo ACT 1118 ANANUM, • Departamento de Matem´aticas, Universidad de La Serena, • Centro de Modelamiento Matem´atico(CMM), Universidad de Chile, and • Centro de Investigaci´onen Ingenier´ıaMatem´atica(CI2MA), Universidad de Concepci´on. In addition, we express our recognition and gratitude to all speakers for making La Serena Num´ericaII possible. Organizing Committee Raimund B¨urger Gabriel N. Gatica Ricardo Oyarz´ua H´ectorTorres La Serena, January 2015 2 2 Wednesday, January 14 8.30-9.15 REGISTRATION 9.15-9.30 WELCOME SPEECH [Chairman: N. HEUER] 9.30-10.15 Alexandre Ern, Martin Vohral´ık: Polynomial-degree-robust a posteri- ori estimates in a unified setting. 10.15-10.45 Dante Kalise: High-order semi-Lagrangian schemes for static Hamilton- Jacobi-Bellman equations.
    [Show full text]
  • AN13 and CT13 Abstracts
    AN13 and CT13 Abstracts Abstracts are printed as submitted by the authors. Society for Industrial and Applied Mathematics 3600 Market Street, 6th Floor Philadelphia, PA 19104-2688 USA Telephone: +1-215-382-9800 Fax: +1-215-386-7999 Conference E-mail: [email protected] Conference web: www.siam.org/meetings/ Membership and Customer Service: 800-447-7426 (US & Canada) or +1-215-382-9800 (worldwide) 2 2013 SIAM Annual Meeting • SIAM Conference on Control & Its Applications Table of Contents Annual Meeting (AN13) Abstracts ...............................................3 Control & Its Applications (CT13) .............................................127 SIAM Presents Since 2008, SIAM has recorded many Invited Lectures, Prize Lectures, and selected Minisymposia from various conferences. These are available by visiting SIAM Presents (http://www.siam.org/meetings/presents.php). 2013 SIAM Annual Meeting • SIAM Conference on Control & Its Applications 3 AN13 Abstracts 4 AN13 Abstracts IC1 system including cars, buses, pedestrians, ants and molecu- Social Networks as Information Filters lar motors, which are considered as ”self-driven particles”. We recently call this interdisciplinary research on jamming Social networks, especially online social networks, are of self-driven particles as ”jamology”. This is based on driven by information sharing. But just how much informa- mathematical physics, and and includes engineering appli- tion sharing is influenced by social networks? A large-scale cations as well. In the talk, starting from the backgroud experiment measured the effect of the social network on the of this research, simple mathematical models, such as the quantity and diversity of information being shared within asymmetric simple exclusion process and the Burgers equa- Facebook. While strong ties were found to be individu- tion, are introduced as basis of all kinds of traffic flow.
    [Show full text]
  • AN17-CT17-GD17 Abstracts
    AN17-CT17-GD17 Abstracts Abstracts are printed as submitted by the authors. SIAM 2017 Events Mobile App Scan the QR code with any QR reader and download the TripBuilder EventMobile™ app to your iPhone, iPad, iTouch or Android mobile device. You can also visit www.tripbuildermedia.com/apps/siamevents Society for Industrial and Applied Mathematics 3600 Market Street, 6th Floor Philadelphia, PA 19104-2688 USA Telephone: +1-215-382-9800 Fax: +1-215-386-7999 Conference E-mail: [email protected] Conference Web: www.siam.org/meetings/ Membership and Customer Service: (800) 447-7426 (USA & Canada) or +1-215-382-9800 (worldwide) www.siam.org/meetings/ 2 2017 SIAM Annual Meeting Table of Contents AN17 Annual Meeting Abstracts ................................................3 CT17 Abstracts .........................................................................130 GD17 Abstracts .........................................................................183 SIAM Presents Since 2008, SIAM has recorded many Invited Lectures, Prize Lectures, and selected Minisymposia from various conferences. These are available by visiting SIAM Presents (http://www.siam.org/meetings/presents.php). 2017 SIAM Annual Meeting AN17 Abstracts 4 AN17 Abstracts IP1 for mapping services this is obtained from mobile phones, Genetic Consequences of Range Expansion Under while for ambulance fleets it is obtained from automatic ve- Climate Change hicle location devices. We demonstrate greatly improved accuracy relative to a system used in Bing Maps, and show Range expansion is a crucial population response to cli- the impact of our methods for improving ambulance fleet mate change. Genetic consequences are coupled to ecolog- management decisions. ical dynamics that, in turn, are driven by shifting climate conditions. We model a population with a reactiondiffu- Dawn Woodard sion system, coupled to a heterogeneous environment that Uber, USA shifts with time due to climate change.
    [Show full text]
  • 4 AN10 Abstracts
    4 AN10 Abstracts IC1 IC4 On Dispersive Equations and Their Importance in Kinematics and Numerical Algebraic Geometry Mathematics Kinematics underlies applications ranging from the design Dispersive equations, like the Schr¨odinger equation for ex- and control of mechanical devices, especially robots, to the ample, have been used to model several wave phenom- biomechanical modelling of human motion. The major- ena with the distinct property that if no boundary con- ity of kinematic problems can be formulated as systems ditions are imposed then in time the wave spreads out of polynomial equations to be solved and so fall within spatially. In the last fifteen years this field has seen an the domain of algebraic geometry. While symbolic meth- incredible amount of new and sophisticated results proved ods from computer algebra have a role to play, numerical with the aid of mathematics coming from different fields: methods such as polynomial continuation that make strong Fourier analysis, differential and symplectic geometry, an- use of algebraic-geometric properties offer advantages in alytic number theory, and now also probability and a bit of efficiency and parallelizability. Although these methods, dynamical systems. In this talk it is my intention to present collectively called Numerical Algebraic Geometry, are ap- few simple, but still representative examples in which one plicable wherever polynomials arise, e.g., chemistry, biol- can see how these different kinds of mathematics are used ogy, statistics, and economics, this talk will concentrate on in this context. applications in mechanical engineering. A brief review of the main algorithms of the field will indicate their broad Gigliola Staffilani applicability.
    [Show full text]
  • High Order ADER Schemes and GLM Curl Cleaning for a First Order Hyperbolic Formulation of Compressible Flow with Surface Tension
    High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension Simone Chiocchetti∗a, Ilya Peshkova, Sergey Gavrilyukb, Michael Dumbsera aLaboratory of Applied Mathematics, Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano 77, 38123 Trento, Italy bIUSTI UMR 7343-CNRS, D´epartement de M´ecanique, Aix-Marseille Universit´e,5 rue Enrico Fermi, 13453 Marseille, France Abstract In this work, we introduce two novel reformulations of the weakly hyperbolic model for two-phase flow with surface tension, recently forwarded by Schmidmayer et al. In the model, the tracking of phase boundaries is achieved by using a vector interface field, rather than a scalar tracer, so that the surface-force stress tensor can be expressed directly as an algebraic function of the state variables, without requiring the computation of gradients of the scalar tracer. An interesting and important feature of the model is that this interface field obeys a curl involution constraint, that is, the vector field is required to be curl-free at all times. The proposed modifications are intended to restore the strong hyperbolicity of the model, and are closely related to divergence-preserving numerical approaches developed in the field of numerical magnetohydrody- namics (MHD). The first strategy is based on the theory of Symmetric Hyperbolic and Thermodynamically Compatible (SHTC) systems forwarded by Godunov in the 60s and 70s and yields a modified system of governing equations which includes some symmetrisation terms, in analogy to the approach adopted later by Powell et al in the 90s for the ideal MHD equations.
    [Show full text]
  • Final Program
    Final Program The SIAM Conference on the Life Sciences is sponsored by the SIAM Activity Group on Life Sciences (SIAG/LS) The SIAM Activity Group on the Life Sciences was established to foster the application of mathematics to the life sciences and research in mathematics that leads to new methods and techniques useful in the life sciences. The life sciences have become quantitative as new technologies facilitate collection and analysis of vast amounts of data ranging from complete genomic sequences of organisms to satellite imagery of forest landscapes on continental scales. Computers enable the study of complex models of biological processes. The activity group brings together researchers who seek to develop and apply mathematical and computational methods in all areas of the life sciences. It provides a forum that cuts across disciplines to catalyze mathematical research relevant to the life sciences and rapid diffusion of advances in mathematical and computational methods. AN16/LS16 Mobile App Scan the QR code with any QR reader and download the TripBulder EventMobileTM app to your iPhone, iPad, iTouch, or Android devices. You can also visit www.tripbuildermedia.com/apps/siam2016events Society for Industrial and Applied Mathematics 3600 Market Street, 6th Floor Philadelphia, PA 19104-2688 USA Telephone: +1-215-382-9800 Fax: +1-215- 386-7999 Conference E-mail: [email protected] Conference Web: www.siam.org/meetings/ Membership and Customer Service: (800) 447-7426 (US & Canada) or +1-215-382-9800 (worldwide) www.siam.org/meetings 2 2016 SIAM Annual Meeting General Information Table of Contents C. David Levermore SIAM Registration Desk University of Maryland, College Park, USA General Information ...............................2 The SIAM registration desk is located on Rachel Levy Exhibitor and Sponsor Information .6 the Concourse Level of the Westin Boston Harvey Mudd College, USA Waterfront.
    [Show full text]
  • KSIAM 2020 Annual Meeting - C O N T EN TS
    KSIAM 2020 Annual Meeting - C O N T EN TS - >> Plenary Talks Stochastic Optimization in Mathematical - Moran Ki ······························································· 2 Finance II - Eun-Jae Park (KSIAM-KUMGOK Award) ······· 3 - yun Jin Jang ······················································ 45 - Hyung Hee Cho ················································· 5 - Hyungbin Park ·················································· 46 - Minsuk Kwak ····················································· 47 - Geonwoo Kim ··················································· 48 >> KSIAM Awards Deep Learning and Image Process - Yongho Choi ····················································· 12 - Han-Soo Choi ··················································· 50 - Junseok Kim ······················································ 13 - Hyomin Ahn ······················································ 51 - Sun Xiang ··························································· 14 - Geonho Hwang ················································· 52 - Yea Chan Park ·················································· 53 >> Special Session Finance·Fishery·Manufacture Industrial Mathematics Center on Big Data Math for Public Safety based on Modeling - Changsin Kim ···················································· 55 and Data analysis - Seong-Uk Nam ················································· 59 - Hyuk Kang ························································· 17 - Yuanmeng Hu ··················································· 60 -
    [Show full text]
  • Implementation and Validation of Semi-Implicit WENO Schemes Using Openfoam®
    computation Article Implementation and Validation of Semi-Implicit WENO Schemes Using OpenFOAM® Tobias Martin 1,* and Ivan Shevchuk 2 1 Department of Civil and Environmental Engineering, Norwegian University of Science and Technology, 7491 Trondheim, Norway 2 Chair of Modeling and Simulation, University of Rostock, 18059 Rostock, Germany; [email protected] * Correspondence: [email protected] Received: 3 December 2017; Accepted: 23 January 2018; Published: 24 January 2018 Abstract: In this article, the development of high-order semi-implicit interpolation schemes for convection terms on unstructured grids is presented. It is based on weighted essentially non-oscillatory (WENO) reconstructions which can be applied to the evaluation of any field in finite volumes using its known cell-averaged values. Here, the algorithm handles convex cells in arbitrary three-dimensional meshes. The implementation is parallelized using the Message Passing Interface. All schemes are embedded in the code structure of OpenFOAM® resulting in the access to a huge open-source community and the applicability to high-level programming. Several verification cases and applications of the scalar advection equation and the incompressible Navier-Stokes equations show the improved accuracy of the WENO approach due to a mapping of the stencil to a reference space without scaling effects. An efficiency analysis indicates an increased computational effort of high-order schemes in comparison to available high-resolution methods. However, the reconstruction time can be efficiently decreased when more processors are used. Keywords: CFD; high-order methods; WENO reconstruction; semi-implicit; unstructured grids 1. Introduction In recent years, open source Computational Fluid Dynamics (CFD) codes experienced an increased influence not only on ongoing science but also on the industry.
    [Show full text]
  • AMMCS2011 a Laurier Centennial Conference JULY 25 – 29 WILFRID LAURIER UNIVERSITY | WATERLOO, ONTARIO, CANADA
    The INTERNATIONAL CONFERENCE on APPLIED MATHEMATICS, MODELING and COMPUTATIONAL SCIENCE AMMCS2011 A Laurier Centennial Conference JULY 25 – 29 WILFRID LAURIER UNIVERSITY | WATERLOO, ONTARIO, CANADA Mathematics and Computation in Biological Sciences and Medicine Partial Differential and Integral Equations in Mathematical Modeling Applications of Dynamical Systems and Differential Equations Computational Physics and Chemistry Computational Algebra, Combinatorics and Optimization Mathematical Models in Social Sciences BOOK OF Computational Mechanics and Engineering Financial Mathematics and Computation ABSTRACTS Statistical Modeling in Envrionmental Sciences Computational Methods for Hyperbolic Problems Technical Design by Cameron Davidson-Pilon Applied Problems and Methods in Research and Education BOOK OF ABSTRACTS: THE INTERNATIONAL CONFERENCE ON APPLIED MATHEMATICS,MODELING AND COMPUTATIONAL SCIENCE WATERLOO,CANADA,JULY 25 - 29, 2011 LAURIER CENTENNIAL CONFERENCE: AMMCS-2011 Disclaimer: This book contains abstracts of the International Conference on Applied Mathematics, Modeling and Computational Science (AMMCS-2011). Authors are responsible for the contents and accuracy. Opinions expressed may not necessarily reflect the position of the AMMCS-2011 Scientific and Organizing Committees. Publisher: AMMCS@2011 ISBN: 2810000004711 Organizing Committee General Co-chairs Ilias Kotsireas Roderick Melnik Program Chair Brian West Treasurer Zilin Wang Student Prize Committee Chair Cristina Stoica Local Organizing Committee Francine Vinette SIAM
    [Show full text]
  • ON the USE of ENO-BASED LIMITERS for DISCONTINUOUS GALERKIN METHODS for AERODYNAMICS Pdfsubject
    Thesis presented to the Instituto Tecnol´ogico de Aeron´autica, in partial fulfillment of the requirements for the Degree of Master in Science in the Program of Aeronautical and Mechanical Engineering, Field Aerodynamics, Propulsion and Energy. Andr´eFernando de Castro da Silva ON THE USE OF ENO-BASED LIMITERS FOR DISCONTINUOUS GALERKIN METHODS FOR AERODYNAMICS Thesis approved in its final version by signatories below: Prof. Dr. Marcos Aur´elio Ortega Advisor Prof. Dr. Celso Massaki Hirata Head of the Faculty of the Department of Graduate Studies Campo Montenegro S~aoJos´edos Campos, SP - Brazil 2012 Cataloging-in Publication Data Documentation and Information Division da Silva, Andr´eFernando de Castro ON THE USE OF ENO-BASED LIMITERS FOR DISCONTINUOUS GALERKIN METHODS FOR AERODYNAMICS / Andr´eFernando de Castro da Silva. S~aoJos´edos Campos, 2012. 146f. Thesis of Master in Science { Program of Aeronautical and Mechanical Engineering. Field of Aerodynamics, Propulsion and Energy. { Aeronautics Institute of Technology, 2012. Advisor: Prof. Dr. Marcos Aur´elio Ortega. 1. Din^amica dos fluidos computacional. 2. M´etodo de Galerkin. 3. Mec^anica dos fluidos. 4. Aerodin^amica. 5. F´ısica. I.Instituto Tecnol´ogico de Aeron´autica. II. ON THE USE OF ENO-BASED LIMITERS FOR DISCONTINUOUS GALERKIN METHODS FOR AERODYNAMICS BIBLIOGRAPHIC REFERENCE DA SILVA, Andr´eFernando de Castro. ON THE USE OF ENO-BASED LIMITERS FOR DISCONTINUOUS GALERKIN METHODS FOR AERODYNAMICS. 2012. 146f. Thesis of Master in Science in Aerodynamics, Propulsion and Energy. { Instituto Tecnol´ogico de Aeron´autica, S~aoJos´edos Campos. CESSION OF RIGHTS AUTHOR NAME: Andr´eFernando de Castro da Silva PUBLICATION TITLE: ON THE USE OF ENO-BASED LIMITERS FOR DISCONTINUOUS GALERKIN METHODS FOR AERODYNAMICS.
    [Show full text]
  • Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes
    Acta Numerica (2020), pp. 1–63 c Cambridge University Press, 2020 doi:10.1017/S09624929XXXXXX Printed in the United Kingdom Essentially non-oscillatory and weighted essentially non-oscillatory schemes Chi-Wang Shu Division of Applied Mathematics Brown University Providence, RI 02912, USA chi-wang [email protected] Essentially non-oscillatory (ENO) and weighted ENO (WENO) schemes were designed for solving hyperbolic and convection-diffusion equations with pos- sibly discontinuous solutions or solutions with sharp gradient regions. The main idea of ENO and WENO schemes is actually an approximation proce- dure, aimed at achieving arbitrarily high order accuracy in smooth regions and resolving shocks or other discontinuities sharply and in an essentially non-oscillatory fashion. Both finite volume and finite difference schemes have been designed using the ENO or WENO procedure, and these schemes are very popular in applications, most noticeably in computational fluid dynam- ics, but also in other areas of computational physics and engineering. Since the main idea of the ENO and WENO schemes is an approximation proce- dure not directly related to partial differential equations (PDEs), ENO and WENO schemes also have non-PDE applications. In this paper we will sur- vey the basic ideas behind ENO and WENO schemes, discuss their properties, and present examples of their applications to different types of PDEs as well as to non-PDE problems. 2 Chi-Wang Shu CONTENTS 1 Introduction 2 2 The ENO interpolation and reconstruction proce- dures 4 3 The WENO interpolation and reconstruction pro- cedures 14 4 ENO and WENO schemes for conservation laws 23 5 ENO and WENO schemes for Hamilton-Jacobi equations 31 6 ENO and WENO schemes for convection-diffusion equations 32 7 Time discretizations 33 8 Boundary conditions 34 9 Bound-preserving WENO schemes 34 10 Applications 35 11 Concluding remarks 45 References 45 1.
    [Show full text]