Computing Aviation Sparing Policies: Solving a Large Nonlinear Integer Program

Deployed US Navy aircraft carriers must stock a large number of spare parts to support the various types of aircraft embarked on the ship. The sparing policy determines the spares that will be stocked on the ship to keep the embarked aircraft ready to fly. Given a fleet of ten or more aircraft carriers and a cost of approximately 50 million dollars per carrier plus the cost of spares maintained in warehouses in the United States, the sparing problem constitutes a significant portion of the Navy’s resources. The objective of this work is to find a minimum-cost sparing policy that meets the readiness requirements of the embarked aircraft. This is a very large, nonlinear, integer optimization problem. The cost function is piecewise linear and convex while the constraint mapping is highly nonlinear. The distinguishing characteristics of this problem from an optimization viewpoint are that a large number of decision variables are required to be integer and that the nonlinear constraint functions are essentially “black box” functions; that is, they are very difficult (and expensive) to evaluate and their derivatives are not available. Moreover, they are not convex. Integer programming problems with a large number of variables are difficult to solve in general and most successful approaches to solving nonlinear integer problems have involved linear approximation and relaxation techniques that, because of the complexity of the constraint functions, are inappropriate for attacking this problem. We instead employ a pattern search method to each iteration of an interior point-type algorithm to solve the relaxed version of the problem. From the solution found by the pattern search on each interior point iteration, we begin another pattern search on the integer lattice to find a good integer solution. The best integer solution found across all interations is returned as the optimal solution. The pattern searches are distributed across a local area network of non-dedicated, heterogeneous computers in an office environment, thus, drastically reducing the time required to find the solution.     

Determination of Surface Functional Groups in Lignocellulosic Materials by Chemical Derivatization and ESCA Analysis

Simple and convenient methods for determining surface chemical composition of lignocellulosic materials are described. The methods are based on vapor phase fluorine surface derivatization with either trifluoro acetic anhydride (TFAA), tri-fluoro ethanol (TFE) or pentafluorophenyl hydrazine (PFPH) and subsequent Electron Spectroscopy for Chemical Analysis (ESCA). Model cellulosic surfaces with well defined functionalities were used to optimize the derivatization reaction conditions. Detection and accessibility of surface hydroxyl functional groups were investigated in cotton and regenerated cellulose as models. Carboxymethyl cellulose (CMC) was used as a model surface for detection and quantification of carboxylic acid groups. Theoretical conversion curves for derivatization reactions were calculated and used to evaluate the extent of the reactions on the model surfaces. It was found that the conversion was higher for the regenerated cellulose and CMC than for cotton. The protocols developed using the model surfaces were applied to a case study on wood fibers with different degrees of complexity, namely dissolving and chemithermomechanical (CTMP) pulp. Untreated and oxygen-plasma modified pulps were compared with respect to the surface composition of functional groups. According to the derivatization reactions, functionalities containing oxygen were significantly increased on the plasma-treated samples. The effect of the treatment was found to be dependent on the type of pulp. Fluorine derivatization is shown to be an unambiguous method for clear assessment of the chemical functionalities of cellulosic surfaces.     

On the Evaluation of Uncertain Courses of Action

We consider the problem of decision making under uncertainty. The fuzzy measure is introduced as a general way of representing available information about the uncertainty. It is noted that generally in uncertain environments the problem of comparing alternative courses of action is difficult because of the multiplicity of possible outcomes for any action. One approach is to convert this multiplicity of possible of outcomes associated with an alternative into a single value using a valuation function. We describe various ways of providing a valuation function when the uncertainty is represented using a fuzzy measure. We then specialize these valuation functions to the cases of probabilistic and possibilistic uncertainty.     

