Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix power”

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

Enhancing ACPF Analysis: Integrating Newton-Raphson Method with Gradient Descent and Computational Graphs

This paper presents a new method for enhancing Alternating Current Power Flow (ACPF) analysis. The method integrates the Newton-Raphson (NR) method with Enhanced-Gradient Descent (GD) and computational graphs. The integration of renewable energy sources in power systems introduces variability and unpredictability, and this method addresses these challenges. It leverages the robustness of NR for accurate approximations and the flexibility of GD for handling variable conditions, all without requiring Jacobian matrix inversion. Furthermore, computational graphs provide a structured and visual framework that simplifies and systematizes the application of these methods. The goal of this fusion is to overcome the limitations of traditional ACPF methods and improve the resilience, adaptability, and efficiency of modern power grid analyses. We validate the effectiveness of our advanced algorithm through comprehensive testing on established IEEE benchmark systems. Furthermore, our findings demonstrate that our approach not only speeds up the convergence process but also ensures consistent performance across diverse system states, representing a significant advancement in power flow computation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A highly stable Cu(OH) 2 -Poly(vinyl alcohol) nanocomposite membrane for dramatically enhanced direct borohydride fuel cell performance

Nano-additive aggregation, limited performance, and unclear modification mechanisms are the main obstacles in developing nanocomposite anion exchange membranes (AEMs). In this work, for the first time, an effective and highly stable Poly(vinyl alcohol) (PVA)-based AEMs with dispersive Cu(OH)$_2$ nanoclusters (Cu-AEMs) are prepared by a simple and eco-friendly three-step method: ‘CuCl$_2$ doping-casting-KOH immersing’. The doped Cu$^{2+}$ ions chemically combine with OH- ions to form anionic conductive Cu(OH)2 nanoclusters intermediated by attaching resins. The PVA skeletons wrap around Cu(OH)$_2$ while the hydroxyl groups expose to bulk water, forming Cu(OH)$_2$-PVA complex, which avoids nano-additive aggregation, increases anionic channels, and strengthens additive-matrix connection. A direct borohydride fuel cell using Cu-AEM with 0.56 wt % CuCl$_2$ possess the highest power density of 403.3 mW cm$^{-2}$ at 60 °C and a life span of over 200 h. The high-performance and durability come from the unique structure of Cu(OH)$_2$ nanoclusters-PVA complex. ‘Vehicle’ theory is considered to be the dominant mechanism for enhancing such nanocomposite Cu-AEMs. This work demonstrates a new concept for preparing the stable AEMs toward high-performance fuel cells. The synthetic chemistry involved can be broadly extended for fabricating versatile AEMs.

25 ENERGY STORAGE↗
DOE OSTI2025

Second-order renormalized Hamiltonian of Yukawa theory

Using the renormalization group procedure for effective particles we calculate the effective Hamiltonians in the theory of a fermion field coupled to a scalar field via the Yukawa interaction. The theory is renormalized by the addition of counterterms. Necessary counterterms are determined by computing matrix elements of the effective Hamiltonian. All calculations are performed up to the second order in the expansion in powers of the coupling constant. Renormalized effective Hamiltonians are well-defined symmetric forms acting in the Fock space as opposed to the renormalized bare Hamiltonian, which is not well defined without regularization. We introduce computational techniques that should streamline higher-order calculations and may be of independent interest.

Ab initio calculations↗
DOE OSTI2024

Developing Platinum-Group-Metal-Free Catalysts for Oxygen Reduction Reaction in Acid: Beyond the Single Metal Site

