Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Programming, Linear”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Autodifferentiable Spectrum Model for High-dispersion Characterization of Exoplanets and Brown Dwarfs

We present an autodifferentiable spectral modeling of exoplanets and brown dwarfs. This model enables a fully Bayesian inference of the high-dispersion data to fit the ab initio line-by-line spectral computation to the observed spectrum by combining it with the Hamiltonian Monte Carlo in recent probabilistic programming languages. An open-source code, ExoJAX, developed in this study, was written in Python using the GPU/TPU compatible package for automatic differentiation and accelerated linear algebra, JAX. We validated the model by comparing it with existing opacity calculators and a radiative transfer code and found reasonable agreements for the output. As a demonstration, we analyzed the high-dispersion spectrum of a nearby brown dwarf, Luhman 16 A, and found that a model including water, carbon monoxide, and H 2 /He collision-induced absorption was well fitted to the observed spectrum (R = 10 5 and 2.28-2.30 μm). As a result, we found that T0 = 1295$^{-32}_{+35}$ K at 1 bar and C/O = 0.62 ± 0.03, which is slightly higher than the solar value. This work demonstrates the potential of a full Bayesian analysis of brown dwarfs and exoplanets as observed by high-dispersion spectrographs and also directly imaged exoplanets as observed by high-dispersion coronagraphy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mini Micropulse Lidar measurements in vertical and scanning modes

The Atmospheric Radiation Measurement (ARM) program operates ten Sigma Space (currently Droplet Measurement Technologies/DMT) micropulse lidar systems (MPL-4B-IDS/POL-FS series), which detect the backscattered light in two orthogonal linear polarization channels allowing to distinguish between spherical (e.g. water) and non-spherical (e.g. ice) scattering particles. The ARM MPLs have been providing valuable measurements to the aerosol and cloud research community for well over a decade. DMT’s Mini-Micropulse lidar (MiniMPL) is a compact version of the standard MPL. It reduces the power-aperture product to minimize cost, size, weight, power requirements and has all the electronics integrated with optics in a single box. Although as a result, it has a lower detection range than the MPL, limiting the measurements to the troposphere, its performance is promised to match that of the standard MPL within the troposphere. The MiniMPL weighs 13 kg, and the complete system, consisting of an optical transceiver and a laptop running the data acquisition and post-processing software, fits in a storm case, allowing its transport on domestic or international flights. This can be especially valuable for short field campaigns or Intensive Operational Periods, where time is of the essence when replacing a faulty instrument. In addition, the scanning polarization lidar capability (MiniMPL-SCAN option) could be of great interest for satisfying ARM measurement needs for certain research topics related to clouds and precipitation. To evaluate its full potential in comparison with co-located remote sensing instruments including the ARM MPL-4B-IDS/POL-FS at SGP, a miniMPL was hosted at the ARM Southern Great Plains (SGP) mega-site in 2020 for approximately 3 months. During this time the miniMPL was tested in both vertical sampling and scanning modes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Comparison of Approach-to-Critical Results in Current and Pulse Mode for Systems with High Starter Neutron Rates

Reactors and critical assemblies use a variety of detection systems to monitor the neutron population. The count rate is proportional to the neutron flux present at the location of the detector. When such systems are placed external to an assembly, it is often assumed that the relative leakage multiplication will be proportional to the detector count rate (assuming that the source term, system geometry, and detector placement have not changed). Such systems are often used in an approach-to-critical during reactor startup to ensure that the critical configuration is well predicted. Various types of detectors have been used during an approach-to-critical. These include 3 He, BF 3 , ion chambers, fission chambers, and fission foils. Any of these types of systems (or others) should work well when adequate counting statistics are available. These detector systems can be operated in either pulse or current mode. The National Criticality Experiments Research Center (NCERC) has two detection systems that are commonly used in critical assembly operations. The start-up (referred to as "SU" in this work) system is made up of 3 He proportional counters in pulse mode and the linear counter system (referred to as "LC" in this work) consists of compensated ion-chambers in current mode. Typically the SU system is used for approach-to-critical operations and the LC system is only used at/above delayed critical ( k eff = 1). This work investigates the use of the LC system for an approach-to-critical. It has been long hypothesized that such an approach would be feasible for systems with high starter neutron rates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

