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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

On the convergence of an IEQ-based first-order semi-discrete scheme for the Beris-Edwards system

We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.

Weber, Franziska↗
DOE OSTI2021

Megavolt bremsstrahlung measurements from linear induction accelerators demonstrate possible use as a FLASH radiotherapy source to reduce acute toxicity

Abstract Recent studies indicate better efficacy and healthy tissue sparing with high dose-rate FLASH radiotherapy (FLASH-RT) cancer treatment. This technique delivers a prompt high radiation dose rather than fractional doses over time. While some suggest thresholds of > 40 Gy s −1 with a maximal effect at > 100 Gy s −1 , accumulated evidence shows that instantaneous dose-rate and irradiation time are critical. Mechanisms are still debated, but toxicity is minimized while inducing apoptosis in malignant tissue. Delivery technologies to date show that a capability gap exists with clinic scale, broad area, deep penetrating, high dose rate systems. Based on these trends, if FLASH-RT is adopted, it may become a dominant approach except in the least technologically advanced countries. The linear induction accelerator (LIA) developed for high instantaneous and high average dose-rate, species independent charged particle acceleration, has yet to be considered for this application. We review the status of LIA technology, explore the physics of bremsstrahlung-converter-target interactions and our work on stabilizing the electron beam. While the gradient of the LIA is low, we present our preliminary work to improve the gradient by an order of magnitude, presenting a point design for a multibeam FLASH-RT system using a single accelerator for application to conformal FLASH-RT.

42 ENGINEERING↗
DOE OSTI2024

A two-dimensional numerical study of ion-acoustic turbulence

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

An in-situ view cell system for investigating swelling behavior of elastomers upon high-pressure hydrogen exposure

The transition to hydrogen as a clean and efficient energy carrier is impeded by challenges in the compatibility of hydrogen with materials used within hydrogen infrastructure. Elastomers, crucial in sealing components, often exhibit premature failures in high-pressure hydrogen environments due to excessive swelling. This study employs an innovative in-situ view cell system to assess the swelling behavior of hydrogenated nitrile butadiene rubber (HNBR) under various hydrogen conditions. The system, designed to withstand pressures up to 96.5 MPa, incorporates Digital Image Correlation (DIC) for strain measurements and volume estimation. Results reveal non-linear volume increases during depressurization, challenging conventional assumptions. Furthermore, investigations into peak hydrogen pressures and pressure-holding scenarios during decompression highlight complex swelling trends. The introduction of a novel computer vision (CV) method enhances precision in volume estimation, overcoming DIC limitations. The study provides insights into mitigating elastomer swelling, crucial for developing robust materials to support future hydrogen-driven energy systems.

Elastomer↗
DOE OSTI2024

Universal Spreading of Conditional Mutual Information in Noisy Random Circuits

For this work, we study the evolution of conditional mutual information (CMI) in generic open quantum systems, focusing on one-dimensional random circuits with interspersed local noise. Unlike in noiseless circuits, where CMI spreads linearly while being bounded by the light cone, we find that noisy random circuits with an error rate 𝑝 exhibit superlinear propagation of CMI, which diverges far beyond the light cone at a critical circuit depth 𝑡 𝑐 ∝ 𝑝 −1 . We demonstrate that the underlying mechanism for such rapid spreading is the combined effect of local noise and a scrambling unitary, which selectively removes short-range correlations while preserving long-range correlations. To analytically capture the dynamics of CMI in noisy random circuits, we introduce a coarse-graining method, and we validate our theoretical results through numerical simulations. Furthermore, we identify a universal scaling law governing the spreading of CMI.

decoherence↗
DOE OSTI2025

Field Test Report Neutron Scintillator Array Dry Storage Cask Scanner FY2024

