Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulator Study”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Investigating the limits of PET/CT imaging at very low true count rates and high random fractions in ion-beam therapy monitoring

Purpose: External beam radiotherapy with protons and heavier ions enables a tighter conformation of the applied dose to arbitrarily shaped tumor volumes with respect to photons, but is more sensitive to uncertainties in the radiotherapeutic treatment chain. Consequently, an independent verification of the applied treatment is highly desirable. For this purpose, the irradiation-induced β{sup +}-emitter distribution within the patient is detected shortly after irradiation by a commercial full-ring positron emission tomography/x-ray computed tomography (PET/CT) scanner installed next to the treatment rooms at the Heidelberg Ion-Beam Therapy Center (HIT). A major challenge to this approach is posed by the small number of detected coincidences. This contribution aims at characterizing the performance of the used PET/CT device and identifying the best-performing reconstruction algorithm under the particular statistical conditions of PET-based treatment monitoring. Moreover, this study addresses the impact of radiation background from the intrinsically radioactive lutetium-oxyorthosilicate (LSO)-based detectors at low counts. Methods: The authors have acquired 30 subsequent PET scans of a cylindrical phantom emulating a patientlike activity pattern and spanning the entire patient counting regime in terms of true coincidences and random fractions (RFs). Accuracy and precision of activity quantification, image noise, and geometrical fidelity of the scanner have been investigated for various reconstruction algorithms and settings in order to identify a practical, well-suited reconstruction scheme for PET-based treatment verification. Truncated listmode data have been utilized for separating the effects of small true count numbers and high RFs on the reconstructed images. A corresponding simulation study enabled extending the results to an even wider range of counting statistics and to additionally investigate the impact of scatter coincidences. Eventually, the recommended reconstruction scheme has been applied to exemplary postirradiation patient data-sets. Results: Among the investigated reconstruction options, the overall best results in terms of image noise, activity quantification, and accurate geometrical recovery were achieved using the ordered subset expectation maximization reconstruction algorithm with time-of-flight (TOF) and point-spread function (PSF) information. For this algorithm, reasonably accurate (better than 5%) and precise (uncertainty of the mean activity below 10%) imaging can be provided down to 80 000 true coincidences at 96% RF. Image noise and geometrical fidelity are generally improved for fewer iterations. The main limitation for PET-based treatment monitoring has been identified in the small number of true coincidences, rather than the high intrinsic random background. Application of the optimized reconstruction scheme to patient data-sets results in a 25% − 50% reduced image noise at a comparable activity quantification accuracy and an improved geometrical performance with respect to the formerly used reconstruction scheme at HIT, adopted from nuclear medicine applications. Conclusions: Under the poor statistical conditions in PET-based treatment monitoring, improved results can be achieved by considering PSF and TOF information during image reconstruction and by applying less iterations than in conventional nuclear medicine imaging. Geometrical fidelity and image noise are mainly limited by the low number of true coincidences, not the high LSO-related random background. The retrieved results might also impact other emerging PET applications at low counting statistics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-202-07: A Methodology for Generating Numerical Phantoms for Radiation Therapy Using Geometric Attribute Distribution Models

Purpose: A persistent challenge for the quality assessment of radiation therapy treatments (e.g. contouring accuracy) is the absence of the known, ground truth for patient data. Moreover, assessment results are often patient-dependent. Computer simulation studies utilizing numerical phantoms can be performed for quality assessment with a known ground truth. However, previously reported numerical phantoms do not include the statistical properties of inter-patient variations, as their models are based on only one patient. In addition, these models do not incorporate tumor data. In this study, a methodology was developed for generating numerical phantoms which encapsulate the statistical variations of patients within radiation therapy, including tumors. Methods: Based on previous work in contouring assessment, geometric attribute distribution (GAD) models were employed to model both the deterministic and stochastic properties of individual organs via principle component analysis. Using pre-existing radiation therapy contour data, the GAD models are trained to model the shape and centroid distributions of each organ. Then, organs with different shapes and positions can be generated by assigning statistically sound weights to the GAD model parameters. Organ contour data from 20 retrospective prostate patient cases were manually extracted and utilized to train the GAD models. As a demonstration, computer-simulated CT images of generated numerical phantoms were calculated and assessed subjectively and objectively for realism. Results: A cohort of numerical phantoms of the male human pelvis was generated. CT images were deemed realistic both subjectively and objectively in terms of image noise power spectrum. Conclusion: A methodology has been developed to generate realistic numerical anthropomorphic phantoms using pre-existing radiation therapy data. The GAD models guarantee that generated organs span the statistical distribution of observed radiation therapy patients, according to the training dataset. The methodology enables radiation therapy treatment assessment with multi-modality imaging and a known ground truth, and without patient-dependent bias.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Effects of unburnt reaction progress on stretch flame dynamics under elevated temperatures

