Engineering Papers⌕ Search

Engineering topics

Zhong, Z

Publications and source records attributed to Zhong, Z.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

TU-AB-202-05: GPU-Based 4D Deformable Image Registration Using Adaptive Tetrahedral Mesh Modeling

Purpose: Deformable image registration (DIR) has been employed today as an automated and effective segmentation method to transfer tumor or organ contours from the planning image to daily images, instead of manual segmentation. However, the computational time and accuracy of current DIR approaches are still insufficient for online adaptive radiation therapy (ART), which requires real-time and high-quality image segmentation, especially in a large datasets of 4D-CT images. The objective of this work is to propose a new DIR algorithm, with fast computational speed and high accuracy, by using adaptive feature-based tetrahedral meshing and GPU-based parallelization. Methods: The first step is to generate the adaptive tetrahedral mesh based on the image features of a reference phase of 4D-CT, so that the deformation can be well captured and accurately diffused from the mesh vertices to voxels of the image volume. Subsequently, the deformation vector fields (DVF) and other phases of 4D-CT can be obtained by matching each phase of the target 4D-CT images with the corresponding deformed reference phase. The proposed 4D DIR method is implemented on GPU, resulting in significantly increasing the computational efficiency due to its parallel computing ability. Results: A 4D NCAT digital phantom was used to test the efficiency and accuracy of our method. Both the image and DVF results show that the fine structures and shapes of lung are well preserved, and the tumor position is well captured, i.e., 3D distance error is 1.14 mm. Compared to the previous voxel-based CPU implementation of DIR, such as demons, the proposed method is about 160x faster for registering a 10-phase 4D-CT with a phase dimension of 256×256×150. Conclusion: The proposed 4D DIR method uses feature-based mesh and GPU-based parallelism, which demonstrates the capability to compute both high-quality image and motion results, with significant improvement on the computational speed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-D-207-04: GPU-Based 4D Cone-Beam CT Reconstruction Using Adaptive Meshing Method

Purpose: Due to the limited number of projections at each phase, the image quality of a four-dimensional cone-beam CT (4D-CBCT) is often degraded, which decreases the accuracy of subsequent motion modeling. One of the promising methods is the simultaneous motion estimation and image reconstruction (SMEIR) approach. The objective of this work is to enhance the computational speed of the SMEIR algorithm using adaptive feature-based tetrahedral meshing and GPU-based parallelization. Methods: The first step is to generate the tetrahedral mesh based on the features of a reference phase 4D-CBCT, so that the deformation can be well captured and accurately diffused from the mesh vertices to voxels of the image volume. After the mesh generation, the updated motion model and other phases of 4D-CBCT can be obtained by matching the 4D-CBCT projection images at each phase with the corresponding forward projections of the deformed reference phase of 4D-CBCT. The entire process of this 4D-CBCT reconstruction method is implemented on GPU, resulting in significantly increasing the computational efficiency due to its tremendous parallel computing ability. Results: A 4D XCAT digital phantom was used to test the proposed mesh-based image reconstruction algorithm. The image Result shows both bone structures and inside of the lung are well-preserved and the tumor position can be well captured. Compared to the previous voxel-based CPU implementation of SMEIR, the proposed method is about 157 times faster for reconstructing a 10 -phase 4D-CBCT with dimension 256×256×150. Conclusion: The GPU-based parallel 4D CBCT reconstruction method uses the feature-based mesh for estimating motion model and demonstrates equivalent image Result with previous voxel-based SMEIR approach, with significantly improved computational speed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

The relationship between coefficient of restitution and state of charge of zinc alkaline primary LR6 batteries [Bouncing alkaline batteries: A basic solution]

The coefficient of restitution of alkaline batteries had been shown to increase as a function of depth of discharge. In this work, using non-destructive mechanical testing, the change in coefficient of restitution is compared to in situ energy-dispersive x-ray diffraction data to determine the cause of the macroscopic change in coefficient of restitution. The increase in coefficient of restitution correlates to the formation of a percolation pathway of ZnO within the anode of the cell, and that the coefficient of restitution saturates at a value of 0.63 ± .05 at 50% state if charge when the anode has densified into porous ZnO solid. Of note is the sensitivity of coefficient of restitution to the amount of ZnO formation that rivals the sensitivity on in situ energy-dispersive x-ray diffraction spectroscopy.

25 ENERGY STORAGE↗
DOE OSTI2014

SU-F-BRF-09: A Non-Rigid Point Matching Method for Accurate Bladder Dose Summation in Cervical Cancer HDR Brachytherapy

Purpose: To propose and validate a deformable point matching scheme for surface deformation to facilitate accurate bladder dose summation for fractionated HDR cervical cancer treatment. Method: A deformable point matching scheme based on the thin plate spline robust point matching (TPSRPM) algorithm is proposed for bladder surface registration. The surface of bladders segmented from fractional CT images is extracted and discretized with triangular surface mesh. Deformation between the two bladder surfaces are obtained by matching the two meshes' vertices via the TPS-RPM algorithm, and the deformation vector fields (DVFs) characteristic of this deformation is estimated by B-spline approximation. Numerically, the algorithm is quantitatively compared with the Demons algorithm using five clinical cervical cancer cases by several metrics: vertex-to-vertex distance (VVD), Hausdorff distance (HD), percent error (PE), and conformity index (CI). Experimentally, the algorithm is validated on a balloon phantom with 12 surface fiducial markers. The balloon is inflated with different amount of water, and the displacement of fiducial markers is benchmarked as ground truth to study TPS-RPM calculated DVFs' accuracy. Results: In numerical evaluation, the mean VVD is 3.7(±2.0) mm after Demons, and 1.3(±0.9) mm after TPS-RPM. The mean HD is 14.4 mm after Demons, and 5.3mm after TPS-RPM. The mean PE is 101.7% after Demons and decreases to 18.7% after TPS-RPM. The mean CI is 0.63 after Demons, and increases to 0.90 after TPS-RPM. In the phantom study, the mean Euclidean distance of the fiducials is 7.4±3.0mm and 4.2±1.8mm after Demons and TPS-RPM, respectively. Conclusions: The bladder wall deformation is more accurate using the feature-based TPS-RPM algorithm than the intensity-based Demons algorithm, indicating that TPS-RPM has the potential for accurate bladder dose deformation and dose summation for multi-fractional cervical HDR brachytherapy. This work is supported in part by the National Natural ScienceFoundation of China (no 30970866 and no 81301940)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-D-9A-01: A Novel Mesh-Based Deformable Surface-Contour Registration

Purpose: Initial guess is vital for 3D-2D deformable image registration (DIR) while dealing with large deformations for adaptive radiation therapy. A fast procedure has been developed to deform body surface to match 2D body contour on projections. This surface-contour DIR will provide an initial deformation for further complete 3D DIR or image reconstruction. Methods: Both planning CT images and come-beam CT (CBCT) projections are preprocessed to create 0–1 binary mask. Then the body surface and CBCT projection body contours are extracted by Canny edge detector. A finite element modeling system was developed to automatically generate adaptive meshes based on the image surface. After that, the projections of the CT surface voxels are computed and compared with corresponding 2D projection contours from CBCT scans. As a result, the displacement vector field (DVF) on mesh vertices around the surface was optimized iteratively until the shortest Euclidean distance between the pixels on the projections of the deformed CT surface and the corresponding CBCT projection contour is minimized. With the help of the tetrahedral meshes, we can smoothly diffuse the deformation from the surface into the interior of the volume. Finally, the deformed CT images are obtained by the optimal DVF applied on the original planning CT images. Results: The accuracy of the surface-contour registration is evaluated by 3D normalized cross correlation increased from 0.9176 to 0.9957 (sphere-ellipsoid phantom) and from 0.7627 to 0.7919 (H and N cancer patient data). Under the GPU-based implementation, our surface-contour-guided method on H and N cancer patient data takes 8 seconds/iteration, about 7.5 times faster than direct 3D method (60 seconds/iteration), and it needs fewer optimization iterations (30 iterations vs 50 iterations). Conclusion: The proposed surface-contour DIR method can substantially improve both the accuracy and the speed of reconstructing volumetric images, which is helpful for applying in adaptive radiotherapy. This research is supported by CPRIT individual investigator award RP110329.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Pressure effects on phase equilibria and solid solubility in MgO-Y2O3 nanocomposites

We study the temperature and pressure dependence of phase evolution in the 0.5MgO-0.5Y{sub 2}O{sub 3} nanocomposite system using a diamond anvil apparatus in conjunction with in situ synchrotron energy dispersive x-ray diffraction at 7 GPa hydrostatic pressure. At (298 K, 7.0 GPa), structural transformations in the Y{sub 2}O{sub 3} phase are observed, giving rise to the co-existence of its cubic, hexagonal, and monoclinic polymorphs together with cubic MgO. An increase in temperature to 1273 K causes the crystallinity of the Y{sub 2}O{sub 3} hexagonal and monoclinic phases to increase. Isothermal and isobaric hold at (1273 K, 7.0 GPa) for 60 min results in yttrium dissolution in cubic MgO, causing {approx}1.0% expansive volumetric lattice strain despite the large differences in the ionic radii of the cations. Cooling the nanocomposite to (298 K, 0 GPa) after a 60 min soak yields four phase co-existence among cubic MgO and cubic, hexagonal, and monoclinic Y{sub 2}O{sub 3}. The residual MgO unit cell volume expansion is 0.69% at 298 K, indicating solid solution formation at room temperature despite large differences in the ionic radii of Mg{sup 2+} and Y{sup 3+}. The macroscopic shrinkage due to densification is 3% by volume. Thermodynamic considerations suggest that the relative molar partial volume of Y{sup 3+} in MgO is a negative quantity, indicating that the partial molar volume of Y{sup 3+} in the solid solution is smaller than its molar volume in the pure state. Aging of the nanocomposites for 240 h does not change the observed 4 phase co-existence. We propose a crystallographic model in which the observed volumetric expansion of the MgO unit cell is primarily attributed to two hydrostatic expansive strain components accompanying solid solution formation: (i) Coulomb repulsion among O{sup 2-} ions in the immediate vicinity of Mg{sup 2+} vacancies, and (ii) misfit strain due to differences in ionic radii upon Y{sup 3+} substitution on Mg{sup 2+} sites.

36 MATERIALS SCIENCE↗
DOE OSTI2011

High-pressure phase transformations in MgO-Y2O3 nanocomposites

Temperature and pressure dependence of phase evolution in 0.5 MgO-0.5 Y{sub 2}O{sub 3} is studied using the diamond anvil method. At 300 K and 5.5 GPa, transformations in Y{sub 2}O{sub 3} are observed, resulting in co-existence of cubic, hexagonal, and monoclinic phases. Heating to 1293 K results in increased crystallinity and increase in hexagonal and monoclinic phase content. Isothermal and isobaric hold at 1273 K and 5.5 GPa for 120 min results in yttrium dissolution in cubic MgO, causing 0.83% volumetric strain. Cooling to 300 K and 0 GPa yields a four phase co-existence among cubic MgO and cubic, hexagonal, and monoclinic Y{sub 2}O{sub 3}. The residual MgO unit cell volume expansion is 0.67% at 300 K, indicating solid solution formation. Aging of the nanocomposites for 240 h does not change the observed 4 phase co-existence.

42 ENGINEERING↗
DOE OSTI2011

Surface Curvatures and Diffraction Profiles of Sagittaly Bent Laue Crystals

The performance of a bent Laue crystal monochromator crucially depends on the sagittal and meridional bending curvatures of the crystal. To optimize the design of monochromator crystals, the surface curvatures and diffraction profiles of a set of sagittally bent Laue crystals with different aspect ratios have been studied experimentally by optical metrology and X-ray measurements. The results were confirmed with finite-element analysis using large-deformation theory. The nonlinear relationship between the curvatures necessitates an experimentally determined parameter in the theoretical modeling of the diffraction profiles. By taking into account the local stress and the aspect ratio of the sagittally bent Laue crystal, the modified analytical approach successfully predicts the rocking-curve width and the integrated reflecting power. The effect of extreme sagittal bending on the rocking curve is also discussed. To retain high reflectivity, the bending curvature should not exceed its critical value for the specified crystal geometry. Furthermore, the uniformity of the bending curvatures across the crystal surface has been examined, which suggests that the minimum crystal dimension should be approximately twice the size of the beam footprint.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