Abstract

We introduce an integration scheme for parabolic problems. Our parallelizable method uses adaptive h p -finite elements in space, and finite differences in time. The strategy can also be combined with a classical finite element method parallelization technique based on domain decomposition. We verified the performance of our method against two different benchmarks, in both two-dimensional (model problem on an L-shaped domain) and three-dimensional (Pennes bioheat equation) settings. Results show a significant speedup in computational time when compared with the sequential version of the algorithm. Moreover, we develop a mathematical framework to analyze similar schemes which include h p spatial adaptivity. Our framework describes error propagation rigorously, and as such allows to analyze convergence properties of these mixed methods.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

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.