Wave Packet Critical Layers in Shear Flows

In the linear inviscid theory of shear flow stability, the eigenvalue problem for a neutral or weakly amplified mode revolves around possible discontinuities in the eigenfunction as the singular critical point is crossed. Extensions of the linear normal mode approach to include nonlinearity and/or wave packets lead to amplitude evolution equations where, again, critical point singularities are an issue because the coefficients of the amplitude equations generally involve singular integrals. In the past, viscosity, nonlinearity, or time dependence has been introduced in a critical layer centered upon the singular point to resolve these integrals. The form of the amplitude evolution equation is greatly influenced by which choice is made. In this paper, a new approach is proposed in which wave packet effects are dominant in the critical layer and it is argued that in many applications this is the appropriate choice. The theory is applied to two-dimensional wave propagation in homogeneous shear flows and also to stratified shear flows. Other generalizations are indicated.     

Uniform asymptotic analysis for waves in an incompressible elastic rod II. Disturbances superimposed on a pre-stressed state

This paper studies the propagation of disturbances superimposedon a pre-stressed incompressible hyperelastic thin rod. Startingfrom the incremental equations given by Haughton and Ogden (1979,J. Mech. Phys. Solids 27, 179-212 and 489-512), we derive athree parameter-dependent one-dimensional rod equation as thegoverning equation. In particular, it is found that one parameterplays a crucial role. Depending on whether it is larger or smallerthan or equal to a critical value, the shear-wave velocity islarger or smaller than or equal to the bar-wave velocity. Inthe case that these two velocities are equal, there exist travelling-wavesolutions of arbitrary form. This implies that for this particularcase the initial disturbance would propagate along the rod withoutdistortion. To see the influence of the pre-stress in detail,we further consider an initial-value problem with an initialsingularity in the shear strain. The solutions are expressedin terms of integrals through the method of Fourier transform.We then conduct an asymptotic analysis for the solutions. Fora material point in a neighbourhood behind the shear-wave front,the phase function of these integrals has a stationary pointat infinity. Here, we use a technique of uniform asymtotic expansionto handle this case. An asymptotic expansion, correct up toorder O(t^-1), for the shear strain, which is uniformly validin a neighbourhood behind the shear-wave front, is derived.For material points in other spatial domains, the method ofstationary point is applicable, and asymptotic expansions (correctup to order O(t^-1)) are obtained. A novelty is that we are ableto deduce precise qualitative information about the waves inthe far field from our analytic results. Wave profiles for twoconcrete examples are also provided.     

Numerical Simulation of Red Blood Cell Suspensions Behind a Moving Interface in a Capillary

Computational modeling and simulation are presented on the motion of red blood cells behind a moving interface in a capillary. The methodology is based on an immersed boundary method and the skeleton structure of the red blood cell (RBC) membrane is modeled as a spring network. As by the nature of the problem, the computational domain is moving with either a designated RBC or an interface in an infinitely long two-dimensional channel with an undisturbed flow field in front of the computational domain. The tanking-treading and the inclination angle of a cell in a simple shear flow are briefly discussed for the validation purpose. We then present and discuss the results of the motion of red blood cells behind a moving interface in a capillary, which show that the RBCs with higher velocity than the interface speed form a concentrated slug behind the moving interface.     

An estimate of the attraction domain with a specified exponential stability index in a wheeled robot control problem

The problem of synthesizing a law for the control of the plane motion of a wheeled robot is investigated. The rear wheels are the drive wheels and the front wheels are responsible for the turning of the platform. The aim of the control is to steer a target point to a specified trajectory and to stabilize the motion along it. The trajectory is assumed to be specified by a smooth curve. The actual curvature of the trajectory of the target point, which is related to the angle of rotation of the front wheels by a simple algebraic relation, is considered as the control. The control is subjected to bilateral constraints by virtue of the fact that the angle of rotation of the front wheels is limited. The attraction domain in the distance to trajectory - orientation space, is investigated for the proposed control law. Arrival at a trajectory with a specified exponential stability index is guaranteed in the case of initial conditions belonging to the given domain. An estimate of the attraction domain in the form of an ellipse is given.     

