Projection-Based NSCH Frameworks

• Projection-Based NSCH frameworks are advanced numerical methods that employ projection operators to decouple complex multiphysics systems like multiphase flows and nonsmooth mechanics.
• They use both Eulerian and Lagrangian discretizations with time-splitting techniques to achieve strict mass conservation and energy stability.
• Benchmark tests demonstrate that these methods ensure accurate enforcement of physical constraints and robust handling of interface dynamics.

Projection-based NSCH (Navier–Stokes–Cahn–Hilliard) frameworks are an advanced class of numerical methods that employ projection operators for constraint enforcement, decoupling, or stabilization in the simulation of multiphase flows, non-smooth mechanics, and related systems governed by the Navier–Stokes and Cahn–Hilliard equations. These frameworks are central to the development of energy-dissipative, mass-conservative, and stable time-integration schemes in both Eulerian and Lagrangian settings. The principal idea is to utilize projections—often oblique, sometimes orthogonal—to split complex coupled systems into tractable sub-problems subject to physical constraints such as incompressibility or nonsmooth dynamics.

1. Core Principles and Mathematical Foundation

Projection-based NSCH methods rely on decomposition strategies, in which the evolution equations (NS, CH, or joint systems) are advanced by time-splitting coupled with projection steps. In multiphase fluid systems, the canonical projection step enforces the incompressibility condition or mass constraint after prediction of an intermediate velocity. In non-smooth mechanics, oblique projectors enforce constraint satisfaction or impact laws across discontinuities.

The construction of these projection operators is highly problem dependent:

• In incompressible multiphase flow, the projection typically matches the Helmholtz–Hodge decomposition, projecting intermediate velocities onto the divergence-free subspace via the solution of a pressure Poisson equation.
• For non-smooth mechanical systems, the projection is onto the tangent space of the constraint manifold, formulated as a closed-form oblique projector:

$$\mathcal P = I - M^{-1} J^T (J M^{-1} J^T)^\dagger J$$

where $M$ is the symmetric positive-definite inertia matrix, $J$ the constraint Jacobian, and $(\cdot)^\dagger$ the Moore–Penrose pseudoinverse (Aghili, 2021).

These projections are constructed to guarantee compatibility with the system's algebraic constraints, ensuring that discrete solutions satisfy incompressibility or contact conditions as required.

2. Discretizations: Eulerian, Lagrangian, and Mesh-Free Approaches

Projection-based NSCH frameworks have emerged in both classical grid-based (Eulerian) and particle-based (Lagrangian) settings.

• Finite Difference Projection-Type Schemes: For staggered-grid Eulerian discretizations of the NSCH system, projection steps decouple the update of velocity and pressure, ensuring both mass conservation and energy stability at the fully discrete level. Scheme variants handle quasi-incompressibility or variable density by incorporating additional projection terms and corrections (Guo et al., 2017).
• Smoothed Particle Hydrodynamics (SPH) Frameworks: In the Lagrangian, mesh-free setting, as in energy-stable SPH-NSCH, anti-symmetric difference operators are used for discretizing gradients and Laplacians. The projection step enforces discrete divergence-free velocity via a Poisson solve, inheriting exact mass, momentum, and energy dissipation properties even for large deformations (Feng et al., 2022).
• Non-smooth Multibody Dynamics: The use of linear projection operators, particularly oblique projectors, enables closed-form updates across discontinuities (e.g., impacts or topology change), providing a geometric foundation for hybrid integration schemes and guaranteeing energetic consistency under constraints (Aghili, 2021).

3. Time-Splitting and Projection Workflow

Projection-based NSCH schemes proceed in sequential operator-split steps, typically structured as follows:

Step	Description	Key Equations (Example)
Phase-field update	Solve Cahn–Hilliard for phase variable	See $d\phi_i/dt$ in (Feng et al., 2022)
Velocity predictor	Advance momentum without enforcing pressure constraint	$\tilde{\mathbf v}_i^{k+1} = ...$
Projection/Poisson	Enforce incompressibility or constraint via pressure update	$\nabla \cdot \mathbf v = 0$
Particle/grid advect	Update positions/cells using latest velocity	$\mathbf x_i^{k+1} = ...$