This project is to develop M (x) -N-C catalysts with dense multiple metal center (MMC) sites to meet the DOE 2025 activity target of 0.044 mA/cm 2 at 0.9 V IR-free ., as well as other goals such as durability. We have made important contributions to both catalyst development and fundamental understandings of the active sites in M-N-C catalysts in this project. We successfully made M (x) -N-C catalysts with some multiple metal center (MMC) sites by combining ionothermal carbonization with chemical vapor deposition (CVD). These catalysts, however, are not as active as the most active single-atom Fe-N-C catalysts made by the similar CVD process, owing likely to the low site density and the presence of inorganic Fe species such as iron carbides and nanoparticles. The most significant accomplishments we achieved in this project are: (1) we unraveled the formation pathway of Fe-N 4 sites during pyrolysis step-by-step and identified the trans-metalation mechanism, in collaboration with Deborah Myers and her colleagues at Argonne National Laboratory (ANL); (2) inspired by this finding, we pioneered the CVD synthesis of M-N-C catalysts (M = Mn, Fe, and Co), in which the Fe-N-C catalyst by CVD demonstrated an ORR activity of 33 mA/cm 2 at 0.9 V in H 2 -O 2 proton exchange membrane fuel cells (PEMFCs), very close to the ultimate goal of 35 mA/cm 2 of our project. This catalyst is the first Fe-N-C catalyst that contains only D1 sites without the D2 sites; whereas D1 and D2 sites have been always identified by Mossbauer in previous Fe-N-C catalysts. This finding helps to understand what the D1 and D2 sites are and their roles in catalyzing the ORR. (3) by improving the mass transport of the carbon matrix prior to the CVD process, the revised Fe-N-C catalyst made by CVD delivered a maximum power density of 0.53 W/cm 2 in H 2 -air PEMFCs. The improvement strategy was partly inspired by the computational modeling work by Adam Weber from LBNL, the Co-PI of this project, by developing, coding, and exercising a continuum level model of transport phenomena within a PGM-free catalyst layer. The model demonstrated that local resistances combined with limited site density of the PGM-free catalyst can result in limiting currents and poor polarization performance. The model also gave design guidance for impact of ECSA and overall catalyst-layer thickness. However, both Fe-N-C and Co-N-C catalysts developed by CVD showed poor durability in PEMFCS, in comparison with the traditional M-N-C catalysts synthesized via regular pyrolysis process. Consequently, we did not achieve the proposed durability targets. Despite so, we believe the FeNC-CVD catalysts with the poor durability and D1 sites only provides an excellent platform to understand the degradation mechanism of Fe-N-C in PEMFCs, the most important challenge in the development of M-N-C catalysts.

08 HYDROGEN↗
DOE OSTI2024

Forced Oscillation Source Location of Bulk Power Systems Using Synchrosqueezing Wavelet Transform

Forced oscillation source location (FOSL) plays a significant role in mitigating forced oscillations (FOs), which threaten the power system stability. Here, this paper proposes a data-driven approach for FOSL in power systems using synchrosqueezing wavelet transform (SWT). The proposed approach conducts SWT on the measured system responses to obtain the SWT matrix. Then, the SWT-based dissipating energy flow (DEF) model in time-frequency domain and dissipating energy spectrum (DES) model in frequency domain are derived from the traditional DEF model. Further, the characteristics of SWT-based DEF and DES are revealed by referring to the traditional DEF, and the FOSL criteria of the SWT-based DEF and DES can be hereby obtained. Using the obtained FOSL criteria, the FO source can be located from the measured responses. The performance of the proposed FOSL method is evaluated by simulation data of the WECC 179-bus test system and field-measurement PMU data of the ISO New England. The results confirm the accuracy and efficiency of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Performing Bayesian Analyses With AZURE2 Using BRICK: An Application to the 7 Be System

Phenomenological R-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical challenge has always been the estimation of the uncertainty on both the cross sections at the energies of interest and the fit parameters, which can take the form of standard level parameters. Frequentist (χ 2 -based) estimation has been the norm. In this work, a Markov Chain Monte Carlo sampler, emcee, has been implemented for the R-matrix code AZURE2, creating the Bayesian R-matrix Inference Code Kit (BRICK). Bayesian uncertainty estimation has then been carried out for a simultaneous R-matrix fit of the 3 He (α,γ) 7 Be and 3 He (α,α) 3 He reactions in order to gain further insight into the fitting of capture and scattering data. Both data sets constrain the values of the bound state α-particle asymptotic normalization coefficients in 7 Be. The analysis highlights the need for low-energy scattering data with well-documented uncertainty information and shows how misleading results can be obtained in its absence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

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

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Forecasting for the Weather Driven Energy System - A New Task under IEA Wind

The energy system needs a range of forecast types for its operation in addition to the narrow wind power forecast that has been the focus of considerable recent attention. Therefore, the group behind the former IEA Wind Task 36 Forecasting for Wind Energy has initiated a new IEA Wind Task with a much broader perspective, which includes prospective interaction with other IEA Technology Collaboration Programmes such as the ones for PV, hydropower, system integration, hydrogen etc. In the new IEA Wind Task 51 (entitled "Foreacsting for the Weather Drive Energy System") the existing Work Packages (WPs) are complemented by work streams in a matrix structure. The Task is divided in three WPs according to the stakeholders: WP1 is mainly aimed at meteorologists, providing the weather forecast basis for the power forecasts. In WP2, the forecast service vendors are the main stakeholders, while the end users populate WP3. The new Task 51 started in January 2022. Planned activities include 4 workshops. The first will focus on the state of the art in forecasting for the energy system plus related research issues and be held during September 2022 in Dublin. The other three workshops will be held later during the 4-year Task period and address (1) seasonal forecasting with emphasis on Dunkelflaute, storage and hydro, (2) minute-scale forecasting, and (3) extreme power system events. The issues and conclusions of each of the workshops will be documented by a published paper. Additionally, the Recommended Practice on Forecast Solution Selection will be updated to reflect the broader perspective.

