Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear systems”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions. Today, tools such as transfer matrix and finite elements enable the calculation of the rotor dynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems. In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered. The purpose of this paper is to show examples of how the linear vs. nonlinear formulations differ. This paper will also show how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

42 ENGINEERING↗
DOE OSTI2025

Impact of position and density of nanoscale voids on fracture initiation in iron from phase field fracture simulation

Understanding the impact of these bubbles on crack propagation, like that of helium bubble-induced cracking in irradiated materials is incredibly complex. A useful first study towards understanding bubble effects on fracture is to examine how voids impact fracture first. In this work, we used phase-field fracture simulations to examine the influence of voids and their distribution on Mode I fracture in Fe. Assuming brittle fracture, two simulation configurations were considered: (1) nanoscale systems with one or two voids, and (2) nanoscale systems with an experimentally relevant distribution of voids, with up to 20% void area. Results from simulations with one and two voids showed that voids within 10 nm of a crack tip reduce the stress required for crack growth, with the magnitude of reduction depending on void-to-crack orientation. Comparisons with linear elastic fracture mechanics and evaluation of one versus two void systems revealed deviations from linear superposition, implying complex interactions between void and crack tip stress fields. In multi-void simulations, as void sizes increase, the nearest void to the crack tip exerts a greater influence on fracture stress than the overall porosity. Furthermore this study provides valuable insights into the relationship between void size and concentration, and the stress necessary for crack growth, marking a step forward towards understanding He bubble-induced fracture in ferrous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Load Mitigating Feedback Control of WECs

This report outlines the development of load-mitigating feedback control for wave energy converters. A simple, self-tuning multi-objective controller is demonstrated in simulation for a 3-DOF (surge, heave, pitch) point absorber. In previous work, the proposed control architecture has been shown to be effective in experiment for a variety of device archetypes for the single objective of the maximization of electrical power capture: here this architecture is extended to reduce device loading as well. In particular, PTO actuation forces and the minimization of fatigue damage (determined from the sum of wave-exerted and PTO forces) are considered as additional objectives for the self-tuning controller. This controller is demonstrated for two similar, but distinct systems: one described by the identified linear models from physical testing of the WaveBot device, and another based upon a WEC-Sim simulation that expands upon boundary element method data from the WaveBot device. In both cases, because the power surface is consistently fairly flat in the vicinity of control parameters that maximize power capture in contrasting sea-states, it is found to be generally possible to mitigate either fatigue damage or PTO load. However, PTO load is found to conflict with fatigue damage in some sea-states, limiting the efficacy of control objectives that attempt to mitigate both simultaneously. Additionally, coupling between the surge and pitch DOFs also limits the extent to which fatigue damage can be mitigated for both DOFs in some sea-states. Because control objectives can be considered a function of the sea-state (e.g., load mitigation may not be a concern until the sea is sufficiently large) a simple transition strategy is proposed and demonstrated. This transition strategy is found to be effective with some caveats: firstly, it cannot circumvent the aforementioned objective contradictions. Secondly, this objective transition is too slow to act as a system constraint, and objective thresholds must thus be considered quite conservatively. Improvement of the adjustment strategy is demonstrated through the addition of an integral term. Selection of well-performing transition parameters can be a function of sea-state. While a simple selection procedure is proposed, it is non-optimal, and a more robust selection procedure is suggested for future work.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Reconstruction of dynamic scenes based on differences between collected view and synthesized view

A system for generating a 4D representation of a scene in motion given a sinogram collected from the scene while in motion. The system generates, based on scene parameters, an initial 3D representation of the scene indicating linear attenuation coefficients (LACs) of voxels of the scene. The system generates, based on motion parameters, a 4D motion field indicating motion of the scene. The system generates, based on the initial 3D representation and the 4D motion field, a 4D representation of the scene that is a sequence of 3D representations having LACs. The system generates a synthesized sinogram of the scene from the generated 4D representation. The system adjusts the scene parameters and the motion parameters based on differences between the collected sinogram and the synthesized sinogram. The processing is repeated until the differences satisfy a termination criterion.

