Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

Ido Lavi,Nicolas Meunier,Olivier Pantz

2023-05-10

Abstract:In this work, we propose and compare three numerical methods to handle the one-phase Hele-Shaw problem with surface tension in dimension two by using three variational approaches in the spirit of the seminal works \cite{Otto, Gia_Otto}.

Numerical Analysis

What problem does this paper attempt to address?

Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

Nicolas Meunier,Ido Lavi,Olivier Pantz

DOI: https://doi.org/10.1051/m2an/2024005

2024-02-21

ESAIM: Mathematical Modelling and Numerical Analysis

Abstract:In this work, we propose and compare three numerical methods to handle the one-phase Hele-Shaw problem with surface tension in dimension two by using three variational approaches in the spirit of the seminal works \cite{Otto, Gia_Otto}.

mathematics, applied
A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

Z. H. Gu,H. L. Wen,Yuan Yao,C. H. Yu

DOI: https://doi.org/10.1080/10407790.2019.1642048

2019-01-01

Abstract:In this article, a two-phase flow algorithm that uses an explicit Adams-Bashforth scheme coupled with volume of fluid (VOF) method in a uniform staggered Cartesian grid for surface tension dominant simulation is presented. The interface reconstruction procedure is implemented using the Weighed Linear Interface Calculation (WLIC) algorithm. The level set (LS) function instead of VOF function is adopted for calculation of not only the interface normal vectors in the WLIC algorithm but also the surface tension terms in the Navier-Stokes equations. This VOF method is compared with LS method and other experiment results in the literature for vortex deforming, milk crown, droplet impacting into deep liquid pool and bubble bursting at a free surface problems to show whether mass can be well conserved and breaking and merging phenomena can be accurately simulated.
Removing the stiffness from interfacial flows with surface tension

Thomas Y. Hou,John S. Lowengrub,Michael J. Shelley

DOI: https://doi.org/10.1006/jcph.1994.1170

IF: 4.645

1994-10-01

Journal of Computational Physics

Abstract:A new formulation and new methods are presented for computing the motion of fluid interfaces with surface tension in two-dimensional, irrotational, and incompressible fluids. Through the Laplace-Young condition at the interface, surface tension introduces high-order terms, both nonlinear and nonlocal, into the dynamics. This leads to severe stability constraints for explicit time integration methods and makes the application of implicit methods difficult. This new formulation has all the nice properties for time integration methods that are associated with having a linear, constant coefficient, highest order term. That is, using this formulation, we give implicit time integration methods that have no high order time step stability constraint associated with surface tension and are explicit in Fourier space. The approach is based on a boundary integral formulation and applies more generally, even to problems beyond the fluid mechanical context. Here they are applied to computing with high resolution the motion of interfaces in Hele-Shaw flows and the motion of free surfaces in inviscid flows governed by the Euler equations. One Hele-Shaw computation shows the behavior of an expanding gas bubble over long-time as the interface undergoes successive tip-splittings and finger competition. A second computation shows the formation of a very ramified interface through the interaction of surface tension with an unstable density stratification. In Euler flows, the computation of a vortex sheet shows its roll-up through the Kelvin-Helmholtz instability. This motion culminates in the late time self-intersection of the interface, creating trapped bubbles of fluid. This is, we believe, a type of singularity formation previously unobserved for such flows in 2D. Finally, computations of falling plumes in an unstably stratified Boussinesq fluid show a very similar behavior.

physics, mathematical,computer science, interdisciplinary applications
Local Smoothing Effect and Existence for the One-Phase Hele–Shaw Problem with Zero Surface Tension

Xuming Xie

DOI: https://doi.org/10.1080/17476933.2011.575463

2012-01-01

Complex Variables and Elliptic Equations

Abstract:We study an initial value problem for the one-phase Hele-Shaw problem with zero surface tension. We establish local well-posedness for the initial value problem in Sobolev space. Furthermore, we obtain that, on average in time, the solution gains 1/2 derivative of smoothness in spatial variable compared to the initial data.
Second Order Decoupled Implicit/Explicit Method Of The Primitive Equations Of The Ocean I: Time Discretization

Yinnian He

2015-01-01

International Journal of Numerical Analysis and Modeling

