Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “validation method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

FAITH – Fast Assembly Inhibitor Test for HIV

Due to the high number of drug-resistant HIV-1 mutants generated by highly active antiretroviral therapy (HAART), there is continuing demand for new types of inhibitors. Both the assembly of the Gag polyprotein into immature and mature HIV-1 particles are attractive candidates for the blocking of the retroviral life cycle. Currently, no therapeutically-used assembly inhibitor is available. One possible explanation is the lack of a reliable and simple assembly inhibitor screening method. To identify compounds potentially inhibiting the formation of both types of HIV-1 particles, we developed a new fluorescent high-throughput screening assay. This assay is based on the quantification of the assembly efficiency in vitro in a 96-well plate format. The key components of the assay are HIV-1 Gag-derived proteins and a dual-labelled oligonucleotide, which emits fluorescence only when the assembly of retroviral particles is inhibited. The method was validated using three (CAI, BM2, PF74) reported assembly inhibitors. - Highlights: • Allows screening of assembly inhibitors of both mature and immature HIV-1 particles. • Based on Gag-derived proteins with CA in mature or immature conformation. • Simple and sensitive method suitable for high-throughput screening of inhibitors. • Unlike in other HIV assembly methods, works under physiological conditions. • No washing steps are necessary.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Fault Clearing Operation in Low-Frequency High-Voltage AC Systems

Abstract—Low-frequency high-voltage ac (LF-HVac) system has been proposed as an alternative to HVdc and HVac systems to reduce transmission losses and improve grid flexibility. LF-HVac system employs back-to-back (BTB) converters to interconnect and control the voltage and power interchange between two power grid operating at different frequencies. While there are multiple studies on LF-HVac system in the literature, the fault clearing method in the LF-HVac system has not been analyzed. Therefore, this paper focuses on fault clearing operation in LFHVac system and proposes a novel method to reduce fault current by changing operating frequency of power source and converter. The performance of the proposed method is validated using detailed models of BTB converters developed in PSCAD/EMTDC. The ability to change operating frequency at the onset of a fault condition allows breakers to rapidly clear the fault and recloses. The results obtained by using the proposed method substantially diminish the fault current.

Kim, Woosung↗
DOE OSTI2022

Characterizing and Mitigating Intraday Variability: Reconstructing Source Structure in Accreting Black Holes with mm-VLBI

The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87* and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse (u, v) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87*. We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87*, finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Technical Note: Cortical thickness and density estimation from clinical CT using a prior thickness-density relationship

Purpose: Cortical thickness and density are critical components in determining the strength of bony structures. Computed tomography (CT) is one possible modality for analyzing the cortex in 3D. In this paper, a model-based approach for measuring the cortical bone thickness and density from clinical CT images is proposed. Methods: Density variations across the cortex were modeled as a function of the cortical thickness and density, location of the cortex, density of surrounding tissues, and imaging blur. High resolution micro-CT data of cadaver proximal femurs were analyzed to determine a relationship between cortical thickness and density. This thickness-density relationship was used as prior information to be incorporated in the model to obtain accurate measurements of cortical thickness and density from clinical CT volumes. The method was validated using micro-CT scans of 23 cadaver proximal femurs. Simulated clinical CT images with different voxel sizes were generated from the micro-CT data. Cortical thickness and density were estimated from the simulated images using the proposed method and compared with measurements obtained using the micro-CT images to evaluate the effect of voxel size on the accuracy of the method. Then, 19 of the 23 specimens were imaged using a clinical CT scanner. Cortical thickness and density were estimated from the clinical CT images using the proposed method and compared with the micro-CT measurements. Finally, a case-control study including 20 patients with osteoporosis and 20 age-matched controls with normal bone density was performed to evaluate the proposed method in a clinical context. Results: Cortical thickness (density) estimation errors were 0.07 ± 0.19 mm (−18 ± 92 mg/cm{sup 3}) using the simulated clinical CT volumes with the smallest voxel size (0.33 × 0.33 × 0.5 mm{sup 3}), and 0.10 ± 0.24 mm (−10 ± 115 mg/cm{sup 3}) using the volumes with the largest voxel size (1.0 × 1.0 × 3.0 mm{sup 3}). A trend for the cortical thickness and density estimation errors to increase with voxel size was observed and was more pronounced for thin cortices. Using clinical CT data for 19 of the 23 samples, mean errors of 0.18 ± 0.24 mm for the cortical thickness and 15 ± 106 mg/cm{sup 3} for the density were found. The case-control study showed that osteoporotic patients had a thinner cortex and a lower cortical density, with average differences of −0.8 mm and −58.6 mg/cm{sup 3} at the proximal femur in comparison with age-matched controls (p-value < 0.001). Conclusions: This method might be a promising approach for the quantification of cortical bone thickness and density using clinical routine imaging techniques. Future work will concentrate on investigating how this approach can improve the estimation of mechanical strength of bony structures, the prevention of fracture, and the management of osteoporosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-201-03: A Real-Time Method to Simultaneously Measure Linear Energy Transfer and Dose for Proton Therapy Using Organic Scintillators

