Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numberical Integration”

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 73 records · Page 4

Automatic Generation of Event Trees and Fault Trees: A Model-Based Approach

In the past few decades, the increasing complexity of modern engineering systems has been driven by the integration of a large number of components whose operations may involve many disciplines (e.g., thermal hydraulics, plant operations, cybersecurity). Most computational tools used by industry and regulators for system safety and reliability assessments are still based on the traditional fault tree (FT) and event tree (ET) approach, which may not be able to capture complex interactions among system constituents. The use of simulation tools has widely increased in the past few decades to improve the fidelity of the reliability and safety analyses. However, the direct use of simulation tools as part of dynamic probabilistic risk assessment (DPRA) methods is not getting traction since (1) modeling the whole system under consideration with DPRA methods may be computationally expensive and unnecessary, and (2) the manual integration of DPRA models into existing state-of-practice probabilistic risk assessment models (i.e., based on FTs and ETs) can be time consuming and prone to errors. Here, in this paper we propose a procedure to overcome this limitation by presenting several algorithms designed to automatically construct subsystem ETs and FTs from DPRA methods for integration into an existing ET/FT system model.

97 MATHEMATICS AND COMPUTING↗

Layer-by-layer growth of graphene oxide multilayers using robust interlayer linking chemistry. 1. Zr-bisphosphates

Graphene oxide (GO) holds substantial interest because of its utility in applications ranging from chemical sensing to catalysis and energy storage. We report on the layer-by-layer formation of graphene oxide phosphate (P-GO) multilayers on silica and oxidized silicon supports, with the layers connected via Zr-bisphosphate (ZP) linking chemistry. Layers were grown either directly by ZP linkages between P-GO sheets, or with the use of 1,4-phenylene bisphosphate (BP1) or [1,1′-biphenyl]4,4′-diyl bisphosphate (BP2). The layers have been characterized using optical null ellipsometry, X-ray photoelectron spectroscopy (XPS), UV-visible absorption spectroscopy and scanning electron microscopy (SEM). There is a linear dependence of assembly thickness on the number of layers deposited and the integrated area of the optical absorbance bands also increases linearly with number of layers deposited. XPS data provides information on interlayer linking stoichiometry. SEM images provide insight into the morphology of the adlayers, suggesting the structure and length of the interlayer linking moieties used in the multilayer assembly play a significant role in the organization of the resulting system.

