Georgia Institute of Technology College of Sciences

Cone Beam Computer tomography (CBCT) has been broadly applied incancer radiation therapy, mainly for positioning patients to align withtreatment radiation beams. As opposed to tomography reconstruction problemsfor diagnostic purposes, CBCT reconstruction in radiotherapy requires a highcomputational efficiency, since it is performed while patient is lying on acouch, waiting for the treatment. Moreover, the excessive radiation dosefrom frequent scans has become a clinical concern. It is therefore desirableto develop new techniques to reconstruct CBCT images from low dose scans. Inthis talk, I will present our recent work on an iterative low dose CBCTreconstruction technique via total variation regularization and tight frameregularization. It is found that 40~60 x-ray projections are sufficient toreconstruct a volumetric image with satisfactory quality in about 2min. Wehave also studied 4 dimensional CBCT (4DCBCT) reconstruction problem viatemporal non-local means (TNLM) and high quality 4DCBCT images can beobtained. Our algorithms have been fully implemented on a graphicsprocessing unit. Detailed implementation techniques will also be addressed.

In this talk, we present a parallel finite element implementation ontetrahedral  grids of the nonlinear three-dimensional nonlinear Stokes model for thedynamics and evolution of ice-sheets. Discretization is based on a high-orderaccurate  scheme using the Taylor-Hood element pair. Both no-slip and sliding boundary conditions at the ice-bedrock boundary are studied. In addition, effective solvers using preconditioning techniques for the saddle-point system resulting fromthe  discretization are discussed and implemented. We demonstrate throughestablished ice-sheet benchmark experiments that our finite element nonlinear Stokesmodel  performs at least as well as other published and established Stokes modelsin the  field, and the parallel solver is shown to be efficient, robust, and scalable.

The inverse problem of optical tomography consists of recovering thespatially-varying absorption of a highly-scattering medium from boundarymeasurements. In this talk we will discuss direct reconstruction methods forthis problem that are based on inversion of the Born series. In previouswork we have utilized such series expansions as tools to develop fast imagereconstruction algorithms. Here we characterize their convergence, stabilityand approximation error. Analogous results for the Calderon problem ofreconstructing the conductivity in electrical impedance tomography will alsobe presented.

Vortex dynamics and solid-fluid interaction are two of the most important and most studied topics in fluid dynamics for their relevance to a wide range of applications from geophysical flows to locomotion in moving fluids. In this talk, we investigate two problems in these two areas: Part I studies the viscous evolution of point vortex equilibria; Part II studies the effects of body elasticity on the passive stability of submerged bodies.In Part I, we describe the viscous evolution of point vortex configurations that, in the absence of viscosity, are in a state of fixed or relative equilibrium. In particular, we examine four cases, three of them correspond to relative equilibria in the inviscid point vortex model and one corresponds to a fixed equilibrium. Our goal is to elucidate some of the main transient dynamical features of the flow. Using a multi-Gaussian ``core growing" type of model, we show that all four configurations immediately begin to rotate unsteadily, while the shapes of vortex configurations remain unchanged. We then examine in detail the qualitative and quantitative evolution of the structures as they evolve, and for each case show the sequence of topological bifurcations that occur both in a fixed reference frame, and in an appropriately chosen rotating reference frame. Comparisons between the cases help to reveal different features of the viscous evolution for short and intermediate time ! scales of vortex structures. The dynamical evolution of passive particles in the viscously evolving flow associated with the initial fixed equilibrium is shown and interpreted in relation to the evolving streamline patterns. In Part II, we examine the effects of body geometry and elasticity on the passive stability of motion in a perfect fluid. Our main motivation is to understand the role of body elasticity on the stability of fish swimming. The fish is modeled as an articulated body made of multiple links (assumed to be identical ellipses in 2D or identical ellipsoids in 3D) interconnected by hinge joints. It can undergo shape changes by varying the relative angles between the links. Body elasticity is accounted for via the torsional springs at the joints. The unsteadiness of the flow is modeled using the added mass effect. Equations of motion for the body-fluid system are derived using Newtonian and Lagrangian approaches for both hydrodynamically decoupled and coupled models in 2D and 3D. We specifically examine the stability associated with a relative equilibrium of the equations, traditionally referred to as the ``coast motion" (proved to be unstable for a rigid elongated body model), and f! ound that body elasticity does stabilize the system. Stable regions are identified based on linear stability analysis in the parameter space spanned by aspect ratio (body geometry) and spring constants (muscle stiffness), and the findings based on the linear analysis are verified by direct numerical simulations of the nonlinear system.