Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained 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

Third-Order Møller–Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller–Plesset (MP) perturbation theory is severely constrained by the use of Hartree–Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). As such, MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N 6 ) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superconvergence of Online Optimization for Model Predictive Control

We develop a one-Newton-step-per-horizon, online, lag-L, model predictive control (MPC) algorithm for solving discrete-time, equality-constrained, nonlinear dynamic programs. Based on recent sensitivity analysis results for the target problems class, we prove that the approach exhibits a behavior that we call superconvergence; that is, the tracking error with respect to the full horizon solution is not only stable for successive horizon shifts, but also decreases with increasing shift order to a minimum value that decays exponentially in the length of the receding horizon. The key analytical step is the decomposition of the one-step error recursion of our algorithm into algorithmic error and perturbation error. We show that the perturbation error decays exponentially with the lag between two consecutive receding horizons, while the algorithmic error, determined by Newton’s method, achieves quadratic convergence instead. Overall this approach induces our local exponential convergence result in terms of the receding horizon length for suitable values of L. In conclusion, numerical experiments validate our theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Budget Constrained Machine Learning for Early Prediction of Adverse Outcomes for COVID-19 Patients

Background: Machine learning (ML) based risk stratification models of Electronic Health records (EHR) data may help to optimize treatment of COVID-19 patients, but are often limited by their lack of clinical interpretability and cost of laboratory tests. We develop a ML based tool for predicting adverse outcomes based on EHR data to optimize clinical utility under a given cost structure. This cohort study was performed using deidentified EHR data from COVID-19 patients from ProMedica Healthcare in northwest Ohio and southeastern Michigan. Methods: We tested performance of various ML approaches for predicting either increasing ventilatory support or mortality and the set of model features under a budget constraint was optimized via exhaustive search across all combinations of features. Results: The optimal sets of features for predicting ventilation under any budget constraint included demographics and comorbidities (DCM), basic metabolic panel (BMP), D-dimer, lactate dehydrogenase (LDH), erythrocyte sedimentation rate (ESR), CRP, brain natriuretic peptide (BNP), and procalcitonin and for mortality included DCM, BMP, complete blood count, D-dimer, LDH, CRP, BNP, procalcitonin and ferritin. Conclusions: This study presents a quick, accurate and cost-effective method to evaluate risk of deterioration for patients with SARS-CoV-2 infection at the time of clinical evaluation.

Cadena Pico, JoseE.↗

Assimilation of multiple datasets results in large differences in regional- to global-scale NEE and GPP budgets simulated by a terrestrial biosphere model

In spite of the importance of land ecosystems in offsetting carbon dioxide emissions released by anthropogenic activities into the atmosphere, the spatiotemporal dynamics of terrestrial carbon fluxes remain largely uncertain at regional to global scales. Over the past decade, data assimilation (DA) techniques have grown in importance for improving these fluxes simulated by terrestrial biosphere models (TBMs), by optimizing model parameter values while also pinpointing possible parameterization deficiencies. Although the joint assimilation of multiple data streams is expected to constrain a wider range of model processes, their actual benefits in terms of reduction in model uncertainty are still under-researched, also given the technical challenges. In this study, we investigated with a consistent DA framework and the ORCHIDEE-LMDz TBM–atmosphere model how the assimilation of different combinations of data streams may result in different regional to global carbon budgets. To do so, we performed comprehensive DA experiments where three datasets (in situ measurements of net carbon exchange and latent heat fluxes, spaceborne estimates of the normalized difference vegetation index, and atmospheric CO 2 concentration data measured at stations) were assimilated alone or simultaneously. We thus evaluated their complementarity and usefulness to constrain net and gross C land fluxes. We found that a major challenge in improving the spatial distribution of the land C sinks and sources with atmospheric CO 2 data relates to the correction of the soil carbon imbalance.

54 ENVIRONMENTAL SCIENCES↗

Cosmic-structure constraints from a one-degree microwave-background anisotropy experiment

A recent anisotropy experiment with a 30-arcmin beam is close to optimal for testing models in which cosmic structure arose from inflation-generated Gaussian density perturbations. Using a Bayesian analysis to constrain the amplitude of the perturbation spectrum, it is shown that adiabatic hot-dark-matter models are convincingly ruled out and cold-dark-matter models have anisotropies near the derived limits. Theories with broken scale invariance such as isocurvature baryon models are strongly constrained, as are models with extra power crafted to give significant large-scale structures and flows.

