Abstract
The numerical integration of time‐dependent PDEs, especially of Advection Diffusion Reaction type, for two and three spatial variables (in short, 2D and 3D problems) in the MoL framework is considered. The spatial discretization is made by using Finite Differences and the time integration is carried out by means of the L‐stable, third order formula known as the two stage Radau IIA method. The main point for the solution of the large dimensional ODEs is not to solve the stage values of the Radau method until convergence (because the convergence is very slow on the stiff components), but only giving a very few iterations and take as advancing solution the latter stage value computed. The iterations are carried out by using the Approximate Matrix Factorization (AMF) coupled to a Newton‐type iteration (SNI) as indicated in [5], which turns out in an acceptably cheap iteration, like Alternating Directions Methods (ADI) of Peaceman and Rachford (1955). Some stability results for the whole process (AMF)‐(SNI) and a local error estimate for an adaptive time‐integration are also given. Numerical results on two standard PDEs are presented and some conclusions about our method and other well‐known solvers are drawn.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.