Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Dynamic Programming”

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.

88 records · Page 5

An EvoDevo Study of Salmonid Visual Opsin Dynamics and Photopigment Spectral Sensitivity

Salmonids are ideal models as many species follow a distinct developmental program from demersal eggs and a large yolk sac to hatching at an advanced developmental stage. Further, these economically important teleosts inhabit both marine- and freshwaters and experience diverse light environments during their life histories. At a genome level, salmonids have undergone a salmonid-specific fourth whole genome duplication event (Ss4R) compared to other teleosts that are already more genetically diverse compared to many non-teleost vertebrates. Thus, salmonids display phenotypically plastic visual systems that appear to be closely related to their anadromous migration patterns. This is most likely due to a complex interplay between their larger, more gene-rich genomes and broad spectrally enriched habitats; however, the molecular basis and functional consequences for such diversity is not fully understood. This study used advances in genome sequencing to identify the repertoire and genome organization of visual opsin genes (those primarily expressed in retinal photoreceptors) from six different salmonids [Atlantic salmon ( Salmo salar ), brown trout ( Salmo trutta ), Chinook salmon ( Oncorhynchus tshawytcha ), coho salmon ( Oncorhynchus kisutch ), rainbow trout ( Oncorhynchus mykiss ), and sockeye salmon ( Oncorhynchus nerka )] compared to the northern pike ( Esox lucius ), a closely related non-salmonid species. Results identified multiple orthologues for all five visual opsin classes, except for presence of a single short-wavelength-sensitive-2 opsin gene. Several visual opsin genes were not retained after the Ss4R duplication event, which is consistent with the concept of salmonid rediploidization. Developmentally, transcriptomic analyzes of Atlantic salmon revealed differential expression within each opsin class, with two of the long-wavelength-sensitive opsins not being expressed before first feeding. Also, early opsin expression in the retina was located centrally, expanding dorsally and ventrally as eye development progressed, with rod opsin being the dominant visual opsin post-hatching. Modeling by spectral tuning analysis and atomistic molecular simulation, predicted the greatest variation in the spectral peak of absorbance to be within the Rh2 class, with a ∼40 nm difference in λ max values between the four medium-wavelength-sensitive photopigments. Overall, it appears that opsin duplication and expression, and their respective spectral tuning profiles, evolved to maximize specialist color vision throughout an anadromous lifecycle, with some visual opsin genes being lost to tailor marine-based vision.

59 BASIC BIOLOGICAL SCIENCES↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

Final Report on CFD and Thermal-Mechanical Stress Analysis of PWR Surge Line under Transient Condition Thermal Stratification and an Evolutionary Cyclic Plasticity Based Transformative Fatigue Evaluation Approach without Using S~N Curve: Rev. 1