Mukherjee, Neelanjana [Michigan State University, ↗

AUTOMATIC GENERATION OF EVENT TREES AND FAULT TREES: A MODEL-BASED APPROACH

In the past few decades, increasing complexity in modern engineering systems has been driven by the integration of a large number of components and by the fact that the system operations involve many disciplines (e.g., thermal-hydraulics, plant operations, cyber-security). Current safety/reliability modeling approaches to such systems are labor intensive, difficult to learn, and rely heavily on simplistic Boolean logic to depict failure propagation and accident progression. While these methods serve well for simple systems (i.e., linear causal systems with limited small inter- and intra-system interactions), their results are difficult to verify when modeling complex systems (typically performed through the extensive use of modeling assumptions). The development of new methods is addressed to meet these challenges through a model-based system engineering (MBSE) lens. Under MBSE philosophy, every aspect of the system (form or function) is represented by a model that completely characterizes its architecture or behavior. MBSE approach greatly improves the management of design, analysis and verification of complex systems. An integration of Dynamic Probabilistic Risk Assessment (DPRA) methods with MBSE models is proposed to perform safety/reliability analyses of engineering systems. In particular, MBSE representation of the system (performed using Systems Modeling Language [SysML]) is coupled with DPRA methods to automatically generate event trees and fault trees.

97 - MATHEMATICS AND COMPUTING↗

Comparison of Thermal Management Approaches for Integrated Traction Drives in Electric Vehicles

The continuous push to increase power densities of electric vehicle (EV) traction drive systems necessitates combining electric motor and power electronics into one unit. A single, compact traction drive unit with fewer interconnecting components also facilitates fast, automated assembly of electric vehicles, driving production costs down and enabling wider adoption of EVs. There are a number of challenges associated with the integration of power electronics with the electric machine, including thermal management of the combined traction drive system. However, one important benefit of integration from the thermal management system perspective is the potential for using a single fluid loop instead of two separate cooling systems for the electric machine and the power electronics/inverter. This paper reviews several integration approaches and, employing finite element analysis (FEA), compares thermal management solutions for the combined electric machine and power electronics systems. Namely, three different scenarios are modeled: (1) independent component (motor and power electronics) cooling, which is compared to the combined cooling system approach for (2) radially and (3) axially integrated power electronics modules into the motor enclosure. Temperature distributions for selected thermal loads and thermal resistances from the key heat-generating components to the cooling fluid are compared for each scenario.

47 OTHER INSTRUMENTATION↗

A Blockchain Based Co-Simulation Framework for Integrating DERs to Wholesale Electricity Markets

As the number of distributed energy resources (DERs) continue to increase across energy-delivery systems, there remains a need for integrating their capabilities into traditional grid operations. In this paper, a blockchain-based solution is proposed to facilitate FERC's Order No. 2222 implementations. The presented use-case enables small-scale DERs to participate in wholesale market operations through DER aggregators, while also enabling local distribution system operators to enforce distribution system constraints in a secure and traceable manner. The presented use case is built around the Unified Testing Platform (UTP) being developed as a part of the Blockchain for Optimized Security and Energy Management (BLOSEM) project. This is a multi-lab effort intended to simplify the deployment of blockchain-powered grid solutions by enabling the integration of simulation tools, and blockchain technologies through the use of system-agnostic interfaces that provide a modular, interoperable, and reusable connectivity layer.

blockchain interoperability↗

Scalable Predictive Control and Optimization for Grid Integration of Large-Scale Distributed Energy Resources: Preprint

Integration of a large number of distributed energy resources (DERs) into the power grid needs a scalable power balancing method. We formulate the power balancing problem as a look-ahead optimization problem to be solved sequentially by a power distribution system aggregator based on a model predictive control (MPC) framework. Solving large-scale look-ahead control problem requires proper configuration of the control steps. In this paper, to solve large-scale control problems, we propose a variable time granularity where control time steps nearby the current control step have finer resolutions. The aggregator objective includes maximization of power production revenue and minimization of power purchasing expense, renewable power curtailment, and mileage costs for energy storage and electric vehicle (EV) charging stations while satisfying system capacity and operational constraints. The control problem is formulated as a mixed-integer linear program (MILP) and solved using the XpressMP solver. We perform simulations considering a copper plate representation of a large distribution network consisting of 2507 devices (controllable DERs) including curtailable photovoltaics (PVs), energy storage batteries, EV charging stations, and buildings with heating, ventilation, and air conditioning units (HVACs). We show the effectiveness of the proposed approach in managing DERs interactively for maximum energy trading profit and local supply-demand power balancing. Finally, we demonstrate that the proposed method outperformed other benchmark controllers regarding computation time without compromising operational performance.

DER↗

Nanophotonic waveguide chip-to-world beam scanning

A seamless chip-to-world photonic interface enables broad advancements in optical ranging, display, communication, computation and quantum information science. The ideal solution enables two-dimensional scanning of a diffraction-limited beam from anywhere on a photonic integrated circuit to a large number of resolvable spots. Current beam-scanning technologies are limited by a fundamental trade-off: photonic-integrated-circuits with diffractive optics offer scalability but have poor mode quality, whereas inertially limited micromechanical scanners provide high-quality beams but lack scalable integration. Here we report a photonic ski-jump—a nanoscale waveguide monolithically integrated on a piezoelectric cantilever—to overcome these limitations. It passively curls ~90° out-of-plane within a less-than-0.1 mm 2 footprint, emits a submicrometre, broadband diffraction-limited beam, and exhibits kilohertz-rate mechanical resonances with quality factors of over 10,000. Fabricated in a volume complementary metal–oxide–semiconductor (CMOS) foundry, our device enables scalable two-dimensional beam scanning. Driven on-resonance at CMOS-level voltages, it achieves a footprint-adjusted spot rate of 68.6 mega spots s –1 mm–², exceeding state-of-the-art micro-electro-mechanical systems mirrors by more than 50-fold, which is sufficient for one million pixels at 100 Hz from an approximately 1.5 mm diameter footprint. We demonstrate full-colour image and video projection, and single-photon initialization and readout from silicon vacancy centres in diamond. Finally, by demonstrating uniformity across a 64 ski-jump array, we establish a pathway to achieving greater than one gigaspot resolution at kilohertz rates within a sub-5-cm-diameter footprint, creating a seamless optical pipeline between integrated photonic processors and the free-space world.

Displays↗

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗

Large-Scale Inverter Integration in Bulk Power Grids Using Heterogeneous Grid-Forming Control Strategies

The conventional generation within the grid is gradually being replaced by inverter based resources. However, the integration of inverters into realistic transmission grid models has not been thoroughly explored. To address this gap, this work develops an automation framework aimed at facilitating the integration of a large number of inverter-based resources into a large-scale grid. Two types of grid-forming controlled inverters—namely droop-controlled inverters and virtual synchronous machine-controlled inverters—are developed for use in positive-sequence simulation packages. Various penetration levels ranging from 23% to 100% of grid-forming inverters are tested in a grid setup in the US Western Interconnection. The dynamic performance of inverters under different control schemes is compared in response to various disturbances.

Lyu, Xue [BATTELLE (PACIFIC NW LAB)]↗

Inverter Intensive Hybrid Power Plant Modeling with Small-Signal Stability Augmentation through Flexible Operation Mode Transition

Hybrid power plants (HPPs) prompt the penetration of inverter-based renewable energy sources (RES) in transmission systems; however, given their low-inertia nature, HPPs are dominated by power electronic inverters, so there are inevitable challenges in system stability when increasing numbers of HPPs are integrated into the modern power grids. To boost the penetration level of HPPs without jeopardizing system stability, it is desirable to equip them with operational characteristics (i.e., grid-forming capabilities) that are comparable to those of conventional power plants dominated by synchronous generators (SGs). In this paper, a holistic model of inverter intensive HPP is derived and a bi-level hierarchical control is developed to allow HPPs to flexibly switch among the designed operation modes (i.e., P-Q, P-V, and isochronous modes). Compared to SGs, which have limited controllability, the operation mode of each HPP could vary as requested. Such flexible mode transitions could be integrated into the secondary plant-level control and be leveraged as an additional control degree to further augment system stability. Further, the system small-signal stability margin is quantified with varying HPP operation modes. More importantly, modal analysis is thereby conducted to quantify the impacts of mode transition on the system oscillatory modes. The effectiveness of the proposed HPP bi-level hierarchical control is verified using extensive case studies based on the simplified real-world island power grid, and the results validate that the system small-signal stability margin can be enhanced with the additional degree of control flexibility enabled by HPP operation mode transition. The real-time hardware-in-the-loop (HIL) results are also provided to verify the proposed method.

grid-forming control↗

Use of Sobol’ Variance-Based Global Sensitivity Analysis and Multidimensional Legendre Polynomial Fitting for Reduced Order Modeling

Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.

97 MATHEMATICS AND COMPUTING↗

SAM Linear Fresnel Model, Project B (CRADA Final Report)

Objective performance and economic modeling of solar thermal plants is of keen interest to many EPRI funders. One solar thermal technology that is not currently available to model in any non-vendor, non-proprietary tool is linear Fresnel. This technology has garnered enough interest from EPRI funders to merit investing in a tool to objectively model its performance. Early in 2010, EPRI performed a comparison of modeling solar thermal power plants using the IPSEPRO, CNRS and Solar Advisor (SAM) tools. After completing this effort, EPRI decided to adopt NREL’s Solar Advisor Model as its default modeling tool based in part on user friendliness, flexibility, number of technologies covered, integrated financial model and ease of running sensitivities. Furthermore, it was recognized that NREL continues to invest considerable time and resources into improving capabilities and functionality of the model.

14 SOLAR ENERGY↗

Integrated Tiger Series (ITS) Code Validation and Uncertainty Quantification for Electron Number Albedo in Various Materials

In the 1970’s and 1980’s, researchers at Sandia National Laboratories produced electron albedo data for a range of materials. Since that time, the electron albedo data has been used for a wide variety of purposes including the validation of Monte Carlo electron transport codes. This report was compiled to examine the electron albedo experiment results in the context of Integrated Tiger Series (ITS) validation. The report presents tables and figures that could provide insight into the underlying model form uncertainty present in the ITS code. Additionally, the report provides data on potential means to reduce these model form errors by highlighting potential refinements in the cross-section generation process.

97 MATHEMATICS AND COMPUTING↗

End-gas autoignition fraction and flame propagation rate in laser-ignited primary reference fuel mixtures at elevated temperature and pressure

Knock in spark-ignited (SI) engines is initiated by autoignition of the unburned gasses upstream of spark-ignited, propagating, turbulent premixed flames. Knock propensity of fuel/air mixtures is typically quantified using research octane number (RON), motor octane number (MON), or methane number (MN; for gaseous fuels), which are measured using single-cylinder, variable compression ratio engines. In this study, knock propensity of SI fuels was quantified via observations of end-gas autoignition (EGAI) in unburned gasses upstream of laser-ignited, premixed flames at elevated pressures and temperatures in a rapid compression machine. Stoichiometric primary reference fuel (PRF; n-heptane/isooctane) blends of varying reactivity (50 ≤ PRF ≤ 100) were ignited using an Nd:YAG laser over a range of temperatures and pressures, all in excess of 545 K and 16.1 bar. Laser ignition produced outwardly-propagating premixed flames. High-speed pressure measurements and schlieren images indicated the presence of EGAI. The fraction of the total heat release attributed to EGAI (i.e., EGAI fraction) varied with fuel reactivity (i.e., octane number) and the time-integrated temperature of the end-gas prior to ignition. Flame propagation rates, which were measured using schlieren images, were only weakly correlated with octane number but were affected by turbulence caused by variation in piston timing. Under conditions of low turbulence, measured flame propagation rates approached one-dimensional premixed laminar flame speed computations performed at the same conditions. Experiments were simulated with a three-dimensional CONVERGE™ model using reduced chemical kinetics (121 species, 538 reactions). The simulations accurately captured the measured flame propagation rates, as well as the variation in EGAI fraction with fuel reactivity and time-integrated end-gas temperature. The simulations also revealed low-temperature heat release as well as formaldehyde and hydrogen peroxide formation in the end-gas upstream of the propagating flame, which increased the temperature and degree of chain branching in the end-gas, ultimately leading to EGAI.

42 ENGINEERING↗

Ranking Biological Features in Soil-Based Microbial Multi-Omics Data with Integration Modeling

Distinguishing the most important features (e.g. proteins, metabolites, etc.) per group (e.g. control and treatment) is a critical challenge in feature-rich multi-omics experiments, especially in soil data. Traditional feature identification and ranking approaches, such as differential expression, are based on single omics and thus not directly translatable to multi-omics experiments. Here, 5 multi-omics integration models (DIABLO, JACA, MOFA, MultiMLP, and SLIDE) that were not explicitly built for soil data applications were tested using a soil-based multi-omics experiment. The data were obtained from an experimental setup of an autoclaved soil system inoculated with 8 bacteria and using chitin as the carbon source and including samples collected at 0- (control), 4-, 8-, and 12-weeks post-inoculation. The omics data included metaproteomics, 16S rRNA sequencing, and LC-MS/MS metabolomics (in positive and negative mode). Each multi-omics integration model was implemented, and top features were compared to differential univariate statistics per omic type, demonstrating that integration approaches cut the potential number of top features from 2957 identified by differential statistics to 13-224 (a 99.6% to 92.4% reduction). Interestingly, most top features across integration models were not shared; though, scaling and averaging ranks across models shared similar patterns. This work highlights the usefulness of multi-omics integration models in soil-based microbial studies and the power of using multiple integration models together to interpret results.

54 ENVIRONMENTAL SCIENCES↗

Sensor Reduction for Diversion Detection in a Realistic Heat Pipe Microreactor Using Supervised Machine Learning

Microreactors are designed as a smaller, cheaper, and safer alternative to traditional nuclear power plants. Their non-traditional characteristics and prospect of mass production and deployment will likely require new approaches to nuclear safeguards. The primary proliferation concern with microreactors is the diversion of fuel material. Such diversion may produce measurable defects in key physical attributes like neutron flux, which may in turn be detectable using machine learning models. Preliminary work has demonstrated this ability for modeled nominal and diversion scenarios using large quantities of energy integrated neutron flux data. In practice, the number of available sensors for such measurements will be limited and energy integrated flux information will not be available. This work explores the ability of tree-based gradient boosted ensemble models to classify a given microreactor core is nominal or diversion, and determine the number of fuel pins diverted in the case of diversion with reduced numbers of sensors and more realistic detector responses. Classification accuracy of greater than 98% and regression errors as low as 5% of the total number of fuel pins were achieved with as few as 15 sensors, compared to 99% and 4.1% with a maximum of 240 sensors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