Flames in practical combustors and engines are inevitably subject to the effects of upstream chemical reaction progress and stretch, induced by elevated thermodynamic conditions and flow non-uniformities, respectively. Recent shock tube experiments and simulation studies on flame propagation have shown that flame speed under engine-relevant conditions can be enhanced with non-negligible upstream chemical reaction progress, especially when low-temperature heat release is involved in the unburnt mixture. On the other hand, depending on the mixture equivalence ratio and diluents, nonequidiffusion (including the non-unity Lewis number effect and the preferential diffusion effect) can couple with the flame stretch to fundamentally affect the flame propagation, which can manifest as either facilitation or suppression. Depending on the transport property of a reacting mixture, there hence can be either inhibition or promotion from unburnt reaction progress and stretch effects on flame propagation. In the current work, through one-dimensional numerical simulations of transient planar and spherical flames of n-heptane/air under elevated thermodynamic conditions, the combined effects of upstream chemical reaction progress and stretch on flame propagation are investigated. Results show that for both lean and rich n-heptane/air mixtures, flame speed can be substantially promoted with reaction progress, while the rich mixtures can exhibit opposite stretch dependence after first-stage ignition. Different definitions of effective Lewis number are adopted to explain the change in Markstein length for spherical flame in lean and rich mixtures reformed by the low temperature chemistry. Further, this work fills an important gap in laminar premixed flame research relevant to practical combustion systems and can provide useful insight into local turbulent flame behaviors and phi-sensitivity in engine combustion.

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

WE-FG-207B-03: Multi-Energy CT Reconstruction with Spatial Spectral Nonlocal Means Regularization

Purpose: Multi-energy computed tomography (MECT) is an emerging application in medical imaging due to its ability of material differentiation and potential for molecular imaging. In MECT, image correlations at different spatial and channels exist. It is desirable to incorporate these correlations in reconstruction to improve image quality. For this purpose, this study proposes a MECT reconstruction technique that employes spatial spectral non-local means (ssNLM) regularization. Methods: We consider a kVp-switching scanning method in which source energy is rapidly switched during data acquisition. For each energy channel, this yields projection data acquired at a number of angles, whereas projection angles among channels are different. We formulate the reconstruction task as an optimziation problem. A least square term enfores data fidelity. A ssNLM term is used as regularization to encourage similarities among image patches at different spatial locations and channels. When comparing image patches at different channels, intensity difference were corrected by a transformation estimated via histogram equalization during the reconstruction process. Results: We tested our method in a simulation study with a NCAT phantom and an experimental study with a Gammex phantom. For comparison purpose, we also performed reconstructions using conjugate-gradient least square (CGLS) method and conventional NLM method that only considers spatial correlation in an image. ssNLM is able to better suppress streak artifacts. The streaks are along different projection directions in images at different channels. ssNLM discourages this dissimilarity and hence removes them. True image structures are preserved in this process. Measurements in regions of interests yield 1.1 to 3.2 and 1.5 to 1.8 times higher contrast to noise ratio than the NLM approach. Improvements over CGLS is even more profound due to lack of regularization in the CGLS method and hence amplified noise. Conclusion: The proposed ssNLM method for kVp-switching MECT reconstruction can achieve high quality MECT images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Differential Die-Away Self-Interrogation for Treaty Verification Applications (Final Report)

In this report, we summarize our experimental and simulation study of Differential Die-Away Self-Interrogation (DDSI) method used to assay 4 PuO 2 samples of almost identical isotopic composition with total Pu content ranging from 60-614g. With measured and calculated traditional neutron coincidence rates such as singles (S), doubles (D) and triples (T) we explore three different proposed DDSI signatures (D/S, T/S and D L /D F ) over a wide range of gating scenarios and analyze the sensitivity of individual observables to fissile content of the sample. We find that only the D L /D F observable is exclusively sensitive to the fissile content of sample, while D/S and T/S yield sensitivity to both fertile and fissile component of the sample. Based on the results of the Monte Carlo simulations, depending on the mass of the PuO 2 sample and considering 900s assay time with the Epithermal Neutron Multiplicity Counter (ENMC) we are able to reconstruct the fissile content with statistical uncertainty of 0.4-3.7% in case of D/S ratio, 1.1-4.3% in case of T/S ratio, and 4.9-34.3% in the case of the D L /D F .

