A two timescale stochastic approximation scheme for simulation-based parametric optimization

A two timescale stochastic approximation scheme which uses coupled iterations is used for simulation-based parametric optimization as an alternative to traditional "infinitesimal perturbation analysis" schemes, It avoids the aggregation of data present in many other schemes. Its convergence is analyzed, and a queueing example is presented.

Multiscale stochastic approximation for parametric optimization of hidden Markov models

A two-time scale stochastic approximation algorithm is proposed for simulation-based parametric optimization of hidden Markov models, as an alternative to the traditional approaches to ''infinitesimal perturbation analysis.'' Its convergence is analyzed, and a queueing example is presented.

Oxide nanotubes prepared using carbon nanotubes as templates

Hollow nanotubes of SiO2, Al2O3, V2O5, and MoO3 have been prepared using carbon nanotubes as templates. The procedure involves coating the carbon nanotubes with tetraethylorthosilicate, aluminum isopropoxide, or vanadium pentoxide gel, followed by calcination and heating at higher temperatures in air to oxidize the carbon. SiO2 nanotubes containing transition metal ions have been prepared by this procedure since such materials may be of use in catalysis.

Strontium and calcium zirconyl citrates as precursors for the low-temperature synthesis of SrZrO3 and CaZrO3 fine powders

Synthesis and the thermal decomposition behavior of new molecular precursors, strontium, and calcium zirconyl citrates are presented. The pathway to the metazirconate formation has been found to proceed through a multistep process. The precursors yield SrZrO3 and CaZrO3 fine powders at temperatures as low as 650 degrees C. Physico-chemical, spectroscopic, thermoanalytical, and microscopic techniques have enabled the identification of the sequence of events leading to the perovskite formation and proposition of a thermolysis scheme. Retention of the molecular level mixing of the metal ions during the course of the precursor decomposition is supported by these techniques. Prior to the formation of MZrO3 (M = Sr and Ca) an ionic oxycarbonate, M2Zr2O5CO3 (M = SI. and Ca), intermediate is produced by the thermal decomposition of the citrate precursors.

Studies on age-hardening characteristics of ceramic particle/matrix interfaces in Al-Cu-SiCp composites using ultra low-load-dynamic microhardness measurements

Ultra low-load-dynamic microhardness testing facilitates the hardness measurements in a very low volume of the material and thus is suited for characterization of the interfaces in MMC's. This paper details the studies on age-hardening behavior of the interfaces in Al-Cu-5SiC(p) composites characterized using this technique. Results of hardness studies have been further substantiated by TEM observations. In the solution-treated condition, hardness is maximum at the particle/matrix interface and decreases with increasing distance from the interface. This could be attributed to the presence of maximum dislocation density at the interface which decreases with increasing distance from the interface. In the case of composites subjected to high temperature aging, hardening at the interface is found to be faster than the bulk matrix and the aging kinetics becomes progressively slower with increasing distance from the interface. This is attributed to the dislocation density gradient at the interface, leading to enhanced nucleation and growth of precipitates at the interface compared to the bulk matrix. TEM observations reveal that the sizes of the precipitates decrease with increasing distance from the interface and thus confirms the retardation in aging kinetics with increasing distance from the interface.

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Alfven resonances and forced reconnection

The correspondence between the forced magnetic reconnection induced by perturbing the boundary of the simple Taylor model and the surface-wave-induced magnetic reconnection given by Alfven resonance theory is pointed out explicitly by showing that the theory of forced magnetic reconnection is actually embedded in the Alfven resonance theory. The advantages of viewing the forced reconnection as surface-wave-induced reconnection are briefly discussed in the context of the formation of small-scale structures at the magnetospheric boundary and solar coronal heating.

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Boundary conditions at a rigid wall for rough granular gases

