Source author record

Lilia Krivodonova

Lilia Krivodonova appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Distributed, Parallel, and Cluster Computing math.NA Mathematical Software physics.comp-ph physics.flu-dyn

Catalog footprint

What is connected

3works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Discontinuous Galerkin methods on graphics processing units for nonlinear hyperbolic conservation laws

We present a novel implementation of the modal discontinuous Galerkin (DG) method for hyperbolic conservation laws in two dimensions on graphics processing units (GPUs) using NVIDIA's Compute Unified Device Architecture (CUDA). Both flexible and highly accurate, DG methods accommodate parallel architectures well as their discontinuous nature produces element-local approximations. High performance scientific computing suits GPUs well, as these powerful, massively parallel, cost-effective devices have recently included support for double-precision floating point numbers. Computed examples for Euler equations over unstructured triangle meshes demonstrate the effectiveness of our implementation on an NVIDIA GTX 580 device. Profiling of our method reveals performance comparable to an existing nodal DG-GPU implementation for linear problems.

preprint2016arXiv

Three-dimensional simulations of sound propagation in a trumpet with accurate mouthpiece shank geometry

The length and bore geometry of musical instruments directly influences the quality of sound that can be produced. In brass instruments, nonlinear effects from finite-amplitude wave propagation can lead to wave distortion giving sounds a brassy timbre [3, 5, 14, 20, 26]. In this paper, we propose a three-dimensional model to describe nonlinear wave propagation in a trumpet and investigate the importance of the mouthpiece shank geometry. Time pressure waveforms corresponding to B_3(b) and B_4(b) notes were recorded at the mouthpiece shank and used as inputs for our model. To describe the motion of compressible inviscid fluid, we numerically solved the compressible Euler equations using the discontinuous Galerkin method. To validate our approach, the numerical results were compared to the recorded musical notes outside the bell of the trumpet. Simulations were performed on computational trumpets where different bore geometries were considered. Our results demonstrate that the shape of the narrow region near mouthpiece greatly influences the wave propagation and accuracy of the trumpet model.

preprint2015arXiv

Spatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations

We apply the discontinuous Galerkin finite element method with a degree $p$ polynomial basis to the linear advection equation and derive a PDE which the numerical solution solves exactly. We use a Fourier approach to derive polynomial solutions to this PDE and show that the polynomials are closely related to the $\frac{p}{p+1}$ Padé approximant of the exponential function. We show that for a uniform mesh of $N$ elements there exist $(p+1)N$ independent polynomial solutions, $N$ of which can be viewed as physical and $pN$ as non-physical. We show that the accumulation error of the physical mode is of order $2p+1$. In contrast, the non-physical modes are damped out exponentially quickly. We use these results to present a simple proof of the superconvergence of the DG method on uniform grids as well as show a connection between spatial superconvergence and the superaccuracies in dissipation and dispersion errors of the scheme. Finally, we show that for a class of initial projections on a uniform mesh, the superconvergent points of the numerical error tend exponentially quickly towards the downwind based Radau points.