Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dashboard for Marine Energy Site Assessment and Monitoring

The marine energy (ME) industry presently relies upon fragmented site assessment solutions that require high resource expenditure for deployment at each site and do not leverage the wealth of readily available tools and information. A wave energy resource assessment dashboard, currently in development, will substantially improve siting, permitting, operations, and maintenance of ME projects by providing an integrated solution that is a one-stop-shop for a developer’s needs. The Site Energy Assessment and MOnitoring Dashboard (SEAMOD) will be of commercial interest to anyone seeking to deploy an ME project and is easily expandable to include tidal and wind energy site assessments. The integrated dashboard is being developed using state-of-the-art database and cloud computing methods and data-assimilative modeling tools that can be coupled with low-cost, rapidly deployable wave buoys and environmental sensing hardware. The combined software and hardware dashboard will reduce wave energy site characterization and wave climate monitoring costs by more than 60 percent and provide assessments that meet international industry standards. To realize a thriving global ME industry, the physical environment at a potential deployment site must be understood, not only for resource characterization, but also for optimization of device and power conversion performance. SEAMOD directly addresses these needs with a commercially marketable product. SEAMOD is a low-cost solution that provides comprehensive ME resource assessments, baseline environmental monitoring, and offshore characterizations required for successful ME development. The key technical objectives for Phase I were a series of software development goals, which when implemented with monitoring solutions, produced an initial proof-of-concept low-cost wave energy resources dashboard. In Phase II, the development of the prototype SEAMOD continued. The basic framework employed was the development of a revised dashboard and monitoring tool customized for ME applications by focusing on IEC site assessment and method requirements. Development was focused on the integration of full hindcast metocean products to provide hindcast resource characterization and environmental information. The final integrated dashboard provides a low-cost solution that delivers comprehensive, scalable, industry-standard energy resource assessments and offshore characterizations required for successful ME development. The integrated dashboard offers visibility of the most recent site modeling, measurements, and historical data. The application and integration of consensus-based standards for wave energy resource assessment, as determined by the International Electrotechnical Commission (IEC), are crucial for the impact and value of SEAMOD. SEAMOD includes monthly, seasonal, and yearly statistics, as well as the total 30-year record, offering temporal resolution of the IEC parameters to aid potential developers in determining the available wave energy resources in their area of interest.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Optimizing Carbon Capture, Transport, and Storage: Overcoming Challenges with Machine Learning and Cost-Benefit Analysis

Carbon Capture, Utilization, and Storage (CCUS) is a critical strategy for reducing CO₂ emissions and mitigating climate change. However, its widespread deployment faces numerous challenges across the capture, transport, and storage phases. These challenges include the technical complexity of predicting subsurface behaviors during CO₂ injection, ensuring long-term storage integrity, optimizing transportation networks, and balancing the economic and environmental trade-offs. Addressing these issues requires an integrated approach combining advanced subsurface modeling with system-level analyses to assess costs, risks, and benefits. This presentation provides an overview of studies conducted by the National Energy Technology Laboratory (NETL) to tackle these challenges. NETL’s efforts encompass cutting-edge research in subsurface fluid behavior machine learning predictions, alongside the development of innovative tools for system optimization and economic evaluation. By bridging technical expertise and strategic analysis, NETL aims to advance the deployment of CCUS technologies to support global decarbonization efforts. Presented at the Carnegie Mellon University CEE IESS Student Seminar October 4, 2024.

Shih, Chung Yan↗
DOE OSTI2020

Bayesian Optimization of a Free-Electron Laser

The Linac Coherent Light Source X-ray free-electron laser is a complex scienti c apparatus which changes con gurations multiple times per day, necessitating fast tuning strategies to reduce setup time for successive experiments. To this end, we employ a Bayesian approach to maximizing X- ray laser pulse energy by controlling groups of quadrupole magnets. A Gaussian process model provides probabilistic predictions for the machine response with respect to control parameters, enabling a balance of exploration and exploitation in the search for the global optimum. We show that the model parameters can be learned from archived scans, and correlations between devices can be extracted from the beam transport. Furthermore, the result is a sample-e cient optimization routine, combining both historical data and knowledge of accelerator physics to signi cantly outperform existing optimizers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗
DOE OSTI2025

A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations

Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.

