Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal sizing method”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Digital memcomputing: From logic to dynamics to topology

Digital memcomputing machines (DMMs) are a class of computational machines designed to solve combinatorial optimization problems. A practical realization of DMMs can be accomplished via electrical circuits of highly non-linear, point-dissipative dynamical systems engineered so that periodic orbits and chaos can be avoided. A given logic problem is first mapped into this type of dynamical system whose point attractors represent the solutions of the original problem. A DMM then finds the solution via a succession of elementary instantons whose role is to eliminate solitonic configurations of logical inconsistency (“logical defects”) from the circuit. By employing a supersymmetric theory of dynamics, a DMM can be described by a cohomological field theory that allows for computation of certain topological matrix elements on instantons that have the mathematical meaning of intersection numbers on instantons. We discuss the “dynamical” meaning of these matrix elements, and argue that the number of elementary instantons needed to reach the solution cannot exceed the number of state variables of DMMs, which in turn can only grow at most polynomially with the size of the problem. These results shed further light on the relation between logic, dynamics and topology in digital memcomputing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Optimized attenuated interaction: Enabling stochastic Bethe–Salpeter spectra for large systems

We develop an improved stochastic formalism for the Bethe–Salpeter equation (BSE), based on an exact separation of the effective-interaction W into two parts, W = (W – vW) + vW, where the latter is formally any translationally invariant interaction, vW(r – r'). When optimizing the fit of the exchange kernel vW to W, using a stochastic sampling W, the difference W – vW becomes quite small. Then, in the main BSE routine, this small difference is stochastically sampled. Furthermore, the number of stochastic samples needed for an accurate spectrum is then largely independent of system size. While the method is formally cubic in scaling, the scaling prefactor is small due to the constant number of stochastic orbitals needed for sampling W.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Variations of the physical parameters of the blazar Mrk 421 based on analysis of the spectral energy distributions

Abstract We report on the variations of the physical parameters of the jet observed in the blazar Mrk 421, and discuss the origin of X-ray flares in the jet, based on analysis of several spectral energy distributions (SEDs). The SEDs are modeled using the one-zone synchrotron self-Compton model, its parameters determined using a Markov chain Monte Carlo method. The lack of data at TeV energies means many of the parameters cannot be uniquely determined and are correlated. These are studied in detail. We find that the optimal solution can be uniquely determined only when we apply a constraint to one of four parameters: the magnetic field (B), the Doppler factor, the size of the emitting region, and the normalization factor of the electron energy distribution. We used 31 sets of SEDs from 2009 to 2014 with optical–UV data observed with UVOT/Swift and the Kanata telescope, X-ray data with XRT/Swift, and γ-ray data with the Fermi Large Area Telescope. The result of our SED analysis suggests that, in the X-ray faint state, the emission occurs in a relatively small area (∼1016 cm) with a relatively strong magnetic field (B ∼ 10−1 G). The X-ray bright state shows a tendency opposite to that of the faint state, that is, a large emitting area (∼1018 cm), probably downstream of the jet, and a weak magnetic field (B ∼ 10−3 G). The high X-ray flux was due to an increase in the maximum energy of electrons. On the other hand, the presence of two kinds of emitting areas implies that the one-zone model is unsuitable for reproducing at least part of the observed SEDs.

Yamada, Yurika↗
DOE OSTI2020

Learning to Solve Large-Scale Security-Constrained Unit Commitment Problems

Security-constrained unit commitment (SCUC) is a fundamental problem in power systems and electricity markets. In practical settings, SCUC is repeatedly solved via mixed-integer linear programming (MIP), sometimes multiple times per day, with only minor changes in input data. In this work, we propose a number of machine learning techniques to effectively extract information from previously solved instances in order to significantly improve the computational performance of MIP solvers when solving similar instances in the future. Based on statistical data, we predict redundant constraints in the formulation, good initial feasible solutions, and affine subspaces where the optimal solution is likely to lie, leading to a significant reduction in problem size. Computational results on a diverse set of realistic and large-scale instances show that using the proposed techniques, SCUC can be solved on average 4.3 times faster with optimality guarantees and 10.2 times faster without optimality guarantees, with no observed reduction in solution quality. Out-of-distribution experiments provide evidence that the method is somewhat robust against data-set shift. Summary of Contribution. The paper describes a novel computational method, based on a combination of mixed-integer linear programming (MILP) and machine learning (ML), to solve a challenging and fundamental optimization problem in the energy sector. The method advances the state-of-the-art, not only for this particular problem, but also, more generally, in solving discrete optimization problems via ML. We expect that the techniques presented can be readily used by practitioners in the energy sector and adapted, by researchers in other fields, to other challenging operations research problems that are solved routinely.

Machine Learning↗
DOE OSTI2022

Predicting Switchgrass Biomass Yields Using a Spectral Vegetation Index Derived from Multispectral Satellite Imagery

Successful scaling of perennial bioenergy crop production requires a landscape design that optimizes the benefits of finite lands for people, communities, and environments. Utilizing marginal areas is the key to sustainable bioenergy crop production (Ssegane et al., 2015, 2016). Marginal areas are often small-sized lands and unevenly distributed across the agricultural landscape (Ssegane et al., 2016); thus, a systematic, semi-automated remote sensing method is needed as an effective means of estimating bioenergy crop yields across landscapes. Argonne National Laboratory (Argonne) is currently developing a tool, Scaling Up Perennial Bioenergy Economics and Ecosystem Services Tool (SUPERBEEST), to identify marginal agricultural lands and quantify environmental and economic effects of perennial bioenergy crop production systems. The tool aims to provide users a path to foster the sustainable and productive integration of bioenergy crops in the Midwestern agricultural landscape. Reliable, cost-effective, and timely estimation of bioenergy crop yields using remote sensing would help calculate and track the success of integrated bioenergy crops in the landscape for those communities. Argonne previously conducted feasibility studies for estimating biomass yields for bioenergy feedstock, corn and perennial grass using spectral vegetation indices (SVIs)1 derived from optical imagery (Hamada et al., 2015, 2021). In both studies, SVIs, more specifically those sensitive to plant chlorophyll or nitrogen contents, showed potential for estimating or predicting biomass yields with a correlation of determination (R2) ranging from 0.54 to 0.96, indicating a value for further investigation as a viable means of quantifying bioenergy feedstock production across large landscapes. Thus, the goal of this study is to evaluate the feasibility of use of SVIs as a means of estimating or predicting switchgrass biomass yields at harvest using publicly available multispectral satellite imagery. The feasibility analysis was performed using four study areas of mature switchgrass located in Virginia. Objectives are to (1) examine the SVIs and establish their relationships with switchgrass biomass yields at harvest, (2) develop a parsimonious image processing model for predicting at-harvest yields by applying the relationships with the most promising spectral index and (3) map switchgrass yields predicted by the image processing model across the study sites. The calibration to field data will rely on switchgrass biomass yields determined by the baling method, representing a potential challenge to the analysis but an important practical aspect for future applications. With this research design, the study aimed to gain insights into enabling remote sensing-based estimation of bioenergy crop yields in a reliable, cost- effective, and timely manner across large, heterogeneous landscapes.

09 BIOMASS FUELS↗
NASA NTRS2002

Very Large Scale Optimization

The purpose of this research under the NASA Small Business Innovative Research program was to develop algorithms and associated software to solve very large nonlinear, constrained optimization tasks. Key issues included efficiency, reliability, memory, and gradient calculation requirements. This report describes the general optimization problem, ten candidate methods, and detailed evaluations of four candidates. The algorithm chosen for final development is a modern recreation of a 1960s external penalty function method that uses very limited computer memory and computational time. Although of lower efficiency, the new method can solve problems orders of magnitude larger than current methods. The resulting BIGDOT software has been demonstrated on problems with 50,000 variables and about 50,000 active constraints. For unconstrained optimization, it has solved a problem in excess of 135,000 variables. The method includes a technique for solving discrete variable problems that finds a "good" design, although a theoretical optimum cannot be guaranteed. It is very scalable in that the number of function and gradient evaluations does not change significantly with increased problem size. Test cases are provided to demonstrate the efficiency and reliability of the methods and software.

Vanderplaats, Garrett↗
NASA NTRS1994

A parametric sensitivity study for single-stage-to-orbit hypersonic vehicles using trajectory optimization

The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) capability reflect highly integrated airframe and propulsion systems. These designs are also known to exhibit a large degree of interaction between the airframe and engine dynamics. Consequently, even simplified hypersonic models are characterized by tightly coupled nonlinear equations of motion. In addition, hypersonic SSTO vehicles present a major system design challenge; the vehicle's overall mission performance is a function of its subsystem efficiencies including structural, aerodynamic, propulsive, and operational. Further, all subsystem efficiencies are interrelated, hence, independent optimization of the subsystems is not likely to lead to an optimum design. Thus, it is desired to know the effect of various subsystem efficiencies on overall mission performance. For the purposes of this analysis, mission performance will be measured in terms of the payload weight inserted into orbit. In this report, a trajectory optimization problem is formulated for a generic hypersonic lifting body for a specified orbit-injection mission. A solution method is outlined, and results are detailed for the generic vehicle, referred to as the baseline model. After evaluating the performance of the baseline model, a sensitivity study is presented to determine the effect of various subsystem efficiencies on mission performance. This consists of performing a parametric analysis of the basic design parameters, generating a matrix of configurations, and determining the mission performance of each configuration. Also, the performance loss due to constraining the total head load experienced by the vehicle is evaluated. The key results from this analysis include the formulation of the sizing problem for this vehicle class using trajectory optimization, characteristics of the optimal trajectories, and the subsystem design sensitivities.

