An approximate factorization/least squares solution method for a mixed finite element approximation of the Cahn-Hilliard equation
Autor: | Roland Glowinski, D. A. Trevas, Edward J. Dean |
---|---|
Rok vydání: | 1996 |
Předmět: | |
Zdroj: | Japan Journal of Industrial and Applied Mathematics. 13:495-517 |
ISSN: | 1868-937X 0916-7005 |
DOI: | 10.1007/bf03167260 |
Popis: | We discuss in this article the numerical solution of the Cahn-Hilliard equation modelling the spinodal decomposition of binary alloys. The numerical methodology combines a second-order finite difference time discretization with a mixed finite element space approximation and a least squares formulation based on an approximate factorization of a fourth-order elliptic operator which appears in the numerical model. The least squares problem—which is linear—is solved by a preconditioned conjugate gradient algorithm. The results of numerical experiments illustrate the possibilities of the methods discussed in this article. |
Databáze: | OpenAIRE |
Externí odkaz: |