Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Resolving local structural motifs across the phase evolution of zinc titanates with computational x-ray absorption spectroscopy

Resolving the local structure motifs that characterize phase evolution as a function of composition is a key challenge in structure characterization of complex materials. Here, in this study, we combine first-principles simulations and x-ray absorption near-edge structures (XANES) analysis to gain insights into the structure evolution revealed by measurements across a combinatorial zinc titanate thin film, which was grown with smoothly varying composition over a wide range of the Ti:Zn ratio. Specifically, we propose a cluster blind-signal-separation (cBSS) method for XANES spectral analysis based on a library of the structures and spectra of representative local motifs. In addition to motifs from zinc titanate crystals, two types of Ti-defect models constructed in this study are key to the understanding of the structure characteristics in the Zn-rich region. The cBSS method makes use of both spectral clustering of the simulated site-XANES spectra library and the BSS procedure to construct high-fidelity spectral basis functions from an experimental spectral sequence. The method provides a rigorous measure of the spectral sensitivity and basis completeness. The results of the XANES analysis are corroborated with other experimental modalities, including x-ray diffraction and spectroscopic ellipsometry, to validate the cBSS method. The calculated motif weights resulting from fitting the XANES spectra with the cBSS basis probe the atomic structure characteristics of both crystalline and amorphous phases as a function of the Ti/Zn composition. The insights of the local structure motif evolution are pivotal to the understanding of the nonmonotonic trend in the optical gap, which may lead to potential applications through tuning the optical properties of zinc titanate. The workflow of the XANES spectral analysis developed in this work can be generalized to construct the structure-property relationship in a broad material space.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Development of the Advanced D1S for Shutdown Dose Rate Calculations in Fusion Reactors

The assessment of the shutdown dose rate due to neutron activation is fundamental for shielding design, materials selection, licensing and planning of maintenance operations in high performances fusion devices. The shutdown dose rate calculations require the combined use of radiation transports and inventory codes. The typical approach consists in a two steps procedure. A pure neutron transport simulation is performed in a first step to calculate the spatial distribution of neutron spectra. The results are used as input for an activation calculation that provides the decay gamma source used in a second step for a pure gamma transport simulation (rigorous-two step, R2S). The Direct 1-Step (D1S) approach has been developed more than ten years ago by ENEA and the ITER team: this technique relies on a coupled neutron-decay gamma transport in a single simulation. It is based on the use of a modified version of the MCNP5 Monte Carlo code with specially tailored nuclear cross-section data for transport and FISPACT inventory code for temporal behaviour. The 'Advanced-D1S' is an improved version of D1S in which additional computational features have been introduced, such as the mesh tally capability and continuous time evolution assessment of the dose rate in single run dose. It has been validated through measurements at JET tokamak and comparison with R2S codes and a future benchmark experiment during next high performance deuterium-tritium campaign is in preparation at JET. It has also been extensively used for ITER calculations and recently applied to DEMO assessment. Present limitations and in progress code improvement strategy are described and discussed in this paper.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

A Continuous-Discontinuous Hybrid Finite Element Method for S{sub N} Transport

