BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Lauri Oksanen

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

ResearcherUnclaimed source record
math.APmath.DGmath.NAmath.OCmath.PRNumerical Analysismath-phmath.MPmath.STStatistics Theory

Catalog footprint

What is connected

27works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

27 published item(s)

preprint2026arXiv

Recovering a (1+1)-dimensional wave equation from a single white noise boundary measurement

We consider the following inverse problem: Suppose a $(1+1)$-dimensional wave equation on $\mathbb{R}_+$ with zero initial conditions is excited with a Neumann boundary data modelled as a white noise process. Given also the Dirichlet data at the same point, determine the unknown first order coefficient function of the system. We first establish that direct problem is well-posed. The inverse problem is then solved by showing that correlations of the boundary data determine the Neumann-to-Dirichlet operator in the sense of distributions, which is known to uniquely identify the coefficient. This approach has applications in acoustic measurements of internal cross-sections of fluid pipes such as pressurised water supply pipes and vocal tract shape determination.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A stabilized finite element method for inverse problems subject to the convection-diffusion equation. II: convection-dominated regime

We consider the numerical approximation of the ill-posed data assimilation problem for stationary convection-diffusion equations and extend our previous analysis in [Numer. Math. 144, 451--477, 2020] to the convection-dominated regime. Slightly adjusting the stabilized finite element method proposed for dominant diffusion, we draw upon a local error analysis to obtain quasi-optimal convergence along the characteristics of the convective field through the data set. The weight function multiplying the discrete solution is taken to be Lipschitz and a corresponding super approximation result (discrete commutator property) is proven. The effect of data perturbations is included in the analysis and we conclude the paper with some numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Counterexamples to inverse problems for the wave equation

We construct counterexamples to inverse problems for the wave operator on domains in $\mathbb{R}^{n+1}$, $n \ge 2$, and on Lorentzian manifolds. We show that non-isometric Lorentzian metrics can lead to same partial data measurements, which are formulated in terms certain restrictions of the Dirichlet-to-Neumann map. The Lorentzian metrics giving counterexamples are time-dependent, but they are smooth and non-degenerate. On $\mathbb{R}^{n+1}$ the metrics are conformal to the Minkowski metric.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Construction of artificial point sources for a linear wave equation in unknown medium

We study the wave equation on a bounded domain of $\mathbb R^m$ and on a compact Riemannian manifold $M$ with boundary. We assume that the coefficients of the wave equation are unknown but that we are given the hyperbolic Neumann-to-Dirichlet map $Λ$ that corresponds to the physical measurements on the boundary. Using the knowledge of $Λ$ we construct a sequence of Neumann boundary values so that at a time $T$ the corresponding waves converge to zero while the time derivative of the waves converge to a delta distribution. Such waves are called an artificial point source. The convergence of the wave takes place in the function spaces naturally related to the energy of the wave. We apply the results for inverse problems and demonstrate the focusing of the waves numerically in the 1-dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inverse problems for real principal type operators

We consider inverse boundary value problems for general real principal type differential operators. The first results state that the Cauchy data set uniquely determines the scattering relation of the operator and bicharacteristic ray transforms of lower order coefficients. We also give two different boundary determination methods for general operators, and prove global uniqueness results for determining coefficients in nonlinear real principal type equations. The article presents a unified approach for treating inverse boundary problems for transport and wave equations, and highlights the role of propagation of singularities in the solution of related inverse problems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A stabilized finite element method for inverse problems subject to the convection-diffusion equation. I: diffusion-dominated regime

The numerical approximation of an inverse problem subject to the convection--diffusion equation when diffusion dominates is studied. We derive Carleman estimates that are on a form suitable for use in numerical analysis and with explicit dependence on the Péclet number. A stabilized finite element method is then proposed and analysed. An upper bound on the condition number is first derived. Combining the stability estimates on the continuous problem with the numerical stability of the method, we then obtain error estimates in local $H^1$- or $L^2$-norms that are optimal with respect to the approximation order, the problem's stability and perturbations in data. The convergence order is the same for both norms, but the $H^1$-estimate requires an additional divergence assumption for the convective field. The theory is illustrated in some computational examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Correlation based passive imaging with a white noise source

