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

Ilaria Perugia

Ilaria Perugia contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
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

9 published item(s)

preprint2026arXiv

Stability, convergence, and geometric properties of second-order-in-time space-time discretizations for linear and semilinear wave equations

We revisit second-order-in-time space-time discretizations of the linear and semilinear wave equations by establishing precise equivalences with first-order-in-time formulations. Focusing on schemes using continuous piecewise-polynomial trial functions in time, we analyze their stability, convergence, and geometric properties. We consider first a weak space-time formulation with test functions projected onto discontinuous polynomials of one degree lower in time, showing that it is equivalent to the scheme proposed in [French, Peterson 1996] in the linear case, and extended in [Karakashian, Makridakis 2005] to the semilinear case. In particular, this equivalence shows that this method conserves energy at mesh nodes but is not symplectic. We then introduce two symplectic variants, obtained through Gauss-Legendre and Gauss-Lobatto quadratures in time, and show that they correspond to specific Runge-Kutta time integrators. These connections clarify the geometric structure of the space-time methods considered.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An adaptive finite element method for high-frequency scattering problems with smoothly varying coefficients

We introduce a new numerical method for solving time-harmonic acoustic scattering problems. The main focus is on plane waves scattered by smoothly varying material inhomogeneities. The proposed method works for any frequency $ω$, but is especially efficient for high-frequency problems. It is based on a time-domain approach and consists of three steps: \emph{i)} computation of a suitable incoming plane wavelet with compact support in the propagation direction; \emph{ii)} solving a scattering problem in the time domain for the incoming plane wavelet; \emph{iii)} reconstruction of the time-harmonic solution from the time-domain solution via a Fourier transform in time. An essential ingredient of the new method is a front-tracking mesh adaptation algorithm for solving the problem in \emph{ii)}. By exploiting the limited support of the wave front, this allows us to make the number of the required degrees of freedom to reach a given accuracy significantly less dependent on the frequency $ω$. We also present a new algorithm for computing the Fourier transform in \emph{iii)} that exploits the reduced number of degrees of freedom corresponding to the adapted meshes. Numerical examples demonstrate the advantages of the proposed method and the fact that the method can also be applied with external source terms such as point sources and sound-soft scatterers. The gained efficiency, however, is limited in the presence of trapping modes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An entropy structure preserving space-time formulation for cross-diffusion systems: Analysis and Galerkin discretization

Cross-diffusion systems are systems of nonlinear parabolic partial differential equations that are used to describe dynamical processes in several application, including chemical concentrations and cell biology. We present a space-time approach to the proof of existence of bounded weak solutions of cross-diffusion systems, making use of the system entropy to examine long-term behavior and to show that the solution is nonnegative, even when a maximum principle is not available. This approach naturally gives rise to a novel space-time Galerkin method for the numerical approximation of cross-diffusion systems that conserves their entropy structure. We prove existence and convergence of the discrete solutions, and present numerical results for the porous medium, the Fisher-KPP, and the Maxwell-Stefan problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the exponential time-decay for the one-dimensional wave equation with variable coefficients

We consider the initial-value problem for the one-dimensional, time-dependent wave equation with positive, Lipschitz continuous coefficients, which are constant outside a bounded region. Under the assumption of compact support of the initial data, we prove that the local energy decays exponentially fast in time, and provide the explicit constant to which the solution converges for large times. We give explicit estimates of the rate of this exponential decay by two different techniques. The first one is based on the definition of a modified, weighted local energy, with suitably constructed weights. The second one is based on the integral formulation of the problem and, under a more restrictive assumption on the variation of the coefficients, allows us to obtain improved decay rates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A numerical study of the dispersion and dissipation properties of virtual element methods for the Helmholtz problem

