The essentials of the finite element method to solve differential equations: an illustrative case in physics
M Chavez-Paez,Enrique Gonzalez-Tovar,Guillermo Ivan Guerrero-Garcia,Guillermo Guevara-Zapata,Jose Marcelo Padron-Zamora
DOI: https://doi.org/10.1088/1402-4896/ad340e
2024-03-15
Physica Scripta
Abstract:Here we review the basic steps of the modern and accurate finite element method in its application to solve ordinary differential equations in physics. To exemplify this well-in-vogue numerical technique, we have chosen the second-order Poisson-Boltzmann equation, which is a classic equation of colloid science. Aiming to formulate a viable, but didactic, implementation of the finite element technique, we have combined a linear basis of functions, the Galerkin weighted residuals method, the Swartz-Wendroff approximation and the Picard iteration algorithm. In summary, the finite element method transforms a differential equation into a simpler system of algebraic equations for the coefficients of the approximate solution in terms of a set of basis functions. We describe the full computational realisation of the finite element procedure and, also, we examine the corresponding Poisson-Boltzmann numerical predictions for various representative conditions.
physics, multidisciplinary