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

Sebastian Acosta

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

ResearcherUnclaimed source record
math.APmath.NAphysics.comp-phNumerical Analysismath-phmath.MPQuantitative Methodsmath.OC

Catalog footprint

What is connected

19works
8topics
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

19 published item(s)

preprint2022arXiv

Nonlinear ultrasound imaging modeled by a Westervelt equation

We consider the ultrasound imaging problem governed by a nonlinear wave equation of Westervelt type with variable wave speed. We show that the coefficient of nonlinearity can be recovered uniquely from knowledge of the Dirichlet-to-Neumann map. Our proof is based on a second order linearization and the use of Gaussian beam solutions to reduce the problem to the inversion of a weighted geodesic ray transform. We propose an inversion algorithm and report the results of a numerical implementation to solve the nonlinear ultrasound imaging problem in a transmission setting in the frequency domain.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Local on-surface radiation condition for multiple scattering of waves from convex obstacles

We propose a novel on-surface radiation condition to approximate the outgoing solution to the Helmholtz equation in the exterior of several impenetrable convex obstacles. Based on a local approximation of the Dirichlet-to-Neumann operator and a local formula for wave propagation, this new method simultaneously accounts for the outgoing behavior of the solution as well as the reflections arising from the multiple obstacles. The method involves tangential derivatives only, avoiding the use of integration over the surfaces of the obstacles. As a consequence, the method leads to sparse matrices and O(N) complexity. Numerical results are presented to illustrate the performance of the proposed method. Possible improvements and extensions are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A weight-adjusted discontinuous Galerkin method for wave propagation in coupled elastic-acoustic media

This paper presents a high-order discontinuous Galerkin (DG) scheme for the simulation of wave propagation through coupled elastic-acoustic media. We use a first-order stress-velocity formulation, and derive a simple upwind-like numerical flux which weakly imposes continuity of the normal velocity and traction at elastic-acoustic interfaces. When combined with easily invertible weight-adjusted mass matrices, the resulting method is efficient, consistent, and energy stable on curvilinear meshes and for arbitrary heterogeneous media, including anisotropy and sub-cell (micro) heterogeneities. We numerically verify the high order accuracy and stability of the proposed method, and investigate its performance for applications in photoacoustic tomography.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High order methods for acoustic scattering: Coupling Farfield Expansions ABC with Deferred-Correction methods

Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs, with the appropriate number of terms, to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the novel method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solvability for Photoacoustic Imaging with Idealized Piezoelectric Sensors

Most reconstruction algorithms for photoacoustic imaging assume that the pressure field is measured by ultrasound sensors placed on a detection surface. However, such sensors do not measure pressure exactly due to their non-uniform directional and frequency responses, and resolution limitations. This is the case for piezoelectric sensors that are commonly employed for photoacoustic imaging. In this paper, using the method of matched asymptotic expansions and the basic constitutive relations for piezoelectricity, we propose a simple mathematical model for piezoelectric transducers. The approach simultaneously models how the pressure waves induce the piezoelectric measurements and how the presence of the sensors affects the pressure waves. Using this model, we analyze whether the data gathered by piezoelectric sensors leads to the mathematical solvability of the photoacoustic imaging problem. We conclude that this imaging problem is well-posed in certain normed spaces and under a geometric assumption. We also propose an iterative reconstruction algorithm that incorporates the model for piezoelectric measurements. A numerical implementation of the reconstruction algorithm is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Recovery of Pressure and Wave Speed for Photoacoustic Imaging under a Condition of Relative Uncertainty

In this paper, we study the photoacoustic tomography problem for which we seek to recover both the initial state of the pressure field and the wave speed of the medium from the knowledge of a single boundary measurement. The goal is to propose practical assumptions to define a set of initial conditions and wave speeds over which uniqueness for this inverse problem is guaranteed. The main result of the paper is that given two sets of wave speeds and pressure profiles, they cannot produce the same acoustic measurements if the relative difference between the wave speeds is much smaller than the relative difference between the pressure profiles. Implications for iterative joint-reconstruction algorithms are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Well-posedness for Photoacoustic Tomography with Fabry-Perot Sensors