We study numerically the dispersion and dissipation properties of the plane wave virtual element method and the nonconforming Trefftz virtual element method for the Helmholtz problem. Whereas the former method is based on a conforming virtual partition of unity approach in the sense that the local (implicitly defined) basis functions are given as modulations of lowest order harmonic virtual element functions with plane waves, the latter one represents a pure Trefftz method with local edge-related basis functions that are eventually glued together in a nonconforming fashion. We will see that the qualitative and quantitative behavior of dissipation and dispersion of the method hinges upon the level of conformity and the use of Trefftz basis functions. To this purpose, we also compare the results to those obtained in [15] for the plane wave discontinuous Galerkin method, and to those for the standard polynomial based finite element method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Space-time discontinuous Galerkin approximation of acoustic waves with point singularities

We develop a convergence theory of space-time discretizations for the linear, 2nd-order wave equation in polygonal domains $Ω\subset\mathbb{R}^2$, possibly occupied by piecewise homogeneous media with different propagation speeds. Building on an unconditionally stable space-time DG formulation developed in [Moiola, Perugia 2018], we (a) prove optimal convergence rates for the space-time scheme with local isotropic corner mesh refinement on the spatial domain, and (b) demonstrate numerically optimal convergence rates of a suitable \emph{sparse} space-time version of the DG scheme. The latter scheme is based on the so-called \emph{combination formula}, in conjunction with a family of anisotropic space-time DG-discretizations. It results in optimal-order convergent schemes, also in domains with corners, with a number of degrees of freedom that scales essentially like the DG solution of one stationary elliptic problem in $Ω$ on the finest spatial grid. Numerical experiments for both smooth and singular solutions support convergence rate optimality on spatially refined meshes of the full and sparse space-time DG schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fast Least-Squares Padé approximation of problems with normal operators and meromorphic structure

In this work, we consider the approximation of Hilbert space-valued meromorphic functions that arise as solution maps of parametric PDEs whose operator is the shift of an operator with normal and compact resolvent, e.g. the Helmholtz equation. In this restrictive setting, we propose a simplified version of the Least-Squares Padé approximation technique introduced in [6] following [11]. In particular, the estimation of the poles of the target function reduces to a low-dimensional eigenproblem for a Gramian matrix, allowing for a robust and efficient numerical implementation (hence the "fast" in the name). Moreover, we prove several theoretical results that improve and extend those in [6], including the exponential decay of the error in the approximation of the poles, and the convergence in measure of the approximant to the target function. The latter result extends the classical one for scalar Padé approximation to our functional framework. We provide numerical results that confirm the improved accuracy of the proposed method with respect to the one introduced in [6] for differential operators with normal and compact resolvent.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A space-time Trefftz discontinuous Galerkin method for the acoustic wave equation in first-order formulation

We introduce a space-time Trefftz discontinuous Galerkin method for the first-order transient acoustic wave equations in arbitrary space dimensions, extending the one dimensional scheme of Kretzschmar et al. (2016, IMA J. Numer. Anal., 36, 1599-1635). Test and trial discrete functions are space-time piecewise polynomial solutions of the wave equations. We prove well-posedness and a priori error bounds in both skeleton-based and mesh-independent norms. The space-time formulation corresponds to an implicit time-stepping scheme, if posed on meshes partitioned in time slabs, or to an explicit scheme, if posed on "tent-pitched" meshes. We describe two Trefftz polynomial discrete spaces, introduce bases for them and prove optimal, high-order $h$-convergence bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A priori error analysis of space-time Trefftz discontinuous Galerkin methods for wave problems

We present and analyse a space-time discontinuous Galerkin method for wave propagation problems. The special feature of the scheme is that it is a Trefftz method, namely that trial and test functions are solution of the partial differential equation to be discretised in each element of the (space-time) mesh. The method considered is a modification of the discontinuous Galerkin schemes of Kretzschmar et al., and of Monk and Richter. For Maxwell's equations in one space dimension, we prove stability of the method, quasi-optimality, best approximation estimates for polynomial Trefftz spaces and (fully explicit) error bounds with high order in the meshwidth and in the polynomial degree. The analysis framework also applies to scalar wave problems and Maxwell's equations in higher space dimensions. Some numerical experiments demonstrate the theoretical results proved and the faster convergence compared to the non-Trefftz version of the scheme.

0 citations0 reviewsNo code linkedOpen access