Particle transport phenomenon for neutral particles, such as neutrons and photons, is governed by the linear Boltzmann equation, which is a first-order partial differential equation (PDE). In order to solve the transport equation, discontinuous finite element method (DFEM) is heavily developed. However, the large memory consumption due to the extra degrees of freedom is problematic in multi-D calculations. Another path is to recast the transport equation to a second-order form, such as the well-known self-adjoint angular flux equation (SAAF) and the even parity equations, and spatially discretize the equation using continuous finite element method (CFEM). Yet, such forms generally have problems in void situation as the streaming operator contains a 1/σ{sub t} term and becomes singular in void, which draws numerous efforts on developing void-treatment methods. An intuition to mitigate the shortcomings of both methods is then to develop a method that utilizes DFEM in part of the problem where CFEM has issues while in the other part, CFEM is employed to avoid DFEM's high cost. In this work, we present a novel continuous-discontinuous finite element method (CDFEM) for solving the discrete-ordinates (S{sub N}) transport equations. We perform the hybridization between 1/σ{sub t}-weighted least-squares finite element methods (1/σ{sub t}-LSFEM), with continuous basis functions, and first-order solver. Different solvers are defined in different regions, or subdomains, of the problem. The resulting scheme is such that CFEM-SAAF is performed in subdomains with continuous basis functions while upwinding based DFEM solver is utilized in subdomains with discontinuous basis functions. The connection between subdomains is through an upwinding operation on the interface, which can be performed via the sweep procedure. A foreseeable benefit is such a scheme is compatible with void as first-order DFEM solver can be used in void. On the other hand, as DFEM is only used in part of the problem, CDFEM can effectively avoid the notorious memory consumption of DFEM in multi-D calculation for the future work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Bulk-solvent and overall scaling revisited: faster calculations, improved results

A fast and robust method for determining the parameters for a flat (mask-based) bulk-solvent model and overall scaling in macromolecular crystallographic structure refinement and other related calculations is described. This method uses analytical expressions for the determination of optimal values for various scale factors. The new approach was tested using nearly all entries in the PDB for which experimental structure factors are available. In general, the resulting R factors are improved compared with previously implemented approaches. In addition, the new procedure is two orders of magnitude faster, which has a significant impact on the overall runtime of refinement and other applications. An alternative function is also proposed for scaling the bulk-solvent model and it is shown that it outperforms the conventional exponential function. Similarly, alternative methods are presented for anisotropic scaling and their performance is analyzed. All methods are implemented in the Computational Crystallography Toolbox (cctbx) and are used in PHENIX programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Tuning chemical precompression: Theoretical design and crystal chemistry of novel hydrides in the quest for warm and light superconductivity at ambient pressures

Over the past decade, a combination of crystal structure prediction techniques and experimental synthetic work has thoroughly explored the phase diagrams of binary hydrides under pressure. The fruitfulness of this dual approach is demonstrated in the recent identification of several superconducting hydrides with T c s approaching room temperature. We start with an overview of the computational procedures for predicting stable structures and estimating their propensity for superconductivity. A survey of phases with high T c reveals some common structural features that appear conducive to the strong coupling of the electronic structure with atomic vibrations that leads to superconductivity. We discuss the stability and superconducting properties of phases containing two of these—molecular H 2 units mixed with atomic H and hydrogenic clathrate-like cages—as well as more unique motifs. Lastly, we argue that ternary hydride phases, whose exploration is still in its infancy, are a promising route to achieve simultaneous superconductivity at high temperatures and stability at low pressures. Several ternary hydrides arise from the addition of a third element to a known binary hydride structure through site mixing or onto a new site, and several more are based on altogether new structural motifs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Modern sensitivity analysis - 15651

