A Generalized Finite Volume Method for Density Driven Flows in Porous Media

In this article, we consider a time evolution equation for solute transport, coupled with a pressure equation in space dimension 2. For the numerical discretization, we combine the generalized finite volume method SUSHI on adaptive meshes with a time semi-implicit scheme. In the first part of this article, we present numerical simulations for two problems: a rotating interface between fresh and salt water and a well-known test case proposed by Henry. In the second part, we also introduce heat transfer and perform simulations for a system from the documentation of the software SEAWAT.

Metrics for Performance Quantification of Adaptive Mesh Refinement

Non-uniform, dynamically adaptive meshes are a useful tool for reducing computational complexities for geophysical simulations that exhibit strongly localised features such as is the case for tsunami, hurricane or typhoon prediction. Using the example of a shallow water solver, this study explores a set of metrics as a tool to distinguish the performance of numerical methods using adaptively refined versus uniform meshes independent of computational architecture or implementation. These metrics allow us to quantify how a numerical simulation benefits from the use of adaptive mesh refinement. The type of meshes we are focusing on are adaptive triangular meshes that are non-uniform and structured. Refinement is controlled by physics-based indicators that capture relevant physical processes and determine the areas of mesh refinement and coarsening. The proposed performance metrics take into account a number of characteristics of numerical simulations such as numerical errors, spatial resolution, as well as computing time. Using a number of test cases we demonstrate that correlating different quantities offers insight into computational overhead, the distribution of numerical error across various mesh resolutions as well as the evolution of numerical error and run-time per degree of freedom.

Ensemble Kalman Filter for non-conservative moving mesh solvers with a joint physics and mesh location update

<p>Numerical solvers using adaptive meshes can focus computational power on important regions of a model domain capturing important or unresolved physics. The adaptation can be informed by the model state, external information, or made to depend on the model physics. <br> In this latter case, one can think of the mesh configuration  as part of the model state. If observational data is to be assimilated into the model, the question of updating the mesh configuration with the physical values arises. Adaptive meshes present significant challenges when using popular ensemble Data Assimilation (DA) methods. We develop a novel strategy for ensemble-based DA for which the adaptive mesh is updated along with the physical values. This involves including the node locations as a part of the model state itself allowing them to be updated automatically at the analysis step. This poses a number of challenges which we resolve to produce an effective approach that promises to apply with some generality. We evaluate our strategy with two testbed models in 1-d comparing to a strategy that we previously developed that does not update the mesh configuration. We find updating the mesh improves the fidelity and convergence of the filter. An extensive analysis on the performance of our scheme beyond just the RMSE error is also presented.</p>

Stable Dynamical Adaptive Mesh Refinement

We consider accurate and stable interpolation procedures for numerical simulations utilizing time dependent adaptive meshes. The interpolation of numerical solution values between meshes is considered as a transmission problem with respect to the underlying semi-discretized equations, and a theoretical framework using inner product preserving operators is developed, which allows for both explicit and implicit implementations. The theory is supplemented with numerical experiments demonstrating practical benefits of the new stable framework. For this purpose, new interpolation operators have been designed to be used with multi-block finite difference schemes involving non-collocated, moving interfaces.