Passive imaging refers to problems where waves generated by unknown sources are recorded and used to image the medium through which they travel. The sources are typically modelled as a random variable and it is assumed that some statistical information is available. In this paper we study the stochastic wave equation $\partial_t^2 u - Δ_g u = χW$, where $W$ is a random variable with the white noise statistics on ${\mathbb R}^{1+n}$, $n \ge 3$, $χ$ is a smooth function vanishing for negative times and outside a compact set in space, and $Δ_g$ is the Laplace-Beltrami operator associated to a smooth non-trapping Riemannian metric tensor $g$ on ${\mathbb R}^n$. The metric tensor $g$ models the medium to be imaged, and we assume that it coincides with the Euclidean metric outside a compact set. We consider the empirical correlations on an open set $\mathcal X \subset {\mathbb R}^n$, $$ C_T(t_1, x_1, t_2, x_2) = \frac 1 T \int_0^T u(t_1+s,x_1) u(t_2+s,x_2) ds, \quad t_1,t_2>0,\ x_1,x_2\in \mathcal X, $$ for $T>0$. Supposing that $χ$ is non-zero on $\mathcal X$ and constant in time after $t > 1$, we show that in the limit $T \to \infty$, the data $C_T$ becomes statistically stable, that is, independent of the realization of $W$. Our main result is that, with probability one, this limit determines the Riemannian manifold $({\mathbb R}^n,g)$ up to an isometry. To our knowledge, this is the first result showing that a medium can be determined in a passive imaging setting, without assuming a separation of scales.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Data assimilation for the heat equation using stabilized finite element methods

We consider data assimilation for the heat equation using a finite element space semi-discretization. The approach is optimization based, but the design of regularization operators and parameters rely on techniques from the theory of stabilized finite elements. The space semi-discretized system is shown to admit a unique solution. Combining sharp estimates of the numerical stability of the discrete scheme and conditional stability estimates of the ill-posed continuous pde-model we then derive error estimates that reflect the approximation order of the finite element space and the stability of the continuous model. Two different data assimilation situations with different stability properties are considered to illustrate the framework. Full detail on how to adapt known stability estimates for the continuous model to work with the numerical analysis framework is given in appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Global uniqueness in an inverse problem for time fractional diffusion equations

Given $(M,g)$, a compact connected Riemannian manifold of dimension $d \geq 2$, with boundary $\partial M$, we consider an initial boundary value problem for a fractional diffusion equation on $(0,T) \times M$, $T>0$, with time-fractional Caputo derivative of order $α\in (0,1) \cup (1,2)$. We prove uniqueness in the inverse problem of determining the smooth manifold $(M,g)$ (up to an isometry), and various time-independent smooth coefficients appearing in this equation, from measurements of the solution on a subset of $\partial M$ at fixed time. In the "flat" case where $M$ is a compact subset of $\mathbb R^d$, two out the three coefficients $ρ$ (weight), $a$ (conductivity) and $q$ (potential) appearing in the equation $ρ\partial_t^αu-\textrm{div}(a \nabla u)+ q u=0$ on $(0,T)\times Ω$ are recovered simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Recovery of time-dependent coefficient on Riemanian manifold for hyperbolic equations

Given $(M,g)$, a compact connected Riemannian manifold of dimension $d \geq 2$, with boundary $\partial M$, we study the inverse boundary value problem of determining a time-dependent potential $q$, appearing in the wave equation $\partial_t^2u-Δ_g u+q(t,x)u=0$ in $\bar M=(0,T)\times M$ with $T>0$. Under suitable geometric assumptions we prove global unique determination of $q\in L^\infty(\bar M)$ given the Cauchy data set on the whole boundary $\partial \bar M$, or on certain subsets of $\partial \bar M$. Our problem can be seen as an analogue of the Calderón problem on the Lorentzian manifold $(\bar M, dt^2 - g)$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Time reversal method with stabilizing boundary conditions for Photoacoustic tomography

