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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Accuracy of Length of Virtual Stents in Treatment of Intracranial Wide-Necked Aneurysms

Background and Purpose: Precise stent deployment is important for successful treatment of intracranial aneurysms by stent-assisted coiling (SAC). We evaluated the accuracy of virtual stents generated using commercial stent planning software by comparing the length of virtual and actually deployed intracranial laser cut stents on three-dimensional digital subtraction angiography (3D-DSA) images. Methods: We retrospectively analyzed the data of 75 consecutive cases of intracranial wide-necked aneurysms treated with the SAC technique using laser cut stents. Based on 3D-DSA images acquired by C-arm CT, stent sizing and placement were intraoperatively simulated by a commercial software application. The difference in length of the stents was estimated by measuring proximal discrepancies between the end points of the virtual and actually deployed stents on fused pre-procedural and post-procedural 3D-DSA images. Discrepancies between distal stent end points were manually minimized. The Kruskal–Wallis test was applied to test whether stent location, type, and length had an effect on difference in length between virtual and real stent. Results: The median difference in length between virtual and real stents was 1.58 mm with interquartile range 1.12–2.12 mm. There was no evidence for an effect of stent location (p = 0.23), stent type (p = 0.33), or stent length (p = 0.53) on difference in length between virtual and real stents. Conclusions: Stent planning software allows 3D simulation of laser cut stents overlain on 3D-DSA images of vessels and may thus be useful for stent selection and deployment of laser cut stents during stent-assisted coiling of intracranial aneurysms.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Analyzing simulation-based PRA data through traditional and topological clustering: A BWR station blackout case study

Here, dynamic probabilistic risk assessment (DPRA) methodologies couple system simulator codes (e.g., RELAP and MELCOR) with simulation controller codes (e.g., RAVEN and ADAPT). Whereas system simulator codes model system dynamics deterministically, simulation controller codes introduce both deterministic (e.g., system control logic and operating procedures) and stochastic (e.g., component failures and parameter uncertainties) elements into the simulation. Typically, a DPRA is performed by sampling values of a set of parameters and simulating the system behavior for that specific set of parameter values. For complex systems, a major challenge in using DPRA methodologies is to analyze the large number of scenarios generated, where clustering techniques are typically employed to better organize and interpret the data. In this paper, we focus on the analysis of two nuclear simulation datasets that are part of the risk-informed safety margin characterization (RISMC) boiling water reactor (BWR) station blackout (SBO) case study. We provide the domain experts a software tool that encodes traditional and topological clustering techniques within an interactive analysis and visualization environment, for understanding the structures of such high-dimensional nuclear simulation datasets. We demonstrate through our case study that both types of clustering techniques complement each other for enhanced structural understanding of the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Quality control for quantitative multicenter whole-body PET/MR studies: A NEMA image quality phantom study with three current PET/MR systems

Purpose: Integrated positron emission tomography/magnetic resonance (PET/MR) systems derive the PET attenuation correction (AC) from dedicated MR sequences. While MR-AC performs reasonably well in clinical patient imaging, it may fail for phantom-based quality control (QC). The authors assess the applicability of different protocols for PET QC in multicenter PET/MR imaging. Methods: The National Electrical Manufacturers Association NU 2 2007 image quality phantom was imaged on three combined PET/MR systems: a Philips Ingenuity TF PET/MR, a Siemens Biograph mMR, and a GE SIGNA PET/MR (prototype) system. The phantom was filled according to the EANM FDG-PET/CT guideline 1.0 and scanned for 5 min over 1 bed. Two MR-AC imaging protocols were tested: standard clinical procedures and a dedicated protocol for phantom tests. Depending on the system, the dedicated phantom protocol employs a two-class (water and air) segmentation of the MR data or a CT-based template. Differences in attenuation- and SUV recovery coefficients (RC) are reported. PET/CT-based simulations were performed to simulate the various artifacts seen in the AC maps (μ-map) and their impact on the accuracy of phantom-based QC. Results: Clinical MR-AC protocols caused substantial errors and artifacts in the AC maps, resulting in underestimations of the reconstructed PET activity of up to 27%, depending on the PET/MR system. Using dedicated phantom MR-AC protocols, PET bias was reduced to −8%. Mean and max SUV RC met EARL multicenter PET performance specifications for most contrast objects, but only when using the dedicated phantom protocol. Simulations confirmed the bias in experimental data to be caused by incorrect AC maps resulting from the use of clinical MR-AC protocols. Conclusions: Phantom-based quality control of PET/MR systems in a multicenter, multivendor setting may be performed with sufficient accuracy, but only when dedicated phantom acquisition and processing protocols are used for attenuation correction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Experimental evaluation of actual delivered dose using mega-voltage cone-beam CT and direct point dose measurement

