Computable scalar fields: A basis for PDE software

Abstract Partial differential equations (PDEs) are fundamental in the formulation of mathematical models of the physical world. Computer simulation of PDEs is an efficient and important tool in science and engineering. Implicit in this is the question of the computability of PDEs. In this context we present the notions of scalar and tensor fields, and discuss why these abstractions are useful for the practical formulation of solvers for PDEs. Given computable scalar fields, the operations on tensor fields will also be computable. As a consequence we get computable solvers for PDEs. The traditional numerical methods for achieving computability by various approximation techniques (e.g., finite difference, finite element or finite volume methods), all have artifacts in the form of numerical inaccuracies and various forms of noise in the solutions. We hope these observations will inspire the development of a theory for computable scalar fields, which either lets us understand why these artefacts are inherent, or provides us with better tools for constructing these basic building blocks.

[1]  Jonathan M. Lees,et al.  Boundary conditions on a finite grid: Applications with pseudospectral wave propagation , 1997 .

[2]  David W. Zingg,et al.  Comparison of High-Accuracy Finite-Difference Methods for Linear Wave Propagation , 2000, SIAM J. Sci. Comput..

[3]  J. V. Tucker,et al.  Computable functions and semicomputable sets on many-sorted algebras , 2001, Logic in Computer Science.

[4]  R. Mittet Free Surface Boundary Conditions for Elastic Staggered-Grid Modelling Schemes , 1998 .

[5]  Stephen Smale,et al.  Some Remarks on the Foundations of Numerical Analysis , 1990, SIAM Rev..

[6]  Klaus Weihrauch,et al.  Computable Analysis: An Introduction , 2014, Texts in Theoretical Computer Science. An EATCS Series.

[7]  Magne Haveraaen,et al.  On the Role of Mathematical Abstractions for Scientific Computing , 2000, The Architecture of Scientific Software.

[8]  Jochen W. Schmidt,et al.  Shape PreservingC3 Data Interpolation andC3 Histopolation with Splines on Threefold Refined Grids , 1996 .

[9]  A. Iserles,et al.  Lie-group methods , 2000, Acta Numerica.

[10]  H. Munthe-Kaas,et al.  Coordinate free numerics : Closing the gap between 'pure' and 'applied' mathematics ? , 1996 .

[11]  Dave Nichols,et al.  Short Note Modeling elastic fields across irregular boundaries , 1992 .

[12]  Richard E. Ewing,et al.  Interface conditions for acoustic and elastic wave propagation , 1991 .

[13]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[14]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[15]  David W. Zingg,et al.  High-Accuracy Finite-Difference Schemes for Linear Wave Propagation , 1996, SIAM J. Sci. Comput..

[16]  Ronald F. Boisvert,et al.  The Architecture of Scientific Software , 2001, IFIP — The International Federation for Information Processing.

[17]  Magne Haveraaen,et al.  Coordinate free programming of computational fluid dynamics problems , 2000, Sci. Program..

[18]  W. A. Mulder,et al.  Experiments with Higdon's absorbing boundary conditions for a number of wave equations , 1997 .

[19]  Rune Mittet,et al.  Free-surface boundary conditions for elastic staggered-grid modeling schemes , 2002 .

[20]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

[21]  Magne Haveraaen,et al.  Case study on algebraic software methodologies for scientific computing , 2000, Sci. Program..

[22]  J. Strikwerda Finite Difference Schemes and Partial Differential Equations , 1989 .

[23]  Kenneth W. Regan,et al.  Computability , 2022, Algorithms and Theory of Computation Handbook.

[24]  Ann K. Stehney,et al.  Geometrical Methods of Mathematical Physics: Connections for Riemannian manifolds and gauge theories , 1980 .

[25]  F. Muir,et al.  Modeling elastic fields across irregular boundaries , 1992 .

[26]  Magne Haveraaen,et al.  Use of coordinate-free numerics in elastic wave simulation , 2001 .

[27]  Magne Haveraaen,et al.  Formal Software Engineering for Computational Modelling , 1999, Nord. J. Comput..

[28]  S. Tsynkov Numerical solution of problems on unbounded domains. a review , 1998 .