This report presents an evolutionary cyclic plasticity based transformative approach for fatigue evaluation of safety critical components. This approach is generic which can be used not only for the present generation light water reactor components (which is our focus) but also for advanced reactor components. The approach is based on the fundamental concept of time-evolution of progressive fatigue damage rather than the end-of-life data based conventional S~N curve approaches. Series of fatigue tests on 316 stainless steel uniaxial specimens were conducted under strain-controlled and stress-controlled cyclic loading in two different environment such as in-air at 300 °C and PWR coolant water condition at 300 °C. Material models are developed to capture the uniaxial test data and are subsequently programmed into commercially available ABAQUS code to transform the uniaxial material behavior to multiaxial loading domain. The developed evolutionary cyclic-plasticity approach is verified by both analytical and 3D-FE modeling of uniaxial fatigue test specimens. Results show that the developed material models can capture the time-dependence of material hardening and softening with great accuracy. In addition, results show that the material models not only can capture the material behavior under constant amplitude loading but also can capture the load-sequence effect under variable and random amplitude loading. The ANL developed approach is one of its kind, which shows first in the fatigue research community that fatigue evaluation of a component (in our case verified for laboratory scale fatigue test specimens with great accuracy) can be performed for its entire fatigue life, including the hardening and softening behavior, without using a conventional S~N curve based approach. Additionally results are presented related to fatigue evaluation of real reactor components such as of a pressurized water reactor (PWR) surge line (SL), under upset transient conditions. SLs of PWR experience complex thermal stratification due to the mixing of differential temperature flow from hot leg and pressurizer. A detailed parametric study of thermal stratification in SL is performed through both standalone computational fluid dynamics (CFD) analysis and fluid-structure-interaction based coupled CFD and heat transfer analysis. Results show that the mass flow rate in SL substantially affects the extent of its thermal stratification. The resulting thermal boundaries conditioned are used as inputs for thermal-mechanical stress analysis and component level fatigue evaluation. Fatigue evaluation of SL pipe under example loading conditions was performed using the above-mentioned evolutionary-plasticity based approach and ASME code (for in-air life) and NUREG-6909 (for PWR environment life) based approaches. Results show that the evolutionary plasticity based approach estimate substantially different life compared to the conventional ASME based approach. To show the entirety of the proposed evolutionary fatigue-modeling framework, this report includes the results presented in semi-annual March report (ANL/LWRS-17/01). Sections 2-5 of March report are directly included as Sections 2-5 in this report, without any significant changes. However, Sections 6-9 of this report present the work performed during March-September of 2017. Whereas, Sections 1 and 10 present the organization of this report and the summary and future work, respectively. This report (ANL/LWRS-17/03 Rev. 1) is a slightly revised version of the earlier published report (ANL/LWRS-17/03). In this revised version, some of the results and associated discussion have been corrected/updated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY 2016 Status Report: Documentation of All CIRFT Data including Hydride Reorientation Tests (Draft M2)

The first portion of this report provides a detailed description of fiscal year (FY) 2015 test result corrections and analysis updates based on FY 2016 updates to the Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) program methodology, which is used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal conditions of transport (NCT). The CIRFT consists of a U-frame test setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages connecting to a universal testing machine. The curvature SNF rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are clamped to the side connecting plates of the U-frame and used to capture deformation of the rod. The second portion of this report provides the latest CIRFT data, including data for the hydride reorientation test. The variations in fatigue life are provided in terms of moment, equivalent stress, curvature, and equivalent strain for the tested SNFs. The equivalent stress plot collapsed the data points from all of the SNF samples into a single zone. A detailed examination revealed that, at the same stress level, fatigue lives display a descending order as follows: H. B. Robinson Nuclear Power Station (HBR), LMK, and mixed uranium-plutonium oxide (MOX). Just looking at the strain, LMK fuel has a slightly longer fatigue life than HBR fuel, but the difference is subtle. The third portion of this report provides finite element analysis (FEA) dynamic deformation simulation of SNF assemblies . In a horizontal layout under NCT, the fuel assembly’s skeleton, which is formed by guide tubes and spacer grids, is the primary load bearing apparatus carrying and transferring vibration loads within an SNF assembly. These vibration loads include interaction forces between the SNF assembly and the canister basket walls. Therefore, the integrity of the guide tubes and spacer grids critically affects the vibration intensity of the fuel assembly during transport and must be considered when developing the multipurpose purpose canister (MPC) design for safe SNF transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigating the Influence of River Geomorphology on Human Presence Using Night Light Data: A Case Study in the Indus Basin

Human settlements have historically thrived near rivers due to enhanced navigation and trade, and the availability of water supply and resources. The use of night light data, representing economic activities, provides a novel approach to studying the interactions between human activity and rivers over time. Here, we use the Defense Meteorological Satellite Program (DMSP) stable night light data from 2000 to 2013 as a proxy for human presence and activities to quantify the statistical relationships between night light presence and intensity in the Indus Basin, Asia. We test how these data are affected by proximity to trunk channels and by channel type (single/multi-thread) in the study area. We find that night light presence is enhanced by 26% within a 0 to 5 km proximity range of the Indus River and its tributaries, relative to the basin as a whole. We interpret this to represent increased human presence and activity within this zone. However, the mean intensity is lower near the river and higher away from the river, signifying denser settlements, such as towns and cities, which are preferentially located away from the Indus and its tributaries. Moreover, the enhancement of lit pixels signifying human presence and activities is increased by 18% near single-thread sections of the Indus River, compared to segments of the Indus displaying multi-thread morphologies. We suggest that this is due to the enhanced stability of single-threaded channels, relative to mobile multi-threaded channel reaches. This study demonstrates how night lights are an important tool in studying the relationship between human presence and river dynamics in large catchments such as the Indus, and we suggest that these data will have an important role in assessing differential flood spatial and social vulnerability at a regional scale.

