Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution multiplicity”

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

Single Grid Error Estimation for Neutron Transport Solvers

The method of nearby problems (MNP) is a solution verification technique that does not require the use of multiple spatial grids. To estimate spatial discretization error without requiring a high-fidelity spatial grid, an analytical curve fit is interpolated from the numerical solution. The residual between the curve fit solution and numerical solution is calculated and added as an additional source term to the governing equation. The nearby solution is estimated using the updated source term and boundary conditions to remain consistent with the curve fit interpolation. The nearby solution can be compared to the curve fit solution as a discretization error estimation while using a single spatial grid. Without the use of higher fidelity spatial grids, the MNP is able to approximate the spatial discretization error, a facet of solution verification. The application of the method of nearby problems is presented for one- and two-dimensional neutron transport problems for both fixed source and criticality problems on the spatial variable. The fixed source results demonstrate the effectiveness of nearby problems for spatial error identification using the discrete ordinates method. Criticality results are shown to identify area of high spatial error for the C5G7 problem as well as for the discrete ordinates solver. A novel approach of combining the capabilities of Monte Carlo with the discrete ordinates nearby problems is presented for one- and two-dimensional fixed source problems. In conclusion, the MNP demonstrates its effectiveness at identifying spatial error on a single structured grid with a wide variety of neutron transport problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reducing the computational burden of a microgrid energy management system

As renewable technology advances and decreases in cost, microgrids are becoming an appealing means of distributed generation both for isolated communities and integrated with existing electrical grid systems. Due to their small size, however, microgrids may have financial limitations which preclude them from using commercial software to optimize control of their assets. Open-source optimization solvers are a viable alternative, but increase computation time. This work expands on a rolling horizon optimization framework for economic dispatch within an existing residential microgrid located in Hoover, Alabama. The microgrid has an open-source solver requirement and a need for quick solution time on a rolling horizon as opposed to a day-ahead commitment. We present a method of reducing integer variables by relaxation which completes two goals: reduction in computation time for real-time operations, and reduction in daily operational cost for the microgrid. Seasonal data for load and photovoltaic (PV) power was also collected from the microgrid to facilitate simulation testing. Computation time was successfully reduced using multiple variations of the relaxation method, while obtaining solution quality with operational cost similar to or better than the original model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Contrasting Trivalent Lanthanide and Actinide Complexation by Polyoxometalates via Solution-State NMR

Deciphering the solution chemistry and speciation of actinides is inherently difficult due to radioactivity, rarity, and cost constraints, especially for transplutonium elements. In this context, the development of new chelating platforms for actinides and associated spectroscopic techniques is particularly important. In this study, we investigate a relatively overlooked class of chelators for actinide binding, namely, polyoxometalates (POMs). We provide the first NMR measurements on americium–POM and curium–POM complexes, using one-dimensional (1D) 31 P NMR, variable-temperature NMR, and spin-lattice relaxation time (T 1 ) experiments. The proposed POM–NMR approach allows for the study of trivalent f-elements even when only microgram amounts are available and in phosphate-containing solutions where f-elements are typically insoluble. The solution-state speciation of trivalent americium, curium, plus multiple lanthanide ions (La 3+ , Nd 3+ , Sm 3+ , Eu 3+ , Yb 3+ , and Lu 3+ ), in the presence of the model POM ligand PW 11 O 39 7– was elucidated and revealed the concurrent formation of two stable complexes, [M III (PW 11 O 39 )(H 2 O) x ] 4– and [M III (PW 11 O 39 ) 2 ] 11– . Interconversion reaction constants, reaction enthalpies, and reaction entropies were derived from the NMR data. The NMR results also provide experimental evidence of the weakly paramagnetic nature of the Am 3+ and Cm 3+ ions in solution. Furthermore, the study reveals a previously unnoticed periodicity break along the f-element series with the reversal of T 1 relaxation times of the 1:1 and 1:2 complexes and the preferential formation of the long T 1 species for the early lanthanides versus the short T 1 species for the late lanthanides, americium, and curium. Furthermore, given the broad variety of POM ligands that exist, with many of them containing NMR-active nuclei, the combined POM–NMR approach reported here opens a new avenue to investigate difficult-to-study elements such as heavy actinides and other radionuclides.

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

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗
DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗
DOE OSTI2022

