A new adaptive local mesh refinement algorithm and its application on fourth order thin film flow problem

Pengtao Sun, Robert D. Russell, Jinchao Xu

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

A new adaptive local mesh refinement method is presented for thin film flow problems containing moving contact lines. Based on adaptation on an optimal interpolation error estimate in the Lp norm (1 < p ≤ ∞) [L. Chen, P. Sun, J. Xu, Multilevel homotopic adaptive finite element methods for convection dominated problems, in: Domain Decomposition Methods in Science and Engineering, Lecture Notes in Computational Science and Engineering 40 (2004) 459-468], we obtain the optimal anisotropic adaptive meshes in terms of the Hessian matrix of the numerical solution. Such an anisotropic mesh is optimal for anisotropic solutions like the solution of thin film equations on moving contact lines. Thin film flow is described by an important type of nonlinear degenerate fourth order parabolic PDE. In this paper, we address the algorithms and implementation of the new adaptive finite element method for solving such fourth order thin film equations. By means of the resulting algorithm, we are able to capture and resolve the moving contact lines very precisely and efficiently without using any regularization method, even for the extreme degenerate cases, but with fewer grid points and degrees of freedom in contrast to methods on a fixed mesh. As well, we compare the method theoretically and computationally to the positivity-preserving finite difference scheme on a fixed uniform mesh which has proven useful for solving the thin film problem.

Original languageEnglish (US)
Pages (from-to)1021-1048
Number of pages28
JournalJournal of Computational Physics
Volume224
Issue number2
DOIs
StatePublished - Jun 10 2007

Fingerprint

mesh
Thin films
thin films
electric contacts
finite element method
engineering
Hessian matrices
Finite element method
Domain decomposition methods
pulse detonation engines
lectures
norms
Sun
Contacts (fluid mechanics)
preserving
interpolation
Interpolation
sun
convection
degrees of freedom

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

@article{17a4c13ec18c4db2b70d33af61a80f68,
title = "A new adaptive local mesh refinement algorithm and its application on fourth order thin film flow problem",
abstract = "A new adaptive local mesh refinement method is presented for thin film flow problems containing moving contact lines. Based on adaptation on an optimal interpolation error estimate in the Lp norm (1 < p ≤ ∞) [L. Chen, P. Sun, J. Xu, Multilevel homotopic adaptive finite element methods for convection dominated problems, in: Domain Decomposition Methods in Science and Engineering, Lecture Notes in Computational Science and Engineering 40 (2004) 459-468], we obtain the optimal anisotropic adaptive meshes in terms of the Hessian matrix of the numerical solution. Such an anisotropic mesh is optimal for anisotropic solutions like the solution of thin film equations on moving contact lines. Thin film flow is described by an important type of nonlinear degenerate fourth order parabolic PDE. In this paper, we address the algorithms and implementation of the new adaptive finite element method for solving such fourth order thin film equations. By means of the resulting algorithm, we are able to capture and resolve the moving contact lines very precisely and efficiently without using any regularization method, even for the extreme degenerate cases, but with fewer grid points and degrees of freedom in contrast to methods on a fixed mesh. As well, we compare the method theoretically and computationally to the positivity-preserving finite difference scheme on a fixed uniform mesh which has proven useful for solving the thin film problem.",
author = "Pengtao Sun and Russell, {Robert D.} and Jinchao Xu",
year = "2007",
month = "6",
day = "10",
doi = "10.1016/j.jcp.2006.11.005",
language = "English (US)",
volume = "224",
pages = "1021--1048",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "2",

}

A new adaptive local mesh refinement algorithm and its application on fourth order thin film flow problem. / Sun, Pengtao; Russell, Robert D.; Xu, Jinchao.

In: Journal of Computational Physics, Vol. 224, No. 2, 10.06.2007, p. 1021-1048.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A new adaptive local mesh refinement algorithm and its application on fourth order thin film flow problem

AU - Sun, Pengtao

AU - Russell, Robert D.

AU - Xu, Jinchao

PY - 2007/6/10

Y1 - 2007/6/10

N2 - A new adaptive local mesh refinement method is presented for thin film flow problems containing moving contact lines. Based on adaptation on an optimal interpolation error estimate in the Lp norm (1 < p ≤ ∞) [L. Chen, P. Sun, J. Xu, Multilevel homotopic adaptive finite element methods for convection dominated problems, in: Domain Decomposition Methods in Science and Engineering, Lecture Notes in Computational Science and Engineering 40 (2004) 459-468], we obtain the optimal anisotropic adaptive meshes in terms of the Hessian matrix of the numerical solution. Such an anisotropic mesh is optimal for anisotropic solutions like the solution of thin film equations on moving contact lines. Thin film flow is described by an important type of nonlinear degenerate fourth order parabolic PDE. In this paper, we address the algorithms and implementation of the new adaptive finite element method for solving such fourth order thin film equations. By means of the resulting algorithm, we are able to capture and resolve the moving contact lines very precisely and efficiently without using any regularization method, even for the extreme degenerate cases, but with fewer grid points and degrees of freedom in contrast to methods on a fixed mesh. As well, we compare the method theoretically and computationally to the positivity-preserving finite difference scheme on a fixed uniform mesh which has proven useful for solving the thin film problem.

AB - A new adaptive local mesh refinement method is presented for thin film flow problems containing moving contact lines. Based on adaptation on an optimal interpolation error estimate in the Lp norm (1 < p ≤ ∞) [L. Chen, P. Sun, J. Xu, Multilevel homotopic adaptive finite element methods for convection dominated problems, in: Domain Decomposition Methods in Science and Engineering, Lecture Notes in Computational Science and Engineering 40 (2004) 459-468], we obtain the optimal anisotropic adaptive meshes in terms of the Hessian matrix of the numerical solution. Such an anisotropic mesh is optimal for anisotropic solutions like the solution of thin film equations on moving contact lines. Thin film flow is described by an important type of nonlinear degenerate fourth order parabolic PDE. In this paper, we address the algorithms and implementation of the new adaptive finite element method for solving such fourth order thin film equations. By means of the resulting algorithm, we are able to capture and resolve the moving contact lines very precisely and efficiently without using any regularization method, even for the extreme degenerate cases, but with fewer grid points and degrees of freedom in contrast to methods on a fixed mesh. As well, we compare the method theoretically and computationally to the positivity-preserving finite difference scheme on a fixed uniform mesh which has proven useful for solving the thin film problem.

UR - http://www.scopus.com/inward/record.url?scp=34248525081&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34248525081&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2006.11.005

DO - 10.1016/j.jcp.2006.11.005

M3 - Article

AN - SCOPUS:34248525081

VL - 224

SP - 1021

EP - 1048

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -