Articles | Volume 3
https://doi.org/10.5194/ars-3-219-2005
https://doi.org/10.5194/ars-3-219-2005
12 May 2005
12 May 2005

Time-Domain Techniques for Computation and Reconstruction of One-Dimensional Profiles

M. Rahman and R. Marklein

Abstract. This paper presents a time-domain technique to compute the electromagnetic fields and to reconstruct the permittivity profile within a one-dimensional medium of finite length. The medium is characterized by a permittivity as well as conductivity profile which vary only with depth. The discussed scattering problem is thus one-dimensional. The modeling tool is divided into two different schemes which are named as the forward solver and the inverse solver. The task of the forward solver is to compute the internal fields of the specimen which is performed by Green’s function approach. When a known electromagnetic wave is incident normally on the media, the resulting electromagnetic field within the media can be calculated by constructing a Green’s operator. This operator maps the incident field on either side of the medium to the field at an arbitrary observation point. It is nothing but a matrix of integral operators with kernels satisfying known partial differential equations. The reflection and transmission behavior of the medium is also determined from the boundary values of the Green's operator. The inverse solver is responsible for solving an inverse scattering problem by reconstructing the permittivity profile of the medium. Though it is possible to use several algorithms to solve this problem, the invariant embedding method, also known as the layer-stripping method, has been implemented here due to the advantage that it requires a finite time trace of reflection data. Here only one round trip of reflection data is used, where one round trip is defined by the time required by the pulse to propagate through the medium and back again. The inversion process begins by retrieving the reflection kernel from the reflected wave data by simply using a deconvolution technique. The rest of the task can easily be performed by applying a numerical approach to determine different profile parameters. Both the solvers have been found to have the ability to deal with different types of slabs and incident electromagnetic pulses. Slabs having continuous and discontinuous relative permittivity have already been tested successfully. The tested electromagnetic pulses are a Dirac, Gaussian and sinusoidal pulse. Due to sampling, the resolution of the system also plays a significant role in obtaining better outputs from this scheme.