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
  • Implementation of the Accurate Conservative Phase Field Method for two-phase incompressible flows in a finite volume framework, Proc. of the 9th Int. and 49th National Conf. on Fluid Mechanics and Fluid Power (FMFP), December 14-16, 2022, IIT Roorkee, Roorkee-247667, Uttarakhand, India. FMFP2022–1300.
Title: Implementation of the Accurate Conservative Phase Field Method for two-phase incompressible flows in a finite volume framework, Proc. of the 9th Int. and 49th National Conf. on Fluid Mechanics and Fluid Power (FMFP), December 14-16, 2022, IIT Roorkee, Roorkee-247667, Uttarakhand, India. FMFP2022–1300.
The phase field method provides a simple mass conserving method for solving two-phase immiscible - incompressible Navier-Stokes Equations. The relative ease in implementing this method compared to other interface reconstruction methods, coupled with its conservativeness and boundedness makes it an attractive alternative. We implement the method in a parallel structured multi-block generalized coordinate finite volume solver using a collocated grid arrangement within the framework of the fractional-step method. The discretization uses a second-order central difference method for both the Navier-Stokes and the phase field equations. A TVD-based averaging technique is used for calculating density at cell faces in the pressure correction step to handle high-density ratios. The simulation framework is verified in standard test cases: Zalesak Disk, a droplet in shear flow, Solitary Wave Runup, Rayleigh Taylor Instability, and the Dam Break Problem. A second-order rate of convergence and excellent phase volume conservation is observed.  more » « less
Award ID(s):
2152694
PAR ID:
10436309
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Proc. of the 9th Int. and 49th National Conf. on Fluid Mechanics and Fluid Power (FMFP)
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 We investigate the temporal accuracy of two generalized‐ schemes for the incompressible Navier‐Stokes equations. In a widely‐adopted approach, the pressure is collocated at the time steptn + 1while the remainder of the Navier‐Stokes equations is discretized following the generalized‐ scheme. That scheme has been claimed to besecond‐order accurate in time. We developed a suite of numerical code using inf‐sup stable higher‐order non‐uniform rational B‐spline (NURBS) elements for spatial discretization. In doing so, we are able to achieve high spatial accuracy and to investigate asymptotic temporal convergence behavior. Numerical evidence suggests that onlyfirst‐order accuracyis achieved, at least for the pressure, in this aforesaid temporal discretization approach. On the other hand, evaluating the pressure at the intermediate time step recovers second‐order accuracy, and the numerical implementation is simplified. We recommend this second approach as the generalized‐ scheme of choice when integrating the incompressible Navier‐Stokes equations. 
    more » « less
  3. ; ; ; (, 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
  4. ; (, Mathematical Models and Methods in Applied Sciences)
    null (Ed.)
    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
  5. ; ; (, International journal of numerical analysis and modeling)
    This report presents a low computational and cognitive complexity, stable, time accurate and adaptive method for the Navier-Stokes equations. The improved method requires a minimally intrusive modification to an existing program based on the fully implicit / backward Euler time discretization, does not add to the computational complexity, and is conceptually simple. The backward Euler approximation is simply post-processed with a two-step, linear time filter. The time filter additionally removes the overdamping of Backward Euler while remaining unconditionally energy stable, proven herein. Even for constant stepsizes, the method does not reduce to a standard / named time stepping method but is related to a known 2-parameter family of A-stable, two step, second order methods. Numerical tests confirm the predicted convergence rates and the improved predictions of flow quantities such as drag and lift. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email