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

Impacts of year-to-year weather variability and inter-panel spacing on agrivoltaic crop yields in Massachusetts

Deployment of utility-scale solar power plants could lead to agricultural land-use changes. Agrivoltaics, a dual land use combining solar and agriculture on the same land, can provide multiple environmental benefits, including improving soil quality and water use efficiency. The body of agrivoltaic field data is still growing, and crop responses to different solar configurations under different local climates are highly varied. We investigate the impact of adding spacing between adjacent solar panels in a fixed-tilt system to improve light diffusion to crops. For four crops (broccoli, peppers, kale, Swiss chard) grown across 3 years in an agrivoltaic system in Massachusetts, we found that only kale had a linearly increasing trend as the inter-panel spacing increased from 0.6 m to 1.5 m (2 ft to 5 ft). However, there were significant year-to-year differences in the yield of agrivoltaic versus control fields. Agrivoltaic and full sun fields produced equivalent yields in a hot, dry year, whereas the full-sun control beds produced more salable yield for all four crops in a warm, wet year. This demonstrates variability of agricultural outcomes and the need for more multi-year studies to ensure agrivoltaic impacts are not under- or overestimated.

14 SOLAR ENERGY↗
DOE OSTI2021

A slip-spring framework to study relaxation dynamics of entangled wormlike micelles with kinetic Monte Carlo algorithm

Hypothesis: Wormlike micelles (WLMs) formed due to the self-assembly of amphiphiles in aqueous solution have similar viscoelastic properties as polymers. Owing to this similarity, in this work, it is postulated that kinetic Monte Carlo (kMC) sampling of slip-springs dynamics, which is able to model the rheology of polymers, can also be extended to capture the relaxation dynamics of WLMs. Theory: The proposed modeling framework considers the following relaxation mechanisms: reptation, union-scission, and constraint release. Specifically, each of these relaxation mechanisms is simulated as separate kMC events that capture the relaxation dynamics while considering the living nature of WLMs within the slip-spring framework. As a case study, the model is implemented to a system of sodium oleate and sodium chloride to predict the linear rheology and the characteristic relaxation times associated with the individual relaxation mechanisms at different pH and salt concentrations. Findings: Linear rheology predictions were found to be in good agreement with experimental data. Furthermore, the calculated relaxation times highlighted that reptation contributed to a continuous increase in viscosity while union-scission contributed to the decrease in viscosity of WLM solutions at a higher salinity and pH. Finally, this manifests the proposed model’s capability to provide insights into the key processes governing WLM’s rheology.

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

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Universal signal scaling in microwave impedance microscopy

Microwave impedance microscopy (MIM) is an emerging scanning probe technique that measures the local complex dielectric function using near-field microwave. Although it has made significant impacts in diverse fields, a systematic, quantitative understanding of the signal's dependence on various important design parameters is lacking. Here, we show that for a wide range of MIM implementations, given a complex tip-sample admittance change ΔΥ, the MIM signal—the amplified change in the reflected microwave amplitude—is –G · ΔΥ/2Υ 0 · η 2 · V in , where η is the ratio of the microwave voltage at the probe to the incident microwave amplitude, Y o is the system admittance, and G is the total voltage gain. For linear circuits, η is determined by the circuit design and does not depend on V in . We show that the maximum achievable signal for different designs scales with η 2 or η when limited by input power or sample perturbation, respectively. Furthermore, this universal scaling provides guidance on diverse design goals, including maximizing narrow-band signal for imaging and balancing bandwidth and signal strength for spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Toward scalable quantum computations of atomic nuclei

We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3 He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2 H and 3 He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Furthermore, based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Collective excitations and quantum incompressibility in electron-hole bilayers

In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Spin polarization engineering in 𝑑-wave altermagnets

Altermagnets host unconventional spin-polarized bands despite zero net magnetization, but controlling their spin structure remains challenging. Here, we theoretically propose a multifield approach to engineer spin polarization in 𝑑-wave altermagnets using gating, optical driving, and in-plane electric fields, which enable tunable and switchable polarizations along multiple directions. Optical driving induces out-of-plane (𝑧) polarization, while gating and in-plane fields generate 𝑥 and 𝑦 polarizations via the Edelstein effect, all of which are experimentally detectable. We further find that spin- and band-selective doping induces chiral optical activity, a feature unique to altermagnets. Our approach offers a versatile means to control the direction of spin polarization in altermagnets.

Floquet systems↗
DOE OSTI2024

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fast and Scalable Sparse Triangular Solver for Multi-GPU Based HPC Architectures

