Robust topology optimization of three-dimensional photonic-crystal band-gap structures

We perform full 3D topology optimization (in which "every voxel" of the unit cell is a degree of freedom) of photonic-crystal structures in order to find optimal omnidirectional band gaps for various symmetry groups, including fcc (including diamond), bcc, and simple-cubic lattices. Even without imposing the constraints of any fabrication process, the resulting optimal gaps are only slightly larger than previous hand designs, suggesting that current photonic crystals are nearly optimal in this respect. However, optimization can discover new structures, e.g. a new fcc structure with the same symmetry but slightly larger gap than the well known inverse opal, which may offer new degrees of freedom to future fabrication technologies. Furthermore, our band-gap optimization is an illustration of a computational approach to 3D dispersion engineering which is applicable to many other problems in optics, based on a novel semidefinite-program formulation for nonconvex eigenvalue optimization combined with other techniques such as a simple approach to impose symmetry constraints. We also demonstrate a technique for \emph{robust} topology optimization, in which some uncertainty is included in each voxel and we optimize the worst-case gap, and we show that the resulting band gaps have increased robustness to systematic fabrication errors.Comment: 17 pages, 9 figures, submitted to Optics Expres

Beam Orientation Optimization for Intensity Modulated Radiation Therapy using Adaptive l1 Minimization

Beam orientation optimization (BOO) is a key component in the process of IMRT treatment planning. It determines to what degree one can achieve a good treatment plan quality in the subsequent plan optimization process. In this paper, we have developed a BOO algorithm via adaptive l_1 minimization. Specifically, we introduce a sparsity energy function term into our model which contains weighting factors for each beam angle adaptively adjusted during the optimization process. Such an energy term favors small number of beam angles. By optimizing a total energy function containing a dosimetric term and the sparsity term, we are able to identify the unimportant beam angles and gradually remove them without largely sacrificing the dosimetric objective. In one typical prostate case, the convergence property of our algorithm, as well as the how the beam angles are selected during the optimization process, is demonstrated. Fluence map optimization (FMO) is then performed based on the optimized beam angles. The resulted plan quality is presented and found to be better than that obtained from unoptimized (equiangular) beam orientations. We have further systematically validated our algorithm in the contexts of 5-9 coplanar beams for 5 prostate cases and 1 head and neck case. For each case, the final FMO objective function value is used to compare the optimized beam orientations and the equiangular ones. It is found that, our BOO algorithm can lead to beam configurations which attain lower FMO objective function values than corresponding equiangular cases, indicating the effectiveness of our BOO algorithm.Comment: 19 pages, 2 tables, and 5 figure

GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy

Online adaptive radiation therapy (ART) has great promise to significantly reduce normal tissue toxicity and/or improve tumor control through real-time treatment adaptations based on the current patient anatomy. However, the major technical obstacle for clinical realization of online ART, namely the inability to achieve real-time efficiency in treatment re-planning, has yet to be solved. To overcome this challenge, this paper presents our work on the implementation of an intensity modulated radiation therapy (IMRT) direct aperture optimization (DAO) algorithm on graphics processing unit (GPU) based on our previous work on CPU. We formulate the DAO problem as a large-scale convex programming problem, and use an exact method called column generation approach to deal with its extremely large dimensionality on GPU. Five 9-field prostate and five 5-field head-and-neck IMRT clinical cases with 5\times5 mm2 beamlet size and 2.5\times2.5\times2.5 mm3 voxel size were used to evaluate our algorithm on GPU. It takes only 0.7~2.5 seconds for our implementation to generate optimal treatment plans using 50 MLC apertures on an NVIDIA Tesla C1060 GPU card. Our work has therefore solved a major problem in developing ultra-fast (re-)planning technologies for online ART

Implementation of a motor control system for electric bus based on DSP

© 2017 IEEE. Motor control system may be the most important part of electric vehicles. To implement the control strategies, a lot of practical problems need to be taken into account. In this paper, an induction motor control system for electric bus is developed based on digital signal processor (DSP). The control strategy is based on field-oriented control and space vector pulse width modulation. Over-modulation, field weakening control, PI controller and fault diagnosis are also applied in this DSP algorithm. As a practical product running on a real electric bus with an 100 kW induction motor, communication with vehicle control unit (VCU) by controller area network (CAN bus), control system safety and PC software designed for experiment at lab are also discussed. The transient and steady-state performances of this motor control system are analyzed by experiments. Its performance is satisfactory when applied to the real electric bus

High-Temporal-Resolution High-Spatial-Resolution Spaceborne SAR Based on Continuously Varying PRF