Abstract:In this article, we propose the time discretization of the second order decoupled implicit/explicit method of the 3D primitive equations of the ocean in the case of the Dirichlet boundary conditions on the side. We deduce the second order optimal error estimates on the L-2 and H-1 norms of the time discrete velocity and density and the L-2 norm of the time discrete pressure under the restriction of the time step 0 < tau <= beta for some positive constant beta. Also, we deduce some stability results on the time discrete solution under the same restriction on the time step.
A Second-Order Decoupled Implicit/Explicit Method of the 3d Primitive Equations of Ocean Ii: Finite Element Spatial Discretization

Yinnian He,Hui Xu,Zhangxin Chen

DOI: https://doi.org/10.1002/nme.5235

IF: 3.021

2016-01-01

International Journal for Numerical Methods in Engineering

Abstract:A fully discrete second-order decoupled implicit/explicit method is proposed for solving 3D primitive equations of ocean in the case of Dirichlet boundary conditions on the side, where a second-order decoupled implicit/explicit scheme is used for time discretization, and a finite element method based on the P-1 (P-1) - P-1 - P-1 (P-1) elements for velocity, pressure and density is used for spatial discretization of these primitive equations. Optimal H-1 - L-2 - H-1 error estimates for numerical solution (u(h)(n), p(h)(n), theta(n)(h)) and an optimal L-2 error estimate for (u(h)(n), theta(n)(h)) are established under the convergence condition of 0 < h <= beta(1), 0 < tau <= beta(2), and tau <= beta(3)h for some positive constants beta(1), beta(2), and beta(3). Furthermore, numerical computations show that the H-1 - L-2 - H-1 convergence rate for numerical solution (u(h)(n), p(h)(n), theta(n)(h)) is of O (h + tau(2)) and an L-2 convergence rate for (u(h)(n), theta(n)(h)) is O (h(2) + tau(2)) with the assumed convergence condition, where h is a mesh size and tau is a time step size. More practical calculations are performed as a further validation of the numerical method. Copyright (C) 2016 John Wiley & Sons, Ltd.
A new family of semi-implicit Finite Volume / Virtual Element methods for incompressible flows on unstructured meshes

Walter Boscheri,Andrea Chiozzi,Michele Giuliano Carlino,Giulia Bertaglia

DOI: https://doi.org/10.1016/j.cma.2023.116140

2023-04-15

Abstract:We introduce a new family of high order accurate semi-implicit schemes for the solution of non-linear hyperbolic partial differential equations on unstructured polygonal meshes. The time discretization is based on a splitting between explicit and implicit terms that may arise either from the multi-scale nature of the governing equations, which involve both slow and fast scales, or in the context of projection methods, where the numerical solution is projected onto the physically meaningful solution manifold. We propose to use a high order finite volume (FV) scheme for the explicit terms, ensuring conservation property and robustness across shock waves, while the virtual element method (VEM) is employed to deal with the discretization of the implicit terms, which typically requires an elliptic problem to be solved. The numerical solution is then transferred via suitable L2 projection operators from the FV to the VEM solution space and vice-versa. High order time accuracy is achieved using the semi-implicit IMEX Runge-Kutta schemes, and the novel schemes are proven to be asymptotic preserving and well-balanced. As representative models, we choose the shallow water equations (SWE), thus handling multiple time scales characterized by a different Froude number, and the incompressible Navier-Stokes equations (INS), which are solved at the aid of a projection method to satisfy the solenoidal constraint of the velocity field. Furthermore, an implicit discretization for the viscous terms is devised for the INS model, which is based on the VEM technique. Consequently, the CFL-type stability condition on the maximum admissible time step is based only on the fluid velocity and not on the celerity nor on the viscous eigenvalues. A large suite of test cases demonstrates the accuracy and the capabilities of the new family of schemes to solve relevant benchmarks in the field of incompressible fluids.

Numerical Analysis,Computational Engineering, Finance, and Science,Fluid Dynamics
High order Finite Difference/Discontinuous Galerkin schemes for the incompressible Navier-Stokes equations with implicit viscosity

Walter Boscheri,Maurizio Tavelli,Nicola Paoluzzi

DOI: https://doi.org/10.2478/caim-2022-0003

2022-01-01

Communications in Applied and Industrial Mathematics