Silvestri, Simone [Massachusetts Institute of Tech↗
DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sub-quantum Fisher information

We report the quantum Fisher information (QFI) plays a crucial role in quantum information theory and in many practical applications such as quantum metrology. However, computing the QFI is generally a computationally demanding task. In this work we analyze a lower bound on the QFI which we call the sub-quantum Fisher information (sub-QFI). The bound can be efficiently estimated on a quantum computer for an n-qubit state using 2n qubits. The sub-QFI is based on the super-fidelity, an upper bound on Uhlmann's fidelity. We analyze the sub-QFI in the context of unitary families, where we derive several crucial properties including its geometrical interpretation. In particular, we prove that the QFI and the sub-QFI are maximized for the same optimal state, which implies that the sub-QFI is faithful to the QFI in the sense that both quantities share the same global extrema. Based on this faithfulness, the sub-QFI acts as an efficiently computable surrogate for the QFI for quantum sensing and quantum metrology applications. Finally, we provide additional meaning to the sub-QFI as a measure of coherence, asymmetry, and purity loss.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Unraveling the Origin of Glassy Thermal Transport in Medium-Entropy Semiconductors: From Nanoscale Phase Segregation to Atomic-Scale Lattice Distortion

Thermoelectric (TE) materials can directly and reversibly convert heat into electricity, offering a promising pathway for reshaping the global energy landscape. The performance of thermoelectric materials is evaluated using the figure of merit, zT = S 2 σT/κ, which can be optimized by improving Seebeck coefficient (S) and electrical conductivity (σ) while suppressing the thermal conductivity (κ) at the operating temperature (T). Here, we introduce a medium-entropy metal chalcogenide that exhibits phonon-glass electron-crystal behavior. The conventional perspective hypothesizes that the suppression of thermal transport arises solely from disordered atomic occupation, which enhances alloy scattering mechanisms for phonon transport. However, emerging evidence suggests the presence of structural and chemical inhomogeneity at multiple length scales in entropy-engineered materials, a phenomenon that has rarely been investigated in detail.

Liu, Yukun [Northwestern Univ]↗
DOE OSTI2023

Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations

The Keller–Segel–Navier–Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L 1 (Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller–Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data optimization for large batch distributed training of deep neural networks

Distributed training in deep learning (DL) is common practice as data and models grow. The current practice for distributed training of deep neural networks faces the challenges of communication bottlenecks when operating at scale, and model accuracy deterioration with an increase in global batch size. Present solutions focus on improving message exchange efficiency as well as implementing techniques to tweak batch sizes and models in the training process. The loss of training accuracy typically happens because the loss function gets trapped in a local minima. We observe that the loss landscape minimization is shaped by both the model and training data and propose a data optimization approach that utilizes machine learning to implicitly smooth out the loss landscape resulting in fewer local minima. Our approach filters out data points which are less important to feature learning, enabling us to speed up the training of models on larger batch sizes to improved accuracy.

Gahlot, Shubhankar↗
DOE OSTI2024

Membrane Contactors for Ammonia Recovery from Anaerobic Digester Centrate: Pretreatment and Process Optimization

Haber-Bosch process allows for the production of modern fertilizers and is crucial for meeting increasing demands for agricultural production. The process requires large amounts of natural gas and contributes to global warming. An alternative to the Haber-Bosch process utilizes membrane contactors to recover nitrogen from the anaerobic digester centrate. In previous studies, high recoveries have been achieved, but membrane fouling decreased the performance and required maintenance to clean the membranes. Furthermore, this study investigated the effect of a settling and ultrafiltration system to pretreat centrate from an anaerobic digester to prevent fouling of the membrane contactor system. The system achieved high recoveries of over 90% for 10 cycles without any performance decline. Tests with increasing distillate concentrations of ammonium sulfate up to 165,000 mg/L-N could not identify a significant decline in membrane performance either. This allows for concentration up to crystallization in a single stage.

09 BIOMASS FUELS↗
DOE OSTI2022

Itaconic acid production is regulated by LaeA in Aspergillus pseudoterreus

The global regulator LaeA controls secondary metabolism in diverse Aspergillus species. Here we explored its role in regulation of itaconic acid production in Aspergillus pseudoterreus. To understand its role in regulating metabolism, we deleted and overexpressed laeA, and assessed the transcriptome, proteome, and secreted metabolome prior to and during initiation of phosphate limitation induced itaconic acid production. We found that secondary metabolite clusters, including the itaconic acid biosynthetic gene cluster, are regulated by laeA and that laeA is required for high yield production of itaconic acid. Overexpression of laeA improves itaconic acid yield at the expense of biomass by increasing the expression of key biosynthetic pathway enzymes and attenuating the expression of genes involved in phosphate acquisition and scavenging. Increased yield was observed in optimized conditions as well as conditions containing excess nutrients that may be present in inexpensive sugar containing feedstocks such as excess phosphate or complex nutrient sources. This suggests that global regulators of metabolism may be useful targets for engineering metabolic flux that is robust to environmental heterogeneity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Optimizing the electron transport chain to sustainably improve photosynthesis

Abstract Genetically improving photosynthesis is a key strategy to boosting crop production to meet the rising demand for food and fuel by a rapidly growing global population in a warming climate. Many components of the photosynthetic apparatus have been targeted for genetic modification for improving photosynthesis. Successful translation of these modifications into increased plant productivity in fluctuating environments will depend on whether the electron transport chain (ETC) can support the increased electron transport rate without risking overreduction and photodamage. At present atmospheric conditions, the ETC appears suboptimal and will likely need to be modified to support proposed photosynthetic improvements and to maintain energy balance. Here, I derive photochemical equations to quantify the transport capacity and the corresponding reduction level based on the kinetics of redox reactions along the ETC. Using these theoretical equations and measurements from diverse C3/C4 species across environments, I identify several strategies that can simultaneously increase the transport capacity and decrease the reduction level of the ETC. These strategies include increasing the abundances of reaction centers, cytochrome b6f complexes, and mobile electron carriers, improving their redox kinetics, and decreasing the fraction of secondary quinone–nonreducing photosystem II reaction centers. I also shed light on several previously unexplained experimental findings regarding the physiological impacts of the abundances of the cytochrome b6f complex and plastoquinone. The model developed, and the insights generated from it facilitate the development of sustainable photosynthetic systems for greater crop yields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗
DOE OSTI2023

Learning Provably Stable Local Volt/Var Controllers for Efficient Network Operation

Here this paper develops a data-driven framework to synthesize local Volt/Var control strategies for distributed energy resources (DERs) in power distribution grids (DGs). Aiming to improve DG operational efficiency, as quantified by a generic optimal reactive power flow (ORPF) problem, we propose a two-stage approach. The first stage involves learning the manifold of optimal operating points determined by an ORPF instance. To synthesize local Volt/Var controllers, the learning task is partitioned into learning local surrogates (one per DER) of the optimal manifold with voltage input and reactive power output. Since these surrogates characterize efficient DG operating points, in the second stage, we develop local control schemes that steer the DG to these operating points. We identify the conditions on the surrogates and control parameters to ensure that the locally acting controllers collectively converge, in a global asymptotic sense, to a DG operating point agreeing with the local surrogates. We use neural networks to model the surrogates and enforce the identified conditions in the training phase. AC power flow simulations on the IEEE 37-bus network empirically bolster the theoretical stability guarantees obtained under linearized power flow assumptions. The tests further highlight the optimality improvement compared to prevalent benchmark methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Harnessing the power of gradient-based simulations for multi-objective optimization in particle accelerators

Abstract Particle accelerator operation requires simultaneous optimization of multiple objectives. Multi-objective optimization (MOO) is particularly challenging due to trade-offs between the objectives. Evolutionary algorithms, such as genetic algorithms (GAs), have been leveraged for many optimization problems, however, they do not apply to complex control problems by design. This paper demonstrates the power of differentiability for solving MOO problems in particle accelerators using a deep differentiable reinforcement learning (DDRL) algorithm. We compare the DDRL algorithm with model-free reinforcement learning (MFRL), GA, and Bayesian optimization (BO) for simultaneous optimization of heat load and trip rates in the continuous electron beam accelerator facility. The underlying problem enforces strict constraints on both individual states and actions as well as cumulative (global) constraints on energy requirements of the beam. Using historical accelerator data, we develop a physics-based surrogate model which is differentiable and allows for back-propagation of gradients. The results are evaluated in the form of a Pareto-front with two objectives. We show that the DDRL outperforms MFRL, BO, and GA on high dimensional problems.

43 PARTICLE ACCELERATORS↗