We study an inverse initial source problem that models Photoacoustic tomography measurements with array detectors, and introduce a method that can be viewed as a modification of the so called back and forth nudging method. We show that the method converges at an exponential rate under a natural visibility condition, with data given only on a part of the boundary of the domain of wave propagation. In this paper we consider the case of noiseless measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Determination of the Spacetime from Local Time Measurements

We consider an inverse problem for a Lorentzian spacetime $(M,g)$, and show that time measurements, that is, the knowledge of the Lorentzian time separation function on a submanifold $Σ\subset M$ determine the $C^\infty$-jet of the metric in the Fermi coordinates associated to $Σ$. We use this result to study the global determination of the spacetime $(M,g)$ when it has a real-analytic structure or is stationary and satisfies the Einstein-scalar field equations. In addition to this, we require that $(M,g)$ is geodesically complete modulo scalar curvature singularities. The results are Lorentzian counterparts of extensively studied inverse problems in Riemannian geometry - the determination of the jet of the metric and the boundary rigidity problem. We give also counterexamples in cases when the assumptions are not valid, and discuss inverse problems in general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Construction of Virtual Interior Point Source Travel Time Distances from the Hyperbolic Neumann-to-Dirichlet Map

We introduce a new algorithm to construct travel time distances between a point in the interior of a Riemannian manifold and points on the boundary of the manifold, and describe a numerical implementation of the algorithm. It is known that the travel time distances for all interior points determine the Riemannian manifold in a stable manner. We do not assume that there are sources or receivers in the interior, and use the hyperbolic Neumann-to-Dirichlet map, or its restriction, as our data. Our algorithm is a variant of the Boundary Control method, and to our knowledge, this is the first numerical implementation of the method in a geometric setting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Regularization strategy for inverse problem for 1+1 dimensional wave equation

An inverse boundary value problem for a 1+1 dimensional wave equation with wave speed $c(x)$ is considered. We give a regularisation strategy for inverting the map $\mathcal A:c\mapsto Λ,$ where $Λ$ is the hyperbolic Neumann-to-Dirichlet map corresponding to the wave speed $c$. More precisely, we consider the case when we are given a perturbation of the Neumann-to-Dirichlet map $\tilde Λ=Λ+\mathcal E $, where $\mathcal E$ corresponds to the measurement errors, and reconstruct an approximate wave speed $\tilde c$. We emphasize that $\tilde Λ$ may not not be in the range of the map $\mathcal A$. We show that the reconstructed wave speed $\tilde c$ satisfies $\| \tilde c-c\|_{L^\infty}<C \|E\|^{1/18}$. Our regularization strategy is based on a new formula to compute $c$ from $Λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Uniqueness for a seismic inverse source problem modeling a subsonic rupture

We consider an inverse problem for an inhomogeneous wave equation with discrete-in-time sources, modeling a seismic rupture. We assume that the sources occur along a path with subsonic velocity, and that data are collected over time on some detection surface. We explore the question of uniqueness for these problems, show how to recover the times and locations of sources microlocally, and then reconstruct the smooth part of the source assuming that it is the same at each source location.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inverse problem for the wave equation with a white noise source

We consider a smooth Riemannian metric tensor $g$ on $\R^n$ and study the stochastic wave equation for the Laplace-Beltrami operator $\p_t^2 u - Δ_g u = F$. Here, $F=F(t,x,ω)$ is a random source that has white noise distribution supported on the boundary of some smooth compact domain $M \subset \R^n$. We study the following formally posed inverse problem with only one measurement. Suppose that $g$ is known only outside of a compact subset of $M^{int}$ and that a solution $u(t,x,ω_0)$ is produced by a single realization of the source $F(t,x,ω_0)$. We ask what information regarding $g$ can be recovered by measuring $u(t,x,ω_0)$ on $\R_+ \times \p M$? We prove that such measurement together with the realization of the source determine the scattering relation of the Riemannian manifold $(M, g)$ with probability one. That is, for all geodesics passing through $M$, the travel times together with the entering and exit points and directions are determined. In particular, if $(M,g)$ is a simple Riemannian manifold and $g$ is conformally Euclidian in $M$, the measurement determines the metric $g$ in $M$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Photoacoustic and thermoacoustic tomography with an uncertain wave speed

