Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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

Sufficient conditions for a local minimum of the Bolza problem with a scalar terminal point constraint

Sufficient conditions for a weak relative minimum for a form of the Bolza problem of variational calculus are derived. Testing of the second-order conditions requires the backward integration of fewer matrix elements than in the case of most previously published sets of conditions. The derivation is felt to be more complete and straightforward than previous derivations. The variational problem considered is relatively simple, with just a scalar constraint to implicitly or explicitly determine the final time, in order to avoid the complexities associated with controllability considerations. Sufficient conditions for a local minimum for more general optimal control problems can be approached by building upon the derivation and results presented here.

Wood, Lincoln J.↗

Effect of a localized minimum in equatorial field strength on resistive tearing instability in the geomagnetotail

A two-dimensional, resistive-MHD computer code is used to investigate the spontaneous reconnection of magnetotaillike configurations. The initial conditions adopted in the simulations are of two types: (1) in which the equatorial normal magnetic field component B(ze) declines monotonically down the tail, and (2) in which B(ze) exhibits a deep minimum in the near-earth plasma sheet. To represent the case where the earthward convection stops before the X line forms, zero-flow boundary conditions are imposed at the edges of the computational box. The initial configurations are in equilibrium and table within ideal MHD. The dynamic evolution of the system starts after the resistivity is turned on. The main results of these simulations basically support the neutral-line model of substorms and confirm Birn's (1980) computer studies. Spontaneous formation of an X-type neutral point and a single O-type plasmoid with strong tailward flow on the tailward side of the X point is found. The time interval from the turning on of the resistivity to the formation of a plasmoid is much shorter in the case where there is an initial deep minimum. A simple analytic calculation is also carried out to demonstrate why the configuration with a deep minimum is more susceptible to the development of the neutral point.

Hau, L.-N.↗

Sufficient conditions for a local minimum of the Bolza problem with multiple terminal point constraints

Sufficient conditions for a weak relative minimum are derived for a form of the Bolza problem of variational calculus. The derivation ties together first-order and second-order conditions in a unified consistent manner, addresses controllability considerations in some detail, and clarifies some small inconsistencies in earlier work. The resulting second-order conditions of optimality involve the integration of fewer backward-sweep matrix elements than the standard conditions in the literature for problems with an unspecified final time. As a result, the backward-sweep matrices have the same general structure and dimensionality whether the final time is specified or free, with the terminal values of two of three sweep matrices being more complicated in the latter case.

Wood, Lincoln J.↗

Sufficient conditions for a local minimum of the Bolza problem with a variable initial point

Sufficient conditions for a weak relative minimum of a form of the Bolza problem of variational calculus are derived. The variational problem considered includes arbitrary numbers of constraints on the initial and terminal points, but assumes no control constraints or interior point constraints. Testing of the second-order conditions requires the backward integration of fewer matrix elements than in the case of previously published sets of conditions. The derivation, though lengthy, is felt to be simpler in concept than previous derivations. The application of these conditions is demonstrated in the context of a simple geometric problem.

Wood, Lincoln J.↗

Thermal Radiator Pointing for International Space Station

In order to provide thermal radiation environments that result in adequate beat rejection, the single-phase, liquid ammonia (NH3) heat rejection system on the International Space Station (ISS) requires that its two thermal radiator wings be dynamically rotated as the ISS travels through its orbit. This paper discusses the closed-loop, thermal radiator pointing system that is used on ISS to ensure adequate heat rejection by the radiators, while preventing freezing of the ammonia under low heat loads and cold-environmental conditions. Although initial designs used an open-loop approach for radiator pointing, concerns about performance robustness, algorithm complexity, memory requirements, and sustaining support drove the development of a more robust, simpler, closed-loop system. Hence, the challenge of the closed-loop system was to utilize existing sensors, actuators and computers to fit into the existing hardware and software architecture of the ISS. Using a proportional-integral (PI) control architecture with limited output and an anti-windup integrator, the temperature of the ammonia coming out of the radiator is measured and controlled by adjusting the radiator wing orientation. The radiator wing orientation for the local minimum environment is fed forward to the control system, and the closed-loop controller is used to generate a bias off of that local minimum environment in order to heat up the ammonia when necessary to avoid freezing. In the earth's shadow, the controller is suspended and the radiator wing is oriented to face the earth, the local maximum thermal environment which further prevents freezing of the ammonia. This control architecture is shown to provide adequate heat rejection and avoid freezing of the ammonia, even though the physical system consists of large transport delays and time-varying dynamics which change dramatically due to orbit motion and variable heat loads.