During two weeks of Field Testing at the Idaho National Laboratory INTEC Cask Farm in July and August 2024, the LLNL Dry Storage Cask Scanner Array was lifted on top of an MC-10 dry storage fuel cask and operated to acquire neutron and gamma-ray data from the 24 fuel bundle positions. Neutron and gamma-ray data acquisition scans across the top of the cask of varying dwell times were performed July 15-18, 2024 and August 19-22, 2024 to evaluate the ability of the scanner data to reveal asymmetries in the fuel positions that reflect asymmetries in the MC-10 cask fuel bundle loading. The MC-10 cask 24 position fuel bundle loading at the INTEC Cask Farm is well documented, including the locations of six empty fuel bundle positions. This loading presents an opportunity to test the ability of the scanner system to detect diversion of spent fuel bundles as well as to validate the MC-10 cask MCNP modeling. The cask scanner array consists of six Stilbene crystal scintillator detectors and a linear actuator frame that moves the six detectors across the MC-10 dry storage cask to obtain data above each of the 24 fuel bundle positions. The detectors are connected to a pulse-shape discrimination data acquisition system capable of generating separate neutron and gamma-ray spectra for each detector and for each scan position. From the prior single detector Field Test in 2021 and iteration with MCNP modeling, the neutron and gamma-ray data were analyzed in multiple energy regions to identify an analysis method that would provide the strongest and most consistent signature of the asymmetric MC-10 cask fuel loading1 . From both the 2021 Field Test and the current Field Test results, the neutron capture gamma-ray count rate around 2.2 MeV provides the strongest signature of the asymmetric MC-10 cask fuel loading and has qualitative agreement with MCNP calculations. Counting all gamma-rays produces a similar signature. Neutrons emerging from the cask top are moderated and captured by the hydrogen in the polyethylene moderator and scintillator detector, producing a 2.2 MeV gamma ray which is seen in the scintillator gamma-ray spectrum. The count rate in the 2.2 MeV gamma-ray region is ~50 c/s, which is ~1000x higher than the ~0.05 n/s rate in the > 4MeV neutron region, and ~50x greater than the ~1 n/s rate in the neutrons > 500 keV region. Analysis of the 2.2 MeV neutron-capture Compton-scattered gamma-rays produces a statistically significant signature of the INTEC Cask Farm MC-10 asymmetric fuel loading. MCNP simulations indicate that the average neutron energy spectrum offers the potential to detect a large asymmetry from several missing bundles as well as individual missing fuel bundles. Testing this feature will require measurements on a cask with single missing elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A Fast Flexibility-Driven Generation Portfolio Planning Method for Sustainable Power Systems

We are witnessing an acceleration in the uptake of renewable energy in power systems. Because of the associated variability and uncertainty of renewables, power systems need to have an adequate supply of flexibility to allow for suitable management of short-term operations. So far most of the work in this area has neglected how flexibility needs associated with renewables are fulfilled as part of dispatchable generation capital investments decisions. To address this challenge, we propose an approach to plan the dispatchable generation mix of a power system as needed to counteract variability and uncertainty associated with significant shares of variable renewable generation. The approach exploits the linear time-invariant feature of variable generation variability using historical phase planes of capacity (in MW) and ramp (in MW/min) to bridge the gap between long-term capacity planning and short-term intra-hour flexibility needs. This approach is much more computationally tractable than other proposals, while also being able to capture adequately short-term operational features like ramping and net load variability. Numerical tests are performed on realistic datasets to substantiate the effectiveness of the approach.

bulk power system planning↗
DOE OSTI2021

Renewable Thermal Energy Systems Designed for Industrial Process Solutions in Multiple Industries

The need for renewable heat in industry is vital for the next decade and beyond. Industrial decarbonization is a key area that must be accelerated, to foster the removal of fossil fuels from the provision of heat, especially at low temperatures. This paper looks at the development and results of case studies for understanding the economics and potential for renewable thermal energy systems (RTES), particularly in hybrid configurations to provide industrial process heat (IPH). For the case studies, these include non-concentrating e.g., heat pumps, and concentrating collectors e.g., parabolic trough collectors and direct steam generation (DSG)-linear Fresnel collectors (LFCs). The results show that a levelized cost of heat (LCOH) of $6-$8 per million British Thermal Units (MMBTU) is possible, depending on the direct normal irradiance (DNI) and the system sizing e.g., to increase the solar fractions. In Arizona for example, with a DNI of 7.36 kWh/m2/day, the base case for the DSG-LFC system with 6hrs of thermal energy storage could potentially meet a 1 MWth load 80% of the year.

