skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, February 13 until 2:00 AM ET on Friday, February 14 due to maintenance. We apologize for the inconvenience.


This content will become publicly available on July 1, 2025

Title: Fully decoupled unconditionally stable Crank–Nicolson leapfrog numerical methods for the Cahn–Hilliard–Darcy system
Abstract

We develop two totally decoupled, linear and second‐order accurate numerical methods that are unconditionally energy stable for solving the Cahn–Hilliard–Darcy equations for two phase flows in porous media or in a Hele‐Shaw cell. The implicit‐explicit Crank–Nicolson leapfrog method is employed for the discretization of the Cahn–Hiliard equation to obtain linear schemes. Furthermore the artificial compression technique and pressure correction methods are utilized, respectively, so that the Cahn–Hiliard equation and the update of the Darcy pressure can be solved independently. We establish unconditionally long time stability of the schemes. Ample numerical experiments are performed to demonstrate the accuracy and robustness of the numerical methods, including simulations of the Rayleigh–Taylor instability, the Saffman–Taylor instability (fingering phenomenon).

 
more » « less
Award ID(s):
2310340
PAR ID:
10521008
Author(s) / Creator(s):
;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Numerical Methods for Partial Differential Equations
Volume:
40
Issue:
4
ISSN:
0749-159X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We propose two mass and heat energy conservative, unconditionally stable, decoupled numerical algorithms for solving the Cahn–Hilliard–Navier–Stokes–Darcy–Boussinesq system that models thermal convection of two‐phase flows in superposed free flow and porous media. The schemes totally decouple the computation of the Cahn–Hilliard equation, the Darcy equations, the heat equation, the Navier–Stokes equations at each time step, and thus significantly reducing the computational cost. We rigorously show that the schemes are conservative and energy‐law preserving. Numerical results are presented to demonstrate the accuracy and stability of the algorithms.

     
    more » « less
  2. Abstract

    In this paper, we establish the fully decoupled numerical methods by utilizing scalar auxiliary variable approach for solving Cahn–Hilliard–Darcy system. We exploit the operator splitting technique to decouple the coupled system and Galerkin finite element method in space to construct the fully discrete formulation. The developed numerical methods have the features of second order accuracy, totally decoupling, linearization, and unconditional energy stability. The unconditionally stability of the two proposed decoupled numerical schemes are rigorously proved. Abundant numerical results are reported to verify the accuracy and effectiveness of proposed numerical methods.

     
    more » « less
  3. Abstract

    This paper concentrates on a priori error estimates of two monolithic schemes for Biot's consolidation model based on the three‐field formulation introduced by Oyarzúa et al. (SIAM J Numer Anal, 2016). The spatial discretizations are based on the Taylor–Hood finite elements combined with Lagrange elements for the three primary variables. We employ two different schemes to discretize the time domain. One uses the backward Euler method, and the other applies the combination of the backward Euler and Crank‐Nicolson methods. A priori error estimates show that both schemes are unconditionally convergent with optimal error orders. Detailed numerical experiments are presented to validate the theoretical analysis.

     
    more » « less
  4. In this paper, we consider numerical approximations for the anisotropic Cahn–Hilliard equation. We develop two linear and second-order schemes that combine the IEQ approach with the stabilization technique, where several extra linear stabilization terms are added in and they can be shown to be crucial to suppress the non-physical spatial oscillations caused by the strong anisotropy. We show the well-posedness of the resulting linear systems and further prove their corresponding unconditional energy stabilities rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy, and efficiency of the proposed schemes. 
    more » « less
  5. We consider numerical approximations for a phase-field dendritic crystal growth model, which is a highly nonlinear system that couples the anisotropic Allen–Cahn type equation and the heat equation. By combining the stabilized-Invariant Energy Quadratization method with a novel decoupling technique, the scheme requires solving only a sequence of linear elliptic equations at each time step, making it the first, to the best of the author’s knowledge, totally decoupled, linear, unconditionally energy stable scheme for the model. We further prove the unconditional energy stability rigorously and present various numerical simulations to demonstrate the stability and accuracy. 
    more » « less