Green, Scott↗

Heliospheric shocks and catastrophe theory

Various configurations of forward and reverse shocks that occur in the outer heliosphere can be classified using catastrophe theory. The existence of a forward shock is associated with a local maximum of a polynomial, and the existence of a reverse shock is associated with a local minimum of a polynomial. A configuration with N forward shocks and N reverse shocks corresponds to a polynomial with N maxima and N minima. The formation of forward and reverse shocks corresponds to the creation of maxima and minima of a polynomial, which is described by the separatrices of the catastrophes. The coalescence of two forward (reverse) shocks corresponds to the situation when two maxima (minima) of a polynomial have equal values, and the interaction of a forward shock with a reverse shock corresponds to a polynomial with a local maximum equal to a local minimum; these situations are described by the Maxwell sets of the appropriate catastrophes.

Burlaga, L. F.↗

Algorithms For Segmentation Of Complex-Amplitude SAR Data

Several algorithms implement improved method of segmenting highly speckled, high-resolution, complex-amplitude synthetic-aperture-radar (SAR) digitized images into regions, within each backscattering characteristics similar or homogeneous from place to place. Method provides for approximate, deterministic solution by two alternative algorithms almost always converging to local minimums: one, Iterative Conditional Modes (ICM) algorithm, which locally maximizes posterior probability density of region labels; other, Maximum Posterior Marginal (MPM) algorithm, which maximizes posterior marginal density of region labels at each pixel location. ICM algorithm optimizes reconstruction of underlying scene. MPM algorithm minimizes expected number of misclassified pixels, possibly better in remote sensing of natural scenes.

Rignot, Eric J. M.↗

Markov Chain Monte Carlo Bayesian Learning for Neural Networks

Conventional training methods for neural networks involve starting al a random location in the solution space of the network weights, navigating an error hyper surface to reach a minimum, and sometime stochastic based techniques (e.g., genetic algorithms) to avoid entrapment in a local minimum. It is further typically necessary to preprocess the data (e.g., normalization) to keep the training algorithm on course. Conversely, Bayesian based learning is an epistemological approach concerned with formally updating the plausibility of competing candidate hypotheses thereby obtaining a posterior distribution for the network weights conditioned on the available data and a prior distribution. In this paper, we developed a powerful methodology for estimating the full residual uncertainty in network weights and therefore network predictions by using a modified Jeffery's prior combined with a Metropolis Markov Chain Monte Carlo method.

Goodrich, Michael S.↗

A Change in the Solar He II EUV Global Network Structure as an Indicator of the Geo-Effectiveness of Solar Minima

Solar activity during 2007 - 2009 was very low, causing anomalously low thermospheric density. A comparison of solar extreme ultraviolet (EUV) irradiance in the He II spectral band (26 to 34 nm) from the Solar Extreme ultraviolet Monitor (SEM), one of instruments on the Charge Element and Isotope Analysis System (CELIAS) on board the Solar and Heliospheric Observatory (SOHO) for the two latest solar minima showed a decrease of the absolute irradiance of about 15 +/- 6 % during the solar minimum between Cycles 23 and 24 compared with the Cycle 22/23 minimum when a yearly running-mean filter was used. We found that some local, shorter-term minima including those with the same absolute EUV flux in the SEM spectral band show a higher concentration of spatial power in the global network structure from the 30.4 nm SOHO/Extreme ultraviolet Imaging Telescope (EIT) images for the local minimum of 1996 compared with the minima of 2008 - 2011.We interpret this higher concentration of spatial power in the transition region's global network structure as a larger number of larger-area features on the solar disk. These changes in the global network structure during solar minima may characterize, in part, the geo-effectiveness of the solar He II EUV irradiance in addition to the estimations based on its absolute levels.

