scholarly journals Adaptive Mesh Refinement Strategies for a Novel Modelof Immiscible Fluid Flow in Fractures

2021 ◽  
Author(s):  
Sobhan Hatami ◽  
Stuart Walsh

In this paper, we consider two Adaptive Mesh Refinement (AMR) methods to simulate flow through fractures using a novel multiphase model. The approach represents the fluid using a two-dimensional parallel-plate model that employs techniques adapted from lattice-Boltzmann simulations to track the fluid interface. Here, we discuss different mesh refinement strategies for the model and compare their performance to that of a uniform grid. Results from the simulations are demonstrated showing excellent agreement between the model and analytical solutions for both unrefined and refined meshes. We also present results from the study that illustrate the behavior of the AMR front-tracking method. The AMR model is able to accurately track the interfacial properties in cases where uniform fine meshes would significantly increase the simulation cost.The ability of the model to dynamically refine the domain is demonstrated by presenting the results from an example with evolving interfaces.

2018 ◽  
Author(s):  
Thiago Dias dos Santos ◽  
Mathieu Morlighem ◽  
Hélène Seroussi ◽  
Philippe Remy Bernard Devloo ◽  
Jefferson Cardia Simões

Abstract. Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of the grounding line indeed plays a major role in ice sheet dynamics, as it is a fundamental control on marine ice sheet stability. Numerical modeling of grounding line requires significant computational resources since the accuracy of its position depends on grid or mesh resolution. A technique that improves accuracy with reduced computational cost is the adaptive mesh refinement approach, AMR. We present here the implementation of the AMR technique in the finite element Ice Sheet System Model (ISSM) to simulate grounding line dynamics under two different benchmarks, MISMIP3d and MISMIP+. We test different refinement criteria: (a) distance around grounding line, (b) a posteriori error estimator, the Zienkiewicz-Zhu (ZZ) error estimator, and (c) different combinations of (a) and (b). We find that for MISMIP3d setup, refining 5 km around the grounding line, both on grounded and floating ice, is sufficient to produce AMR results similar to the ones obtained with uniformly refined meshes. However, for the MISMIP+ setup, we note that there is a minimum distance of 30 km around the grounding line required to produce accurate results. We find this AMR mesh-dependency is linked to the complex bedrock topography of MISMIP+. In both benchmarks, the ZZ error estimator presents high values around the grounding line. Particularly for MISMIP+ setup, the estimator also presents high values in the grounded part of the ice sheet, following the complex shape of the bedrock geometry. This estimator helps guide the refinement procedure such that AMR performance is improved. Our results show that computational time with AMR depends on the required accuracy, but in all cases, it is significantly shorter than for uniformly refined meshes. We conclude that AMR without an associated error estimator should be avoided, especially for real glaciers that have a complex bed geometry.


2018 ◽  
Vol 50 (04) ◽  
pp. 561-570
Author(s):  
I. A. QAZI ◽  
A. F. ABBASI ◽  
M. S. JAMALI ◽  
INTIZAR INTIZAR ◽  
A. TUNIO ◽  
...  

Author(s):  
Alexander Haberl ◽  
Dirk Praetorius ◽  
Stefan Schimanko ◽  
Martin Vohralík

AbstractWe consider a second-order elliptic boundary value problem with strongly monotone and Lipschitz-continuous nonlinearity. We design and study its adaptive numerical approximation interconnecting a finite element discretization, the Banach–Picard linearization, and a contractive linear algebraic solver. In particular, we identify stopping criteria for the algebraic solver that on the one hand do not request an overly tight tolerance but on the other hand are sufficient for the inexact (perturbed) Banach–Picard linearization to remain contractive. Similarly, we identify suitable stopping criteria for the Banach–Picard iteration that leave an amount of linearization error that is not harmful for the residual a posteriori error estimate to steer reliably the adaptive mesh-refinement. For the resulting algorithm, we prove a contraction of the (doubly) inexact iterates after some amount of steps of mesh-refinement/linearization/algebraic solver, leading to its linear convergence. Moreover, for usual mesh-refinement rules, we also prove that the overall error decays at the optimal rate with respect to the number of elements (degrees of freedom) added with respect to the initial mesh. Finally, we prove that our fully adaptive algorithm drives the overall error down with the same optimal rate also with respect to the overall algorithmic cost expressed as the cumulated sum of the number of mesh elements over all mesh-refinement, linearization, and algebraic solver steps. Numerical experiments support these theoretical findings and illustrate the optimal overall algorithmic cost of the fully adaptive algorithm on several test cases.


Author(s):  
Weiqun Zhang ◽  
Andrew Myers ◽  
Kevin Gott ◽  
Ann Almgren ◽  
John Bell

Block-structured adaptive mesh refinement (AMR) provides the basis for the temporal and spatial discretization strategy for a number of Exascale Computing Project applications in the areas of accelerator design, additive manufacturing, astrophysics, combustion, cosmology, multiphase flow, and wind plant modeling. AMReX is a software framework that provides a unified infrastructure with the functionality needed for these and other AMR applications to be able to effectively and efficiently utilize machines from laptops to exascale architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving accurate descriptions of different physical processes in complex multiphysics algorithms. AMReX supports algorithms that solve systems of partial differential equations in simple or complex geometries and those that use particles and/or particle–mesh operations to represent component physical processes. In this article, we will discuss the core elements of the AMReX framework such as data containers and iterators as well as several specialized operations to meet the needs of the application projects. In addition, we will highlight the strategy that the AMReX team is pursuing to achieve highly performant code across a range of accelerator-based architectures for a variety of different applications.


Sign in / Sign up

Export Citation Format

Share Document