Purpose: The light generated in organic scintillators depends on both the radiation dose and the linear energy transfer (LET). The LET dependence leads to an under-response of the detector in the Bragg peak of proton beams. This phenomenon, called ionization quenching, must be corrected to obtain accurate dose measurements of proton beams. This work exploits the ionization quenching phenomenon to provide a method of measuring LET and auto correcting quenching. Methods: We exposed simultaneously four different organic scintillators (BCF-12, PMMA, PVT, and LSD; 1mm in diameter) and a plane parallel ionization chamber in passively scattered proton beams to doses between 32 and 43 cGy and fluence averaged LET values from 0.47 to 1.26 keV/µm. The LET values for each irradiation condition were determined using a validated Monte Carlo model of the beam line. We determined the quenching parameter in the Birk’s equation for scintillation in BCF-12 for dose measurements. One set of irradiation conditions was used to correlate the scintillation response ratio to the LET values and plot a scintillation response ratio versus LET calibration curve. Irradiation conditions independent from the calibration ones were used to validate this method. Comparisons to the expected values were made on both the basis of dose and LET. Results: Among all the scintillators investigated, the ratio of PMMA to BCF-12 provided the best correlation to LET values and was used as the LET calibration curve. The expected LET values in the validation set were within 2%±6%, which resulted in dose accuracy of 1.5%±5.8% for the range of LET values investigated in this work. Conclusion: We have demonstrated the feasibility of using the ratio between the light output of two organic scintillators to simultaneously measure LET and dose of therapeutic proton beams. Further studies are needed to verify the response in higher LET values.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

An Evaluation of Actinide Reactivity with CO 2 , O 2 , and O 2 /He Gases using Inductively Coupled Plasma Tandem Mass Spectrometry: Application to Simultaneous Measurement of 241 Am/ 241 Pu Ratios in Unseparated Complex Matrices

Accurate actinide measurements are critical within the field of nuclear science. Traditional methods for actinide quantification require time-consuming sample processing prior to analysis. There is a need for rapid analytical techniques that still maintain a high degree of accuracy. In this work, actinide reactivity was assessed for multiple oxygen-containing reaction gases using quadrupole inductively coupled plasma tandem mass spectrometry (Q-ICP-MS/MS) to evaluate actinide analysis in complex sample matrices without analyte-matrix separation. A novel method was developed to measure 241 Am/ 241 Pu in complex sample matrices using O 2 /He reaction gas with no matrix removal or analyte pre-concentration. This inline method reduces matrix-derived polyatomic interferences that complicate traditional ICP-MS analyses by mass-shifting to 241 Am 16 O + and 241 Pu 16 O 2 + , allowing Am and Pu to be mass separated for simultaneous analysis. While mass shifting is efficient, a small portion of Am + (<1.3%) and Pu + (<1.4%) react to from AmO 2 + and PuO + , respectively. Therefore, a mass balance approach was used, in combination with reactivity determined from 242 Pu and 243 Am standard solutions, to correct for residual 241 PuO + and 241 AmO 2 + . The method was validated by measuring 241 Am/ 241 Pu in Pu isotope standards CRM-136 and CRM-137 (separated in March/April 1970 and February 2022, respectively) in both neat solutions and complex matrices containing diluted soil (NIST SRM 2711a, >1000 µg·g -1 ). Method detection limits of 15.9 and 9.6 fg·g -1 were determined for 241 Am and 241 Pu, respectively, and 241 Am/ 241 Pu ratios were measured with accuracies within <3.5%. In conclusion, this work presents the first direct analysis of 241 Am/ 241 Pu in unseparated complex matrices, advancing capabilities for rapid actinide measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗
DOE OSTI2018

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: a tomographic measurement of cosmic structure growth and expansion rate based on optimal redshift weights

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the Extended Zel’dovich mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of 0.8 < z < 2.2. We report a joint measurement of fσ8 and two-dimensional BAO parameters D A and H at four effective redshifts of z eff = 0.98, 1.23, 1.52, and 1.94, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, Main Galaxy Sample (MGS), and 6 degree Field Galaxy Survey (6dFGS) BAO measurements, we find that the existence of dark energy is supported by observations at a 7.4σ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be γ = 0.580 ± 0.082, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On state-of-charge determination for lithium-ion batteries

