Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x greater-or-equal, slanted 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as var epsilon → 0 is obtained.

Pulsatile flow of a viscous fluid in elliptical pipe of variable cross-section

The pulsatile flow of an incompressible viscous fluid in an elliptical pipe of slowly varying cross-section is considered. Asymptotic series solutions for the velocity distribution and pressure gradient are obtained in terms of Mathieu functions for a low Reynold number flow in which the volume flux is prescribed. An expression for shear stress on the boundary is derived. The physically significant quantities governing the flow are computed numerically and analysed for different types of constrictions. The effect of eccentricity and Womerslay parameter on the flow is discussed.

Non-linearity in piezo-fiber reinforced composites: an asymptotic approach

An asymptotically correct analysis is developed for Macro Fiber Composite unit cell using Variational Asymptotic Method (VAM). VAM splits the 3D nonlinear problem into two parts: A 1D nonlinear problem along the length of the fiber and a linear 2D cross-sectional problem. Closed form solutions are obtained for the 2D problem which are in terms of 1D parameters.

The scratch hardness and friction of a soft rigid-plastic solid

The paper describes an experimental and analytical study of the normal and scratch hardnesses of a model soft rigid-plastic solid. The material known as ‘Plasticine’, a mixture of dry particles and a mineral oil, has been deformed with a range of rigid conical indentors with included angles of between 30° and 170°. The sliding velocity dependence of the computed scratch hardness and friction has been examined in the velocity range 0.19 mm/s to 7.3 m/s. Data are also described for the time dependence of the normal hardness and also the estimated rate dependence of the intrinsic flow stress. The latter values were estimated from data obtained during the upsetting of right cylinders. Three major conclusions are drawn from these data and the associated analysis. (1) A first-order account of the scratching force may be provided by adopting a model which sums the computed plastic deformation and interfacial sliding contributions to the total sliding work. This is tantamount to the adoption of the two-term non-interacting model of friction. (2) For this system during sliding, at high sliding velocities at least, the interface shear stress which defines the boundary condition is not directly related to the bulk shear stress. The interface rheological characteristics indicate an appreciable dependence on the imposed strain or strain rate. In particular, the relative contributions of the slip and stick boundary conditions appear to be a function of the imposed sliding velocity. (3) The computed normal and scratch hardness values are not simply interrelated primarily because of the evolving boundary conditions which appear to exist in the scratching experiments.

An actor-critic algorithm with function approximation for discounted cost constrained Markov decision processes

We develop in this article the first actor-critic reinforcement learning algorithm with function approximation for a problem of control under multiple inequality constraints. We consider the infinite horizon discounted cost framework in which both the objective and the constraint functions are suitable expected policy-dependent discounted sums of certain sample path functions. We apply the Lagrange multiplier method to handle the inequality constraints. Our algorithm makes use of multi-timescale stochastic approximation and incorporates a temporal difference (TD) critic and an actor that makes a gradient search in the space of policy parameters using efficient simultaneous perturbation stochastic approximation (SPSA) gradient estimates. We prove the asymptotic almost sure convergence of our algorithm to a locally optimal policy. (C) 2010 Elsevier B.V. All rights reserved.

Saturation throughput analysis of an input queueing ATM switch with multiclass bursty traffic

In this paper we consider an N x N non-blocking, space division ATM switch with input cell queueing. At each input, the cell arrival process comprises geometrically distributed bursts of consecutive cells for the various outputs. Motivated by the fact that some input links may be connected to metropolitan area networks, and others directly to B-ISDN terminals, we study the situation where there are two classes of inputs with different values of mean burst length. We show that when inputs contend for an output, giving priority to an input with smaller expected burst length yields a saturation throughput larger than if the reverse priority is given. Further, giving priority to less bursty traffic can give better throughput than if all the inputs were occupied by this less bursty traffic. We derive the asymptotic (as N --> infinity) saturation throughputs for each priority class.

Shear-induced enhancement of self-diffusion in interacting colloidal suspensions

We set up the generalized Langevin equations describing coupled single-particle and collective motion in a suspension of interacting colloidal particles in a shear how and use these to show that the measured self-diffusion coefficients in these systems should be strongly dependent on shear rate epsilon. Three regimes are found: (i) an initial const+epsilon(.2), followed by (ii) a large regime of epsilon(.1/2) behavior, crossing over to an asymptotic power-law approach (iii) D-o - const x epsilon(.-1/2) to the Stokes-Einstein value D-o. The shear dependence is isotropic up to very large shear rates and increases with the interparticle interaction strength. Our results provide a straightforward explanation of recent experiments and simulations on sheared colloids.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Recovery of tropical rainforest avifauna in relation to vegetation succession following shifting cultivation in Mizoram, north-east India

