Pith. sign in

REVIEW 2 major objections 6 minor 28 references

A Note on Hyperbolic Relaxation of the Navier-Stokes-Cahn-Hilliard system for incompressible two-phase flow

T0 review · 2 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper proves that a three-parameter relaxed approximation of the Navier-Stokes-Cahn-Hilliard system has a hyperbolic first-order sub-system in one dimension, so its inviscid part can be solved with standard conservation-law methods.

desk verdict A clean, correct hyperbolicity proof for a new formally derived relaxation of NSCH, but the approximation claim is only formal and the numerics are preliminary. read the letter →

arxiv 2412.11904 v2 pith:OE55YHXZ submitted 2024-12-16 math.AP physics.comp-phphysics.flu-dyn

Jens Keim , Hasel-Cicek Konan , Christian Rohde This is my paper · ORCID
classification math.APphysics.comp-phphysics.flu-dyn MSC 35L6535Q3576T06
keywords two-phaseflowdiffuse-interfacemodelingNavier-Stokes-Cahn-Hilliardsystemfirst-orderhyperbolicrelaxationhyperbolicityentropy-entropyfluxpairartificialcompressibilityCahn-Hilliardapproximation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to establish that the Navier-Stokes-Cahn-Hilliard system for two incompressible immiscible fluids, which is burdened by a divergence constraint and a fourth-order Cahn-Hilliard operator, can be approximated by a first-order system that is hyperbolic in its inviscid part. The proposed relaxed friction-type model combines artificial compressibility for the velocity-pressure block, a friction variable that formally replaces the chemical-potential gradient, and a relaxation variable $\omega$ that formally replaces the Laplacian-of-phase term. In one space dimension the main theorem shows that, whenever the relaxation parameter $\beta$ lies below a threshold set by the most negative curvature of the free energy, the homogeneous first-order subsystem admits a convex entropy-entropy flux pair and hence has real wave speeds and a complete eigenbasis. If this is right, the inviscid two-phase dynamics can be attacked with the toolbox of hyperbolic conservation laws, including shock-capturing finite-volume schemes, while the full system still dissipates a natural energy. The paper's own caveat is that the convergence of relaxed solutions to the original NSCH solutions as the three parameters vanish is not proven.

What carries the argument

The load-bearing object is the entropy/entropy-flux pair $(\eta,q)$ from (4.12), with $\eta(Q)=\tfrac{\alpha}{2}p^2+\tfrac12 u^2+W(c)+\tfrac{1}{2\beta}c^2+\tfrac{\delta}{2}v^2$ and $q(Q)=pu+\tfrac12 u^3+c(W'(c)+\beta^{-1}c)u+(W'(c)+\beta^{-1}c)v$. The entropy is convex because its Hessian is $\operatorname{diag}(\alpha,1,W''(c)+\beta^{-1},\delta)$, and the condition on $\beta$ keeps the phase-field entry positive for every $c\in[-1,1]$. The flux $f$ is engineered so that $\nabla\eta^T Df=\nabla q^T$, the compatibility relation that makes $q$ the entropy flux; the classical convex-extension fact then certifies hyperbolicity without solving the quartic characteristic polynomial. Structurally, the relaxed system (4.1) is assembled from three devices: artificial compressibility adds $p_t$ to enforce incompressibility in the $\alpha\to0$ limit, the friction variable $v$ is designed to approach the negative chemical-potential gradient as $\delta\to0$, and the elliptic relation $-\gamma\Delta\omega+\beta^{-1}\omega=\beta^{-1}c$ is designed so that $\beta^{-1}(\omega-c)$ approaches $\gamma\Delta c$ as $\beta\to0$.

What would settle it

Compute the $\beta\to0$ limit of $\beta^{-1}(\omega^\beta-c^\beta)$ in the one-dimensional model problem (4.3)-(4.4); if it is not $\gamma c_{xx}$, the relaxation is not reproducing the Cahn-Hilliard capillarity and the model is not an approximation of two-phase flow, even though the hyperbolicity theorem itself could remain true.

Watch

Extended reading notes

Core claim