geophysics computing↗
DOE OSTI2020

NLOX, a one-loop provider for Standard Model processes

NLOX is a computer program for calculations in high-energy particle physics. It provides fully renormalized scattering matrix elements in the Standard Model of particle physics, up to one-loop accuracy for all possible coupling-power combinations in the strong and electroweak couplings, and for processes with up to six external particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Gauge invariance of radiative jet functions in the position-space formulation of SCET

In subleading powers of soft-collinear effective theory (SCET), the Lagrangian contains couplings between soft quarks and hard-collinear quarks. Matrix elements of the hard-collinear parts of these couplings are radiative jet functions. In the position-space formulation of SCET, the Lagrangians are constructed from operators that appear to be gauge invariant. Nevertheless, we find violations of gauge invariance arise in the hard-collinear sector because gauge transformations can shift the momentum of a hard-collinear quark field from the hard-collinear sector to the soft sector, where the hard-collinear fields, by definition, have no support. The violations of gauge invariance are manifested in perturbation theory in the hard-collinear sector through the absence of certain Feynman diagrams that would be present in full QCD. A consequence of the absence of these diagrams is that the radiative jet functions that follow directly from the position-space Lagrangians are not gauge invariant, and we demonstrate this through explicit calculations in lower-order perturbation theory. We obtain gauge-invariant Lagrangians by adding to existing position-space Lagrangians terms that are proportional to the soft-quark equation of motion. These gauge-invariant Lagrangians are valid for nonzero, as well as zero, quark masses. We also remark briefly on the gauge invariance of certain Lagrangians that have been constructed in the label-momentum formulation of SCET. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Simulation of enhanced light extraction from periodic, disordered, and quasi-periodic OLED structures

Here we have developed a rigorous scattering matrix theory of light emission from periodically structured media using a Green’s function approach. We computationally simulate the spectral power inside the structure, incorporating Purcell factor enhancements, and find internal waveguiding modes and plasmonic losses. We simulate the light-outcoupling factor η out and describe how corrugations and structured media can enhance η out . We have extended our framework to describe non-periodic disordered arrays as well as aperiodic quasi-crystalline arrays. Flat organic light-emitting diodes (OLEDs) have low outcoupling with η out ~20% since most of the light is trapped in waveguided modes in higher index layers or lost to plasmonic excitations at the metal cathode. Periodically corrugated OLEDs can achieve highly enhanced η out ~65%-70% for pitch values between 1000 and 2000 nm, representing an enhancement factor >3 over planar structures. Periodic corrugations strongly diffract trapped waveguided and plasmonic modes to the emissive air cone. Disordered templates also can cause significant enhanced outcoupling with η out ~50%-55% for smaller nearest-neighbor separations of 300–400 nm. Quasi-crystalline tilings can also lead to enhancements of η out ~50%-55%. This framework can be utilized to design novel structured media that can generate high light extraction.

42 ENGINEERING↗
DOE OSTI2021

Fabrication of MiniFuel Compacts for High-Power Irradiation Testing of TRISO Fuel

MiniFuel compacts containing 20 individual tristructural isotropic (TRISO)-coated fuel particles were characterized to provide supporting preirradiation characterization data. The preirradiated characterization data supports irradiation capsule design and safety analysis, as well as downstream postirradiation examination (PIE) for a planned MiniFuel irradiation to explore high particle powers. The analysis included dimensional inspection, x-ray radiography and tomography, as-fabricated defect fraction analysis, and a matrix impurity analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Online Data-Enabled Predictive Control

We develop an online data-enabled predictive (ODeePC) control method for trajectory tracking of unknown systems, building upon the recently proposed DeePC. Our proposed ODeePC method leverages a primal-dual algorithm with real-time measurement feedback to iteratively compute the corresponding real-time optimal control policy as system conditions change. Specifically, our developed ODeePC: a) records data from the unknown system and updates the underlying primal-dual algorithm dynamically, b) can track changes in the system's operating point and adjust the control inputs, and c) is computationally efficient as it deploys a Fast Fourier Transform-based algorithm enabling the fast computation of the product of a non-square Hankel matrix with a vector. We provide theoretical guarantees regarding the asymptotic behavior of ODeePC and demonstrate its performance through a power system application.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