Newton's iteration is studied for the numerical solution of an elliptic PDE with nonlinear boundary conditions. At each iteration of Newton's method, a conjugate gradient based decomposition method is applied to the matrix of the linearized system. The decomposition is such that all the remaining linear systems have the same constant matrix. Numerical results confirm the savings with respect to the computational cost, compared with the classical Newton method with factorization at each step.
The conjugate gradient method of Liu and Storey is an efficient minimization algorithm which uses second derivatives information, without saving matrices, by finite difference approximation. It is shown that the finite difference scheme can be removed by using a quasi-Newton approximation for computing a search direction, without loss of convergence. A conjugate gradient method based on BFGS approximation is proposed and compared with existing methods of the same class.
We present a numerical simulation of two coupled Navier-Stokes flows, using ope-rator-split-ting and optimization-based non-overlapping domain decomposition methods. The model problem consists of two Navier-Stokes fluids coupled, through a common interface, by a nonlinear transmission condition. Numerical experiments are carried out with two coupled fluids; one with an initial linear profile and the other in rest. As expected, the transmission condition generates a recirculation within the fluid...
We present a mesh adaptation method by node movement for two-dimensional linear elasticity problems with unilateral contact. The adaptation is based on a hierarchical estimator on finite element edges and the node displacement techniques use an analogy of the mesh topology with a spring network. We show, through numerical examples, the efficiency of the present adaptation method.
Download Results (CSV)