Stochastic Algorithms for Discrete Parameter Simulation Optimization

We present two efficient discrete parameter simulation optimization (DPSO) algorithms for the long-run average cost objective. One of these algorithms uses the smoothed functional approximation (SFA) procedure, while the other is based on simultaneous perturbation stochastic approximation (SPSA). The use of SFA for DPSO had not been proposed previously in the literature. Further, both algorithms adopt an interesting technique of random projections that we present here for the first time. We give a proof of convergence of our algorithms. Next, we present detailed numerical experiments on a problem of admission control with dependent service times. We consider two different settings involving parameter sets that have moderate and large sizes, respectively. On the first setting, we also show performance comparisons with the well-studied optimal computing budget allocation (OCBA) algorithm and also the equal allocation algorithm. Note to Practitioners-Even though SPSA and SFA have been devised in the literature for continuous optimization problems, our results indicate that they can be powerful techniques even when they are adapted to discrete optimization settings. OCBA is widely recognized as one of the most powerful methods for discrete optimization when the parameter sets are of small or moderate size. On a setting involving a parameter set of size 100, we observe that when the computing budget is small, both SPSA and OCBA show similar performance and are better in comparison to SFA, however, as the computing budget is increased, SPSA and SFA show better performance than OCBA. Both our algorithms also show good performance when the parameter set has a size of 10(8). SFA is seen to show the best overall performance. Unlike most other DPSO algorithms in the literature, an advantage with our algorithms is that they are easily implementable regardless of the size of the parameter sets and show good performance in both scenarios.

KCL, Ray Theories and Application to Sonic Boom

We first review a general formulation of ray theory and write down the conservation forms of the equations of a weakly nonlinear ray theory (WNLRT) and a shock ray theory (SRT) for a weak shock in a polytropic gas. Then we present a formulation of the problem of sonic boom by a maneuvering aerofoil as a one parameter family of Cauchy problems. The system of equations in conservation form is hyperbolic for a range of values of the parameter and has elliptic nature else where, showing that unlike the leading shock, the trailing shock is always smooth.

Calculation of the front part of the sonic boom signature for a Maneuvering aerofopil

The sonic boom at a large distance from its source consists of a leading shock, a trailing shock and a one parameter family of nonlinear wavefronts in between these shocks. A new ray theoretical method using a shock ray theory and a weakly nonlinear lay theory has been used to obtain the shock fronts and wavefronts respectively, for a maneuvering aerofoil in a homogeneous medium. This method introduces a one parameter family of Cauchy problems to calculate the shock and wave fronts emerging from the surface of the aerofoil. These problems are solved numerically to obtain the leading shock front and the nonlinear wavefronts emerging from the front portion of the aerofoil.

Laser damage threshold studies on urea L-malic acid: A nonlinear optical crystal

A detailed study of surface laser damage performed on a nonlinear optical crystal, urea L-malic acid, using 7 ns laser pulses at 10 Hz repetition rate from a Q-switched Nd:YAG laser at wavelengths of 532 and 1064 nm is reported. The single shot and multiple shot surface laser damage threshold values are determined to be 26.64±0.19 and 20.60±0.36 GW cm−2 at 1064 nm and 18.44±0.31 and 7.52±0.22 GW cm−2 at 532 nm laser radiation, respectively. The laser damage anisotropy is consistent with the Vickers mechanical hardness measurement performed along three crystallographic directions. The Knoop polar plot also reflects the damage morphology. Our investigation reveals a direct correlation between the laser damage profile and hardness anisotropy. Thermal breakdown of the crystal is identified as the possible mechanism of laser induced surface damage.

Mode-coupling glass transition in a fluid confined by a periodic potential

