Journal of Computational Physics
The multidimensional positive definite advection transport algorithm: nonoscillatory option
Journal of Computational Physics
A synchronous and iterative flux-correction formalism for coupled transport equations
Journal of Computational Physics
MPDATA: a finite-difference solver for geophysical flows
Journal of Computational Physics
Antidiffusive Velocities for Multipass Donor Cell Advection
SIAM Journal on Scientific Computing
Second-order sign-preserving conservative interpolation (remapping) on general grids
Journal of Computational Physics
Numerical solution of the reaction-advection-diffusion equation on the sphere
Journal of Computational Physics
Building resolving large-eddy simulations and comparison with wind tunnel experiments
Journal of Computational Physics
Iterated upwind schemes for gas dynamics
Journal of Computational Physics
An edge-based unstructured mesh discretisation in geospherical framework
Journal of Computational Physics
SIAM Journal on Scientific Computing
Using blue gene/p and GPUs to accelerate computations in the EULAG model
LSSC'11 Proceedings of the 8th international conference on Large-Scale Scientific Computing
Parallelization of EULAG model on multicore architectures with GPU accelerators
PPAM'11 Proceedings of the 9th international conference on Parallel Processing and Applied Mathematics - Volume Part II
EULAG, a computational model for multiscale flows: An MHD extension
Journal of Computational Physics
An unstructured-mesh atmospheric model for nonhydrostatic dynamics
Journal of Computational Physics
Journal of Computational Physics
Hi-index | 31.48 |
We present an advancement in the evolution of MPDATA (multidimensional positive definite advection transport algorithm). Over the last two decades, MPDATA has proven successful in applications using single-block structured cuboidal meshes (viz. Cartesian meshes), while employing homeomorphic mappings to accommodate time-dependent curvilinear domains. Motivated by the strengths of the Cartesian-mesh MPDATA, we develop a new formulation in an arbitrary finite-volume framework with a fully unstructured polyhedral hybrid mesh. In MPDATA, as in any Taylor-series based integration method for PDE, the choice of data structure has a pronounced impact on the technical details of the algorithm. Aiming at a broad range of applications with a large number of control-volume cells, we select a general, compact and computationally efficient, edge-based data structure. This facilitates the use of MPDATA for problems involving complex geometries and/or inhomogeneous anisotropic flows where mesh adaptivity is advantageous. In this paper, we describe the theory and implementation of the basic finite-volume MPDATA, and document extensions important for applications: a fully monotone scheme, diffusion scheme, and generalization to complete flow solvers. Theoretical discussions are illustrated with benchmark results in two and three spatial dimensions.