Accurate estimation of state-of-charge (SOC) of a battery through its life remains challenging in battery research. Although improved precisions continue to be reported at times, almost all are based on regression methods empirically, while accuracy is often not properly addressed. Here, a comprehensive review is set to address such issues, from fundamental principles that are supposed to define SOC to methodologies to estimate SOC for practical use. It covers topics from calibration, regression (including modeling methods) to validation in terms of precision and accuracy. At the end, we intend to answer the following questions: 1) can SOC estimation be self-adaptive without bias? 2) Why Ah-counting is a necessity in almost all battery-model-assisted regression methods? 3) How to establish a consistent framework of coupling in multi-physics battery models? 4) To assess the accuracy in the SOC estimation, statistical methods should be employed to analyze factors that contribute to the uncertainty. Finally, we hope, through this proper discussion of the principles, accurate SOC estimation can be widely achieved.

25 ENERGY STORAGE↗
DOE OSTI2018

In situ coherent diffractive imaging

Coherent diffractive imaging (CDI) has been widely applied in the physical and biological sciences using synchrotron radiation, X-ray free-electron laser, high harmonic generation, electrons, and optical lasers. One of CDI’s important applications is to probe dynamic phenomena with high spatiotemporal resolution. Here, we report the development of a general in situ CDI method for real-time imaging of dynamic processes in solution. By introducing a time-invariant overlapping region as real-space constraint, we simultaneously reconstructed a time series of complex exit wave of dynamic processes with robust and fast convergence. We validated this method using optical laser experiments and numerical simulations with coherent X-rays. Our numerical simulations further indicated that in situ CDI can potentially reduce radiation dose by more than an order of magnitude relative to conventional CDI. With further development, we envision in situ CDI could be applied to probe a range of dynamic phenomena in the future.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Extraction of the electron excess temperature in terahertz quantum cascade lasers from laser characteristics

Abstract We propose a method to extract the upper laser level’s (ULL’s) excess electronic temperature from the analysis of the maximum light output power ( P max ) and current dynamic range Δ J d = ( J max − J th ) of terahertz quantum cascade lasers (THz QCLs). We validated this method, both through simulation and experiment, by applying it on THz QCLs supporting a clean three-level system. Detailed knowledge of electronic excess temperatures is of utmost importance in order to achieve high temperature performance of THz QCLs. Our method is simple and can be easily implemented, meaning an extraction of the excess electron temperature can be achieved without intensive experimental effort. This knowledge should pave the way toward improvement of the temperature performance of THz QCLs beyond the state-of-the-art.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Using high-throughput amplicon sequencing to determine diet of generalist lady beetles in agricultural landscapes

Determining feeding relationships is central to understanding biological control potential in the field. However, methods to differentiate actual (or realized) feeding relationships from potential feeding relationships is lacking especially for small, generalist predators such as lady beetles. In this study, we used DNA metabarcoding approaches to characterize actual feeding relationships of lady beetles (Coccinellidae) in the field and validated our methods with a lab study. We first asked whether high-throughput amplicon sequencing (HTS) can characterize diets of lady beetles ranging from monotypic diets to diverse diet mixtures in the lab. We then examined whether diet composition and breadth of lady beetles collected from different habitat types in southern WI varied between monocultures of soybean and corn, diverse tallgrass prairie, and urban habitats. Lastly, we asked whether different body or tissue types (partial-body versus whole-body specimens) would change the likelihood of prey detection for both studies. In our controlled lab study, we found that HTS can accurately assess diet composition and diet breadth for lady beetle populations, but at the individual level, HTS has limitation for individuals feeding on more than three species of prey at any given time. In our field study, we documented lower prey richness and diet breadth in corn compared to soybean and grassland, and greater prey DNA in regurgitants than in whole-body specimens. Finally, we found that reduced diet diversity was associated with an increased prevalence of intraguild predation, but that habitat and prey diversity are not necessarily correlated. The prey detection rates (proportion of specimens with prey DNA) in our study were comparable to other studies (55% in the lab study, 33% in field study) and varied with diet composition and habitat type. Furthermore, the relatively low detection rates suggest that many consumer individuals would need to be assayed to fully assess diet diversity, especially in diverse systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Electromagnetic Shower from Cosmic Muons in NOvA Detector