Instrumentation Related to Nuclear Science & Techn↗
DOE OSTI2014

Scientific Discovery through Advanced Computing

My work during the past year has involved application and development of the M3D code. M3D (Multi-level 3D) carries out simulation studies of plasmas using multiple levels of physics, geometry, and grid models. M3D combines a two dimensional unstructured mesh with finite element discretization in poloidal planes, with a pseudo spectral representation in the toroidal direction. Two critical MHD issues have emerged with the ITER effort: the suppression of edge localized modes (ELMs) by resonant magnetic perturbations (RMPs) and disruptions. Large ELMs can cause an unacceptable heat and particle flux to the wall and divertor. The application of RMPs to ITER is now being planned, although there are still some questions about how RMPs work. Disruptions can produce a much larger heat and particle than ELMs, and in addition can cause unacceptable mechanical stress on the wall and blanket. I have been carrying out extended MHD simulations which address these issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

MCP-PMT studies at the High-B test facility at Jefferson Lab

Here we present preliminary results for the gain performance of commercially available 3- mum and 6- mum pore-size single-anode microchannel-plate photomultipliers (MCP PMTs) in magnetic fields up to 5 T and for various orientations of the sensor relative to the field direction. The measurements were performed at Thomas Jefferson National Accelerator Facility in Newport News, VA. Our results show that smaller-pore-size PMTs have better gain performance in magnetic fields. At various angles, the shape of the gain dependence on the strength of the magnetic field strongly depends on the type of the sensor. Also, for each sensor, the azimuthal dependence is strongly correlated with the polar angle. Overall, the sensors exhibit a reasonable performance up to 2 T, although that upper limit depends on the sensor, the applied high voltage, and the orientation of the sensor relative to the field. To optimize the operational and design parameters of MCP PMTs for performance in high magnetic fields, further measurements and simulation studies will be pursued. Furthermore, our studies are part of an R&D for development of a Detector of Internally Reflected Cherenkov Light for the central detector of a future U.S. Electron Ion Collider.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Investigation of powertrain system decarbonization using electrically assisted turbocharging and hybridization in off-road vehicles

In recent years, the causes, effects, and existential threat of a global anthropogenic climate shift have drawn significant attention and stimulated mitigation efforts from all genres of scientific, political, and industrial bodies. The greenhouse effect and the detrimental environmental impact of excess greenhouse gas (GHG) emissions, like carbon dioxide, is well studied and targeted as the primary culprit for aversive action. However, some fields exist where reducing GHG emissions is met by considerable challenges. One such field is the production and operations of off-road vehicles. Applications of these vehicles are highly diverse and are often characterized by rugged, especially transient, and power intensive duty cycles that make one-fits-all vehicle configuration solutions impractical. This research surveys existing literature to identify the modern technologies packages and challenges that face development and configuration of powertrain systems which currently navigate the regulated emissions and unique duty cycle requirements of this market space. A novel powertrain concept is proposed and evaluated with respect to the pinnacle objectives of load performance improvement, criteria pollutant reduction, and reducing the life cycle GHG intensity of its operation. A sure-fire path to vehicular GHG reduction is through improving the fuel efficiency of internal combustion engines (ICEs), an entrenched component of the off-road vehicle sector. Unfortunately, this is more easily said than done. An approach that has proven successful in this endeavor is engine downsizing and turbocharging, where a larger engine is replaced with a smaller one with added air system boosting via turbocharger to reduce frictional and pumping losses while also enabling access to additional fuel energy. However, practical realization of these potential benefits is often impeded by the transient response capability of the smaller engine across the operating space. For this reason, downsized ICE powertrains have turned to electrified forced induction systems (EFISs) for a decoupling of exhaust energy and engine speed from boost capability. Platformed on 48V hybrid technology, these systems introduce the need for more sophisticated controls around the engine gas-exchange process for the management of boost performance and exhaust gas emissions. On this account, simulation studies utilizing a GT-POWER model of a turbocharged 4.5 L engine outfit with an EFIS using an electrically driven compressor (eBooster®) are conducted to provide insight into the performance of this air handling architecture on an off-road engine. The results show that improvement in transient torque response time in sync with reduction in engine-out soot and NOx emissions are possible with an engine recalibration that leverages the transient air-fuel ratio authority of the eBooster®. Benefits are further demonstrated when duty cycle simulations of the 48V mild-hybrid engine concept are exercised, showing an acute decrease in cumulative fuel usage and soot production. Powertrain hybridization is another technological pathway achieving pronounced GHG reduction successes in modern on-road vehicle applications through integration of Li-ion battery technology. In the off-road vehicle segment, a review of available literature concludes that hybridized architectures are present but generally lack the depth of technologies that have both high specific energy and power capabilities, and thus are limited in their inclusion of Li-ion batteries for ICE assistance and enhanced energy storage capability. Therefore, building on the mild-hybrid engine results, the downsized and eBoosted engine concept was integrated into a larger high-voltage battery-hybrid series-electric powertrain system. Hybrid powertrain parameter sensitivity studies were carried out in a numerical charge-sustaining framework, providing novel insights into power flows between the battery and the engine and how their respective capabilities and operation contribute to GHG and criteria pollutant emissions of diverse duty cycles. Application of supervisory power management introduced robustness into the power sourcing and battery SOC control process and showed that optimum specification of battery properties can yield synergies between GHG emissions and battery energy capacity. Furthermore, examination of recent literature on Li-ion batteries has shown that pack manufacturing is a highly energy intensive process, thereby producing considerable quantities of GHGs that scale with energy storage capacity. Also scaling with a battery’s energy storage capacity is its investment cost. In consideration of these factors, an inclusive technoeconomic and GHG life cycle analysis is conducted. This analysis systematically compares the carbon footprint and total cost of ownership associated with the proposed hybrid powertrain concept to reference and alternative powertrain configurations, facilitating a thorough evaluation of the decarbonization effectiveness and economic viability.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2013