Proximity Portability and in Transit , M-to-N Data Partitioning and Movement in SENSEI [Book Chapter]

In high-performance parallel in situ processing, the term in transit processing refers to those configurations where data must move from a producer to a consumer that runs on separate resources. In the context of parallel and distributed computing on an HPC platform one of the central challenges is to determine a mapping of data from producer ranks to consumer ranks. This problem is complicated by the heterogeneity that arises in producer-consumer pairs, such as when producer and consumer codes have different levels of concurrency, different scaling characteristics, or different data models. The resulting mapping and movement of data from M producer to N consumer ranks can have a significant impact on aggregate application performance, particularly when the data consumer requires only a subset of the overall data for its task. This chapter focuses on the design considerations that underlie SENSEI’s implementation to this challenging problem. These design considerations extend the core SENSEI architecture and include ideas like the need to accommodate flexibility in the choice of different partitioning methods, the ability for a data consumer to request and receive only the subset of data needed for its particular operation, and the ability to leverage any of several different data transport tools. The idea of proximity portability, being able to use different data transport methods as part of an in transit workflow, is illustrated through the use of three different transport layers where switching from one transport tool to another is accomplished with only a configuration file change. Here, the chapter also includes a performance analysis summary showing the performance gains that are possible in terms of multiple metrics, such as memory footprint, time to solution, and amount of data moved, when using optimized partitioners in an in transit setting, gains that are made possible by the implementation shaped by specific design considerations.

Bethel, E. Wes↗
DOE OSTI2023

Systematic mapping of the conformational landscape and dynamism of soluble fibrinogen

Background: Fibrinogen is a soluble, multisubunit, and multidomain dimeric protein, which, upon its proteolytic cleavage by thrombin, is converted to insoluble fibrin, initiating polymerization that substantially contributes to clot growth. Fibrinogen contains numerous, transiently accessible “cryptic” epitopes for hemostatic and immunologic proteins, suggesting that fibrinogen exhibits conformational flexibility, which may play functional roles in its temporal and spatial interactions. Hitherto, there have been limited integrative approaches characterizing the solution structure and internal flexibility of fibrinogen. Methods: Here, utilizing a multipronged, biophysical approach involving 2 solution-based techniques, temperature-dependent hydrogen-deuterium exchange mass spectrometry and small angle X-ray scattering, corroborated by negative stain electron microscopy, we present a holistic, conformationally dynamic model of human fibrinogen in solution. Results: Our data reveal 4 major and distinct conformations of fibrinogen accommodated by a high degree of internal protein flexibility along its central scaffold. We propose that the fibrinogen structure in the solution consists of a complex, conformational landscape with multiple local minima. This is further supported by the location of numerous point mutations that are linked to dysfibrinogenemia and posttranslational modifications, residing near the identified fibrinogen flexions. Conclusion: Further, this work provides a molecular basis for the structural “dynamism” of fibrinogen that is expected to influence the broad swath of its functionally diverse macromolecular interactions and fine-tune the structural and mechanical properties of blood clots.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Charge-density based convolutional neural networks for stacking fault energy prediction in concentrated alloys

A descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. Here, we illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R 2 of 6.18 mJ/m 2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MULTI-MODAL global surveillance methodology for predictive and on-demand characterization of localized processes using cube satellite platforms and deep learning techniques