Environmental Sciences & Ecology↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FY 2016 Status Report: CIRFT Testing Data Analyses and Updated Curvature Measurements

This report provides a detailed description of FY15 test result corrections/analysis based on the FY16 Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) test program methodology update used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal transportation conditions. The CIRFT consists of a U-frame testing setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages to a universal testing machine. The curvature of rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are used and clamped to the side connecting plates of the U-frame to capture the deformation of the rod. The contact-based measurement, or three-LVDT-based curvature measurement system, on SNF rods has been proven to be quite reliable in CIRFT testing. However, how the LVDT head contacts the SNF rod may have a significant effect on the curvature measurement, depending on the magnitude and direction of rod curvature. It has been demonstrated that the contact/curvature issues can be corrected by using a correction on the sensor spacing. The sensor spacing defines the separation of the three LVDT probes and is a critical quantity in calculating the rod curvature once the deflections are obtained. The sensor spacing correction can be determined by using chisel-type probes. The method has been critically examined this year and has been shown to be difficult to implement in a hot cell environment, and thus cannot be implemented effectively. A correction based on the proposed equivalent gauge-length has the required flexibility and accuracy and can be appropriately used as a correction factor. The correction method based on the equivalent gauge length has been successfully demonstrated in CIRFT data analysis for the dynamic tests conducted on Limerick (LMK) (17 tests), North Anna (NA) (6 tests), and Catawba mixed oxide (MOX) (10 tests) SNF samples. These CIRFT tests were completed in FY14 and FY15. Specifically, the data sets obtained from measurement and monitoring were processed and analyzed. The fatigue life of rods has been characterized in terms of moment, curvature, and equivalent stress and strain..

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulator for an Accelerator-Driven Subcritical Fissile Solution System

LANL has developed a process to generate a progressive family of system models for a fissile solution system. This family includes a dynamic system simulation comprised of coupled nonlinear differential equations describing the time evolution of the system. Neutron kinetics, radiolytic gas generation and transport, and core thermal hydraulics are included in the DSS. Extensions to explicit operation of cooling loops and radiolytic gas handling are embedded in these systems as is a stability model. The DSS may then be converted to an implementation in Visual Studio to provide a design team the ability to rapidly estimate system performance impacts from a variety of design decisions. This provides a method to assist in optimization of the system design. Once design has been generated in some detail the C++ version of the system model may then be implemented in a LabVIEW user interface to evaluate operator controls and instrumentation and operator recognition and response to off-normal events. Taken as a set of system models the DSS, Visual Studio, and LabVIEW progression provides a comprehensive set of design support tools.

42 ENGINEERING↗

Laboratory Experiments and Instrument Intercomparison Studies of Carbonaceous Aerosol Particles

