Researcher profile

Irene Kyza

Irene Kyza contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

math.NA math.AP

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

2topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2016arXiv

$hp$-Adaptive Galerkin Time Stepping Methods for Nonlinear Initial Value Problems

This work is concerned with the derivation of an a posteriori error estimator for Galerkin approximations to nonlinear initial value problems with an emphasis on finite-time existence in the context of blow-up. The stucture of the derived estimator leads naturally to the development of both h and hp versions of an adaptive algorithm designed to approximate the blow-up time. The adaptive algorithms are then applied in a series of numerical experiments, and the rate of convergence to the blow-up time is investigated.

preprint2016arXiv

A posteriori error analysis for evolution nonlinear Schrödinger equations up to the critical exponent

We provide a posteriori error estimates in the $L^\infty(L^2)-$norm for relaxation time discrete and fully discrete schemes for a class of evolution nonlinear Schrödinger equations up to the critical exponent. In particular for the discretization in time we use the relaxation Crank-Nicolson-type scheme introduced by Besse in \cite{Besse}. For the discretization in space we use finite element spaces that are allowed to change between time steps. The estimates are obtained using the reconstruction technique. Through this technique the problem is converted to a perturbation of the original partial differential equation and this makes it possible to use nonlinear stability arguments as in the continuous problem. In particular, main ingredients we use in our analysis are the Gagliardo-Nirenberg inequality and the two conservation laws (mass and energy conservation) of the continuous problem. Numerical results illustrate that the estimates are indeed of optimal order of convergence.

preprint2015arXiv

Adaptivity and blow-up detection for nonlinear evolution problems

This work is concerned with the development of a space-time adaptive numerical method, based on a rigorous a posteriori error bound, for a semilinear convection-diffusion problem which may exhibit blow-up in finite time. More specifically, a posteriori error bounds are derived in the $L^{\infty}(L^2)+L^2(H^1)$-type norm for a first order in time implicit-explicit (IMEX) interior penalty discontinuous Galerkin (dG) in space discretization of the problem, although the theory presented is directly applicable to the case of conforming finite element approximations in space. The choice of the discretization in time is made based on a careful analysis of adaptive time stepping methods for ODEs that exhibit finite time blow-up. The new adaptive algorithm is shown to accurately estimate the blow-up time of a number of problems, including one which exhibits regional blow-up.

preprint2015arXiv

Regularized semiclassical limits: linear flows with infinite Lyapunov exponents

Semiclassical asymptotics for linear Schrödinger equations with non-smooth potentials give rise to ill-posed formal semiclassical limits. These problems have attracted a lot of attention in the last few years, as a proxy for the treatment of eigenvalue crossings, i.e. general systems. It has recently been shown that the semiclassical limit for conical singularities is in fact well-posed, as long as the Wigner measure (WM) stays away from singular saddle points. In this work we develop a family of refined semiclassical estimates, and use them to derive regularized transport equations for saddle points with infinite Lyapunov exponents, extending the aforementioned recent results. In the process we answer a related question posed by P. L. Lions and T. Paul in 1993. If we consider more singular potentials, our rigorous estimates break down. To investigate whether conical saddle points, such as $-|x|$, admit a regularized transport asymptotic approximation, we employ a numerical solver based on posteriori error control. Thus rigorous upper bounds for the asymptotic error in concrete problems are generated. In particular, specific phenomena which render invalid any regularized transport for $-|x|$ are identified and quantified. In that sense our rigorous results are sharp. Finally, we use our findings to formulate a precise conjecture for the condition under which conical saddle points admit a regularized transport solution for the WM.

preprint2013arXiv

A posteriori error control & adaptivity for Crank-Nicolson finite element approximations for the linear Schrödinger equation

We derive optimal order a posteriori error estimates for fully discrete approximations of linear Schrödinger-type equations, in the $L^\infty(L^2)-$norm. For the discretization in time we use the Crank-Nicolson method, while for the space discretization we use finite element spaces that are allowed to change in time. The derivation of the estimators is based on a novel elliptic reconstruction that leads to estimates which reflect the physical properties of Schrödinger equations. The final estimates are obtained using energy techniques and residual-type estimators. Various numerical experiments for the one-dimensional linear Schrödinger equation in the semiclassical regime, verify and complement our theoretical results. The numerical implementations are performed with both uniform partitions and adaptivity in time and space. For adaptivity, we further develop and analyze an existing time-space adaptive algorithm to the cases of Schrödinger equations. The adaptive algorithm reduces the computational cost substantially and provides efficient error control for the solution and the observables of the problem, especially for small values of the Planck constant.