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

Sebastian Acosta

Sebastian Acosta contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NAphysics.comp-phNumerical Analysismath-phmath.MPmath.OC

Trust snapshot

Quick read

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

12 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
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