A global Newton method II: Analytic centers

This paper modifies the convergence conditions of a back-tracking global Newton method announced in Goldstein (1991), making them sharper and easier to apply. A new version of the Kantorovich inequalities is presented that is simple to state and prove. An application is made to the centering problem for polytopes. Based on an idea of Ye (1989), an algorithm is given for the feasibility problem of linear inequalities.This paper is dedicated to Phil Wolfe on the occasion of his 65th birthday.     

Inequalities for linear combinations of monomials in p-Newton sequences

The partial order on monomials that corresponds to domination when evaluated at positive Newton sequences is fully understood. Here we take up the corresponding partial order on linear combinations of monomials. In part using analysis based upon the cone structure of the exponents in p-Newton sequences, an array of conditions is given for this new partial order. It appears that a characterization in general will be difficult. Within the case in which all coefficients are 1, the situation in which, for general sequence length, there are two monomials, each of length two and nonnegative integer exponents, the partial order is fully characterized. The characterization is combinatorial, in terms of indices in the monomials, and, already here there is much more than term-wise domination.     

The Spin-Coating Process: Analysis of the Free Boundary Value Problem

In this paper, an accurate model of the spin-coating process is presented and investigated from the analytical point of view. More precisely, the spin-coating process is being described as a one-phase free boundary value problem for Newtonian fluids subject to surface tension and rotational effects. It is proved that for T > 0 there exists a unique, strong solution to this problem in (0, T) belonging to a certain regularity class provided the data and the speed of rotation are small enough in suitable norms. The strategy of the proof is based on a transformation of the free boundary value problem to a quasilinear evolution equation on a fixed domain. The keypoint for solving the latter equation is the so-called maximal regularity approach. In order to pursue in this direction one needs to determine the precise regularity classes for the associated inhomogeneous linearized equations. This is being achieved by applying the Newton polygon method to the boundary symbol.     

Modified two-grid method for solving coupled Navier-Stokes/Darcy model based on Newton iteration

A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the err...     

Advancements in subsea riser analysis using quasi-rotations and the Newton–Raphson method

Beam structures undergoing finite deflections and rotations in space have extensive application in the subsea industry particularly for the analysis of holistic systems with larger numbers of mooring and riser components. In using the finite element analysis approach, there is an increasing requirement for large element sizes which preserve accuracy with regard to the coupling of axial, bending and torsion response.The authors outline a method for improving the current state of practice for the analysis of riser systems. The approach draws on the convected coordinates method, Euler–Bernoulli beam theory, the principle of virtual work and the finite element method. Two quasi-rotation measures are developed including a quasi-material rotation definition for rotational deformation relative to the convected axis of a beam and a quasi-space rotation definition to deal with the path dependent nature of rotations in three dimensions.The novel aspect of this work is to relate the rate of change of the quasi-material rotation vector along the beam axis to a linear transformation of the beam axis rate-of-rotation vector through utilising the convected coordinates axes system. In this way, incremental values of quasi-material rotation are directly linked to incremental values of nodal quasi-space rotation and a global Newton–Raphson solution technique for interconnecting beam elements is straightforward to assemble.Furthermore, this leads to accurate definitions of coupled axial, bending and torque response for beams with significant deflection. The approach has particular advantages in the analysis of subsea riser sections. Also, the accuracy of the solution is preserved for a fewer number of elements compared to alternative solutions for computationally sensitive load cases with highly non-linear loading regimes.     

A computational method for geometric optimization of enhanced heat transfer devices based upon entropy generation minimization

A computational fluid dynamics‐based optimization methodology is developed, appropriate for the geometric optimization of enhanced heat transfer devices based upon the principle of entropy generation minimization, in which the objective function is evaluated from a flow field obtained by computational simulation. A quasi‐Newton optimization procedure is employed, with computation of the objective function gradients based upon a finite difference approach. The optimization procedure is developed to be general with regard to the choice of objective function, the details of the problem under consideration, and the computational methodology employed in solving the fluid flow and heat transfer problems. A novel implementation of a Taylor series‐based procedure for the fast solution of nearby problems is presented, which is found to greatly benefit the efficiency of the present methodology. Finally, a numerical experiment is presented, illustrating the use of the present method in the geometric optimization of a practical enhanced heat transfer device on the basis of the criterion of entropy generation minimization. The optimization of the fin spacing of a simple plate fin heat sink is considered, and a comparison of the computational results with results obtained by analytical optimization based upon empirical friction factor and Nusselt number correlations is given. Copyright © 2012 John Wiley & Sons, Ltd.     

An operator splitting strategy for fluid–structure interaction problems with thin elastic structures in an incompressible Newtonian flow

We present a computational framework based on the use of the Newton and level set methods to model fluid–structure interaction problems involving elastic membranes freely suspended in an incompressible Newtonian flow. The Mooney–Rivlin constitutive model is used to model the structure. We consider an extension to a more general case of the method described in Laadhari (2017) to model the elasticity of the membrane. We develop a predictor–corrector finite element method where an operator splitting scheme separates different physical phenomena. The method features an affordable computational burden with respect to the fully implicit methods. An exact Newton method is described to solve the problem, and the quadratic convergence is numerically achieved. Sample numerical examples are reported and illustrate the accuracy and robustness of the method.     

An equivalence between two algorithms for a class of quadratic programming problems with mmatrices *

In this paper the behaviour of the gradient projection method for the Quadratic Programming problems with M-Matrices and lower bounds on variables is analysed. It is shown that the Chandbasekaran method for solving such problems is simply a realization of the Newton projection method if for the latter one the starting point has components equal to lower bounds on variables.     

Path-following barrier and penalty methods for linearly constrained problems

In the present paper some barrier and penalty methods (e.g. logarithmic barriers, SUMT, exponential penalties), which define a continuously differentiable primal and dual path, applied to linearly constrained convex problems are studied, in particular, the radius of convergence of Newton’s method depending on the barrier and penalty para-meter is estimated, Unlike using self-concordance properties the convergence bounds are derived by direct estimations of the solutions of the Newton equations. The obtained results establish parameter selection rules which guarantee the overall convergence of the considered barrier and penalty techniques with only a finite number of Newton steps at each parameter level. Moreover, the obtained estimates support scaling method which uses approximate dual multipliers as available in barrier and penalty methods