In the mathematical analysis of photoacoustic imaging, it is usually assumed that the acoustic pressure (Dirichlet data) is measured on a detection surface. However, actual ultrasound detectors gather data of a different type. In this paper, we propose a more realistic mathematical model of ultrasound measurements acquired by the Fabry--Perot sensor. This modeling incorporates directional response of such sensors. We study the solvability of the resulting photoacoustic tomography problem, concluding that the problem is well-posed under certain assumptions. Numerical reconstructions are implemented using the Landweber iterations, after discretization of the governing equations using the finite element method.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

High order surface radiation conditions for time-harmonic waves in exterior domains

We formulate a new family of high order on-surface radiation conditions to approximate the outgoing solution to the Helmholtz equation in exterior domains. Motivated by the pseudo-differential expansion of the Dirichlet-to-Neumann operator developed by Antoine et al. (J. Math. Anal. Appl. 229:184-211, 1999), we design a systematic procedure to apply pseudo-differential symbols of arbitrarily high order. Numerical results are presented to illustrate the performance of the proposed method for solving both the Dirichlet and the Neumann boundary value problems. Possible improvements and extensions are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High order local absorbing boundary conditions for acoustic waves in terms of farfield expansions

We devise a new high order local absorbing boundary condition (ABC) for radiating problems and scattering of time-harmonic acoustic waves from obstacles of arbitrary shape. By introducing an artificial boundary $S$ enclosing the scatterer, the original unbounded domain $Ω$ is decomposed into a bounded computational domain $Ω^{-}$ and an exterior unbounded domain $Ω^{+}$. Then, we define interface conditions at the artificial boundary $S$, from truncated versions of the well-known Wilcox and Karp farfield expansion representations of the exact solution in the exterior region $Ω^{+}$. As a result, we obtain a new local absorbing boundary condition (ABC) for a bounded problem on $Ω^{-}$, which effectively accounts for the outgoing behavior of the scattered field. Contrary to the low order absorbing conditions previously defined, the order of the error induced by this ABC can easily match the order of the numerical method in $Ω^{-}$. We accomplish this by simply adding as many terms as needed to the truncated farfield expansions of Wilcox or Karp. The convergence of these expansions guarantees that the order of approximation of the new ABC can be increased arbitrarily without having to enlarge the radius of the artificial boundary. We include numerical results in two and three dimensions which demonstrate the improved accuracy and simplicity of this new formulation when compared to other absorbing boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photoacoustic imaging taking into account thermodynamic attenuation

In this paper we consider a mathematical model for photoacoustic imaging which takes into account attenuation due to thermodynamic dissipation. The propagation of acoustic (compressional) waves is governed by a scalar wave equation coupled to the heat equation for the excess temperature. We seek to recover the initial acoustic profile from knowledge of acoustic measurements at the boundary. We recognize that this inverse problem is a special case of boundary observability for a thermoelastic system. This leads to the use of control/observability tools to prove the unique and stable recovery of the initial acoustic profile in the weak thermoelastic coupling regime. This approach is constructive, yielding a solvable equation for the unknown acoustic profile. Moreover, the solution to this reconstruction equation can be approximated numerically using the conjugate gradient method. If certain geometrical conditions for the wave speed are satisfied, this approach is well--suited for variable media and for measurements on a subset of the boundary. We also present a numerical implementation of the proposed reconstruction algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A control approach to recover the wave speed (conformal factor) from one measurement

In this paper we consider the problem of recovering the conformal factor in a conformal class of Riemannian metrics from the boundary measurement of one wave field. More precisely, using boundary control operators, we derive an explicit equation satisfied by the contrast between two conformal factors (or wave speeds). This equation is Fredholm and generically invertible provided that the domain of interest is properly illuminated at an initial time. We also show locally Lipschitz stability estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiwave imaging in an enclosure with variable wave speed