We consider the mathematical model of photoacoustic and thermoacoustic tomography in media with a variable sound speed. When the sound speed is known, the explicit reconstruction formula by P. Stefanov and G. Uhlmann (Inverse Problems, 25(7):075011, 16, 2009) can be used. We study how a modelling error in the sound speed affects the reconstruction formula and quantify the effect in terms of a stability estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Lipschitz stable reconstruction formula for the inverse problem for the wave equation

We consider the problem to reconstruct a wave speed $c \in C^\infty(M)$ in a domain $M \subset \R^n$ from acoustic boundary measurements modelled by the hyperbolic Dirichlet-to-Neumann map $Λ$. We introduce a reconstruction formula for $c$ that is based on the Boundary Control method and incorporates features also from the complex geometric optics solutions approach. Moreover, we show that the reconstruction formula is locally Lipschitz stable for a low frequency component of $c^{-2}$ under the assumption that the Riemannian manifold $(M, c^{-2} dx^2)$ has a strictly convex function with no critical points. That is, we show that for all bounded $C^2$ neighborhoods $U$ of $c$, there is a $C^1$ neighborhood $V$ of $c$ and constants $C, R > 0$ such that |\F\ll(\tilde c^{-2} - c^{-2}\rr)(ξ)| \le C e^{2R |ξ|} \norm{\tilde Λ- Λ}_*, \quad ξ\in \R^n, for all $\tilde c \in U \cap V$, where $\tilde Λ$ is the Dirichlet-to-Neumann map corresponding to the wave speed $\tilde c$ and $\norm{\cdot}_*$ is a norm capturing certain regularity properties of the Dirichlet-to-Neumann maps.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inverse problem for the Riemannian wave equation with Dirichlet data and Neumann data on disjoint sets

We consider the inverse problem to determine a smooth compact Riemannian manifold with boundary $(M, g)$ from a restriction $Λ_{\Src, \Rec}$ of the Dirichlet-to-Neumann operator for the wave equation on the manifold. Here $\Src$ and $\Rec$ are open sets in $\p M$ and the restriction $Λ_{\Src, \Rec}$ corresponds to the case where the Dirichlet data is supported on $\R_+\times \Src$ and the Neumann data is measured on $\R_+\times \Rec$. In the novel case where $\bar \Src \cap \bar \Rec = \emptyset$, we show that $Λ_{\Src, \Rec}$ determines the manifold $(M,g)$ uniquely, assuming that the wave equation is exactly controllable from the set of sources $\Src$. Moreover, we show that the exact controllability can be replaced by the Hassell-Tao condition for eigenvalues and eigenfunctions, that is, λ_j \le C \norm{\p_νϕ_j}_{L^2(\Src)}^2, \quad j =1, 2, ..., where $λ_j$ are the Dirichlet eigenvalues and $(ϕ_j)_{j=1}^\infty$ is an orthonormal basis of the corresponding eigenfunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving an inverse obstacle problem for the wave equation by using the boundary control method

We introduced in [arXiv:1106.3204] a method to locate discontinuities of a wave speed in dimension two from acoustic boundary measuments modelled by the hyperbolic Neumann-to-Dirichlet operator. Here we extend the method for sound hard obstacles in arbitrary dimension. We present numerical experiments with simulated noisy data suggesting that the method is robust against measurement noise.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Inverse obstacle problem for the non-stationary wave equation with an unknown background

