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 1,027 records · Page 57

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↗

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

The described investigation is concerned with he development of a finite element scheme which can be used in a study of the acoustics of aircraft-engine ducts. In the absence of suitable variational principles for acoustic fluctuations within an aircraft fan engine, an acoustic analysis must proceed directly from the differential equations which describe compressible flow. The derived equations cannot be solved algebraically. The numerical technique used for solving them makes use of a linear rectangular element of a type considered by Zienkiewicz (1971). Attention is given to aspects of element derivation, the global matrix assembly, the solution of the matrix equation, questions of acoustic attenuation, and illustrations of the potential of the current model in duct optimization.

Abrahamson, A. L.↗

GCM Hindcasts for SST Forced Climate Variability over Agriculturally Intensive Regions

The ability to forecast seasonal climate is of great practical interest. One of the most obvious benefits would be agriculture, for which various preparations (planting, machinery, irrigation, manpower) would be enabled. The expectation of being able to make such forecasts far enough in advance (on the order of 9 months) hinges on components of the system with the longest persistence or predictability. The mixed results of El Nino forecasts has raised the hope that tropical Pacific sea surface temperatures (SST) fall into this category. For agriculturally-relevant forecasts to be made, and utilized, requires several conditions. The SST in the regions that affect agricultural areas must be forecast successfully, many months in advance. The climate response to such sea surface temperatures must then be ascertained, either through the use of historical empirical studies or models (e.g., GCMS). For practical applications, the agricultural production must be strongly influenced by climate, and farmers on either the local level or through commercial concerns must be able to adjust to using such forecasts. In a continuing series of papers, we will explore each of these components. This article concerns the question of utilizing SST to forecast the climate in several regions of agricultural production. We optimize the possibility of doing so successfully by using observed SST in a hindcast mode (i.e., a perfect forecast), and we also use the globally observed values (rather than just those from the tropical Pacific, for which predictability has been shown). This then is the ideal situation; in subsequent papers we will explore degrading the results by using only tropical Pacific SSTs, and then using only

Druyan, Leonard M.↗

Space Weather Investigation Frontier (SWIFT) Mission Concept: Continuous Distributed Observations of Geo-Effective, Heliospheric Structures from the Vantage Points of Sun-Earth L1 and Sub- L1

Continuous, in-situ, multi-point observations along the Sun-Earth line at and inside the Lagrange point L1 (sub-L1) will enable a better understanding of the three-dimensional structure and temporal evolution of heliospheric structures that drive terrestrial space weather. The proposed SWIFT (Space Weather Investigation Frontier) mission will use a new solar sail propulsion system developed by NASA to enable a suite of science instruments onboard a smallsat to maintain observations along the Sun-Earth line, sub-L1, for extended periods. Three identically instrumented small spacecraft at L1 will fly in concert with the sailcraft at sub-L1 in an optimized tetrahedron constellation, covering distances between 10 to 100’s of Earth radii. This viewing geometry will enable scientists to distinguish between local and global processes driving space weather by revealing the spatial characteristics, temporal evolution, and geo-effectiveness of small- to meso-scale solar wind structures and substructures of macro-scale structures, such as interplanetary coronal mass ejections (ICMEs) and stream interaction regions (SIRs). In addition, real time measurements of earth-bound heliospheric structures from sub-L1 will improve our current forecasting lead-times by up to 35 percent. This paper will provide an overview of the proposed SWIFT science and mission concept.

Spacecraft Propulsion↗

Space Weather Investigation Frontier (SWIFT)

Continuous, in-situ, multi-point observations along the Sun-Earth line at and inside the Lagrange point L1 (sub-L1) will enable a better understanding of the three-dimensional structure and temporal evolution of heliospheric structures that drive terrestrial space weather. The proposed SWIFT (Space Weather Investigation Frontier) mission will use a new solar sail propulsion system developed by NASA to enable a suite of science instruments onboard a smallsat to maintain observations along the Sun-Earth line, sub-L1, for extended periods. Three identically instrumented small spacecraft at L1 will fly in concert with the sailcraft at sub-L1in an optimized tetrahedron constellation, covering distances between 10 to 100’s of Earth radii. This viewing geometry will enable scientists to distinguish between local and global processes driving space weather by revealing the spatial characteristics, temporal evolution, and geo-effectiveness of small-to meso-scale solar wind structures and substructures of macro-scale structures, such as interplanetary coronal mass ejections (ICMEs) and stream interaction regions (SIRs). In addition, real time measurements of earth-bound heliospheric structures from sub-L1 will improve our current forecasting lead-times by up to 35 percent. This paper will provide an overview of the proposed SWIFT science and mission concept.

Heliophysics↗

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↗

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↗

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↗

Combined control-structure optimization

