Premium
Toward a reduction of mesh imprinting
Author(s) -
Lung T. B.,
Roe P. L.
Publication year - 2014
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.3941
Subject(s) - vorticity , stencil , mathematics , polygon mesh , limiter , laplacian smoothing , algorithm , robustness (evolution) , isotropy , mesh generation , computer science , geometry , vortex , finite element method , mechanics , physics , computational science , optics , gene , thermodynamics , telecommunications , biochemistry , chemistry
SUMMARY This paper is the initial investigation into a new Lagrangian cell‐centered hydrodynamic scheme that is motivated by the desire for an algorithm that resists mesh imprinting and has reduced complexity. Key attributes of the new approach include multidimensional construction, the use of flux‐corrected transport (FCT) to achieve second order accuracy, automatic determination of the mesh motion through vertex fluxes, and vorticity control. Toward this end, vorticity preserving Lax–Wendroff (VPLW) type schemes for the two‐dimensional acoustic equations were analyzed and then implemented in a FCT algorithm. Here, mesh imprinting takes the form of anisotropic dispersion relationships. If the stencil for the LW methods is limited to nine points, four free parameters exist. Two parameters were fixed by insisting that no spurious vorticity be created. Dispersion analysis was used to understand how the remaining two parameters could be chosen to increase isotropy. This led to new VPLW schemes that suffer less mesh imprinting than the rotated Richtmyer method. A multidimensional, vorticity preserving FCT implementation was then sought using the most promising VPLW scheme to address the problem of spurious extrema. A well‐behaved first order scheme and a new flux limiter were devised in the process. The flux limiter is unique in that it acts on temporal changes and does not place a priori bounds on the solution. Numerical results have demonstrated that the vorticity preserving FCT scheme has comparable performance to an unsplit MUSCL‐H algorithm at high Courant numbers but with reduced mesh imprinting and superior symmetry preservation. Copyright © 2014 John Wiley & Sons, Ltd.