Practical fully three-dimensional reconstruction algorithms for diffuse optical tomography

We have developed an efficient fully three-dimensional (3D) reconstruction algorithm for diffuse optical tomography (DOT). The 3D DOT, a severely ill-posed problem, is tackled through a pseudodynamic (PD) approach wherein an ordinary differential equation representing the evolution of the solution on pseudotime is integrated that bypasses an explicit inversion of the associated, ill-conditioned system matrix. One of the most computationally expensive parts of the iterative DOT algorithm, the reevaluation of the Jacobian in each of the iterations, is avoided by using the adjoint-Broyden update formula to provide low rank updates to the Jacobian. In addition, wherever feasible, we have also made the algorithm efficient by integrating along the quadratic path provided by the perturbation equation containing the Hessian. These algorithms are then proven by reconstruction, using simulated and experimental data and verifying the PD results with those from the popular Gauss-Newton scheme. The major findings of this work are as follows: (i) the PD reconstructions are comparatively artifact free, providing superior absorption coefficient maps in terms of quantitative accuracy and contrast recovery; (ii) the scaling of computation time with the dimension of the measurement set is much less steep with the Jacobian update formula in place than without it; and (iii) an increase in the data dimension, even though it renders the reconstruction problem less ill conditioned and thus provides relatively artifact-free reconstructions, does not necessarily provide better contrast property recovery. For the latter, one should also take care to uniformly distribute the measurement points, avoiding regions close to the source so that the relative strength of the derivatives for measurements away from the source does not become insignificant. (c) 2012 Optical Society of America

A PRACTICAL PROCEDURE TO ESTIMATE THE BEARING CAPACITY OF FOOTINGS ON SAND-APPLICATION TO 87 CASE STUDIES

Following the recent work of the authors in development and numerical verification of a new kinematic approach of the limit analysis for surface footings on non-associative materials, a practical procedure is proposed to utilize the theory. It is known that both the peak friction angle and dilation angle depend on the sand density as well as the stress level, which was not the concern of the former work. In the current work, a practical procedure is established to provide a better estimate of the bearing capacity of surface footings on sand which is often non-associative. This practical procedure is based on the results obtained theoretically and requires the density index and the critical state friction angle of the sand. The proposed practical procedure is a simple iterative computational procedure which relates the density index of the sand, stress level, dilation angle, peak friction angle and eventually the bearing capacity. The procedure is described and verified among available footing load test data.