Sensitivity analysis has been applied to complex models, such as performance assessment models, since at least the advent of computational tools to support such an approach to model evaluation. Initially, sensitivity analysis methods were applied to deterministic models and focused on changing one input variable at a time. Probabilistic modeling opened the door to evaluating changes in all input variables simultaneously. Standardized rank regression methods were primarily used in the last century to evaluate probabilistic models. These methods assume linearity in the rank response of the model to the inputs. At the very least this assumes monotonicity in the response, but for the most part the response needs to be near-linear. This is highly unlikely for performance assessment models that are based on differential equations, multiplication and division of variables that are usually not correlated, and conservatism that is not well understood. There are better statistical methods for evaluating highly non-linear and dynamic probabilistic models. Some of these methods were founded in the data mining industry. This century has seen a rapid development of the internet that led directly to the development of data mining techniques for massive databases. Although the goals were different, these techniques can be adapted to sensitivity analysis for complex probabilistic models, which, also, form large, or even massive, databases. These modern sensitivity analysis procedures assess the importance of every explanatory variable (input parameters) used in a probabilistic model. In practice, this means that for each endpoint of interest, every explanatory variable in the model is associated with a quantitative measure of importance. For a given explanatory variable, the importance value depends on the endpoint of interest. This global approach to sensitivity analysis allows all levels of interactions to be evaluated, so that the variation in the output can be measured at every level of possible interaction. The effects are collected together during sensitivity analysis processing, and then separated out to attribute the overall contribution of each input parameter to the model output. These global sensitivity analysis procedures have the ability to characterize non-linear and non-monotonic relationships between explanatory variables and endpoints of interest. This is critically important because of the need to characterize complex interactions among multiple explanatory variables in PA models that often have non-linear and non-monotonic effects. The options for sensitivity analysis are fully explored, including an evaluation of several methods when applying them to a well-known complex non-linear non-monotonic model. The results clearly demonstrate the relative benefit of using data mining methods over simple statistical approaches such as standardized rank regression. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Nanoindentation mapping defects filtration for heterogeneous materials using generative adversarial networks

Advanced composite materials with multiple phases and heterogeneous microstructure necessitate spatial mapping characterization of elastic modulus to develop constitutive relations and overall mechanical response. Such modulus mapping can be obtained using the nanoindentation technique, where the indenter tip raster over the selected microstructure region. Typically, a surface preparation procedure is done in the specimens to ensure proper contact between the indenter tip and sample surface. However, a near-perfect surface finish is unachievable in heterogeneous materials, primarily with ceramic reinforcements, due to the differential material removal rate during polishing. Thus, the nanoindenter records localized erroneous measurements due to differences in surface roughness and corresponding force response. This study establishes a novel deep learning-based strategy to rectify incorrect experimental spatial measurements acquire during nanoindentation modulus mapping. Here, the integrated bicubic interpolation and generative adversarial networks (GANs) model was trained using 14 ceramic and 18 metallic data sets, each comprising 65,536 measurements. The developed algorithm was validated against experimental measurements on four unknown specimens. The standard deviation in measured elastic modulus reduces by ~50% in ceramics and ~72% in metallic samples. This computational framework proposes a novel approach to reducing uncertainty in materials’ properties using state-of-the-art computer vision techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A physics-based solver to optimize the illumination of cylindrical targets in spherically distributed high power laser systems

In recent years, our understanding of high energy density plasmas has played an important role in improving inertial fusion confinement and in emerging new fields of physics, such as laboratory astrophysics. Every new idea required developing innovative experimental platforms at high power laser facilities, such as OMEGA or NIF. These facilities, designed to focus all their beams onto spherical targets or hohlraum windows, are now required to shine them on more complex targets. While the pointing on planar geometries is relatively straightforward, it becomes problematic for cylindrical targets or target with more complex geometries. This publication describes how the distribution of laser beams on a cylindrical target can be done simply by using a set of physical laws as a pointing procedure. The advantage of the method is threefold. First, it is straightforward, requiring no mathematical enterprise besides solving ordinary differential equations. Second, it will converge if a local optimum exists. Finally, it is computationally inexpensive. Experimental results show that this approach produces a geometrical beam distribution that yields cylindrically symmetric implosions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

The coyote universe extended: Precision emulation of the matter power spectrum