Designing efficient and scalable sparse linear algebra kernels on modern multi-GPU based HPC systems is a daunting task due to significant irregular memory references and workload imbalance across the GPUs. This is particularly the case for \textit{Sparse Triangular Solver (SpTRSV)} which introduces additional two-dimensional computation dependencies among subsequent computation steps. Dependency information is exchanged and shared among GPUs, thus warrant for efficient memory allocation, data partitioning, and workload distribution as well as fine-grained communication and synchronization support. In this work, we demonstrate that directly adopting unified memory can adversely affect the performance of SpTRSV on multi-GPU architectures, despite linking via fast interconnect like NVLinks and NVSwitches. Alternatively, we employ the latest NVSHMEM technology based on Partitioned Global Address Space programming model to enable efficient fine-grained communication and drastic synchronization overhead reduction. Furthermore, to handle workload imbalance, we propose a malleable task-pool execution model which can further enhance the utilization of GPUs. By applying these techniques, our experiments on the NVIDIA multi-GPU supernode V100-DGX-1 and DGX-2 systems demonstrate that our design can achieve on average 3.53x (up to 9.86x) speedup on a DGX-1 system and 3.66x (up to 9.64x) speedup on a DGX-2 system with 4-GPUs over the Unified-Memory design. The comprehensive sensitivity and scalability studies also show that the proposed zero-copy SpTRSV is able to fully utilize the computing and communication resources of the multi-GPU system.

Xie, Chenhao↗
DOE OSTI2021

District EnErgy Planner (DEEP) v2

The District EnErgy Planner - DEEP is an Mixed Integer Linear Programming (MILP) optimization software tool to solve district energy system, technology configuration, operation strategies, and energy retail business schemes. DEEP takes user input energy demand on heating, cooling and electricity demand. Energy loss through pipeline is calculated for heating and cooling from central plants. The optimal system configuration to minimize the district energy system life cycle cost or minimize CO2 emission is calculated based on energy balance equations and various constraints. The recent update on DEEP include a feature of waste heat utilization energy balance based optimization to take into accounts of different temperature bands of heat into account to maximize system efficiency. Another feature is to optimize district energy retail by provide optimal pricing scheme for heating, cooling and electricity to end users to incorporate district to grid demand response.

Feng, Wei↗
DOE OSTI2024

Modeling In-Reactor Temperature Gradients in Lithium Aluminate Ceramic Pellets

In reactor thermal performance of lithium aluminate ceramic pellets are calculated to determine operating temperatures if internal components shift during irradiation. A 2D conduction model was setup to model this scenario with additional physics, including thermal expansion of components to reactor temperatures, creep down of the cladding due to differential pressure between the inside of the rod and reactor coolant system, and net radiation exchange between components. The non-linear coupling of the governing equations presented convergence challenges due to the very different effects from a small temperature change on the radiation vs radial and circumferential conduction equations. Initial temperature guesses for the multidimensional solution were provided from conservative 1D models simulating bounding scenarios. Convergence was greatly improved by taking small steps towards the root solutions and averaging with initial guesses. The final method robustly converged all 132 points at an axial plane after a few dozen iterations. The temperature results show that the additional physics have a secondary effect as the dominant heat transfer mechanism is radial conduction across the gas gap followed by circumferential conduction around components. Circumferential conduction in the lithium aluminate pellets is enhanced by the relatively high thermal conductivity of this ceramic.

Carstens, Nathan A.↗
DOE OSTI2022

Strengthening the Security of Operational Technology: Understanding Contemporary Bill of Materials

