Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adaptive Sampling”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

AN OPTICAL CATALOG OF GALAXY CLUSTERS OBTAINED FROM AN ADAPTIVE MATCHED FILTER FINDER APPLIED TO SLOAN DIGITAL SKY SURVEY DATA RELEASE 6

We present a new cluster catalog extracted from the Sloan Digital Sky Survey Data Release 6 (SDSS DR6) using an adaptive matched filter (AMF) cluster finder. We identify 69,173 galaxy clusters in the redshift range 0.045 {<=} z < 0.78 in 8420 deg{sup 2} of the sky. We provide angular position, redshift, richness, core, and virial radii estimates for these clusters, as well as an error analysis for each of these quantities. We also provide a catalog of more than 205,000 galaxies representing the three brightest galaxies in the r band which are possible brightest cluster galaxy (BCG) candidates. We show basic properties of the BCG candidates and study how their luminosity scales in redshift and cluster richness. We compare our catalog with the maxBCG and GMBCG catalogs, as well as with that of Wen et al. We match between 30% and 50% of clusters between catalogs over all overlapping redshift ranges. We find that the percentage of matches increases with the richness for all catalogs. We cross match the AMF catalog with available X-ray data in the same area of the sky and find 539 matches, 119 of which with temperature measurements. We present scaling relations between optical and X-ray properties and cluster center comparison. We find that both {Lambda}{sub 200} and R{sub 200} correlate well with both L{sub X} and T{sub X} , with no significant difference in trend if we restrict the matches to flux-limited X-ray samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Adaptive Cybersecurity for Distributed Energy Resources (AdCyDER): Online Reinforcement Learning with Stackelberg-Optimized Defenses — Pipeline Architecture, Evaluation Methodology, and Findings from a Synthetic-Data Evaluation

This report documents the design and evaluation of an integrated online-learning pipeline developed within the AdCyDER project for Distributed Energy Resource (DER) cybersecurity. The pipeline couples a Reinforcement Learning (RL) attack classifier — which produces an attack-type probability distribution — with a Stackelberg game-theoretic (GT) defense selector that consumes those distributions alongside SME-encoded priors over (defense, attack) effectiveness pairings and perdefense costs to choose grid-health-preserving defenses. The objective is not attack classification per se but production of distributions that drive effective defense selection through the Stackelberg layer, learned from delayed grid-health feedback rather than labeled attack data. AdCyDER as a whole is broader than the work presented here; this report covers the specific RL/GT loop integration and its evaluation. We present the integrated pipeline (SCADA telemetry with Fronius inverter physics, Suricata IDS, time-windowed aggregation, per-facility LSTM classifier, Stackelberg optimizer, OpenC2 actuators), an experimental campaign of 28 eight-hour iterations across three baseline modes, and a pipeline-ordered diagnostic protocol. The protocol identifies two distinct failure modes within the loop: paired supervised ceilings on the same features establish that the deployed online RL classifier (macro F1 ≈ 0.07) sits at least 4.7× below a same-architecture supervised LSTM (≈ 0.34) and 10–11× below a linear feature-signal ceiling (≈ 0.70–0.79 depending on per-facility isolation), localizing the dominant failure to the training procedure; and the reward signal driving online updates carries weak directional coupling with classifier correctness in the methodology-expected direction (multi-lens convergent: top-decile P(true) records produce more frequent state changes and slightly larger improvements, top-vs-bot Cohen’s 𝑑 ≈ −0.19), but at effect magnitudes too small to drive gradient-based learning at the campaign sample size. The original learning hypothesis is not supported by the data. The primary contributions are the diagnostic methodology — proposed as a transferable falsification protocol for online RL/GT defense pipelines learning from delayed environmental reward — and the open, reproducible experimental infrastructure. We outline reward reformulation as the highest-priority aspirational next step given the underpowered-but-aligned Q6 reading, with hardware-in-the-loop evaluation as the broadest scope-expansion option.