Modern sky surveys are returning precision measurements of cosmological statistics such as weak lensing shear correlations, the distribution of galaxies, and cluster abundance. To fully exploit these observations, theorists must provide predictions that are at least as accurate as the measurements, as well as robust estimates of systematic errors that are inherent to the modeling process. In the nonlinear regime of structure formation, this challenge can only be overcome by developing a large-scale, multi-physics simulation capability covering a range of cosmological models and astrophysical processes. As a first step to achieving this goal, we have recently developed a prediction scheme for the matter power spectrum (a so-called emulator), accurate at the 1% level out to k ∼ 1 Mpc{sup –1} and z = 1 for wCDM cosmologies based on a set of high-accuracy N-body simulations. It is highly desirable to increase the range in both redshift and wavenumber and to extend the reach in cosmological parameter space. To make progress in this direction, while minimizing computational cost, we present a strategy that maximally reuses the original simulations. We demonstrate improvement over the original spatial dynamic range by an order of magnitude, reaching k ∼ 10 h Mpc{sup –1}, a four-fold increase in redshift coverage, to z = 4, and now include the Hubble parameter as a new independent variable. To further the range in k and z, a new set of nested simulations run at modest cost is added to the original set. The extension in h is performed by including perturbation theory results within a multi-scale procedure for building the emulator. This economical methodology still gives excellent error control, ∼5% near the edges of the domain of applicability of the emulator. A public domain code for the new emulator is released as part of the work presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-T-173: Evaluation of the MicroSTAR Reader for the OSLD System for Remote Dosimetry Audits at the IROC HQAC

Purpose: Analysis of the performance of the microSTARii reader for optically stimulated luminescence dosimeters (OSLD) used by the IROC Houston Quality Assurance Center (IROC HQAC) for external audits compare to the characteristics of the Microstar reader. Methods: The IROC HQAC uses the Microstar reader for the OSLD program for verification of output of photon, electron and proton beams. The calculation of dose from the OSLD system is based on a group of factors defined at time of the commissioning of a batch of detectors. Factors like system sensitivity (SS), depletion (KD), element correction factor (ECF), linearity (KL), energy correction (KE). The new microSTARii unit presents many hardware and software upgrades that were considered useful for this program. Based on these changes many factors, that were considered reader dependent, were revised in order to analyze the effect of the new reading process. The SS, KD, ECF and KL were evaluated and compared with data defined based on reading done on Microstar reader. Results: The SS is reader specific and specified at 100 cGy dose level. This value is define per reading session and monitored over time. The KD factor was found to be different because of reading procedure are different. The ECF values changed for a group of nano dots compare to values defined based on reading done on the Microstar reader. The KL was defined for the reader. Conclusion: The new microSTARii reader presents new features that improve the efficiency of the OSLD program at the IROC HQAC. New characterization is needed before final implementation is done to match the requirements of the existing OSLD system defined for the Microstar reader. Changes in uncertainty of the results has not been analized.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Performance Testing of a Moving-Bed Gasifier Using Coal, Biomass, and Waste Plastic Blends with Washed and Unwashed Legacy Coals and Other Waste Fuels to Generate White Hydrogen

