An adaptive moving mesh method is developed for the numerical solution of an enthalpy formulation of heat conduction problems with a phase change. The algorithm is based on a very simple mesh modification strategy that allows the smooth evolution of mesh nodes to track interfaces. At each time step
On The Numerical Solution of One-Dimensional PDEs Using Adaptive Methods Based on Equidistribution
โ Scribed by G. Beckett; J.A. Mackenzie; A. Ramage; D.M. Sloan
- Publisher
- Elsevier Science
- Year
- 2001
- Tongue
- English
- Weight
- 149 KB
- Volume
- 167
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
โฆ Synopsis
Numerical experiments are described that illustrate some important features of the performance of moving mesh methods for solving one-dimensional partial differential equations (PDEs). The particular method considered here is an adaptive finite difference method based on the equidistribution of a monitor function and it is one of the moving mesh methods proposed by W. Huang, Y. Ren, and R. D. Russell (1994, SIAM J. Numer. Anal. 31 709). We show how the accuracy of the computations is strongly dependent on the choice of monitor function, and we present a monitor function that yields an optimal rate of convergence. Motivated by efficiency considerations for problems in two or more space dimensions, we demonstrate a robust and efficient algorithm in which the mesh equations are uncoupled from the physical PDE. The accuracy and efficiency of the various formulations of the algorithm are considered and a novel automatic time-step control mechanism is integrated into the scheme.
๐ SIMILAR VOLUMES
This paper presents a parallel algorithm and its systolic array architecture for the BFGS (Broyden, Fletcher, Goldfarb, and Shanno) quasi-Newton method of minimizing an n-vector function. The calculation of search direction vectors and the update of approximation to Hessian matrices by the BFGS upda
The application of the uncoupling of boundary integral and finite element methods to solve exterior boundary value problems in R yields a weak formulation that contains only one boundary term. This is the so-called uncoupling term, which is determined by the boundary integral operator of the single-