Recent advances in deterministic and stochastic numerics for evolutive problems
Dajana Conte (Department of Mathematics, University of Salerno)
This contribution regards the numerical solution of evolutionary problems related to natural phenomena and physical processes modelled by functional equations of various type, with specific characteristics such as: stiffness, oscillations, presence of memory together with stochastic terms, by focusing especially on stability properties of the methods. We will describe the construction of adapted numerical methods, that is, strongly oriented to the problem and with excellent stability properties, where the stability will be also interpreted as preservation of the intrinsic qualitative characteristics of the problem itself.
For the numerical solution of stiff problems, originating for example from convection and diffusion terms in Partial Differential Equations (PDEs), methods with good stability properties can be obtained by exploiting the usage of the Jacobian matrix of the problem [1, 5]. Recent works on Time-Accurate and Highly-Stable (TASE) operators [2] show as the employment of these operators, together with explicit Runge-Kutta methods, can lead to A-stable schemes with a considerable gain in efficiency. We will focus our attention on the class of explicit and parallel peer methods [8, 9] combined with TASE operators, describing the derivation of highly stable numerical methods.
For problems with oscillating solution, we will show how the use of non-polynomial bases reveals a powerful tool for adapting numerical methods to the known behavior of the solution [3, 4, 7].
For problems with memory and influenced by stochastic phenomena, modeled by stochastic Volterra integral equations, the focus will be on the construction of numerical methods that can inherit good stability properties from the methods for stochastic differential equations [6].
The presented results have been obtained in collaboration with Beatrice Paternoster, Leila Moradi and Giovanni Pagano (University of Salerno), Raffaele D’Ambrosio (University of L’Aquila), Liviu Ixaru (Institute of Physics and Nuclear Engineering, Bucharest), Fakhrodin Mohamadi (University of Hormozgan, Iran), Ali Abdi (University of Tabriz, Iran).
A. Abdi, D. Conte, Implementation of second derivative general linear methods, Calcolo 57(3), 20 (2020).
M. Calvo, J.I. Montijano, L. Rández, A note on the stability of time–accurate and highly–stable explicit operators for stiff differential equations, J. Comput. Phys. 436, 110316 (2021).
D. Conte, R. D’Ambrosio, M.P. D’Arienzo, B. Paternoster, Multivalue mixed collocation methods, Appl. Math. Comput., 409, 126346 (2021).
D. Conte, R. D’Ambrosio, G. Giordano, L.Gr. Ixaru, B. Paternoster, User-friendly expressions of the coefficients of some exponentially fitted methods, LNCS 12249, 47–62, (2020).
D. Conte, R. D’Ambrosio, G. Pagano, B. Paternoster, Jacobian-dependent vs Jacobian-free discretizations for nonlinear differential problems, Comput. Appl. Math. 39 (3), 171 (2020).
D. Conte, R. D’Ambrosio, D. Paternoster, Improved \(\vartheta\)-methods for stochastic Volterra integral equations, Commun. Nonlinear Sci. Numer. Simul., 93, 105528 (2021).
D. Conte, L. Moradi, B. Paternoster, F. Mohamadi, Exponentially fitted two-step peer methods for oscillatory problems, Comput. Appl. Math. 39 (3), 174 (2020).
B.A. Schmitt, R. Weiner, Parallel two-step W-methods with peer variables, SIAM J. Numer. Anal. 42, 265–282 (2004).
R. Weiner, K. Biermann, B. Schmitt, H. Podhaisky: Explicit two-step peer methods, Computh. Math. Appl. 55 (4), 609-619 (2008).