Aerosols containing black carbon (and some specific types of organic particulate matter) directly absorb incoming light, heating the atmosphere. In addition, all aerosol particles backscatter solar light, leading to a net-cooling effect. Indirect effects involve hydrophilic aerosols, which serve as cloud condensation nuclei (CCN) that affect cloud cover and cloud stability, impacting both atmospheric radiation balance and precipitation patterns. At night, all clouds produce local warming, but overall clouds exert a net-cooling effect on the Earth. The effect of aerosol radiative forcing on climate may be as large as that of the greenhouse gases, but predominantly opposite in sign and much more uncertain. The uncertainties in the representation of aerosol interactions in climate models makes it problematic to use model projections to guide energy policy. The objective of our program is to reduce the uncertainties in the aerosol radiative forcing in the two areas highlighted in the ASR Science and Program Plan. That is, (1) addressing the direct effect by correlating particle chemistry and morphology with particle optical properties (i.e. absorption, scattering, extinction), and (2) addressing the indirect effect by correlating particle hygroscopicity and CCN activity with particle size, chemistry, and morphology. In this connection we are systematically studying particle formation, oxidation, and the effects of particle coating. The work is specifically focused on carbonaceous particles where the uncertainties in the climate relevant properties are the highest. The ongoing work consists of laboratory experiments and related instrument inter-comparison studies both coordinated with field and modeling studies, with the aim of providing reliable data to represent aerosol processes in climate models. The work is performed in the aerosol laboratory at Boston College. At the center of our laboratory setup are two main sources for the production of aerosol particles: (a) two well-characterized source of soot particles and (b) a flow reactor for controlled OH and/or O3 oxidation of relevant gas phase species to produce well-characterized SOA particles. After formation, the aerosol particles are subjected to physical and chemical processes that simulate aerosol growth and aging. A suite of instruments in our laboratory is used to characterize the physical and chemical properties of aerosol particles before and after processing. The Time of Flight Aerosol Mass Spectrometer (ToF-AMS) together with a Scanning Mobility Particle Sizer (SMPS) measures particle mass, volume, density, composition (including black carbon content), dynamic shape factor, and fractal dimension. The–ToF-AMS was developed at ARI with Boston College participation. About 120 AMS instruments are now in service (including 5 built for DOE laboratories) performing field and laboratory studies world-wide. Other major instruments include a thermal denuder, two Differential Mobility Analyzers (DMA), a Cloud Condensation Nuclei Counter (CCN), a Thermal desorption Aerosol GC/MS (TAG) and the new Soot Particle Aerosol Mass Spectrometer (SP-AMS). Optical instrumentation required for the studies have been brought to our laboratory as part of ongoing and planned collaborative projects with colleagues from DOE, NOAA and university laboratories. Optical instruments that will be utilized include a Photoacoustic Spectrometer (PAS), a Cavity Ring Down Aerosol Extinction Spectrometer (CRD-AES), a Photo Thermal Interferometer (PTI), a new 7-wavelength Aethalometer and a Cavity Attenuated Phase Shift Extinction Monitor (CAPS). These instruments are providing aerosol absorption, extinction and scattering coefficients at a range of atmospherically relevant wavelengths. During the past two years our work has continued along the lines of our original proposal. We report on 12 completed and/or continuing projects conducted during the period 08/14 to 0814/2015. These projects are described in 17 manuscripts published in refereed journals.

54 ENVIRONMENTAL SCIENCES↗

Exploring the interplay of resilience and energy consumption for a task-based partial differential equations preconditioner

Here, we discuss algorithm-based resilience to silent data corruptions (SDCs) in a task-based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDCs. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ~51K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single and double bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamic voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

97 MATHEMATICS AND COMPUTING↗

Dynamics of an Adiabatic Point Reactor

Here, the problem of adiabatic excursions in a reactor is studied in general. We let the prompt temperature reactivity feedback be an unspecified function of temperature, ρ = ρ 0 = ρ 0 + f(T), where ρ is total reactivity, ρ 0 initial step reactivity and f(T) the feedback function. The similarity of the behavior of the reactor for different f(T) is established by means of a topological (qualitative) analysis. A quantitative asymptotic solution of the non-linear system of DE describing the reactor is presented. In delayed critical excursions, the delayed neutrons play a determining role. In the first part of a prompt excursion, the delayed-neutron source is nil; however this is not so in the second part, where it contributes appreciably to the excursion. These conclusions are shown to be valid in general, and allow us to write down almost directly the (approximate) quantitative solution of the non-linear system for any f(T). These results are correlated with the experimental data for the adiabatic excursions of a UO 2 core in SPERT I; in this case the (prompt) dependence of the reactivity on energy is of the form ρ = ρ 0 - 4.588 × 10 -4 E 0.74 .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Disruption of histone acetylation homeostasis reveals multilayered chromatin regulation for transcriptional resiliency