R-matrix Resolved Resonance Region Evaluation of 140,142 Ce

Oak Ridge National Laboratory completed the resolved resonance region (RRR) evaluation of the two most abundant cerium isotopes, 140 Ce (88.45%) and 142 Ce (11.11%), as requested by the US Nuclear Criticality Safety Program. These evaluations are based on recent high-resolution transmission and capture measurements performed on nat Ce and highly enriched 142 Ce samples at the JRC-Geel Linear Accelerator facility, as well as measured thermal constants available from the EXFOR database. Starting from the resonance parameters of the ENDF/B-VIII.0 library followed by a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information were derived by fitting these measured data using the Reich–Moore approximation of the R -matrix theory, as implemented in the SAMMY code system. The 140 Ce RRR upper energy limit was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. Updated statistical properties were obtained for the new evaluations and compared to those derived from the ENDF/B-VIII.0 nuclear data library. The new evaluation work improved some of the discrepancies found in previous work, such as the capture resonance integral and stellar Maxwellian-averaged cross sections. These integral quantities were mainly derived from the fit of the latest measured data, especially the neutron capture yield data for 142 Ce isotope.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Strategies for Constructing and Operating DNA Origami Linear Actuators

Abstract Linear actuators are ubiquitous components at all scales of engineering. DNA nanotechnology offers a unique opportunity for bottom‐up assembly at the molecular scale, providing nanoscale precision with multiple methods for constructing and operating devices. In this paper, DNA origami linear actuators with up to 200 nm travel, based on a rail threading a topologically locked slider, are demonstrated. Two strategies, one‐ and two‐pot assembly, are demonstrated whereby the two components are folded from one or two DNA scaffold strands, respectively. In order to control the position of the slider on the rail, the rail and the inside of the slider are decorated with single‐stranded oligonucleotides with distinct sequences. Two positioning strategies, based on diffusion and capture of signaling strands, are used to link the slider reversibly to determined positions on the rail with high yield and precision. These machine components provide a basis for applications in molecular machinery and nanoscale manufacture including programmed chemical synthesis.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Differential response of the photosynthetic machinery to dehydration in older and younger resurrection plants

Abstract A group of vascular plants called homoiochlorophyllous resurrection plants evolved unique capabilities to protect their photosynthetic machinery against desiccation-induced damage. This study examined whether the ontogenetic status of the resurrection plant Craterostigma pumilum has an impact on how the plant responds to dehydration at the thylakoid membrane level to prepare cells for the desiccated state. Thus, younger plants (<4 months) were compared with their older (>6 months) counterparts. Ultrastructural analysis provided evidence that younger plants suppressed senescence-like programs that are realized in older plants. During dehydration, older plants degrade specific subunits of the photosynthetic apparatus such as the D1 subunit of PSII and subunits of the cytochrome b6f complex. The latter leads to a controlled down-regulation of linear electron transport. In contrast, younger plants increased photoprotective high-energy quenching mechanisms and maintained a high capability to replace damaged D1 subunits. It follows that depending on the ontogenetic state, either more degradation-based or more photoprotective mechanisms are employed during dehydration of Craterostigma pumilum.

Plant Sciences↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Near-Optimal Performance of Stochastic Model Predictive Control

Here, this article presents a regret analysis for stochastic model predictive control (SMPC) in linear systems with quadratic performance index and additive and multiplicative uncertainties. Under a finite support assumption, the problem can be cast as a finite-dimensional quadratic program, but the problem becomes quickly intractable as the problem size grows exponentially in the horizon length. SMPC aims to compute approximate solutions by solving a sequence of problems with truncated prediction horizons and committing the solution in a receding-horizon fashion. Although this approach is widely used in practice, its performance relative to the optimal solution is not well understood. This article reports for the first time a rigorous near-optimal performance guarantee of SMPC: under stabilizability and detectability conditions, the regret of SMPC is exponentially small in the prediction horizon length, allowing SMPC to achieve near-optimal performance at a substantially reduced computational expense.

93E20, 93B45↗
DOE OSTI2025

An Agile/XP software development process for modernizing the accelerator control system at Fermilab

