Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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.

318 records · Page 18

On the Stability of Power Transmission Systems Under Persistent Inverter Attacks: A Bi-Linear Matrix Approach

We investigate the stability and robustness properties of a power transmission system under persistent deceiving attacks on inverter-interfaced energy resources. The attacks can corrupt the damping coefficients in the inverters' controllers and measurements of the frequency at the points of coupling. Leveraging tools from hybrid dynamical systems theory, we characterize a broad family of persistent (and not necessarily periodic) attacks acting on the inverters, under which the stability properties of the transmission system can be shown to not be compromised. To address potentially conservative conditions identified through conventional bounding techniques, sufficient conditions on the average activation time of the attacks are identified via Lyapunov theory, as well as the formulation and solution of a class of bilinear matrix inequalities (BMI). The results are obtained for constant and slowly time-varying loads via input-to-state stability (ISS) tools. Numerical simulations on the IEEE 39-bus test system are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear analysis and vibro-impact characteristics of a shaft-bearing assembly

Here this study investigates the nonlinear frequency response of a shaft-bearing assembly with vibro-impacts occurring at the bearing clearances. The formation of nonlinear behavior as system parameters change is examined, along with the effects of asymmetries in the nominal, inherently symmetric system. The primary effect of increasing the forcing magnitude or decreasing the contact gap sizes is the formation of grazing-induced chaotic solution branches occurring over a wide frequency range near each system resonance. The system's nominal setup has very hard contact stiffness and shows no evidence of isolas or superharmonic resonances over the frequency ranges of interest. Moderate contact stiffnesses cause symmetry breaking and introduce superharmonic resonance branches of primary resonances. Even if some primary resonances are not present due to the system's inherent symmetry, their superharmonic resonances still manifest. Branches of quasiperiodic isolas (isolated resonance branches) are also discovered, along with a cloud of isolas near a high-frequency resonance. Parameter asymmetries are found to produce a few significant changes in behavior: asymmetric linear stiffness, contact stiffness, and gap size could affect the behavior of primary resonant frequencies and isolas.

42 ENGINEERING↗

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Analysis of the Effects of Renewable Energy Intermittency on the 2030 Korean Electricity Market

Republic of Korea has unique geographical characteristics similar to those of an island, resulting in an isolated power system. For this reason, securing sufficient operating reserves for the system’s stability and reliability in the face of the intermittency of increasing variable renewable energy (VRE) is paramount, and this will pave the way to achieving the nation’s decarbonization target and carbon neutrality. However, the current reserve-operation method in Republic of Korea does not take into account energy-system conditions, such as the intermittency of the VRE. Therefore, this paper presents an analysis of the impact of changes in reserve-operation methods on the electricity market in the future Republic of Korean power system, with the increased levels of VRE that are currently envisioned. Specifically, three reserve-operation methods, including Korea’s current reserve-power-operation standards, were applied to the two power-system plans announced by the Korean government to analyze the annual generator operation and costs. The analysis results show that securing reserves proportional to the VRE would exert negative effects, such as increased power-generation costs and the curtailment of nuclear and VRE generation. These results can contribute to the estimation of operational reserves needed for high levels of VRE and to the design of new the Korean reserve market, to be introduced in 2025.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Understanding Self-Assembly and the Stabilization of Liquid/Liquid Interfaces: The Importance of Ligand Tail Branching and Oil-Phase Solvation

Organophosphorus-based ligands represent a versatile set of solvent extraction reagents whose chemical makeup plays an important role in extraction mechanism. We hypothesize that the branching of the extractant hydrophobic tail and its oil-phase solvation affect the liquid/liquid interfacial structure. Understanding the structure mediated adsorption and interfacial ordering becomes key in designing ligands with enhanced selectivity and efficiency for targeted extractions. We employed vibrational sum frequency generation spectroscopy and interfacial tension measurements to extract thermodynamic adsorption energies, map interfacial ordering, and rationalize disparate behaviors of model di-(2-ethylhexyl) phosphoric acid and dioctyl phosphoric acid ligands at the hexadecane water interface. With increased surface loading, ligands with branched hydrophobic tails formed stable interfaces at much lower concentrations than those observed for ligands with linear alkyl tails. The lack of an oil phase and associated solvation results in markedly different interfacial properties, and thus measurements made at air/liquid surfaces cannot be assumed to correlate with the processes occurring at buried liquid/liquid interfaces. Here, we attribute these differences in the surface mediated self-assembly to key variations in hydrophobic interactions and tail solvation taking place in the oil phase demonstrating that interactions in both the polar and nonpolar phases are essential to understand self-assembly and function.

36 MATERIALS SCIENCE↗

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗

High-Temperature Linear Receiver Enabled by Multicomponent Aerogels