Expected runtimes of a simple evolutionary algorithm for the multi-objective minimum spanning tree problem

Evolutionary algorithms are applied to problems that are not well understood as well as to problems in combinatorial optimization. The analysis of these search heuristics has been started for some well-known polynomial solvable problems. Such analyses are starting points for the analysis of evolutionary algorithms on difficult problems. We present the first runtime analysis of a multi-objective evolutionary algorithm on a NP-hard problem. The subject of our analysis is the multi-objective minimum spanning tree problem for which we give upper bounds on the expected time until a simple evolutionary algorithm has produced a population including for each extremal point of the Pareto front a corresponding spanning tree. These points are of particular interest as they give a 2-approximation of the Pareto front. We show that in expected pseudopolynomial time a population is produced that includes for each extremal point a corresponding spanning tree.     

The harmonic oscillations of a viscoelastic rod of triangular cross-section

Two exact solutions of the plane strain problem of the harmonic oscillations of a viscoelastic rod, the cross-section of which is a right triangle, are proposed. Either the normal displacement and the shear stress or the shear displacement and the normal stress of the side surface of the rod are given. Six dimensionless parameters which affect the dynamic deformation process are derived. Two parameters characterize the contribution of the viscous properties with respect to the elastic properties, two others define the logarithmic decrement of the longitudinal and shear harmonic waves, and two other parameters affect the wavelength of the corresponding wave and the velocity of motion of the wave front of these waves. The velocities of both types of waves and their wavelengths turn out to be greater than the velocities and wavelengths of the corresponding elastic waves. It is shown that, for certain values of the viscosity and the oscillation frequency, pseudo-resonance frequencies are possible which are higher than the resonance frequencies for an elastic medium.     

Propagation of Gabor singularities for Schrödinger equations with quadratic Hamiltonians

We study propagation of the Gabor wave front set for a Schrödinger equation with a Hamiltonian that is the Weyl quantization of a quadratic form with nonnegative real part. We point out that the singular space associated with the quadratic form plays a crucial role for the understanding of this propagation. We show that the Gabor singularities of the solution to the equation for positive times are always contained in the singular space, and that they propagate in this set along the flow of the Hamilton vector field associated with the imaginary part of the quadratic form. As an application we obtain for the heat equation a sufficient condition on the Gabor wave front set of the initial datum tempered distribution that implies regularization to Schwartz regularity for positive times.     

A restart algorithm for computing fixed points without an extra dimension

An algorithm to compute a fixed point of an upper semicontinuous point to set mapping using a simplicial subdivision is introduced. The new element of the algorithm is that for a given grid it does not start with a subsimplex but with one (arbitrary) point only; the algorithm will terminate always with a subsimplex. This subsimplex yields an approximation of a fixed point and provides the starting point for a finer grid. Some numerical results suggest that this algorithm converges more rapidly than the known algorithms. Moreover, it is very simple to implement the algorithm on the computer.     

Stress and fracture analysis in delaminated orthotropic composite plates using third-order shear deformation theory

The third-order shear deformable plate theory is applied in this work to calculate the stresses and energy release rates in delaminated orthotropic composite plates with straight crack front. The delaminated parts are modeled by the general third-order plate theory, while a double-plate model with interface constraint is developed for the uncracked portion of the plate. The governing equations of the uncracked part are formulated by considering the equilibrium and the displacement continuity along the interface. As an example, a simply-supported delaminated orthotropic plate subjected to a point force is solved adopting Lévy plate formulation and the state-space approach. The mode-II and mode-III energy release rate distributions along the crack front were calculated by the J-integral. To verify the analytical results the 3D finite element model of the plate was constructed and the energy release rates were calculated by the virtual crack-closure technique. A previous second-order plate theory solution was also utilized in the course of the comparison. The results indicate a good agreement between analysis and numerical computation and that third-order theory is better in some cases than the second-order approximation.     

