Source author record

Loc H. Nguyen

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2022arXiv

A Carleman-based numerical method for quasilinear elliptic equations with over-determined boundary data and applications

We propose a new iterative scheme to compute the numerical solution to an over-determined boundary value problem for a general quasilinear elliptic PDE. The main idea is to repeatedly solve its linearization by using the quasi-reversibility method with a suitable Carleman weight function. The presence of the Carleman weight function allows us to employ a Carleman estimate to prove the convergence of the sequence generated by the iterative scheme above to the desired solution. The convergence of the iteration is fast at an exponential rate without the need of an initial good guess. We apply this method to compute solutions to some general quasilinear elliptic equations and a large class of first-order Hamilton-Jacobi equations. Numerical results are presented.

preprint2022arXiv

Carleman estimates and the contraction principle for an inverse source problem for nonlinear hyperbolic equation

The main aim of this paper is to solve an inverse source problem for a general nonlinear hyperbolic equation. Combining the quasi-reversibility method and a suitable Carleman weight function, we define a map of which fixed point is the solution to the inverse problem. To find this fixed point, we define a recursive sequence with an arbitrary initial term by the same manner as in the classical proof of the contraction principle. Applying a Carleman estimate, we show that the sequence above converges to the desired solution with the exponential rate. Therefore, our new method can be considered as an analog of the contraction principle. We rigorously study the stability of our method with respect to noise. Numerical examples are presented.

preprint2022arXiv

