A second order finite element method for the one‐dimensional Stefan problem
暂无分享,去创建一个
We describe a finite element method for the one-dimensional Stefan problem. The elements are quadrilaterals of the space-time plane which are determined at each time-step in relation with the position of the free boundary. The method appears as a generalization of the classical Crank-Nicolson scheme, since it is identical to this scheme in the case of rectangular elements; it has the advantage of providing a simple and accurate determination of the free boundary. Numerical experiments show that the order of accuracy is equal to 2.
[1] Gunter H. Meyer,et al. On a free interface problem for linear ordinary differential equations and the one-phase Stefan problem , 1970 .
[2] Philippe G. Ciarlet,et al. THE COMBINED EFFECT OF CURVED BOUNDARIES AND NUMERICAL INTEGRATION IN ISOPARAMETRIC FINITE ELEMENT METHODS , 1972 .
[3] Jim Douglas,et al. On the numerical integration of a parabolic differential equation subject to a moving boundary condition , 1955 .