Synthetic Aperture Radar (SAR) is a well-established and powerful imaging technique for acquiring high-spatial-resolution images of the Earth's surface. With the development of beam steering techniques, sliding spotlight and staring spotlight modes have been employed to support high-spatial-resolution applications. In addition to this strengthened high-spatial-resolution and wide-swath capability, high-temporal-resolution (short repeat-observation interval) represents a key capability for numerous applications. However, conventional SAR systems are limited in that the same patch can only be illuminated for several seconds within a single pass. This paper considers a novel high-squint-angle system intended to acquire high-spatial-resolution spaceborne SAR images with repeat-observation intervals varying from tens of seconds to several minutes within a single pass. However, an exponentially increased range cell migration would arise and lead to a conflict between the receive window and 'blind ranges'. An efficient data acquisition technique for high-temporal-resolution, high-spatial-resolution and high-squint-angle spaceborne SAR, in which the pulse repetition frequency (PRF) is continuously varied according to the changing slant range, is presented in this paper. This technique allows echo data to remain in the receive window instead of conflicting with the transmitted pulse or nadir echo. Considering the precision of hardware, a compromise and practical strategy is also proposed. Furthermore, a detailed performance analysis of range ambiguities is provided with respect to parameters of TerraSAR-X. For strong point-like targets, the range ambiguity of this technique would be better than that of uniform PRF technique. For this innovative technique, a resampling strategy and modified imaging algorithm have been developed to handle the non-uniformly sampled echo data. Simulations are performed to validate the efficiency of the proposed technique and the associated imaging algorithm

GPU-based Iterative Cone Beam CT Reconstruction Using Tight Frame Regularization

X-ray imaging dose from serial cone-beam CT (CBCT) scans raises a clinical concern in most image guided radiation therapy procedures. It is the goal of this paper to develop a fast GPU-based algorithm to reconstruct high quality CBCT images from undersampled and noisy projection data so as to lower the imaging dose. For this purpose, we have developed an iterative tight frame (TF) based CBCT reconstruction algorithm. A condition that a real CBCT image has a sparse representation under a TF basis is imposed in the iteration process as regularization to the solution. To speed up the computation, a multi-grid method is employed. Our GPU implementation has achieved high computational efficiency and a CBCT image of resolution 512\times512\times70 can be reconstructed in ~5 min. We have tested our algorithm on a digital NCAT phantom and a physical Catphan phantom. It is found that our TF-based algorithm is able to reconstrct CBCT in the context of undersampling and low mAs levels. We have also quantitatively analyzed the reconstructed CBCT image quality in terms of modulation-transfer-function and contrast-to-noise ratio under various scanning conditions. The results confirm the high CBCT image quality obtained from our TF algorithm. Moreover, our algorithm has also been validated in a real clinical context using a head-and-neck patient case. Comparisons of the developed TF algorithm and the current state-of-the-art TV algorithm have also been made in various cases studied in terms of reconstructed image quality and computation efficiency.Comment: 24 pages, 8 figures, accepted by Phys. Med. Bio

Mechanisms of enhanced heterogeneous nucleation during solidification in binary Al-Mg alloys

This is the post-print version of the Article. The official published version can be accessed from the link below - Copyright @ 2012 ElsevierThe mechanisms involved in the grain refinement of Al–Mg alloys through varying the Mg content and applying intensive melt shearing were investigated. It was found that the oxide formed in Al–Mg alloys under normal melting conditions is MgAl2O4, which displays an equiaxed and faceted morphology with {1 1 1} planes exposed as its natural surfaces. Depending on the Mg content, MgAl2O4 particles exist either as oxide films in dilute Al–Mg alloys (Mg 1 wt.%). Such MgAl2O4 particles can act as potent sites for nucleation of α-Al grains, which is evidenced by the well-defined cube-on-cube orientation relationship between MgAl2O4 and α-Al. Enhanced heterogeneous nucleation in Al–Mg alloys can be attributed to the high potency of MgAl2O4 particles with a lattice misfit of 1.4% and the increased number density of MgAl2O4 particles due to either natural dispersion by the increased Mg content or forced dispersion through intensive melt shearing. It was also found that intensive melt shearing leads to significant grain refinement of dilute Al–Mg alloys by effective dispersion of the MgAl2O4 particles entrapped in oxide films, but it has marginal effect on the grain refinement of concentrated Al–Mg alloys, where MgAl2O4 particles have been naturally dispersed into individual particles by the increased Mg content.This study is funded from the EPSRC Grant EP/H026177/1