Lanczos iterated time-reversal.

A new iterative time-reversal algorithm capable of identifying and focusing on multiple scatterers in a relatively small number of iterations is developed. It is recognized that the traditional iterated time-reversal method is based on utilizing power iterations to determine the dominant eigenpairs of the time-reversal operator. The convergence properties of these iterations are known to be suboptimal. Motivated by this, a new method based on Lanczos iterations is developed. In several illustrative examples it is demonstrated that for the same number of transmitted and received signals, the Lanczos iterations based approach is substantially more accurate.

Reconstruction of periodic structures from optical scattering measurements using adjoint equations.

A new numerical approach to efficiently reconstruct the profile of a grating from measurements of reflection coefficients is demonstrated. The problem is posed in the mathematical framework of an inverse scattering problem and solved using gradient-based algorithms. The gradient is computed efficiently using adjoint equations, which amounts to an extra scattering computation per iteration. For symmetric profiles it is shown that only knowledge of the scattered field is sufficient to compute the gradient. As a result, complex profiles can be reconstructed rapidly, and the method can be potentially used in metrology applications in semiconductors. The technique is demonstrated for the case of TE polarization.

A one-dimensional finite element method for simulation-based medical planning for cardiovascular disease.

We have previously described a new approach to planning treatments for cardiovascular disease, Simulation-Based Medical Planning, whereby a physician utilizes computational tools to construct and evaluate a combined anatomic/physiologic model to predict the outcome of alternative treatment plans for an individual patient. Current systems for Simulation-Based Medical Planning utilize finite element methods to solve the time-dependent, three-dimensional equations governing blood flow and provide detailed data on blood flow distribution, pressure gradients and locations of flow recirculation, low wall shear stress and high particle residence. However, these methods are computationally expensive and often require hours of time on parallel computers. This level of computation is necessary for obtaining detailed information about blood flow, but likely is unnecessary for obtaining information about mean flow rates and pressure losses. We describe, herein, a space-time finite element method for solving the one-dimensional equations of blood flow. This method is applied to compute flow rate and pressure in a single segment model, a bifurcation, an idealized model of the abdominal aorta, in three alternate treatment plans for a case of aorto-iliac occlusive disease and in a vascular bypass graft. All of these solutions were obtained in less than 5 min of computation time on a personal computer.