skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • On a SAV-MAC scheme for the Cahn–Hilliard–Navier–Stokes phase-field model and its error analysis for the corresponding Cahn–Hilliard–Stokes case
Title: On a SAV-MAC scheme for the Cahn–Hilliard–Navier–Stokes phase-field model and its error analysis for the corresponding Cahn–Hilliard–Stokes case
We construct a numerical scheme based on the scalar auxiliary variable (SAV) approach in time and the MAC discretization in space for the Cahn–Hilliard–Navier–Stokes phase- field model, prove its energy stability, and carry out error analysis for the corresponding Cahn–Hilliard–Stokes model only. The scheme is linear, second-order, unconditionally energy stable and can be implemented very efficiently. We establish second-order error estimates both in time and space for phase-field variable, chemical potential, velocity and pressure in different discrete norms for the Cahn–Hilliard–Stokes phase-field model. We also provide numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy of our scheme.  more » « less
Award ID(s):
2012585
PAR ID:
10225298
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Mathematical Models and Methods in Applied Sciences
Volume:
30
Issue:
12
ISSN:
0218-2025
Page Range / eLocation ID:
2263 to 2297
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Journal of Computational Physics)
    In this paper we construct a novel discretization of the Cahn-Hilliard equation coupled with the Navier-Stokes equations. The Cahn-Hilliard equation models the separation of a binary mixture. We construct a very simple time integration scheme for simulating the Cahn-Hilliard equation, which is based on splitting the fourth-order equation into two second-order Helmholtz equations. We combine the Cahn-Hilliard equation with the Navier-Stokes equations to simulate phase separation in a two-phase fluid flow in two dimensions. The scheme conserves mass and momentum and exhibits consistency between mass and momentum, allowing it to be used with large density ratios. We introduce a novel discretization of the surface tension force from the phase-field variable that has finite support around the transition region. The model has a parameter that allows it to transition from a smoothed continuum surface force to a fully sharp interface formulation. We show that our method achieves second-order accuracy, and we compare our method to previous work in a variety of experiments. 
    more » « less
  2. ; ; ; (, International Journal for Numerical Methods in Engineering)
    ABSTRACT This paper is concerned with efficient and accurate numerical schemes for the Cahn‐Hilliard‐Navier‐Stokes phase field model of binary immiscible fluids. By introducing two Lagrange multipliers for each of the Cahn‐Hilliard and Navier‐Stokes parts, we reformulate the original model problem into an equivalent system that incorporates the energy evolution process. Such a nonlinear, coupled system is then discretized in time using first‐ and second‐order backward differentiation formulas, in which all nonlinear terms are treated explicitly and no extra stabilization term is imposed. The proposed dynamically regularized Lagrange multiplier (DRLM) schemes are mass‐conserving and unconditionally energy‐stable with respect to the original variables. In addition, the schemes are fully decoupled: Each time step involves solving two biharmonic‐type equations and two generalized linear Stokes systems, together with two nonlinear algebraic equations for the Lagrange multipliers. A key feature of the DRLM schemes is the introduction of the regularization parameters which ensure the unique determination of the Lagrange multipliers and mitigate the time step size constraint without affecting the accuracy of the numerical solution, especially when the interfacial width is small. Various numerical experiments are presented to illustrate the accuracy and robustness of the proposed DRLM schemes in terms of convergence, mass conservation, and energy stability. 
    more » « less
  3. ; ; ; ; ; ; (, Computer physics communications)
    We present a fully-coupled, implicit-in-time framework for solving a thermodynamically-consistent Cahn-Hilliard Navier-Stokes system that models two-phase flows. In this work, we extend the block iterative method presented in Khanwale et al. [Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes, J. Comput. Phys. (2020)], to a fully-coupled, provably second-order accurate scheme in time, while maintaining energy-stability. The new method requires fewer matrix assemblies in each Newton iteration resulting in faster solution time. The method is based on a fully-implicit Crank-Nicolson scheme in time and a pressure stabilization for an equal order Galerkin formulation. That is, we use a conforming continuous Galerkin (cG) finite element method in space equipped with a residual-based variational multiscale (RBVMS) procedure to stabilize the pressure. We deploy this approach on a massively parallel numerical implementation using parallel octree-based adaptive meshes. We present comprehensive numerical experiments showing detailed comparisons with results from the literature for canonical cases, including the single bubble rise, Rayleigh-Taylor instability, and lid-driven cavity flow problems. We analyze in detail the scaling of our numerical implementation. 
    more » « less
  4. ; ; ; (, SIAM journal on scientific computing)
    In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions. 
    more » « less
  5. ; ; ; (, Mathematics of Computation)
    A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment for the singular logarithmic terms and the surface diffusion term, and an explicit update for the expansive concave term. The convective term for the phase variable, as well as the coupled term in the Stokes equation, is approximated in a semi-implicit manner. In the spatial discretization, the marker and cell difference method is applied, which evaluates the velocity components, the pressure and the phase variable at different cell locations. Such an approach ensures the divergence-free feature of the discrete velocity, and this property plays an important role in the analysis. The positivity-preserving property and the unique solvability of the proposed numerical scheme are theoretically justified, utilizing the singular nature of the logarithmic term as the phase variable approaches the singular limit values. An unconditional energy stability analysis is standard, as an outcome of the convex-concave decomposition technique. A convergence analysis with accompanying error estimate is provided for the proposed numerical scheme. In particular, a higher order consistency analysis, accomplished by supplementary functions, is performed to ensure the separation properties of numerical solution. In turn, using the approach of rough and refined error estimates, we are able to derive an optimal rate convergence. To conclude, several numerical experiments are presented to validate the theoretical analysis. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email