We derive boundary conditions at a rigid wall for a granular material comprising rough, inelastic particles. Our analysis is confined to the rapid flow, or granular gas, regime in which grains interact by impulsive collisions. We use the Chapman-Enskog expansion in the kinetic theory of dense gases, extended for inelastic and rough particles, to determine the relevant fluxes to the wall. As in previous studies, we assume that the particles are spheres, and that the wall is corrugated by hemispheres rigidly attached to it. Collisions between the particles and the wall hemispheres are characterized by coefficients of restitution and roughness. We derive boundary conditions for the two limiting cases of nearly smooth and nearly perfectly rough spheres, as a hydrodynamic description of granular gases comprising rough spheres is appropriate only in these limits. The results are illustrated by applying the equations of motion and boundary conditions to the problem of plane Couette flow.

Woody plant seedling distribution under invasive Lantana camara thickets in a dry-forest plot in Mudumalai, southern India

Lantana camara, a shrub of Central and South American origin, has become invasive across dry forests worldwide. The effect of the thicket-forming habit of L. camara as a dispersal and recruitment barrier in a community of native woody seedlings was examined in a 50-ha permanent plot located in the seasonally dry forest of Mudumalai, southern India. Sixty 100-m(2) plots were enumerated for native woody seedlings between 10-100 cm in height. Of these, 30 plots had no L. camara thickets, while the other 30 had dense thickets. The frequency of occurrence and abundance of seedlings were modelled as a function of dispersal mode (mammal, bird or mechanical) and affinities to forest habitats (dry forest, moist forest or ubiquitous) as well as presence or absence of dense L. camara thickets. Furthermore, frequency of occurrence and abundance of individual species were also compared between thickets and no L. camara. At the community level, L. camara density, dispersal mode and forest habitat affinities of species determined both frequency of occurrence and abundance of seedlings, with the abundance of dry-forest mammal-dispersed species and ubiquitous mechanically dispersed species being significantly lower under L. camara thickets. Phyllanthus emblica and Kydia calycina were found to be significantly less abundant under L. camara, whereas most other species were not affected by the presence of thickets. It was inferred that, by affecting the establishment of native tree seedlings, L. camara thickets could eventually alter the community composition of such forests.

Morphology and thermal characteristics of nano-sized Pb-Sn inclusions in Al

Nanoembedded aluminum alloys with bimetallic dispersoids of Sn and Pb of compositions Sn-82-Pb-18,Pb- Sn-64-Pb-36, and Sn-54-Pb-46 were synthesized by rapid solidification. The two phases, face-centered-cubic Pb and tetragonal Sn solid-solution, coexist in all the particles. The crystallographic relation between the two phases and the matrix depends upon the solidification pathways adopted by the particles. For Al-(Sn-82-Pb-18), we report a new orientation relation given by [011]Al//[010]Sn and (011)Al//(101)Sn. Pb exhibits a cube-on-cube orientation with Al in few particles, while in others no orientation relationship could be observed. In contrast, Pb in Sn-64-Pb-36 and Sn-54-Pb-46 particles always exhibits cube-on-cube orientation with the matrix. Sn does not show any orientation relationship with Al or Pb in these cases. Differential scanning calorimetry studies revealed melting at eutectic temperature for all compositions, although solidification pathways are different. Attempts were made to correlate these with the melting and heterogeneous nucleation. characteristics.

Performance analysis of the random early detection algorithm

In this article we consider a finite queue with its arrivals controlled by the random early detection algorithm. This is one of the most prominent congestion avoidance schemes in the Internet routers. The aggregate arrival stream from the population of transmission control protocol sources is locally considered stationary renewal or Markov modulated Poisson process with general packet length distribution. We study the exact dynamics of this queue and provide the stability and the rates of convergence to the stationary distribution and obtain the packet loss probability and the waiting time distribution. Then we extend these results to a two traffic class case with each arrival stream renewal. However, computing the performance indices for this system becomes computationally prohibitive. Thus, in the latter half of the article, we approximate the dynamics of the average queue length process asymptotically via an ordinary differential equation. We estimate the error term via a diffusion approximation. We use these results to obtain approximate transient and stationary performance of the system. Finally, we provide some computational examples to show the accuracy of these approximations.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.