Kim, Hyojin↗
DOE OSTI2024

Evaluating Heat Extraction Performance of Closed-Loop Geothermal Systems with Thermally Conductive Enhancements in Conduction-Only Reservoirs: Preprint

We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.

advanced geothermal systems↗
DOE OSTI2020

2020 International Conference on Smart Grids and Energy Systems (SGES)

This paper is concerned about phasor-based linear state estimation of three-phase power grids which contain non transposed lines and unsymmetrical branches. It builds upon our previous work where the linear measurement equations were decoupled into three independent sets via the symmetrical component transformation. Such a transformation may however not fully decouple the measurement equations if there are nontransposed transmission lines and unsymmetrical branches in the system. This paper addresses this limitation by deriving additive terms that compensate for the remaining coupling between sequence measurements after the application of the conventional transformation. It will be shown that general unbalanced threephase power grids can still be transformed into three decoupled subsystems and their state estimation solution can be obtained in parallel by the proposed approach.

Khalili, Ramtin↗
DOE OSTI2023

Agriculture, bioenergy, and water implications of constrained cereal trade and climate change impacts

International trade increases connections and dependencies between countries, weaving a network of global supply chains. Agricultural commodity trade has implications for crop producers, consumers, crop prices, water and land uses, and other human systems. Interconnections among these systems are not always easy to observe when external impacts penetrate across multiple sectors. To better understand the interactions of non-linear and globally coupled agricultural-bioenergy-water systems under the broader economy, we introduce systematic perturbations in two dimensions, one human (restrictions on agricultural trade) and the other physical (climate impacts on crop yields). We explore these independently and in combination to distinguish the consequences of individual perturbation and interactive effects in long-term projections. We show that most regions experience larger changes in cereal consumption due to cereal import dependency constraints than due to the impacts of climate change on agricultural yields. In the scenario where all regions ensure an import dependency ratio of zero, the global trade of cereals decreases ~50% in 2050 compared to the baseline, with smaller decreases in cereal production and consumption (4%). The changes in trade also impact water and bioenergy: global irrigation water consumption increases 3% and corn ethanol production decreases 7% in 2050. Climate change results in rising domestic prices and declining consumption of cereal crops in general, while the import dependency constraint exacerbates the situation in regions which import more cereals in the baseline. The individual and interactive effects of trade perturbations and climate change vary greatly across regions, which are also affected by the regional ability to increase agricultural production through intensification or extensification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

LCLS-II Helium Refrigeration System Commissioning Results

SLAC National Accelerator Laboratory has upgraded to LCLS-II, featuring a 4 GeV superconducting linear accelerator composed of 37 cryomodules and two large helium refrigeration systems with a cooling capacity of 4 kW at 2.0 K. The LCLS-II Helium Refrigeration System (HRS) consists of two compressor stations, each with a power of approximately 4.5 MW, two 4.5 K cold boxes, each with a power of 18 kW equivalent at 4.5 K, and two sets of cold compressors that can each produce a flow of 230 g/s at 31 mbar (2.0K). Performance tests of the HRS were meticulously planned and successfully carried out, with results demonstrating that it exceeded the process requirements for LCLS-II operations. This paper provides a detailed presentation of the LCLS-II HRS performance and the challenges encountered during the commissioning phase.

42 ENGINEERING↗
DOE OSTI2025

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Enhancing Resilience of Urban Systems Against Climate-Induced Floods Using Advanced Data-Driven and Computing Techniques: A Driver-Pressure-State-Impact-Response (DPSIR) Framework