Solar extreme ultraviolet irradiance↗

The density minimum at the earth's magnetic equator

Observations of the density structure in the plasmapause region reveal the existence of a local minimum in the total electron density at the magnetic equator. Data from the plasma wave instrument and ion mass spectrometer on the DE-1 satellite are used to study this phenomenon. The density depletion typically extends from +/- 5 to +/- 20 deg in latitude and is found at altitudes from 2 to 5 RE. Density depletions of 10-70 percent are found in regions where the off-equator density ranges from 10 to 1000/cu cm. This density structure is associated with equator crossings where the thermal plasma has been heated over normal plasmasphere values. The heated plasma is the equatorially trapped plasma previously reported from DE 1 and the SCATHA satellite. Within the plasmasphere, the drop in total (electron) density corresponds to a decrease in the cold-ion density, in both H(+) and He(+). There is a rough pressure balance provided by the warm tail of the distribution, which is a few percent by density but 1-2 orders of magnitude higher in temperature.

Olsen, R. C.↗

New Evidence for Equatorially Trapped Thermal Plasma During Early Post-Storm Recovery

Almost 20 years ago Olsen et al. [1987] reported Dynamics Explorer 1 Retarding Ion Mass Spectrometer observations of equatorially trapped, cold ions in the vicinity of the plasmapause. In that study the trapped population corresponded to a local minimum in density at the magnetic equator. During that time period observations were uncovered of a local maximum in plasma density at the equator. Until IMAGE there has been no good opportunity to experimentally revisit this topic, however until now no direct evidence of a relevant equatorial process has been recognized near the plasmapause during early recovery conditions. It appears that evidence has now been found in both the Extreme Ultraviolet Imager and Radio Plasma Imager observations. The observations, conditions, and properties of what appears to be an equatorially trapped and enhanced density near the magnetic equator will be presented and discussed.

Gallagher, D. L.↗

The faintest stars - The luminosity and mass functions at the bottom of the main sequence

We present IR K-band photometry of complete samples of VLM candidates constructed from IIIaF and IVN plates in 10 fields taken as part of the POSSII and UKSRC surveys. Using the I-K colors constructed for these stars we estimate a bolometric luminosity function which extends to M(Bol) = 13.75. We find significant evidence for a luminosity function decreasing toward these luminosities. We also find that our results are consistent with those of studies based on the Nearby Star sample, when those data are presented as a bolometric luminosity function. We convert our observed luminosity function into a mass function, which extends with reasonable statistics to 0.08 solar masses - the H-burning minimum mass. We find significant evidence for features in the mass function at these masses. Specifically, the mass function 'turns over' at 0.25 solar mass, goes through a local minimum at about 0.15 solar mass, and may increase again below 0.1 solar mass - none of these features are predicted by any of the current theories of star formation. Lastly, the mass density we observe just above the H-burning minimum mass makes it difficult to envisage brown dwarfs contributing significant quantities of missing mass without invoking either a mass function in this region significantly steeper than that seen for main-sequence stars, or an extremely low cutoff mass to the mass function.

Tinney, Christopher G.↗

When Gravity Fails: Local Search Topology