Magnetohydrodynamic shear flow along a flat plate with uniform suction or blowing

An analysis is made of the steady shear flow of an incompressible viscous electrically conducting fluid past an electrically insulating porous flat plate in the presence of an applied uniform transverse magnetic field. It is shown that steady shear flow exists for suction at the plate only when the square of the suction parameter S is less than the magnetic parameter Q. In this case the velocity at a given point increases with increase in either the magnetic field or suction velocity. The shear stress at the plate increases with increase in either S or the free-stream shear-rate parameter σ₁ or Q. The analysis further reveals that solution exists for steady shear flow past a porous flat plate subject to blowing only when the square of the blowing parameter S₁ is less than Q. It is found that the induced magnetic field at a given location decreases with increase in Q. Further the wall shear stress decreases with increase in S₁. No steady shear flow is possible for blowing at the plate when S₁² > Q. Received: June 16, 2004; revised: October 24, 2004     

Magnetic field effects on Newtonian and non-Newtonian ferrofluid flow past a circular cylinder

The changes in the flow properties under the action of electromagnetic body forces are investigated numerically for ferrofluid flow past a circular cylinder. Ferrofluid is modeled as both a Newtonian and a non-Newtonian Power-Law fluid. Magnetic forces are applied by placing magnets at different locations on the surface of the cylinder. The magnetostatic effects on the structure of the wake region, on drag reduction and on vortex formation length and frequency are shown and compared in terms of Reynolds number, interaction parameter, Power-Law index and magnet location. It is shown that the increase in the interaction parameter reduces drag for both Newtonian and non-Newtonian model. This decrease is observed to be higher for shear thinning and lower for shear thickening fluid compared to Newtonian case. It is also shown that vortex street formation in the wake region behind the cylinder may be delayed under high magnetic effects. The Strouhal number is higher for shear thinning case at both low and high Reynolds numbers, and lower for shear thickening case at high Reynolds numbers, compared to Newtonian fluid. The vortex formation frequency also decreases under the action of the magnetic field in all cases, however the vortex formation length increases. Placing the magnet towards the front region of the cylinder increases considerably the drag coefficient for both Newtonian and non-Newtonian model. This increase in drag coefficient is higher in the shear thinning fluid and lower in the shear thickening fluid compared to the Newtonian case.     

基于非负矩阵理论的同步网络AIMD算法分析

利用非负矩阵理论提出了一个多源竞争共享带宽的网络模型.研究发现,满足模型假设的网络通常都有一个唯一的、全局指数收敛的稳态点.利用这个模型,分析了全局网络收敛到稳态点的收敛速度和稳定性.     

Cost Functions in the Location of Depots for Multiple-Delivery Journeys

Most wholesale distribution is performed during multiple-delivery journeys. Mathematical methods of locating depots utilize simple functions of delivery data, e.g. weight and distance from the depot, to measure the delivery "cost"; the total "cost" is minimized to find the depot location.It is pointed out that the cost of a delivery is influenced by the occurrence of other deliveries. It is shown that, in a few examples, simple functions of the delivery data are not always good measures of variable cost, as measured by the length of journeys planned to carry out the same deliveries, and that the minimum points of the simple functions rarely coincide with the point of minimum variable cost. It is concluded that, subject to reservations, which are discussed, about the experiments, the use of simple cost functions to locate wholesale distribution depots will probably give misleading results.     

Front propagation in periodic excitable media