Focal Area: This white paper proposes advanced data-driven and computing techniques to improve the predictability of extreme climate events (precipitation) and focuses on comprehending their impact dynamics to enhance resilience of urban systems. Science Challenge: In urban systems, climate change has presented significant challenges related to accurate spatiotemporal prediction of extreme climate events (heavy precipitation) and their impacts, i.e. flooding. The occurrence of extreme climate events and their compounding repercussions on the built environment make urban systems increasingly more vulnerable. The interactions between extreme climate events, urban systems, and the corresponding impacts are irregular, non-linear, and dynamic, making it difficult to standardize resilience. As it stands, the management of complex urban systems is affected by uncertainties related to the nature of short- and long-term consequences caused by extreme events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

pvplr-python: Python package implementation of PVplr for Performance Loss Rate (PLR) analysis

Due to software fragmentation, PV system modeling teams can be limited to language specific packages, preventing cross-sectional analysis of different modeling techniques and workflows. To this end, PVplr, a popular PV performance modeling R software package, has been ported to the Python programming language. To verify and test the robustness of the port, NSRDB data has been used to simulated PV installations at native resolution (~2 million Sites), with a variety of degradation rates, degradation patterns, and modules. Performance Ratios were calculated using the ported functions from pvplr-python and compared against Rdtools YoY values. Due to the complicated nature of degradation, a new metric has been proposed to quantify the performance loss of a system. The cumulative production loss, is the total amount of energy lost due to the degrading performance of the system. Cumulative production loss alleviates the problems with fitting linear functions to non-linear degradation. Cumulative Production loss was shown to better estimate the total lost revenue for non-linear degradation patterns. $XbX + UTC$ was found to most accurately predict the total lost revenue in simulated systems.

Kumar, Suraj↗
DOE OSTI2024

PIP-II Beam Current Monitor Fault Case & Beam Position Monitor Linearity Studies in CST Studio

The use of non-invasive sensors & systems to measure particle beam characteristics is a crucial part of modern accelerator control systems. As such, simulations predicting the behavior of these sensors are essential for guiding beamline control system design. This poster details the results returned by two beam sensor models created using CST Studio software: the signal linearity of an elliptical beam position monitor (BPM), and the fault cases of an AC current transformer (ACCT) beam current monitor (BCM).

Rouzky, Adam↗
DOE OSTI2021

Measurement placement in electric power transmission and distribution grids: Review of concepts, methods, and research needs

Sensing and measurement systems are quintessential to the safe and reliable operation of electric power grids. Their strategic placement is of ultimate importance because it is not economically viable to install measurement systems on every node and branch of a power grid, though they need to be monitored. An overwhelming number of strategies have been developed to meet oftentimes multiple conflicting objectives. The prime challenge in formulating the problem lies in developing a heuristic or an optimisation model that, though mathematically tractable and constrained in cost, leads to trustworthy technical solutions. Further, large-scale, long-term deployments pose additional challenges because the boundary conditions change as technologies evolve. For instance, the advent of new technologies in sensing and measurement, as well as in communications and networking, might impact the cost and performance of available solutions and shift initially set conditions. Also, the placement strategies developed for transmission grids might not be suitable for distribution grids, and vice versa, because of unique characteristics; therefore, the strategies need to be flexible, to a certain extent, because no two power grids are alike. Despite the extensive literature on the present topic, the focus of published works tends to be on a specific subject, such as the optimal placement of measurements to ensure observability in transmission grids. There is a dearth of work providing a comprehensive picture for developing optimal placement strategies. Because of the ongoing efforts on the modernisation of electric power grids, there is a need to consolidate the status quo while exposing its limitations to inform policymakers, industry stakeholders, and researchers on the research-and-development needs to push the boundaries for innovation. Accordingly, this paper first reviews the state-of-the-art considering both transmission and distribution grids. Then, it consolidates the key factors to be considered in the problem formulation. Finally, it provides a set of perspectives on the measurement placement problem, and it concludes with future research directions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2022

Effect of ion interactions on the Raman spectrum of NO 3 − : Toward monitoring of low-activity nuclear waste at Hanford

