Tang, Jet Hoe. Solving forward and inverse Helmholtz equations via controllability methods. 2019, Doctoral Thesis, University of Basel, Faculty of Science.

PDF
54Mb 
Official URL: http://edoc.unibas.ch/diss/DissB_13528
Downloads: Statistics Overview
Abstract
Waves are useful for probing an unknown medium by illuminating it with a source.
To infer the characteristics of the medium from (boundary) measurements,
for instance, one typically formulates inverse scattering problems
in frequency domain as a PDEconstrained optimization problem.
Finding the medium, where the simulated wave field
matches the measured (real) wave field, the inverse problem
requires the repeated solutions of forward (Helmholtz) problems.
Typically, standard numerical methods, e.g. direct solvers or iterative methods,
are used to solve the forward problem.
However, largescaled (or highfrequent) scattering problems
are known being competitive in computation and storage for standard methods.
Moreover, since the optimization problem is severely illposed
and has a large number of
local minima, the inverse problem requires additional regularization
akin to minimizing the total variation.
Finding a suitable regularization for the inverse problem is critical
to tackle the illposedness and to reduce the computational cost and storage requirement.
In my thesis, we first apply standard methods to forward problems.
Then, we consider the controllability method (CM)
for solving the forward problem: it
instead reformulates the problem in the time domain
and seeks the timeharmonic solution of the corresponding wave equation.
By iteratively reducing the mismatch between the solution at
initial time and after one period with the conjugate gradient (CG) method,
the CMCG method greatly speeds up the convergence to the timeharmonic
asymptotic limit. Moreover, each conjugate gradient iteration
solely relies on standard numerical algorithms,
which are inherently parallel and robust against higher frequencies.
Based on the original CM, introduced in 1994 by Bristeau et al.,
for soundsoft scattering problems, we extend the CMCG method to
general boundaryvalue problems governed by the Helmholtz equation.
Numerical results not only show the usefulness, robustness, and efficiency
of the CMCG method for solving the forward problem,
but also demonstrate remarkably accurate solutions.
Second, we formulate the PDEconstrained optimization
problem governed by the inverse scattering problem
to reconstruct the unknown medium.
Instead of a gridbased discrete representation combined with
standard Tikhonovtype regularization, the unknown medium is
projected to a small finitedimensional subspace,
which is iteratively adapted using dynamic thresholding.
The adaptive (spectral) space is governed by solving
several Poissontype eigenvalue problems.
To tackle the illposedness that the Newtontype optimization
method converges to a false local minimum,
we combine the adaptive spectral inversion (ASI) method with the frequency stepping strategy.
Numerical examples illustrate the usefulness of the ASI approach,
which not only efficiently and remarkably reduces the dimension of the
solution space, but also yields an accurate and robust method.
To infer the characteristics of the medium from (boundary) measurements,
for instance, one typically formulates inverse scattering problems
in frequency domain as a PDEconstrained optimization problem.
Finding the medium, where the simulated wave field
matches the measured (real) wave field, the inverse problem
requires the repeated solutions of forward (Helmholtz) problems.
Typically, standard numerical methods, e.g. direct solvers or iterative methods,
are used to solve the forward problem.
However, largescaled (or highfrequent) scattering problems
are known being competitive in computation and storage for standard methods.
Moreover, since the optimization problem is severely illposed
and has a large number of
local minima, the inverse problem requires additional regularization
akin to minimizing the total variation.
Finding a suitable regularization for the inverse problem is critical
to tackle the illposedness and to reduce the computational cost and storage requirement.
In my thesis, we first apply standard methods to forward problems.
Then, we consider the controllability method (CM)
for solving the forward problem: it
instead reformulates the problem in the time domain
and seeks the timeharmonic solution of the corresponding wave equation.
By iteratively reducing the mismatch between the solution at
initial time and after one period with the conjugate gradient (CG) method,
the CMCG method greatly speeds up the convergence to the timeharmonic
asymptotic limit. Moreover, each conjugate gradient iteration
solely relies on standard numerical algorithms,
which are inherently parallel and robust against higher frequencies.
Based on the original CM, introduced in 1994 by Bristeau et al.,
for soundsoft scattering problems, we extend the CMCG method to
general boundaryvalue problems governed by the Helmholtz equation.
Numerical results not only show the usefulness, robustness, and efficiency
of the CMCG method for solving the forward problem,
but also demonstrate remarkably accurate solutions.
Second, we formulate the PDEconstrained optimization
problem governed by the inverse scattering problem
to reconstruct the unknown medium.
Instead of a gridbased discrete representation combined with
standard Tikhonovtype regularization, the unknown medium is
projected to a small finitedimensional subspace,
which is iteratively adapted using dynamic thresholding.
The adaptive (spectral) space is governed by solving
several Poissontype eigenvalue problems.
To tackle the illposedness that the Newtontype optimization
method converges to a false local minimum,
we combine the adaptive spectral inversion (ASI) method with the frequency stepping strategy.
Numerical examples illustrate the usefulness of the ASI approach,
which not only efficiently and remarkably reduces the dimension of the
solution space, but also yields an accurate and robust method.
Advisors:  Grote, Marcus J. and Gander, Martin J. 

Faculties and Departments:  05 Faculty of Science > Departement Mathematik und Informatik > Mathematik > Numerik (Grote) 
UniBasel Contributors:  Grote, Marcus J. 
Item Type:  Thesis 
Thesis Subtype:  Doctoral Thesis 
Thesis no:  13528 
Thesis status:  Complete 
Number of Pages:  1 OnlineRessource (139 Seiten) 
Language:  English 
Identification Number: 

edoc DOI:  
Last Modified:  25 Feb 2020 05:30 
Deposited On:  24 Feb 2020 15:39 
Repository Staff Only: item control page