SEARCH · Engineering Papers
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.
A method is presented for determining if the connector adaptor assembly (CAA) for a particular unit, as determined by Coordinate Measurement Machine (CMM), is suitable for use in an Alternate Configuration (AC).
Abstract not provided.
Journal article based on previous SDRD work is to be submitted to Springer Nature
Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.
Abstract We develop a new method for analytical inversion of binned exoplanet transit spectra and for retrieval of planet parameters. The method has a geometrical interpretation and treats each observed spectrum as a single vector r → in the multidimensional spectral space of observed bin values. We decompose the observed r → into two orthogonal components: a wavelength-independent component r → ∥ corresponding to the spectral mean across all observed bins, and a transverse component r → ⊥ that is wavelength dependent and contains the relevant information about the atmospheric chemistry. The method allows us to extract, without any prior assumptions or additional information, the relative mass (or volume) mixing ratios of the absorbers in the atmosphere, the scale height to stellar radius ratio, H / R S , and the atmospheric temperature. The method is illustrated and validated with several examples of increasing complexity.
The Savannah River National Laboratory Nuclear Measurements Group was tasked by the Mark 18A Program team with developing three analytical characterization methods. These methods will support process and waste characterization of the Mark 18A program material. While methods already existed for most analytes of interest to the Mark 18A program team, the Nuclear Measurements Group developed and refined methods to quantify the Cf isotopes, 107 Pd, and 121m Sn. The results of that effort are presented in this work. Methods were successfully developed to characterize the Mark 18A samples for the needed isotopes and the NMG is ready to receive Mark 18A program samples.
The analysis of particulates from environmental sampling is routinely performed for nuclear forensics applications. In order to increase the tools available for nuclear forensics applications, capability development materials (CDMs), or reference particulates, are necessary for developing and benchmarking new analytical methods. Historically, these CDM particulates have been produced with highly controlled and characterized isotopic and size parameters for applications in developing and benchmarking particle sizers and mass spectrometry analytical methods [1- 4]. However, the development of CDMs with highly characterized particle morphology and phase of the particles is vital for aiding in the development and benchmarking of particle analytical methods for those parameters. In many cases, key properties such as crystallographic phase, morphology, and microstructure, may be correlated to processing history of environmental sampling particulates [5]. To that end, this work investigates determining the structural fingerprint of produced Pu-doped uranium oxide CDMs using transmission electron microscopy (TEM). The particulates, one of which shown in Figure 1, were synthesized to Fig. 1. Single particulate of uranium oxide fabricated using the THESEUS technique. develop capabilities for environmental sampling investigations. These particles are monodisperse at diameter of 1 μm and a range of plutonium concentrations from 0- 1,000 ppm. Given their small diameter, TEM analysis was ideal for studying the particulate structural fingerprint in detail. Previous investigations confirmed the external homogeneity of the particulates at different plutonium concentrations, however this analysis focuses on determining the grain structure, phase distribution, and porosity of the particulates.
Network macroscopic fundamental diagrams (MFDs) have recently been shown to exist in real-world urban traffic networks. When present, MFDs can be used to model traffic dynamics within an urban network by dividing the network into a set of spatially compact homogeneous regions and tracking the average level of congestion in each region. Existing analytical methods to estimate MFD mostly focus on the behavior of a single type of vehicle and do not capture the patterns of mixed traffic (e.g., cars and buses). The existence of buses matters since a bus will block the movements of other vehicles when it dwells at the bus stop. This paper proposes an analytical method to estimate the impact of bus dwelling on a network’s MFD based on the network’s geometric features, traffic control strategies, and bus operation parameters, and validates the performance of the proposed method using simulations based on microscopic traffic models. Comparisons of the analytical and simulation results show that the proposed analytical method can generally provide a good estimate of the lower bound and upper bound of the network’s MFD.
With continuing advances in high-performance parallel computing platforms, parallel algorithms have become powerful tools for development of faster than real-time power system dynamic simulations. In particular, it has been demonstrated in recent years that parallel-in-time (Parareal) algorithms have the potential to achieve such an ambitious goal. Here, the selection of a fast and reasonably accurate coarse operator of the Parareal algorithm is crucial for its effective utilization and performance. This paper examines semi-analytical solution (SAS) methods as the coarse operators of the Parareal algorithm and explores performance of the SAS methods to the standard numerical time integration methods. Two promising time-power series-based SAS methods were considered; Adomian decomposition method and Homotopy analysis method with a windowing approach for improving the convergence. Numerical performance case studies on 10-generator 39-bus system and 327-generator 2383-bus system were performed for these coarse operators over different disturbances, evaluating the number of Parareal iterations, computational time, and stability of convergence. All the coarse operators tested with different scenarios have converged to the same corresponding true solution (if they are convergent) and the SAS methods provide comparable computational speed, while having more stable convergence to the true solution in many cases.
This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.
COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.
Creep crack growth is a phenomenon which arises in damaged metallic structures under combined primary and secondary loads in the creep regime. The High Temperature Flaw Evaluation Code Committee of the American Society of Mechanical Engineers Boiler and Pressure Vessel Code (BPVC) is evaluating methods in extending Code Case N-934 to capture transient creep crack growth. Here, this paper provides context and basic examples on the selected approach to transient creep crack growth methods based on the C(t)-integral. The basis for selected C(t)-integral solution as analytical method is established, including an overview of its derivation. Practical assessments of a crack growing under creep conditions in a realistic component are conducted to illustrate the analytical approach. Additional considerations in the application of the analytical methods and limitations are discussed.
Magnetic integration is a hot topic in power electronics that concerns the use of a transformer’s leakage and magnetizing inductances purposefully in isolated power electronic converters, thereby giving the opportunity to save the cost and footprint of any additional inductor. This is of prime interest, especially in CLLLC resonant converters which require up to three inductors. For a complete integration of these inductances, the concept of a variable inductance transformer (VIT) is introduced in this thesis. A VIT is an adaptive magnetic structure that facilitates an easy adjustment of both magnetizing and leakage inductances to meet their desired values. However, for a more promising design, an accurate estimation of these inductances is necessary. While the evaluation of magnetizing inductance is quite straightforward, the calculation of leakage inductance is rather convoluted, because the leakage inductance is influenced by both the winding layout and the operating frequency. In this thesis, three new semi-analytical methods for calculating the frequency-independent leakage inductance, and a novel semi-analytical method for evaluating the frequency-dependent leakage inductance are proposed. These methods can calculate the respective leakage inductances of a VIT within an outstanding ±5% uncertainty. Finally, a bidirectional CLLLC resonant dc-dc converter is investigated for the constant current constant voltage (CCCV) charging of the next-generation 900 V traction battery of an electric vehicle. A new voltage gain equation is derived for designing the CLLLC resonant tank, and a small-signal model is presented for designing the variable-frequency feedback controller. Furthermore, a new methodology to design a VIT is developed to overcome the challenges associated with small coupling coefficients and guarantee a complete magnetic integration of the tank inductances. All theoretical results presented herein are verified through simulations and experiments performed on hardware prototypes designed in the lab.
An islanded microgrid can be used to serve critcal load as a resiliency source when a severe outage occurs. In an islanded mode, control of a microgrid relies on the communication system significantly. Hence, microgrids are cyberphysical systems and, therefore, the cyber system plays a crucial role in the performance of the cyber-power system. A microgrid control scheme is proposed for power dispatch and regulation based on the droop and proportional-integral (PI) feedback control. The proposed control strategy is validated for transient stability following dynamic events. To evaluate the impact of a networked control system on control performance, a cyber model is developed to represent data acquisition periods and communication delays. An analytical method is proposed to determine the critical values for the data reporting periods and communication delays. The analytical method based on a state space model of a networked control system is applicable for large-scale systems. A 2-dimensional stability region of a microgrid in the space of cyber parameters can be obtained by the proposed method, and the critical values of cyber parameters are determined based on the stability region. Simulation results validate that the design of a microgrid as a cyber-physical system needs to be guided by critical values for the data reporting period and communication delay to prevent system instability.
Electromagnetic Pulse (EMP) poses a significant threat to the normal operation of power systems, especially with the increasing penetration of renewable energy. Without adequate protection, EMP can severely damage equipment or result in circuit breakdowns or short circuits. Solar photovoltaic (PV) facilities are particularly susceptible to EMP since PV systems are outdoors and exposed to EMP radiation. To assess and mitigate this threat, this paper summarizes various models and tests used to study the effects of EMP on PV systems, assesses the nature of the threat, and identifies measures to mitigate it. The paper provides a comprehensive review by categorizing analytical methods based on their pros and cons, introducing and describing EMP-related events, and providing a taxonomy of analytical methods for studying EMP. This paper compares the processes of modeling, testing, and mitigating EMP at both the component and system levels of PV systems. Furthermore, it also presents a case study that reveals the vulnerability of solar inverters to EMP using the pulse current injection method. Finally, the paper discusses some of the remaining challenges that should be considered in future solar PV system design and deployment.
The purpose of this project is to develop analytical methods that can be deployed for analysis of lithium isotope enrichment products. One method involves rapid, high-throughput, relatively low precision analysis for direct monitoring of enrichment processes. A second method involves slower, but more precise analysis of final products of the enrichment process. We will use commercially available analytical equipment (quadrupole base inductively coupled plasma mass spectrometry and multi-collector magnetic sector inductively coupled plasma mass spectrometry) so that the analytical methods developed can be applied broadly.
The Savannah River Site stores approximately 36 million gallons of radioactive and hazardous waste that contains approximately 245 million curies. The waste is sent through various chemical processes to reduce its volume and to separate various components. The facility plans to replace formic acid (a chemical used to reduce soluble mercury) with glycolic acid. Recycle solution with glycolate may flow back to the tank farm, where the glycolate can generate hydrogen gas by thermal and radiolytic mechanisms. The current analytical method for detecting glycolate (ion chromatography) in supernatant requires a large dilution to reduce interference from the nitrate anions. Hydrogen nuclear magnetic resonance is an analytical method that requires less sample dilution. It takes advantage of the CH 2 group in glycolate. Liquid samples were spiked with four different levels of glycolate to build a calibration line, as it is recommended in the standard addition method. The detection and quantitation limits determined were 1 and 5 ppm, respectively, for 32 scans, which is well below the process limit of 10 ppm. In one test, 800 scans of a supernatant spiked with 1 ppm glycolate resulted in a -CH 2 peak with a signal-to-noise ratio of 36.
We develop an analytic method of inverting the Tolman–Oppenheimer–Volkoff relations to high accuracy. In principle, a specified energy density–pressure relation gives a unique mass–radius (M–R) relation and vice versa. Our method is developed from the strong correlations that are shown to exist between the neutron star mass–radius curve and the equation of state (EOS) or pressure–energy density relation. Selecting points that have masses equal to fixed fractions of the maximum mass, we find a semi-universal power-law relation between the central energy densities, pressures, sound speeds, chemical potentials, and number densities of those stars, with the maximum mass and the radii of one or more fractional maximum mass points. Rms fitting accuracies, for EOSs without large first-order phase transitions, are typically 0.5% for all quantities at all mass points. The method also works well, although less accurately, in reconstructing the EOS of hybrid stars with first-order phase transitions. These results permit, in effect, an analytic method of inverting an arbitrary M–R curve to yield its underlying EOS. We discuss applications of this inversion technique to the inference of the dense matter EOS from measurements of neutron star masses and radii as a possible alternative to traditional Bayesian approaches.