Local search algorithms for combinatorial search problems frequently encounter a sequence of states in which it is impossible to improve the value of the objective function; moves through these regions, called {\em plateau moves), dominate the time spent in local search. We analyze and characterize {\em plateaus) for three different classes of randomly generated Boolean Satisfiability problems. We identify several interesting features of plateaus that impact the performance of local search algorithms. We show that local minima tend to be small but occasionally may be very large. We also show that local minima can be escaped without unsatisfying a large number of clauses, but that systematically searching for an escape route may be computationally expensive if the local minimum is large. We show that plateaus with exits, called benches, tend to be much larger than minima, and that some benches have very few exit states which local search can use to escape. We show that the solutions (i.e. global minima) of randomly generated problem instances form clusters, which behave similarly to local minima. We revisit several enhancements of local search algorithms and explain their performance in light of our results. Finally we discuss strategies for creating the next generation of local search algorithms.

Frank, Jeremy↗

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗

A finite state, finite memory minimum principle, part 2

The development of the theory of the finite - state, finite - memory (FSFM) stochastic control problem is discussed. The sufficiency of the FSFM minimum principle (which is in general only a necessary condition) was investigated. By introducing the notion of a signaling strategy as defined in the literature on games, conditions under which the FSFM minimum principle is sufficient were determined. This result explicitly interconnects the information structure of the FSFM problem with its optimality conditions. The min-H algorithm for the FSFM problem was studied. It is demonstrated that a version of the algorithm always converges to a particular type of local minimum termed a person - by - person extremal.

Sandell, N. R., Jr.↗

Wavefront Control Toolbox for James Webb Space Telescope Testbed

We have developed a Matlab toolbox for wavefront control of optical systems. We have applied this toolbox to the optical models of James Webb Space Telescope (JWST) in general and to the JWST Testbed Telescope (TBT) in particular, implementing both unconstrained and constrained wavefront optimization to correct for possible misalignments present on the segmented primary mirror or the monolithic secondary mirror. The optical models implemented in Zemax optical design program and information is exchanged between Matlab and Zemax via the Dynamic Data Exchange (DDE) interface. The model configuration is managed using the XML protocol. The optimization algorithm uses influence functions for each adjustable degree of freedom of the optical mode. The iterative and non-iterative algorithms have been developed to converge to a local minimum of the root-mean-square (rms) of wavefront error using singular value decomposition technique of the control matrix of influence functions. The toolkit is highly modular and allows the user to choose control strategies for the degrees of freedom to be adjusted on a given iteration and wavefront convergence criterion. As the influence functions are nonlinear over the control parameter space, the toolkit also allows for trade-offs between frequency of updating the local influence functions and execution speed. The functionality of the toolbox and the validity of the underlying algorithms have been verified through extensive simulations.

Shiri, Ron↗

Path planning for assembly of strut-based structures

A path planning method with collision avoidance for a general single chain nonredundant or redundant robot is proposed. Joint range boundary overruns are also avoided. The result is a sequence of joint vectors which are passed to a trajectory planner. A potential field algorithm in joint space computes incremental joint vectors delta-q = delta-q(sub a) + delta-q(sub c) + delta-q(sub r). Adding delta-q to the robot's current joint vector leads to the next step in the path. Delta-q(sub a) is obtained by computing the minimum norm solution of the underdetermined linear system J delta-q(sub a) = x(sub a) where x(sub a) is a translational and rotational force vector that attracts the robot to its goal position and orientation. J is the manipulator Jacobian. Delta-q(sub c) is a collision avoidance term encompassing collisions between the robot (links and payload) and obstacles in the environment as well as collisions among links and payload of the robot themselves. It is obtained in joint space directly. Delta-q(sub r) is a function of the current joint vector and avoids joint range overruns. A higher level discrete search over candidate safe positions is used to provide alternatives in case the potential field algorithm encounters a local minimum and thus fails to reach the goal. The best first search algorithm A* is used for graph search. Symmetry properties of the payload and equivalent rotations are exploited to further enlarge the number of alternatives passed to the potential field algorithm.

Muenger, Rolf↗