Investigating Short Circuit Models for Wind Turbine Generators

Short circuit study needs to be performed at the planning stage before connecting large gen erators with the electric grid. This study will reveal the impact of new generation on the fault levels of the grid. The fault levels will determine the rating of nearby equipment and setting of relays. However, a short circuit model of the new generation is required for this study. This report investigates and quantifies the short circuit models of different wind turbine generators (WTGs) and of wind farms, and using these models, explores the effects of the wind farms on the nearby switchgear. WTGs, commonly classified as Type 1, Type 2, Type 3, and Type 4, exhibit different short circuit behavior depending upon their architecture and controls. The short circuit model (synonymous to Thevenin model) of Type 1 WTG is well known; it comprises of a voltage source behind a transient reactance, X'. The transient reactance X' is derived using the reactance of stator, rotor and the magnetizing reactance. In case of Type 2 WTG, an external rotor resistance is used to increase the range of operating slips. The value of the external rotor resistance increases at higher operating slips. The effect of the external rotor resistance on the fault current response is analyzed by simulating a validated model of Type 2 WTG. The validation is performed using the field data from an actual 67-machine Type 2 wind farm located in Wyoming, USA. The results shows that the external rotor resistance should be considered in its Thevenin impedance for better accuracy. A wind farm contains a large number of WTGs, pad mounted transformers, and cables typically connected in a ‘daisy chain” configuration. The effect of these cable impedances on the short circuit model of Type 2 wind farm is examined using the validated models. It is found that the collector circuit cables inside the farm do not have a significant effect on the short circuit model of the farm. Based on this finding, a simple equivalencing procedure for the wind farm is proposed and validated. Faults inside and outside the wind farm are analyzed and recommendations to select main circuit breakers and generator circuit breakers are made. Type 3 WTG, also known as Doubly Fed Asynchronous Generator (DFAG), uses back to back converters in the rotor circuit to inject slip frequency voltage into the rotor circuit. Until recently, these machines used crowbar to short the rotor circuit during faults to protect the converter com ponents form high rotor currents resulting from faults. This type of WTG is simulated and the model is validated using field data from a 66-machine Type 3 wind farm in Wyoming, USA. The simulation study reveals that the crowbar feature is activated almost immediately after the fault takes place, and remains activated throughout the fault-period. This confirms that the short circuit model for Type 3 WTG with crowbar is similar to the model for Type 2 WTG. In the past, wind farms were allowed to disconnect during faults. However, recent grid codes around the world and in the United States demand that the wind farms be kept connected to the grid during faults, and provide certain reactive power quantified based on the voltage dip at the point of interconnection. This requirement, also called low voltage ride through (LVRT), can be accomplished by using Type 4 WTGs that employ full back to back converter interface between generators (permanent magnet synchronous generators are typically used) and the grid. However, the controls are proprietary, and hence the short circuit models are not well defined. This report proposes controls that enable the WTG to comply with the strictest grid codes, and based on this model, proposes a short circuit model of the Type 4 WTG. This model is not validated due to lack of validation data, but the use of grid code requirements as control objectives lends good credibility to the model. PSCAD/EMTDC ® , a powerful and time-tested time-domain simulation software is used to implement all models. This software is capable of modeling in detail all components of wind farms, and all components and controls of converters in time domain.

