65 research outputs found
Nonsmooth Newton methods for set-valued saddle point problems
We present a new class of iterative schemes for large scale set-valued saddle point problems as arising, e.g., from optimization problems in the presence of linear and inequality constraints. Our algorithms can be regarded either as nonsmooth Newton-type methods for the nonlinear Schur complement or as Uzawa-type iterations with active set preconditioners. Numerical experiments with a control constrained optimal control problem and a discretized Cahn–Hilliard equation with obstacle potential illustrate the reliability and efficiency of the new approach
Composing Scalable Nonlinear Algebraic Solvers
Most efficient linear solvers use composable algorithmic components, with the
most common model being the combination of a Krylov accelerator and one or more
preconditioners. A similar set of concepts may be used for nonlinear algebraic
systems, where nonlinear composition of different nonlinear solvers may
significantly improve the time to solution. We describe the basic concepts of
nonlinear composition and preconditioning and present a number of solvers
applicable to nonlinear partial differential equations. We have developed a
software framework in order to easily explore the possible combinations of
solvers. We show that the performance gains from using composed solvers can be
substantial compared with gains from standard Newton-Krylov methods.Comment: 29 pages, 14 figures, 13 table
Improving Pseudo-Time Stepping Convergence for CFD Simulations With Neural Networks
Computational fluid dynamics (CFD) simulations of viscous fluids described by
the Navier-Stokes equations are considered. Depending on the Reynolds number of
the flow, the Navier-Stokes equations may exhibit a highly nonlinear behavior.
The system of nonlinear equations resulting from the discretization of the
Navier-Stokes equations can be solved using nonlinear iteration methods, such
as Newton's method. However, fast quadratic convergence is typically only
obtained in a local neighborhood of the solution, and for many configurations,
the classical Newton iteration does not converge at all. In such cases,
so-called globalization techniques may help to improve convergence.
In this paper, pseudo-transient continuation is employed in order to improve
nonlinear convergence. The classical algorithm is enhanced by a neural network
model that is trained to predict a local pseudo-time step. Generalization of
the novel approach is facilitated by predicting the local pseudo-time step
separately on each element using only local information on a patch of adjacent
elements as input. Numerical results for standard benchmark problems, including
flow through a backward facing step geometry and Couette flow, show the
performance of the machine learning-enhanced globalization approach; as the
software for the simulations, the CFD module of COMSOL Multiphysics is
employed
- …