This paper is devoted to the study of pulsating travelling fronts for reaction‐diffusion‐advection equations in a general class of periodic domains with underlying periodic diffusion and velocity fields. Such fronts move in some arbitrarily given direction with an unknown effective speed. The notion of pulsating travelling fronts generalizes that of travelling fronts for planar or shear flows. Various existence, uniqueness and monotonicity results are proved for two classes of reaction terms. Firstly, for a combustion‐type nonlinearity, it is proved that the pulsating travelling front exists and that its speed is unique. Moreover, the front is increasing with respect to the time variable and unique up to translation in time. We also consider one class of monostable nonlinearity which arises either in combustion or biological models. Then, the set of possible speeds is a semi‐infinite interval, closed and bounded from below. For each possible speed, there exists a pulsating travelling front which is increasing in time. This result extends the classical Kolmogorov‐Petrovsky‐Piskunov case. Our study covers in particular the case of flows in all of space with periodic advections such as periodic shear flows or a periodic array of vortical cells. These results are also obtained for cylinders with oscillating boundaries or domains with a periodic array of holes. © 2002 Wiley Periodicals, Inc.     

An Erratum in the Paper “A Plane Shear Wave Excited in an Elastic Half-Space by Oscillating Normal Pressure Distributed Nonuniformly”

Excitation of an elastic half-space by an oscillating normal load nonuniformly distributed over its surface is considered. The medium is assumed to be linear, isotropic, and homogeneous. Exterior load distribution is smooth, and the scale of its variation is small compared to wavelengths. A simple expression that describes the shear wave with plane front set, running away from the surface, is found. This wave vanishes if the load is uniform. Bibliography: 1 title.     

The Stability of Attached Shock Waves in Magnetogasdynamic Flow Past a Wedge

The question of the physical significance of the new phenomenaindicated by Cabannes' work on magnetogasdynamic flow past awedge is considered from the point of view of the stabilityof the shock waves. Analytical and heuristic reasons are givensuggesting that downstream-facing shocks are stable if the upstreamflow is supersonic and unstable if it is subsonic, while upstreamfacing shocks are always to be considered unstable.     

A time integration procedure for plastic deformation in elastic-viscoplastic metals

A simple unconditionally stable numerical procedure for time integration of the flow rule for large plastic deformation of an elastic-viscoplastic metal is developed. Specific attention is focused on a unified set of constitutive equations which represents a generalization (for large deformation and thermomechanical response) of the Bodner-Partom model [6, 7]. An analytical solution is obtained for large deformation simple shear at constant shear rate. Numerical examples of simple shear, a corner test exhibiting the transition from uniaxial compression to shear, and simple tension are considered which demonstrate the stability and accuracy of the procedure. It is shown that the same procedure can be used for a rate insensitive metal characterized by a yield function as well as for a rate sensitive metal characterized by an overstress model. Finally, an appendix is provided which records the basic equations associated with the small deformation theory.     

Properties of simple model problems for reacting gas flows

The compressible Navier–Stokes equations for reacting gases are extremely complex. Simpler models have been considered, and for these completely non-physical propagation speeds have been observed. These model problems are stiff, meaning that several different scales are present in the solution. Numerical solution of non-reacting flows almost always involves addition of extra dissipation. It will be shown that this action will render a totally wrong propagation speed for a simple model equation of reacting flows. This problem will be accentuated by increasing stiffness of the problem. Existence and uniqueness of a solution to this model equation is proved. The dependence of the propagation speed on the viscosity and a term governing the stiffness (comparable to the reaction rate for a more complete model) is investigated. A remedy for the wrong propagation speed for this simple model equation is proposed such that the speed is correct although the front is smeared out.     

A simplified kinematic wave model at a merge bottleneck

A merge is a point of a highway where two or more streams of traffic flow into one. It is always easy to solve the demand–supply problem at a merge when the merge is operating under uncongested condition. However, when congestion backs up exceeding the merging point where multiple streams of traffic meet, one is typically faced with splitting downstream supply among the merging branches. Solutions of this situation are multiple and several merge queuing models have been proposed in the past. To address the drawbacks of the past models, this paper proposes a capacity-based weighted fair queuing (CBWFQ) model that is characterized by its fidelity (approximation to real situation), simplicity, and extensibility. Based on the CBWFQ merge queuing model, a simplified kinematic wave model is formulated to model traffic operation at a merge bottleneck.