A Recursive Algorithm for Order Cycle-time that Minimizes Logistics Cost

A firm's logistics cost, including shipping and inventory-carrying, is a substantial percentage of its sales. Nevertheless, typical inventory-control methods ignore or insufficiently represent the shipping cost. This paper describes a recursive algorithm that determines the reorder cycle-time that minimizes total logistics cost. It allows for a realistic accounting of shipping cost, which is modelled here as a function of shipping distance and weight. The algorithm uses a relaxation procedure to identify a suitable initial approximation to the optimal order cycle-time and then, through a series of recursive steps, moves to the optimal result. We demonstrate the algorithm with a single item, with a group of items that share a common order cycle, and with multi-items when item demands are random variables. Experience with this algorithm indicates that it converges to the optimal result in a very few steps.     

Formation of hydrogenated amorphous carbon films from polymer pyrolysis

Free-standing and supported hydrogenated amorphous carbon films (a:C–H) were prepared upon pyrolysis of the polymer formed by ethanolamine (EA) and citric acid (CA), under an ambient atmosphere at 300 °C. EA facilitates the formation of the macroscopic films, while CA is essential for obtaining the a:C–H microstructure, which comprises a mixture of sp² and sp³ carbon. Received: 29 May 2001 / Accepted: 17 August 2001 / Published online: 20 December 2001     

Structures and stability of Al<Subscript><Emphasis Type="Bold">7</Emphasis></Subscript>C and Al<Subscript><Emphasis Type="Bold">7</Emphasis></Subscript>N clusters

We report results of the atomic and electronic structures of Al₇C cluster using ab initio molecular dynamics with ultrasoft pseudopotentials and generalized gradient approximation. The lowest energy structure is found to be the one in which carbon atom occupies an interstitial position in Al₇ cluster. The electronic structure shows that the recent observation [Chem. Phys. Lett. 316, 31 (2000)] of magic behavior of Al₇C^- cluster is due to a large highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) gap which makes Al₇C^- chemically inert. These results have further led us to the finding of a new neutral magic cluster Al₇N which has the same number of valence electrons as in Al₇C^- and a large HOMO-LUMO gap of 1.99 eV. Further, calculations have been carried out on (Al₇N)₂ to study interaction between magic clusters. Received 28 July 2001     

Nonlinear optical parameters and optical limitation in cobalt-doped polyvinylpyrrolidone solutions

Nonlinear optical parameters (nonlinear refractive indices and nonlinear absorption coefficients) of solutions of polyvinylpyrrolidone doped by cobalt to various concentrations are measured at the lasing and second-harmonic wavelengths of a picosecond Nd:YAG laser (λ=1064 and 532 nm, respectively, and τ=35 ps). Data on optical limitation in these solutions are presented. The absence of nonlinear absorption in the IR spectral range and its significant effect in the visible range are demonstrated. Optical limitation at a wavelength of 1064 nm is related to defocusing, whereas at 532 nm, this effect is caused by two-photon absorption and partially by inverse saturated absorption and defocusing. Nonlinear optical parameters of metal-polymer complexes are reported.     

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow vs. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque G = T/ν² L, where L is the cylinder length, scales with Reynolds number R and gap width η, G = 1.46η^3/2(1 - η)^-7/4 R ^3/2. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: G = 0.50 . These predictions are found to be in excellent agreement with avail-able experimental data. Predictions for scaling of velocity fluctuations are also provided. Received 7 June 2001 and Received in final form 7 December 2001     

A Stochastic Analysis of a Staffing Problem

Consider an organization with several categories of staff, each having a specified number of members. The age of entry into service and the age of retirement are both specified. Withdrawals occur at time-independent rates. Under this set-up, the authors have computed the maximum ages at which promotions should take place, assuming that the staff members continuing at various ages are random. Some illustrative calculations have been carried out.     

Three Electrons in a Two-Dimensional Parabolic Trap: The Relative Motion Solution

In agreement with the Kohn theorem the relative motion (rel) of three electrons in a two-dimensional parabolic trap separates from the centre-of-mass (CM) motion. By introducing new coordinates the Hamiltonian for relative motion in the approximation of non-interacting electrons can be taken to the normal form. The eigenstates of the normalized Hamiltonian are products of the Fock-Darwin states for normal modes. The energy levels for relative motion are obtained by diagonalizing the exact Hamiltonian in the eigenbasis for the non-interacting case. In this basis the interaction matrix elements can be obtained in the analytical form. Since the rank of the Hamiltonian matrix is significantly reduced, the calculations are faster and more accurate than those for the full (CM + rel) motion. This advantage is especially important for the calculations of excited states and the analysis of energy spectra.     

Primal cutting plane algorithms revisited

Dual fractional cutting plane algorithms, in which cutting planes are used to iteratively tighten a linear relaxation of an integer program, are well-known and form the basis of the highly successful branch-and-cut method. It is rather less well-known that various primal cutting plane algorithms were developed in the 1960s, for example by Young. In a primal algorithm, the main role of the cutting planes is to enable a feasible solution to the original problem to be improved. Research on these algorithms has been almost non-existent.  In this paper we argue for a re-examination of these primal methods. We describe a new primal algorithm for pure 0-1 problems based on strong valid inequalities and give some encouraging computational results. Possible extensions to the case of general mixed-integer programs are also discussed.     

Rank-perfect and" weakly rank-perfect graphs

An edge e of a perfect graph G is critical if G−e is imperfect. We would like to decide whether G−e is still “almost perfect” or already “very imperfect”. Via relaxations of the stable set polytope of a graph, we define two superclasses of perfect graphs: rank-perfect and weakly rank-perfect graphs. Membership in those two classes indicates how far an imperfect graph is away from being perfect. We study the cases, when a critical edge is removed from the line graph of a bipartite graph or from the complement of such a graph.     

Possible Loss and Recovery of Gibbsianness¶During the Stochastic Evolution of Gibbs Measures

We consider Ising-spin systems starting from an initial Gibbs measure ν and evolving under a spin-flip dynamics towards a reversible Gibbs measure μ≠ν. Both ν and μ are assumed to have a translation-invariant finite-range interaction. We study the Gibbsian character of the measure νS(t) at time t and show the following: (1) For all ν and μ, νS(t) is Gibbs for small t. (2) If both ν and μ have a high or infinite temperature, then νS(t) is Gibbs for all t > 0. (3) If ν has a low non-zero temperature and a zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t and non-Gibbs for large t. (4) If ν has a low non-zero temperature and a non-zero magnetic field and μ has a high or infinite temperature, then νS(t) is Gibbs for small t, non-Gibbs for intermediate t, and Gibbs for large t. The regime where μ has a low or zero temperature and t is not small remains open. This regime presumably allows for many different scenarios. Received: 26 April 2001 / Accepted: 10 October 2001