An approach for combined control-structure optimization keyed to enhancing early design trade-offs is outlined and illustrated by numerical examples. The approach employs a homotopic strategy and appears to be effective for generating families of designs that can be used in these early trade studies. Analytical results were obtained for classes of structure/control objectives with linear quadratic Gaussian (LQG) and linear quadratic regulator (LQR) costs. For these, researchers demonstrated that global optima can be computed for small values of the homotopy parameter. Conditions for local optima along the homotopy path were also given. Details of two numerical examples employing the LQR control cost were given showing variations of the optimal design variables along the homotopy path. The results of the second example suggest that introducing a second homotopy parameter relating the two parts of the control index in the LQG/LQR formulation might serve to enlarge the family of Pareto optima, but its effect on modifying the optimal structural shapes may be analogous to the original parameter lambda.

Salama, M.↗

Comparing Sensor Fusion and Multimodal Chemometric Models for Monitoring U(VI) in Complex Environments Representative of Irradiated Nuclear Fuel

Optical sensors and chemometric models were leveraged for the quantification of uranium(VI) (0–100 μg mL –1 ), europium (0–150 μg mL –1 ), samarium (0–250 μg mL –1 ), praseodymium (0–350 μg mL –1 ), neodymium (0–1000 μg mL –1 ), and HNO 3 (2–4 M) with varying corrosion product (iron, nickel, and chromium) levels using laser fluorescence, Raman scattering, and ultraviolet–visible–near-infrared absorption spectra. In this paper, an efficient approach to developing and evaluating tens of thousands of partial least-squares regression (PLSR) models, built from fused optical spectra or multimodal acquisitions, is discussed. Each PLSR model was optimized with unique preprocessing combinations, and features were selected using genetic algorithm filters. The 7-factor D-optimal design training set contained just 55 samples to minimize the number of samples. The performance of PLSR models was evaluated by using an automated latent variable selection script. PLS1 regression models tailored to each species outperformed a global PLS2 model. PLS1 models built using fused spectra data and a multimodal (i.e., analyzed separately) approach yielded similar information, resulting in percent root-mean-square error of prediction values of 0.9–5.7% for the seven factors. Further, the optical techniques and data processing strategies established in this study allow for the direct analysis of numerous species without measuring luminescence lifetimes or relying on a standard addition approach, making it optimal for near-real-time, in situ measurements. Nuclear reactor modeling helped bound training set conditions and identified elemental ratios of lanthanide fission products to characterize the burnup of irradiated nuclear fuel. Leveraging fluorescence, spectrophotometry, experimental design, and chemometrics can enable the remote quantification and characterization of complex systems with numerous species, monitor system performance, help identify the source of materials, and enable rapid high-throughput experiments in a variety of industrial processes and fundamental studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Achieving Breakthroughs in Global Hydrologic Science by Unlocking the Power of Multisensor, Multidisciplinary Earth Observations

Over the last half century, remote sensing has transformed hydrologic science. Whereas early efforts were devoted to observation of discrete variables, we now consider spaceborne missions dedicated to interlinked global hydrologic processes.Furthermore, cloud computing and computational techniquesare accelerating analyses of these data. How will the hydrologic community use these new resources to better understand the world’s water and relatedchallenges facing society? In this Commentary, we suggest that optimizing the benefits of remote sensing for advancing hydrologic research will happen byintegratingmultidisciplinary and multisensor data, leveraging commercial satellite measurements, and employingdata assimilation, cloud computing, and machine learning.We provide several recommendations to these ends. Plain Language Summary Observations from satellites have transformed hydrologic science. Early efforts, five decades ago, mapped attributes like snow cover, rainfall, topography, and vegetation, but now we consider new missions specifically designed to study global hydrologic processes. We also take advantageof new technologies like cloud computing and artificial intelligence. We describe strategiesfor maximizing the benefits of remote sensing for hydrology, encouraging research across disciplines using multiple sensors, using new commercially available satellites, and combining remote sensing measurements with hydrologic models.

Michael Durand↗

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗

Charged lepton flavor violation at the EIC

We present a comprehensive analysis of the potential sensitivity of the Electron-Ion Collider (EIC) to charged lepton flavor violation (CLFV) in the channel ep → τX, within the model-independent framework of the Standard Model Effective Field Theory (SMEFT). We compute the relevant cross sections to leading order in QCD and electroweak corrections and perform simulations of signal and SM background events in various τ decay channels, suggesting simple cuts to enhance the associated estimated efficiencies. To assess the discovery potential of the EIC in τ-e transitions, we study the sensitivity of other probes of this physics across a broad range of energy scales, from pp → eτX at the Large Hadron Collider to decays of B mesons and τ leptons, such as τ → eγ, τ → eℓ+ℓ–, and crucially the hadronic modes τ → eY with Yϵ {π, K, ππ, Kπ, …}. We find that electroweak dipole and four-fermion semi-leptonic operators involving light quarks are already strongly constrained by τ decays, while operators involving the c and b quarks present more promising discovery potential for the EIC. An analysis of three models of leptoquarks confirms the expectations based on the SMEFT results. We also identify future directions needed to maximize the reach of the EIC in CLFV searches: these include an optimization of the τ tagger in hadronic channels, an exploration of background suppression through tagging b and c jets in the final state, and a global fit by turning on all SMEFT couplings, which will likely reveal new discovery windows for the EIC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Robust cooperative control strategy for a platoon of connected and autonomous vehicles against sensor errors and control errors simultaneously in a real-world driving environment