We show that a fluid under strong spatially periodic confinement displays a glass transition within mode-coupling theory at a much lower density than the corresponding bulk system. We use fluctuating hydrodynamics, with confinement imposed through a periodic potential whose wavelength plays an important role in our treatment. To make the calculation tractable we implement a detailed calculation in one dimension. Although we do not expect simple 1d fluids to show a glass transition, our results are indicative of the behavior expected in higher dimensions. In a certain region of parameter space we observe a three-step relaxation reported recently in computer simulations [S. H. Krishnan, Ph.D. thesis, Indian Institute of Science (2005); Kim et al., Eur. Phys. J. Special Topics 189, 135 (2010)] and a glass-glass transition. We compare our results to those of Krakoviack [Phys. Rev. E 75, 031503 (2007)] and Lang et al. [Phys. Rev. Lett. 105, 125701 (2010)].

Nonlinear dielectric behavior in three-component ferroelectric superlattices

Three-component ferroelectric superlattices consisting of alternating layers of SrTiO3, BaTiO3, and CaTiO3 (SBC) with variable interlayer thickness were fabricated on Pt(111)/TiO2/SiO2/Si (100) substrates by pulsed laser deposition. The presence of satellite reflections in x-ray-diffraction analysis and a periodic concentration of Sr, Ba, and Ca throughout the film in depth profile of secondary ion mass spectrometry analysis confirm the fabrication of superlattice structures. The Pr (remnant polarization) and Ps (saturation polarization) of SBC superlattice with 16.4-nm individual layer thickness (SBC16.4) were found to be around 4.96 and 34 μC/cm2, respectively. The dependence of polarization on individual layer thickness and lattice strain were studied in order to investigate the size dependence of the dielectric properties. The dielectric constant of these superlattices was found to be much higher than the individual component layers present in the superlattice configuration. The relatively higher tunability ( ∼ 55%) obtained around 300 K indicates that the superlattice is a potential electrically tunable material for microwave applications at room temperature. The enhanced dielectric properties were thus discussed in terms of the interfacial strain driven polar region due to high lattice mismatch and electrostatic coupling due to polarization mismatch between individual layers.

The R-x, R2-x Planes; Some New Ideas for the Study of Nonlinear Systems

Use of some new planes such as the R-x, R2-x (where R represents in the n-dimensional phase space, the radius vector from the origin to any point on the trajectory described by the system) is suggested for analysis of nonlinear systems of any kind. The stability conditions in these planes are given. For easy understanding of the method, the transformation from the phase plane to the R-x, R2-x planes is brought out for second-order systems. In general, while these planes serve as useful as the phase plane, they have proved to be simpler in determining quickly the general behavior of certain classes of second-order nonlinear systems. A chart and a simple formula are suggested to evaluate time easily from the R-x and R2-x trajectories, respectively. A means of solving higher-order nonlinear systems is also illustrated. Finally, a comparative study of the trajectories near singular points on the phase plane and on the new planes is made.

The x¿n-x Plane for Analysis of Certain Second-Order Nonlinear Systems

Analysis of certain second-order nonlinear systems, not easily amenable to the phase-plane methods, and described by either of the following differential equations xÿn-2ÿ+ f(x)xÿ2n+g(x)xÿn+h(x)=0 ÿ+f(x)xÿn+h(x)=0 n≫0 can be effected easily by drawing the entire portrait of trajectories on a new plane; that is, on one of the xÿnÿx planes. Simple equations are given to evaluate time from a trajectory on any of these n planes. Poincaré's fundamental phase plane xÿÿx is conceived of as the simplest case of the general xÿnÿx plane.

Techniques for Analysis of Certain Nonlinear Systems

This paper suggests the use of simple transformations like ÿ=kx, kx2 for second-order nonlinear differential equations to effect rapid plotting of the phase-plane trajectories. The method is particularly helpful in determining quickly the trajectory slopes along simple curves in any desired region of the phase plane. New planes such as the tÿ-x, tÿ2-x are considered for the study of some groups of nonlinear time-varying systems. Suggestions for solving certain higher-order nonlinear systems are also made.