Basins of attraction and critical curves for Newton-type methods in a phase equilibrium problem

Many engineering problems are described by systems of nonlinear equations, which may exhibit multiple solutions, in a challenging situation for root-finding algorithms. The existence of several solutions may give rise to complex basins of attraction for the solutions in the algorithms, with severe influence in their convergence behavior. In this work, we explore the relationship of the basins of attractions with the critical curves (the locus of the singular points of the Jacobian of the system of equations) in a phase equilibrium problem in the plane with two solutions, namely the calculation of a double azeotrope in a binary mixture. The results indicate that the conjoint use of the basins of attraction and critical curves can be a useful tool to select the most suitable algorithm for a specific problem.

[1]  G. Platt,et al.  Multiple solution algorithm with applications to robot kinematics , 2017 .

[2]  Miquel Grau-Sánchez,et al.  On the computational efficiency index and some iterative methods for solving systems of nonlinear equations , 2011, J. Comput. Appl. Math..

[3]  P. Deuflhard Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms , 2011 .

[4]  Angelo Lucia,et al.  Process simulation in the complex domain , 1993 .

[5]  Miquel Grau-Sánchez,et al.  On some computational orders of convergence , 2010, Appl. Math. Lett..

[6]  The fractal response of robust solution techniques to the stationary point problem , 1993 .

[7]  Changbum Chun,et al.  COMPARATIVE STUDY OF METHODS OF VARIOUS ORDERS FOR FINDING SIMPLE ROOTS OF NONLINEAR EQUATIONS , 2019, Journal of Applied Analysis & Computation.

[8]  Yanxing Zhao,et al.  Azeotropic refrigerants and its application in vapor compression refrigeration cycle , 2019 .

[9]  E. Zotos Determining the Newton-Raphson basins of attraction in the electromagnetic Copenhagen problem , 2017, 1704.02273.

[10]  Angelo Lucia,et al.  Simple Process Equations, Fixed-Point Methods, and Chaos , 1990 .

[11]  Alicia Cordero,et al.  On a Novel Fourth-Order Algorithm for Solving Systems of Nonlinear Equations , 2012, J. Appl. Math..

[12]  Thomas P. Wihler,et al.  THE NEWTON–RAPHSON METHOD AND ADAPTIVE ODE SOLVERS , 2011 .

[13]  W. Świȩtosławski,et al.  Azeotropy and polyazeotropy , 1963 .

[14]  Junfeng Chen,et al.  Population diversity of particle swarm optimisation algorithms for solving multimodal optimisation problems , 2018, Int. J. Comput. Sci. Eng..

[15]  Qin Ni,et al.  Incomplete Jacobian Newton method for nonlinear equations , 2008, Comput. Math. Appl..

[16]  Changbum Chun,et al.  Basin attractors for various methods , 2011, Appl. Math. Comput..

[17]  J. Wisniak,et al.  Comments on computing all the azeotropes in refrigerant mixtures through equations of state by naveed aslam and aydin k. Sunol [fluid phase equilib. 224 (2004), 97-109] , 2005 .

[18]  Zhixin Ma,et al.  IFOA: an improved forest algorithm for continuous nonlinear optimisation , 2019 .

[19]  Gustavo Mendes Platt,et al.  Azeotropy in a refrigerant system: a useful scenario to test and compare metaheuristics , 2018, Int. J. Metaheuristics.

[20]  Thomas P. Wihler,et al.  An adaptive Newton-method based on a dynamical systems approach , 2014, Commun. Nonlinear Sci. Numer. Simul..

[21]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[22]  Higinio Ramos,et al.  The application of Newton's method in vector form for solving nonlinear scalar equations where the classical Newton method fails , 2015, J. Comput. Appl. Math..

[23]  Yang Chen,et al.  Novel fruit fly algorithm for global optimisation and its application to short-term wind forecasting , 2019, Connect. Sci..

[24]  G. Nedzhibov A family of multi-point iterative methods for solving systems of nonlinear equations , 2008 .

[25]  Aline De Lima Guedes,et al.  Prediction of azeotropic behaviour by the inversion of functions from the plane to the plane , 2015 .

[26]  Ivan Petkovic,et al.  Computer methodologies for comparison of computational efficiency of simultaneous methods for finding polynomial zeros , 2020, J. Comput. Appl. Math..

[27]  L. Gardini,et al.  Contact bifurcations related to critical sets and focal points in iterated maps of the plane , 2011 .

[28]  Gustavo Mendes Platt,et al.  Double Azeotropy: calculations with Newton-like methods and continuous GRASP (C-GRASP) , 2011, Int. J. Math. Model. Numer. Optimisation.

[29]  Alicia Cordero,et al.  Pseudocomposition: A technique to design predictor-corrector methods for systems of nonlinear equations , 2012, Appl. Math. Comput..

[30]  Katherine J. Evans,et al.  Implementation of the Jacobian-free Newton-Krylov method for solving the first-order ice sheet momentum balance , 2011, J. Comput. Phys..

[31]  Haihua Zhao,et al.  Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with Jacobian-free Newton–Krylov method , 2016 .

[32]  Alicia Cordero,et al.  A fast algorithm to solve systems of nonlinear equations , 2019, J. Comput. Appl. Math..