Premium
A mixed‐variable continuously deforming finite element method for parabolic evolution problems. Part III: Numerical implementation and computational results
Author(s) -
Kececioglu Ifiyenia,
Rubinsky Boris
Publication year - 1989
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620281203
Subject(s) - finite element method , discretization , partial differential equation , ordinary differential equation , mathematics , mixed finite element method , flow (mathematics) , system of linear equations , extended finite element method , variable (mathematics) , mathematical analysis , mathematical optimization , differential equation , geometry , physics , thermodynamics
In Part I of this paper, 1 the conceptual framework of a rate variational least squares formulation of a continuously deforming mixed‐variable finite element method was presented for solving a single evolution equation. In Part II 2 a system of ordinary differential equations with respect to time was derived for solving a system of three coupled evolution equations by the deforming grid mixed‐variable least squares rate variational finite element method. The system of evolution equations describes the coupled heat flow, fluid flow and trace species transport in porous media under conditions when the flow velocities and constituent phase transitions induce sharp fronts in the solution domain. In this paper, we present the method we have adopted to integrate with respect to time the resulting spatially discretized system of non‐linear ordinary differential equations. Next, we present computational results obtained using the code in which this deforming mixed finite element method was implemented. Because several features of the formulation are novel and have not been previously attempted, the problems were selected to exercise these features with the objective of demonstrating that the formulation is correct and that the numerical procedures adopted converge to the correct solutions.