Bond, J. R.↗

Analytical solutions to constrained hypersonic flight trajectories

The flight trajectory of aerospace vehicles subject to a class of path constraints is considered. The constrained dynamics is shown to be a natural two-time-scale system. Asymptotic analytical solutions are obtained. Problems of trajectory optimization and guidance can be dramatically simplified with these solutions. Applications in trajectory design for an aerospace plane strongly support the theoretical development.

Lu, Ping↗

Analytical solutions to constrained hypersonic flight trajectories

The flight trajectory of aerospace vehicles subject to a class of path constraints is considered. The constrained dynamics is shown to be a natural two-time-scale system. Asymptotic analytical solutions are obtained. Problems of trajectory optimization and guidance can be dramatically simplified with these solutions. Applications in trajectory design for an aerospace plane strongly support the theoretical development.

Lu, Ping↗

A Globally Convergent Augmented Lagrangian Pattern Search Algorithm for Optimization with General Constraints and Simple Bounds

We give a pattern search adaptation of an augmented Lagrangian method due to Conn, Gould, and Toint. The algorithm proceeds by successive bound constrained minimization of an augmented Lagrangian. In the pattern search adaptation we solve this subproblem approximately using a bound constrained pattern search method. The stopping criterion proposed by Conn, Gould, and Toint for the solution of this subproblem requires explicit knowledge of derivatives. Such information is presumed absent in pattern search methods; however, we show how we can replace this with a stopping criterion based on the pattern size in a way that preserves the convergence properties of the original algorithm. In this way we proceed by successive, inexact, bound constrained minimization without knowing exactly how inexact the minimization is. So far as we know, this is the first provably convergent direct search method for general nonlinear programming.

Lewis, Robert Michael↗

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗

Optical and Mechanical Designs of the Multi-band SWIR Receiver for the Lunar Flashlight CubeSat Mission

Mapping and quantifying lunar water ice addresses one of NASA’s Strategic Knowledge Gaps to understand the lunar resource potential for future human exploration of the Moon. Lunar Flashlight is an innovative NASA CubeSat mission dedicated to mapping water ice in the permanently-shadowed and occasionally-sunlit regions in the vicinity of the lunar South Pole. Lunar Flashlight will acquire these measurements from lunar orbit using a multi-band laser reflectometer composed of an optical receiver aligned with four lasers emitting different wavelengths in the shortwave infrared spectral region between 1 µm and 2 µm. The receiver measures the laser radiance reflected from the lunar surface in each spectral band and continuum/absorption reflectance band ratios are then analyzed to quantify water ice concentration in the illuminated spot. The receiver utilizes a 70×70-mm, aluminum, off-axis paraboloidal mirror with a focal length of 70 mm, which collects the incoming light onto a single, 2 mm diameter InGaAs detector with a cutoff wavelength of 2.4 µm. We present the optical and mechanical designs of the receiver, including its optimization for rejection of solar stray-light from outside its intended field of view. This highly mass- and volume-constrained instrument payload will demonstrate several firsts, including being one of the first instruments onboard a CubeSat performing science measurements beyond low Earth orbit and the first planetary mission to use multi-band active reflectometry from orbit.

Sellar, R. Glenn↗

A Functional Interpolation Approach to Compute Period Orbits in the Circular Restricted Three-body Problem

In this paper, we develop a method to solve for periodic orbits, i.e., Lyapunov and Halo orbits, using a functional interpolation scheme called the Theory of Func- tional Connections (TFC). Using this technique, a periodic constraint is analyti- cally embedded into the TFC constrained expression. By doing this, the system of differential equations governing the three-body problem is transformed into an unconstrained optimization problem where simple numerical schemes can be used to find a solution, e.g., nonlinear least-squares is used. This allows for a simpler numerical implementation with comparable accuracy and speed to the traditional differential corrector method.

Mortari, Daniele↗

Broadband Predictions of Optimized Proprotors in Axial Forward Flight

