Abstract. Time-dependent simulations of ice sheets require two equations to be solved:the mass transport equation, derived from the conservation of mass, and thestress balance equation, derived from the conservation of momentum. The masstransport equation controls the advection of ice from the interior of the icesheet towards its periphery, thereby changing its geometry. Because it isbased on an advection equation, a stabilization scheme needs to beemployed when solved using the finite-element method. Several stabilizationschemes exist in the finite-element method framework, but their respectiveaccuracy and robustness have not yet been systematically assessed forglaciological applications. Here, we compare classical schemes used in thecontext of the finite-element method: (i) artificial diffusion, (ii)streamline upwinding, (iii) streamline upwind Petrov–Galerkin, (iv)discontinuous Galerkin, and (v) flux-corrected transport. We also look at thestress balance equation, which is responsible for computing the ice velocitythat “advects” the ice downstream. To improve the velocity computationaccuracy, the ice-sheet modeling community employs several sub-elementparameterizations of physical processes at the grounding line, the point wherethe grounded ice starts to float onto the ocean. Here, we introduce a newsub-element parameterization for the driving stress, the force that drives theice-sheet flow. We analyze the response of each stabilization scheme byrunning transient simulations forced by ice-shelf basal melt. The simulationsare based on an idealized ice-sheet geometry for which there is no influenceof bedrock topography. We also perform transient simulations of the AmundsenSea Embayment, West Antarctica, where real bedrock and surface elevations areemployed. In both idealized and real ice-sheet experiments, stabilizationschemes based on artificial diffusion lead systematically to a bias towardsmore mass loss in comparison to the other schemes and therefore should beavoided or employed with a sufficiently high mesh resolution in the vicinityof the grounding line. We also run diagnostic simulations to assess theaccuracy of the driving stress parameterization, which, in combination with anadequate parameterization for basal stress, provides improved numericalconvergence in ice speed computations and more accurate results.
more »
« less
A macroelement stabilization for mixed finite element/finite volume discretizations of multiphase poromechanics
Strong coupling between geomechanical deformation and multiphase fluid flow appears in a variety of geoscience applications. A common discretization strategy for these problems is a continuous Galerkin finite element scheme for the momentum balance equation and a finite volume scheme for the mass balance equations. When applied within a fully implicit solution strategy, however, this discretization is not intrinsically stable. In the limit of small time steps or low permeabilities, spurious oscillations in the piecewise-constant pressure field, i.e., checkerboarding, may be observed. Further, eigenvalues associated with the spurious modes will control the conditioning of the matrices and can dramatically degrade the convergence rate of iterative linear solvers. Here, we propose a stabilization technique in which the mass balance equations are supplemented with stabilizing flux terms on a macroelement basis. The additional stabilization terms are dependent on a stabilization parameter. We identify an optimal value for this parameter using an analysis of the eigenvalue distribution of the macroelement Schur complement matrix. The resulting method is simple to implement and preserves the underlying sparsity pattern of the original discretization. Another appealing feature of the method is that mass is exactly conserved on macroelements, despite the addition of artificial fluxes. The efficacy of the proposed technique is demonstrated with several numerical examples.
more »
« less
- Award ID(s):
- 1914780
- NSF-PAR ID:
- 10162828
- Date Published:
- Journal Name:
- Computational Geosciences
- ISSN:
- 1420-0597
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Summary In this paper, we propose and analyze two stabilized mixed finite element methods for the dual‐porosity‐Stokes model, which couples the free flow region and microfracture‐matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual‐porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods.
-
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
-
Abstract We present a new class of discontinuous Galerkin methods for the space discretization of the time-dependent Maxwell equations whose main feature is the use of time derivatives and/or time integrals in the stabilization part of their numerical traces.These numerical traces are chosen in such a way that the resulting semidiscrete schemes exactly conserve a discrete version of the energy.We introduce four model ways of achieving this and show that, when using the mid-point rule to march in time, the fully discrete schemes also conserve the discrete energy.Moreover, we propose a new three-step technique to devise fully discrete schemes of arbitrary order of accuracy which conserve the energy in time.The first step consists in transforming the semidiscrete scheme into a Hamiltonian dynamical system.The second step consists in applying a symplectic time-marching method to this dynamical system in order to guarantee that the resulting fully discrete method conserves the discrete energy in time.The third and last step consists in reversing the above-mentioned transformation to rewrite the fully discrete scheme in terms of the original variables.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s10596-020-09964-3
