High-resolution scheme
{{short description|Scheme used in the numerical solution of partial differential equations}}
High-resolution schemes are used in the numerical solution of partial differential equations where high accuracy is required in the presence of shocks or discontinuities. They have the following properties:
- Second- or higher-order spatial accuracy is obtained in smooth parts of the solution.
- Solutions are free from spurious oscillations or wiggles.
- High accuracy is obtained around shocks and discontinuities.
- The number of mesh points containing the wave is small compared with a first-order scheme with similar accuracy.
General methods are often not adequate for accurate resolution of steep gradient phenomena; they usually introduce non-physical effects such as smearing of the solution or spurious oscillations. Since publication of Godunov's order barrier theorem, which proved that linear methods cannot provide non-oscillatory solutions higher than first order ({{harvnb|Godunov|1954}}, {{harvnb|Godunov|1959}}), these difficulties have attracted much attention and a number of techniques have been developed that largely overcome these problems. To avoid spurious or non-physical oscillations where shocks are present, schemes that exhibit a Total Variation Diminishing (TVD) characteristic are especially attractive. Two techniques that are proving to be particularly effective are MUSCL (Monotone Upstream-Centered Schemes for Conservation Laws), a flux/slope limiter method ({{harvnb|van Leer|1979}}, {{harvnb|Hirsch|1991}}, {{harvnb|Anderson|Tannehill|Pletcher|2016}}, {{harvnb|Laney|1998}}, {{harvnb|Toro|1999}}) and the WENO (Weighted Essentially Non-Oscillatory) method ({{harvnb|Shu|1998}}, {{harvnb|Shu|2009}}). Both methods are usually referred to as high resolution schemes (see diagram).
MUSCL methods are generally second-order accurate in smooth regions (although they can be formulated for higher orders) and provide good resolution, monotonic solutions around discontinuities. They are straightforward to implement and are computationally efficient.
For problems comprising both shocks and complex smooth solution structure, WENO schemes can provide higher accuracy than second-order schemes along with good resolution around discontinuities. Most applications tend to use a fifth order accurate WENO scheme, whilst higher order schemes can be used where the problem demands improved accuracy in smooth regions.
The method of holistic discretisation systematically analyses subgrid scale dynamics to algebraically construct closures for numerical discretisations that are both accurate to any specified order of error in smooth regions, and automatically adapt to cater for rapid grid variations through the algebraic learning of subgrid structures ({{harvnb|Roberts|2003}}). [https://profajroberts.github.io/holistic1.html A web service analyses any PDE in a class that may be submitted].
See also
References
{{refbegin}}
- {{cite thesis | last = Godunov | first = Sergei K. | date = 1954 | title = Ph.D. Dissertation: Different Methods for Shock Waves | publisher = Moscow State University}}
- {{cite journal |first=Sergei K. |last=Godunov |title=A Difference Scheme for Numerical Solution of Discontinuous Solution of Hydrodynamic Equations |journal=Mat. Sbornik |volume=47 |pages=271–306 |year=1959 }} translated US Joint Publ. Res. Service, JPRS 7226, 1969.
- {{cite journal |first=A. |last=Harten |title=High Resolution Schemes for Hyperbolic Conservation Laws |journal=J. Comput. Phys. |volume=49 |issue=3 |pages=357–393 |year=1983 |doi=10.1016/0021-9991(83)90136-5 |bibcode=1983JCoPh..49..357H |hdl=2060/19830002586 |hdl-access=free }}
- {{cite book |first=Charles |last=Hirsch |series=Numerical Computation of Internal and External Flows |volume=2 |title=Computational Methods for Inviscid and Viscous Flows |date=1991 |publisher=Wiley |isbn=978-0-471-92452-4}}
- {{cite book |first=Culbert B. |last=Laney |title=Computational Gasdynamics |url=https://books.google.com/books?id=pcwLAQAAQBAJ |date=1998 |publisher=Cambridge University Press |isbn=978-1-107-39360-8}}
- {{cite journal |first=A.J. |last=Roberts |title=A holistic finite difference approach models linear dynamics consistently |journal=Mathematics of Computation |volume=72 |issue= 241|pages=247–262 |year=2003 |doi=10.1090/S0025-5718-02-01448-5 |arxiv=math/0003135 |s2cid=11525980 }}
- {{cite book |first=C-W. |last=Shu |chapter=Essentially Non-oscillatory and Weighted Essential Non-oscillatory Schemes for Hyperbolic Conservation Laws. In: Cockburn |editor-link= Alfio Quarteroni|editor-first=Alfio |editor-last=Quarteroni |title=Advanced Numerical Approximation of Nonlinear Hyperbolic Equations |publisher=Springer |series=Lecture Notes in Mathematics |volume=1697 |year=1998 |isbn=978-3-540-49804-9 |pages=325–432 |doi=10.1007/BFb0096355|hdl=2060/19980007543 }}
- {{cite journal |first=C-W. |last=Shu |title=High Order Weighted Essentially Non-oscillatory Schemes for Convection Dominated Problems |journal=SIAM Review |volume=51 |issue=1 |pages=82–126 |year=2009 |doi=10.1137/070679065 |bibcode=2009SIAMR..51...82S }}
- {{cite book |first1=Dale |last1=Anderson |first2=John C. |last2=Tannehill |first3=Richard H. |last3=Pletcher |title=Computational Fluid Mechanics and Heat Transfer |url=https://books.google.com/books?id=-GTOBQAAQBAJ |date=2016 |publisher=Taylor & Francis |isbn=978-1-4665-7830-2 |edition=3rd}}
- {{cite book |first=Eleuterio F. | last = Toro |title=Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction |url=https://books.google.com/books?id=BSD4CAAAQBAJ |date=1999 |publisher=Springer |edition=2nd | isbn=978-3-662-03915-1 }}
- {{cite journal |first=B. |last=van Leer |title=Towards the ultimate conservative difference scheme V. A second order sequel to Godunov's method |journal=J. Comput. Phys. |volume=32 |issue=1 |pages=101–136 |year=1979 |doi=10.1016/0021-9991(79)90145-1 |bibcode=1979JCoPh..32..101V }}
{{refend}}
{{Numerical PDE}}
{{DEFAULTSORT:High-Resolution Scheme}}