Blakely, Benjamin [Argonne National Laboratory (AN↗
DOE OSTI2016

Determination of krypton diffusion coefficients in uranium dioxide using atomic scale calculations

We present a study of the diffusion of krypton in UO 2 using atomic scale calculations combined with diffusion models adapted to the system studied. The migration barriers of the elementary mechanisms for interstitial or vacancy assisted migration are calculated in the DFT + U framework using the nudged elastic band method. The attempt frequencies are obtained from the phonon modes of the defect at the initial and saddle points using empirical potential methods. The diffusion coefficients of Kr in UO 2 are then calculated by combining this data with diffusion models accounting for the concentration of vacancies and the interaction of vacancies with Kr atoms. We determined the preferred mechanism for Kr migration and the corresponding diffusion coefficient as a function of the oxygen chemical potential μ O or nonstoichiometry. For very hypostoichiometric (or U-rich) conditions, the most favorable mechanism is interstitial migration. For hypostoichiometric UO 2 , migration is assisted by the bound Schottky defect and the charged uranium vacancy, V U 4– . Around stoichiometry, migration assisted by the charged uranium–oxygen divacancy (V UO 2– ) and V U 4– is the favored mechanism. Finally, for hyperstoichiometric or O-rich conditions, the migration assisted by two V U 4– dominates. Kr migration is enhanced at higher μ O , and in this regime, the activation energy will be between 4.09 and 0.73 eV depending on nonstoichiometry. The experimental values available are in the latter interval. Since it is very probable that these values were obtained for at least slightly hyperstoichiometric samples, our activation energies are consistent with the experimental data, even if further experiments with precisely controlled stoichiometry are needed to confirm these results. Finally, the mechanisms and trends with nonstoichiometry established for Kr are similar to those found in previous studies of Xe.

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

MAPPING THE CLUMPY STRUCTURES WITHIN SUBMILLIMETER GALAXIES USING LASER-GUIDE STAR ADAPTIVE OPTICS SPECTROSCOPY

We present the first integral-field spectroscopic observations of high-redshift submillimeter-selected galaxies (SMGs) using Laser-Guide Star Adaptive Optics. We target H{alpha} emission of three SMGs at redshifts z {approx} 1.4-2.4 with the OH-Suppressing Infrared Imaging Spectrograph on Keck. The spatially resolved spectroscopy of these galaxies reveals unresolved broad-H{alpha} line regions (FWHM >1000 km s{sup -1}) likely associated with an active galactic nucleus (AGN) and regions of diffuse star formation traced by narrow-line H{alpha} emission (FWHM {approx}< 500 km s{sup -1}) dominated by multiple H{alpha}-bright stellar clumps, each contributing 1%-30% of the total clump-integrated H{alpha} emission. We find that these SMGs host high star formation rate surface densities, similar to local extreme sources, such as circumnuclear starbursts and luminous infrared galaxies. However, in contrast to these local environments, SMGs appear to be undergoing such intense activity on significantly larger spatial scales as revealed by extended H{alpha} emission over 4-16 kpc. H{alpha} kinematics show no evidence of ordered global motion as would be found in a disk, but rather large velocity offsets ({approx}few Multiplication-Sign 100 km s{sup -1}) between the distinct stellar clumps. Together with the asymmetric distribution of the stellar clumps around the AGN in these objects, it is unlikely that we are unveiling a clumpy disk structure as has been suggested in other high-redshift populations of star-forming galaxies. The SMG clumps in this sample may correspond to remnants of originally independent gas-rich systems that are in the process of merging, hence triggering the ultraluminous SMG phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Geochemistry and metagenomics analyses of bacterial community structure in selected waste dumpsites in Lagos Metropolis, Nigeria

Dumpsites are reservoirs of persistent organic pollutants (POPs) and heavy metals (HMs), constituting environmental hazards to humanity. Autochthonous microorganisms in dumpsites exhibit various degrees of responses to contaminants. Unfortunately, there is a dearth of information on the types and concentration of pollutants and the array of microorganisms in these dumpsites which may play important roles in the metabolism of such pollutants or other community processes. Therefore, determining the microbial community structure in such contaminated sites across a municipality is essential for profiling the taxa that would serve as consensus degraders of the pollutants. In this study, soil samples from three dumpsites (Cele, CS; Solous, SS; and Computer Village, CVS) were characterized for geochemical properties using GC-MS, MP-AES, and other analytical protocols, while the dynamics of bacterial communities were evaluated based on their 16S rRNA gene barcodes. A significant difference in the bacterial communities was observed among the dumpsites in relation to the extent of pollution caused by POPs and HMs. CVS, with the highest HM contamination, was rich in Actinobacteria (41.7%) and Acidobacteria (10.2%), in contrast to CS and SS. Proteobacteria (34.1%) and Firmicutes (20%) were the dominant phyla in CS (highest POP contamination), while Bacteroidetes (45.5%) and Proteobacteria (39.9%) were dominant in SS soil.Bacilluswas the dominant genus in the most polluted dumpsite. Canonical correspondence analysis revealed that polycyclic aromatic hydrocarbons (PAHs) and HMs shaped the structure of the bacterial operational taxonomic units (OTUs) in the most polluted dumpsite. Out of a total of 706 OTUs, 628 OTUs exhibited a significant correlation (>50%) with benzo(b)fluoranthene, azobenzene, dibenzofurans, pyrene, dibenzo(a,l)pyrene, Cu, and Zn. In particular, Proteobacteria (Achromobactersp. andSerratiasp.), Bacteroidetes (Zhouiasp.), and Firmicutes (Bacillussp.) were suggested to be pivotal to the ecophysiology of dumpsite soils contaminated with POPs and HMs. The results generally underscored the importance of metagenomic and physicochemical analyses of polluted systems in enabling correlations for useful prediction of drivers of such ecosystems. This will further improve our understanding of the metabolic potential and adaptation of organisms in such systems.

Environmental Sciences & Ecology↗
DOE OSTI2012

THE FORMATION OF THE FIRST COSMIC STRUCTURES AND THE PHYSICS OF THE z {approx} 20 UNIVERSE

We perform a suite of cosmological simulations in the {Lambda}CDM paradigm of the formation of the first structures in the universe prior to astrophysical reheating and reionization (15 {approx}< z < 200). These are the first simulations initialized in a manner that self-consistently accounts for the impact of pressure on the rate of growth of modes, temperature fluctuations in the gas, and the dark matter-baryon supersonic velocity difference. Even with these improvements, these are still difficult times to simulate accurately as the Jeans length of the cold intergalactic gas must be resolved while also capturing a representative sample of the universe. We explore the box size and resolution requirements to meet these competing objectives. Our simulations support the finding of recent studies that the dark matter-baryon velocity difference has a surprisingly large impact on the accretion of gas onto the first star-forming minihalos (which have masses of {approx}10{sup 6} M {sub Sun }). In fact, the halo gas is often significantly downwind of such halos and with lower densities in the simulations in which the baryons have a bulk flow with respect to the dark matter, modulating the formation of the first stars by the local value of this velocity difference. We also show that dynamical friction plays an important role in the nonlinear evolution of the dark matter-baryon differential velocity, acting to erase this velocity difference quickly in overdense gas, as well as sourcing visually apparent bow shocks and Mach cones throughout the universe. We use simulations with both the GADGET and Enzo cosmological codes to test the robustness of these conclusions. The comparison of these codes' simulations also provides a relatively controlled test of these codes themselves, allowing us to quantify some of the tradeoffs between the algorithms. For example, we find that particle coupling in GADGET between the gas and dark matter particles results in spurious growth that mimics nonlinear growth in the matter power spectrum for standard initial setups. This coupling is alleviated by using adaptive gravitational softening for the gas. In a companion paper, we use the simulations presented here to make detailed estimates for the impact of the dark matter-baryon velocity differential on redshifted 21 cm radiation. The initial conditions generator used in this study, CICSASS, can be publicly downloaded.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Alternative High-Performance Ceramic Waste Forms

This final report (M5NU-12-NY-AU # 0202-0410) summarizes the results of the project titled “Alternative High-Performance Ceramic Waste Forms,” funded in FY12 by the Nuclear Energy University Program (NEUP Project # 12-3809) being led by Alfred University in collaboration with Savannah River National Laboratory (SRNL). The overall focus of the project is to advance fundamental understanding of crystalline ceramic waste forms and to demonstrate their viability as alternative waste forms to borosilicate glasses. We processed single- and multiphase hollandite waste forms based on simulated waste streams compositions provided by SRNL based on the advanced fuel cycle initiative (AFCI) aqueous separation process developed in the Fuel Cycle Research and Development (FCR&D). For multiphase simulated waste forms, oxide and carbonate precursors were mixed together via ball milling with deionized water using zirconia media in a polyethylene jar for 2 h. The slurry was dried overnight and then separated from the media. The blended powders were then subjected to melting or spark plasma sintering (SPS) processes. Microstructural evolution and phase assemblages of these samples were studied using x-ray diffraction (XRD), scanning electron microscopy (SEM), energy dispersion analysis of x-rays (EDAX), wavelength dispersive spectrometry (WDS), transmission electron spectroscopy (TEM), selective area x-ray diffraction (SAXD), and electron backscatter diffraction (EBSD). These results showed that the processing methods have significant effect on the microstructure and thus the performance of these waste forms. The Ce substitution into zirconolite and pyrochlore materials was investigated using a combination of experimental (in situ XRD and x-ray absorption near edge structure (XANES)) and modeling techniques to study these single phases independently. In zirconolite materials, a transition from the 2M to the 4M polymorph was observed with increasing Ce content. The resulting powders were consolidated via SPS. Ce was reduced to the trivalent oxidation state and the zirconolite was converted into undesirable perovskite. The zirconolite polymorphs found in the synthesized powders were recovered after a post-SPS heat treatment in air. These results demonstrated the potential of processing in controlling the phase assemblage in these waste forms. Hollandites with Cr 3+ trivalent cations were identified as potential hosts for Cs immobilization and are being investigated for Cs retention properties. Series of compositions Ba 1.15-x Cs 2x Cr 2.3 Ti 5.7 O 16 , with increasing Cs loadings, were prepared by sol-gel process and characterized for structural parameters. Structural characterization was performed by a combination of powder XRD and neutron powder diffraction. Phase pure hollandite adapting monoclinic symmetry (I2/m) was observed for 0 ≤ x ≤ 0.55. These results were used to develop a new structural model to interpret Cs immobilization in these hollandites. Performance of these waste forms were evaluated for chemical durability and radiation resistance. Product consistency testing (PCT) and vapor hydration testing (VHT) were used for testing of chemical durability. Radiation resistance was tested using He + ions to simulate $\alpha$ particles and heavy ions such as Au 3+ to simulate a recoil. These results showed that these waste forms were chemically durable. The waste forms also amorphized to various degrees on exposure to simulated radiation.

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

MOA 2010-BLG-477Lb: CONSTRAINING THE MASS OF A MICROLENSING PLANET FROM MICROLENSING PARALLAX, ORBITAL MOTION, AND DETECTION OF BLENDED LIGHT

Microlensing detections of cool planets are important for the construction of an unbiased sample to estimate the frequency of planets beyond the snow line, which is where giant planets are thought to form according to the core accretion theory of planet formation. In this paper, we report the discovery of a giant planet detected from the analysis of the light curve of a high-magnification microlensing event MOA 2010-BLG-477. The measured planet-star mass ratio is q = (2.181 {+-} 0.004) Multiplication-Sign 10{sup -3} and the projected separation is s = 1.1228 {+-} 0.0006 in units of the Einstein radius. The angular Einstein radius is unusually large {theta}{sub E} = 1.38 {+-} 0.11 mas. Combining this measurement with constraints on the 'microlens parallax' and the lens flux, we can only limit the host mass to the range 0.13 < M/M{sub Sun} < 1.0. In this particular case, the strong degeneracy between microlensing parallax and planet orbital motion prevents us from measuring more accurate host and planet masses. However, we find that adding Bayesian priors from two effects (Galactic model and Keplerian orbit) each independently favors the upper end of this mass range, yielding star and planet masses of M{sub *} = 0.67{sup +0.33}{sub -0.13} M{sub Sun} and m{sub p} = 1.5{sup +0.8}{sub -0.3} M{sub JUP} at a distance of D = 2.3 {+-} 0.6 kpc, and with a semi-major axis of a = 2{sup +3}{sub -1} AU. Finally, we show that the lens mass can be determined from future high-resolution near-IR adaptive optics observations independently from two effects, photometric and astrometric.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-672: A Novel Kernel-Based Dose Engine for KeV Photon Beams

Purpose: Mimicking state-of-the-art patient radiotherapy with high precision irradiators for small animals allows advanced dose-effect studies and radiobiological investigations. One example is the implementation of pre-clinical IMRT-like irradiations, which requires the development of inverse planning for keV photon beams. As a first step, we present a novel kernel-based dose calculation engine for keV x-rays with explicit consideration of energy and material dependencies. Methods: We follow a superposition-convolution approach adapted to keV x-rays, based on previously published work on micro-beam therapy. In small animal radiotherapy, we assume local energy deposition at the photon interaction point, since the electron ranges in tissue are of the same order of magnitude as the voxel size. This allows us to use photon-only kernel sets generated by MC simulations, which are pre-calculated for six energy windows and ten base materials. We validate our stand-alone dose engine against Geant4 MC simulations for various beam configurations in water, slab phantoms with bone and lung inserts, and on a mouse CT with (0.275mm)3 voxels. Results: We observe good agreement for all cases. For field sizes of 1mm{sup 2} to 1cm{sup 2} in water, the depth dose curves agree within 1% (mean), with the largest deviations in the first voxel (4%) and at depths>5cm (<2.5%). The out-of-field doses at 1cm depth agree within 8% (mean) for all but the smallest field size. In slab geometries, the mean agreement was within 3%, with maximum deviations of 8% at water-bone interfaces. The γ-index (1mm/1%) passing rate for a single-field mouse irradiation is 71%. Conclusion: The presented dose engine yields an accurate representation of keV-photon doses suitable for inverse treatment planning for IMRT. It has the potential to become a significantly faster yet sufficiently accurate alternative to full MC simulations. Further investigations will focus on energy sampling as well as calculation times. Research at ICR is also supported by Cancer Research UK under Programme C33589/A19727 and NHS funding to the NIHR Biomedical Research Centre at RMH and ICR. MFF is supported by Cancer Research UK under Programme C33589/A19908.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Adaptation of the CVT algorithm for catheter optimization in high dose rate brachytherapy

Purpose: An innovative, simple, and fast method to optimize the number and position of catheters is presented for prostate and breast high dose rate (HDR) brachytherapy, both for arbitrary templates or template-free implants (such as robotic templates).Methods: Eight clinical cases were chosen randomly from a bank of patients, previously treated in our clinic to test our method. The 2D Centroidal Voronoi Tessellations (CVT) algorithm was adapted to distribute catheters uniformly in space, within the maximum external contour of the planning target volume. The catheters optimization procedure includes the inverse planning simulated annealing algorithm (IPSA). Complete treatment plans can then be generated from the algorithm for different number of catheters. The best plan is chosen from different dosimetry criteria and will automatically provide the number of catheters and their positions. After the CVT algorithm parameters were optimized for speed and dosimetric results, it was validated against prostate clinical cases, using clinically relevant dose parameters. The robustness to implantation error was also evaluated. Finally, the efficiency of the method was tested in breast interstitial HDR brachytherapy cases.Results: The effect of the number and locations of the catheters on prostate cancer patients was studied. Treatment plans with a better or equivalent dose distributions could be obtained with fewer catheters. A better or equal prostate V100 was obtained down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of prostate V100 and D90. Implantation errors up to 3 mm were acceptable since no statistical difference was found when compared to 0 mm error (p > 0.05). No significant difference in dosimetric indices was observed for the different combination of parameters within the CVT algorithm. A linear relation was found between the number of random points and the optimization time of the CVT algorithm. Because the computation time decrease with the number of points and that no effects were observed on the dosimetric indices when varying the number of sampling points and the number of iterations, they were respectively fixed to 2500 and to 100. The computation time to obtain ten complete treatments plans ranging from 9 to 18 catheters, with the corresponding dosimetric indices, was 90 s. However, 93% of the computation time is used by a research version of IPSA. For the breast, on average, the Radiation Therapy Oncology Group recommendations would be satisfied down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of V100, dose homogeneity index, and D90.Conclusions: The authors have devised a simple, fast and efficient method to optimize the number and position of catheters in interstitial HDR brachytherapy. The method was shown to be robust for both prostate and breast HDR brachytherapy. More importantly, the computation time of the algorithm is acceptable for clinical use. Ultimately, this catheter optimization algorithm could be coupled with a 3D ultrasound system to allow real-time guidance and planning in HDR brachytherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Genetic monitoring of steelhead in the Klickitat River to estimate productivity, straying, and migration timing

Abstract Objective Salmonids with a complex life history variation present challenges for conservation management, but genetic approaches alongside fisheries monitoring can address questions regarding the viability of the natural populations. Methods We genotyped adult (n = 3108) and juvenile (n = 2624) samples of steelhead Oncorhynchus mykiss that were collected in the Klickitat River, Washington, USA at traps in the lower drainage to examine tributary level productivity, straying from outside sources, and variation in adult migration timing. Result Genetic assignment of steelhead from this system indicated that the majority were produced within or near tributaries of the middle Klickitat River (juvenile mean = 72.8%; adult mean = 87.3%). Analyses with parentage-based tagging identified that most hatchery-origin adults were assigned to the Skamania Hatchery (80.8%), as expected, since this has been the release stock for decades within the Klickitat River drainage. Hatchery-origin adults were also identified from programs operating outside the Klickitat River, which were primarily strays from Snake River hatcheries. Most natural-origin steelhead were assigned to the Klickitat River, but there were also natural-origin fish identified as strays from other regions of the Columbia River (22.3% of natural returns). We also examined genes known to be associated with migration timing in adult steelhead observed at the trap and observed a strong relationship between migration date and alleles for early and late migration, but individual outliers were detected across seasons. Conclusion Our results indicate that genetic variation of steelhead in the Klickitat River has been influenced by hatchery programs as well as natural-origin straying from other subbasins, but genetic diversity remains high throughout the subbasin, and both early and late migration alleles are maintained. The genetic diversity present in Klickitat River steelhead may enable this Endangered Species Act listed (threatened) species to better adapt to stochastic environmental conditions compared to less diverse populations.

Collins, Erin E. (ORCID:0000000260978479)↗
DOE OSTI2014

Validation of deformable image registration algorithms on CT images of ex vivo porcine bladders with fiducial markers

Purpose: The spatial accuracy of deformable image registration (DIR) is important in the implementation of image guided adaptive radiotherapy techniques for cancer in the pelvic region. Validation of algorithms is best performed on phantoms with fiducial markers undergoing controlled large deformations. Excised porcine bladders, exhibiting similar filling and voiding behavior as human bladders, provide such an environment. The aim of this study was to determine the spatial accuracy of different DIR algorithms on CT images ofex vivo porcine bladders with radiopaque fiducial markers applied to the outer surface, for a range of bladder volumes, using various accuracy metrics. Methods: Five excised porcine bladders with a grid of 30–40 radiopaque fiducial markers attached to the outer wall were suspended inside a water-filled phantom. The bladder was filled with a controlled amount of water with added contrast medium for a range of filling volumes (100–400 ml in steps of 50 ml) using a luer lock syringe, and CT scans were acquired at each filling volume. DIR was performed for each data set, with the 100 ml bladder as the reference image. Six intensity-based algorithms (optical flow or demons-based) implemented in theMATLAB platform DIRART, a b-spline algorithm implemented in the commercial software package VelocityAI, and a structure-based algorithm (Symmetric Thin Plate Spline Robust Point Matching) were validated, using adequate parameter settings according to values previously published. The resulting deformation vector field from each registration was applied to the contoured bladder structures and to the marker coordinates for spatial error calculation. The quality of the algorithms was assessed by comparing the different error metrics across the different algorithms, and by comparing the effect of deformation magnitude (bladder volume difference) per algorithm, using the Independent Samples Kruskal-Wallis test. Results: The authors found good structure accuracy without dependency on bladder volume difference for all but one algorithm, and with the best result for the structure-based algorithm. Spatial accuracy as assessed from marker errors was disappointing for all algorithms, especially for large volume differences, implying that the deformations described by the registration did not represent anatomically correct deformations. The structure-based algorithm performed the best in terms of marker error for the large volume difference (100–400 ml). In general, for the small volume difference (100–150 ml) the algorithms performed relatively similarly. The structure-based algorithm exhibited the best balance in performance between small and large volume differences, and among the intensity-based algorithms, the algorithm implemented in VelocityAI exhibited the best balance. Conclusions: Validation of multiple DIR algorithms on a novel physiological bladder phantom revealed that the structure accuracy was good for most algorithms, but that the spatial accuracy as assessed from markers was low for all algorithms, especially for large deformations. Hence, many of the available algorithms exhibit sufficient accuracy for contour propagation purposes, but possibly not for accurate dose accumulation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-BRA-01: 4D Monte Carlo Simulations for Verification of Dose Delivered to a Moving Anatomy

Purpose: To validate 4D Monte Carlo (MC) simulations of dose delivery by an Elekta Agility linear accelerator to a moving phantom. Methods: Monte Carlo simulations were performed using the 4DdefDOSXYZnrc/EGSnrc user code which samples a new geometry for each incident particle and calculates the dose in a continuously moving anatomy. A Quasar respiratory motion phantom with a lung insert containing a 3 cm diameter tumor was used for dose measurements on an Elekta Agility linac with the phantom in stationary and moving states. Dose to the center of tumor was measured using calibrated EBT3 film and the RADPOS 4D dosimetry system. A VMAT plan covering the tumor was created on the static CT scan of the phantom using Monaco V.5.10.02. A validated BEAMnrc model of our Elekta Agility linac was used for Monte Carlo simulations on stationary and moving anatomies. To compare the planned and delivered doses, linac log files recorded during measurements were used for the simulations. For 4D simulations, deformation vectors that modeled the rigid translation of the lung insert were generated as input to the 4DdefDOSXYZnrc code as well as the phantom motion trace recorded with RADPOS during the measurements. Results: Monte Carlo simulations and film measurements were found to agree within 2mm/2% for 97.7% of points in the film in the static phantom and 95.5% in the moving phantom. Dose values based on film and RADPOS measurements are within 2% of each other and within 2σ of experimental uncertainties with respect to simulations. Conclusion: Our 4D Monte Carlo simulation using the defDOSXYZnrc code accurately calculates dose delivered to a moving anatomy. Future work will focus on more investigation of VMAT delivery on a moving phantom to improve the agreement between simulation and measurements, as well as establishing the accuracy of our method in a deforming anatomy. This work was supported by the Ontario Consortium of Adaptive Interventions in Radiation Oncology (OCAIRO), funded by the Ontario Research Fund Research Excellence program.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Performance of 12 DIR algorithms in low-contrast regions for mass and density conserving deformation

Purpose: Deformable image registration (DIR) has become a key tool for adaptive radiotherapy to account for inter- and intrafraction organ deformation. Of contemporary interest, the application to deformable dose accumulation requires accurate deformation even in low contrast regions where dose gradients may exist within near-uniform tissues. One expects high-contrast features to generally be deformed more accurately by DIR algorithms. The authors systematically assess the accuracy of 12 DIR algorithms and quantitatively examine, in particular, low-contrast regions, where accuracy has not previously been established.Methods: This work investigates DIR algorithms in three dimensions using deformable gel (DEFGEL) [U. J. Yeo, M. L. Taylor, L. Dunn, R. L. Smith, T. Kron, and R. D. Franich, “A novel methodology for 3D deformable dosimetry,” Med. Phys. 39, 2203–2213 (2012)], for application to mass- and density-conserving deformations. CT images of DEFGEL phantoms with 16 fiducial markers (FMs) implanted were acquired in deformed and undeformed states for three different representative deformation geometries. Nonrigid image registration was performed using 12 common algorithms in the public domain. The optimum parameter setup was identified for each algorithm and each was tested for deformation accuracy in three scenarios: (I) original images of the DEFGEL with 16 FMs; (II) images with eight of the FMs mathematically erased; and (III) images with all FMs mathematically erased. The deformation vector fields obtained for scenarios II and III were then applied to the original images containing all 16 FMs. The locations of the FMs estimated by the algorithms were compared to actual locations determined by CT imaging. The accuracy of the algorithms was assessed by evaluation of three-dimensional vectors between true marker locations and predicted marker locations.Results: The mean magnitude of 16 error vectors per sample ranged from 0.3 to 3.7, 1.0 to 6.3, and 1.3 to 7.5 mm across algorithms for scenarios I to III, respectively. The greatest accuracy was exhibited by the original Horn and Schunck optical flow algorithm. In this case, for scenario III (erased FMs not contributing to driving the DIR calculation), the mean error was half that of the modified demons algorithm (which exhibited the greatest error), across all deformations. Some algorithms failed to reproduce the geometry at all, while others accurately deformed high contrast features but not low-contrast regions—indicating poor interpolation between landmarks.Conclusions: The accuracy of DIR algorithms was quantitatively evaluated using a tissue equivalent, mass, and density conserving DEFGEL phantom. For the model studied, optical flow algorithms performed better than demons algorithms, with the original Horn and Schunck performing best. The degree of error is influenced more by the magnitude of displacement than the geometric complexity of the deformation. As might be expected, deformation is estimated less accurately for low-contrast regions than for high-contrast features, and the method presented here allows quantitative analysis of the differences. The evaluation of registration accuracy through observation of the same high contrast features that drive the DIR calculation is shown to be circular and hence misleading.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Miniaturized Magnetoelastic Sensor System

This article describes the design, assembly, and implementation of a hand-held, magnetic-field-based sensor system that can be adapted for a variety of sensing applications. The miniaturized system is based on Chemical Identification by Magneto-Elastic Sensing (ChIMES) technology, which uses three concentric solenoid coils to wirelessly interrogate a sensor body comprised of a response material coupled to a magnetoelastic wire. The response material expands when it encounters a target, imposing mechanical stress on the wire and altering its magnetic permeability. The sensor bodies are passive, requiring no external power source, and they are small, measuring about 15 mm in length and 3.0 mm in diameter. Up to four sensor bodies can be configured as an evenly-spaced linear array. The sensor system operates by applying a low-frequency, current-stabilized, filtered triangle wave to a uniform-density excitation coil to switch the magnetic domains within the wire. Further, the responses from the sensors are picked up by a detection coil as stress-induced changes in the Faraday voltage, and the strong magnetic field induced by the excitation coil in the detection coil is nullified by a cancellation coil reverse-wound in series with the detection coil. The responses of the sensors in an array are separated in time by a linear gradient dc biasing coil. The sensors can be interrogated through metallic and nonmetallic barriers. The signals from the detection coil and the excitation coil are digitized by a pair of bipolar analog-to-digital converters (ADCs). A Raspberry Pi single-board computer (SBC) and associated software perform data acquisition and control all aspects of the sensor system hardware. The program allows the user to select the number of sensors in the array, the type of signal that is being collected, and the number of samples to take. The program also allows for signal processing of the sensor data, such as baseline correction. The program can differentiate sensor peaks from each other and calculate the magnitude of each sensor response with less than 1% error. The data are then displayed along with a graph of the signal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Mechanisms of interfacial reactivity in near surface and extreme environments

The local water structure surrounding ions in aqueous solutions greatly affects their chemical properties such as reaction rates, ion association, and proton and electron transport. These properties result in the behavior of ions in natural aqueous environments. For example ore transport is facilitated by chloride ion pair formation and the reaction of ions in an interface is strongly dependent on the dehydration of the ion hydration shell. We are developing the use of high-­resolution XAFS observations and 1st principles based MD-­XAFS analysis (spectra simulated using 1st principle methods with no adjustable parameters, AIMD) to interpret the solution properties of strongly interacting aqueous solutes under arbitrary pressure and temperature conditions. In the 1st principle MD-­XAFS method density functional theory (DFT) based MD simulations(Car and Parrinello, 1985) are used to generate a large ensemble of structural snap shots of the hydration region. These are then used to generate scattering intensities. I emphasize three points about this novel approach to analyzing XAFS data. 1st: As illustrated in Figure 1, the level of agreement between the calculated and observed intensities is considerably higher than has been obtained by any XAFS analysis to date (note 2nd shell region, R> 2 Å). 2nd: This result was obtained from a parameter free simulation with no fitting of the interaction potentials to any data. This supports the use of these methods for more difficult environments and more complex solutes (polyions). 3rd: New information about the shell structure (Figure 1) is now available because of this more detailed agreement. We note also that both multiple scattering and second shell features are well represented in the analysis. As far as we know this is the 1st analysis of second shell structure and multiple scattering. Excellent agreement has been obtained for most of the third row metal ions: Ca 2+ , Zn 2+ , Cu 2+ , Ni 2+ , Co 2+ , Mn 2+ , Fe 3+ , Cr 3+ . Calculations on these systems are demanding because of their open electronic shells, and high ionic charge. Principal Investigator: Professor John Weare (University of California, San Diego) The prediction of the interactions of geochemical fluids with minerals, nanoparticles, and colloids under extreme near surface conditions of temperature (T) and pressure (P) is a grand challenge research need in geosciences (U.S. DOE 2007, Basic Research Needs for Geosciences: Facilitating the 21st Energy Systems.). To evaluate the impact of these processes on energy production and management strategies it is necessary to have a high level of understanding of the interaction between complex natural fluids and mineral formations. This program emphasizes 1st principle parameter free simulations of complex chemical processes in solutions, in the mineral phase, and in the interfaces between these phases The development of new computational tools (with emphasis on oxide materials and reaction dynamics) tailored to treat wide range of conditions and time scales experienced in such geochemical applications is has been developed. Because of the sensitivity of the interaction in these systems to electronic structure and local bonding environments, and of the need to describe bond breaking/formation, our simulations are based on interactions calculated at the electronic structure level (ab-initio molecular dynamics, AIMD). The progress in the computational aspects of program may be summarized in terms of the following themes (objectives); Development of efficient parameter free dynamical simulation technology based on 1st principles force and energy calculations especially adapted for geochemical applications (e.g., mineral, interfaces and aqueous solutions) (continuing program); Calculation of the dynamics of water structure of in the surface-water interface of transition metal oxides and oxihydroxides; and Development of improved (beyond DFT+GGA) electronic structure calculations for minerals and the interface region that more accurately calculate electron correlation, spin density, and localization. The focus of the program is also on the iron oxide and oxihydroxide minerals and Fe 2+ (aq)/Fe 3+ (cr) oxidation in the mineral solution interface region. These methods included the development of model Hamilitonian methods that can be solved to near convergence for single site models (DMFT) and many-body perturbation methods (MP2, GW); Development of time decomposition methods to extend time scales of molecular dynamics (MD) simulations and support the use of high complexity electronic structure calculations (MP2, CCSD(T)) of forces for use in dynamical simulations where very high chemical accuracy is required (microsolvated reactions in absorbed surface layers); and The development of a new linear scaling finite element solver for eigenvalue problem that supports solution of quantum problems with unusual potential and boundary values. Application progress of the above new simulation technology to problems of geochemical interests includes; The prediction of metal oxide surface structure and the reduction/oxidation of Fe 3+ (cr)/Fe 2+ (aq) in metal oxide (hematite, goethite)/solution interfaces. Result: water interacts strongly with the 001 Hematite surface; interaction of water with the 100 goethite is weak; The study of ion solvation and the composition of ion hydration shells under extreme conditions (focus on Fe 3+/2+ , Al 3+ and Mg 2+ and their hydroxide speciation). Result: Ion association in water solutions can be calculated from 1st principle methods. Efficient sampling of the free energy requires more development; The continued development of new high resolution analysis of XAFS scattering of disordered systems (particularly Al, Mg) and of XANES calculations for aqueous ions. Result: EXAFS spectra can be calculated to high accuracy with DFT level dynamic simulations; The exploration of electron localization and electron transport in metal oxides (highly correlated materials). Result: Proper description of electron localization requires levels of calculation beyond DFT; and Localization of electrons in DFT type Hamiltonians was studied. Result: For very Dirac high exchange new solutions (New unphysical bifrucations) to the eigenvalue problem are found. The program was highly collaborative involving faculty and students in mathematics, physics and computer science departments as well as coworkers at the Pacific Northwest National Laboratories (PNNL). The students in this program had the opportunity to develop skills in the development of methods, the implementation of method on high performance parallel computers and the application of these methods to problem in geochemical science. Much of the software that was developed was incorporated in the NWchem software package maintained by PNNL.

58 GEOSCIENCES↗
DOE OSTI2025

Degradable Biocomposite Thermoplastic Polyurethanes

In this project, the team developed tough and degradable biocomposite thermoplastic polyurethanes (TPUs) by incorporating bacterial spores into TPUs as a biofunctional living filler. The team screened various bacteria and selected the Bacillus subtilis ATCC 6633 strain as the final candidate, primarily due to its genomic availability, sporulation ability and TPU assimilation activity. The heat-shock tolerance of ATCC 6633 spores was further improved through evolutionary engineering via Adaptive Laboratory Evolution (ALE), demonstrating a 17.7-fold enhanced germination efficiency post heat-shock treatment compared to the wild-type strain (WT). The team fabricated biocomposite TPUs by incorporating lyophilized powder of heat-shock tolerized (HST) spores during the hot melt extrusion (HME) of TPU at 135 °C. The baseline TPU used in this project is a commercially available soft-grade TPU (BCF45) manufactured by BASF. Colony forming unit (CFU) assays quantified that WT and HST spores in the TPU matrix retained approximately 20% and 100% survivability, respectively, after HME. Tensile testing demonstrated that the spores behaved as a polymer-reinforcing filler, positively affecting the overall tensile properties of the biocomposite TPU. For example, biocomposite TPU with WT and HST spores (BC TPU WT and BC TPU HST , respectively) exhibited up to 25% and 37% improved toughness, respectively, compared to TPU without spores. BC TPU HST showed remarkably improved disintegration in autoclaved compost (92% mass loss in 5 months), which simulated a microbially poor environment for TPU degradation. When compared to TPU without spores (44% mass loss in 5 months) the acceleration of degradation is marked. Respirometry confirmed that 72% of BC TPU HST was biomineralized into CO2 within 6 months, indicating that spores in the biocomposite TPU were germinated by utilizing nutrients in the autoclaved compost, facilitating TPU degradation at the end of the material's life. The team demonstrated the scale-up of biocomposite TPU fabrication using continuous extrusion and injection molding techniques. Processing conditions optimized in a lab-scale microcompounder were successfully transferred to a continuous extruder with a 30-fold increased throughput. Biocomposite TPUs prepared using these industry-relevant processes showed comparable toughness improvements to samples prepared in the lab-scale extruder. Excitingly, following compounding in the pilot-extruder the composite material could be injection molded, while retaining high spore viability and similar toughness improvements. The team also found that spores in biocomposite TPU served as antioxidants, preventing toughness decay during the recycled extrusion of BC TPU HST . Long-term storage tests over one year showed that the addition of spores had no negative effect on the longevity of the TPU. Furthermore, the team demonstrated the fabrication of spore-bearing biocomposite polymers with other polyesters such as PBAT, PLA, and PCL. We obtained promising preliminary data that showed overall toughness improvements for all polymers with spore addition. Finally, life cycle assessment (LCA) and techno-economic analysis (TEA) were carried out, which indicated minimal additional cost of fabrication. Overall, a tough and degradable biocomposite thermoplastic was successfully developed through this project, with all tasks completed successfully, achieving >100% of the objectives.

36 MATERIALS SCIENCE↗