In a real-world driving environment, a platoon of connected and autonomous vehicles (CAVs) is subject to many internal and external disturbances, resulting in uncertain vehicle dynamics. In general, the disturbances can be categorized into two types: disturbances due to vehicle sensor errors (e.g., GPS error) and disturbances due to vehicle control errors (e.g., actuator delay). In the literature, many control strategies have been proposed to improve the robustness of the CAV platoon against uncertain vehicle dynamics induced by these disturbances. However, most of these strategies only consider one type of disturbance and cannot tackle both types of disturbances simultaneously. Furthermore, they are designed to maximize the benefits of each vehicle in the platoon independently, which can deteriorate the performance of the platoon. Here, to address these problems, this study proposes a robust cooperative control (RCC) strategy to maneuver the vehicles in the platoon cooperatively to counteract the impacts of both types of disturbances. The RCC strategy is developed based on a minimax problem, where the maximization subproblem seeks to find the worst inputs for the uncertainty terms in the vehicle dynamics equation to minimize the platoon performance, while the minimization subproblem seeks to find the optimal control decisions for all subsequent vehicles to maximize the platoon performance in the worst case. To solve the minimax problem, this study proposes a globally convergent solution algorithm. It can solve the minimax problem very efficiently to enable real time deployment of the RCC strategy. Numerical application indicates that compared to the existing methods, the RCC strategy can dramatically improve the robustness of the CAV platoon against the uncertain vehicle dynamics induced by both vehicle state detection errors and vehicle control errors. Therefore, it can maneuver the CAV platoon safely and efficiently in a real-world driving environment.

33 ADVANCED PROPULSION SYSTEMS↗

A circuit-generated quantum subspace algorithm for the variational quantum eigensolver

Recent research has shown that wavefunction evolution in real and imaginary time can generate quantum subspaces with significant utility for obtaining accurate ground state energies. Inspired by these methods, we propose combining quantum subspace techniques with the variational quantum eigensolver (VQE). In our approach, the parameterized quantum circuit is divided into a series of smaller subcircuits. The sequential application of these subcircuits to an initial state generates a set of wavefunctions that we use as a quantum subspace to obtain high-accuracy groundstate energies. We call this technique the circuit subspace variational quantum eigensolver (CSVQE) algorithm. By benchmarking CSVQE on a range of quantum chemistry problems, we show that it can achieve significant error reduction in the best case compared to conventional VQE, particularly for poorly optimized circuits, greatly improving convergence rates. Furthermore, we demonstrate that when applied to circuits trapped at local minima, CSVQE can produce energies close to the global minimum of the energy landscape, making it a potentially powerful tool for diagnosing local minima.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local chromatic correction optics for Future Circular Collider e + e −

Local chromatic correction optics are proposed for the Future Circular e + − e − Collider. These new optics assume an identical layout of the magnets at all operation energies and include a new design for the arcs, for the straight sections, and for the final focus. The arcs design is a step forward from the classic FODO lattice that achieves near cancellation up to the fourth order of chromatic and geometric aberrations. Straight sections and the final focus benefit from the application of transparency conditions and are thus minimally impacting the global beam dynamics. The final focus design includes both vertical and horizontal chromaticity correction sections, crab sextupoles, and additional specific sextupoles and decapoles for the optimization of the dynamics off-energy. Decapoles are very effective in mitigating the reduction of dynamic aperture due to synchrotron radiation. Nonlinear magnet settings are further optimized with multiobjective algorithms. The overall properties of the new optics proposed are analyzed in the presence of errors and compared to other optics designs.

43 PARTICLE ACCELERATORS↗

Multi-Level Adaptive Techniques (MLAT) for singular-perturbation problems

The multilevel (multigrid) adaptive technique, a general strategy of solving continuous problems by cycling between coarser and finer levels of discretization is described. It provides very fast general solvers, together with adaptive, nearly optimal discretization schemes. In the process, boundary layers are automatically either resolved or skipped, depending on a control function which expresses the computational goal. The global error decreases exponentially as a function of the overall computational work, in a uniform rate independent of the magnitude of the singular-perturbation terms. The key is high-order uniformly stable difference equations, and uniformly smoothing relaxation schemes.

Brandt, A.↗