Abstract:Abstract In this work we propose a novel numerical method for the solution of the incompressible Navier-Stokes equations on Cartesian meshes in 3D. The semi-discrete scheme is based on an explicit discretization of the nonlinear convective flux tensor and an implicit treatment of the pressure gradient and viscous terms. In this way, the momentum equation is formally substituted into the divergence-free constraint, thus obtaining an elliptic equation on the pressure which eventually maintains at the discrete level the involution on the divergence of the velocity field imposed by the governing equations. This makes our method belonging to the class of so-called structure-preserving schemes. High order of accuracy in space is achieved using an efficient CWENO reconstruction operator that is exploited to devise a conservative finite difference scheme for the convective terms. Implicit central finite differences are used to remove the numerical dissipation in the pressure gradient discretization. To avoid the severe time step limitation induced by the viscous eigenvalues related to the parabolic terms in the governing equations, we propose to devise an implicit local discontinuous Galerkin (DG) solver. The resulting viscous sub-system is symmetric and positive definite, therefore it can be efficiently solved at the aid of a matrix-free conjugate gradient method. High order in time is granted by a semi-implicit IMEX time stepping technique. Convergence rates up to third order of accuracy in space and time are proven, and a suite of academic benchmarks is shown in order to demonstrate the robustness and the validity of the novel schemes, especially in the context of high viscosity coefficients.
A fully-coupled algorithm with implicit surface tension treatment for interfacial flows with large density ratios

Romain Janodet,Berend van Wachem,Fabian Denner

DOI: https://doi.org/10.1016/j.jcp.2024.113520

2024-10-23

Abstract:The stability of most surface-tension-driven interfacial flow simulations is governed by the capillary time-step constraint. This concerns particularly small-scale flows and, more generally, highly-resolved liquid-gas simulations with moderate inertia. To date, the majority of interfacial-flow simulations are performed using an explicit surface-tension treatment, which restrains the performance of such simulations. Recently, an implicit treatment of surface tension able to breach the capillary time-step constraint using the volume-of-fluid (VOF) method was proposed, based on a fully-coupled pressure-based finite-volume algorithm. To this end, the interface-advection equation is incorporated implicitly into the linear flow solver, resulting in a tight coupling between all implicit solution variables (colour function, pressure, velocity). However, this algorithm is limited to uniform density and viscosity fields. Here, we present a fully-coupled algorithm for interfacial flows with implicit surface tension applicable to interfacial flows with large density and viscosity ratios. This is achieved by solving the continuity and momentum equations in conservative form, whereby the density is treated implicitly with respect to the colour function, and the advection term of the interface-advection equation is discretised using the THINC/QQ algebraic VOF scheme, yielding a consistent discretisation of the advective terms. This new algorithm is tested by considering representative surface-tension-dominated interfacial flows, including the Laplace equilibrium of a stationary droplet and the three-dimensional Rayleigh-Plateau instability of a liquid filament. The presented results demonstrate that interfacial flows with large density and viscosity ratios can be simulated and energy conservation is ensured, even with a time step larger than the capillary time-step constraint.

Fluid Dynamics,Computational Physics
Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. II. Numerical study

R. Folch,J. Casademunt,A. Hernandez-Machado,L. Ramirez-Piscina

DOI: https://doi.org/10.1103/PhysRevE.60.1734

1999-09-13

Abstract:We implement a phase-field simulation of the dynamics of two fluids with arbitrary viscosity contrast in a rectangular Hele-Shaw cell. We demonstrate the use of this technique in different situations including the linear regime, the stationary Saffman-Taylor fingers and the multifinger competition dynamics, for different viscosity contrasts. The method is quantitatively tested against analytical predictions and other numerical results. A detailed analysis of convergence to the sharp interface limit is performed for the linear dispersion results. We show that the method may be a useful alternative to more traditional methods.

Soft Condensed Matter,Materials Science,Statistical Mechanics,Fluid Dynamics
Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

Jason Goulding,Mehrnaz Ayazi,Tamar Shinar,Craig Schroeder

DOI: https://doi.org/10.1016/j.jcp.2024.113382

IF: 4.645

2024-09-03

Journal of Computational Physics

Abstract: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.

physics, mathematical,computer science, interdisciplinary applications
Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness

Helmut Abels,Harald Garcke,Andrea Poiatti

2024-07-21

Abstract:A Cahn-Hilliard-Navier-Stokes system for two-phase flow on an evolving surface with non-matched densities is derived using methods from rational thermodynamics. For a Cahn-Hilliard energy with a singular (logarithmic) potential short time well-posedness of strong solutions together with a separation property is shown, under the assumption of a priori prescribed surface evolution. The problem is reformulated with the help of a pullback to the initial surface. Then a suitable linearization and a contraction mapping argument for the pullback system are used. In order to deal with the linearized system, it is necessary to show maximal $L^2$-regularity for the surface Stokes operator in the case of variable viscosity and to obtain maximal $L^p$-regularity for the linearized Cahn-Hilliard system.