The evolution of cyber-physical infrastructure has made its security more challenging. The last few years have witnessed a convergence of hardware and software segments in various domains, including operational technology (OT) which is responsible for carrying out critical tasks such as monitoring and controlling power grids, nuclear plants, transportation, and emergency services. Both hardware and software encapsulate numerous open source and proprietary subcomponents, making it crucial for end-users to understand the composition of the products they are using. For example, wind turbines incorporate thousands of lines of code (software) used for the turbine's design, planning, operation, and analytics in addition to the numerous hardware subcomponents that construct it. Due to the highly complex nature of software and hardware, knowledge of the components and subcomponents is required to mitigate cyber vulnerabilities and defend against cyberattacks. There has also been a transformation from a traditional linear supply chain into a global, dynamic, diverse, and interconnected system. The digitization of the supply chain makes it easier to find and exploit vulnerabilities. Critical infrastructures (e.g., power grids, oil, natural gas, water, and wastewater) rely on OT to function, and if the OT is compromised, equipment damage and potential interruption of services could result. A significant security measure to protect OT systems from disruption is to develop a supply chain bill of materials (BoM) corresponding to the software and hardware used in OT, along with attestations amongst vendors and asset owners. A supply chain BoM is a proactive way to understand the inherent vulnerabilities in the system and mitigate them in advance of being exploited. BoMs bolster the trust placed in the digital infrastructure and enhance software supply chain security by sustaining the management of component obsolescence and compliance, along with the seclusion of unsafe segments of a specific product. Adopting BoM tools is becoming increasingly important across various government sectors, as evidenced by the recent U.S. executive order on cybersecurity (NIST 2021). This paper aims to classify BoMs based on structure, functionality, component type, and architecture. The work also discusses case studies to further highlight the benefits of BoMs. In addition, it identifies missing pieces in existing BoM implementations so that future research may identify bounds on where it could expect to make improvements and directly enable researchers to identify promising areas for exploration. Further, the authors provide valuable recommendations to tool developers, researchers, and standardizing organizations (policymakers), additionally benefitting critical infrastructure owners and government executives. This aids in paving a path for future work, thereby, providing suggestions to determine a tool for consumers that best suit their needs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Efficient continuous-wave normal conducting accelerator for industrial applications

A normal conducting, high power, high efficiency copper linear accelerator prototype is being developed for industrial applications. The system will be powered by low-cost high-efficiency magnetron RF sources and will use a gridded thermionic cathode electron gun. Leveraging the significant accelerator expertise at JLab and industry partners, these technologies will be combined to deliver high-power (>100 kW) electron beams with energies of 1 MeV or higher that are cost-effective to produce and operate. The design is modular such that energy and power can be increased by adding additional sections as required. The status of the design, prototype fabrication and plans for a beam demonstration at JLab are described.

Accelerator Physics↗
DOE OSTI2024

Bench-Scale Testing of Monolithic PPI Structured Contactors for Direct Air Capture of CO 2

The overall project objective was to develop, optimize and bench-scale test the integrated embodiment of a leading direct air capture (DAC) sorbent composition, linear-poly(propylenimine) (l-PPI), in a structured material system, specifically a monolithic contactor for achieving low pressure drop, to justify its further scale-up in a subsequent program. The project team was led by CORMETECH, a leading monolithic gas/solid contactor manufacturer for environmental applications, and included Global Thermostat, a well-known DAC start-up, and Georgia Tech, a prominent US academic institution in carbon capture and DAC. The monolithic contactor for l-PPI in this project employed an advanced, cost-effective fabrication technique from CORMETECH, namely a porous monolith substrate directly impregnated with the amine, an approach contrasting to the traditional wash-coating methods used for activating monoliths. The l-PPI monolithic contractor is targeted for use in DAC systems such as the Global Thermostat DAC process, which involves the cyclic operation of ambient air flow over a monolithic amine contactor followed by steam-mediated thermal desorption and CO 2 collection, maximizing volumetric productivity while reducing the auxiliary power required to capture CO 2 from air. The Global Thermostat DAC process currently utilizes poly(ethyleneimine) (PEI) materials as baseline sorbents, which are subject to oxidative degradation and slow capacity fade at the elevated temperature required for efficient CO 2 removal. Georgia Tech, with Global Thermostat, had discovered in lab-scale efforts that l-PPI had superior resistance to oxidative degradation compared to PEI, with similar CO 2 adsorption capacity. While the novel l-PPI sorbent offers significant advantages over PEI for DAC, including potential design simplification and increased process efficiency, it is not commercially available, and prior to this project, had not yet been evaluated on the bench-scale in a structured contactor.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impacts of Year-to-Year Weather Variability and Inter-Panel Spacing on Crop Yields in a Massachusetts Agrivoltaics System

This presentation summarizes recently published work on agrivoltaic field data and irradiance modeling. The body of agrivoltaic field data is still growing, and crop responses to different solar configurations under different local climates are highly varied. We investigate the impact of adding spacing between adjacent solar panels in a fixed-tilt system to improve light diffusion to crops. For four crops (broccoli, peppers, kale, Swiss chard) grown across 3 years in an agrivoltaic system in Massachusetts, we found that only kale had a linearly increasing trend as the inter-panel spacing increased from 0.6 m to 1.5 m (2 ft to 5 ft). However, there were significant year-to-year differences in the yield of agrivoltaic versus control fields. Agrivoltaic and full sun fields produced equivalent yields in a hot, dry year, whereas the full-sun control beds produced more salable yield for all four crops in a warm, wet year. This demonstrates variability of agricultural outcomes and the need for more multi-year studies to ensure agrivoltaic impacts are not under- or overestimated.

14 SOLAR ENERGY↗