pith. sign in

arxiv: 2604.21697 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

Structure-preserving approximation for non-isothermal phase-field models in melt flow

Aaron Brunk , Dennis H\"ohn This is my paper

Pith reviewed 2026-05-09 21:15 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords finite element methodphase-field modelAllen-Cahn-Navier-Stokesentropy preservationstructure-preserving schemenon-isothermal melt flowthermodynamic consistencynumerical dissipation
0
0 comments X

The pith

A finite-element scheme for the non-isothermal Allen-Cahn-Navier-Stokes system preserves entropy production exactly.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a conforming finite-element discretization for the coupled incompressible Navier-Stokes equations and non-isothermal Allen-Cahn phase-field model that handles periodic, closed, and thermal boundary conditions. It proves that this scheme inherits exact entropy production from the continuous system while conserving total energy up to a controlled negative numerical dissipation term. Convergence in space and time is verified through tests, and the method is applied to representative melt-flow examples. A sympathetic reader cares because the structure preservation prevents unphysical entropy decrease or energy creation, enabling reliable long-time simulations of thermodynamic processes such as melting and solidification without ad-hoc fixes.

Core claim

The central claim is that the proposed conforming finite-element scheme for the non-isothermal Allen-Cahn-Navier-Stokes system exactly preserves entropy production and maintains total energy conservation up to a negative numerical dissipation, for the given boundary conditions, without requiring additional stabilization terms.

What carries the argument

The structure-preserving conforming finite-element discretization that directly inherits the thermodynamic consistency (non-negative entropy production) of the continuous non-isothermal Allen-Cahn-Navier-Stokes system.

If this is right

  • The scheme supports stable long-time integrations of melt-flow problems without artificial thermodynamic violations.
  • Spatial and temporal convergence rates can be verified numerically for the coupled system.
  • The method applies directly to periodic, closed, and thermal boundary conditions without extra stabilization.
  • Representative numerical examples confirm practical effectiveness for non-isothermal phase transitions in fluid flow.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The exact entropy preservation may extend the applicability of this discretization to other coupled phase-field systems with similar energy-dissipation structures.
  • Avoiding extra stabilization terms could simplify implementation for related non-isothermal models in materials science.
  • The negative numerical dissipation bound suggests the scheme remains dissipative overall, which may aid unconditional stability proofs in future analyses.

Load-bearing premise

The continuous non-isothermal Allen-Cahn-Navier-Stokes system obeys thermodynamic consistency with non-negative entropy production, and the chosen boundary conditions allow the discrete scheme to inherit these properties directly.

What would settle it

A computed solution in which the discrete entropy production rate becomes negative at any time step, or in which total energy increases beyond the documented negative dissipation bound, would falsify the preservation property.

read the original abstract

This work presents a conforming finite-element scheme for the non-isothermal Allen-Cahn-Navier-Stokes system, incorporating periodic, closed, and thermal boundary conditions. The system comprises the incompressible Navier-Stokes equations coupled with the non-isothermal Allen-Cahn equation, which includes a non-conserved phase-field equation and a temperature equation. The proposed numerical scheme preserves entropy production exactly and maintains total energy conservation up to a negative numerical dissipation. Convergence tests in both space and time are conducted, and representative examples are provided to demonstrate the scheme's effectiveness.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a conforming finite-element discretization of the non-isothermal Allen-Cahn-Navier-Stokes system subject to periodic, closed, and thermal boundary conditions. The central claim is that the scheme exactly preserves discrete entropy production while conserving total energy up to a negative numerical dissipation term. Convergence tests in space and time are reported together with representative numerical examples.

Significance. If the exact preservation properties are rigorously established, the work would provide a useful contribution to structure-preserving methods for coupled phase-field and incompressible flow problems. Thermodynamically consistent discretizations are valuable for long-time simulations of melt flow and non-isothermal phase transitions, where artificial entropy production or energy drift can distort physical behavior. The reported convergence tests and examples add practical support for the approach.