The objective of this effort, primarily funded by the United States Department of Energy (DOE), and led by the Electric Power Research Institute, Inc. (EPRI), with support by Hamilton Maurer International (HMI) and Sotacarbo S.p.A. (Sotacarbo), has been to qualify coal, biomass, and plastic waste blends based on performance testing of selected fuel pellet compositions in a pilot-scale updraft moving-bed (UDMB) gasifier. The testing provided relevant data to advance the commercial-scale design of the moving-bed gasifier to be able to successfully use these feedstocks to produce hydrogen. In particular, the effects of waste plastics on feedstock development (i.e., blending and pelletizing) and the resulting products (i.e., syngas compositions, organic condensate production, and ash characteristics) are the focus. The gasifier used for testing is HMI’s moving-bed gasifier, which has been proven capable of gasifying nearly all coal ranks. It has also shown the ability in prior testing work to gasify wood chips (biomass). However, mixtures of these fuels with plastic wastes have not been prepared and gasified together. The three feedstocks were densified and pelletized by California Pellet Mill (CPM) to meet the feedstock size required by Sotacarbo’s 30mm ID UDMB gasifier, under contract to HMI. The technical tasks and results from this two-year research project included: (1) Feed Procurement and Preparation: Nine different tri-fuel pellets were prepared from varying compositions of fresh mined PRB coal, corn stover biomass, and car fluff waste plastics. Tri-fuel pellets were produced by CPM and shipped to Sotacarbo’s test facility in Carbonia, Sardinia, Italy. (2) Test Plan Development: A test plan was created to define the test runs to be performed. The test plan detailed the different UDMB gasification tests to be performed in Sotacarbo’s 12-inch ID pilot scale gasifier, the process monitoring instrumentation used, and the extractive samples recovered for analysis of the total gasification process mass and energy balance. (3) Gasifier Testing: Nine different gasification runs were performed in the pilot-scale gasifier at Sotacarbo using nine different fuel feedstock compositions generated from varying mixtures of PRB coal, biomass, and plastic wastes. The testing generated performance data on gasification reaction efficiency and performance, yielding relevant data for models used to scale up the gasifier design. This task also included work to refurbish and reassemble the pilot gasifier at Sotacarbo and perform a baseline 100% PRB coal run. (4) Data Analysis and Reporting: Review of the data, determination of figures of merit, and interpretation of the results are reported in the project’s final report, published in March 2024. The results show that all tri-fuel pellets gasified well and maintained structural integrity throughout the gasification process. The syngas generated can be shifted to hydrogen by using commercial syngas shifting technologies. (5) High Fidelity computational fluid dynamics (CFD) Simulation: The National Energy Technology Laboratory (NETL) team performed CFD simulations of the UDMB gasifier for two of the tri-fuel pellets gasified in Sotacarbo’s pilot scale gasifier. The kinetic mechanisms for the pyrolysis of each constituent, PRB coal, corn stover biomass, and waste plastics are based on thermogravimetric analysis performed by Sotacarbo. The gasification model was validated by comparing the predicted syngas composition at the exit of the gasifier with the measured syngas composition. In addition, the reactor’s measured internal temperature profile agreed well with the predicted internal reactor temperature profile. These results validate that the model can be used to predict the performance of the updraft moving bed gasifier for different feedstocks and operating conditions. This paper summarizes the results of the completed work in which the pelletizing procedure was validated to ensure the viability of the tri-fuel pellets for the gasification runs performed at Sotacarbo’s 30 mm UDMB gasifier. The gasification performance data from this series of nine runs will enable modeling of a full-scale HMI industrial scale gasifier supporting both combined heat and power, and Hydrogen production from coal (both fresh mined and legacy) combined with various biomass and waste plastics. Additionally, plans and progress on a follow-up project, being executed by the same project team, will be presented. In this project, a total of twenty (20) different feedstocks are being prepared from varying compositions of biomass (both woody biomass and corn stover) with a mixture of legacy coal waste, plastic waste, and refuse-derived fuel (RDF). The testing will provide information on gasification reaction efficiency/performance, yielding relevant data for models used to scale up the gasifier design to 50 megawatt electric (MWe) (equivalent hydrogen production). Tests will also be performed on a bench-scale fluidized-bed gasifier for comparison purposes. The results of this testing will be used to specify the range of feedstock blends that can be successfully gasified as well as quantify gasifier outputs based on specific blends.

08 HYDROGEN↗
DOE OSTI2014

TH-C-BRD-05: Reducing Proton Beam Range Uncertainty with Patient-Specific CT HU to RSP Calibrations Based On Single-Detector Proton Radiography