Raman spectroscopy is a valuable in-situ technique for many applications. The concentrations of species in complex ionic mixtures of nuclear waste can be estimated using Raman measurements. However, it has been experimentally observed that ion interactions can cause a modification of Raman peak intensities and positions. The present work explores the nonlinear behavior associated with the nitrate anion, NO 3 − , which is present in abundance in low activity nuclear waste. We examine changes in the main Raman peak of the nitrate anion in the presence of other ions. A wide range of concentrations are covered, including those expected during direct-feed low-activity waste (DFLAW) processing at the Hanford site in the State of Washington. The experiments showed that the ions interact and associate to form ion pairs, which results in a blue shift (i.e., a shift towards higher wavenumbers) in the main Raman peak of nitrate. The results indicate that cation concentration is a better predictor of the peak shift, compared to ionic strength, both for binary and multicomponent mixtures. These findings have direct implications on the development of spectra-to-composition models for the DFLAW system, since they show deviations from the linearity assumptions used in common chemometric models.

42 ENGINEERING↗
DOE OSTI2022

Optimal control to handle variations in moisture content and reactor in-feed rate

The variations in feedstock characteristics, such as moisture and particle size distribution, lead to an inconsistent flow of feedstock from the biomass pre-processing system to the reactor in-feed system. These inconsistencies result in low on-stream times at the reactor in-feed equipment. This research develops an optimal process control method for a biomass pre-processing system comprised of milling and densification operations to provide the consistent flow of feedstock to a reactor's throat. This method uses a mixed-integer optimization model to identify optimal bale sequencing, equipment in-feed rate, and buffer location and size in the biomass pre-processing system. This method, referred to as the hybrid process control (HPC), aims to maximize throughput over time. We compare HPC with a baseline feed forward process control. Our case study based on switchgrass finds that HPC reduces the variation of a reactor's feeding rate by up to 100% without increasing the operating cost of the biomass pre-processing system for biomass with moisture ranging from 10 to 25%. Additionally, HPC reduces the cost of processing biomass by 0.36%–2.22%, and reduces processing time by 0.35%-2.24%. Furthermore, a biorefinery can adapt HPC to achieve its design capacity.

09 BIOMASS FUELS↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2022

Anatomy of Protein Electrospray Ionization Mass Spectra by Superconducting Tunnel Junction Mass and Energy Spectrometry

Cryogenic superconducting tunnel junction (STJ) detectors have the advantage of single-particle sensitivity, high quantum efficiency, low noise, and the ability to detect the time and relative impact energy of deposited ions. This makes them attractive for use in mass spectrometry (MS) and as a form of energy spectrometry. STJ cryodetectors have been coupled to time-offlight (TOF) mass spectrometers equipped with a matrix-assisted laser desorption ionization (MALDI) source and to an electrospray ionization (ESI) TOF mass spectrometer. Here, a lab-made linear ion trap (LIT) mass spectrometer system was coupled to an ESI source and a 16-channel Nb-STJ array with improved readout electronics. The goal was to investigate fundamentals of ESI-generated protein ions by further exploiting the advantage of resolving these ions in a third dimension of the relative energy deposited into the STJs. The proteins equine cytochrome c, bovine carbonic anhydrase, bovine serum albumin, and murine immunoglobulin G were studied using this ESI-LIT-STJ-MS instrument. Multiply charged monomers, multimers, and fragments from metastable ions were resolved from monomer peaks by differences in ion deposition energy even when these ions have the same mass-to-charge ratio as the corresponding monomer. The determination of fragment mass from metastable decomposition is accomplished without knowing the charge state of the fragment. The average charge state of the multimers is reduced with each addition of a protein which is a direct reflection of the surface area available for charging. Multiply charged in-source fragments have also been observed and distinguished in the mass spectrum of carbonic anhydrase by using the differences in the energy deposited in the STJs.

42 ENGINEERING↗