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 469 records · Page 26

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗

On relaxations of the max k -cut problem formulations

Here, a tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.

97 MATHEMATICS AND COMPUTING↗

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

Optimal Control via Self-Generated Stochasticity

The problem of global maxima of functionals has been examined. Mathematical roots of local maxima are the same as those for a much simpler problem of finding global maximum of a multi-dimensional function. The second problem is instability even if an optimal trajectory is found, there is no guarantee that it is stable. As a result, a fundamentally new approach is introduced to optimal control based upon two new ideas. The first idea is to represent the functional to be maximized as a limit of a probability density governed by the appropriately selected Liouville equation. Then, the corresponding ordinary differential equations (ODEs) become stochastic, and that sample of the solution that has the largest value will have the highest probability to appear in ODE simulation. The main advantages of the stochastic approach are that it is not sensitive to local maxima, the function to be maximized must be only integrable but not necessarily differentiable, and global equality and inequality constraints do not cause any significant obstacles. The second idea is to remove possible instability of the optimal solution by equipping the control system with a self-stabilizing device. The applications of the proposed methodology will optimize the performance of NASA spacecraft, as well as robot performance.

Zak, Michail↗

Ozone Profiles and Tropospheric Ozone from Global Ozone Monitoring Experiment

Ozone profiles are derived from backscattered radiances in the ultraviolet spectra (290-340 nm) measured by the nadir-viewing Global Ozone Monitoring Experiment using optimal estimation. Tropospheric O3 is directly retrieved with the tropopause as one of the retrieval levels. To optimize the retrieval and improve the fitting precision needed for tropospheric O3, we perform extensive wavelength and radiometric calibrations and improve forward model inputs. Retrieved O3 profiles and tropospheric O3 agree well with coincident ozonesonde measurements, and the integrated total O3 agrees very well with Earth Probe TOMS and Dobson/Brewer total O3. The global distribution of tropospheric O3 clearly shows the influences of biomass burning, convection, and air pollution, and is generally consistent with our current understanding.

Liu, X.↗

A constraint on historic growth in global photosynthesis due to rising CO 2

Theory predicts that rising CO 2 increases global photosynthesis, a process known as CO 2 fertilization, and that this is responsible for much of the current terrestrial carbon sink. The estimated magnitude of the historic CO 2 fertilization, however, differs by an order of magnitude between long-term proxies, remote sensing-based estimates and terrestrial biosphere models. Here we constrain the likely historic effect of CO 2 on global photosynthesis by combining terrestrial biosphere models, ecological optimality theory, remote sensing approaches and an emergent constraint based on global carbon budget estimates. Our analysis suggests that CO 2 fertilization increased global annual terrestrial photosynthesis by 13.5 ± 3.5% or 15.9 ± 2.9 PgC (mean ± s.d.) between 1981 and 2020. Our results help resolve conflicting estimates of the historic sensitivity of global terrestrial photosynthesis to CO 2 and highlight the large impact anthropogenic emissions have had on ecosystems worldwide.

54 ENVIRONMENTAL SCIENCES↗

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗

Assimilation of Precipitation Information Retrieved from Passive Microwave Sensors

Estimates of surface precipitation and the associated vertical latent heating structure can be retrieved from space-borne passive microwave sensors such as the TRMM Microwave Imager (TMI) and Special Sensor Microwave/Imager (SSM/I). Studies have shown that assimilation of microwave rain rates can improve the quality of global analysis and forecast. Numerical weather prediction (NWP) centers are beginning to make operational use of these data. In the next few years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission to be launched around 2007. Realizing the full potential of these observations in data assimilation will require continued advances in retrieval algorithms, assimilation techniques, and model physics. I present results from two research efforts at NASA GSFC. The first is on the assimilation of tropical precipitation using temperature/moisture tendency corrections within a continuous variational assimilation framework to improve analyses and forecasts produced by the Goddard Earth Observing System (GEOS) global data assimilation system. Results show that assimilation of TMI and SSM/I rain rates improves not only precipitation and moisture, but also the related clouds, radiation energy fluxes, and large-scale circulations in GEOS analyses. The improved analyses also yield better short-range forecasts, quantitative precipitation forecast (QPF) threat scores, and storm track predictions. The second part consists of results from exploratory experiments on variational assimilation of convective and stratiform latent heating profiles within the general framework of model parameter estimation as a way to identify systematic model errors and optimize physical paremeterization schemes in global forecast and climate models.

