Search NASA⌕ Search

SEARCH · Search NASA

Results for “outer approximation”

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 19 records

A Convexification-Based Outer-Approximation Method for Convex and Nonconvex MINLP

The advancement of domain reduction techniques has significantly enhanced the performance of solvers in mathematical programming. This paper delves into the impact of integrating convexification and domain reduction techniques within the Outer-Approximation method. We propose a refined convexification-based Outer-Approximation method alongside a Branch-and-Bound method for both convex and nonconvex Mixed-Integer Nonlinear Programming problems. These methods have been developed and incorporated into the open-source Mixed-Integer Nonlinear Decomposition Toolbox for Pyomo-MindtPy. Comprehensive benchmark tests were conducted, validating the effectiveness and reliability of our proposed algorithms. These tests highlight the improvements achieved by incorporating convexification and domain reduction techniques into the Outer-Approximation and Branch-and-Bound methods.

Optimization↗

Alternative regularizations for Outer-Approximation algorithms for convex MINLP

In this work, we extend the regularization framework from Kronqvist et al. (Math Program 180(1):285–310, 2020) by incorporating several new regularization functions and develop a regularized single-tree search method for solving convex mixed-integer nonlinear programming (MINLP) problems. We propose a set of regularization functions based on distance metrics and Lagrangean approximations, used in the projection problem for finding new integer combinations to be used within the Outer-Approximation (OA) method. The new approach, called Regularized Outer-Approximation (ROA), has been implemented as part of the open-source Mixed-integer nonlinear decomposition toolbox for Pyomo—MindtPy. We compare the OA method with seven regularization function alternatives for ROA. Moreover, we extend the LP/NLP Branch and Bound method proposed by Quesada and Grossmann (Comput Chem Eng 16(10–11):937–947, 1992) to include regularization in an algorithm denoted RLP/NLP. We provide convergence guarantees for both ROA and RLP/NLP. Finally, we perform an extensive computational experiment considering all convex MINLP problems in the benchmark library MINLPLib. The computational results show clear advantages of using regularization combined with the OA method.

Convex Mixed-integer nonlinear programming↗

Optimal Network Topology for Node-Breaker Representations With AC Power Flow Constraints

It has been demonstrated that network topology optimization (NTO) may change the topology of power system networks, and consequently, provide additional flexibility to reduce network congestion and violations. Most NTO problems are formulated based on the bus-branch model in which it is challenging to represent a realistic picture of all substation configurations. In this paper, we explore advantages of substation reconfiguration modeling based on node-breaker representations for NTO problem with full nonlinear alternating current power flow. It also proposes a tailored solution algorithm to solve this nonconvex mixed-integer nonlinear programming through the outer approximation method. The proposed solution approach iterates between a mixed-integer linear programming and a nonlinear subproblem. Additional enhancements to further accelerate the iteration process are illustrated. Numerical case studies demonstrate the relative economic and operational impact of optimal network topology with node-breaker representations.

42 ENGINEERING↗

Practical algorithms for multivariate rational approximation

We present two approaches for computing rational approximations to multivariate functions, motivated by their effectiveness as surrogate models for high-energy physics (HEP) applications. Our first approach builds on the Stieltjes process to efficiently and robustly compute the coefficients of the rational approximation. Our second approach is based on an optimization formulation that allows us to include structural constraints on the rational approximation (in particular, constraints demanding the absence of singularities), resulting in a semi-infinite optimization problem that we solve using an outer approximation approach. We present results for synthetic and real-life HEP data, and we compare the approximation quality of our approaches with that of traditional polynomial approximations.

97 MATHEMATICS AND COMPUTING↗

Deep Interior Mission: Imaging the Interior of Near-Earth Asteroids Using Radio Reflection Tomography