NO$\nu$A (NuMI Off-Axis $\nu_e$ Appearance) experiment is a neutrino $\nu_\mu \rightarrow \nu_e$ oscillation experiment with a baseline of 810 km. The signature for a $\nu_\mu \rightarrow \nu_e$ oscillation in the far detector is the production of an electromagnetic (EM) shower due to the charge current (CC) interaction of $\nu_e$ with the detector material. Cosmic muon induced electromagnetic showers can serve as a source of EM shower sample that can be used as a data-driven method for validating the EM shower modelling and EM shower identification algorithms at NO$\nu$A. The work reported in the thesis is the first ever implementation of cosmic muon induced EM shower for such a purpose in the NO$\nu$A experiment. We extracted EM shower sample from two sources - muon undergoing bremsstrahlung (Brem) and muon decay in flight to electron (DiF). Brem sample, extracted from cosmic muon, turns out to be a very useful data-driven method for testing EM shower modelling and the part icle identification (PID) algorithms for the detection of $\nu_e$-CC interactions. A new muon removal (MR) technique is developed and applied to remove the muon hits inside the bremsstrahlung shower. However, a complete removal of the muon energy in the event is not possible with high efficiency. DiF induced EM shower, on the other hand, is exactly similar to $\nu_e$-CC EM shower because DiF is a pure electron initiated EM shower. Therefore, DiF EM shower, once extracted, mimicks the shower of a $\nu_e$-CC interaction event. We present the extraction methods for both brem and DiF events in MC and cosmic data in far and near detector of NO$\nu$A and its use in data driven validation of $\nu_e$ EM shower modelling and $\nu_e$ signal selection efficiency across the detector. A systematic uncertainty is also assigned on $\nu_e$ appearance signal selection based on this study to take into account the difference in EM shower modelling in data and MC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

A Feature-Based Stochastic Permeability of Shale: Part 1—Validation and Two-Phase Permeability in a Utica Shale Sample

Estimate of permeability plays a crucial role in flow-based studies of fractured tight-rocks. It is well known that most of the flow through tight-rocks (e.g., shales) is controlled by permeable features (e.g., fractures, laminations, etc.), and there is negligible flow through the matrix. However, current approaches in the literature to model permeability of tight-rocks do not account for such features present within the rock ranging from micro-scale to field-scale. Current permeability modeling approach assumes a single continuum without considering the presence of permeable features within the matrix (e.g., micro-fractures) or outside the matrix (e.g., natural fractures). Although the laboratory-measured permeability implicitly captures discrete features present in that sample (e.g., fractures, laminations, micro-fractures), most of the permeability models proposed for shale do not account for these features. Fracture permeability in the literature is typically modeled using an ideal slit assumption; however, this highly overestimates its permeability because fractures in real medium are non-ideal in terms of their porosity and tortuosity, which affect their permeability. Additionally, the transition zone between fracture and matrix also affects the permeability of fracture. Here in this study, part of a two-part series, a new method to predict permeability of fractured shale by discretizing the medium into matrix (inorganic and organic) and fractures is presented. New analytical expressions of permeability are derived to account for non-ideal nature of porous medium and two-phase flow in fractures. Rock feature in each cell of the grid is identified as one of the three elements (organic matter, inorganic matter, or fracture), and permeability of that cell is estimated using a suitable analytical expression. This method allows estimating permeability at any scale of interest and more robustly than by a pure analytical approach. The proposed method is validated against local and global-scale measurements on three fractured samples from laboratory. Finally, the method is used to predict two-phase flow permeability of supercritical CO 2 displacing water within a fracture in a Utica shale sample. The proposed two-phase flow permeability equations can be used as a quick analytical tool to predict relative permeability estimates of two-phase flow in fractured shale samples. In Part 2, the proposed method is used to estimate field-scale permeability through an optimization process that uses field-scale production and other readily available information.

58 GEOSCIENCES↗
DOE OSTI2024

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

TU-H-CAMPUS-IeP3-01: Simultaneous PET Restoration and PET/CT Co-Segmentation Using a Variational Method

Purpose: PET images are usually blurred due to the finite spatial resolution, while CT images suffer from low contrast. Segment a tumor from either a single PET or CT image is thus challenging. To make full use of the complementary information between PET and CT, we propose a novel variational method for simultaneous PET image restoration and PET/CT images co-segmentation. Methods: The proposed model was constructed based on the Γ-convergence approximation of Mumford-Shah (MS) segmentation model for PET/CT co-segmentation. Moreover, a PET de-blur process was integrated into the MS model to improve the segmentation accuracy. An interaction edge constraint term over the two modalities were specially designed to share the complementary information. The energy functional was iteratively optimized using an alternate minimization (AM) algorithm. The performance of the proposed method was validated on ten lung cancer cases and five esophageal cancer cases. The ground truth were manually delineated by an experienced radiation oncologist using the complementary visual features of PET and CT. The segmentation accuracy was evaluated by Dice similarity index (DSI) and volume error (VE). Results: The proposed method achieved an expected restoration result for PET image and satisfactory segmentation results for both PET and CT images. For lung cancer dataset, the average DSI (0.72) increased by 0.17 and 0.40 than single PET and CT segmentation. For esophageal cancer dataset, the average DSI (0.85) increased by 0.07 and 0.43 than single PET and CT segmentation. Conclusion: The proposed method took full advantage of the complementary information from PET and CT images. This work was supported in part by the National Cancer Institute Grants R01CA172638. Shan Tan and Laquan Li were supported in part by the National Natural Science Foundation of China, under Grant Nos. 60971112 and 61375018.

60 APPLIED LIFE SCIENCES↗