Analysis of the Picard-Newton iteration for the incompressible Boussinesq equations
<|reference_start|>Analysis of the Picard-Newton iteration for the incompressible Boussinesq equations: We study the Picard-Newton iteration for the incompressible Boussinesq equations, which is a two-step iteration resulting from the composition of the Picard and Newton iterations. We prove that this iterative method retains Newton's quadratic convergence but has less restrictive sufficient conditions for convergence than Newton and also is unconditionally stable under a small data condition. In this sense, Picard-Newton can be considered as a Newton iteration that is nonlinearly preconditioned with Picard. Our numerical tests illustrate this quadratic convergence and stability on benchmark problems. Furthermore, the tests show convergence for significantly higher Rayleigh number than both Picard and Newton, which illustrates the larger convergence basin of Picard-Newton that the theory predicts. We also introduce Anderson acceleration into the Picard step in our Picard-Newton numerical tests, and this enables convergence for even higher Rayleigh number.<|reference_end|>
@article{hawkins2024analysis,
    title={Analysis of the Picard-Newton iteration for the incompressible
  Boussinesq equations},
    author={Elizabeth Hawkins},
    journal={arXiv preprint arXiv:2408.16872},
    year={2024},
    archivePrefix={arXiv},
    eprint={2408.16872},
    primaryClass={math.NA cs.NA}
}