Newton linearization of the Navier-Stokes equations for flow computations using a fully coupled finite volume method

Abstract Newton linearization is implemented for the discretized advection terms in the Navier-Stokes momentum equations. A modified Newton linearization algorithm is developed by analyzing how to properly account for mass conservation implicitly in the linearization. The numerical performance of the modified Newton linearization algorithm is investigated in terms of solution stability, convergence behaviour, and accuracy. A set of standard fluid flow test problems is solved using a pressure-based cell-centered fully coupled finite volume solver linearized by Picard, standard Newton, and modified Newton methods. The accuracy of numerical solutions is assessed by comparisons with provided benchmark solutions and those obtained by the present coupled solver linearized by the Picard method. The numerical results showed the modified Newton linearization algorithm converged to a solution with a rate up to 22 times higher than that of Picard linearization algorithm for the steady flow fields for lid-driven cavity test problems. The modified Newton linearization technique is also evaluated by computations of steady flow over a backward facing step and unsteady flow around a circular cylinder. The developed algorithm reduces the iteration number of the linearization cycle required up to 10 times for the backward facing step flow problem. Also, the unsteady flow computations are performed by modified Newton linearization algorithm with a lower number of iterations in comparison with the Picard linearization algorithm. The numerical experiments indicate that the modified Newton linearization algorithm maintains the solution accuracy with respect to the Picard linearization algorithm.

[1]  Vincenzo Casulli,et al.  A Nested Newton-Type Algorithm for Finite Volume Methods Solving Richards' Equation in Mixed Form , 2010, SIAM J. Sci. Comput..

[2]  Ionut Danaila,et al.  A Newton method with adaptive finite elements for solving phase-change problems with natural convection , 2014, J. Comput. Phys..

[3]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[4]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[5]  Sang Dong Kim,et al.  Newton's method for the Navier-Stokes equations with finite-element initial guess of stokes equations , 2006, Comput. Math. Appl..

[6]  S. Vakilipour,et al.  A Coupled Pressure-Based Co-Located Finite-Volume Solution Method for Natural-Convection Flows , 2012 .

[7]  F. W. Schmidt,et al.  USE OF A PRESSURE-WEIGHTED INTERPOLATION METHOD FOR THE SOLUTION OF THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON A NONSTAGGERED GRID SYSTEM , 1988 .

[8]  Dongjoo Kim,et al.  A Second-Order Time-Accurate Finite Volume Method for Unsteady Incompressible Flow on Hybrid Unstructured Grids , 2000 .

[9]  Gerry E. Schneider,et al.  Control Volume Finite-Element Method for Heat Transfer and Fluid Flow Using Colocated Variables— 1. Computational Procedure , 1987 .

[10]  Denz Lee,et al.  A hybrid adaptive gridding procedure for recirculating fluid flow problems , 1993 .

[12]  D. Spalding.  A novel finite difference formulation for differential expressions involving both first and second derivatives , 1972 .

[13]  S. Majumdar.  Role of underrelaxation in momentum interpolation for calculation of flow with nonstaggered grids , 1988 .

[14]  Chang-Ho Lee,et al.  USE OF THE MOMENTUM INTERPOLATION METHOD FOR FLOWS WITH A LARGE BODY FORCE , 2003 .

[15]  G. D. Raithby,et al.  A critical evaluation of upstream differencing applied to problems involving fluid flow , 1976 .

[16]  Laszlo Fuchs,et al.  Numerical and experimental study of driven flow in a polar cavity , 1985 .

[17]  G. D. Raithby,et al.  Skew upstream differencing schemes for problems involving fluid flow , 1976 .

[18]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[19]  Raytcho D. Lazarov,et al.  A finite volume element method for a non-linear elliptic problem , 2005, Numer. Linear Algebra Appl..

[20]  D. Gartling.  A test problem for outflow boundary conditions—flow over a backward-facing step , 1990 .

[21]  S. Vakilipour,et al.  Developing implicit pressure‐weighted upwinding scheme to calculate steady and unsteady flows on unstructured grids , 2008 .

[22]  Iuliu Sorin Pop,et al.  A modified L-scheme to solve nonlinear diffusion problems , 2019, Comput. Math. Appl..