We consider boundary measurements for the wave equation on a bounded domain $M \subset \R^2$ or on a compact Riemannian surface, and introduce a method to locate a discontinuity in the wave speed. Assuming that the wave speed consist of an inclusion in a known smooth background, the method can determine the distance from any boundary point to the inclusion. In the case of a known constant background wave speed, the method reconstructs a set contained in the convex hull of the inclusion and containing the inclusion. Even if the background wave speed is unknown, the method can reconstruct the distance from each boundary point to the inclusion assuming that the Riemannian metric tensor determined by the wave speed gives simple geometry in $M$. The method is based on reconstruction of volumes of domains of influence by solving a sequence of linear equations. For $τ\in C(\p M)$ the domain of influence $M(τ)$ is the set of those points on the manifold from which the distance to some boundary point $x$ is less than $τ(x)$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Solving an inverse problem for the wave equation by using a minimization algorithm and time-reversed measurements

We consider the inverse problem for the wave equation on a compact Riemannian manifold or on a bounded domain of $\R^n$, and generalize the concept of {\em domain of influence}. We present an efficient minimization algorithm to compute the volume of a domain of influence using boundary measurements and time-reversed boundary measurements. Moreover, we show that if the manifold is simple, then the volumes of the domains of influence determine the manifold. For a continuous real valued function $τ$ on the boundary of the manifold, the domain of influence is the set of those points on the manifold from which the travel time to some boundary point $y$ is less than $τ(y)$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An inverse problem for the wave equation with one measurement and the pseudorandom noise

We consider the wave equation $(\p_t^2-Δ_g)u(t,x)=f(t,x)$, in $\R^n$, $u|_{\R_-\times \R^n}=0$, where the metric $g=(g_{jk}(x))_{j,k=1}^n$ is known outside an open and bounded set $M\subset \R^n$ with smooth boundary $\p M$. We define a deterministic source $f(t,x)$ called the pseudorandom noise as a sum of point sources, $f(t,x)=\sum_{j=1}^\infty a_jδ_{x_j}(x)δ(t)$, where the points $x_j,\ j\in\Z_+$, form a dense set on $\p M$. We show that when the weights $a_j$ are chosen appropriately, $u|_{\R\times \p M}$ determines the scattering relation on $\p M$, that is, it determines for all geodesics which pass through $M$ the travel times together with the entering and exit points and directions. The wave $u(t,x)$ contains the singularities produced by all point sources, but when $a_j=λ^{-λ^{j}}$ for some $λ>1$, we can trace back the point source that produced a given singularity in the data. This gives us the distance in $(\R^n, g)$ between a source point $x_j$ and an arbitrary point $y \in \p M$. In particular, if $(\bar M,g)$ is a simple Riemannian manifold and $g$ is conformally Euclidian in $\bar M$, these distances are known to determine the metric $g$ in $M$. In the case when $(\bar M,g)$ is non-simple we present a more detailed analysis of the wave fronts yielding the scattering relation on $\p M$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inverse problem for wave equation with sources and observations on disjoint sets

We consider an inverse problem for a hyperbolic partial differential equation on a compact Riemannian manifold. Assuming that $Γ_1$ and $Γ_2$ are two disjoint open subsets of the boundary of the manifold we define the restricted Dirichlet-to-Neumann operator $Λ_{Γ_1,Γ_2}$. This operator corresponds the boundary measurements when we have smooth sources supported on $Γ_1$ and the fields produced by these sources are observed on $Γ_2$. We show that when $Γ_1$ and $Γ_2$ are disjoint but their closures intersect at least at one point, then the restricted Dirichlet-to-Neumann operator $Λ_{Γ_1,Γ_2}$ determines the Riemannian manifold and the metric on it up to an isometry. In the Euclidian space, the result yields that an anisotropic wave speed inside a compact body is determined, up to a natural coordinate transformations, by measurements on the boundary of the body even when wave sources are kept away from receivers. Moreover, we show that if we have three arbitrary non-empty open subsets $Γ_1,Γ_2$, and $Γ_3$ of the boundary, then the restricted Dirichlet-to-Neumann operators $Λ_{Γ_j,Γ_k}$ for $1\leq j<k\leq 3$ determine the Riemannian manifold to an isometry. Similar result is proven also for the finite-time boundary measurements when the hyperbolic equation satisfies an exact controllability condition.

0 citations0 reviewsNo code linkedOpen access