Fermilab is undergoing the most ambitious upgrade to its accelerator control system of the 21st century. As part of the ACORN project, hundreds of legacy control system applications written in C/C++ will be re-imagined and developed from the ground up. In addition, applications to support Fermilab’s new super-conducting linear accelerator are already under construction. To manage the development of modern controls applications, the Controls department has adopted an Agile software development process based on eXtreme Programming. In this paper we will describe our process and detail our experience applying it to the development of two case studies.

Diamond, John [Fermilab]↗
DOE OSTI2022

Purification and Concentration of 99 Mo from a UREX raffinate

The medical isotope 99m Tc ( t ½ = 6.0 h) accounts for over 80% of isotopes used in diagnostic nuclear medicine today and is normally derived from its transient equilibrium parent 99 Mo ( t ½ = 66 h). To reduce U.S. dependence on 99 Mo derived from HEU by foreign suppliers, the NNSA/M3 program—under the American Medical Isotope Production Act—has been tasked with facilitating the work of domestic 99 Mo suppliers that do not utilize HEU. Superconducting electron linear accelerators with high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that can induce photonuclear reactions and LEU fission. After sufficient production intervals, targets can be rotated out and processed while another batch is irradiated. The process flow chemistry is shown in Figure 1. The irradiated U 3 O 8 targets are retrieved and dissolved in HNO 3 ; the volatile fission products are expelled and captured during this step. The HNO 3 liquor bearing the uranyl, fission and activation products is injected into the UREX suite of liquid-liquid extraction banks. The uranyl is partitioned in typical fashion using tri-$\textit{n}$-butyl phosphate (TBP) in a hydrocarbon diluent. The raffinate—comprising mainly 99 Mo and fission products—is then injected into the MoLLE (Molybdenum Liquid-Liquid Extraction) flowsheet, where the Mo is selectively extracted by an organophosphorous acid extractant such as di-(2-ethylhexyl phosphoric acid) (HDEHP). Trace amounts of Nb, Zr, Np, I, and Te are co-extracted. The Mo is stripped using acetohydroxamic acid (AHA) and fed onto an anion exchange column. Following a series of hydroxide, HCl, and oxalic acid wash steps, the Mo can be recovered in NaOH/NaCl. A manuscript describing this process chemistry in more detail was recently published.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

140,142 Ce Neutron Cross Section Resolved Resonance Region Evaluation [Abstract]

A resolved resonance region evaluation of 140,142 Ce has been carried out by Oak Ridge National Laboratory. Requested by the US Nuclear Criticality Safety Program, this evaluation is based on recent high-resolution transmission and capture high-resolution measurements of nat Ce and 142 Ce conducted at JRC-GEEL at the Geel Linear Accelerator facility, as well as recently measured thermal constants available from the EXFOR database. Starting from the resonance parameters from the ENDF/B-VIII.0 library and following a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information was derived by the fit of these experimental datasets using the Reich-Moore approximation of the R-matrix theory as implemented in the SAMMY code system. The resolved resonance region upper energy limit for 140 Ce was kept at 200 keV while the 142 Ce resonance region was extended from 13 to 26 keV. This new evaluation was found to be in good agreement not only with several integral quantities of interest to the reactor physics community, but also with the stellar Maxwellian-averaged cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Computer Vision on Edge Devices for the Short Term Prediction of Cloud Cover

Edge Computing and IoT are important pieces of today's technological landscape. Here, we build a low-cost IoT sensor for sky imaging and program it using AWS GreenGrass, one of the leading IoT platforms. We demonstrate remote reprogramming of this device to load software that predicts sun shading events through the linear advection method, which is a baseline algorithm that can be used to benchmark algorithmic improvements in future work. Some future directions for sky imaging research are enumerated.

14 SOLAR ENERGY↗
DOE OSTI2022

Koopman Model Predictive Control for Eco-Driving of Automated Vehicles

In this paper, we develop a data-driven process for building a model predictive control (MPC) for eco-driving of automated vehicles. The process involves performing system identification in which the non-linear vehicle dynamics model is approximated by the Koopman operator, a linear predictor of higher state-dimension, in a data-driven framework. This approach allows us to formulate the eco-driving problem in a constrained quadratic program that leads to a computationally fast MPC. The MPC is then implemented as a closed-loop control of an electric vehicle in numerical simulations for demonstration.