concentrated solar thermal↗
DOE OSTI2020

Inspection Tools for Hanford Tanks and Waste Transport Systems - 20444

Recent integrity issues in the double-shell tanks (DST) at Hanford have motivated the need for developing innovative tools that can provide information regarding the health of the tanks. These issues include the primary linear failure of AY-102 and recent concerns of thinning in the DST secondary liners. Other concerns include erosion or corrosion on transfer lines and processing pipes. In recent years, Florida International University (FIU) has supported DOE-EM by developing tools that can assist in understanding the health of tanks and waste transport system at Hanford. More specifically, FIU has developed a miniature rover that can be deployed through small risers and gain access to refractory slots in the Hanford DSTs. It traverses through the slots on the primary liner upside down via magnets to avoid debris in the slots. Recent modifications include a spring loaded magnetic arm that allows the unit to traverse over large weld seams but does not reduce the overall functionality of the system. The system has also incorporated a control capsule that will sit within the annulus of the tank and improves on issues related to signal and voltage loss over long tether lines. Lastly, the unit is also being augmented with a second module that will include a means to prep the tank surface and obtain a thickness measurement using a dual element ultrasonic sensor. FIU has also developed a 6-inch peristaltic pipe crawler that operates similar to the previous pipe crawlers developed at FIU. This crawler is a marsupial type crawler that will navigate through the 6-inch drain lines, and deploy a small rover at the entrance of the drain slots. This rover will also traverse upside down on the secondary liner and provide information regarding the health of the liner. An initial prototype for the crawler unit has been developed, assembled and tested. The system is pneumatic and consists of three modules - two grippers and one extender. The validation tests are conducted in a mockup of the drain lines and demonstrate the system's ability to navigate through long pipe lengths and elbows. In this paper, a general overview of each of the aforementioned systems is presented along with results from testing and demonstrations that shows the applicability of each system. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods

Here, we present an energy-conserving discontinuous Galerkin scheme for the full-f electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for ∂A II /∂t, the inductive component of the parallel electric field, using a generalized Ohm’s law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms

The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that we recently proposed. In particular, this application relies on the near-identity MT (NIMT); however, the NIMT approximation used so far is not exactly unitary and leads to numerical instability. Here, we develop a discrete MT that is exactly unitary, and approximate it to obtain a discrete NIMT that is also unitary and can be computed in linear time. We prove that the discrete NIMT converges to the discrete MT when iterated, thereby allowing the NIMT to compute MTs that are not necessarily near-identity. Finally, we then demonstrate the new algorithms with a series of examples.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Floquet Spin Splitting and Spin Generation in Antiferromagnets

In antiferromagnetic spintronics, accessing the spin degrees of freedom is essential for generating spin currents and manipulating magnetic order, which generally requires lifting spin degeneracy. This is typically achieved through relativistic spin-orbit coupling or nonrelativistic spin splitting in altermagnets. Here, we propose an alternative approach: a dynamical spin splitting induced by an optical field in antiferromagnets. By coupling the driven system to a thermal bath, we demonstrate the emergence of steady-state pure spin currents as well as linear-response longitudinal and transverse spin currents. Crucially, thermal bath engineering enables a nonrelativistic Edelstein effect—the generation of a net spin accumulation—without relying on spin-orbit coupling. Furthermore, our results provide a broadly applicable and experimentally tunable route to control spins in antiferromagnets, offering new opportunities for spin generation and manipulation in antiferromagnetic spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Energy amplification in plasma crystals due to multiple torsions