Radiation therapy in patients is planned by using computed tomography (CT) images acquired before start of the treatment course. Here, tumor shrinkage or weight loss or both, which are common during the treatment course for patients with head-and-neck (H and N) cancer, causes unexpected differences from the plan, as well as dose uncertainty with the daily positional error of patients. For accurate clinical evaluation, it is essential to identify these anatomical changes and daily positional errors, as well as consequent dosimetric changes. To evaluate the actual delivered dose, the authors proposed direct dose measurement and dose calculation with mega-voltage cone-beam CT (MVCBCT). The purpose of the present study was to experimentally evaluate dose calculation by MVCBCT. Furthermore, actual delivered dose was evaluated directly with accurate phantom setup. Because MVCBCT has CT-number variation, even when the analyzed object has a uniform density, a specific and simple CT-number correction method was developed and applied for the H and N site of a RANDO phantom. Dose distributions were calculated with the corrected MVCBCT images of a cylindrical polymethyl methacrylate phantom. Treatment processes from planning to beam delivery were performed for the H and N site of the RANDO phantom. The image-guided radiation therapy procedure was utilized for the phantom setup to improve measurement reliability. The calculated dose in the RANDO phantom was compared to the measured dose obtained by metal-oxide-semiconductor field-effect transistor detectors. In the polymethyl methacrylate phantom, the calculated and measured doses agreed within about +3%. In the RANDO phantom, the dose difference was less than +5%. The calculated dose based on simulation-CT agreed with the measured dose within±3%, even in the region with a high dose gradient. The actual delivered dose was successfully determined by dose calculation with MVCBCT, and the point dose measurement with the image-guided radiation therapy procedure.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

A plastic miniature x-ray emission spectrometer based on the cylindrical von Hamos geometry

We present a short working distance miniature x-ray emission spectrometer (miniXES) based on the cylindrical von Hamos geometry. We describe the general design principles for the spectrometer and detail a specific implementation that covers K{beta} and valence level emission from Fe. Large spatial and angular access to the sample region provides compatibility with environmental chambers, microprobe, and pump/probe measurements. The primary spectrometer structure and optic is plastic, printed using a 3-dimensional rapid-prototype machine. The spectrometer is inexpensive to construct and is portable; it can be quickly set up at any focused beamline with a tunable narrow bandwidth monochromator. The sample clearance is over 27 mm, providing compatibility with a variety of environment chambers. An overview is also given of the calibration and data processing procedures, which are implemented by a multiplatform user-friendly software package. Finally, representative measurements are presented. Background levels are below the level of the K{beta}{sub 2,5} valence emission, the weakest diagram line in the system, and photometric analysis of count rates finds that the instrument is performing at the theoretical limit.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part II