This paper presents the work completed towards the development of a multi-modal global surveillance methodology using cube satellite (CubeSat) platforms and novel data analysis techniques. A CubeSat system equipped with adequate sensors and data analytics capabilities can autonomously characterize various phenomena of interest on the Earth’s surface. CubeSats are advantageous over conventional satellites in certain remote monitoring applications because of their reduced construction costs (due to the availability of commercially-off-the-shelf components) and are easier to launch. The CubeSat surveillance system developed in this paper focused on phenomena of interest surrounding the nuclear fuel cycle in support of nuclear non-proliferation and emergency response. To observe the phenomena, a constellation of 3U and 6U CubeSats deployed from the ISS with adequate components was chosen. Four different sensor configurations were identified for remote sensing: panchromatic/multispectral in the visible and near-infrared spectrum, multispectral in infrared spectrum, hyperspectral in infrared spectrum, and multispectral in ultraviolet spectrum. While a panchromatic/multispectral sensor configuration has CubeSat flight heritage at the required spatial resolutions, the other three sensor types need future 3 development to meet signature and system requirements. Once each sensor onboard the CubeSat system collects data on a target of interest, the onboard computers would then apply the deep learning-based characterization methodology developed in this paper to identify phenomena. Four surrogate datasets containing representative simplified “images” were created for each sensor type to train the characterization methodology. A convolutional neural network was applied to each dataset and produced recall rates for the phenomena between 89.7% - 99.3% and precision rates between 92.3% - 99.9%. Each phenomenon’s presence probability from each network is then combined into a final characterization solution for a target area. This paper covers multiple interdisciplinary areas to develop the foundation for a CubeSat surveillance system focused on phenomena surrounding the nuclear fuel cycle.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Ion Distribution and Cation Exchange at Mica–Electrolyte Interfaces Probed with Deep Potential Molecular Dynamics

Here, we investigate the Stern layer structure and cation exchange mechanism at muscovite mica-electrolyte interfaces using nanosecond timescale molecular dynamics simulations based on deep neural network interatomic potentials trained on Density Functional Theory (DFT) data. Focusing on mica with exposed surface K + interfaced with aqueous NaCl and mica with surface Na + interfaced with KCl solution, we find that K + remains predominantly in inner-sphere configurations, while Na + exhibits notable populations in outer-sphere states. Most importantly, our simulations show that contact with an electrolyte solution results in the co-adsorption of multiple cation species, making the mica surface locally overcharged and thus reshaping the cation speciation in a manner that enhances the tendency of neighboring surface cations to desorb. These findings are consistent with recent experimental observations that co-adsorption of different cation species induces changes in cation speciation and slow kinetics of cation exchange at the muscovite-water interface, providing a basis for their detailed understanding.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2024

3D Integrated Sensing Solutions

We will present the proposed research project to develop 3D-integrated sensors using advanced manufacturing capability for novel sensors, heterogeneously integrated with energy efficient readout circuits. The lack of precision timing in particle tracking detectors and the absence of low power, high throughput communications solutions to read them out limits progress in multiple fields of fundamental science. The aim of the project is to develop technology to enable large-scale particle detectors with 3D-integrated ASIC designs to simultaneously achieve 10 μm position resolution and 10 ps precision timing, with low-power consumption and high throughput rates. This research program leverages the unique combination of facilities and cross-disciplinary expertise of scientists and engineers at SLAC, FNAL, and LLNL and industrial partners.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Understanding Twinning and Deformation in High Entropy Alloys

A combination of high strength and high ductility has been observed in multi-principal element alloys due to twin formation attributed to low stacking fault energy (SFE). In the pursuit of low SFE alloys, a key bottleneck is the lack of understanding of the composition–SFE cor- relations that would guide tailoring SFE via alloy composition. Using density functional theory (DFT), we show that dopant radius, which have been postulated as a key descriptor for SFE in dilute alloys, does not fully explain SFE trends across different host metals. Instead, charge density is a much more central descriptor. It allows us to (1) explain contrasting SFE trends in Ni and Cu host metals due to various dopants in dilute concentrations, (2) explain the large SFE variations observed in the literature even within a given alloy composition due to the nearest neighbor environments in “model” concentrated alloys, and (3) develop a machine learning model that can be used to predict SFEs in multi-elemental alloys. This model opens a possibility to use charge density as a descriptor for predicting SFE in alloys. Furthermore, a descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. We illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R2 of 6.18 mJ/m2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

iNovo479: Metabolic Modeling Provides a Roadmap to Optimize Bioproduct Yield from Deconstructed Lignin Aromatics by Novosphingobium aromaticivorans