Hou, Arthur Y.↗

Global sensitivity analysis in control-augmented structural synthesis

In this paper, an integrated approach to structural/control design is proposed in which variables in both the passive (structural) and active (control) disciplines of an optimization process are changed simultaneously. The global sensitivity equation (GSE) method of Sobieszczanski-Sobieski (1988) is used to obtain the behavior sensitivity derivatives necessary for the linear approximations used in the parallel multidisciplinary synthesis problem. The GSE allows for the decoupling of large systems into smaller subsystems and thus makes it possible to determine the local sensitivities of each subsystem's outputs to its inputs and parameters. The advantages in using the GSE method are demonstrated using a finite-element representation of a truss structure equipped with active lateral displacement controllers, which is undergoing forced vibration.

Bloebaum, Christina L.↗

Aeroelastic Wingbox Stiffener Topology Optimization

This work considers an aeroelastic wingbox model seeded with run-out blade stiffeners along the skins. Topology optimization is conducted within the shell webs of the stiffeners, in order to add cutouts and holes for mass reduction. This optimization is done with a global-local approach in order to moderate the computational cost: aeroelastic loads are computed at the wing-level, but the topology and sizing optimization is conducted at the panel-level. Each panel is optimized separately under stress, buckling, and adjacency constraints, and periodically reassembled to update the trimmed aeroelastic loads. The resulting topology is baselined against a design with standard full-depth solid stiffener blades, and found to weigh 7.43% less.

Stanford, Bret K.↗

Accelerating high-order mesh optimization using finite element partial assembly on GPUs

In this paper we present a new GPU-oriented mesh optimization method based on high order finite elements. Our approach relies on node movement with fixed topology, through the Target-Matrix Optimization Paradigm (TMOP) and uses a global nonlinear solve over the whole computational mesh, i.e., all mesh nodes are moved together. A key property of the method is that the mesh optimization process is recast in terms of finite element operations, which allows us to utilize recent advances in the field of GPU-accelerated high order finite element algorithms. For example, we reduce data motion by using tensor factorization and matrix-free methods, which have superior performance characteristics compared to traditional full finite element matrix assembly and offer advantages for GPU based HPC hardware. Furthermore, we describe the major mathematical components of the method along with their efficient GPU-oriented implementation. In addition, we propose an easily reproducible mesh optimization test that can serve as a performance benchmark for the mesh optimization community.

97 MATHEMATICS AND COMPUTING↗

Guidance strategies for near-optimum take-off performance in a windshear

This paper is concerned with guidance strategies for near-optimum performance in a windshear. This is a wind characterized by sharp change in intensity and direction over a relatively small region of space. The take-off problem is considered with reference to flight in a vertical plane. First, trajectories for optimum performance in a windshear are determined for different windshear models and different windshear intensities. Use is made of the methods of optimal control theory in conjunction with the dual sequential gradient-restoration algorithm (DSGRA) for optimal control problems. In this approach, global information on the wind flow field is needed. Then, guidance strategies for near-optimum performance in a windshear are developed, starting from the optimal trajectories. Specifically, three guidance schemes are presented: (1) gamma guidance, based on the relative path inclination; (2) theta guidance, based on the pitch attitude angle; and (3) acceleration guidance, based on the relative acceleration. In this approach, local information on the wind flow field is needed.

Miele, A.↗

Optimal Design of Grid-Stiffened Panels and Shells With Variable Curvature