Background Epigenetic modifications, nucleosome occupancy, and three-dimensional chromatin architecture collectively create a multi-layered, highly interactive regulatory system for controlling genomic functionality. Dysregulation of epigenetic processes leads to a plethora of abnormalities including disease states. Therapies focused on epigenetic modulation can alter gene expression to correct dysfunction, though the perpetuation of these states and the relationships among chromatin regulatory layers is not well understood. Results Here, we investigated global and local chromatin structural and functional responses after acute histone deacetylase inhibitor treatment (suberoylanilide hydroxamic acid) in lung cancer cells across time. Treatment substantially increased global histone acetylation resulting in a pervasive but not distinctive signature. The spread of acetylation did not significantly impact global chromatin accessibility, and nucleosome remodeling largely occurred at finer scales in functionally relevant genomic regions. Indeed, both H 3 K 4 trimethylation, a mark of active transcription, and gene expression changes were altered in a controlled locus-specific manner, suggesting aberrant acetylation indirectly leads to balanced and bidirectional gene expression profiles from tighter regulation of other chromatin features. HDACi treatment induced (13%) genomic rearrangement in chromatin compartmentalization and moderate weakening of topologically associating domains. Conclusions Continuous wavelet analysis of these features demonstrates that scale-dependent, locus-specific factors influence the relationship between chromatin architecture and functional output, suggesting that regulation of transcription and nucleosome remodeling is not entirely (nor linearly) dependent upon large scale compartment exchange. Structural and functional responses are most pronounced early after treatment with partial persistence of differential local chromatin features and expression later in time; this highlights the plasticity of chromatin regulation, which may have implications for the efficacy of epigenetic treatments. These results demonstrate the effectiveness of multi-layered regulation of transcription: in resilient systems, disruption of one chromatin feature does not distort the regulation of other features in supporting a transcriptional program that allows for survival.

59 BASIC BIOLOGICAL SCIENCES↗

Programming interfacial energetic offsets and charge transfer in β-Pb 0.33 V 2 O 5 /quantum-dot heterostructures: Tuning valence-band edges to overlap with midgap states

Here, semiconductor heterostructures for solar energy conversion interface light-harvesting semiconductor nanoparticles with wide-band-gap semiconductors that serve as charge acceptors. In such heterostructures, the kinetics of charge separation depend on the thermodynamic driving force, which is dictated by energetic offsets across the interface. A recently developed promising platform interfaces semiconductor quantum dots (QDs) with ternary vanadium oxides that have characteristic midgap states situated between the valence and conduction bands. In this work, we have prepared CdS/β-Pb 0.33 V 2 O 5 heterostructures by both linker-assisted assembly and surface precipitation and contrasted these materials with CdSe/β-Pb 0.33 V 2 O 5 heterostructures prepared by the same methods. Increased valence-band (VB) edge onsets in X-ray photoelectron spectra for CdS/β-Pb 0.33 V 2 O 5 heterostructures relative to CdSe/β-Pb 0.33 V 2 O 5 heterostructures suggest a positive shift in the VB edge potential and, therefore, an increased driving force for the photoinduced transfer of holes to the midgap state of β-Pb 0.33 V 2 O 5 . This approach facilitates a ca. 0.40 eV decrease in the thermodynamic barrier for hole injection from the VB edge of QDs suggesting an important design parameter. Transient absorption spectroscopy experiments provide direct evidence of hole transfer from photoexcited CdS QDs to the midgap states of β-Pb 0.33 V 2 O 5 NWs, along with electron transfer into the conduction band of the β-Pb 0.33 V 2 O 5 NWs. Hole transfer is substantially faster and occurs at <1-ps time scales, whereas completion of electron transfer requires 5—30 ps depending on the nature of the interface. The differentiated time scales of electron and hole transfer, which are furthermore tunable as a function of the mode of attachment of QDs to NWs, provide a vital design tool for designing architectures for solar energy conversion. More generally, the approach developed here suggests that interfacing semiconductor QDs with transition-metal oxide NWs exhibiting intercalative midgap states yields a versatile platform wherein the thermodynamics and kinetics of charge transfer can be systematically modulated to improve the efficiency of charge separation across interfaces.