Near-Earth asteroids are important exploration targets since they provide clues to the evolution of the solar system. They are also of interest since they present a clear danger to Earth. Our mission objective is to image the internal structure of two NEOs using radio reflection tomography (RRT) in order to explore the record of asteroid origin and impact evolution, and to test the fundamental hypothesis that some NEOs are rubble piles rather than consolidated bodies. Our mission s RRT technique is analogous to doing a CAT scan of the asteroid from orbit. Closely sampled radar echoes are processed to yield volumetric maps of mechanical and compositional boundaries, and to measure interior material dielectric properties. The RRT instrument is a radar that operates at 5 and 15 MHz with two 30-m (tip-to-tip) dipole antennas that are used in a cross-dipole configuration. The radar transmitter and receiver electronics have heritage from JPL's MARSIS contribution to Mars Express, and the antenna is similar to systems used in IMAGE and LACE missions. The 5-MHz channel is designed to penetrate greater than 1 km of basaltic rock, and 15-MHz penetrates a few hundred meters or more. In addition to RRT volumetric imaging, we use redundant color cameras to explore the surface expressions of unit boundaries, in order to relate interior radar imaging to what is observable from spacecraft imaging and from Earth. The camera also yields stereo color imaging for geology and RRT-related compositional analysis. Gravity and high fidelity geodesy are used to explore how interior structure is expressed in shape, density, mass distribution and spin. Ion thruster propulsion is utilized by Deep Interior to enable tomographic radar mapping of multiple asteroids. Within the Discovery AO scheduling parameters we identify two targets, S-type 1999 ND43 (approximately 500 m diameter) and V-type 3908 Nyx (approximately 1 km), asteroids whose compositions bracket the diversity of solar system materials that we are likely to encounter, from undifferentiated to highly evolved. The 5-15 MHz radar is capable of probing more primitive bodies (e.g. comets or C-types) that may be available given other launch schedules. 5 MHz radar easily penetrates, with the required SNR , greater than 1 km of basalt (a good analog for Nyx). Basalt has a greater loss tangent than expected for most asteroids, although iron-rich M-types are probably not appropriate targets. 15 MHz radar penetrates the outer approximately 100 m of rocky 1 km asteroids and the deep interiors of comets. Laboratory studies of the most common NE0 materials expected (S-, C- and V-type meteorite analogs) will commence in 2005.

Safaeinili, A.↗

Reducing Conservativeness of Polytopic Linear-Parameter-Varying Robust Vehicle Sideslip Angle Observer Through Minimum-Area Convex Quadrilateral Design

The polytopic linear-parameter-varying (LPV) method has become increasingly popular for designing intelligent control and estimation systems for ground vehicles, particularly for the vehicle sideslip angle observer. In vehicle lateral dynamics, the vehicle longitudinal velocity-induced nonlinearities are conventionally outer-approximated with polytopes such as rectangles or triangles to obtain a polytopic LPV system. Yet, such polytopic approximations tend to be conservative and may lead to inadequate observer performance. To address this issue, a minimum-area convex quadrilateral construction is proposed in this paper to reduce design conservatism. The suggested design is demonstrated through the synthesis of an LPV H ∞ robust vehicle sideslip angle observer. Furthermore, a dSPACE-ASM simulation study is conducted to demonstrate the effectiveness and the advantage of the proposed polytopic construction over a baseline approach.

Zhou, Xingyu↗

Application of Hybrid Optimization-Expert System for Optimal Power Management on Board Space Power Station

The space power system has two sources of energy: photo-voltaic blankets and batteries. The optimal power management problem on-board has two broad operations: off-line power scheduling to determine the load allocation schedule of the next several hours based on the forecast of load and solar power availability. The nature of this study puts less emphasis on speed requirement for computation and more importance on the optimality of the solution. The second category problem, on-line power rescheduling, is needed in the event of occurrence of a contingency to optimally reschedule the loads to minimize the 'unused' or 'wasted' energy while keeping the priority on certain type of load and minimum disturbance of the original optimal schedule determined in the first-stage off-line study. The computational performance of the on-line 'rescheduler' is an important criterion and plays a critical role in the selection of the appropriate tool. The Howard University Center for Energy Systems and Control has developed a hybrid optimization-expert systems based power management program. The pre-scheduler has been developed using a non-linear multi-objective optimization technique called the Outer Approximation method and implemented using the General Algebraic Modeling System (GAMS). The optimization model has the capability of dealing with multiple conflicting objectives viz. maximizing energy utilization, minimizing the variation of load over a day, etc. and incorporates several complex interaction between the loads in a space system. The rescheduling is performed using an expert system developed in PROLOG which utilizes a rule-base for reallocation of the loads in an emergency condition viz. shortage of power due to solar array failure, increase of base load, addition of new activity, repetition of old activity etc. Both the modules handle decision making on battery charging and discharging and allocation of loads over a time-horizon of a day divided into intervals of 10 minutes. The models have been extensively tested using a case study for the Space Station Freedom and the results for the case study will be presented. Several future enhancements of the pre-scheduler and the 'rescheduler' have been outlined which include graphic analyzer for the on-line module, incorporating probabilistic considerations, including spatial location of the loads and the connectivity using a direct current (DC) load flow model.