A design strategy for optimal design of composite grid-stiffened structures with variable curvature subjected to global and local buckling constraints is developed using a discrete optimizer. An improved smeared stiffener theory is used for the global buckling analysis. Local buckling of skin segments is assessed using a Rayleigh-Ritz method that accounts for material anisotropy and transverse shear flexibility. The local buckling of stiffener segments is also assessed. Design variables are the axial and transverse stiffener spacing, stiffener height and thickness, skin laminate, and stiffening configuration. Stiffening configuration is herein defined as a design variable that indicates the combination of axial, transverse and diagonal stiffeners in the stiffened panel. The design optimization process is adapted to identify the lightest-weight stiffening configuration and stiffener spacing for grid-stiffened composite panels given the overall panel dimensions. in-plane design loads, material properties. and boundary conditions of the grid-stiffened panel or shell.

Ambur, Damodar R.↗

Beyond Component Optimization: Systems Level Biodesign for Lanthanide Recovery

Global demand for lanthanides (Ln) is projected to rise sharply over the next decade, while geographically concentrated supply chains and the low concentrations and matrix complexity of secondary feedstocks limit the reach of conventional hydro- and pyrometallurgical separation. Engineered biological systems offer a selective, low-energy alternative, and component-level advances in Ln-binding proteins, AI-designed selective scaffolds, and cell-surface display platforms now rival synthetic chelators in affinity and selectivity. These components, however, remain functionally isolated. Currently, there are no engineered chassis coupling recognition, intracellular trafficking, accumulation, and controlled release into an end-to-end pipeline. Here, we outline how new biodesign strategies and chassis selection must move beyond bioleaching to encompass the full recovery pathway. Achieving this requires integrating AI/ML-guided design, genome-scale build tools, high-throughput phenotyping, and biophysical transport modeling within a Design–Build–Test–Learn cycle tuned to recognition, trafficking, accumulation, and release.

Biodesign↗

One Stomatal Model to Rule Them All? Toward Improved Representation of Carbon and Water Exchange in Global Models

Abstract Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land‐surface model (LSM) simulations during drought. However, systematic evaluations of both optimality‐based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf‐level observations of gas‐exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr −1 . We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf‐gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses.

54 ENVIRONMENTAL SCIENCES↗

Optimization of COMPASS-U Support Structure

The global support structure for Compass Upgrade (COMPASS-U) employs many bolted and pinned components to improve the manufacturability of the whole structures. This is a departure from other designs used for highly loaded tokamaks, like Alcator C-Mod which used large forgings. Modeling of the numerous connections and interfaces was a challenge. This article will describe how to optimize the C-frame support structure with the goal to eliminate differential displacements at those interfaces and, importantly, minimize toroidal mutual displacement at the sliding joint of the toroidal field (TF) coils. The global support structure consists of 16 C-frames, and it can be simplified as a 1/8th cyclic symmetry model due to the eight vertical rods mounting that connect bottom CS and PF1-3 assemblies to the C-frames. Structure variants with different connection webs between two neighboring C-frames have been analyzed. The preloading of bolts can be split into several preload substeps to simulate the real assembly process of the COMPASS-U support structure. Additionally, the TF crown and the wedge plate can be further improved to become a solid part to better hold the TF coils.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A multilevel control system for the large space telescope

A multilevel scheme was proposed for control of Large Space Telescope (LST) modeled by a three-axis-six-order nonlinear equation. Local controllers were used on the subsystem level to stabilize motions corresponding to the three axes. Global controllers were applied to reduce (and sometimes nullify) the interactions among the subsystems. A multilevel optimization method was developed whereby local quadratic optimizations were performed on the subsystem level, and global control was again used to reduce (nullify) the effect of interactions. The multilevel stabilization and optimization methods are presented as general tools for design and then used in the design of the LST Control System. The methods are entirely computerized, so that they can accommodate higher order LST models with both conceptual and numerical advantages over standard straightforward design techniques.

D D Siljak↗