In this paper we consider the mathematical model of thermo- and photo-acoustic tomography for the recovery of the initial condition of a wave field from knowledge of its boundary values. Unlike the free-space setting, we consider the wave problem in a region enclosed by a surface where an impedance boundary condition is imposed. This condition models the presence of physical boundaries such as interfaces or acoustic mirrors which reflect some of the wave energy back into the enclosed domain. By recognizing that the inverse problem is equivalent to a statement of boundary observability, we use control operators to prove the unique and stable recovery of the initial wave profile from knowledge of boundary measurements. Since our proof is constructive, we explicitly derive a solvable equation for the unknown initial condition. This equation can be solved numerically using the conjugate gradient method. We also propose an alternative approach based on the stabilization of waves. This leads to an exponentially and uniformly convergent Neumann series reconstruction when the impedance coefficient is not identically zero. In both cases, if well-known geometrical conditions are satisfied, our approaches are naturally suited for variable wave speed and for measurements on a subset of the boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical method of characteristics for one-dimensional blood flow

Mathematical modeling at the level of the full cardiovascular system requires the numerical approximation of solutions to a one-dimensional nonlinear hyperbolic system describing flow in a single vessel. This model is often simulated by computationally intensive methods like finite elements and discontinuous Galerkin, while some recent applications require more efficient approaches (e.g. for real-time clinical decision support, phenomena occurring over multiple cardiac cycles, iterative solutions to optimization/inverse problems, and uncertainty quantification). Further, the high speed of pressure waves in blood vessels greatly restricts the time step needed for stability in explicit schemes. We address both cost and stability by presenting an efficient and unconditionally stable method for approximating solutions to diagonal nonlinear hyperbolic systems. Theoretical analysis of the algorithm is given along with a comparison of our method to a discontinuous Galerkin implementation. Lastly, we demonstrate the utility of the proposed method by implementing it on small and large arterial networks of vessels whose elastic and geometrical parameters are physiologically relevant.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Recovery of the absorption coefficient in radiative transport from a single measurement

In this paper, we investigate the recovery of the absorption coefficient from boundary data assuming that the region of interest is illuminated at an initial time. We consider a sufficiently strong and isotropic, but otherwise unknown initial state of radiation. This work is part of an effort to reconstruct optical properties using unknown illumination embedded in the unknown medium. We break the problem into two steps. First, in a linear framework, we seek the simultaneous recovery of a forcing term of the form $σ(t,x,θ) f(x)$ (with $σ$ known) and an isotropic initial condition $u_{0}(x)$ using the single measurement induced by these data. Based on exact boundary controllability, we derive a system of equations for the unknown terms $f$ and $u_{0}$. The system is shown to be Fredholm if $σ$ satisfies a certain positivity condition. We show that for generic term $σ$ and weakly absorbing media, this linear inverse problem is uniquely solvable with a stability estimate. In the second step, we use the stability results from the linear problem to address the nonlinearity in the recovery of a weak absorbing coefficient. We obtain a locally Lipschitz stability estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Source estimation with incoherent waves in random waveguides

We study an inverse source problem for the acoustic wave equation in a random waveguide. The goal is to estimate the source of waves from measurements of the acoustic pressure at a remote array of sensors. The waveguide effect is due to boundaries that trap the waves and guide them in a preferred (range) direction, the waveguide axis, along which the medium is unbounded. The random waveguide is a model of perturbed ideal waveguides which have flat boundaries and are filled with known media that do not change with range. The perturbation consists of fluctuations of the boundary and of the wave speed due to numerous small inhomogeneities in the medium. The fluctuations are uncertain in applications, which is why we model them with random processes, and they cause significant cumulative scattering at long ranges from the source. The scattering effect manifests mathematically as an exponential decay of the expectation of the acoustic pressure, the coherent part of the wave. The incoherent wave is modeled by the random fluctuations of the acoustic pressure, which dominate the expectation at long ranges from the source. We use the existing theory of wave propagation in random waveguides to analyze the inverse problem of estimating the source from incoherent wave recordings at remote arrays. We show how to obtain from the incoherent measurements high fidelity estimates of the time resolved energy carried by the waveguide modes, and study the invertibility of the system of transport equations that model energy propagation in order to estimate the source.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On-surface radiation condition for multiple scattering of waves