At each time level, these schemes enforce discrete energy laws and mass conservation, with rigorous proofs that the projection/correction steps do not increase kinetic or total energy (Guo et al., 2017, Feng et al., 2022).

4. Conservation, Stability, and Energy Laws

A defining characteristic of projection-based NSCH frameworks is their ability to enforce physical conservation laws and monotonic energy dissipation at the discrete level:

• Mass conservation: Achieved via symmetric discretization and flux cancellation in both Eulerian and Lagrangian approaches, validating $\sum \rho c$ (or $\sum m_i \phi_i$) is preserved at machine precision under periodic/no-flux boundaries.
• Momentum conservation: Guaranteed through pairwise symmetry in dissipation and pressure projection steps, in particular for mesh-free methods (Feng et al., 2022).
• Energy stability: Discrete analogues of the continuum free-energy law are established by telescoping summations and energetic cancellation, leading to

$M$0

regardless of time step, with dissipation terms corresponding to viscosity, mobility, and numerical stabilization (Guo et al., 2017, Feng et al., 2022).

• Energetic consistency in nonsmooth mechanics: For projections enforcing constraints post-impact, kinetic energy is guaranteed not to increase, with analytically bounded energy loss ratios according to restitution coefficients (Aghili, 2021).

5. Algorithmic Realizations and Computational Strategies

Projection-based NSCH frameworks employ a combination of sparse linear (or nonlinear) solves, multigrid accelerations, and stabilization mechanisms:

• SPH projection-based NSCH: Requires the solution of SPD systems for the CH and Poisson problems at each step, amenable to mesh-free multigrid or AMG (algebraic multigrid) preconditioners (Feng et al., 2022).
• Finite difference schemes: Multigrid solvers with Vanka-type smoothers on staggered grids handle the coupled CH and NS sub-systems efficiently, exploiting local block structure for robust convergence (Guo et al., 2017).
• Nonsmooth multibody mechanics: Projection formulas admit closed-form implementation, with conditioning of the constraint-inertia matrix optimized by regularization. Condition number minimization is achieved by selecting the regularization parameter within the spectrum of projected inertia (Aghili, 2021).

Numerical tests across frameworks (e.g., droplet deformation, Rayleigh–Taylor, capillary waves) demonstrate exact constraint satisfaction, monotone energy decay, and robustness to strong interface dynamics.

6. Extensions: Multiphysics, Nonsmooth Events, and Energetic Constraints

Projection-based NSCH approaches extend naturally to systems with variable density, viscosity, or nontrivial constraint architectures:

• Variable-parameter NSCH systems: Projection-based schemes accommodate density and viscosity dependence in both convective and diffusive terms. Mass and energy conservation are rigorously maintained even for large density ratios (Guo et al., 2017).
• Locally varying restitution in impacts: Oblique projectors can be parameterized to handle codependent restitution, with local coefficients constrained by LMIs to guarantee energetic admissibility (Aghili, 2021). The set of all restitution matrices that avoid artificial energy creation is characterized exactly.
• Coupling to other multiphysics: The framework is compatible with additional forces, topological change, and phase transitions by extending the basis for projection to new constraints or conservation laws.

7. Representative Results and Performance Considerations

Projection-based NSCH methods show favorable computational properties:

• Time-stepping stability: Unconditional energy-stability allows larger $M$1 than classical weakly-compressible or explicit methods. Practical time steps up to $M$2 are reported (Feng et al., 2022).
• Accuracy and benchmarking: Quantitative validation includes mass conservation to $M$3 or better, accurate interface and velocity dynamics, and robust handling of vortical and merging-droplet phenomena against analytical and benchmark solutions (Guo et al., 2017).
• Efficiency: The dominant cost lies in linear/quasi-linear solv es, with high scalability on modern architectures using multigrid or local block smoothers.

A plausible implication is that the integration of projection operators across discretization paradigms, from mesh-based, mesh-free, to hybrid or adaptive formulations, will continue to support the development of robust, physically faithful simulation capabilities for complex multiphysics and nonsmooth systems.

---

References:

• Energy-stable SPH-NSCH: (Feng et al., 2022)
• Nonsmooth mechanics with linear projection: (Aghili, 2021)
• Mass conservative, energy stable projection NSCH: (Guo et al., 2017)