1. Recovery of rainforest bird community structure and composition, in relation to forest succession after slash-and-burn shifting cultivation or jhum was studied in Mizoram, north-east India. Replicate fallow sites abandoned after shifting cultivation 1, 5, 10, 25 and approximate to 100 years ago, were compared with primary evergreen and semi-evergreen forest using transect and quadrat sampling. 2. Vegetation variables such as woody plant species richness, tree density and vertical stratification increased with fallow age in a rapid. nun-linear, asymptotic manner. Principal components analysis of vegetation variables summarized 92.8% of the variation into two axes: PC1 reflecting forest development and woody plant succession (variables such as tree density, woody plant species richness), and PC2 depicting bamboo density, which increased from 1 to 25 years and declined thereafter. 3. Bird species richness, abundance and diversity, increased rapidly and asymptotically during succession paralleling vegetation recovery as shown by positive correlations with fallow age and PC1 scores of sites. Bamboo density reflected by PC2 had a negative effect on bird species richness and abundance. 4. The bird community similarity (Morisita index) of sites with primary forest also increased asymptotically with fallow age indicating sequential species turnover during succession. Bird community similarity of sites with primary forest (or between sites) was positively correlated with both physiognomic and floristic similarities with primary forest (or between sites). 5. The number of bird species in guilds associated with forest development and woody plants (canopy insectivores, frugivores: bark feeders) was correlated with PCI scores of the sites. Species in other guilds (e. g. granivores, understorey insectivores) appeared to dominate during early and mid-succession. 6. The non-linear relationships imply that fallow periods less than a threshold of 25 years for birds, and about 50-75 years for woody plants, are likely to cause substantial community alteration. 7. As 5-10-year rotation periods or jhum cycles prevail in many parts of north-east India. there is a need to protect and conserve tracts of late-successional and primary forest.

First-principle description of magnonic PdnFem multilayers

Ab-initio calculations are used to determine the parameters that determine magnonic band structure of PdnFem multilayers (n = 2, m <= 8). We obtain the layer-resolved magnetization, the exchange coupling, and the magnetic anisotropy of the Pd-Fe structures. The Fe moment is 3.0 mu(B) close to the Pd layers and 2.2 mu(B) in the middle of the Fe layers. An intriguing but not usually considered aspect is that the elemental Pd is nonmagnetic, similar to Cu spacer layers in other multilayer systems. This leads to a pre-asymptotic ferromagnetic coupling through the Pd (about 40 mJ/m(2)). Furthermore, the Pd acquires a small moment due to spin polarization by neighboring Fe atoms, which translates into magnetic anisotropy. The anisotropies are large, in the range typical for L1(0) structures, which is beneficial for high-frequency applications. (C) 2011 American Institute of Physics. doi:10.1063/1.3556763]

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.

Quasi-Orthogonal Design and Performance Analysis of Amplify-And-Forward Relay Networks with Multiple-Antennas

This paper is on the design and performance analysis of practical distributed space-time codes for wireless relay networks with multiple antennas terminals. The amplify-andforward scheme is used in a way that each relay transmits a scaled version of the linear combination of the received symbols. We propose distributed generalized quasi-orthogonal space-time codes which are distributed among the source antennas and relays, and valid for any number of relays. Assuming M-PSK and M-QAM signals, we derive a formula for the symbol error probability of the investigated scheme over Rayleigh fading channels. For sufficiently large SNR, this paper derives closed-form average SER expression. The simplicity of the asymptotic results provides valuable insights into the performance of cooperative networks and suggests means of optimizing them. Our analytical results have been confirmed by simulation results, using full-rate full-diversity distributed codes.

Effect of longitudinal surface curvature on boundary layers

We consider here the higher order effect of moderate longitudinal surface curvature on steady, two-dimensional, incompressible laminar boundary layers. The basic partial differential equations for the problem, derived by the method of matched asymptotic expansions, are found to possess similarity solutions for a family of surface curvatures and pressure gradients. The similarity equations obtained by this anaylsis have been solved numerically on a computer, and show a definite decrease in skin friction when the surface has convex curvature in all cases including zero pressure gradient. Typical velocity profiles and some relevant boundary-layer characteristics are tabulated, and a critical comparison with previous work is given.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Steady compressible flow of cohesionless granular materials through a wedge-shaped bunker

A continuum model based on the critical-state theory of soil mechanics is used to generate stress, density, and velocity profiles, and to compute discharge rates for the flow of granular material in a mass flow bunker. The bin–hopper transition region is idealized as a shock across which all the variables change discontinuously. Comparison with the work of Michalowski (1987) shows that his experimentally determined rupture layer lies between his prediction and that of the present theory. However, it resembles the former more closely. The conventional condition involving a traction-free surface at the hopper exit is abandoned in favour of an exit shock below which the material falls vertically with zero frictional stress. The basic equations, which are not classifiable under any of the standard types, require excessive computational time. This problem is alleviated by the introduction of the Mohr–Coulomb approximation (MCA). The stress, density, and velocity profiles obtained by integration of the MCA converge to asymptotic fields on moving down the hopper. Expressions for these fields are derived by a perturbation method. Computational difficulties are encountered for bunkers with wall angles θw [gt-or-equal, slanted] 15° these are overcome by altering the initial conditions. Predicted discharge rates lie significantly below the measured values of Nguyen et al. (1980), ranging from 38% at θw = 15° to 59% at θw = 32°. The poor prediction appears to be largely due to the exit condition used here. Paradoxically, incompressible discharge rates lie closer to the measured values. An approximate semi-analytical expression for the discharge rate is obtained, which predicts values within 9% of the exact (numerical) ones in the compressible case, and 11% in the incompressible case. The approximate analysis also suggests that inclusion of density variation decreases the discharge rate. This is borne out by the exact (numerical) results – for the parameter values investigated, the compressible discharge rate is about 10% lower than the incompressible value. A preliminary comparison of the predicted density profiles with the measurements of Fickie et al. (1989) shows that the material within the hopper dilates more strongly than predicted. Surprisingly, just below the exit slot, there is good agreement between theory and experiment.