[23]  Reza Hosseini,et al.  High-order pole-treatment in cylindrical coordinates for incompressible flow simulations with finite-difference collocated schemes , 2015, J. Comput. Phys..

[24]  G. Raithby,et al.  TREATMENT OF NON-LINEARITIES IN THE NUMERICAL SOLUTION OF THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS , 1986 .

[25]  K. E. Torrance,et al.  Upstream-weighted differencing schemes and their application to elliptic problems involving fluid flow , 1974 .

[26]  Bo Yu,et al.  DISCUSSION ON MOMENTUM INTERPOLATION METHOD FOR COLLOCATED GRIDS OF INCOMPRESSIBLE FLOW , 2002 .

[27]  Zhen F. Tian,et al.  A compact scheme for the streamfunction-velocity formulation of the 2D steady incompressible Navier-Stokes equations in polar coordinaes , 2013, J. Sci. Comput..

[28]  F. Radu,et al.  A study on iterative methods for solving Richards’ equation , 2015, Computational Geosciences.

[29]  A finite-difference method of high-order accuracy for the solution of transient nonlinear diffusive-convective problem in three dimensions , 2014 .

[30]  Christian Masson,et al.  Development of a fully coupled control‐volume finite element method for the incompressible Navier–Stokes equations , 2004 .

[31]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[32]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[33]  Tony W. H. Sheu,et al.  Newton linearization of the incompressible Navier–Stokes equations , 2004 .

[34]  Frédéric Hecht,et al.  A full discretisation of the time-dependent Boussinesq (buoyancy) model with nonlinear viscosity , 2018, Calcolo.

[35]  O. Posdziech,et al.  A systematic approach to the numerical calculation of fundamental quantities of the two-dimensional flow over a circular cylinder , 2007 .

[36]  Graeme J. Pettet,et al.  A finite volume method with linearisation in time for the solution of advection-reaction-diffusion systems , 2014, Appl. Math. Comput..

[37]  F. Thiele,et al.  ON THE CALCULATION OF MOMENTUM, HEAT, AND MASS TRANSFER IN LAMINAR AND TURBULENT BOUNDARY LAYER FLOWS ALONG A VAPORIZING LIQUID FILM , 1983 .

[38]  M. Mohammadi,et al.  Developing a physical influence upwind scheme for pressure‐based cell‐centered finite volume methods , 2018, International Journal for Numerical Methods in Fluids.

[39]  D. Balsara,et al.  A divergence‐free semi‐implicit finite volume scheme for ideal, viscous, and resistive magnetohydrodynamics , 2018, International journal for numerical methods in fluids.

[40]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[41]  George D. Raithby,et al.  PREDICTION OF HEAT AND FLUID FLOW IN COMPLEX GEOMETRIES USING GENERAL ORTHOGONAL COORDINATES , 1986 .

[42]  Seok-Ki Choi,et al.  Note on the use of momentum interpolation method for unsteady flows , 1999 .

[43]  S. M. H. Karimian,et al.  Advances in control-volume-based finite-element methods for compressible flows , 1994 .

[44]  F. Moukalled,et al.  A coupled finite volume solver for the solution of incompressible flows on unstructured grids , 2009, J. Comput. Phys..

[45]  R. Kessler,et al.  Comparison of finite-volume numerical methods with staggered and colocated grids , 1988 .

[46]  Masoud Darbandi,et al.  Conceptual linearization of Euler governing equations to solve high speed compressible flow using a pressure‐based method , 2008 .

[47]  Homer F. Walker,et al.  NITSOL: A Newton Iterative Solver for Nonlinear Systems , 1998, SIAM J. Sci. Comput..

[48]  Jiten C. Kalita,et al.  A 4OEC scheme for the biharmonic steady Navier-Stokes equations in non-rectangular domains , 2015, Comput. Phys. Commun..

[49]  A. Mamun,et al.  Nonlinear Effects on the Convergence of Picard and Newton Iteration Methods in the Numerical Solution of One-Dimensional Variably Saturated–Unsaturated Flow Problems , 2017 .

[50]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .