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
  • Fully decoupled unconditionally stable Crank–Nicolson leapfrog numerical methods for the Cahn–Hilliard–Darcy system
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. ; ; ; (, Numerical Methods for Partial Differential Equations)
    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. ; ; (, Numerical Methods for Partial Differential Equations)
    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. ; ; ; (, Journal of Scientific Computing)
    Abstract We consider a family of variable time-stepping Dahlquist-Liniger-Nevanlinna (DLN) schemes, which is unconditionally non-linear stable and second order accurate, for the Allen-Cahn equation. The finite element methods are used for the spatial discretization. For the non-linear term, we combine the DLN scheme with two efficient temporal algorithms: partially implicit modified algorithm and scalar auxiliary variable algorithm. For both approaches, we prove the unconditional, long-term stability of the model energy under any arbitrary time step sequence. Moreover, we provide rigorous error analysis for the partially implicit modified algorithm with variable time-stepping. Efficient time-adaptive algorithms based on these schemes are also proposed. Several one- and two-dimensional numerical tests are presented to verify the properties of the proposed time-adaptive DLN methods. 
    more » « less
  4. ; ; (, Journal of Scientific Computing)
    In contrast to the classical Allen-Cahn equation, the conservative Allen-Cahn equation with a nonlocal Lagrange multiplier not only satisfies the maximum bound principle (MBP) and energy dissipation law but also ensures mass conservation. Many existing schemes often fail to preserve all these properties at the discrete level or require high regularity in time on the exact solution for convergence analysis. In this paper, we construct a new class of low regularity integrators (LRIs) for time discretization of the conservative Allen-Cahn equation by repeatedly using Duhamel's formula. The proposed first- and second-order LRI schemes are shown to conserve mass unconditionally and satisfy the MBP under some time step size constraints. Temporal error estimates for these schemes are derived under a low regularity assumption that the exact solution is only Lipschitz continuous in time, followed by a rigorous proof for energy stability of the corresponding time-discrete solutions. Various numerical experiments and comparisons in two and three dimensions are presented to verify the theoretical results and illustrate the performance of the LRI schemes, especially when the interfacial parameter approaches zero. 
    more » « less
  5. ; ; ; (, 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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email