Analysis of PDEs
Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Global Well-Posedness

Helmut Abels,Harald Garcke,Andrea Poiatti

2024-08-14

Abstract:We show existence and uniqueness of strong solutions to a Navier-Stokes/Cahn-Hilliard type system on a given two-dimensional evolving surface in the case of different densities and a singular (logarithmic) potential. The system describes a diffuse interface model for a two-phase flow of viscous incompressible fluids on an evolving surface. We also establish the validity of the instantaneous strict separation property from the pure phases. To show these results we use our previous achievements on local well-posedness together with suitable novel regularity results for the convective Cahn-Hilliard equation. The latter allows to obtain higher-order energy estimates to extend the local solution globally in time. To this aim the time evolution of energy type quantities has to be calculated and estimated carefully.

Analysis of PDEs
The Euler Implicit/explicit Scheme for the 2D Time-Dependent Navier-Stokes Equations with Smooth or Non-Smooth Initial Data

Yinnian He

DOI: https://doi.org/10.1090/s0025-5718-08-02127-3

IF: 2.118

2008-01-01

Mathematics of Computation

Abstract:This paper considers the stability and convergence results for the Euler implicit/explicit scheme applied to the spatially discretized two-dimensional (2D) time-dependent Navier-Stokes equations. A Galerkin finite element spatial discretization is assumed, and the temporal treatment is implicit/explict scheme, which is implicit for the linear terms and explicit for the nonlinear term. Here the stability condition depends on the smoothness of the initial data u(0) is an element of H-alpha, i.e., the time step condition is tau <= C-0 in the case of alpha = 2, tau|log h| <= C-0 in the case of alpha = 1 and tau h(-2) = C-0 in the case of alpha = 0 for mesh size h and some positive constant C-0. We provide the H-2-stability of the scheme under the stability condition with alpha = 0, 1, 2 and obtain the optimal H-1 - L-2 error estimate of the numerical velocity and the optimal L-2 error estimate of the numerical pressure under the stability condition with alpha = 1, 2.
A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow

Harald Garcke,Michael Hinze,Christian Kahle

DOI: https://doi.org/10.48550/arXiv.1402.6524

2014-02-26

Abstract:A new time discretization scheme for the numerical simulation of two-phase flow governed by a thermodynamically consistent diffuse interface model is presented. The scheme is consistent in the sense that it allows for a discrete in time energy inequality. An adaptive spatial discretization is proposed that conserves the energy inequality in the fully discrete setting by applying a suitable post processing step to the adaptive cycle. For the fully discrete scheme a quasi-reliable error estimator is derived which estimates the error both of the flow velocity, and of the phase field. The validity of the energy inequality in the fully discrete setting is numerically investigated.

Analysis of PDEs,Numerical Analysis
Quinpi: Integrating stiff hyperbolic systems with implicit high order finite volume schemes

Gabriella Puppo,Matteo Semplice,Giuseppe Visconti

2024-07-16

Abstract:Many interesting physical problems described by systems of hyperbolic conservation laws are stiff, and thus impose a very small time-step because of the restrictive CFL stability condition. In this case, one can exploit the superior stability properties of implicit time integration which allows to choose the time-step only from accuracy requirements, and thus avoid the use of small time-steps. We discuss an efficient framework to devise high order implicit schemes for stiff hyperbolic systems without tailoring it to a specific problem. The nonlinearity of high order schemes, due to space- and time-limiting procedures which control nonphysical oscillations, makes the implicit time integration difficult, e.g.~because the discrete system is nonlinear also on linear problems. This nonlinearity of the scheme is circumvented as proposed in (Puppo et al., Comm.~Appl.~Math.~\& Comput., 2023) for scalar conservation laws, where a first order implicit predictor is computed to freeze the nonlinear coefficients of the essentially non-oscillatory space reconstruction, and also to assist limiting in time. In addition, we propose a novel conservative flux-centered a-posteriori time-limiting procedure using numerical entropy indicators to detect troubled cells. The numerical tests involve classical and artificially devised stiff problems using the Euler's system of gas-dynamics.