The central discovery is Theorem 4.3: for the one-dimensional reduction $Q=(p,u,c,v)$, the first-order conservation system $Q_t+f(Q)_x=0$ with flux $f(Q)=(u/\alpha,\tfrac34 u^2+p+G(c),cu+v,\delta^{-1}(W'(c)+\beta^{-1}c))$ and $G'(c)=cW''(c)+\beta^{-1}c$ is hyperbolic whenever $\alpha,\delta>0$ and $\beta < -\left(\min_{c\in[-1,1]}\min\{W''(c),0\}\right)^{-1}$. The proof exhibits the convex entropy $\eta(Q)=\tfrac{\alpha}{2}p^2+\tfrac12 u^2+W(c)+\tfrac{1}{2\beta}c^2+\tfrac{\delta}{2}v^2$ and the matching entropy flux $q(Q)=pu+\tfrac12 u^3+c(W'(c)+\beta^{-1}c)u+(W'(c)+\beta^{-1}c)v$. Convexity of $\eta$ is exactly the condition on $\beta$, and the compatibility relation $\nabla\eta^T Df=\nabla q^T$ is verified algebraically. By the classical convex-extension principle, this yields four real eigenvalues and a complete eigenbasis at every admissible state. The paper also proves an energy-dissipation law for the full relaxed system and presents numerical characteristic plots showing the four wave speeds arranged in an outer pair controlled by $\delta$ and an inner pair controlled by $\alpha$, shifted by the mean velocity $u$.

Load-bearing premise

The whole approach rests on unproven formal limits: as $\alpha$, $\delta$, and $\beta$ tend to zero, the relaxed solutions are expected to converge to the original Navier-Stokes-Cahn-Hilliard solutions, and the authors explicitly state that this convergence question remains open.

Editorial extensions

If this is right

  • The inviscid part of the relaxed model can be discretized with standard finite-volume schemes for hyperbolic conservation laws; the paper demonstrates this in one dimension with a MUSCL-Hancock/Rusanov solver.
  • The four real wave speeds form an outer pair controlled by $\delta$ and an inner pair controlled by $\alpha$, so any accurate scheme has to resolve the fast scales that appear as either parameter goes to zero.
  • The full relaxed system dissipates the energy $E_\epsilon$ at the rate $-\int_\Omega(|v^\epsilon|^2+\nu|\nabla u^\epsilon|^2)\,dx$, matching the NSCH dissipation structure in the formal limit.
  • Convexity of the entropy makes entropy-stable discretizations of the subsystem available, giving a stability handle for two-phase flow computations.
  • The one-dimensional theorem is the stated first step toward the multi-dimensional and asymptotic-preserving time-integration schemes announced as forthcoming work.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the theorem's threshold on $\beta$ involves only the free energy's curvature, the proof covers both quartic double-well potentials and logarithmic potentials; in practice $\beta$ should be chosen below the reciprocal of the steepest downward curvature of $W''$.
  • The characteristic polynomial in (4.15) is independent of $p$ and $v$, which suggests those two fields may be linearly degenerate or at least wave-speed-neutral; a tailored Riemann solver could exploit this to reduce the four-wave problem to the $(c,u)$ plane.
  • A testable route to make the approximation rigorous is to prove the formal relaxation limit $\beta^{-1}(\omega^\epsilon-c^\epsilon)\to\gamma\Delta c$ in the full NSCH setting; if it holds, the same convex entropy could support multi-dimensional stability estimates and turn the relaxed model into a convergent numerical method for two-phase flow.
  • The fast-slow wave splitting seen in the numerical spectra suggests an implicit-explicit time integrator with the fast waves treated implicitly; this is a concrete recipe for the stiffness the paper identifies as the main obstacle.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. The paper proposes a first-order hyperbolic relaxation of the inviscid Navier-Stokes-Cahn-Hilliard (NSCH) system. It introduces a three-parameter (α, δ, β) system, equation (4.1), that combines artificial compressibility for the incompressibility constraint, a friction-type approximation for the phase-field flux, and a relaxation of the third-order capillarity term through an auxiliary elliptic variable ω. The main analytical results are energy-dissipation identities for the friction and relaxed systems (Theorems 3.2 and 4.2) and, in one space dimension, a proof that the first-order subsystem (4.8) is hyperbolic via an explicit entropy/entropy-flux pair under condition (4.11) (Theorem 4.3). A characteristic polynomial is derived and studied numerically. The final sections contain 1D numerical experiments comparing the relaxed model with Cahn-Hilliard reference solutions for droplet relaxation, spinodal decomposition, Ostwald ripening, and Riemann-type shock-tube setups. The paper explicitly states in Section 6 that the question of convergence of solutions of the approximate system to solutions of the NSCH system remains open.

Significance. If Theorem 4.3 is accepted, the paper provides a rigorous and self-contained basis for applying standard hyperbolic conservation-law solvers to the first-order subsystem (4.8). The explicit entropy/entropy-flux pair and the energy-dissipation identities are clean, correct, and potentially useful for future numerical analysis. The numerical experiments give encouraging qualitative evidence that the relaxation may approximate the Cahn-Hilliard/NSCH dynamics in one dimension. However, the paper's stated purpose is to provide an approximative system for NSCH, and that approximation property rests entirely on formal singular limits (Remarks 3.1 and 4.1) plus numerical evidence; no convergence theorem or quantitative asymptotic-preserving estimate is supplied. The hyperbolicity result itself appears correct, but its interpretation as a foundation for a numerical method for two-phase flow depends on the open convergence question.

major comments (2)
  1. [Sections 3.1, 4.1, and 6] The central claim that (4.1) approximates the NSCH system (2.1) rests on three unproven singular limits: α→0 (artificial compressibility), δ→0 (friction), and β→0 with β^{-1}(ω-c)→γΔc. Remark 4.1 only states formal expectations, and Section 6 explicitly says that the convergence of U^ε to U remains open. No compactness, relative-entropy, or asymptotic-preserving estimate is supplied that would justify the reduction of the u- and v-equations in (4.1) to the NSCH force and chemical-potential flux. Since the experiments in Section 5 are presented as validation of convergence to the CH/NSCH solution, the manuscript should either provide a convergence proof in a simplified setting or explicitly recalibrate its claims and label the full-system experiments as heuristic.
  2. [Theorem 4.3 and Sections 5.1, 5.3.4] The hyperbolicity statement is proved only for c∈[-1,1] and under condition (4.11), but several numerical experiments leave that regime. In Section 5.1, initial data (5.2) set c=2, outside Q. In Section 5.3.4, the potential in (5.6) has W''(c)=12c^2-36c+26, which is positive on [-1,1]; hence min_{c∈[-1,1]} min{W''(c),0}=0 and the right-hand side of (4.11) is -∞, so the condition as written is not satisfied by β=10 (or by any positive β). The intended sufficient condition is evidently W''(c)+β^{-1}>0 for all c in the relevant interval, which would cover (5.6); the theorem should be reformulated accordingly, and the numerical claims should be restricted to the regime in which hyperbolicity is actually proved.
minor comments (6)
  1. [Notation] The paper switches inconsistently between ε and ϵ for the parameter vector and for superscripts of approximate solutions; for example, (3.1) uses ϵ while (4.1) uses ε, and the proof of Theorem 4.2 contains vϵ. The notation should be unified.
  2. [Section 3.1] The text refers to 'Theorem 2 of Section 3.1', but the theorem is numbered Theorem 3.2; the cross-reference should be corrected.
  3. [Figure captions 4.3-4.10] Several captions appear to have reversed interval endpoints or duplicated values: for example, Figures 4.3 and 4.4 list α∈[0.01,0.001] instead of an increasing interval, and Figures 4.9 and 4.10 state α∈[0.01,0.01]. These should be corrected.
  4. [Section 5.1] The convex free energy is written as W(c)=c−1; if W(c)=c^{-1} is intended, the domain of c and the meaning of the condition on [-1,1] should be stated explicitly, especially because the initial data in (5.2) use c=2.
  5. [Section 5.3.4] There are typographical errors: 'friciton' should be 'friction', and 'only surfs the purpose' should be 'only serves the purpose'.
  6. [Figure 5.14] The bottom-left plot is described as showing the 'expected first-order convergence' of ω to c, but the plot only displays ||cε-ωε||_2 against β; since the formal limit in Remark 4.1 concerns β^{-1}(ω-c)→γΔc, the presentation should either plot the relevant scaled quantity or describe the observed rate as numerical evidence rather than an expected result.

Circularity Check

0 steps flagged · score 1.0 of 10

No load-bearing circularity; hyperbolicity theorem is self-contained, and the approximation claims are explicitly conditional.

full rationale

The central result (Theorem 4.3) is proved by direct computation rather than by importing the conclusion. The authors define the entropy η in (4.12), compute its Hessian, impose (4.11) to make it positive definite, verify the compatibility relation (∇η)^T Df = (∇q)^T using the definition G' = cW'' + β^{-1}c, and then invoke the classical Friedrichs–Lax/Godunov criterion. No quantity is fitted to data and no prior theorem of the same authors is used to establish hyperbolicity. The choice of the quadratic term 1/(2β)c^2 in η is an ansatz designed to make the Hessian diagonally dominant under (4.11); this is a legitimate sufficient-condition proof, not a circular one, because the condition is verified algebraically for the given system rather than assumed. The approximative system (4.1) is motivated by a relaxation technique from [24] and formal limits referenced to [9,12], but the paper itself marks the convergence as expected and open (Remark 4.1 and Section 6), so the approximation claim is explicitly conditional and not disguised as a prediction. Numerical tests are benchmarked against the CH/NSCH limit solution, providing external comparison; even the β=10 shock-tube experiment outside the proved regime (β<1 for the modified potential) is a parameter-domain issue, not a circularity. There is minor self-citation in the derivation of the model form, but it is not load-bearing for the validity of Theorem 4.3.

Assumptions & free parameters 3 free parameters · 4 assumptions · 1 invented entities

The central hyperbolicity theorem is self-contained and relies only on standard assumptions about the free energy and the parameter regime. The approximative nature of the system, however, rests on unproven formal limits for α, δ, β, which are listed as axioms. No data fitting is involved.

free parameters (3)
  • α (artificial compressibility) = varied in experiments (e.g., 10⁻⁴ to 0.5)
    Small parameter scaling the pressure time derivative; not fitted to data, chosen by hand.
  • δ (friction parameter) = varied (10⁻⁸ to 0.0625)
    Controls the friction-type relaxation of v; chosen by hand.
  • β (relaxation parameter) = varied (0.001 to 10)
    Controls the elliptic relaxation of the capillarity term; chosen to satisfy (4.11) in the hyperbolicity theorem, but set to 10 in §5.3.4 outside the theorem's regime.
assumptions (4)
  • domain assumption Quartic double-well free energy W(c) = (c²-1)²/4
    Adopted from [17]; hyperbolicity condition (4.11) is evaluated for this W.
  • ad hoc to paper Formal limit β→0 with β⁻¹(ω-c) → γ∆c
    Used in Remark 4.1 to justify the relaxation; not proven, convergence left open.
  • ad hoc to paper Formal limit δ→0 with v → -∇(W'(c)-γ∆c)
    Used in Remark 3.1 to recover the convective Cahn-Hilliard equation; singular perturbation limit not justified rigorously.
  • domain assumption Artificial compressibility ansatz pt + (1/α)∇·u = 0
    Standard approach [4]; expected to converge to the incompressible solution as α→0.
invented entities (1)
  • Relaxation variable ω
    purpose: Auxiliary scalar field governed by the elliptic equation -γ∆ω + (1/β)ω = (1/β)c, used to replace the third-order capillarity term by a first-order flux plus a constraint.
    ω is a purely mathematical device with no physical interpretation; no independent measurement is proposed, and its only role is the formal relaxation limit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Note on Hyperbolic Relaxation of the Navier-Stokes-Cahn-Hilliard system for incompressible two-phase flow." pith.science (2026). https://pith.science/paper/OE55YHXZ

@misc{pith2026241211904,
  author       = {Pith},
  title        = {Pith review of: A Note on Hyperbolic Relaxation of the Navier-Stokes-Cahn-Hilliard system for incompressible two-phase flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OE55YHXZ}},
  note         = {Machine review of arXiv:2412.11904}
}
read the original abstract

We consider the two-phase dynamics of two incompressible and immiscible fluids. As a mathematical model we rely on the Navier-Stokes-Cahn-Hilliard system that belongs to the class of diffuse-interface models. Solutions of the Navier-Stokes-Cahn-Hilliard system exhibit strong non-local effects due to the velocity divergence constraint and the fourth-order Cahn-Hilliard operator. We suggest a new first-order approximative system for the inviscid sub-system. It relies on the artificial-compressibility ansatz for the Navier-Stokes equations, a friction-type approximation for the Cahn-Hilliard equation and a relaxation of a third-order capillarity term. We show under reasonable assumptions that the first-order operator within the approximative system is hyperbolic; precisely we prove for the spatially one-dimensional case that it is equipped with an entropy-entropy flux pair with convex (mathematical) entropy. For specific states we present a numerical characteristic analysis. Thanks to the hyperbolicity of the system, we can employ all standard numerical methods from the field of hyperbolic conservation laws. We conclude the paper with preliminary numerical results in one spatial dimension.

Figures

Figures reproduced from arXiv: 2412.11904 by the authors.

Figure 2.1
Figure 2.1. Sketch of a diffuse-interface configuration with phase boundary of width γ. the parameter γ tends to zero, one expects that an appropriate sharp-interface model is attained in the limit. Here, we are not interested in the sharp-interface limit and therefore fix the parameter γ throughout the entire note. As a meanwhile well-known diffuse-interface approach, we start from the classical Navier–Stokes–Cahn–Hilliard (NS… view at source ↗
Figure 2.2
Figure 2.2. Graphs of the logarithmic free energy W (red) and the quartic free energy W (blue). To close the NSCH system (2.1), we assume homogeneous Dirichlet boundary condition for u and homoge￾neous Neumann boundary conditions for c and µ(c), which are u = 0, ∇c · n = 0, ∇µ(c) · n = 0 on ∂Ω × (0, T). (2.4) Appropriate initial conditions for u and c have to be added. 2.2. Energy dissipation for the NSCH system (2.1) A major s… view at source ↗
Figure 4.3
Figure 4.3. Color plot of sgn(p(.; Q)) for fixed δ = 0.09 and α ∈ [0.01, 0.001] on the vertical axis [PITH_FULL_IMAGE:figures/full_fig_p012_4_3.png] view at source ↗
Figures from the paper (13 more)
Figure 4.5
Figure 4.5. Figure 4.5: Color plot of sgn(p(.; Q)) for fixed α = 0.01 and δ ∈ [0.01, 0.001] on the vertical axis [PITH_FULL_IMAGE:figures/full_fig_p013_4_5.png]
Figure 4.7
Figure 4.7. Figure 4.7: Color plot of sgn(p(.; Q)) for fixed δ = 0.01 and α ∈ [0.01, 0.001] on the vertical axis [PITH_FULL_IMAGE:figures/full_fig_p013_4_7.png]
Figure 4.9
Figure 4.9. Figure 4.9: Color plot of sgn(p(.; Q)) for fixed δ = 0.01 and α ∈ [0.01, 0.01] on the vertical axis [PITH_FULL_IMAGE:figures/full_fig_p014_4_9.png]
Figure 5.11
Figure 5.11. Figure 5.11: Comparison of the pressure p, the velocity u, the flux v and the phase-field variable c for the compression-like test case and different choices of the relaxation parameters, α ∈ {0.5, 0.1, 0.05} and δ = 0.0625 at t = 0.15. both initial conditions, different choices…
Figure 5.12
Figure 5.12. Figure 5.12: Comparison of the pressure p, the velocity u, the flux v and the phase-field variable c for the expansion-like test case and different choices of the relaxation parameters, α ∈ {0.5, 0.1, 0.05} and δ = 0.0625 at t = 0.15. 5.2. The relaxed friction-type model with fr…
Figure 5.13
Figure 5.13. Figure 5.13: Comparison of the pressure p, the velocity u, the flux v and the phase-field variable c for compression-like test case and different choices of the relaxation parameters, α = 0.05 and δ ∈ {0.0625, 0.01, 10−4 , 10−6 } at t = 0.15. one-dimensional Ostwald ripening. Fi…
Figure 5.14
Figure 5.14. Figure 5.14: Comparison of the phase-field variable c ε , the artificial velocity v ε , the integral L2 difference between the phase-field variable and the order parameter and the evolution of the total energy per unit volume E ε over time, for α = 10−4 and δ = 10−8 , γ = 0.001 …
Figure 5.15
Figure 5.15. Figure 5.15: Comparison of the phase-field variable c ε at different time instances for α = β = 10−4 , γ = 0.001 and δ ∈ {10−2 , 10−4 , 10−6 }. t = 0.1 t = 0.2 t = 0.3 0 0.2 0.4 0.6 0.8 1 −1 −0.5 0 0.5 1 x [−] c, c ε [−] CH α = 10−1 α = 10−2 α = 10−3 0 0.2 0.4 0.6 0.8 1 −1 −0.5 …
Figure 5.16
Figure 5.16. Figure 5.16: Comparison of the phase-field variable c ε at different time instances for α ∈ {10−1 , 10−2 , 10−3 }, β = 10−4 , γ = 0.001 and δ = 10−8 . velocity field. Finally, similar to the previous investigations, a monotonic decrease of the total energy per unit volume can be…
Figure 5.17
Figure 5.17. Figure 5.17: Comparison of the velocity u ε , the artificial velocity v ε at t = 0.2 and the evolution of the total energy per unit volume over time for α ∈ {10−1 , 10−2 , 10−3 }, β = 10−4 , γ = 0.001 and δ = 10−8 . t = 0.1 t = 0.2 t = 0.3 0 0.2 0.4 0.6 0.8 1 −1 −0.5 0 0.5 1 x […
Figure 5.18
Figure 5.18. Figure 5.18: Evolution of the phase field c ε over time for α ∈ {10−1 , 10−2 , 10−3 }, β = 10−4 , γ = 0.001 and δ = 10−8 and the modified Ostwald ripening setup. choice of the velocity, the final profile of the phase-field variable at t = 0.3 should coincide with the one of the …
Figure 5.19
Figure 5.19. Figure 5.19: Evolution of the phase field c ε , the velocity u ε and the pressure p ε at t = 0.15 for α = 0.5, β = 0.1, γ = 10−4 and δ = 0.0625 for the compression-like setup. −1 −0.5 0 0.5 1 1 1.5 2 x [−] c ε [−] −1 −0.5 0 0.5 1 −1 0 1 x [−] u ε [−] −1 −0.5 0 0.5 1 −1 −0.5 0 x …
Figure 5.20
Figure 5.20. Figure 5.20: Evolution of the phase field c ε , the velocity u ε and the pressure p ε at t = 0.15 for α = 0.5, β = 0.1, γ = 10−4 and δ = 0.0625 for the expansion-like setup. 2000 cells and for the time step restriction CF L = 0.9 is chosen. The relaxation parameters and the capi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 27 canonical work pages

  1. [1]

    Abels, H

    H. Abels, H. Garcke, and G. Gr¨ un. Thermodynamically consistent, frame indifferent diffuse interface models for incompressible two-phase flows with different densities. Math. Models Methods Appl. Sci. , 22(3):1150013, 40, 2012

    2012

  2. [2]

    D. M. Anderson, G. B. McFadden, and A. A. Wheeler. Diffuse-interface methods in fluid mechanics. Ann. Rev. Fluid Mech. , 30(1):139–165, 1998

    work page 1998

  3. [3]

    Boyer and C

    F. Boyer and C. Lapuerta. Study of a three component Cahn-Hilliard flow model. M2AN Math. Model. Numer. Anal. , 40(4):653–687, 2006

    work page 2006

  4. [4]

    A. J. Chorin. A numerical method for solving incompressible viscous flow problems. J. Comput. Phys. , 2(1):12–26, 1967

    work page 1967

  5. [5]

    A. J. Chorin. The numerical solution of the Navier-Stokes equations for an incompressible fluid. B. Am. Math. Soc., 73(6):928– 931, 1967

    work page 1967

  6. [6]

    A. J. Chorin. Numerical solution of the Navier-Stokes equations. Math. Comput. , 22(104):745–762, 1968

    work page 1968

  7. [7]

    Dhaouadi and M

    F. Dhaouadi and M. Dumbser. A first order hyperbolic reformulation of the Navier-Stokes-Korteweg system based on the GPR model and an augmented Lagrangian approach. J. Comput. Phys. , 470:Paper No. 111544, 30, 2022

    work page 2022

  8. [8]

    A first-order hyperbolic reformulation of the Cahn-Hilliard equation

    F. Dhaouadi, M. Dumbser, and S. Gavrilyuk. A first-order hyperbolic reformulation of the Cahn-Hilliard equation. ArXiv preprint: https://arxiv.org/abs/2408.03862, 2024

    work page Pith review arXiv 2024

Show all 28 references
  1. [9]

    Engel, A

    P. Engel, A. Viorel, and C. Rohde. A low-order approximation for viscous-capillary phase transition dynamics. Port. Math. , 70(4):319–344, 2013

    2013

  2. [10]

    K. O. Friedrichs and P. D. Lax. Systems of conservation equations with a convex extension. Proc. Nat. Acad. Sci. U.S.A. , 68:1686–1688, 1971

    1971

  3. [11]

    Gatti, M

    S. Gatti, M. Grasselli, V. Pata, and A. Miranville. Hyperbolic relaxation of the viscous Cahn-Hilliard equation in 3-D. Math. Models Methods Appl. Sci. , 15(2):165–198, 2005

    2005

  4. [12]

    Giesselmann

    J. Giesselmann. A relative entropy approach to convergence of a low order approximation to a nonlinear elasticity model with viscosity and capillarity. SIAM J. Math. Anal. , 46(5):3518–3539, 2014

    2014

  5. [13]

    Giorgini and R

    A. Giorgini and R. Temam. Weak and strong solutions to the nonhomogeneous incompressible Navier-Stokes-Cahn-Hilliard system. J. Math. Pures Appl. , 144:194–249, 2020

    2020

  6. [14]

    S. K. Godunov. An interesting class of quasi-linear systems. Dokl. Akad. Nauk SSSR , 139:521–523, 1961

    1961

  7. [15]

    P. M. Gresho. Some current CFD issues relevant to the incompressible Navier-Stokes equations. Comput. Methods Appl. Mech. Eng., 87(2-3):201–252, 1991

    1991

  8. [16]

    T. Hitz, J. Keim, C.-D. Munz, and C. Rohde. A parabolic relaxation model for the Navier-Stokes-Korteweg equations. J. Comput. Phys. , 421:109714, 2020

    2020

  9. [17]

    P. C. Hohenberg and B. I. Halperin. Theory of dynamic critical phenomena. Rev. Mod. Phys. , 49(3):435, 1977

    1977

  10. [18]

    Huang, C

    Q. Huang, C. Rohde, W.-A. Yong, and R. Zhang. A hyperbolic relaxation system of the incompressible navier-stokes equations with artificial compressibility. ArXiv preprint: https://arxiv.org/abs/2411.15575, 2024

    2024 arXiv

  11. [19]

    Keim, H.-C

    J. Keim, H.-C. Konan, and C. Rohde. Robust time-discretization methods for hyperbolic approximations of the Navier-Stokes- Cahn-Hilliard system. In preparation

  12. [20]

    Keim, C.-D

    J. Keim, C.-D. Munz, and C. Rohde. A relaxation model for the non-isothermal Navier-Stokes-Korteweg equations in confined domains. J. Comput. Phys. , 474:Paper No. 111830, 28, 2023

    2023

  13. [21]

    Lattanzio and A

    C. Lattanzio and A. E. Tzavaras. From gas dynamics with large friction to gradient flows describing diffusion theories.Commun. Partial Differ. Equ.. , 42(2):261–290, 2017

    2017

  14. [22]

    Massa, L

    F. Massa, L. Ostrowski, F. Bassi, and C. Rohde. An artificial equation of state based Riemann solver for a discontinuous Galerkin discretization of the incompressible Navier-Stokes equations. J. Comput. Phys. , 448:Paper No. 110705, 25, 2022. 23

    2022

  15. [23]

    Nithiarasu

    P. Nithiarasu. An efficient artificial compressibility (AC) scheme based on the characteristic based split (CBS) method for incompressible flows. Int. J. Numer. Methods Eng. , 56(13):1815–1845, 2003

    2003

  16. [24]

    C. Rohde. A local and low-order Navier-Stokes-Korteweg system. In Nonlinear partial differential equations and hyperbolic wave phenomena , volume 526 of Contemp. Math. , pages 315–337. Amer. Math. Soc., Providence, RI, 2010

    2010

  17. [25]

    Rohde and L

    C. Rohde and L. von Wolff. Homogenization of nonlocal Navier-Stokes-Korteweg equations for compressible liquid-vapor flow in porous media. SIAM J. Math. Anal. , 52(6):6155–6179, 2020

    2020

  18. [26]

    Rohde and L

    C. Rohde and L. von Wolff. A ternary Cahn-Hilliard–Navier-Stokes model for two-phase flow with precipitation and dissolution. Math. Models Methods Appl. Sci. , 31(1), 2021

    2021

  19. [27]

    R. Temam. Sur l’approximation de la solution des ´ equations de Navier-Stokes par la m´ ethode des pas fractionnaires (I).Arch. Ration. Mech. Anal., 32:135–153, 1969

    1969

  20. [28]

    R. Temam. Sur l’approximation de la solution des ´ equations de Navier-Stokes par la m´ ethode des pas fractionnaires (II).Arch. Ration. Mech. Anal., 33:377–385, 1969. 24

    1969

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.