X-ray Bragg coherent diffraction imaging (BCDI) has been demonstrated as a powerful 3D microscopy approach for the investigation of sub-micrometre-scale crystalline particles. The approach is based on the measurement of a series of coherent Bragg diffraction intensity patterns that are numerically inverted to retrieve an image of the spatial distribution of the relative phase and amplitude of the Bragg structure factor of the diffracting sample. This 3D information, which is collected through an angular rotation of the sample, is necessarily obtained in a non-orthogonal frame in Fourier space that must be eventually reconciled. To deal with this, the approach currently favored by practitioners (detailed in Part I) is to perform the entire inversion in conjugate non-orthogonal real- and Fourier-space frames, and to transform the 3D sample image into an orthogonal frame as a post-processing step for result analysis. In this article, which is a direct follow-up of Part I, two different transformation strategies are demonstrated, which enable the entire inversion procedure of the measured data set to be performed in an orthogonal frame. The new approaches described here build mathematical and numerical frameworks that apply to the cases of evenly and non-evenly sampled data along the direction of sample rotation ( i.e. the rocking curve). The value of these methods is that they rely on the experimental geometry, and they incorporate significantly more information about that geometry into the design of the phase-retrieval Fourier transformation than the strategy presented in Part I. Finally, two important outcomes are (1) that the resulting sample image is correctly interpreted in a shear-free frame and (2) physically realistic constraints of BCDI phase retrieval that are difficult to implement with current methods are easily incorporated. Computing scripts are also given to aid readers in the implementation of the proposed formalisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Optimization of oncological {sup 18}F-FDG PET/CT imaging based on a multiparameter analysis

Purpose: This paper describes a method to achieve consistent clinical image quality in {sup 18}F-FDG scans accounting for patient habitus, dose regimen, image acquisition, and processing techniques. Methods: Oncological PET/CT scan data for 58 subjects were evaluated retrospectively to derive analytical curves that predict image quality. Patient noise equivalent count rate and coefficient of variation (CV) were used as metrics in their analysis. Optimized acquisition protocols were identified and prospectively applied to 179 subjects. Results: The adoption of different schemes for three body mass ranges (<60 kg, 60–90 kg, >90 kg) allows improved image quality with both point spread function and ordered-subsets expectation maximization-3D reconstruction methods. The application of this methodology showed that CV improved significantly (p < 0.0001) in clinical practice. Conclusions: Consistent oncological PET/CT image quality on a high-performance scanner was achieved from an analysis of the relations existing between dose regimen, patient habitus, acquisition, and processing techniques. The proposed methodology may be used by PET/CT centers to develop protocols to standardize PET/CT imaging procedures and achieve better patient management and cost-effective operations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Status of Electronics Upgrades to the LANL Green is Clean Phoswich Detector Systems - 16419

Los Alamos National Laboratory radiological facilities produce low-density room trash that, in many cases, is not contaminated with radioactivity. It has been estimated that 50 to 90% of low-density room trash is free of radioactive contamination and eligible for inclusion in LANL's Green is Clean (GIC) program. The GIC program is a verification program for non-regulated waste from radiation controlled areas that has been actively segregated as nonradioactive through the use of the waste generator's acceptable knowledge. The GIC program has employed two large area Phoswich gamma ray detector systems (HERCULES and ZEUS) to assay and verify that GIC candidate waste meets free release disposal requirements. The large area (127 cm{sup 2}) Phoswich detectors consist of a thin NaI front crystal (3 mm thick) that is optically coupled to a thick CsI back crystal (50 mm thick). The crystals are set in an oxygen-free high-conductivity copper housing with a very thin (0.025 mm) aluminum entrance window. The scintillation properties of the two crystals are different enough to allow the system's electronic components to identify the origin of any pulse. The original systems utilize Nuclear Instrumentation Module (NIM) electronics to process the detector signals. These electronic modules provide the buffering, amplification, and pulse shape analyzer needed for the proper presentation of signals to the custom Phosmux router module. The router module separates the sodium iodide and cesium iodide signals by their rise times and transfers this information into the computer-mounted multichannel analyzer, a Canberra Model S-100 system. While these NIM Bin counting electronics configurations have been effective for approximately 20 years of operation, they are bulky requiring numerous modules per detector and the important Phosmux module is no longer manufactured by LANL. Modern advancement of commercial counting electronics now offer more compact options for equivalent detector signal processing. A pulse shape discriminator module, the DGF Pixie-4, made by XIA LLC, has been identified that can replace the majority of the original NIM bin electronics in a much more compact chassis - one Pixie-4 replaces four NIM amplifiers and four NIM PSAs and fits in a small National Instruments{sup TM} chassis. Thus, LANL has initiated a system upgrade project to replace the historic NIM electronics with the XIA Pixie4 and associated configuration software. The hardware for one Phoswich systems has been procured for the project and the configuration software, customized for the LANL Phoswich detector application, has been written by the manufacturer. Initial setup and optimization of the Pixie-4 configuration has been performed. Qualitative and quantitative baseline comparisons to the NIM electronics has been performed. Good agreement of key low energy gamma lines between the systems has allowed 'Proof of Concept' to be concluded. Further optimization and performance evaluations are in progress. Modifications to the LabVIEW{sup TM} based operating software (originally written by VI Control Systems Ltd.) are being performed by a LANL LabVIEW{sup TM} programmer to accommodate the electronics upgrade. A new LabVIEW{sup TM} interface to control the Pixie-4 has already been developed. Final optimization and testing, formal software QA, and revisions to operating procedures are in progress to achieve completion of the GIC system upgrade project. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