Numerical Analysis,Computational Physics
Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

Nathan Touroux,Masakiyo Kitazawa,Koichi Murase,Marlene Nahrgang

DOI: https://doi.org/10.1093/ptep/ptae058

2024-04-23

Progress of Theoretical and Experimental Physics

Abstract:Abstract We propose a new method to solve the relativistic hydrodynamic equations based on implicit Runge-Kutta methods with a locally optimized fixed-point iterative solver. For numerical demonstration, we implement our idea for ideal hydrodynamics using the one-stage Gauss-Legendre method as an implicit method. The accuracy and computational cost of our new method are compared with those of explicit ones for the (1+1)-dimensional Riemann problem, as well as the (2+1)-dimensional Gubser flow and event-by-event initial conditions for heavy-ion collisions generated by TRENTo. We demonstrate that the solver converges with only one iteration in most cases, and as a result, the implicit method requires a smaller computational cost than the explicit one at the same accuracy in these cases, while it may not converge with an unrealistically large Δt. By showing a relationship between the one-stage Gauss-Legendre method with the iterative solver and the two-step Adams-Bashforth method, we argue that our method benefits from both the stability of the former and the efficiency of the latter.

physics, multidisciplinary, particles & fields
Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models with dynamic boundary conditions

Andrea Giorgini,Patrik Knopf

DOI: https://doi.org/10.1007/s00021-023-00811-w

2023-10-24

Abstract:We derive a novel thermodynamically consistent Navier--Stokes--Cahn--Hilliard system with dynamic boundary conditions. This model describes the motion of viscous incompressible binary fluids with different densities. In contrast to previous models in the literature, our new model allows for surface diffusion, a variable contact angle between the diffuse interface and the boundary, and mass transfer between bulk and surface. In particular, this transfer of material is subject to a mass conservation law including both a bulk and a surface contribution. The derivation is carried out by means of local energy dissipation laws and the Lagrange multiplier approach. Next, in the case of fluids with matched densities, we show the existence of global weak solutions in two and three dimensions as well as the uniqueness of weak solutions in two dimensions.

Analysis of PDEs,Mathematical Physics
Implicit and semi-implicit second-order time stepping methods for the Richards equation

Sana Keita,Abdelaziz Beljadid,Yves Bourgault

DOI: https://doi.org/10.1016/j.advwatres.2020.103841

IF: 5.361

2021-02-01

Advances in Water Resources

Abstract:This study concerns numerical methods for efficiently solving the Richards equation where different weak formulations and computational techniques are analyzed. The spatial discretizations are based on standard or mixed finite element methods. Different implicit and semi-implicit temporal discretization techniques of second-order accuracy are studied. To obtain a linear system for the semi-implicit schemes, we propose second-order techniques using extrapolation formulas and/or semi-implicit Taylor approximations for the temporal discretization of nonlinear terms. A numerical convergence study and a series of numerical tests are performed to analyze efficiency and robustness of the different schemes. The developed scheme, based on the proposed temporal extrapolation techniques and the mixed formulation involving the saturation and pressure head and using the standard linear Lagrange element, performs better than other schemes based on the saturation and the flux and using the Raviart-Thomas elements. The proposed semi-implicit scheme is a good alternative when implicit schemes meet convergence issues.

water resources

Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

Removing the stiffness from interfacial flows with surface tension

Local Smoothing Effect and Existence for the One-Phase Hele–Shaw Problem with Zero Surface Tension

Second Order Decoupled Implicit/Explicit Method Of The Primitive Equations Of The Ocean I: Time Discretization

A Second-Order Decoupled Implicit/Explicit Method of the 3d Primitive Equations of Ocean Ii: Finite Element Spatial Discretization

A new family of semi-implicit Finite Volume / Virtual Element methods for incompressible flows on unstructured meshes

High order Finite Difference/Discontinuous Galerkin schemes for the incompressible Navier-Stokes equations with implicit viscosity

A fully-coupled algorithm with implicit surface tension treatment for interfacial flows with large density ratios

Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. II. Numerical study

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Global Well-Posedness

The Euler Implicit/explicit Scheme for the 2D Time-Dependent Navier-Stokes Equations with Smooth or Non-Smooth Initial Data

A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow

Quinpi: Integrating stiff hyperbolic systems with implicit high order finite volume schemes

Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models with dynamic boundary conditions

Implicit and semi-implicit second-order time stepping methods for the Richards equation