17 WIND ENERGY↗
DOE OSTI2018

Multicoordination Control Strategy Performance in Hybrid Power Systems

This paper evaluates a state-space methodology of a multi-input multi-output (MIMO) control strategy using a 2 × 2 tightly coupled scenario applied to a physical gas turbine fuel cell hybrid power system. A centralized MIMO controller was preferred compared to a decentralized control approach because previous simulation studies showed that the coupling effect identified during the simultaneous control of the turbine speed and cathode airflow was better minimized. The MIMO controller was developed using a state-space dynamic model of the system that was derived using first-order transfer functions empirically obtained through experimental tests. The controller performance was evaluated in terms of disturbance rejection through perturbations in the gas turbine operation, and setpoint tracking maneuver through turbine speed and cathode airflow steps. The experimental results illustrate that a multicoordination control strategy was able to mitigate the coupling of each actuator to each output during the simultaneous control of the system, and improved the overall system performance during transient conditions. On the other hand, the controller showed different performance during validation in simulation environment compared to validation in the physical facility, which will require a better dynamic modeling of the system for the implementation of future multivariable control strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

High-dimensional multivariate autoregressive model estimation of human electrophysiological data using fMRI priors

Multivariate autoregressive (MVAR) model estimation enables assessment of causal interactions in brain networks. However, accurately estimating MVAR models for high-dimensional electrophysiological recordings is challenging due to the extensive data requirements. Hence, the applicability of MVAR models for study of brain behavior over hundreds of recording sites has been very limited. Prior work has focused on different strategies for selecting a subset of important MVAR coefficients in the model to reduce the data requirements of conventional least-squares estimation algorithms. Here we propose incorporating prior information, such as resting state functional connectivity derived from functional magnetic resonance imaging, into MVAR model estimation using a weighted group least absolute shrinkage and selection operator (LASSO) regularization strategy. The proposed approach is shown to reduce data requirements by a factor of two relative to the recently proposed group LASSO method of Endemann et al (Neuroimage 254:119057, 2022) while resulting in models that are both more parsimonious and more accurate. The effectiveness of the method is demonstrated using simulation studies of physiologically realistic MVAR models derived from intracranial electroencephalography (iEEG) data. The robustness of the approach to deviations between the conditions under which the prior information and iEEG data is obtained is illustrated using models from data collected in different sleep stages. This approach allows accurate effective connectivity analyses over short time scales, facilitating investigations of causal interactions in the brain underlying perception and cognition during rapid transitions in behavioral state.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Shielded radiography with a laser-driven MeV-energy X-ray source

We report the results of experimental and numerical-simulation studies of shielded radiography using narrowband MeV-energy X-rays from a compact all-laser-driven inverse-Compton-scattering X-ray light source. This recently developed X-ray light source is based on a laser-wakefield accelerator with ultra-high-field gradient (GeV/cm). We demonstrate experimentally high-quality radiographic imaging (image contrast of 0.4 and signal-to-noise ratio of 2:1) of a target composed of 8-mm thick depleted uranium shielded by 80-mm thick steel, using a 6-MeV X-ray beam with a spread of 45% (FWHM) and 10 7 photons in a single shot. The corresponding dose of the X-ray pulse measured in front of the target is ~100 nGy/pulse. Simulations performed using the Monte-Carlo code MCNPX accurately reproduce the experimental results. These simulations also demonstrate that the narrow bandwidth of the Compton X-ray source operating at 6 and 9 MeV leads to a reduction of deposited dose as compared to broadband bremsstrahlung sources with the same end-point energy. The X-ray beam’s inherently low-divergence angle (~mrad) is advantageous and effective for interrogation at standoff distance. Furthermore, these results demonstrate significant benefits of all-laser driven Compton X-rays for shielded radiography.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
DOE OSTI2016

