Two-dimensional meshless solution of the non-linear convection diffusion reaction equation by the Local Hermitian Interpolation method

A meshless numerical scheme is developed for solving a generic version of the non-linear convection-diffusion-reaction equation in two-dimensional domains. The Local Hermitian Interpolation (LHI) method is employed for the spatial discretization and several strategies are implemented for the solution...

Full description

Autores:
Bustamante Chaverra, Carlos A
Power, Henry
Florez Escobar, Whady
Hill Betancourt, Alan F
Tipo de recurso:
Fecha de publicación:
2013
Institución:
Universidad EAFIT
Repositorio:
Repositorio EAFIT
Idioma:
eng
OAI Identifier:
oai:repository.eafit.edu.co:10784/14409
Acceso en línea:
http://hdl.handle.net/10784/14409
Palabra clave:
Radial Basis Functions
Meshless Methods
Symmetric Method
Newton Raphson
Homotopy Analysis Method
Funciones De Base Radial
Métodos Sin Malla
Método Simétrico
Newton Raphson
Método De Análisis De Homotopía
Rights
License
Acceso abierto
Description
Summary:A meshless numerical scheme is developed for solving a generic version of the non-linear convection-diffusion-reaction equation in two-dimensional domains. The Local Hermitian Interpolation (LHI) method is employed for the spatial discretization and several strategies are implemented for the solution of the resulting non-linear equation system, among them the Picard iteration, the Newton Raphson method and a truncated version of the Homotopy Analysis Method (HAM). The LHI method is a local collocation strategy in which Radial Basis Functions (RBFs) are employed to build the interpolation function. Unlike the original Kansa’s Method, the LHI is applied locally and the boundary and governing equation differential operators are used to obtain the interpolation function, giving a symmetric and non-singular collocation matrix. Analytical and Numerical Jacobian matrices are tested for the Newton-Raphson method and the derivatives of the governing equation with respect to the homotopy parameter are obtained analytically. The numerical scheme is verified by comparing the obtained results to the one-dimensional Burgers’ and two-dimensional Richards’ analytical solutions. The same results are obtained for all the non-linear solvers tested, but better convergence rates are attained with the Newton Raphson method in a double iteration scheme.