Purpose: Beam range uncertainty in proton treatment comes primarily from converting the patient's X-ray CT (xCT) dataset to relative stopping power (RSP). Current practices use a single curve for this conversion, produced by a stoichiometric calibration based on tissue composition data for average, healthy, adult humans, but not for the individual in question. Proton radiographs produce water-equivalent path length (WEPL) maps, dependent on the RSP of tissues within the specific patient. This work investigates the use of such WEPL maps to optimize patient-specific calibration curves for reducing beam range uncertainty. Methods: The optimization procedure works on the principle of minimizing the difference between the known WEPL map, obtained from a proton radiograph, and a digitally-reconstructed WEPL map (DRWM) through an RSP dataset, by altering the calibration curve that is used to convert the xCT into an RSP dataset. DRWMs were produced with Plastimatch, an in-house developed software, and an optimization procedure was implemented in Matlab. Tests were made on a range of systems including simulated datasets with computed WEPL maps and phantoms (anthropomorphic and real biological tissue) with WEPL maps measured by single detector proton radiography. Results: For the simulated datasets, the optimizer showed excellent results. It was able to either completely eradicate or significantly reduce the root-mean-square-error (RMSE) in the WEPL for the homogeneous phantoms (to zero for individual materials or from 1.5% to 0.2% for the simultaneous optimization of multiple materials). For the heterogeneous phantom the RMSE was reduced from 1.9% to 0.3%. Conclusion: An optimization procedure has been designed to produce patient-specific calibration curves. Test results on a range of systems with different complexities and sizes have been promising for accurate beam range control in patients. This project was funded equally by the Engineering and Physical Sciences Research Council (UK) and Ion Beam Applications (Louvain-La-Neuve, Belgium)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

Finding a maximal element of a non-negative convex set through its characteristic cone: an application to finding a strictly complementary solution

In order to express a polyhedron as the Minkowski sum of a polytope and a polyhedral cone, Motzkin (Beiträge zur Theorie der linearen Ungleichungen. Dissertation, University of Basel, 1936) devised a homogenization technique that translates the polyhedron to a polyhedral cone in one higher dimension. Refining his technique, we present a conical representation of a set in the Euclidean space. Then, we use this representation to reach four main results: First, we establish a convex programming based framework for determining a maximal element—an element with the maximum number of positive components—of a non-negative convex set—a convex set in the non-negative Euclidean orthant. Second, we develop a linear programming problem for finding a relative interior point of a polyhedron. Third, we propose two procedures for identifying a strictly complementary solution in linear programming. Finally, we generalize Motzkin’s (Beiträge zur Theorie der linearen Ungleichungen. Dissertation, University of Basel, 1936) representation theorem for a class of closed convex sets in the Euclidean space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Dead-Time Correction of Pulsed Neutron Source Experiments with the Backward Extrapolation Method

This paper introduces a novel procedure to correct pulsed neutron source experimental data for dead-time effects. The procedure is based on the application of the backward extrapolation method to the reaction rate value in each time interval (bin) of the neutron counts distribution as a function of time. The theoretical framework has been successfully applied to the experimental data (timestamps) from the YALINA Thermal facility driven by a pulsed D-D neutron source. Dead-time effects may have more than 500 pcm (1 pcm=1 E-5) impact on the effective multiplication factor calculated by the area method. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

SU-F-T-47: MRI T2 Exclusive Based Planning Using the Endocavitary/interstitial Gynecological Benidorm Applicator: A Proposed TPS Library and Preplan Efficient Methodology

Purpose: ABS and GEC-ESTRO have recommended MRI T2 for image guided brachytherapy. Recently, a new applicator (Benidorm Template, TB) has been developed in our Department (Rodriguez et al 2015). TB is fully MRI compatible because the Titanium needles and it allows the use of intrauterine tandem. Currently, TPS applicators library are not currently available for non-rigid applicators in case of interstitial component as the TB.The purpose of this work is to present the development of a library for the TB, together with its use on a pre-planning technique. Both new goals allow a very efficient and exclusive T2 MRI based planning clinical TB implementation. Methods: The developed library has been implemented in Oncentra Brachytherapy TPS, version 4.3.0 (Elekta) and now is being implemented on Sagiplan v 2.0 TPS (Eckert&Ziegler BEBIG). To model the TB, free and open software named FreeCAD and MeshLab have been used. The reconstruction process is based on three inserted A-vitamin pellets together with the data provided by the free length. The implemented preplanning procedure is as follow: 1) A MRI T2 acquisition is performed with the template in place just with the vaginal cylinder (no uterine tube nor needles). 2) The CTV is drawn and the required needles are selected using a developed Java based application and 3) A post-implant MRI T2 is performed. Results: This library procedure has been successfully applied by now in 25 patients. In this work the use of the developed library will be illustrated with clinical examples. The preplanning procedure has been applied by now in 6 patients, having significant advantages: needle depth estimation, needle positions and number are optimized a priori, time saving, etc Conclusion: TB library and pre-plan techniques are feasible and very efficient and their use will be illustrated in this work.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗
DOE OSTI2016