Concentrating solar thermal (CST) technology has significant potential mainly due to its dispatchability and low cost of storage. However, to compete with other sources, including utility-scale solar PV, its final cost (cents/kWh) still needs to be lowered. Cost reduction can be achieved by improving the system level efficiency of the CST plants through the deployment of advanced power cycles, which operate at high temperatures of ~700°C. However, optical, and thermal losses pose a major challenge to the efficiency of such CST systems. The overall aim of this project is to investigate and de-risk a linear solar receiver concept called an Aerogel Insulated Receiver (AIR) that generates high temperatures (up to 700°C) at a low solar concentration ratio (<100) and with a high collection efficiency (optical × receiver). Our prior work has demonstrated the thermal stability1 and optical and heat-insulating properties2,3 of transparent aerogel insulation at a one-inch scale. The focus of this work is on (1) co-optimization of the geometry of the aerogel tile and receiver enclosure to fit a standard parabolic collector (PTC), (2) scale-up of aerogels into 4-inch tiles while preserving key properties, (3) experimental measurement of receiver heat loss (W/m) in a >70-cm test stand and validation of anticipated receiver performance at high temperatures. Regarding (1), appropriate optical and thermal models for a parabolic trough receiver (PTR) are developed and validated. The geometry of the aerogels and the receiver enclosure are co-optimized to maximize the collection efficiency. The model predicts a 54% collection efficiency at 700°C for an AIR design based on flat aerogels. By combining the collection efficiency with the power block efficiency of supercritical CO 2 cycles, we predict >10% improvements in peak plant efficiency relative to existing line-focusing CST systems. The application of curved plasmonic aerogels is predicted to further increase the collection efficiency to 64%. Regarding (2), we demonstrate the successful development of 6-inch-long refractory aerogel tiles with optical, thermal, and stability characteristics consistent with our prior work. This scale-up requires a transition to a larger ALD station and modifying the ALD process variables such as exposure time and the number of precursor doses. Regarding (3), we design and develop an AIR test stand measuring 3 feet in length. Heat loss performance analysis is carried out using the test stand. The results indicate that aerogel insulation can significantly reduce receiver thermal losses at the high operating temperatures required for next-generation PTRs. The experimental results agree with the heat loss performance predicted by our receiver model. Lastly, we conducted preliminary failure mode and effects (FMEA) and techno-economic (TEA) analyses to identify failure mitigation strategies and commercial opportunities, respectively. Overall, this project identifies key opportunities and challenges in deploying aerogel insulating receivers in next-generation line-focusing CST technologies.

14 SOLAR ENERGY↗

Investigations into the Ternary NaF-KF-UF4 Salt System – Phase A

A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stable electrospray signal on a microfabricated glass chip with three-dimensional open edge and tiered depth geometries

This paper presents the microfabrication and performance of a three-dimensional electrospray ionization (ESI) emitter tip made from glass, which achieves stable current signals important for chemical analysis. Our novel microfabrication process and custom-built signal conditioning hardware provides the advantage of providing accurate features and steady signals. The fabrication process relies on standard microfabrication techniques (i.e., deposition, photolithography, and wet etching). This fabrication method involves the novel application of two layers of positive and negative photoresists in addition to Parafilm® wax tape. Open edge and tiered depth details were successfully created from a multilayer planar mask. This is a benefit for integrated miniaturized and microfluidic systems that often require micro features for their functionality but relatively large millimeter size features for their physical periphery. We demonstrate the fundamental performance of electrospray with this microfluidic chip. The emitter tip was fixed on a linear axis stage with high resolution (10 μm) to finely control the tip distance from a metal counter electrode plate. Further, a custom printed circuit board system was built to safely control four voltages applied to the microchip ports from a single high voltage power supply. To readily form the electrospray, non-aqueous solvents were used for their low viscosity and a constant voltage of +2.7 kV was applied to the sheath electrospray microchannel. The liquid being sprayed was 80/20 (v/v) methanol/acetonitrile with 0.1% acetic acid in the sheath microchannel and with ammonium acetate (10–40 mM) in its remaining microchannels. The electrospray signal was measured in response to varying the distance (1.4 to 1.6mm) between the electrospray emitter tip and a metal counter electrode plate in addition to the varying concentration of the background electrolyte, ammonium acetate. Stable and repeatable electrospray signal showed linear relationships with emitter tip distance and concentration (r 2 ≥ 0.95).

42 ENGINEERING↗

CMB-S4: Foreground-cleaning Pipeline Comparison for Measuring Primordial Gravitational Waves

We compare multiple foreground-cleaning pipelines for estimating the tensor-to-scalar ratio, r, using simulated maps of the planned CMB-S4 experiment within the context of the South Pole Deep Patch. To evaluate robustness, we analyze bias and uncertainty on r across various foreground suites using map-based simulations. The foreground-cleaning methods include: a parametric maximum likelihood approach applied to auto- and cross-power spectra between frequency maps; a map-based parametric maximum-likelihood method; and a harmonic-space internal linear combination using frequency maps. We summarize the conceptual basis of each method to highlight their similarities and differences. To better probe the impact of foreground residuals, we implement an iterative internal delensing step, leveraging a map-based pipeline to generate a lensing B-mode template from the large aperture telescope frequency maps. Our results show that the performance of the three approaches is comparable for simple and intermediate-complexity foregrounds, with σ(r) ranging from 3–5 ×10 −4 . However, biases at the 1σ–2σ level appear when analyzing more complex forms of foreground emission. By extending the baseline pipelines to marginalize over foreground residuals, we demonstrate that contamination can be reduced to within statistical uncertainties, albeit with a pipeline-dependent impact on σ(r), which translates to a detection significance between 2σ and 4σ for an input value of r = 0.003. These findings suggest varying levels of maturity among the tested pipelines, with the auto- and cross-spectra-based approach demonstrating the best stability and overall performance. Moreover, given the extremely low noise levels, mutual validation of independent foreground-cleaning pipelines is essential to ensure the robustness of any potential detection.

astronomy data analysis↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