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 MATLAB Code for Three Dimensional Linear Elastostatics using Constant Boundary Elements

Present work presents a code written in the very simple programming language MATLAB, for three dimensional linear elastostatics, using constant boundary elements. The code, in full or in part, is not a translation or a copy of any of the existing codes. Present paper explains how the code is written, and lists all the formulae used. Code is verified by using the code to solve a simple problem which has the well known approximate analytical solution. Of course, present work does not make any contribution to research on boundary elements, in terms of theory. But the work is justified by the fact that, to the best of author’s knowledge, as of now, one cannot find an open access MATLAB code for three dimensional linear elastostatics using constant boundary elements. Author hopes this paper to be of help to beginners who wish to understand how a simple but complete boundary element code works, so that they can build upon and modify the present open access code to solve complex engineering problems quickly and easily. The code is available online for open access (as supplementary file for the present paper), and may be downloaded from the website for the present journal.

Extracting Three Dimensional Surface Model of Human Kidney from the Visible Human Data Set using Free Software

Three dimensional digital model of a representative human kidney is needed for a surgical simulator that is capable of simulating a laparoscopic surgery involving kidney. Buying a three dimensional computer model of a representative human kidney, or reconstructing a human kidney from an image sequence using commercial software, both involve (sometimes significant amount of) money. In this paper, author has shown that one can obtain a three dimensional surface model of human kidney by making use of images from the Visible Human Data Set and a few free software packages (ImageJ, ITK-SNAP, and MeshLab in particular). Images from the Visible Human Data Set, and the software packages used here, both do not cost anything. Hence, the practice of extracting the geometry of a representative human kidney for free, as illustrated in the present work, could be a free alternative to the use of expensive commercial software or to the purchase of a digital model.

Porous Three Dimensional Arrays of Plasmonic Nanoparticles

Plasmonic interactions in a well-defined array of metallic nanoparticles can lead to interesting optical effects, such as local electric field enhancement and shifts in the extinction spectra, which are of interest in diverse technological applications, including those pertaining to biochemical sensing and photonic circuitry. Here, we report on a single-step wafer scale fabrication of a three-dimensional array of metallic nanoparticles whose sizes and separations can be easily controlled to be anywhere between fifty to a few hundred nanometers, allowing the optical response of the system to be tailored with great control in the visible region of the spectrum. The substrates, apart from having a large surface area, are inherently porous and therefore suitable for optical sensing applications, such as surface enhanced Raman scattering, containing a high density of spots with enhanced local electric fields arising from plasmonic couplings.

Recapitulating tumour microenvironment in chitosan-gelatin three-dimensional scaffolds: an improved in vitro tumour model

Owing to the reduced co-relationship between conventional flat Petri dish culture (two-dimensional) and the tumour microenvironment, there has been a shift towards three-dimensional culture systems that show an improved analogy to the same. In this work, an extracellular matrix (ECM)-mimicking three-dimensional scaffold based on chitosan and gelatin was fabricated and explored for its potential as a tumour model for lung cancer. It was demonstrated that the chitosan-gelatin (CG) scaffolds supported the formation of tumoroids that were similar to tumours grown in vivo for factors involved in tumour-cell-ECM interaction, invasion and metastasis, and response to anti-cancer drugs. On the other hand, the two-dimensional Petri dish surfaces did not demonstrate gene-expression profiles similar to tumours grown in vivo. Further, the three-dimensional CG scaffolds supported the formation of tumoroids, using other types of cancer cells such as breast, cervix and bone, indicating a possible wider potential for in vitro tumoroid generation. Overall, the results demonstrated that CG scaffolds can be an improved in vitro tool to study cancer progression and drug screening for solid tumours.

Demonstrating the Usefulness of CAELinux for Computer Aided Engineering using an Example of the Three Dimensional Reconstruction of a Pig Liver