Lovell, T. Alan↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1984

A nonlinear programming method for system design with results that have been implemented

A general nonlinear programming algorithm (NICO) is discussed. An academic optimization example is given. The NICO multi-input, multi-output control system design is discussed. NICO applications relative to launch vehicle autopilot design, space shuttle static balance, transient response criteria in the design of a reentry vehicle control system, and waterjet propulsion and lift system components sized to a large surface effect ship are noted.

Hauser, F.↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Control co-design under uncertainty for offshore wind farms: Optimizing grid integration, energy storage, and market participation

Offshore wind farms (OWFs) are set to significantly contribute to global decarbonization efforts. Developers often use a sequential approach to optimize design variables and market participation for grid-integrated offshore wind farms. However, this method can lead to sub-optimal system performance, and uncertainties associated with renewable resources are often overlooked in decision-making. Here, this paper proposes a control co-design approach, optimizing design and control decisions for integrating OWFs into the power grid while considering energy market and primary frequency market participation. Additionally, we introduce optimal sizing solutions for energy storage systems deployed onshore to enhance revenue for OWF developers over time. This framework addresses uncertainties related to wind resources and energy prices. We analyze five U.S. west-coast offshore wind farm locations and potential interconnection points, as identified by the Bureau of Ocean Energy Management (BOEM). Results show that optimized control co-design solutions can increase market revenue by 3.2% and provide flexibility in managing wind resource uncertainties.