Momoh, James↗

An approximation to the outer planet ephemeris errors in JPL's DE 200

The outer planet positions of JPL's recent planetary ephemeris, DE 202, have begun to show significant differences from DE 200, the basic ephemeris of the major national almanacs. The differences DE 202 - DE 200 are plotted and are assumed to approximate the errors of DE 200. For Jupiter, the difference in right ascension varies over the planet's 12 yr period between -0.1 and -0.2 arcsec throughout the century; for Saturn, the right ascension also varies over the 30 yr period, but in addition, shows a drift which reaches -0.25 arcsec at present; for Uranus, the difference is small through the first half of the century, but now has reached -0.4 arcsec; for Neptune, the error was +0.6 arcsec at 1900 and is near -1.0 arcsec by the year 2000; the error for Pluto exceeds +2.0 arcsec by the end of the century and is rapidly increasing. The declination errors are generally periodic and smaller than the right ascensions.

Standish, E. M., Jr.↗

Hydrogen ionization in the shock layer for entry into the outer planets

Approximate assessments are presented for the chemical state of the shock layer with respect to the hydrogen ionization reaction for representative entries into Saturn, Uranus, and Jupiter. It is shown that, except for a steep entry into the Saturn nominal atmosphere, the computation of all entry probe gascaps for the outer planets can probably be considered on either a frozen or chemical equilibrium basis.

Howe, J. T.↗

An approximate inviscid radiating flow field analysis for outer planet entry probes

An approximate computational technique has been developed for predicting inviscid, radiating flows about blunt probes entering atmospheres consisting of hydrogen and helium. The technique is rapid and versatile and is well suited for performing parametric trade studies for outer planet entries. Details of the computational technique, the thermodynamic correlations, the 58-step absorption coefficient model and the analytic shock shape equations are discussed. Good comparisons of the radiative heating computed by the approximate method and by detailed calculations are obtained.

Zoby, E. V.↗

Impact penetration experiments in teflon targets of variable thickness

Approximately 20.4 sq m of Teflon thermal blankets on the nonspinning Long Duration Exposure Facility (LDEF) were exposed to the orbital debris and micrometeoroid environment in low-Earth orbit (LEO) for approximately 5.7 years. Each blanket consisted of an outer layer (approximately 125 micron thick) of FEP Teflon that was backed by a vapor-deposited metal mirror (Inconel; less than 1 micron thick). The inner surface consisted of organic binders and Chemglaze thermal protective paint (approximately 50 micron thick) resulting in a somewhat variable, total blanket thickness of approximately 180 to 200 microns. There was at least one of these blankets, each exposing approximately 1.2 sq m of surface area, on nine of LDEF's 12 principal pointing directions, the exceptions being Rows 3, 9, and 12. As a consequence, these blankets represent a significant opportunity for micrometeoroid and debris studies, in general, and specifically they provide an opportunity to address those issues that require information about pointing direction (i.e., spatial density of impact events as a function of instrument orientation). During deintegration of the LDEF spacecraft at KSC, all penetration holes greater than or equal to 300 micron in diameter were documented and were recently synthesized in terms of spatial density as a function of LDEF viewing direction by. The present report describes ongoing cratering and penetration experiments in pure Teflon targets, which are intended to establish the relationships between crater or penetration-hole diameters and the associated projectile dimensions at laboratory velocities (i.e., 6 km/s). The ultimate objective of these efforts is to extract reliable mass-frequencies and associated fluxes of hypervelocity particles in LEO.

Hoerz, F.↗

First measurement of the characteristic depletion radius of dark matter haloes from weak lensing

