2024/10/16 by Golo A. Wimmer, Wimmer, Golo A., Ben S. Southworth +3
Computer Science · Engineering · Materials Science · #Advanced Mathematical Modeling in Engineering #Advanced Numerical Methods in Computational Mathematics #FOS: Mathematics #Numerical Analysis (math.NA) #Solidification and crystal growth phenomena
paper · pdf · doi:10.48550/arxiv.2410.13087
openalex publication_date 2024/10/16 · openalex created_date 2024/10/20 · openalex updated_date 2026/08/01
We present a novel spatial discretization for the Cahn-Hilliard equation including transport. The method is given by a mixed discretization for the two elliptic operators, with the phase field and chemical potential discretized in discontinuous Galerkin spaces, and two auxiliary flux variables discretized in a divergence-conforming space. This allows for the use of an upwind-stabilized discretization for the transport term, while still ensuring a consistent treatment of structural properties including mass conservation and energy dissipation. Further, we couple the novel spatial discretization to an adaptive time stepping method in view of the Cahn-Hilliard equation's distinct slow and fast time scale dynamics. The resulting implicit stages are solved with a robust preconditioning strategy, which is derived for our novel spatial discretization based on an existing one for continuous Galerkin based discretizations. Our overall scheme's accuracy, robustness, efficient time adaptivity as well as structure preservation and stability with respect to advection dominated scenarios are demonstrated in a series of numerical tests.