minor comments (3)
  1. [Abstract] Abstract: the statement that entropy production is preserved 'exactly' would benefit from a brief parenthetical reference to the specific discrete identity (e.g., the summation-by-parts or integration-by-parts relation used) so that readers can immediately locate the supporting argument.
  2. [Section 4] Section 4 (or wherever the discrete energy/entropy identities are derived): the proof that the dissipation term is non-positive should be written out explicitly rather than summarized, including the precise sign of the numerical dissipation contribution.
  3. [Numerical experiments] Numerical experiments: the convergence tables would be clearer if they reported both the observed rates and the norms employed (e.g., L2 or H1 for the phase field and velocity); the current description leaves the precise error measure ambiguous.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary and recommendation of minor revision. The referee's description accurately captures the central contribution: a conforming finite-element scheme for the non-isothermal Allen-Cahn-Navier-Stokes system that exactly preserves discrete entropy production and conserves energy up to a negative dissipation term. No major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper presents a conforming finite-element discretization explicitly constructed so that discrete analogs of the continuous thermodynamic identities (non-negative entropy production and energy balance) hold exactly or up to controlled dissipation, given the chosen boundary conditions. This structure-preserving design is verified by direct analysis of the scheme rather than by fitting parameters or re-labeling inputs as outputs. No load-bearing step reduces to a self-citation chain, self-definition, or renaming of a known result; the central preservation claims follow from the discretization choices and are independently checked via convergence tests and examples. The derivation chain remains self-contained against the continuous model and external thermodynamic consistency assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard thermodynamic consistency of the continuous model and properties of conforming finite-element spaces for incompressible flow; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The continuous non-isothermal Allen-Cahn-Navier-Stokes system satisfies thermodynamic consistency with non-negative entropy production.
    Required for the discrete scheme to inherit exact entropy preservation.

pith-pipeline@v0.9.0 · 5378 in / 1167 out tokens · 47716 ms · 2026-05-09T21:15:03.309395+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

