This paper extends the author’s previous two-dimensional work with Ou and LeVeque to high-resolution finite volume modeling of systems of fluids and poroelastic media in three dimensions, using logically rectangular mapped grids. A method is described for calculating consistent cell face areas and normal vectors for a finite volume method on a general non-rectilinear hexahedral grid. A novel li...

We investigate computing models that are presented as families of finite computing devices with a uniformity condition on the entire family. Examples of such models include Boolean circuits, membrane systems, DNA computers, chemical reaction networks and tile assembly systems, and there are many others. However, in such models there are actually two distinct kinds of uniformity condition. The f...

The Arnold–Winther element successfully discretizes the Hellinger–Reissner variational formulation of linear elasticity; its development was one key early breakthroughs finite exterior calculus. Despite great utility, it is not available in standard software, because degrees freedom are preserved under Piola push-forward. In this work we apply novel transformation theory recently developed by K...

In this work we develop a high-resolution mapped-grid finite volume method code to model wave propagation in two dimensions in systems of multiple orthotropic poroelastic media and/or fluids, with curved interfaces between different media. We use a unified formulation to simplify modeling of the various interface conditions — open pores, imperfect hydraulic contact, or sealed pores — that may e...

Efficient numerical seismic wavefield modelling is a key component of modern seismic imaging techniques, such as reverse-time migration and full-waveform inversion. Finite difference methods are perhaps the most widely used numerical approach for forward modelling, and here we introduce a novel scheme for implementing finite difference by introducing a time-to-space wavelet mapping. Finite diff...

Memory-mapped transactions combine the advantages of both memory mapping and transactions to provide a programming interface for concurrently accessing data on disk without explicit I/O or locking operations. This interface enables a programmer to design a complex serial program that accesses only main memory, and with little to no modification, convert the program into correct code with multip...