ABSTRACT We use weak lensing observations to make the first measurement of the characteristic depletion radius, one of the three radii that characterize the region where matter is being depleted by growing haloes. The lenses are taken from the halo catalogue produced by the extended halo-based group/cluster finder applied to DESI Legacy Imaging Surveys DR9, while the sources are extracted from the DECaLS DR8 imaging data with the fourier_quad pipeline. We study halo masses $12 \lt \log (M_{\rm grp} \,[{\rm M_{\odot }}\, h^{-1}]) \le 15.3$ within redshifts 0.2 ≤ z ≤ 0.3. The virial and splashback radii are also measured and used to test the original findings on the depletion region. When binning haloes by mass, we find consistency between most of our measurements and predictions from the cosmicgrowth simulation, with exceptions to the lowest mass bins. The characteristic depletion radius is found to be roughly 2.5 times the virial radius and 1.7–3 times the splashback radius, in line with an approximately universal outer density profile, and the average enclosed density within the characteristic depletion radius is found to be roughly 29 times the mean matter density of the Universe in our sample. When binning haloes by both mass and a proxy for halo concentration, we do not detect a significant variation of the depletion radius with concentration, on which the simulation prediction is also sensitive to the choice of concentration proxy. We also confirm that the measured splashback radius varies with concentration differently from simulation predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

Passive seismic experiment

The network of seismometers installed by the Apollo 17 and other Apollo missions is described. The effects of the impacts of lunar modules and S-4B stages on the lunar surfaces are discussed. The information concerning lunar composition which is obtained by analyzing the seismic signals generated by moonquakes and meteoroid impacts are analyzed. It is concluded that the seismic activity within the moon is extremely low compared to that with the earth. The moon is characterized by a rigid, dynamically inactive outer shell, approximately 1000 kilometers thick, surrounding a core that has markedly different elastic properties.

Latham, G. V.↗

Improving turbine engine compressor performance retention through airfoil coatings

In order to evaluate the potential effectiveness of coatings in limiting erosive damage to compressor airfoils, an effort was initiated to evaluate candidate coatings for substrate alloys typically used in commercial engine high compressor blades. Laboratory and rig erosion testing of plasma deposited and diffusion coatings described in this paper have shown the potential of a two to four fold improvement in erosion life. The selective application of these coatings to approximately the outer third of the airfoil avoids coating the fatigue critical region of the blade, thus providing erosion resistance potentially without compromising the fatigue strength of the blade. Both the plasma and the diffusion coatings also offer the advantage of low initial cost and a multi-source production base.

Friedrich, L. A.↗

Timescales for planetary accretion and the structure of the protoplanetary disk

Models of planetary accretion which assume the mass of condensable matter in the protoplanetary disk was equal to that present in the planets today predict accretion timescales for the outer planets approximately or less than 10 to the 8th years. Such timescales are inconsistent with observations of star forming regions, which suggest that most of the gas in disks around one solar mass is removed in a few x 10 to the 6th years. A unified scenario was outlined for solar system formation consistent with astrophysical constraints. Jupiter's core could have grown by runaway accretion of planetesimals to a mass sufficient to initiate rapid accretion of gas in times of order of 500,000 to 5,000,000 years, provided the surface density of solids in its accretion zone was at least 5 to 10 times greater than that required by minimum mass models of the protoplanetary disk. The inner planets and the asteroids can be accounted for in this picture if the surface density of the solar nebula was relatively uniform out to Jupiter's orbit. The formation of such a protoplanetary disk requires significant transport of mass and angular momentum, and is consistent with viscous accretion disk models of the solar nebula.

Lissauer, Jack J.↗

Laboratory and theoretical studies of baroclinic processes

An understanding is being developed for processes which may be important in the atmosphere, and the definition and analysis of baroclinic experiments utilizing the geophysical fluid flow cells (GFFC) apparatus in microgravity space flights. Included are studies using numerical codes, theoretical models, and terrestrial laboratory experiments. The numerical modeling is performed in three stages: calculation of steady axisymmetric flow, calculation of fastest-growing linear eigenmodes, and nonlinear effects (first, wave-mean flow interactions, then wave-wave interactions). The code can accommodate cylindrical, spherical, or channel geometry. It uses finite differences in the vertical and meridional directions, and is spectral in the azimuthal. The theoretical work was mostly in the area of effects of topography upon the baroclinic instability problem. The laboratory experiments are performed in a cylindrical annulus which has a temperture gradient imposed upon the lower surface and an approximately isothermal outer wall, with the upper and inner surfaces being nominally thermally insulating.

Miller, Timothy↗