A promising limited angular computed tomography reconstruction via segmentation based regional enhancement and total variation minimization

X-ray computed tomography (CT) is a powerful and common inspection technique used for the industrial non-destructive testing. However, large-sized and heavily absorbing objects cause the formation of artifacts because of either the lack of specimen penetration in specific directions or the acquisition of data from only a limited angular range of views. Although the sparse optimization-based methods, such as the total variation (TV) minimization method, can suppress artifacts to some extent, reconstructing the images such that they converge to accurate values remains difficult because of the deficiency in continuous angular data and inconsistency in the projections. To address this problem, we use the idea of regional enhancement of the true values and suppression of the illusory artifacts outside the region to develop an efficient iterative algorithm. This algorithm is based on the combination of regional enhancement of the true values and TV minimization for the limited angular reconstruction. In this algorithm, the segmentation approach is introduced to distinguish the regions of different image knowledge and generate the support mask of the image. A new regularization term, which contains the support knowledge to enhance the true values of the image, is incorporated into the objective function. Then, the proposed optimization model is solved by variable splitting and the alternating direction method efficiently. A compensation approach is also designed to extract useful information from the initial projections and thus reduce false segmentation result and correct the segmentation support and the segmented image. The results obtained from comparing both simulation studies and real CT data set reconstructions indicate that the proposed algorithm generates a more accurate image than do the other reconstruction methods. The experimental results show that this algorithm can produce high-quality reconstructed images for the limited angular reconstruction and suppress the illusory artifacts caused by the deficiency in valid data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Iterative reconstruction using a Monte Carlo based system transfer matrix for dedicated breast positron emission tomography

To maximize sensitivity, it is desirable that ring Positron Emission Tomography (PET) systems dedicated for imaging the breast have a small bore. Unfortunately, due to parallax error this causes substantial degradation in spatial resolution for objects near the periphery of the breast. In this work, a framework for computing and incorporating an accurate system matrix into iterative reconstruction is presented in an effort to reduce spatial resolution degradation towards the periphery of the breast. The GATE Monte Carlo Simulation software was utilized to accurately model the system matrix for a breast PET system. A strategy for increasing the count statistics in the system matrix computation and for reducing the system element storage space was used by calculating only a subset of matrix elements and then estimating the rest of the elements by using the geometric symmetry of the cylindrical scanner. To implement this strategy, polar voxel basis functions were used to represent the object, resulting in a block-circulant system matrix. Simulation studies using a breast PET scanner model with ring geometry demonstrated improved contrast at 45% reduced noise level and 1.5 to 3 times resolution performance improvement when compared to MLEM reconstruction using a simple line-integral model. The GATE based system matrix reconstruction technique promises to improve resolution and noise performance and reduce image distortion at FOV periphery compared to line-integral based system matrix reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Cerebral perfusion computed tomography deconvolution via structure tensor total variation regularization

Purpose: Cerebral perfusion computed tomography (PCT) imaging as an accurate and fast acute ischemic stroke examination has been widely used in clinic. Meanwhile, a major drawback of PCT imaging is the high radiation dose due to its dynamic scan protocol. The purpose of this work is to develop a robust perfusion deconvolution approach via structure tensor total variation (STV) regularization (PD-STV) for estimating an accurate residue function in PCT imaging with the low-milliampere-seconds (low-mAs) data acquisition. Methods: Besides modeling the spatio-temporal structure information of PCT data, the STV regularization of the present PD-STV approach can utilize the higher order derivatives of the residue function to enhance denoising performance. To minimize the objective function, the authors propose an effective iterative algorithm with a shrinkage/thresholding scheme. A simulation study on a digital brain perfusion phantom and a clinical study on an old infarction patient were conducted to validate and evaluate the performance of the present PD-STV approach. Results: In the digital phantom study, visual inspection and quantitative metrics (i.e., the normalized mean square error, the peak signal-to-noise ratio, and the universal quality index) assessments demonstrated that the PD-STV approach outperformed other existing approaches in terms of the performance of noise-induced artifacts reduction and accurate perfusion hemodynamic maps (PHM) estimation. In the patient data study, the present PD-STV approach could yield accurate PHM estimation with several noticeable gains over other existing approaches in terms of visual inspection and correlation analysis. Conclusions: This study demonstrated the feasibility and efficacy of the present PD-STV approach in utilizing STV regularization to improve the accuracy of residue function estimation of cerebral PCT imaging in the case of low-mAs.

60 APPLIED LIFE SCIENCES↗