Control Co-design↗
DOE OSTI2026

Sequential surface synthesis of dispersed sub-nanometer iridium on titanium nitride for acidic water oxidation

Maximizing iridium utilization while maintaining high oxygen evolution reaction (OER) performance remains a persistent challenge in acidic water electrolysis. Immobilizing Ir on conductive, acid-stable supports is promising, yet simultaneously achieving sub-nanometer size, high area coverage, and strong electronic coupling is difficult. Here, we report a sequential surface synthesis on titanium nitride (TiN) that yields uniformly distributed sub-nanometer Ir arrays (∼0.7 nm). Our method uses ethylenediaminetetraacetic acid (EDTA) as a temporal scaffold: it chemisorbs to TiN to install dense chelating sites, captures Ir 3+ ions, and confines Ir cluster growth. A subsequent thermal treatment at 500 °C in a reducing atmosphere removes the ligand shell, while preserving ultrasmall particle size and establishing direct Ir–TiN electronic coupling. The optimized catalyst exhibits mixed Ir 0 /Ir x+ coordination with low charge-transfer resistance (R ct = 19.2 Ω), delivering a mass activity of 342 A g Ir −1 at 1.54 V in acidic electrolyte. In situ X-ray absorption spectroscopy reveals irreversible surface oxidation as the primary stability-limiting factor. This stepwise strategy provides a general framework for supported catalysts that maximize precious metal utilization via sub-nanometer dispersion.

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

Effects of Ink Formulation on the Structure and Performance of PGM-Free Catalyst Layer in PEMFCs

Platinum group metal (PGM) catalysts are the major electrocatalysts for oxygen reduction reaction (ORR) in the polymer electrolyte membrane fuel cells (PEMFCs). The cost becomes unaffordable if the PEMFC is in massive application. The PGM-Free catalyst shows very promising activity in rotation disk electrode (RDE) testing. The half-wave potential could reach 0.91 V versus standard hydrogen electrode (SHE). However, in a membrane electrode assembly (MEA), the performance of PGM-Free catalysts is not good enough to replace the PGM catalysts. Since the PGM-free catalysts are so different from the PGM catalysts in terms of catalytic activity, stability, surface conditions, particle size, etc., the fabrication of PGM-Free catalyst MEA cannot simply copy the method of making PGM MEA. Here we proposed a novel method of fabricating PGM-Free catalyst MEA, so that the intrinsic catalyst activity from RDE can be translated into MEA performance. The method is based on the catalyst coated membrane (CCM) method using optimized ionomer to carbon (I/C) ratio and solvent mixture of catalyst ink. Using this method, the PGM-free catalyst MEA achieved the current density 44.9 mA cm -2 at 0.9 V iR-free in H 2 /O 2 and 150 mA cm -2 at 0.8 V in H 2 /air, which surpassed the performance targets of US Department of Energy (DOE)for PGM-Free catalyst MEA. The property (solvent composition, dispersion of catalyst and ionomer in an ink), structure (pore structure) and the MEA performance have been characterized using, mercury intrusion porosimetry (MIP), MEA testing. A property-structure-performance relationship has been established.

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

Advanced Modeling of Plasma-based Particle Accelerators