Mixing mechanisms of vortex ring formed by gravity slumping motion

 This paper is concerned with an experimental investigation of the mixing inside the vortex ring formed by the gravity slumping motion of a dense cloud in a less dense atmosphere. The dynamics of the spreading and instantaneous structures of the turbulent flow were examined by visualization, single and multi-point measurements of velocity and concentration for two heavy gases, carbondioxide (CO₂) and dichlorodifluoromethane (CCl₂F₂), in a configuration in that heavy gas, initially trapped in a reservoir, was released with the rise of a shutter into calm air of a sector-shaped dispersion channel. Visualization of the cloud as a whole showed a spreading motion in which an advancing frontal structure was followed by a stratified flow with a layer of dense fluid of higher velocities near the wall and, on top of it, a layer of dilute fluid whose concentration is controlled by the mixing mechanisms within the head. During the course of spreading, there was always a phase in which the head attained to a constant speed of advance, which occurred as 0.13 m/s for CO₂ and 0.48 m/s for CCl₂F₂. It was interesting to observe for CO₂ that the phase of constant speed took place in between two acceleration phases; the former was due to the initial slumping of the cloud at the exit of the reservoir, and the latter was attributed to the collapse of the head on the transition to the passive dispersion phase. Instantaneous two-dimensional velocity field, measured with particle image velocimeter (PIV), showed that the cloud overran the ambient air which caused the approaching dense fluid deflected away from the wall with significant vertical velocities and downstream-moving separation, and the air trapped under the head resulted in the density inversion which introduced further intricacy to the turbulent structure of the head. Instabilities at the upper free shear layer due to density and velocity discontinuity rolled into periodic array of vortices which engulfed a considerable amount of air as they were convected backwards over the head, but the incorporation of heavy and light fluids was completed with the appearance of microscales after the collapse on the stratified layer. Analyses of the cloud head at different downstream locations also revealed that its size remained unchanged when the speed of advance was constant, allowing the rate of change of the cloud volume being modeled with the rate of spreading. Contours of concentration obtained from digitized PIV pictures confirmed the kinematic features of the mixing revealed by the velocity field and that the concentration values within the large structures were higher than those at the upper part of the stratified layer. Motivated by the experimental observations, a semi-empirical analysis was presented to describe the results and based on local values of the Richardson and Reynolds numbers. Received: 4 October 1995 / Accepted: 4 July 1996     

Linear dependence of stationary distributions in ergodic Markov decision processes

In ergodic MDPs we consider stationary distributions of policies that coincide in all but n states, in which one of two possible actions is chosen. We give conditions and formulas for linear dependence of the stationary distributions of n+2 such policies, and show some results about combinations and mixtures of policies.     

Oscillations and concentrations in weak solutions of the incompressible fluid equations

The authors introduce a new concept of measure-valued solution for the 3-D incompressible Euler equations in order to incorporate the complex phenomena present in limits of approximate solutions of these equations. One application of the concepts developed here is the following important result: a sequence of Leray-Hopf weak solutions of the Navier-Stokes equations converges in the high Reynolds number limit to a measure-valued solution of 3-D Euler defined for all positive times. The authors present several explicit examples of solution sequences for 3-D incompressible Euler with uniformly bounded local kinetic energy which exhibit complex phenomena involving both persistence of oscillations and development of concentrations. An extensions of the concept of Young measure is developed to incorporate these complex phenomena in the measure-valued solutions constructed here.Partially supported by N.S.F. GrantPartially supported by N.S.F. Grant 84-0223 and 86-11110     

Solubility of carbon dioxide in aqueous solutions of sodium chloride: Experimental results and correlation

The solubility of carbon dioxide in aqueous solutions of sodium chloride was measured in the temperature range from 40 to 160°C, up to 6 mol-kg salt solutions and total pressures up to 10 MPa. Pitzer's⁽¹⁾ equations as well as the Chen and Evans⁽²⁾ model were used to correlate the new data. Results are reported and compared to literature data and correlations.     

Neighborhood unions and the cycle cover number of a graph

For several years, the study of neighborhood unions of graphs has given rise to important structural consequences of graphs. In particular, neighborhood conditions that give rise to hamiltonian cycles have been considered in depth. In this paper we generalize these approaches to give a bound on the smallest number of cycles in G containing all the vertices of G. We show that if for all x, y ? V(G), |N(x) ∩ N(y)| ≧ 2n/5 + 1, then V(G) is coverable by at most two cycles. Several related results and extensions to t cycles are also given.