A low-fidelity toolchain was used to predict broadband self-noise trends for three proprotors in axial flight: a baseline (C24ND) and two acoustically constrained proprotors (OPT-III and COPR-3). Rotor loads were predicted with the ANOPP Propeller Analysis System (PAS) and self-noise was predicted with the semiempirical method of Brooks, Pope, and Marcolini (BPM) implemented in the ANOPP2 Self-Noise Internal Functional Module (ASNIFM). Comparisons to experimental data revealed that trends for turbulent boundary layer trailing edge (TBL-TE) noise could be modeled across several flight conditions by increasing the boundary layer thicknesses via the trip setting. Since the Mach number range of the BPM method is exceeded in these predictions, a dependence of boundary layer displacement thickness on blade station Mach number was suggested as a possible reason for needing to model thicker boundary layers, suggesting that the TBL-TE model needs to be developed further. Bluntness vortex shedding noise (BVS) predictions required tuning the trailing edge thickness and trailing edge closure angle for each flight condition to match experimental trends, demonstrating that the BVS noise model is incomplete and that BVS noise may vary with the angle of attack. This study indicates that the BPM self-noise method needs to be improved, which will lead to more accurate broadband predictions.

Broadband self-noise↗

Broadband Predictions of Optimized Proprotors in Axial Forward Flight

A low-fidelity toolchain was used to predict broadband self-noise trends for three proprotors in axial flight: a baseline (C24ND) and two acoustically constrained proprotors (OPT-III and COPR-3). Rotor loads were predicted with the ANOPP Propeller Analysis System (PAS) and self-noise was predicted with the semiempirical method of Brooks, Pope, and Marcolini (BPM) implemented in the ANOPP2 Self-Noise Internal Functional Module (ASNIFM). Comparisons to experimental data revealed that trends for turbulent boundary layer trailing edge (TBL-TE) noise could be modeled across several flight conditions by increasing the boundary layer thicknesses via the trip setting. Since the Mach number range of the BPM method is exceeded in these predictions, a dependence of boundary layer displacement thickness on blade station Mach number was suggested as a possible reason for needing to model thicker boundary layers, suggesting that the TBL-TE model needs to be developed further. Bluntness vortex shedding noise (BVS) predictions required tuning the trailing edge thickness and trailing edge closure angle for each flight condition to match experimental trends, demonstrating that the BVS noise model is incomplete and that BVS noise may vary with the angle of attack. This study indicates that the BPM self-noise method needs to be improved, which will lead to more accurate broadband predictions.

Broadband self-noise↗

Tula: Optimizing Time, Cost, and Generalization in Distributed Large-Batch Training

Distributed training increases the number of batches processed per iteration either by scaling-out (adding more nodes) or scaling-up (increasing the batch-size). However, the largest configuration does not necessarily yield the best performance. Horizontal scaling introduces additional communication overhead, while vertical scaling is constrained by computation cost and device memory limits. Thus, simply increasing the batch-size leads to diminishing returns: training time and cost decrease initially but eventually plateaus, creating a knee-point in the time/cost vs. batch-size pareto curve. The optimal batch-size therefore depends on the underlying model, data and available compute resources. Large batches also suffer from worse model quality due to the well-known “generalization gap”. In this paper, we present Tula, an online service that automatically optimizes time, cost, and convergence quality for large-batch training of convolutional models. It combines parallel-systems modeling with statistical performance prediction to identify the optimal batchsize. Tula predicts training time and cost within 7.5−14% error across multiple models, and achieves up to 20× overall speedup and improves test accuracy by ≈9% on average over standard large-batch training on various vision tasks, thus successfully mitigating the generalization gap and accelerating training at the same time.

Tyagi, Sahil [ORNL] (ORCID:0009000783144745)↗

Leading-edge slat optimization for maximum airfoil lift

A numerical procedure for determining the position (horizontal location, vertical location, and deflection) of a leading edge slat that maximizes the lift of multielement airfoils is presented. The structure of the flow field is calculated by iteratively coupling potential flow and boundary layer analysis. This aerodynamic calculation is combined with a constrained function minimization analysis to determine the position of a leading edge slat so that the suction peak on the nose of the main airfoil is minized. The slat position is constrained by the numerical procedure to ensure an attached boundary layer on the upper surface of the slat and to ensure negligible interaction between the slat wake and the boundary layer on the upper surface of the main airfoil. The highest angle attack at which this optimized slat position can maintain attached flow on the main airfoil defines the optimum slat position for maximum lift. The design method is demonstrated for an airfoil equipped with a leading-edge slat and a trailing edge, single-slotted flap. The theoretical results are compared with experimental data, obtained in the Ames 40 by 80 Foot Wind Tunnel, to verify experimentally the predicted slat position for maximum lift. The experimentally optimized slat position is in good agreement with the theoretical prediction, indicating that the theoretical procedure is a feasible design method.