MO-DE-207B-05: Predicting Gene Mutations in Renal Cell Carcinoma Based On CT Imaging Features: Validation Using TCGA-TCIA Datasets

Purpose: The goal of this work is to investigate the use of contrast enhanced computed tomographic (CT) features for the prediction of mutations of BAP1, PBRM1, and VHL genes in renal cell carcinoma (RCC). Methods: For this study, we used two patient databases with renal cell carcinoma (RCC). The first one consisted of 33 patients from our institution (UT Southwestern Medical Center, UTSW). The second one consisted of 24 patients from the Cancer Imaging Archive (TCIA), where each patient is connected by a unique identi?er to the tissue samples from the Cancer Genome Atlas (TCGA). From the contrast enhanced CT image of each patient, tumor contour was first delineated by a physician. Geometry, intensity, and texture features were extracted from the delineated tumor. Based on UTSW dataset, we completed feature selection and trained a support vector machine (SVM) classifier to predict mutations of BAP1, PBRM1 and VHL genes. We then used TCIA-TCGA dataset to validate the predictive model build upon UTSW dataset. Results: The prediction accuracy of gene expression of TCIA-TCGA patients was 0.83 (20 of 24), 0.83 (20 of 24), and 0.75 (18 of 24) for BAP1, PBRM1, and VHL respectively. For BAP1 gene, texture feature was the most prominent feature type. For PBRM1 gene, intensity feature was the most prominent. For VHL gene, geometry, intensity, and texture features were all important. Conclusion: Using our feature selection strategy and models, we achieved predictive accuracy over 0.75 for all three genes under the condition of using patient data from one institution for training and data from other institutions for testing. These results suggest that radiogenomics can be used to aid in prognosis and used as convenient surrogates for expensive and time consuming gene assay procedures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Using an operator training simulator in the undergraduate chemical engineering curriculim

