Simulation and Control of a Nonsmooth Cahn-Hilliard Navier-Stokes System
Pith reviewed 2026-05-25 00:14 UTC · model grok-4.3
The pith
Optimal solutions exist for controlling a nonsmooth two-phase flow model, with C- and strong stationarity conditions derived for the bilevel problem.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We establish the existence of optimal solutions and present two distinct approaches to derive suitable stationarity conditions for the bilevel problem, namely C- and strong stationarity. The numerical realization relies on two adaptive solution algorithms that use a specifically developed goal-oriented error estimator. A model order reduction approach using proper orthogonal decomposition replaces high-fidelity models by low-order surrogates, combining POD with space-adapted snapshots to handle different spatial resolutions while conserving the solenoidal property.
What carries the argument
Bilevel optimal control problem for the nonsmooth Cahn-Hilliard Navier-Stokes system, with derivation of C-stationarity and strong stationarity conditions as the core mechanism for characterizing optimality.
If this is right
- Existence of optimal solutions makes the bilevel control problem well-posed for two-phase flows.
- C- and strong stationarity conditions supply necessary optimality criteria that support numerical solvers.
- Goal-oriented adaptive algorithms become applicable for accurate simulation and control.
- POD-MOR with space-adapted snapshots yields low-order surrogates that maintain the divergence-free velocity constraint.
Where Pith is reading between the lines
- The stationarity approaches might extend to other coupled nonsmooth fluid models with similar bilevel structure.
- The handling of varying snapshot resolutions could inform adaptive reduced-order methods in related multiphase problems.
- Numerical tests could check whether strong stationarity yields tighter control bounds than C-stationarity in practice.
Load-bearing premise
The nonsmooth energy potential and the coupled PDE system are assumed to be sufficiently regular to admit optimal solutions and to permit the derivation of C- and strong stationarity conditions.
What would settle it
A concrete example of a nonsmooth Cahn-Hilliard Navier-Stokes system where no optimal control exists, or where computed solutions violate the derived C- or strong stationarity conditions.
Figures
read the original abstract
We are concerned with the simulation and control of a two phase flow model governed by a coupled Cahn-Hilliard Navier-Stokes system involving a nonsmooth energy potential. We establish the existence of optimal solutions and present two distinct approaches to derive suitable stationarity conditions for the bilevel problem, namely C- and strong stationarity. Moreover, we demonstrate the numerical realization of these concepts at the hands of two adaptive solution algorithms relying on a specifically developed goal-oriented error estimator. In addition, we present a model order reduction approach using proper orthogonal decomposition (POD-MOR) in order to replace high-fidelity models by low order surrogates. In particular, we combine POD with space-adapted snapshots and address the challenges which are the consideration of snapshots with different spatial resolutions and the conservation of a solenoidal property.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the simulation and optimal control of a coupled Cahn-Hilliard Navier-Stokes system with a nonsmooth energy potential modeling two-phase flows. It proves existence of optimal solutions to the associated bilevel control problem and derives C-stationarity and strong stationarity conditions via two distinct approaches. Numerically, the work develops two adaptive solution algorithms that employ a goal-oriented error estimator and introduces a POD-based model-order-reduction technique that accommodates snapshots of differing spatial resolutions while preserving the solenoidal constraint.
Significance. If the existence and stationarity results are valid, the paper supplies rigorous first-order conditions for a class of nonsmooth, coupled PDE-constrained optimization problems that are otherwise difficult to treat. The combination of limiting subdifferential calculus with regularization arguments and the practical demonstration of adaptive and reduced-order methods constitute a concrete advance for control of two-phase flow models.
minor comments (2)
- [Introduction] The abstract states that two distinct approaches are used to obtain C- and strong stationarity, yet the introduction does not explicitly contrast the two routes (e.g., which assumptions each route relaxes). Adding a short comparative paragraph would improve readability.
- [Section on POD-MOR] In the POD-MOR section the preservation of the divergence-free property after projection onto the reduced basis is asserted, but the precise construction of the solenoidal reduced space (e.g., via a discrete Helmholtz decomposition or constrained POD) is not detailed; a short algorithmic box or equation would clarify the implementation.
Simulated Author's Rebuttal
We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. The report correctly identifies the core contributions on existence of optimal solutions, C- and strong stationarity conditions for the bilevel problem, goal-oriented adaptive algorithms, and the POD-MOR approach handling space-adapted snapshots while preserving the divergence-free constraint.
Circularity Check
No significant circularity
full rationale
The paper establishes existence of optimal solutions for the bilevel control problem and derives C- and strong stationarity conditions via standard limiting subdifferential arguments, regularization, and function-space compactness results for the nonsmooth Cahn-Hilliard-Navier-Stokes system. No load-bearing step reduces by definition or self-citation to a fitted input or prior result whose validity depends on the present work; the stationarity systems are obtained from external PDE theory rather than internal renaming or ansatz smuggling. The derivation chain is therefore self-contained.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The nonsmooth Cahn-Hilliard Navier-Stokes system admits solutions that allow existence of optimal controls and derivation of stationarity conditions.
Forward citations
Cited by 1 Pith paper
-
An optimal control problem for Stokes-Cahn-Hilliard-Oono equations with regular potential
Existence of an optimal control is established for the Stokes-Cahn-Hilliard-Oono equations, and first-order optimality conditions are derived via the corresponding adjoint system.
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
Aland, S. and Voigt, A. (2012). Benchmark computations of diffuse interface models for two-dimensional bubble dynamics. International Journal for Numerical Methods in Fluids , 69:747–761
work page 2012
-
[4]
J. O. Alff, Modellordnungsreduktion f¨ ur das Cahn-Hilliard System, Bachelorarbeit, Univer- sit¨ at Hamburg (2015)
work page 2015
- [5]
- [6]
-
[7]
A. Alla, C. Gr¨ aßle, M. Hinze, A-posteriori snapshot location for POD in optimal control of linear parabolic equations, ESAIM: M2AN 52(5) (2018), 1847–1873
work page 2018
-
[8]
A. Alla, N. Kutz, Nonlinear model order reduction via dynamic mode decomposition , SIAM J. Sci. Comput., 39(5) (2017), B778–B796
work page 2017
-
[9]
J. O. Alff, Trust Region POD for Optimal Control of Cahn-Hilliard Systems, Master’s Thesis, Universit¨ at Hamburg (2018)
work page 2018
-
[10]
Barbu, Optimal control of variational inequalities , vol
V. Barbu, Optimal control of variational inequalities , vol. 100 of Research Notes in Mathe- matics, Pitman (Advanced Publishing Program), Boston, MA, 1984
work page 1984
-
[11]
T. Borrvall and J. Petersson. Topology optimization of fluids in Stokes flow . Internat. J. Numer. Methods Fluids 41(1) (2003), 77–107
work page 2003
-
[12]
F. Boyer, A theoretical and numerical model for the study of incompressible mixture flows , Computers & fluids, 31 (2002), pp. 41–68
work page 2002
- [13]
-
[14]
Brett, C., Elliott, C.M., Hinterm¨ uller, M., L¨ obhard, C.: Mesh adaptivity in optimal control of elliptic variational inequalities with point-tracking of the state. Interfaces Free Bound. 17(1), 21–53 (2015)
work page 2015
- [15]
-
[16]
O. Burkovska, B. Haasdonk, J. Salomon, B. Wohlmuth, Reduced Basis Methods for Pricing Options with the Black-Scholes and Heston Models , SIAM J. Financial Math., 6(1) (2015), 685–712
work page 2015
-
[17]
M. Barrault, Y. Maday, N. C. Nguyen, A. T. Patera, An “empirical interpolation” method: application to efficient reduced-basis discretization of partial differential equations , C. R. Acad. Sci. Paris, 339(9) (2004), 667–672
work page 2004
-
[18]
F. Ballarin, A. Manzoni, A. Quarteroni, G. Rozza, Supremizer stabilization of POD-Galerkin approximation of parametrized steady incompressible Navier-Stokes equations, Int. J. Numer. Methods Eng. 102(5) (2015), 1136–1161
work page 2015
-
[19]
P. Colli, M. H. Farshbaf-Shaker, G. Gilardi, and J. Sprekels, Optimal boundary control of a viscous Cahn-Hilliard system with dynamic boundary condition and double obstacle potentials, SIAM J. Control Optim., 53 (2015), pp. 2696–2721. SIMULATION AND CONTROL OF A NONSMOOTH CHNS SYSTEM 23
work page 2015
-
[20]
N. Cagniart, Y. Maday, B. Stamm, Model Order Reduction with large Convection Effects , Contributions to Partial Differential Equations (2018), 131–150
work page 2018
-
[21]
S. Chaturantabut, D. C. Sorensen, Nonlinear model order reduction via discrete empirical interpolation, SIAM J. Sci. Comput., 32(5) (2010), 2737–2764
work page 2010
-
[22]
H. Ding, P. DM Spelt, and C. Shu , Diffuse interface model for incompressible two-phase flows with large density ratios, Journal of Computational Physics, 226 (2007), pp. 2078–2095
work page 2007
-
[23]
S. Eckert, P. A. Nikrityuk, B. Willers, D. R ¨abiger, N. Shevchenko, H. Neumann- Heyme, V. Travnikov, S. Odenbach, A. Voigt, and K. Eckert , Electromagnetic melt flow control during solidification of metallic alloys , The European Physical Journal Special Topics, 220 (2013), pp. 123–137
work page 2013
-
[24]
C. M. Elliott, A. Stuart, The global dynamics of discrete semilinear parabolic euqations , SIAM J. Numer. Anal., 30(6) (1993), 1622–1663
work page 1993
-
[25]
D. J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation , MRS Online Proceedings Library Archive 529 (1998)
work page 1998
-
[26]
C. M. Elliott and Z. Songmu, On the Cahn-Hilliard equation, Arch. Rational Mech. Anal., 96 (1986), pp. 339–357
work page 1986
-
[27]
S. Frigeri, E. Rocca, and J. Sprekels , Optimal distributed control of a nonlocal Cahn- Hilliard/Navier-Stokes system in two dimensions , SIAM J. Control Optim., 54 (2016), pp. 221–250
work page 2016
-
[28]
C. G. Gal and M. Grasselli, Asymptotic behavior of a Cahn-Hilliard-Navier-Stokes system in 2D, Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire, 27 (2010), pp. 401–436
work page 2010
-
[29]
C. Gr¨ aßle, M. Hinze,POD reduced-order modeling for evolution equations utilizing arbitrary finite element discretizations , Adv. Comput. Math. 44(6) (2018), 1941–1978
work page 2018
- [30]
- [31]
-
[32]
C. Gr¨ aßle, M. Hinze,The combination of POD model reduction with adaptive finite element methods in the context of phase field models , PAMM, 17(1) (2017), 47–50
work page 2017
-
[33]
C. Gr¨ aßle, M. Hinze, J. Lang, S. Ullmann, POD model order reduction with space-adapted snapshots for incompressible flows , accepted for publication in Adv. Comput. Math. (2018), preprint available https://arxiv.org/abs/1810.03892
-
[34]
C. Gr¨ aßle, M. Hinze, N. Scharmacher,POD for optimal control of the Cahn-Hilliard system using spatially adapted snapshots , in Numerical Mathematics and Advanced Applications ENUMATH 2017 (2019), 703–711
work page 2017
- [35]
-
[36]
M. Hinterm¨ uller, M. Hinze, C. Kahle,An adaptive finite element Moreau-Yosida-based solver for a coupled Cahn-Hilliard/Navier-Stokes system , J. Comput. Phys., 235 (2013), 810–827
work page 2013
-
[37]
M. Hinterm¨ uller, M. Hinze, C. Kahle, T. Keil, A goal-oriented dual-weighted adaptive fi- nite element approach for the optimal control of a nonsmooth Cahn-Hilliard Navier-Stokes system, Optim. Eng. 19(3) (2018), 629–662
work page 2018
-
[38]
M. Hinterm¨ uller, M. Hinze, M. Tber,An adaptive finite element Moreau-Yosida-based solver for a nonsmooth Cahn-Hilliard problem , Optim. Method. Softw., 25 (2011), 777-811
work page 2011
-
[39]
Mathematical Programming 160(1), 271–305 (2016)
Hinterm¨ uller, M., Surowiec, T.: A bundle-free implicit programming approach for a class of elliptic MPECs in function space. Mathematical Programming 160(1), 271–305 (2016)
work page 2016
-
[40]
Hinterm¨ uller, M., Keil, T., Wegner, D.: Optimal control of a semidiscrete Cahn-Hilliard- Navier-Stokes system with nonmatched fluid densities. SIAM J. Control Optim. 55(3), 1954– 1989 (2017)
work page 1954
-
[41]
M. Hinterm ¨uller and I. Kopacka , Mathematical programs with complementarity con- straints in function space: C- and strong stationarity and a path-following algorithm , SIAM J. Optim., 20 (2009), pp. 868–902
work page 2009
-
[42]
M. Hinterm¨uller, B. S. Mordukhovich, and T. M. Surowiec, Several approaches for the derivation of stationarity conditions for elliptic MPECs with upper-level control constraints , Math. Program., 146 (2014), pp. 555–582. 24 GR ¨ASSLE, HINTERM ¨ULLER, HINZE, AND KEIL
work page 2014
-
[43]
M. Hinterm¨uller and D. Wegner, Distributed optimal control of the Cahn-Hilliard system including the case of a double-obstacle homogeneous free energy density , SIAM J. Control Optim., 50 (2012), pp. 388–418
work page 2012
-
[44]
M. Hinterm ¨uller and D. Wegner , Distributed and boundary control problems for the semidiscrete Cahn-Hilliard/Navier-Stokes system with nonsmooth Ginzburg-Landau ener- gies, Isaac Newton Institute preprint:NI14042-FRB, (2014)
work page 2014
-
[45]
M. Hinterm¨uller and D. Wegner, Optimal control of a semidiscrete Cahn-Hilliard-Navier- Stokes system, SIAM J. Control Optim., 52 (2014), pp. 747–772
work page 2014
-
[46]
P. C. Hohenberg and B. I. Halperin , Theory of dynamic critical phenomena , Reviews of Modern Physics, 49 (1977), p. 435
work page 1977
-
[47]
Hysing, S., Turek, S., Kuzmin, D., Parolini, N., Burman, E., Ganesan, S., and Tobiska, L. (2009). Quantitative benchmark computations of two-dimensional bubble dynamics. Inter- national Journal for Numerical Methods in Fluids , 60(11):1259–1288
work page 2009
-
[48]
J. Jaruˇsek, M. Krbec, M. Rao, and J. Soko lowski, Conical differentiability for evolution variational inequalities, J. Differential Equations, 193 (2003), pp. 131–146
work page 2003
-
[49]
J. Kim, K. Kang, and J. Lowengrub , Conservative multigrid methods for Cahn-Hilliard fluids, J. Comput. Phys., 193 (2004), pp. 511–543
work page 2004
- [50]
-
[51]
K. Kunisch, S. Volkwein, Galerkin proper orthogonal decomposition methods for a general equation in fluid dynamics , SIAM J. Numer. Anal., 40(2) (2002), 492–515
work page 2002
-
[52]
J. Lowengrub and L. Truskinovsky, Quasi-incompressible Cahn-Hilliard fluids and topo- logical transitions, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci., 454 (1998), pp. 2617– 2654
work page 1998
- [53]
-
[54]
Mordukhovich, Variational analysis and generalized differentiation
Boris S. Mordukhovich, Variational analysis and generalized differentiation. II, vol. 331 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences], Springer-Verlag, Berlin, 2006. Applications
work page 2006
-
[55]
J. Outrata, M. Ko ˇcvara, and J. Zowe , Nonsmooth approach to optimization problems with equilibrium constraints, vol. 28 of Nonconvex Optimization and its Applications, Kluwer Academic Publishers, Dordrecht, 1998. Theory, applications and numerical results
work page 1998
-
[56]
A Navier Stokes Phase Field Crystal Model for Colloidal Suspensions
S. Praetorius and A. Voigt , A phase field crystal model for colloidal suspensions with hydrodynamic interactions, arXiv preprint arXiv:1310.5495, (2013)
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [57]
- [58]
-
[59]
Sirovich, Turbulence and the dynamics of coherent structures I-III , Q
L. Sirovich, Turbulence and the dynamics of coherent structures I-III , Q. Appl. Math., 45(3) (1987), 561–590
work page 1987
- [60]
-
[61]
H. Scheel and S. Scholtes, Mathematical programs with complementarity constraints: sta- tionarity, optimality, and sensitivity , Math. Oper. Res., 25 (2000), pp. 1–22
work page 2000
-
[62]
Tachim Medjo, T.: Optimal control of a Cahn-Hilliard-Navier-Stokes model with state con- straints. J. Convex Anal. 22(4), 1135–1172 (2015)
work page 2015
- [63]
-
[64]
S. Ullmann, M. Rotkvic, J. Lang, POD-Galerkin reduced-order modeling with adaptive finite element snapshots, J. Comput. Phys. 325 (2016), 244–258
work page 2016
- [65]
-
[66]
J. M. Yong and S. M. Zheng , Feedback stabilization and optimal control for the Cahn- Hilliard equation, Nonlinear Anal., 17 (1991), pp. 431–444
work page 1991
-
[67]
B. Zhou et al. , Simulations of polymeric membrane formation in 2D and 3D , PhD thesis, Massachusetts Institute of Technology, 2006. SIMULATION AND CONTROL OF A NONSMOOTH CHNS SYSTEM 25
work page 2006
-
[68]
J. Zowe and S. Kurcyusz, Regularity and stability for the mathematical programming prob- lem in Banach spaces , Appl. Math. Optim., 5 (1979), pp. 49–62. University of Hamburg, Bundesstr. 55, 20146 Hamburg, Germany E-mail address: carmen.graessle@uni-hamburg.de Weierstraß-Institut, Mohrenstr. 39, 10117 Berlin, Germany E-mail address: michael.hintermueller@...
work page 1979
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.