Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear optimization”

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,387 records · Page 77

Stochastic optimization of cellulosic biofuel supply chain incorporating feedstock yield uncertainty

The global goal to reduce dependence on fossil fuels and to mitigate greenhouse gas emissions has resulted in research focused on environment friendly and socio-economically sustainable renewable energy sources. However, commercial production of bio-energy is constrained by biomass supply uncertainty and associated costs. This study presents an integrated approach to determining the optimal biofuel supply chain considering biomass yield uncertainty. A two-stage stochastic mixed integer linear programming is utilized to minimize the expected system cost while incorporating yield uncertainty in the strategic level decisions related to biomass production and biorefinery investment. Applicability of the stochastic model is illustrated through a case study of switchgrass-based biofuel in west Tennessee.

09 BIOMASS FUELS↗

CHARACTERIZATION OF DISCONTINUITIES IN HIGH-DIMENSIONAL STOCHASTIC PROBLEMS ON ADAPTIVE SPARSE GRIDS

In this paper we present a set of efficient algorithms for detection and identification of discontinuities in high dimensional space. The method is based on extension of polynomial annihilation for edge detection in low dimensions. Compared to the earlier work, the present method poses significant improvements for high dimensional problems. The core of the algorithms relies on adaptive refinement of sparse grids. It is demonstrated that in the commonly encountered cases where a discontinuity resides on a small subset of the dimensions, the present method becomes optimal , in the sense that the total number of points required for function evaluations depends linearly on the dimensionality of the space. The details of the algorithms will be presented and various numerical examples are utilized to demonstrate the efficacy of the method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermoelectric, electronic, optical and chemical bonding properties of Ba{sub 2}PrRuO{sub 6}: At temperature 7 K and 150 K

Highlights: • DFT-FPLAPW method used for calculating the electronic structure. • The Fermi surface of BPRO (7 K and 150 K) is also calculated. • The complex dielectric function has been calculated. • Thermoelectric properties were also calculated using BoltzTraP code. • Power factor shows that both compounds are good thermoelectric materials at 600 K. - Abstract: We present first principles calculations of the band structure, density of states, electronic charge density, Fermi surface and optical properties of Ba{sub 2}PrRuO{sub 6} single crystals at two different temperatures. The atomic positions were optimized by minimizing the forces acting on the atoms. We have employed the full potential linear augmented plane wave method within local density approximation, generalized gradient approximation and Engel–Vosko generalized gradient approximation to treat the exchange correlation potential. The calculation shows that the compound is superconductor with strong hybridization near the Fermi energy level. Fermi surface is composed of two sheets. The calculated electronic specific heat capacities indicate, very close agreement with the experimental one. The bonding features of the compounds are analyzed using the electronic charge density in the (1 0 0) and (0–10) crystallographic planes. The dispersion of the optical constants was calculated and discussed. The thermoelectric properties are also calculated using the BoltzTrap code.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modern deep neural networks for Direct Normal Irradiance forecasting: A classification approach

The escalating energy demand and the adverse environmental impacts of fossil-fuel use necessitate a shift towards cleaner and renewable alternatives. Concentrated Solar Power (CSP) technology emerges as a promising solution, offering a carbon-free alternative for power generation. The efficiency and profitability of CSP depend on the Direct Normal Irradiance (DNI) component of solar radiation; hence, accurate DNI forecasting can help optimize CSP plants’ operations and performance. The unpredictable nature of weather phenomena, particularly cloud cover, introduces uncertainty into DNI projections. Existing DNI forecasting models use meteorological factors, which are both challenging to estimate numerically over short prediction windows and expensive to model through data at a sufficiently high spatial and temporal resolution. This research addresses the challenge by presenting a novel approach that formulates DNI prediction as a multi-class classification problem, departing from conventional regression-based methods. The primary objective of this classification framework is to identify optimal periods aligning with specific operational thresholds for CSP plants, contributing to enhanced dispatch optimization strategies. We model the DNI classification problem using four advanced deep neural networks – rectified linear unit (ReLU) networks, 1D residual networks (ResNets), bidirectional long short-term memory (BiLSTM) networks, and transformers – achieving accuracies up to 93.5% without requiring meteorological parameters.

14 SOLAR ENERGY↗

Hydrogen Analysis by Gas Chromatography–Mass Spectrometry

