Computational Mathematics (CM): Recent submissions
Now showing items 120 of 165

Fast inversion of loggingwhiledrilling resistivity measurements acquired in multiple wells
(Geophysics, 201610)This paper introduces a new method for the fast inversion of borehole resistivity measurements acquired in multiple wells using loggingwhiledrilling (LWD) instruments. There are two key novel contributions. First, we ... 
Finite element approximation of electromagnetic fields using nonfitting meshes for Geophysics
(SIAM Journal on Numerical Analysis, 201807)We analyze the use of nonfitting meshes for simulating the propagation of electromagnetic waves inside the earth with applications to borehole logging. We avoid the use of parameter homogenization and employ standard edge ... 
TimeDomain GoalOriented Adaptivity Using PseudoDual Error Representations
(Computer Methods in Applied Mechanics and Engineering, 201712)Goaloriented adaptive algorithms produce optimal grids to solve challenging engineering problems. Recently, a novel error representation using (unconventional) pseudodual problems for goaloriented adaptivity in the ... 
ExplicitinTime GoalOriented Adaptivity
(Computer Methods in Applied Mechanics and Engineering, 20190415)Goaloriented adaptivity is a powerful tool to accurately approximate physically relevant solution features for partial differential equations. In time dependent problems, we seek to represent the error in the quantity of ... 
Normal lubrication force between spherical particles immersed in a shearthickening fluid
(Physics of Fluids, 20181207)In this work, the inverse biviscous model [Physics of Fluids 29, 103104 (2017)] is used to describe a shearthickening fluid. An analytical velocity profile in a planar Poiseuille flow is utilized to calculate an approximate ... 
Apparent slip mechanism between two spheres based on solvent rheology: Theory and implication for the shear thinning of nonBrownian suspensions
(Physical Review Fluids, 20181210)Analytical results for the apparent slip between two spheres in a simple biviscous model of a shear thinning fluid are presented. Velocity profiles and apparent slip lengths along the surfaces are analyzed in order to ... 
Source time reversal (STR) method for linear elasticity
(Computers & Mathematics with Applications, 2018)We study the problem of source reconstruction for a linear elasticity problem applied to seismicity induced by mining. We assume the source is written as a variable separable function $\mathbf{f(x)}\>g(t)$ . We first present ... 
Nonlinear computations of heave motions for a generic Wave Energy Converter
(In Vittorio Bucci Alberto Marinò (Eds.), Technology and Science for the Ships of the Future, Proceedings of NAV 2018: 19th International Conference on Ship and Maritime Research, IOS Press BV, 2018)A benchmarking activity of numerical methods for analysis of Wave Energy Converters (WEC) was proposed under the Ocean Energy Systems (OES) International Energy Agency (IEA) Task 10 in 2015. The purpose of the benchmark ... 
Finite Element Simulations of Twophase Flow and Floating Bodies Using FEniCSHPC
(13th World Congress in Computational Mechanics WCCM, New York, USA, 201807)We present a variational multiscale stabilized finite element method to solve the variable density incompressible NavierStokes equations for the simulation of twophase flow. We introduce a levelset method based on the ... 
Adaptive Direct FEM Simulation with Unicorn/FEniCSHPC for CS1
(5th International Workshop on High Order CFD Methods, 20180101) 
Fast Onedimensional Finite Element Approximation of Geophysical Measurements
(2018)There exist a wide variety of geophysical prospection methods. In this work, we focus on resistivity methods. We categorize these resistivity prospection methods according to their acquisition location as (a) on the surface, ... 
Gaussian quadrature for $C^1$ cubic CloughTocher macrotriangles
(Journal of Computational and Applied Mathematics, 20181031)A numerical integration rule for multivariate cubic polynomials over ndimensional simplices was designed by Hammer and Stroud [14]. The quadrature rule requires n + 2 quadrature points: the barycentre of the simplex and ... 
Generalization of the Pythagorean Eigenvalue Error Theorem and its Application to Isogeometric Analysis
(Numerical Methods for PDEs, 20181013)This chapter studies the effect of the quadrature on the isogeometric analysis of the wave propagation and structural vibration problems. The dispersion error of the isogeometric elements is minimized by optimally blending ... 
International energy agency ocean energy systems task 10 wave energy converter modeling verification and validation
(Proceedings of European Wave and Tidal Energy Conference, 2017)This is the first joint reference paper for the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling verification and validation group. The group is established under the OES Energy Technology Network program ... 
A partition of unity finite element method for computational diffusion MRI
(Journal of Computational Physics, 2018)The Bloch–Torrey equation describes the evolution of the spin (usually water proton) magnetization under the influence of applied magnetic field gradients and is commonly used in numerical simulations for diffusion MRI ... 
Simulation of the HarshLab floating platform for offshore experimentation using FEniCSHPC
(Poster EUSKAMPUS IV Marine Energy Conference, 2017) 
Direct finite element simulation of turbulent flow for marine based renewable energy
(Companion paper to F. Wendt, et. al. International energy agency ocean energy systems task 10 wave energy converter modeling verification and validation. In 12th European Wave and Tidal Energy Conference, 2017., 2017)In this article we present a computational framework for simulation of turbulent flow in marine based renewable energy applications. In particular, we focus on floating structures and rotating turbines. This work is an ... 
A Numerical 1.5D Method for the Rapid Simulation of Geophysical Resistivity Measurements
(Geosciences, 20180614)In some geological formations, borehole resistivity measurements can be simulated using a sequence of 1D models. By considering a 1D layered media, we can reduce the dimensionality of the problem from 3D to 1.5D via a ... 
Numerical simulation of the von Kármán sodium dynamo experiment
(Journal of Fluid Mechanics, 20180903)We present hydrodynamic and magnetohydrodynamic (MHD) simulations of liquid sodium flows in the von Kármán sodium (VKS) setup. The counterrotating impellers made of soft iron that were used in the successful 2006 experiment ... 
Goaloriented adaptivity using unconventional error representations for the multidimensional Helmholtz equation
(International Journal for Numerical Methods in Engineering, 20170627)In goal‐oriented adaptivity, the error in the quantity of interest is represented using the error functions of the direct and adjoint problems. This error representation is subsequently bounded above by element‐wise error ...