Lignin is an abundant renewable source of aromatics and precursors for the production of other organic chemicals. However, lignin is a heterogeneous polymer, so the mixture of aromatics released during its depolymerization can make its conversion to chemicals challenging. Microbes are a potential solution to this challenge, as some can catabolize multiple aromatic substrates into one product. Novosphingobium aromaticivorans has this ability, and its use as a bacterial chassis for lignin valorization could be improved by the ability to predict product yields based on thermodynamic and metabolic inputs. In this work, we built a genome-scale metabolic model of N. aromaticivorans, iNovo479, to guide the engineering of strains for aromatic conversion into products. iNovo479 predicted product yields from single or multiple aromatics, and the impact of combinations of aromatic and non-aromatic substrates on product yields. We show that enzyme reactions from other organisms can be added to iNovo479 to predict the feasibility and profitability of producing additional products by engineered strains. Thus, we conclude that iNovo479 can help guide the design of bacteria to convert lignin aromatics into valuable chemicals.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

LTE Electrolyzer Data Collection

The goal for NREL is to collect, develop and publish performance metrics relative to low temperature electrolyzer installations. This will be done through the development of: Secure storage solution to house the collection of data from multiple projects Standardization of data to be collected and analyzed. This will be done using data templates developed with the help of partners involved with electrolyzer installations. Analysis that produces metrics of interest for all stakeholders Aggregation of results from multiple projects to view industry progress as a whole Publication of aggregated results in the form of composite data products (CDPs) Collaboration with Idaho National Lab and their work with high temperature electrolyzer installations will enable efficient use of storage and analysis tools.

data↗
DOE OSTI2024

Practical challenges of model predictive control (MPC) for grid interactive small and medium commercial buildings

To the urgent call for mitigating climate change, substantial initiatives have been undertaken to deploy grid-interactive heating, ventilation, and air-conditioning (HVAC) controls, such as model predictive control (MPC) for buildings. These efforts typically aim to curtail peak energy demand, shift load and enhance overall energy efficiency. With the recent development of low-cost MPC technologies that don’t require extensive instrumentation or manual modeling, small and medium commercial buildings (SMCBs), which rarely utilize advanced HVAC control systems, have become candidates for grid-interactive efficient buildings (GEBs). However, despite the potential benefits and maturity of the technology itself, several practical challenges remain in real-world implementation. In this paper, we share the practical challenges that we have encountered in implementing and testing three types of MPC solutions (ON/OFF unit, dualfuel, and VRF systems) on multiple SMCB sites. We describe the MPC deployment process and discuss the lessons learned. The site selection, eligibility, and retrofit availability (e.g., utility price structure, thermostat communications, etc.) are the main discussion points at the beginning of the project. Also, the modeling automation and the best practices for interacting with endusers and handling erroneous situations are presented for successful operations.

woo Ham, Sang↗
DOE OSTI2022

Land‐based climate solutions for the United States

Abstract Meeting end‐of‐century global warming targets requires aggressive action on multiple fronts. Recent reports note the futility of addressing mitigation goals without fully engaging the agricultural sector, yet no available assessments combine both nature‐based solutions (reforestation, grassland and wetland protection, and agricultural practice change) and cellulosic bioenergy for a single geographic region. Collectively, these solutions might offer a suite of climate, biodiversity, and other benefits greater than either alone. Nature‐based solutions are largely constrained by the duration of carbon accrual in soils and forest biomass; each of these carbon pools will eventually saturate. Bioenergy solutions can last indefinitely but carry significant environmental risk if carelessly deployed. We detail a simplified scenario for the United States that illustrates the benefits of combining approaches. We assign a portion of non‐forested former cropland to bioenergy sufficient to meet projected mid‐century transportation needs, with the remainder assigned to nature‐based solutions such as reforestation. Bottom‐up mitigation potentials for the aggregate contributions of crop, grazing, forest, and bioenergy lands are assessed by including in a Monte Carlo model conservative ranges for cost‐effective local mitigation capacities, together with ranges for (a) areal extents that avoid double counting and include realistic adoption rates and (b) the projected duration of different carbon sinks. The projected duration illustrates the net effect of eventually saturating soil carbon pools in the case of most strategies, and additionally saturating biomass carbon pools in the case of forest management. Results show a conservative end‐of‐century mitigation capacity of 110 (57–178) Gt CO 2 e for the U.S., ~50% higher than existing estimates that prioritize nature‐based or bioenergy solutions separately. Further research is needed to shrink uncertainties, but there is sufficient confidence in the general magnitude and direction of a combined approach to plan for deployment now.

54 ENVIRONMENTAL SCIENCES↗