The detection of hydrogen in complex gas mixtures is essential for many applications. Conventional approaches such as gas chromatography (GC) with thermal conductivity detection (TCD) and residual gas analyzers (RGAs) face significant limitations: TCD exhibits poor response when helium is used as a carrier gas, whereas RGAs lack chromatographic separation, preventing reliable quantification of hydrogen because of interference from other species. Traditional GC methods rely on dual-column configurations with packed or molecular-sieve porous layer open tubular (PLOT) columns to separate hydrogen from O2, N2, CO, CO2, CH4, and other hydrocarbons, increasing system complexity and limiting compatibility with mass spectrometry (MS) detection. In this work, we developed and validated a robust GC–MS method capable of directly detecting and quantifying hydrogen using electron ionization (EI) without dopants, reagent gases, or ion–molecule reaction schemes. By integrating a modified EI source and a cryogenically cooled single capillary column configuration, we achieved baseline separation of hydrogen from all major permanent gases and hydrocarbons in a refinery gas mixture using helium as the carrier gas. The method demonstrated excellent linearity, high sensitivity, and exceptional reproducibility. Adjustable sample-loop volumes and split ratios enabled optimization of peak shape and signal-to-noise performance for trace-level and percent-level hydrogen concentrations. Beyond hydrogen quantitation, the method provides simultaneous compositional profiling of other gases in a mixture in a single run, making it valuable for a wide range of tasks.

Lobodin, Vlad [ORNL]↗

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗

Optimized symmetry functions for machine-learning interatomic potentials of multicomponent systems

Current machine-learning methods to reproduce ab initio potential energy landscapes suffer from an unfavorable computational scaling with respect to the number of chemical species. In this work, we propose a new approach by using optimized symmetry functions to explore similarities of structures in multicomponent systems in order to yield linear complexity. Here, we combine these symmetry functions with the charge equilibration via neural network technique, a reliable artificial neural network potential for ionic materials, and apply this method to study alkali-halide materials MX with 6 chemical species (M = {Li, Na, K} and X = {F, Cl, Br}). Our results show that our approach provides good agreement both with experimental and DFT reference data of many physical and structural properties for any chemical combination.

74 ATOMIC AND MOLECULAR PHYSICS↗

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗

Universal control in bosonic systems with weak Kerr nonlinearities

Resonators with weak single-photon self-Kerr nonlinearities can theoretically be used to prepare Fock states in the presence of a loss much larger than their nonlinearities. Two necessary ingredients are large displacements and a two-photon (parametric) drive. Here, in this study, we find that these systems can be controlled to achieve any desired gate operation in a finite-dimensional subspace (whose dimensionality can be chosen at will). Moreover, we show that the two-photon driving requirement can be relaxed and that full controllability is achievable with only one-photon (linear) drives. We make use of both Trotter-Suzuki decompositions and gradient-based optimization to find control pulses for a desired gate, which reduces the computational overhead by using a small blockaded subspace. We also discuss the infidelity arising from input power limitations in realistic settings, as well as from corrections to the rotating-wave approximation. Our universal control protocol opens the possibility for quantum information processing using a wide range of lossy systems with weak nonlinearities.

Yuan, Ming [Univ. of Chicago, IL (United States)] ↗

Semidefinite programming algorithm for the quantum mechanical bootstrap

Here, we present a semidefinite programming algorithm to find eigenvalues of Schrödinger operators within the bootstrap approach to quantum mechanics. The bootstrap approach involves two ingredients: a nonlinear set of constraints on the variables (expectation values of operators in an energy eigenstate), plus positivity constraints (unitarity) that need to be satisfied. By fixing the energy we linearize all the constraints and show that the feasibility problem can be presented as an optimization problem for the variables that are not fixed by the constraints and one additional slack variable that measures the failure of positivity. To illustrate the method we are able to obtain high-precision, sharp bounds on eigenenergies for arbitrary confining polynomial potentials in one dimension.

97 MATHEMATICS AND COMPUTING↗

Homotopy Solver

This software implements parallel versions of an interior-point solver, based on the publicly available ipopt solver. Here we have full control over the linear solver and our algorithm is fully parallel thus enabling scalability to large-scale optimization problems. This package also has a parallel implementation of a homotopy solver developed under the scalable methods for contact LDRD project 23-ERD-017. This solver is an mfem-based implementation of algorithm described in ``A filter trust-region Newton continuation method for nonlinear complementarity problems''. Cosmin G. Petra, Nai-Yuan Chiang, Jingyi Wang, Tucker Hartland, and Michael Puso (submitted), LLNL-JRNL-869761.