47 extracted references · 34 canonical work pages

  1. [1]

    American Mathematical Society, Provi- dence, RI (1971)

    Rubinstein, L.I.: The Stefan Problem. American Mathematical Society, Provi- dence, RI (1971)

    1971

  2. [2]

    Hemisphere Publishing, Washington, DC (1992)

    Alexiades, V., Solomon, A.D.: Mathematical Modeling of Melting and Freezing Processes. Hemisphere Publishing, Washington, DC (1992)

    1992

  3. [3]

    Caginalp, G.: An analysis of a phase field model of a free boundary. Arch. Ration. Mech. Anal.92(3), 205–245 (1986) https://doi.org/10.1007/bf00254827 30

    work page doi:10.1007/bf00254827 1986

  4. [4]

    Caginalp, G.: Stefan and Hele-Shaw type models as asymptotic limits of the phase-field equations. Phys. Rev. A39, 5887–5896 (1989) https://doi.org/10. 1103/PhysRevA.39.5887

    1989

  5. [5]

    Provatas, K

    Provatas, N., Elder, K.: Phase-Field Methods in Materials Science and Engineer- ing. Wiley-VCH, Weinheim (2010). https://doi.org/10.1002/9783527631520

    work page doi:10.1002/9783527631520 2010

  6. [6]

    Steinbach, I.: Phase-field models in materials science. Model. Simul. Mater. Sci. Eng.17(7), 073001 (2009) https://doi.org/10.1088/0965-0393/17/7/073001

    work page doi:10.1088/0965-0393/17/7/073001 2009

  7. [7]

    Karma, A., Rappel, W.-J.: Quantitative phase-field modeling of dendritic growth in two and three dimensions. Phys. Rev. E57(4), 4323–4349 (1998) https://doi. org/10.1103/PhysRevE.57.4323

    work page doi:10.1103/physreve.57.4323 1998

  8. [8]

    Beckermann, C., Diepers, H.-J., Steinbach, I., Karma, A., Tong, X.: Modeling melt convection in phase-field simulations of solidification. JJ. Comput. Phys. 154(2), 468–496 (1999) https://doi.org/10.1006/jcph.1999.6323

    work page doi:10.1006/jcph.1999.6323 1999

  9. [9]

    Liang, C., Shen, Z., Yi, M.: Non-isothermal phase-field model of molten pool dynamics coupled with thermocapillary effect and phase change in metal addi- tive manufacturing. Phys. Fluids38(1), 013309 (2026) https://doi.org/10.1063/ 5.0310906

    2026

  10. [10]

    King, W.E., Barth, H.D., Castillo, V.M., Gallegos, G.F., Gibbs, J.W., Hahn, D.E., Kamath, C., Rubenchik, A.M.: Observation of keyhole-mode laser melting in laser powder-bed fusion additive manufacturing. J. Mater. Process. Technol. 214(12), 2915–2925 (2014) https://doi.org/10.1016/j.jmatprotec.2014.06.005

    work page doi:10.1016/j.jmatprotec.2014.06.005 2014

  11. [11]

    Acta Mater.108, 36–45 (2016) https://doi.org/10.1016/j.actamat.2016.02.014

    Khairallah, S.A., Anderson, A., Rubenchik, J., King, W.E.: Laser powder-bed fusion additive manufacturing: Physics of melting, vapor depression, and spatter. Acta Mater.108, 36–45 (2016) https://doi.org/10.1016/j.actamat.2016.02.014

    work page doi:10.1016/j.actamat.2016.02.014 2016

  12. [12]

    Cambridge University Press, Cambridge (2010)

    Gurtin, M.E., Fried, E., Anand, L.: The Mechanics and Thermodynamics of Con- tinua. Cambridge University Press, Cambridge (2010). https://doi.org/10.1017/ CBO9780511762956

    2010

  13. [13]

    Kurz, W., Fisher, D.J.: Fundamentals of Solidification, 4th ed. edn. Trans Tech Publications, Aedermannsdorf, Switzerland (1998)

    1998

  14. [14]

    Physica D135(1), 175–194 (2000) https://doi.org/10

    Anderson, D.M., McFadden, G.B., Wheeler, A.A.: A phase-field model of solid- ification with convection. Physica D135(1), 175–194 (2000) https://doi.org/10. 1016/S0167-2789(99)00109-8

    2000

  15. [15]

    Freist¨ uhler, H., Kotschote, M.: Phase-field and Korteweg-type models for the time-dependent flow of compressible two-phase fluids. Arch. Ration. Mech. Anal. 224(1), 1–20 (2016) https://doi.org/10.1007/s00205-016-1065-0 31

    work page doi:10.1007/s00205-016-1065-0 2016

  16. [16]

    ZAMP63(1), 145–169 (2011) https://doi.org/10.1007/s00033-011-0139-y

    Heida, M., M´ alek, J., Rajagopal, K.R.: On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework. ZAMP63(1), 145–169 (2011) https://doi.org/10.1007/s00033-011-0139-y

    work page doi:10.1007/s00033-011-0139-y 2011

  17. [17]

    Anderson, D.M., McFadden, G.B., Wheeler, A.A.: Diffuse-interface methods in fluid mechanics. Annu. Rev. Fluid Mech.30, 139–165 (1998) https://doi.org/10. 1146/annurev.fluid.30.1.139

    1998

  18. [18]

    ZAMP66(1), 135–147 (2014) https://doi.org/10.1007/s00033-013-0395-0

    Alessia, B., Bochicchio, I., Fabrizio, M.: Phase separation in quasi-incompressible fluids: Cahn–Hilliard model in the Cattaneo–Maxwell framework. ZAMP66(1), 135–147 (2014) https://doi.org/10.1007/s00033-013-0395-0

    work page doi:10.1007/s00033-013-0395-0 2014

  19. [19]

    Springer, New York (1999)

    Eringen, A.C.: Microcontinuum Field Theories I: Foundations and Solids. Springer, New York (1999). https://doi.org/10.1007/978-1-4612-0555-5

    work page doi:10.1007/978-1-4612-0555-5 1999

  20. [20]

    Forest, S.: Micromorphic approach for gradient elasticity, viscoplasticity, and damage. J. Eng. Mech.135(3), 117–131 (2009) https://doi.org/10.1061/(ASCE) 0733-9399(2009)135:3(117)

    work page doi:10.1061/(asce 2009

  21. [21]

    Elliott, C.M.: The Cahn–Hilliard model for the kinetics of phase separation. Math. Surf., 35–44 (1996)

    1996

  22. [22]

    Novick-Cohen, A., Pego, R.L.: Stable patterns in a viscous diffusion equation. Trans. Amer. Math. Soc.324(1), 331–351 (1991). Accessed 2026-02-13

    1991

  23. [23]

    Colli, P., Visintin, A.: On a class of doubly nonlinear evolution equations. Commun. Partial Differ. Equ.15(5), 737–756 (1990) https://doi.org/10.1080/ 03605309908820706 https://doi.org/10.1080/03605309908820706

    work page doi:10.1080/03605309908820706 1990

  24. [24]

    MRS Proc.529(1998) https://doi.org/10.1557/proc-529-39

    Eyre, D.J.: Unconditionally gradient stable time marching the Cahn-Hilliard equation. MRS Proc.529(1998) https://doi.org/10.1557/proc-529-39

    work page doi:10.1557/proc-529-39 1998

  25. [25]

    Numerical approximations of Allen-Cahn and Cahn-Hilliard equations

    Shen, J., Yang, X.: Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Syst. A28(4), 1669–1691 (2010) https://doi. org/10.3934/dcds.2010.28.1669

    work page doi:10.3934/dcds.2010.28.1669 2010

  26. [26]

    Gr¨ un, G., Klingbeil, F.: Two-phase flow with mass density contrast: Stable schemes for a thermodynamic consistent and frame-indifferent diffuse-interface model. J. Comput. Phys.257, 708–725 (2014) https://doi.org/10.1016/j.jcp.2013. 10.028

    work page doi:10.1016/j.jcp.2013 2014

  27. [27]

    Alt, W.: Existence of solutions for non-isothermal phase separation. Adv. Math. Sci. Appl.1, 319–409 (1992)

    1992

  28. [28]

    Nonlinear Anal

    Kenmochi, N., Niezg´ odka, M.: Evolution systems of nonlinear variational inequal- ities arising from phase change problems. Nonlinear Anal. Theory Methods Appl. 22(9), 1163–1180 (1994) https://doi.org/10.1016/0362-546X(94)90235-6 32

    work page doi:10.1016/0362-546x(94)90235-6 1994

  29. [29]

    Bonetti, E., Colli, P., Fremond, M.: A phase field model with ther- mal memory governed by the entropy balance. Math. Models Methods Appl. Sci.13(11), 1565–1588 (2003) https://doi.org/10.1142/S0218202503003033 https://doi.org/10.1142/S0218202503003033

    work page doi:10.1142/s0218202503003033 2003

  30. [30]

    Colli, P., Gilardi, G., Rocca, E., Schimperna, G.: On a Penrose-Fife phase-field model with nonhomogeneous Neumann boundary conditions for the tempera- ture. Differ. Integral Equ.17(5-6), 511–534 (2004) https://doi.org/10.57262/die/ 1356060345

    work page doi:10.57262/die/ 2004

  31. [31]

    Nonlinear Anal.240, 113461 (2024) https: //doi.org/10.1016/j.na.2023.113461

    Colli, P., Gilardi, G., Signori, A., Sprekels, J.: On a Cahn–Hilliard system with source term and thermal memory. Nonlinear Anal.240, 113461 (2024) https: //doi.org/10.1016/j.na.2023.113461

    work page doi:10.1016/j.na.2023.113461 2024

  32. [32]

    Discrete Contin

    Eleuteri, M., Rocca, E., Schimperna, G.: On a non-isothermal diffuse interface model for two-phase flows of incompressible fluids. Discrete Contin. Dyn. Syst. A 35(6), 2497–2522 (2015) https://doi.org/10.3934/dcds.2015.35.2497

    work page doi:10.3934/dcds.2015.35.2497 2015

  33. [33]

    Nonlinear Anal.213, 112526 (2021) https://doi.org/10.1016/j.na.2021.112526

    Lasarzik, R.: Analysis of a thermodynamically consistent Navier–Stokes–Cahn–Hilliard model. Nonlinear Anal.213, 112526 (2021) https://doi.org/10.1016/j.na.2021.112526

    work page doi:10.1016/j.na.2021.112526 2021

  34. [34]

    Electron

    Lopes, J.H., Planas, G.: Existence of solutions for a non-isothermal Navier- Stokes-Allen-Cahn system with thermo-induced coefficients. Electron. J. Differ. Equations2022(01–87), 72 (2022) https://doi.org/10.58997/ejde.2022.72

    work page doi:10.58997/ejde.2022.72 2022

  35. [35]

    Wu, H.: Well-posedness of a diffuse-interface model for two-phase incompressible flows with thermo-induced marangoni effect. Eur. J. Appl. Math.28(3), 380–434 (2017) https://doi.org/10.1017/S0956792516000322

    work page doi:10.1017/s0956792516000322 2017

  36. [36]

    Hossain, M.A., Ma, L.: Global weak solutions to a phase-field model for seawater solidification with melt convection. Part. Differ. Equ. Appl.5(4) (2024) https: //doi.org/10.1007/s42985-024-00290-2

    work page doi:10.1007/s42985-024-00290-2 2024

  37. [37]

    Paw low, I.: A thermodynamic approach to nonisothermal phase-field models. Appl. Math., 1–63 (2016) https://doi.org/10.4064/am2282-12-2015

    work page doi:10.4064/am2282-12-2015 2016

  38. [38]

    Gonzalez-Ferreiro, B., Gomez, H., Romero, I.: A thermodynamically consistent numerical method for a phase field model of solidification. Commun. Nonlin- ear Sci. Numer. Simul.19(7), 2309–2323 (2014) https://doi.org/10.1016/j.cnsns. 2013.11.016

    work page doi:10.1016/j.cnsns 2014

  39. [39]

    Brunk, A., Habrich, O., Oyedeji, T.D., Yang, Y., Xu, B.-X.: Variational approx- imation for a non-isothermal coupled phase-field system: Structure-preservation & Nonlinear stability. Comput. Methods Appl. Math.25(2), 373–396 (2025) https://doi.org/10.1515/cmam-2023-0274 33

    work page doi:10.1515/cmam-2023-0274 2025

  40. [40]

    PAMM24(2), 202400060 (2024) https://doi.org/10.1002/pamm.202400060

    Brunk, A., Schumann, D.: Nonisothermal Cahn–Hilliard Navier–Stokes system. PAMM24(2), 202400060 (2024) https://doi.org/10.1002/pamm.202400060

    work page doi:10.1002/pamm.202400060 2024

  41. [41]

    Brunk, A., Schumann, D.: Structure-Preserving Approximation for the Non- isothermal Cahn-Hilliard-Navier-Stokes System, pp. 188–197. Springer, ??? (2025). https://doi.org/10.1007/978-3-031-86173-4 19

    work page doi:10.1007/978-3-031-86173-4 2025

  42. [42]

    Rueda-G´ omez, D.A., Rueda-Fern´ andez, E.E., Villamizar-Roa, E.J.: Numerical analysis for a non-isothermal incompressible Navier–Stokes–Allen–Cahn system. J. Math. Fluid Mech.26(4) (2024) https://doi.org/10.1007/s00021-024-00898-9

    work page doi:10.1007/s00021-024-00898-9 2024

  43. [43]

    Sun, S., Li, J., Zhao, J., Wang, Q.: Structure-preserving numerical approximations to a non-isothermal hydrodynamic model of binary fluid flows. J. Sci. Comput. 83(3) (2020) https://doi.org/10.1007/s10915-020-01229-6

    work page doi:10.1007/s10915-020-01229-6 2020

  44. [44]

    Sun, Y., Wu, J., Jiang, M., Wise, S.M., Guo, Z.: A thermodynamically consis- tent phase-field model and an entropy stable numerical method for simulating two-phase flows with thermocapillary effects. Appl. Numer. Math.206, 161–189 (2024) https://doi.org/10.1016/j.apnum.2024.08.010

    work page doi:10.1016/j.apnum.2024.08.010 2024

  45. [45]

    Brunk, A., H¨ ohn, D., Luk´ aˇ cov´ a-Medvid ’ov´ a, M.: Structure-preserving approxima- tion of the non-isothermal Cahn-Hilliard system based on the entropy equation. Appl. Math. Comput.522, 129987 (2026) https://doi.org/10.1016/j.amc.2026. 129987

    work page doi:10.1016/j.amc.2026 2026

  46. [46]

    JOM72(4), 1719–1733 (2020)

    Yang, Y., K¨ uhn, P., Yi, M., Egger, H., Xu, B.-X.: Non-isothermal phase-field mod- eling of heat–melt–microstructure-coupled processes during powder bed fusion. JOM72(4), 1719–1733 (2020)

    2020

  47. [47]

    Gallou¨ et, T., Gastaldo, L., Herbin, R., Latch´ e, J.-C.: An unconditionally sta- ble pressure correction scheme for the compressible barotropic Navier-Stokes equations. ESAIM. Math. Model. Numer. Anal.42(2), 303–331 (2008) 34

    2008