CAELinux is a Linux distribution which is bundled with free software packages related to Computer Aided Engineering (CAE). The free software packages include software that can build a three dimensional solid model, programs that can mesh a geometry, software for carrying out Finite Element Analysis (FEA), programs that can carry out image processing etc. Present work has two goals: 1) To give a brief description of CAELinux 2) To demonstrate that CAELinux could be useful for Computer Aided Engineering, using an example of the three dimensional reconstruction of a pig liver from a stack of CT-scan images. One can note that instead of using CAELinux, using commercial software for reconstructing the liver would cost a lot of money. One can also note that CAELinux is a free and open source operating system and all software packages that are included in the operating system are also free. Hence one can conclude that CAELinux could be a very useful tool in application areas like surgical simulation which require three dimensional reconstructions of biological organs. Also, one can see that CAELinux could be a very useful tool for Computer Aided Engineering, in general.

Three-dimensional acoustic analysis of an elliptical chamber muffler with a side inlet and a side outlet

The acoustical behaviour of an elliptical chamber muffler having a side inlet and side outlet port is analyzed in this paper, wherein a uniform velocity piston source is assumed to model the 3-D acoustic field in the elliptical chamber cavity. Towards this end, we consider the modal expansion of the acoustic pressure field in the elliptical cavity in terms of the angular and radial Mathieu func-tions, subjected to the rigid wall condition. Then, the Green's function due to the point source lo-cated on the side (curved) surface of the elliptical chamber is obtained. On integrating this function over the elliptical piston area on the curved surface of the elliptical chamber and subsequent divi-sion by the area of the elliptic piston, one obtains the acoustic pressure field due to the piston driven source which is equivalent to considering plane wave propagation in the side ports. Thus, one can obtain the acoustic pressure response functions, i.e., the impedance matrix (Z) parameters due to the sources (ports) located on the side surface, from which one may also obtain a progressive wave rep-resentation in terms of the scattering matrix (S). Finally, the acoustic performance of the muffler is evaluated in terms of the Transmission loss (TL) which is computed in terms of the scattering pa-rameters. The effect of the axial length of the muffler and the angular location of the ports on the TL characteristics is studied in detail. The acoustically long chambers show dominant axial plane wave propagation while the TL spectrum of short chambers indicates the dominance of the trans-versal modes. The 3-D analytical results are compared with the 3-D FEM simulations carried on a commercial software and are shown to be in an excellent agreement, thereby validating the analyti-cal procedure suggested in this work.

Wafer scale fabrication of porous three-dimensional plasmonic metamaterials for the visible region: chiral and beyond

We report on a wafer scale fabrication method of a three-dimensional plasmonic metamaterial with strong chiroptical response in the visible region of the electromagnetic spectrum. The system was comprised of metallic nanoparticles arranged in a helical fashion, with high degree of flexibility over the choice of the underlying material, as well as their geometrical parameters. This resulted in exquisite control over the chiroptical properties, most importantly the spectral signature of the circular dichroism. In spite of the large variability in the arrangement, as well as the size and shape of the constituent nanoparticles, the average chiro-optical response of the material remained uniform across the wafer, thus confirming the suitability of this system as a large area chiral metamaterial. By simply heating the substrate for a few minutes, the geometrical properties of the nanoparticles could be altered, thus providing an additional handle towards tailoring the spectral response of this novel material.

Three dimensional PN based impact angle control for higher speed nonmaneuvering targets

Terminal impact angle control is crucial for enhancement of warhead effectiveness. In the literature, this problem has been addressed in the context of targets with lower speeds than the interceptor. However, in the current defence scenario, targets of much higher speed than the interceptor is a reality. This paper presents a generic proportional navigation (PN) based guidance law, that uses the standard PN and novel Retro-PN guidance laws based on the initial engagement geometry and terminal engagement requirements, for three dimensional engagement scenario against higher speed nonmaneuvering targets to control terminal impact angle. Results are obtained on the set of achievable impact angles and conditions on the navigation constant to achieve them. Simulation results are given to support the theoretical findings.

Codes for solving three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces

Three codes, that can solve three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces, are provided here. The file 'bemconst.m' contains a MATLAB code for solving three dimensional linear elastostatic problems using constant boundary elements while ignoring body forces. The file 'bemconst.f90' is a Fortran translation of the MATLAB code contained in the file 'bemconst.m'. The file 'bemconstp.f90' is a parallelized version of the Fortran code contained in the file 'bemconst.f90'. The file 'inbem96.txt' is the input file for the Fortran codes contained in the files 'bemconst.f90' and 'bemconstp.f90'. Author hereby declares that the present codes are the original works of the author. Further, author hereby declares that any of the present codes, in full or in part, is not a translation or a copy of any of the existing codes written by someone else. Author's institution (Indian Institute of Science) has informed the author in writing that the institution is not interested in claiming any copyright on the present codes. Author is hereby distributing the present codes under the MIT License; full text of the license is included in each of the files that contain the codes.

Hypersonic interference heating on flat plate with short three-dimensional protuberances

Experiments were conducted to measure the heat flux in the vicinity of a three-dimensional protuberance placed on a flat plate facing a hypersonic flow at zero angle of attack. The effects of flow enthalpy and height of the protuberance on the interference heating in its vicinity were studied. Evidence of disturbed flow with highly three-dimensional characteristics and heightened vorticity was observed near the protrusion. A parametric study by changing the deflection angle of the protuberance was also made. Correlations exist in the open literature for enthalpy values lower than 2  MJ/kg. This effort has yielded a new correlation that is valid for enthalpies up to 6  MJ/kg. The Z-type schlieren technique was used to visualize the flow features qualitatively for one of the flow conditions studied.

Optimal trajectory planning for path convergence in three-dimensional space

This article addresses the problem of determining the shortest path that connects a given initial configuration (position, heading angle, and flight path angle) to a given rectilinear or a circular path in three-dimensional space for a constant speed and turn-rate constrained aerial vehicle. The final path is assumed to be located relatively far from the starting point. Due to its simplicity and low computational requirements the algorithm can be implemented on a fixed-wing type unmanned air vehicle in real time in missions where the final path may change dynamically. As wind has a very significant effect on the flight of small aerial vehicles, the method of optimal path planning is extended to meet the same objective in the presence of wind comparable to the speed of the aerial vehicles. But, if the path to be followed is closer to the initial point, an off-line method based on multiple shooting, in combination with a direct transcription technique, is used to obtain the optimal solution. Optimal paths are generated for a variety of cases to show the efficiency of the algorithm. Simulations are presented to demonstrate tracking results using a 6-degrees-of-freedom model of an unmanned air vehicle.

Three-Dimensional Study of Natural Convection in a Horizontal Channel With Discrete Heaters on One of Its Vertical Walls

Three-dimensional natural convection in a horizontal channel with an array of discrete flush-mounted heaters on one of its vertical walls is numerically studied. Effects of thermal conductivities of substrate and heaters and convection on outer sides of the channel walls on heat transfer are examined. The substrate affects heat transfer in a wider range of thermal conductivities than do the heaters. At lower heater thermal conductivities a higher heat portion is transferred by direct convection from the heaters to the adjacent coolant. However, higher substrate conductivity is associated with higher heat portion transferred through the substrate. The innermost heater column is found to become the hottest heater column due to the lower coolant accessibility. The heat transfer in the channel is strongly influenced by convection on the outer sides of the channel walls. Correlations are presented for dimensionless temperature maximum and average Nusselt number.

Three Dimensional Fluorescence Imaging Using Multiple Light-Sheet Microscopy

We developed a multiple light-sheet microscopy (MLSM) system capable of 3D fluorescence imaging. Employing spatial filter in the excitation arm of a SPIM system, we successfully generated multiple light-sheets. This improves upon the existing SPIM system and is capable of 3D volume imaging by simultaneously illuminating multiple planes in the sample. Theta detection geometry is employed for data acquisition from multiple specimen layers. This detection scheme inherits many advantages including, background reduction, cross-talk free fluorescence detection and high-resolution at long working distance. Using this technique, we generated 5 equi-intense light-sheets of thickness approximately 7: 5 mm with an inter-sheet separation of 15 mm. Moreover, the light-sheets generated by MLSM is found to be 2 times thinner than the state-of-art SPIM system. Imaging of fluorescently coated yeast cells of size 4 +/- 1 mm (encaged in Agarose gel-matrix) is achieved. Proposed imaging technique may accelerate the field of fluorescence microscopy, cell biology and biophotonics.