36 MATERIALS SCIENCE↗

Digitalization mapping and assessment process supporting ION strategic transformation activities

The existing fleet of commercial nuclear power plants (NPPs) are an important asset in the nation’s portfolio of electrical generating resources. Their continued safe and reliable operation are critical to providing a large source of carbon-free electricity to power the nation’s economy. The United States Department of Energy’s (DOE) Light Water Reactor Sustainability (LWRS) Program develops the scientific bases, methods, and tools, for the continued safe and economical operation of the nation's commercial NPPs. The Plant Modernization Pathway within LWRS Program focuses on providing guidance to industry on the full-scale implementation of modernization solutions for NPPs that significantly reduce the technical and financial risks associated with modernization. This research is focused on helping the nuclear industry understand how to digitize and digitalize their NPPs so that they can design their modernization solutions to be scalable, sustainable, and integrated both laterally and horizontally within their organization. That is, this research creates a digital transformation in NPPs by reshaping cultural mindsets and by identifying business efficiencies. In partnership with industry, and using four previously established guiding principles for digitalization, this research supported NPP modernization through assessing readiness for digitalization as a means to achieve integrated operations for nuclear. Specifically, this research created an assessment to review an entire organization’s work processes to gather information about the digitalization health of the plant. The assessment tools were administered to plant employees, and the results were used to develop a digitalization plan. The survey assessment identified the optimal candidate processes that would most benefit from a digitalization initiative which were revealed through analytical frameworks. One analysis calculated mean digitalization health indicator scores for all endorsed activities which allowed the researchers to rank and color code the results for easy identification. Individual health indicator scores are also provided, should our industry partner wish to understand these findings according to their own organizational priorities, business considerations and desired end-state. The results were also analyzed from the perspective that organizations are comprised of different types of innovators (e.g., generators, optimizers, conceptualizers, and implementers), which differentially affects the organization’s ability to comprehend and adapt to change (i.e., opportunities to innovate). Understanding the relative composition of innovator types at an NPP allows them to gather insights into the strengths and weaknesses they have in innovating how work is performed. For the utility that partnered with this research team, the results showed that implementers make up the largest portion of respondents and conceptualizers the smallest portion. Knowing the proportion of innovator types gave this organization insights on how they can effectively implement their innovation solutions. Additionally, the results were analyzed from a technical, economic, and risk perspective to identify and quantify work reduction opportunities (WROs). Recognizing that not all cost-saving opportunities are the same, a Technical, Economic and Risk Assessment (TERA) was performed to evaluate WROs to identify areas of greatest potential and lowest risk. The key results from TERA included a digitalization opportunity score for each activity, and a calculation of potential cost savings. These two outputs formed the bases for calculating a priority index and rank for the activities/processes assessed. From the prioritization calculations, TERA can then help the utility 1) decide what digitalization priorities to invest money in implementing and then 2) calculates how much should be invested in the digitalization initiatives selected to achieve cost savings and/or an acceptable return on investment. Last, onsite interviews revealed several inefficiencies in the standard work processes that occur cross-departmentally that are due to the absence of digitized and digitalized processes. Examples of these include time spent scanning paper documents and then uploading the documents electronically, obtaining signatures, and searching for desired information. This represents a digital but not digitalized process. Over 15 opportunities to improve work processes were identified through this multi-method digitalization assessment. The various analytical assessments used (e.g., TERA, digitalization health indicator scores), as well as discussions with the utility partner, corroborated that all the opportunities identified had a strong potential to make work processes more efficient and to improve overall performance of the NPP.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