Abstract
In this paper, we propose new adaptive local refinement (ALR) strategies for first-order system least-squares finite elements in conjunction with algebraic multigrid methods in the context of nested iteration. The goal is to reach a certain error tolerance with the least amount of computational cost and nearly uniform distribution of the error over all elements. To accomplish this, the refinement decision at each refinement level is determined based on optimizing efficiency measures that take into account both error reduction and computational cost. Two efficiency measures are discussed: predicted error reduction and predicted computational cost. These methods are first applied to a two-dimensional (2D) Poisson problem with steep gradients, and the results are compared with the threshold-based methods described in [W. Dörfler, SIAM J. Numer. Anal., 33 (1996), pp. 1106–1124]. Next, these methods are applied to a 2D reduced model of the incompressible, resistive magnetohydrodynamic equations. These equations are used to simulate instabilities in a large aspect-ratio tokamak. We show that, by using the new ALR strategies on this system, we are able to resolve the physics using only 10 percent of the computational cost used to approximate the solutions on a uniformly refined mesh within the same error tolerance.
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.