autonomous vehicle↗
DOE OSTI2022

Economic dispatch for electricity merchant with energy storage and wind plant: State of charge based decision making considering market impact and uncertainties

Here this paper investigates how the market impact of electricity merchants and uncertainty of wind generation affect their co-optimized scheduling policy, specifically for merchants who have both energy storage and wind plants. In the existing literature, merchants' trading actions are usually assumed not to affect market prices; however, a large-scale energy storage merchant's actions can affect market prices. To this end, we approximate the electricity price by a linear function of the quantity of power traded by the merchant in the reward function to achieve decision-making incorporating the market impact. This paper utilizes the dynamic programming approach to analyze merchants' optimal multi-period decision-making incorporating market impact, uncertain wind generation, and energy storage constraints. First, our results demonstrate that for a merchant with co-located energy storage facilities and wind power plants, the energy storage's feasible state of charge (SOC) range can be segmented into four possible sub-ranges by three analytically developed SOC reference points. The unique optimal trading decision can be achieved by comparing the current energy inventory and the SOC references of the next period. Second, our results show that market impact and uncertainties substantially change the optimal storage scheduling policy by impacting the values of the reference points. To mitigate the negative effect of the merchant's market impact on buying and selling actions, the merchant may reduce the amount of generating or pumping electricity each period to maximize profit. Moreover, we identify and investigate the trade-off between market price and transaction quantity. Our findings provide co-optimized scheduling guidance for electricity merchants with co-located energy storage and renewable power plants systems.

17 WIND ENERGY↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2025

Practical Implementation of GPU-based Computing at the Grid Edge for Resilience Scenarios

This paper presents a practical implementation of GPU-accelerated computing at the grid edge to enhance power system resilience through next-generation smart meters. Advanced Metering Infrastructure (AMI) systems rely predominantly on centralized processing architectures, which limit real-time response capabilities during grid disturbances. This work proposes the integration of GPU-enabled computational platforms directly within smart meter to enable local execution support for power system analytics, fault detection algorithms, and optimization routines. The proposed framework uses the Julia programming language to leverage highperformance parallel computing capabilities while maintaining code portability and development efficiency. We use two experimental scenarios to benchmark the computational feasibility of this approach: sparse linear system solutions representative of power flow analyses, and multi-stage production cost simulations incorporating unit commitment and economic dispatch operations. Results demonstrate that computationally intensive power system algorithms, such as those supporting resilience scenario calculations, can be effectively executed at the distribution edge using commercially available embedded GPU hardware. Keywords—GPU acceleration, edge computing, smart meters, grid resilience, AMI, resilience.

De Souza, Reubun [School of Electrical Engineering↗
DOE OSTI2025

Challenges in Tracking Waste Reduction Performance Improvement in Manufacturing Plants

The recently released Circularity Gap Report 2023 by the Circle Economy Foundation states that the circularity score for the global economy is declining. The US Environmental Protection Agency (EPA) tracked municipal solid waste from 1960 to 2018 and found that 50% of the waste was destined for landfills. EPA estimates that US industry is responsible for 2.7 Gt of solid nonhazardous waste annually in the US mostly linear economy model. The circular economy framework aims to decouple economic value generation from the extraction of virgin materials from nature. The linear model of material extraction and disposal at the end of life is highly unsustainable. Manufacturing companies are adopting ambitious waste reduction targets to achieve sustainability. Through the Better Plants program, US DOE has established the Waste Reduction Network, which offers technical assistance to partners to achieve their ambitious waste reduction goals. Basic requirements of establishing a target include identifying a baseline, quantifying waste performance, and measuring progress over time. One problem faced by industry is unstandardized metrics for quantifying waste performance that may not be well suited to demonstrate progress. This paper studies traditional methods used to measure waste performance and highlights advantages, disadvantages, and limitations of each method. It also examines the suitability of applying the methods in different manufacturing circumstances. Finally, the paper presents a case study of a large manufacturer that faced inconsistencies in its tracked measurement metric. A solution was proposed and implemented to alter the methodology to enable more accurate waste performance tracking against a baseline.

circular economy↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