Hartland, Tucker [Lawrence Livermore National Labo↗

Subcontract Report: Diffusion Mechanisms and Bond Dynamics in Solid Electrolyte Ion-Conductors

We employ first-principles molecular dynamics simulations and Maximally Localized Wannier Function (MLWF) analysis to explore how halide substitution and nano-phase microstructures affect diffusivity, through the activation energy barrier - E a and D 0 , in the solid electrolyte Li 3 InBr 6-x Cl x . We find that nano-phase microstructures with x=3 (50-50 Br-Cl) mixed composition have a higher diffusivity compared to x=2 and x=3 solid solutions. There is a positive linear relationship between ln(D 0 .) and E a , which suggests that for superionic conductivity optimizing both the activation energy and the D 0 is important. Bond frustration due to mismatch in crystal geometry and ideal coordination number leads to especially high diffusivity through a high D 0 in the x=3 composition.

25 ENERGY STORAGE↗

A Sensitive Electrochemical Cholinesterase-Inhibiting Biosensor for Organophosphorus Pesticides Based on Ti 3 C 2 TX MXene Quantum Dots

Organophosphorus pesticides (OPs) pose significant environmental and health risks due to their widespread use and toxicity, primarily by inhibiting acetylcholinesterase. Traditional detection methods are often slow and costly, highlighting the urgent need for advanced, sensitive, and accessible technologies. This study developed a highly sensitive electrochemical cholinesterase-inhibiting biosensor for OP pesticides, utilizing Ti 3 C 2 T x MXene Quantum Dots (MQDs), which was synthesized via a hydrothermal method. The biosensor’s performance was characterized using electrochemical impedance spectroscopy, differential pulse voltammetry (DPV), and cyclic voltammetry. DPV proved to be the optimal technique, exhibiting an ultralow detection limit of 1 × 10 −17 M and a wide linear range (10 −14 –10 −8 M) for chlorpyrifos (a model OP) with an estimated inhibition constant of 62 nM. The biosensor demonstrated high selectivity for OPs (chlorpyrifos, acephate, glyphosate) over a non-target pyrethroid (permethrin), confirmed by distinct electrochemical signatures and compared to in vitro cholinergic activity assays in bean beetle homogenates. The enhanced performance is attributed to the high surface-to-volume ratio, quantum confinement effects, and superior conductivity of the MQDs, as well as the robust enzyme immobilization facilitated by glutaraldehyde cross-linking and a chitosan matrix. This work presents a promising platform for rapid, sensitive, and selective detection of OP pesticides, with potential applications in environmental monitoring and public health protection.

Chemistry↗

Dipole corrector magnets for the LBNE beam line

The conceptual design of a new dipole corrector magnet has been thoroughly studied. The planned Long-Baseline Neutrino Experiment (LBNE) beam line will require correctors capable of greater range and linearity than existing correctors, so a new design is proposed based on the horizontal trim dipole correctors built for the Main Injector synchrotron at Fermilab. The gap, pole shape, length, and number of conductor turns remain the same. To allow operation over a wider range of excitations without overheating, the conductor size is increased, and to maintain better linearity, the back leg thickness is increased. The magnetic simulation was done using ANSYS to optimize the shape and the size of the yoke. The thermal performance was also modeled and analyzed.

43 PARTICLE ACCELERATORS↗

Quasi-Stochastic Approximation and Off-Policy Reinforcement Learning: Preprint

The Robbins-Monro stochastic approximation algorithm is a foundation of many algorithmic frameworks for reinforcement learning and often an efficient approach to solving (or approximating the solution to) complex optimal control problems. In many cases, however, practitioners are unable to apply these techniques because of an inherent high variance. This paper aims to provide a general foundation for 'quasi-stochastic approximation,' in which all of the processes under consideration are deterministic, much like quasi-Monte Carlo for variance reduction in simulation. The variance reduction can be substantial, subject to tuning of pertinent parameters in the algorithm. This paper introduces a new coupling argument to establish the optimal rate of convergence provided the gain is sufficiently large. These results are established for linear models and tested also in nonideal settings. A major application of these general results is a new class of reinforcement learning algorithms for deterministic state space models. In this setting, the main contribution is a class of algorithms for approximating the value function for a given policy, using a different policy designed to introduce exploration.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Probabilistic Forecast of All-Sky Solar Radiation Using Enhanced WRF-Solar: Preprint

This study presents enhancements of the Weather Research and Forecasting model with solar extensions (WRF-Solar) to provide probabilistic forecasts of solar radiation. Our approach builds ensemble WRF-Solar runs by introducing stochastic perturbations of variables that produce the largest uncertainties in predicting surface irradiance and clouds. The key variables are identified using tangent linear sensitivity analysis of six physics packages responsible for all-sky irradiance variability. An optimal strategy to stochastically perturb the selected variables is developed and applied to WRF-Solar to generate ensemble members for day-ahead solar prediction. The National Solar Radiation Database (NSRDB) is used to validate the ensemble forecast at arbitrary locations on the model grid. Preliminary results indicate that the proposed technique can potentially produce WRF-Solar ensembles providing reliable information of solar prediction uncertainty. This study describes the implemented methodology and initial results as well as future research to improve ensemble-based probabilistic forecasts with WRF-Solar.

all-sky↗

Thermal Analysis for SRS Saltstone Disposal Unit Facility

Saltstone Disposal Units (SDU) are long term disposal units for low activity grout produced from the solidification of decontaminated non-hazardous salt waste at SRS. The unit is a cylindrical disposal vault of 150 ft. in diameter and 23 ft. in height. The total volume the SDU can hold is 3 million gallons of grout. The hydration process of the grout material produces heat from an exothermic chemical reaction. If the grout reaches temperatures higher than the maximum safety temperature of 95 deg. C, cracks can start forming leading to contamination of the environment. The purpose of this study is to develop the optimal baseline domain boundary conditions and to perform a sensitivity analysis with respect to the baseline modeling conditions. The following assumptions were made for the modeling calculations: - The system geometry is axisymmetric; - The system was at steady state; - Ambient temperature was set to 100 deg. F; - Air was assumed to be an incompressible ideal gas; - Radiation was neglected for a conservative estimation; - Natural convection for the SDU vapor space is assumed to be turbulent because Rayleigh number based on typical operating conditions is higher than 10{sup 8}; - Water tables were not taken into consideration. The following sensitivity analyses were performed to find the optimal domain boundary: - Mesh size sensitivity; - Soil depth sensitivity - Soil depth varied from 30 ft to 250 ft with a soil thickness of 90 ft; - Soil thickness sensitivity - Soil thickness varied from 30 ft to 270 ft with a soil depth of 180 ft; - Soil conductivity sensitivity - Soil thermal conductivity varied from 0.25 W/m-K to 1 W/m-K; - Ambient Temperature Sensitivity - Ambient temperature varied from 70 deg. F to 110 deg. F. For the thermal calculations, the temperature of the bottom of the soil was set to 17.4 deg. C. The heat flux on the roof of the building and surface of the earth is 163.2 W/m{sup 2}. The soil sides are set to be adiabatic. The sensitivity analysis shows that the mesh size has no effect on the numerical results. This is mainly due to a small convection zone and large conduction zones in the modeling domain. The optimal domain boundary was found to be a soil thickness of 90 ft and a soil depth of 150 ft. The different sensitivity analysis performed here shows the optimal boundary domain to be used for thermal calculations is a soil depth of 150 ft beneath the SDU and a soil thickness of 90 ft. Anything larger would produce similar solutions but at a longer computational time to solve. The maximum temperature was more sensitive to soil properties as shown here. While the ambient temperature affected the max temperature linearly. Most of the heat transfer takes place downward vertically between the SDU and the soil beneath it. Using the optimal domain boundary determined by the sensitivity analyses, the maximum thermal loading to keep the maximum grout temperature below 95 deg. C is 1.646 W/m{sup 3}.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Laser-Induced Spectrochemical Assay for Uranium Enrichment (LISA-UE)

Uranium hexafluoride (UF6) is the uranium compound typically involved in uranium enrichment process. As the first line of defense against nuclear proliferation, accurate determinations of the uranium enrichment ratio in UF6 are critical for materials verification, accounting and safeguards. Shipping gaseous UF6 samples off-site for analysis with mass spectrometry is cumbersome and costly, and results are not available for some time (months). In-field UF6 enrichment assay has the potential to substantially reduce the time, logistics and expense of sample handling. At present, COMPUCEA is the only accepted method for UF6 enrichment assay in the field. Laser-Induced Spectrochemical Assay for Uranium Enrichment (LISA-UE) is an all-optical (based on laser induced plasma emission) analytical technique intended for fieldable, accurate, precise and rapid UF6 enrichment assay. In its operation, laser induced plasma is created directly in the gaseous UF6 sample. Because different U isotopes emit at slightly different wavelengths, the isotopic information of the UF6 sample is inherently encoded in the atomic emission from the plasma. Isotopic emissions from 235U and 238U are measured simultaneously, which eliminate correlated noise from the laser induced plasma. Isotopic information of the UF6 sample can be extracted from the acquired spectrum with theoretical multi-variable non-linear spectral fitting. To date, advances made by the LISA-UE research team include optimization of the spectral window for direct gaseous UF6 enrichment assay with laser induced plasma, development of data reduction algorithms, and demonstrations of the LISA-UE technique with gaseous UF6 samples. In this presentation, the technical aspect of LISA-UE will be overviewed, the data reduction algorithm will be described, and performance of the technique will be discussed.

Chan, George↗