Plasma-based acceleration (PBA) driven by an intense laser (LWFA) or particle beam (PWFA) can produce ultra-high accelerating fields in excess of a GV/cm. PBA could substantially reduce the size and cost of future linear collider facilities if deployed successfully. PBA enables compact tabletop accelerators that can provide lower energy GeV-class beams in a laboratory setting. PBA enables high quality beam generation suitable for x-ray free electron lasers (XFEL) via controllable methods of self-injection. Challenges in modeling and optimization of multi-stage PBA motivate the need for exascale computing and state-of-the-art PIC codes.

43 PARTICLE ACCELERATORS↗
NASA NTRS1974

Optimal aeroelastic design of an oblique wing structure

A procedure is presented for determining the optimal cover panel thickness of a wing structure to meet specified strength and static aeroelastic divergence requirements for minimum weight. Efficient reanalysis techniques using discrete structural and aerodynamic methods are used in conjunction with redesign algorithms driven by optimality criteria. The optimality conditions for the divergence constraint are established, and expressions are obtained for derivatives of the dynamic pressure at divergence with respect to design variables. The procedure is applied to an oblique wing aircraft where strength and stiffness are critical design considerations for sizing the cover thickness of the wing structure.

Gwin, L. B.↗
DOE OSTI2024

Optimized In Situ Doping Strategy Stabling Single-Crystal Ultrahigh-Nickel Layered Cathode Materials

Single-crystal Ni-rich cathodes offer promising prospects in mitigating intergranular microcracks and side reaction issues commonly encountered in conventional polycrystalline cathodes. However, the utilization of micrometer-sized single-crystal particles has raised concerns about sluggish Li + diffusion kinetics and unfavorable structural degradation, particularly in high Ni content cathodes. Herein, we present an innovative in situ doping strategy to regulate the dominant growth of characteristic planes in the singlecrystal precursor, leading to enhanced mechanical properties and effectively tackling the challenges posed by ultrahigh-nickel layered cathodes. Compared with the traditional dry-doping method, our in situ doping approach possesses a more homogeneous and consistent modifying effect from the inside out, ensuring the uniform distribution of doping ions with large radius (Nb, Zr, W, etc). Further, this mitigates the generally unsatisfactory substitution effect, thereby minimizing undesirable coating layers induced by different solubilities during the calcination process. Additionally, the uniformly dispersed ions from this in situ doping are beneficial for alleviating the two-phase coexistence of H2/H3 and optimizing the Li + concentration gradient during cycling, thus inhibiting the formation of intragranular cracks and interfacial deterioration. Consequently, the in situ doped cathodes demonstrate exceptional cycle retention and rate performance under various harsh testing conditions. Our optimized in situ doping strategy not only expands the application prospects of elemental doping but also offers a promising research direction for developing high-energy-density single-crystal cathodes with extended lifetime.

25 ENERGY STORAGE↗
DOE OSTI2026

Multi-modal characterization of nitrate reduction nano-catalysts with periodic strain distribution

Strain engineering serves as a pivotal strategy to optimize catalytic activity in electrocatalysis. However, the catalyst sizes under industrial conditions are usually large and even beyond nanometer regime. The critical methodological limitations on strain imaging of such catalysts with both large field of view and high spatial resolution obscure the mechanistic understanding of strain-performance correlations. Here, we present an optimized four-dimensional scanning transmission electron microscopy (4D-STEM) method to acquire strain mapping of both bulk and surface across particles up to 500 nm with 0.6 nm spatial resolution and 0.55% precision. We observe the ripple-like periodic strain coupled with elemental fluctuations inside a perovskite-type hydroxide CuCoSn(OH) 6 and find it correlated to electrocatalytic nitrate reduction (NO 3 – RR) absorption energy to achieve the 92.6% Faradaic efficiency and long-term test over 1000 h at membrane electrode assembly (MEA) for ammonia electrosynthesis. This universal framework design offers a practical method that not only develops an advanced measurement combining multi-modal characterization techniques but also reveals the intrinsic structure-property constitutive law of industry-level catalysts.

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

Benchmarking Quantum Simulators Using Ergodic Quantum Dynamics

Here we propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced spatiotemporal control. Our protocol relies on universal fluctuations emerging from generic Hamiltonian dynamics, which we discover in the present work. It does not require fine-tuned control over state preparation, quantum evolution, or readout capability, while achieving near optimal sample complexity: a percent-level precision is obtained with ~103 measurements, independent of system size. Furthermore, the accuracy of our fidelity estimation improves exponentially with increasing system size. We numerically demonstrate our protocol in a variety of quantum simulator platforms, including quantum gas microscopes, trapped ions, and Rydberg atom arrays. We discuss applications of our method for tasks such as multiparameter estimation of quantum states and processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