BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Kenneth Duru

Kenneth Duru contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.geo-phphysics.comp-phComputational Engineering, Finance, and ScienceDistributed, Parallel, and Cluster ComputingMathematical Softwarephysics.ao-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2026arXiv

On well-posed energy/entropy stable boundary conditions for the rotating shallow water equations

We derive and analyze well-posed, energy- and entropy-stable boundary conditions (BCs) for the two-dimensional linear and nonlinear rotating shallow water equations (RSWE) in vector invariant form. The focus of the study is on subcritical flows, which are commonly observed in atmospheric, oceanic, and geostrophic flow applications. We consider spatial domains with smooth boundaries and formulate both linear and nonlinear BCs using mass flux, Riemann's invariants, and Bernoulli's potential, ensuring that the resulting initial boundary value problem (IBVP) is provably entropy- and energy-stable. The linear analysis is comprehensive, providing sufficient conditions to establish the existence, uniqueness, and energy stability of solutions to the linear IBVP. For the nonlinear IBVP, which admits more general solutions, our goal is to develop nonlinear BCs that guarantee entropy stability. We introduce the concepts of linear consistency and linear stability for nonlinear IBVPs, demonstrating that if a nonlinear IBVP is both linearly consistent and linearly stable, then, for sufficiently regular initial and boundary data over a finite time interval, a unique smooth solution exists. Both the linear and nonlinear IBVPs can be efficiently solved using high-order accurate numerical methods. By employing high-order summation-by-parts operators to discretize spatial derivatives and implementing weak enforcement of BCs via penalty techniques, we develop provably energy- and entropy-stable numerical schemes on curvilinear meshes. Extensive numerical experiments are presented to verify the accuracy of the methods and to demonstrate the robustness of the proposed BCs and numerical schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accurate simulations of nonlinear dynamic shear ruptures on pre-existing faults in 3D elastic solids with dual-pairing SBP methods

In this paper we derive and analyse efficient and stable numerical methods for accurate numerical simulations of nonlinear dynamic shear ruptures on non-planar faults embedded in 3D elastic solids using dual-paring (DP) summation by parts (SBP) finite difference (FD) methods. Specifically, for nonlinear dynamic earthquake ruptures, we demonstrate that the DP SBP FD operators [K. Mattsson. J. Comput. Phys., 335:283-310, 2017] generate spurious catastrophic high frequency wave modes that do not diminish with mesh refinement. Meanwhile our new dispersion relation preserving (DRP) SBP FD operators [C. Williams and K Duru, arXiv:2110.04957, 2021] have more accurate numerical dispersion relation properties and do not support poisonous spurious high frequency wave modes. Numerical simulations are performed in 3D with geometrically complex fault surfaces verifying the efficacy of the method. Our method accurately reproduces community developed dynamic rupture benchmark problems, proposed by Southern California Earthquake Center, with less computational effort than standard methods based on traditional SBP FD operators.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On energy-stable and high order finite element methods for the wave equation in heterogeneous media with perfectly matched layers

This paper presents a stable finite element approximation for the acoustic wave equation on second-order form, with perfectly matched layers (PML) at the boundaries. Energy estimates are derived for varying PML damping for both the discrete and the continuous case. Moreover, a priori error estimates are derived for constant PML damping. Most of the analysis is performed in Laplace space. Numerical experiments in physical space validate the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The perfectly matched layer (PML) for hyperbolic wave propagation problems: A review

It is well-known that reliable and efficient domain truncation is crucial to accurate numerical solution of most wave propagation problems. The perfectly matched layer (PML) is a method which, when stable, can provide a domain truncation scheme which is convergent with increasing layer width/damping. The difficulties in using the PML are primarily associated with stability, which can be present at the continuous level or be triggered by numerical approximations. The mathematical and numerical analysis of the PML for hyperbolic wave propagation problems has been an area of active research. It is now possible to construct stable and high order accurate numerical wave solvers by augmenting wave equations with the PML and approximating the equations using summation-by-parts finite difference methods, continuous and discontinuous Galerkin finite element methods. In this review we summarise the progress made, from mathematical, numerical and practical perspectives, point out some open problems and set the stage for future work. We also present numerical experiments of model problems corroborating the theoretical analysis, and numerical simulations of real-world wave propagation demonstrating impact. Stable and parallel implementations of the PML in the high performance computing software packages WaveQLab3D and ExaHyPE allow to sufficiently limit the computational domain of seismological problems with only a few grid points/elements around the computational boundaries where the PML is active, thus saving as much as $96\%$ of the required computational resources for a three space dimensional seismological benchmark problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form

We present a stable discontinuous Galerkin (DG) method with a perfectly matched layer (PML) for three and two space dimensional linear elastodynamics, in velocity-stress formulation, subject to well-posed linear boundary conditions. First, we consider the elastodynamics equation, in a cuboidal domain, and derive an unsplit PML truncating the domain using complex coordinate stretching. Leveraging the hyperbolic structure of the underlying system, we construct continuous energy estimates, in the time domain for the elastic wave equation, and in the Laplace space for a sequence of PML model problems, with variations in one, two and three space dimensions, respectively. They correspond to PMLs normal to boundary faces, along edges and in corners. Second, we develop a DG numerical method for the linear elastodynamics equation using physically motivated numerical flux and penalty parameters, which are compatible with all well-posed, internal and external, boundary conditions. When the PML damping vanishes, by construction, our choice of penalty parameters yield an upwind scheme and a discrete energy estimate analogous to the continuous energy estimate. Third, to ensure numerical stability of the discretization when PML damping is present, it is necessary to extend the numerical DG fluxes, and the numerical inter-element and boundary procedures, to the PML auxiliary differential equations. This is crucial for deriving discrete energy estimates analogous to the continuous energy estimates. By combining the DG spatial approximation with the high order ADER time stepping scheme and the accuracy of the PML we obtain an arbitrarily accurate wave propagation solver in the time domain. Numerical experiments are presented in two and three space dimensions corroborating the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are useful in a wide range of application areas. Applications powered by ExaHyPE can be run on a student's laptop, but are also able to exploit thousands of processor cores on state-of-the-art supercomputers. The engine is able to dynamically increase the accuracy of the simulation using adaptive mesh refinement where required. Due to the robustness and shock capturing abilities of ExaHyPE's numerical methods, users of the engine can simulate linear and non-linear hyperbolic PDEs with very high accuracy. Users can tailor the engine to their particular PDE by specifying evolved quantities, fluxes, and source terms. A complete simulation code for a new hyperbolic PDE can often be realised within a few hours - a task that, traditionally, can take weeks, months, often years for researchers starting from scratch. In this paper, we showcase ExaHyPE's workflow and capabilities through real-world scenarios from our two main application areas: seismology and astrophysics.

0 citations0 reviewsNo code linkedOpen access