Convexification for a CIP for the RTE]{Convexification Numerical Method for a Coefficient Inverse Problem for the Radiative Transport Equation

An $\left( n+1\right) -$D coefficient inverse problem for the radiative stationary transport equation is considered for the first time. A globally convergent so-called convexification numerical \ method is developed and its convergence analysis is provided. The analysis is based on a Carleman estimate. In particular, convergence analysis implies a certain uniqueness theorem. Extensive numerical studies in the 2-D case are presented.

preprint2022arXiv

The Carleman convexification method for Hamilton-Jacobi equations on the whole space

We propose a new globally convergent numerical method to solve Hamilton-Jacobi equations in $\mathbb{R}^d$, $d \geq 1$. This method is named as the Carleman convexification method. By Carleman convexification, we mean that we use a Carleman weight function to convexify the conventional least squares mismatch functional. We will prove a new version of the convexification theorem guaranteeing that the mismatch functional involving the Carleman weight function is strictly convex and, therefore, has a unique minimizer. Moreover, a consequence of our convexification theorem guarantees that the minimizer of the Carleman weighted mismatch functional is an approximation of the viscosity solution we want to compute. Some numerical results in 1D and 2D will be presented.

preprint2021arXiv

Carleman contraction mapping for a 1D inverse scattering problem with experimental time-dependent data

It is shown that the contraction mapping principle with the involvement of a Carleman Weight Function works for a Coefficient Inverse Problem for a 1D hyperbolic equation. Using a Carleman estimate, the global convergence of the corresponding numerical method is established. Numerical studies for both computationally simulated and experimentally collected data are presented. The experimental part is concerned with the problem of computing dielectric constants of explosive-like targets in the standoff mode using severely underdetermined data.

preprint2021arXiv

Numerical viscosity solutions to Hamilton-Jacobi equations via a Carleman estimate and the convexification method

We propose a globally convergent numerical method, called the convexification, to numerically compute the viscosity solution to first-order Hamilton-Jacobi equations through the vanishing viscosity process where the viscosity parameter is a fixed small number. By convexification, we mean that we employ a suitable Carleman weight function to convexify the cost functional defined directly from the form of the Hamilton-Jacobi equation under consideration. The strict convexity of this functional is rigorously proved using a new Carleman estimate. We also prove that the unique minimizer of the this strictly convex functional can be reached by the gradient descent method. Moreover, we show that the minimizer well approximates the viscosity solution of the Hamilton-Jacobi equation as the noise contained in the boundary data tends to zero. Some interesting numerical illustrations are presented.

preprint2021arXiv

The quasi-reversibility method to numerically solve an inverse source problem for hyperbolic equations

We propose a numerical method to solve an inverse source problem of computing the initial condition of hyperbolic equations from the measurements of Cauchy data. This problem arises in thermo- and photo- acoustic tomography in a bounded cavity, in which the reflection of the wave makes the widely-used approaches, such as the time reversal method, not applicable. In order to solve this inverse source problem, we approximate the solution to the hyperbolic equation by its Fourier series with respect to a special orthogonal basis of $L^2$. Then, we derive a coupled system of elliptic equations for the corresponding Fourier coefficients. We solve it by the quasi-reversibility method. The desired initial condition follows. We rigorously prove the convergence of the quasi-reversibility method as the noise level tends to 0. Some numerical examples are provided. In addition, we numerically prove that the use of the special basic above is significant.

preprint2020arXiv

An inverse problem of a simultaneous reconstruction of the dielectric constant and conductivity from experimental backscattering data

This report extends our recent progress in tackling a challenging 3D inverse scattering problem governed by the Helmholtz equation. Our target application is to reconstruct dielectric constants, electric conductivities and shapes of front surfaces of objects buried very closely under the ground. These objects mimic explosives, like, e.g., antipersonnel land mines and improvised explosive devices. We solve a coefficient inverse problem with the backscattering data generated by a moving source at a fixed frequency. This scenario has been studied so far by our newly developed convexification method that consists in a new derivation of a boundary value problem for a coupled quasilinear elliptic system. However, in our previous work only the unknown dielectric constants of objects and shapes of their front surfaces were calculated. Unlike this, in the current work performance of our numerical convexification algorithm is verified for the case when the dielectric constants, the electric conductivities and those shapes of objects are unknown. By running several tests with experimentally collected backscattering data, we find that we can accurately image both the dielectric constants and shapes of targets of interests including a challenging case of targets with voids. The computed electrical conductivity serves for reliably distinguishing conductive and non-conductive objects. The global convergence of our numerical procedure is shortly revisited.

preprint2020arXiv

Convexification and experimental data for a 3D inverse scattering problem with the moving point source

Inverse scattering problems of the reconstructions of physical properties of a medium from boundary measurements are substantially challenging ones. This work aims to verify the performance on experimental data of a newly developed convexification method for a 3D coefficient inverse problem for the case of objects buried in a sandbox a fixed frequency and the point source moving along an interval of a straight line. Using a special Fourier basis, the method of this work strongly relies on a new derivation of a boundary value problem for a system of coupled quasilinear elliptic equations. This problem, in turn, is solved via the minimization of a Tikhonov-like functional weighted by a Carleman Weight Function. The global convergence of the numerical procedure is established analytically. The numerical verification is performed using experimental data, which are raw backscatter data of the electric field. These data were collected using a microwave scattering facility at The University of North Carolina at Charlotte.

preprint2020arXiv

Convexification for a 1D Hyperbolic Coefficient Inverse Problem with Single Measurement Data

A version of the convexification numerical method for a Coefficient Inverse Problem for a 1D hyperbolic PDE is presented. The data for this problem are generated by a single measurement event. This method converges globally. The most important element of the construction is the presence of the Carleman Weight Function in a weighted Tikhonov-like functional. This functional is strictly convex on a certain bounded set in a Hilbert space, and the diameter of this set is an arbitrary positive number. The global convergence of the gradient projection method is established. Computational results demonstrate a good performance of the numerical method for noisy data.

preprint2016arXiv

A globally convergent method for a 3-D inverse medium problem for the generalized Helmholtz equation

A 3-D inverse medium problem in the frequency domain is considered. Another name for this problem is Coefficient Inverse Problem. The goal is to reconstruct spatially distributed dielectric constants from scattering data. Potential applications are in detection and identification of explosive-like targets. A single incident plane wave and multiple frequencies are used. A new numerical method is proposed. A theorem is proved, which claims that a small neigborhood of the exact solution of that problem is reached by this method without any advanced knowledge of that neighborhood. We call this property of that numerical method "global convergence". Results of numerical experiments for the case of the backscattering data are presented.

preprint2016arXiv

A globally convergent numerical method for a 1-d inverse medium problem with experimental data

In this paper, a reconstruction method for the spatially distributed dielectric constant of a medium from the back scattering wave field in the frequency domain is considered. Our approach is to propose a globally convergent algorithm, which does not require any knowledge of a small neighborhood of the solution of the inverse problem in advance. The Quasi-Reversibility Method (QRM) is used in the algorithm. The convergence of the QRM is proved via a Carleman estimate. The method is tested on both computationally simulated and experimental data.

preprint2016arXiv

A globally convergent numerical method for a 3D coefficient inverse problem with a single measurement of multi-frequency data

The goal of this paper is to reconstruct spatially distributed dielectric constants from complex-valued scattered wave field by solving a 3D coefficient inverse problem for the Helmholtz equation at multi-frequencies. The data are generated by only a single direction of the incident plane wave. To solve this inverse problem, a globally convergent algorithm is analytically developed. We prove that this algorithm provides a good approximation for the exact coefficient without any \textit{a priori} knowledge of any point in a small neighborhood of that coefficient. This is the main advantage of our method, compared with classical approaches using optimization schemes. Numerical results are presented for both computationally simulated data and experimental data. Potential applications of this problem are in detection and identification of explosive-like targets.

preprint2016arXiv

Single measurement experimental data for an inverse medium problem inverted by a multi-frequency globally convergent numerical method

The recently developed globally convergent numerical method for an inverse medium problem for the Helmholtz equation is tested on experimental data. The data were originally collected in the time domain, whereas the method works in the frequency domain with the multi-frequency data. Due to a huge discrepancy between the collected and computationally simulated data, the straightforward Fourier transform of the experimental data does not work. Hence, it is necessary to develop a heuristic data preprocessing procedure. This procedure is described. The preprocessed data are used as the input for the inversion algorithm. Numerical results demonstrate good accuracy in the reconstruction of both refracive indices and locations of targets. Furthermore, the reconstruction errors for refractive indices of dielectric targets are significantly less than errors of a posteriori direct measurements.

preprint2015arXiv

Nanostructures imaging via numerical solution of a 3-d inverse scattering problem without the phase information

Inverse scattering problems without the phase information arise in imaging of nanostructures whose sizes are hundreds of nanometers as well as in imaging of biological cells. The governing equation is the 3-d generalized Helmholtz equation with the unknown coefficient, which represents the spatially distributed dielectric constant. It is assumed that only the modulus of the complex valued wave field is measured on a frequency range. The phase is not measured. Two rigorous numerical methods are tested and their performances are compared for realistic ranges of parameters.