Search papers, labs, and topics across Lattice.
This paper introduces a two-stage Newton initial guess strategy that leverages learned features from parameter-space sampling and a database of precomputed solutions to accelerate the convergence of nonlinear parametrized PDE solvers. By constructing two complementary reduced spaces鈥攐ne for converged states and another for intermediate Newton increments鈥攖he method effectively predicts surrogate solutions for unseen parameters and computes a residual-minimizing correction. Numerical experiments demonstrate significant reductions in Newton iterations and overall CPU time, showcasing the method's applicability to a wide range of large-scale nonlinear problems.
Learning features from past solutions can slash Newton's method iterations and drastically speed up convergence for nonlinear PDEs.
It is well known that Newton's method converges faster when the initial guess is closer to a root of a system of nonlinear equations. In this paper, a two-stage Newton initial guess strategy is proposed by learning features from a parameter-space sampling and a database of precomputed solutions. The method uses discrete Newton trajectories to construct two complementary reduced spaces: a solution feature space, built from converged states, and a corrective search direction feature space, built from intermediate Newton increments. For an unseen parameter, a regression model is used to predict a surrogate solution approximation. Then, in a second step, a residual-minimizing correction is computed using a dedicated GMRES-based approach. The resulting state is then used as an initial guess for the high-fidelity Newton method, which completes convergence. The corrective step is computationally inexpensive since it only requires residual evaluations and the solution of a small least-squares problem. The methodology is weakly intrusive once the high-fidelity residual fields and a script-based programming interface are available. This strategy reduces the number of Newton iterations and decreases the overall CPU time. Numerical experiments on representative PDE problems show quantifiable speedups compared with standalone surrogate initialization. Significant speedups are observed. This generic approach can be applied to a broad class of large-scale nonlinear problems.