Poster - 32: Atlas Selection for Automated Segmentation of Pelvic CT for Prostate Radiotherapy

Atlas based-segmentation has recently been evaluated for use in prostate radiotherapy. In a typical approach, the essential step is the selection of an atlas from a database that the best matches of the target image. This work proposes an atlas selection strategy and evaluate it impacts on final segmentation accuracy. Several anatomical parameters were measured to indicate the overall prostate and body shape, all of these measurements obtained on CT images. A brute force procedure was first performed for a training dataset of 20 patients using image registration to pair subject with similar contours; each subject was served as a target image to which all reaming 19 images were affinity registered. The overlap between the prostate and femoral heads was quantified for each pair using the Dice Similarity Coefficient (DSC). Finally, an atlas selection procedure was designed; relying on the computation of a similarity score defined as a weighted sum of differences between the target and atlas subject anatomical measurement. The algorithm ability to predict the most similar atlas was excellent, achieving mean DSCs of 0.78 ± 0.07 and 0.90 ± 0.02 for the CTV and either femoral head. The proposed atlas selection yielded 0.72 ± 0.11 and 0.87 ± 0.03 for CTV and either femoral head. The DSC obtained with the proposed selection method were slightly lower than the maximum established using brute force, but this does not include potential improvements expected with deformable registration. The proposed atlas selection method provides reasonable segmentation accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Metabolic alterations in bladder cancer: applications for cancer imaging

Treatment planning, outcome and prognosis are strongly related to the adequate tumor staging for bladder cancer (BC). Unfortunately, a large discrepancy exists between the preoperative clinical and final pathologic staging. Therefore, an advanced imaging-based technique is crucial for adequate staging. Although Magnetic Resonance Imaging (MRI) is currently the best in vivo imaging technique for BC staging because of its excellent soft-tissue contrast and absence of ionizing radiation it lacks cancer-specificity. Tumor-specific positron emission tomography (PET), which is based on the Warburg effect (preferential uptake of glucose by cancer cells), exploits the radioactively–labeled glucose analogs, i.e., FDG. Although FDG-PET is highly cancer specific, it lacks resolution and contrast quality comparable with MRI. Chemical Exchange Saturation Transfer (CEST) MRI enables the detection of low concentrations of metabolites containing protons. BC is an attractive target for glucose CEST MRI because, in addition to the typical systemic administration, glucose might also be directly applied into the bladder to reduce toxicity-related complications. As a first stage of the development of a contrast-specific BC imaging technique we have studied glucose uptake by bladder epithelial cells and have observed that glucose is, indeed, consumed by BC cells with higher intensity than by non-transformed urothelial cells. This effect might be partly explained by increased expression of glucose transporters GLUT1 and GLUT3 in transformed cells as compared to normal urothelium. We also detected higher lactate production by BC cells which is another cancer-specific manifestation of the Warburg effect. In addition, we have observed other metabolic alterations in BC cells as compared to non-transformed cells: in particular, increased pyruvate synthesis. When glucose was substituted by glutamine in culture media, preferential uptake of glutamine by BC cells was observed. The preferential uptake of glucose by BC cells gives an opportunity to develop NMR based imaging procedures where glucose or its derivatives can serve as a contrasting agent. In addition, metabolic alterations observed in BC cells could provide the basis for development of new anti-cancer therapeutics.

60 APPLIED LIFE SCIENCES↗