Interacting torsions are examined within a two-dimensional monolayer crystal suspended in an argon complex plasma for 1–10 W discharge powers and pressures of 135–155 mTorr. Two torsions embedded in a lattice are shown to amplify the kinetic energy and range of motion of particles located between the torsions to nearly three times that observed in single torsion systems. It is also shown that multiple torsions can interact via amplified particle energy when separated by up to 14 interparticle distances (Δ). The torsion separation distance also showed a positive linear trend with power and a slightly positive correlation with the pressure. This amplification of energy is possible due to the fact that multiple torsions in a lattice increase the interparticle distance of the lattice by 16% more than single torsion systems, leading to additional freedom of motion in the lattice plane. These combined findings show that multiple torsions heat the lattice differently depending on their separation from the other torsion. The midpoint particles between torsions absorb the majority of energy from the two torsions, and energy addition at the midpoint is nonlinear. The addition of more torsions to the lattice may lead to melting of the plasma crystal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Upgrading Fermilab’s accelerator control system with ACORN

The Fermilab Accelerator Complex is the largest national user facility in the Office of High Energy Physics (DOE/HEP) program and the only national user facility operating at Fermilab. Fermilab serves as the host to the Long Baseline Neutrino Facility/Deep Underground Neutrino Experiment (LBNF/DUNE), the laboratory’s flagship project for neutrino science that is under construction. LBNF/DUNE will be powered by megawatt beams from an upgraded accelerator, the Proton Improvement Plan II (PIP-II) that will replace the laboratory’s aging linear accelerator with a new one based on superconducting radio-frequency cavities. The Accelerator Controls Operations Research Network (ACORN) Project will support LBNF/DUNE and PIP-II by modernizing the accelerator control system. The project is at the conceptual design phase and looking to achieve Critical Decision 1 (CD-1) later this year. The scope and structure of the project will be presented, along with an overview of how that has changed in the past year. Current design and technology choices will be shared. Specific challenges facing the project will be addressed, along with current thinking on solutions.

Roehrig, Christian [Fermilab]↗
DOE OSTI2020

Dynamic Optimal Dispatch of Energy Systems with Intermittent Renewables and Damage Model

With the increasing penetration of intermittent renewable energy sources into the grid, there is a growing need for process systems-based strategies that integrate dispatchable and variable energy systems for supplying the demand while maintaining grid reliability. The proposed framework corresponds to a dynamic mixed-integer linear programming optimization approach that integrates coal-fired and natural gas-fired power plants, NaS batteries for energy storage, and solar/wind energy to supply the demand. This optimization approach considers an economic goal and constraints to provide power balance while maintaining the overall damage of the natural gas combined cycle (NGCC) power plant drum under a maximum stress as well as avoiding the overheating of the NGCC superheater and reheater. Renewable curtailment levels are also retained at minimum levels. Case studies are analyzed considering different loads and renewable penetration levels. The results show that the demand was met for all cases. Grid flexibility was mostly provided by the NGCC, while the batteries were used sparingly. In addition, considering a CO 2 equivalent analysis, the environmental performance was intrinsically connected to grid flexibility and the level of renewable penetration. Stress analysis results reinforced the necessity for an equipment health-related constraint.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Managing Power Systems-Induced Wildfire Risks Using Optimal Scheduled Shutoffs

The increasing demands for electricity and the increase in extreme weather conditions are putting unprecedented pressure on our electric grids. Often, this pressure leads to electrical component failures, which might ignite wildfires. This work develops a novel model to balance the reliability of power network operations and the risk of wildfire ignition by opti- mizing the operational schedule of power transmission networks considering time-varying risk measures that include exogenous and operational factors. Energy storage systems are considered to deliver power during peak wildfire hours and enable temporal load shifting. The problem is formulated as a mixed-integer linear program that maximizes a weighted sum of the served power demand and the reduction in grid-induced wildfire risk. The results demonstrate the ability of the model to significantly reduce wildfire risk without considerable load shedding.

power systems operations↗