The formulation of the on-surface radiation condition (OSRC) is extended to handle wave scattering problems in the presence of multiple obstacles. The new multiple-OSRC simultaneously accounts for the outgoing behavior of the wave fields, as well as, the multiple wave reflections between the obstacles. Like boundary integral equations (BIE), this method leads to a reduction in dimensionality (from volume to surface) of the discretization region. However, as opposed to BIE, the proposed technique leads to boundary integral equations with smooth kernels. Hence, these Fredholm integral equations can be handled accurately and robustly with standard numerical approaches without the need to remove singularities. Moreover, under weak scattering conditions, this approach renders a convergent iterative method which bypasses the need to solve single scattering problems at each iteration. Inherited from the original OSRC, the proposed multiple-OSRC is generally a crude approximate method. If accuracy is not satisfactory, this approach may serve as a good initial guess or as an inexpensive pre-conditioner for Krylov iterative solutions of BIE.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the multi-frequency inverse source problem in heterogeneous media

The inverse source problem where an unknown source is to be identified from the knowledge of its radiated wave is studied. The focus is placed on the effect that multi-frequency data has on establishing uniqueness. In particular, it is shown that data obtained from finitely many frequencies is not sufficient. On the other hand, if the frequency varies within a set with an accumulation point, then the source is determined uniquely, even in the presence of highly heterogeneous media. In addition, an algorithm for the reconstruction of the source using multi-frequency data is proposed. The algorithm, based on a subspace projection method, approximates the minimum-norm solution given the available multi-frequency measurements. A few numerical examples are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The DtN nonreflecting boundary condition for multiple scattering problems in the half-plane

The multiple-Dirichlet-to-Neumann (multiple-DtN) non-reflecting boundary condition is adapted to acoustic scattering from obstacles embedded in the half-plane. The multiple-DtN map is coupled with the method of images as an alternative model for multiple acoustic scattering in the presence of acoustically soft and hard plane boundaries. As opposed to the current practice of enclosing all obstacles with a large semicircular artificial boundary that contains portion of the plane boundary, the proposed technique uses small artificial circular boundaries that only enclose the immediate vicinity of each obstacle in the half-plane. The adapted multiple-DtN condition is simultaneously imposed in each of the artificial circular boundaries. As a result the computational effort is significantly reduced. A computationally advantageous boundary value problem is numerically solved with a finite difference method supported on boundary-fitted grids. Approximate solutions to problems involving two scatterers of arbitrary geometry are presented. The proposed numerical method is validated by comparing the approximate and exact far-field patterns for the scattering from a single and from two circular obstacles in the half-plane.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time reversal for radiative transport with applications to inverse and control problems

In this paper we develop a time reversal method for the radiative transport equation to solve two problems: an inverse problem for the recovery of an initial condition from boundary measurements, and the exact boundary controllability of the transport field with finite steering time. Absorbing and scattering effects, modeled by coefficients with low regularity, are incorporated in the formulation of these problems. This time reversal approach leads to a convergent iterative procedure to reconstruct the initial condition provided that the scattering coefficient is sufficiently small in the $L^{\infty}$ norm. Then, using duality arguments, we show that the solvability of the inverse problem leads to exact controllability of the transport field with minimum-norm control obtained constructively. The solution approach to both of these problems may have medical applications in areas such as optical imaging and optimization of radiation delivery for cancer therapy.

0 citations0 reviewsNo code linkedOpen access