An operator training simulator (OTS) is to the chemical engineer what a flight simulator is to the aerospace engineer. The basis of an OTS is a high-fidelity dynamic model of a chemical process that allows an engineer to simulate start-up, shut-down, and normal operation. It can also be used to test the skill and ability of an engineer or operator to respond and control some unforeseen situation(s) through the use of programmed malfunctions. West Virginia University (WVU) is a member of the National Energy Technology Laboratory’s Regional University Alliance (NETL-RUA). Working through the NETL-RUA, the authors have spent the last four years collaborating on the development of a high-fidelity OTS for an Integrated Gasification Combined Cycle (IGCC) power plant with CO{sub 2} capture that is the cornerstone of the AVESTARTM (Advanced Virtual Energy Simulation Training And Research) Center with sister facilities at NETL and WVU in Morgantown, WV. This OTS is capable of real-time dynamic simulation of IGCC plant operation, including start-up, shut-down, and power demand load following. The dynamic simulator and its human machine interfaces (HMIs) are based on the DYNSIM and InTouch software, respectively, from Invensys Operations Management. The purpose of this presentation is to discuss the authors’ experiences in using this sophisticated dynamic simulation-based OTS as a hands-on teaching tool in the undergraduate chemical engineering curriculum. At present, the OTS has been used in two separate courses: a new process simulation course and a traditional process control course. In the process simulation course, concepts of steady-state and dynamic simulations were covered prior to exposing the students to the OTS. Moreover, digital logic and the concept of equipment requiring one or more permissive states to be enabled prior to successful operation were also covered. Students were briefed about start-up procedures and the importance of following a predetermined sequence of actions in order to start-up the plant successfully. Student experience with the dynamic simulator consisted of a six-hour training session in which the Claus sulfur capture unit of the IGCC plant was started up. The students were able to operate the simulator through the InTouch-based HMI displays and study and understand the underlying dynamic modeling approach used in the DYNSIM-based simulator. The concepts learned during the training sessions were further reinforced when students developed their own DYNSIM models for a chemical process and wrote a detailed start-up procedure. In the process control course, students learned how the plant responds dynamically to changes in the manipulated inputs, as well as how the control system impacts plant performance, stability, robustness and disturbance rejection characteristics. The OTS provided the opportunity to study the dynamics of complicated, “real-life” process plants consisting of hundreds of pieces of equipment. Students implemented ideal forcing functions, tracked the time-delay through the entire plant, studied the response of open-loop unstable systems, and learned “good practices” in control system design by taking into account the real-world events where significant deviations from the “ideal” or “expected” response can occur. The theory of closed-loop stability was reinforced by implementing limiting proportional gain for stability limits of real plants. Finally, students were divided into several groups where each group was tasked to control a section of the plant within a set of operating limits in the face of disturbances and simulated process faults. At the end of this test, they suggested ways to improve the control system performance based on the theory they learned in class and the hands-on experience they earned while working on the OTS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Scalable multilevel Monte Carlo methods exploiting parallel redistribution on coarse levels

Here, we study an element agglomeration coarsening strategy that requires data redistribution at coarse levels when the number of coarse elements becomes smaller than the number of MPI processes used on the finest level. The overall procedure generates coarse elements (general unstructured unions of fine grid elements) within the framework of element-based algebraic multigrid methods (or AMGe) studied previously. The AMGe-generated coarse spaces have the ability to exhibit approximation properties of the same order as the fine-level spaces since by construction they contain the piecewise polynomials of the same order as on the fine level. These approximation properties are key for the successful use of AMGe in multilevel solvers for nonlinear partial differential equations as well as for multilevel Monte Carlo (MLMC) simulations. The ability to coarsen without being constrained by the number of MPI processes, as described in the present paper, allows to improve the scalability of these solvers as well as the overall MLMC method. The paper illustrates this latter fact with detailed scalability study of MLMC simulations applied to model Darcy equations with a stochastic log-normal permeability field.

AMGe↗
DOE OSTI2019

A scalable mixed-integer decomposition approach for optimal power system restoration

The optimal restoration problem lies at the foundation of the evaluation and improvement of resilience in power systems. In this paper we present a scalable decomposition algorithm, based on the integer L-shaped method, for solving this problem for realistic power systems. The algorithm works by partitioning the problem into a master problem and a slave problem. The master problem optimizes over energization sequences of generators, buses and branches, for the entire restoration horizon, while respecting stylized energization requirements. The slave problem verifies that there exist power flow solutions supporting the energization sequences proposed by the master problem, adding cuts to the master if an infeasible sequence is detected. The power flow verification is carried out using a piece-wise linear approximation of the power flow equations, capable of capturing overvoltage situations common in restoration procedures. We develop specialized binary cuts and hybrid Benders-binary cuts for separating infeasible islands, which are stronger than regular no-good cuts, as well as heuristics for obtaining feasible solutions. We present numerical results for the IEEE 39- and 118-bus systems, and for the Chilean 1500-bus grid, where our algorithm finds near-optimal restoration sequences orders of magnitude faster than the state-of-the-art, demonstrating the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

A Strong Stability Preserving Analysis for Explicit Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions

High-order strong stability preserving (SSP) time discretizations are usually needed to ensure the nonlinear (and sometimes non-inner-product) strong stability properties of spatial discretizations specially designed for the solution of hyperbolic PDEs. Multi-derivative time-stepping methods have recently been increasingly used for evolving hyperbolic PDEs, and the strong stability properties of these methods are of interest. In our prior work we explored time discretizations that preserve the strong stability properties of spatial discretizations coupled with forward Euler and a second-derivative formulation. Yet, many spatial discretizations do not satisfy strong stability properties when coupled with this second-derivative formulation, but rather with a more natural Taylor series formulation. Here, we demonstrate sufficient conditions for an explicit two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and Taylor series formulation. We call these strong stability preserving Taylor series (SSP-TS) methods. We also prove that the maximal order of SSP-TS methods is p=6, and define an optimization procedure that allows us to find such SSP methods. Several types of these methods are presented and their efficiency compared. Finally, these methods are tested on several PDEs to demonstrate the benefit of SSP-TS methods, the need for the SSP property, and the sharpness of the SSP time-step in many cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Respiration Induced Heart Motion and Indications of Gated Delivery for Left-Sided Breast Irradiation

Purpose: To investigate respiration-induced heart motion for left-sided breast irradiation using a four-dimensional computed tomography (4DCT) technique and to determine novel indications to assess heart motion and identify breast patients who may benefit from a gated treatment. Methods and Materials: Images of 4DCT acquired during free breathing for 20 left-sided breast cancer patients, who underwent whole breast irradiation with or without regional nodal irradiation, were analyzed retrospectively. Dose distributions were reconstructed in the phases of 0%, 20%, and 50%. The intrafractional heart displacement was measured in three selected transverse CT slices using D{sub LAD} (the distance from left ascending aorta to a fixed line [connecting middle point of sternum and the body] drawn on each slice) and maximum heart depth (MHD, the distance of the forefront of the heart to the line). Linear regression analysis was used to correlate these indices with mean heart dose and heart dose volume at different breathing phases. Results: Respiration-induced heart displacement resulted in observable variations in dose delivered to the heart. During a normal free-breathing cycle, heart-induced motion D{sub LAD} and MHD changed up to 9 and 11 mm respectively, resulting in up to 38% and 39% increases of mean doses and V{sub 25.2} for the heart. MHD and D{sub LAD} were positively correlated with mean heart dose and heart dose volume. Respiratory-adapted gated treatment may better spare heart and ipsilateral-lung compared with the conventional non-gated plan in a subset of patients with large D{sub LAD} or MHD variations. Conclusion: Proposed indices offer novel assessment of heart displacement based on 4DCT images. MHD and D{sub LAD} can be used independently or jointly as selection criteria for respiratory gating procedure before treatment planning. Patients with great intrafractional MHD variations or tumor(s) close to the diaphragm may particularly benefit from the gated treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Learning Latent Representations to Bridge Coarse-Grained and Atomistic Resolutions in Polymer Simulations

We present a machine-learning-based framework for learning reduced-order representations of polymer chain conformations across coarse-grained (CG) and united-atom (UA) fidelities. By employing linear singular value decomposition and nonlinear autoencoders, we compress high-dimensional polymer configurations into latent spaces with minimal loss of structural accuracy. Crucially, we demonstrate a near-perfect linear mapping between CG and UA latent spaces, enabling an efficient super-resolution back-mapping procedure that reconstructs high-fidelity UA configurations from CG simulations. While minor structural inaccuracies occur, they are effectively corrected through a brief molecular dynamics relaxation, forming a practical hybrid machine learning−physics scheme. This approach establishes the key structural prerequisites for accelerated polymer dynamics simulations: a compact and accurate latent encoding of polymer chain conformations and a validated multi-fidelity mapping that permits reconstruction of UA structures from CG configurations. The extension of this framework to explicit time evolution within the latent space, enabling dynamics to be propagated at CG fidelity and decoded to UA resolution only when required, represents a natural and well-motivated direction for future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