Olson, L. E.↗

Augmenting Parametric Optimal Ascent Trajectory Modeling with Graph Theory

It has been well documented that decisions made in the early stages of Conceptual and Pre-Conceptual design commit up to 80% of total Life-Cycle Cost (LCC) while engineers know the least about the product they are designing [1]. Once within Preliminary and Detailed design however, making changes to the design becomes far more difficult to enact in both cost and schedule. Primarily this has been due to a lack of detailed data usually uncovered later during the Preliminary and Detailed design phases. In our current budget-constrained environment, making decisions within Conceptual and Pre-Conceptual design which minimize LCC while meeting requirements is paramount to a program's success. Within the arena of launch vehicle design, optimizing the ascent trajectory is critical for minimizing the costs present within such concerns as propellant, aerodynamic, aeroheating, and acceleration loads while meeting requirements such as payload delivered to a desired orbit. In order to optimize the vehicle design its constraints and requirements must be known, however as the design cycle proceeds it is all but inevitable that the conditions will change. Upon that change, the previously optimized trajectory may no longer be optimal, or meet design requirements. The current paradigm for adjusting to these updates is generating point solutions for every change in the design's requirements [2]. This can be a tedious, time-consuming task as changes in virtually any piece of a launch vehicle's design can have a disproportionately large effect on the ascent trajectory, as the solution space of the trajectory optimization problem is both non-linear and multimodal [3]. In addition, an industry standard tool, Program to Optimize Simulated Trajectories (POST), requires an expert analyst to produce simulated trajectories that are feasible and optimal [4]. In a previous publication the authors presented a method for combatting these challenges [5]. In order to bring more detailed information into Conceptual and Pre-Conceptual design, knowledge of the effects originating from changes to the vehicle must be calculated. In order to do this, a model capable of quantitatively describing any vehicle within the entire design space under consideration must be constructed. This model must be based upon analysis of acceptable fidelity, which in this work comes from POST. Design space interrogation can be achieved with surrogate modeling, a parametric, polynomial equation representing a tool. A surrogate model must be informed by data from the tool with enough points to represent the solution space for the chosen number of variables with an acceptable level of error. Therefore, Design Of Experiments (DOE) is used to select points within the design space to maximize information gained on the design space while minimizing number of data points required. To represent a design space with a non-trivial number of variable parameters the number of points required still represent an amount of work which would take an inordinate amount of time via the current paradigm of manual analysis, and so an automated method was developed. The best practices of expert trajectory analysts working within NASA Marshall's Advanced Concepts Office (ACO) were implemented within a tool called multiPOST. These practices include how to use the output data from a previous run of POST to inform the next, determining whether a trajectory solution is feasible from a real-world perspective, and how to handle program execution errors. The tool was then augmented with multiprocessing capability to enable analysis on multiple trajectories simultaneously, allowing throughput to scale with available computational resources. In this update to the previous work the authors discuss issues with the method and solutions.

Patrick D Dees↗

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

An integrated high-throughput robotic platform and active learning approach for accelerated discovery of optimal electrolyte formulations

Solubility of redox-active molecules is an important determining factor of the energy density in redox flow batteries. However, the advancement of electrolyte materials discovery has been constrained by the absence of extensive experimental solubility datasets, which are crucial for leveraging data-driven methodologies. In this study, we design and investigate a highly automated workflow that synergizes a high-throughput experimentation platform with a state-of-the-art active learning algorithm to significantly enhance the solubility of redox-active molecules in organic solvents. Our platform identifies multiple solvents that achieve a remarkable solubility threshold exceeding 6.20 M for the archetype redox-active molecule, 2,1,3-benzothiadiazole, from a comprehensive library of more than 2000 potential solvents. Significantly, our integrated strategy necessitates solubility assessments for fewer than 10% of these candidates, underscoring the efficiency of our approach. Our results also show that binary solvent mixtures, particularly those incorporating 1,4-dioxane, are instrumental in boosting the solubility of 2,1,3-benzothiadiazole. Beyond designing an efficient workflow for developing high-performance redox flow batteries, our machine learning-guided high-throughput robotic platform presents a robust and general approach for expedited discovery of functional materials.

25 ENERGY STORAGE↗