Obliquely incident water waves against a vertical cliff

Utilizing the commutativity property of the Cartesian coordinate differential operators arising in the boundary conditions associated with the propagation of surface water waves against a vertical cliff, under the assumptions of linearized theory, the problem of obliquely incident surface waves is considered for solution. The case of normal incidence, handled by previous workers follow as a particular limiting case of the present problem, which exhibits a source/sink type behavior of the velocity potential at the shore-line. An independent method of attack is also presented to handle the case of normal incidence.

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

Atomic vibration in the Carbon Nanotubes (CNTs) gives rise to non-local interactions. In this paper, an expression for the non-local scaling parameter is derived as a function of the geometric and electronic properties of the rolled graphene sheet in single-walled CNTs. A self-consistent method is developed for the linearization of the problem of ultrasonic wave propagation in CNTs. We show that (i) the general three-dimensional elastic problem leads to a single non-local scaling parameter (e(0)), (ii) e(0) is almost constant irrespective of chirality of CNT in the case of longitudinal wave propagation, (iii) e(0) is a linear function of diameter of CNT for the case of torsional mode of wave propagation, (iv) e(0) in the case of coupled longitudinal-torsional modes of wave propagation, is a function which exponentially converges to that of axial mode at large diameters and to torsional mode at smaller diameters. These results are valid in the long-wavelength limit. (C) 2011 Elsevier Ltd. All rights reserved.

Estimation of low-velocity impact damage in laminated composite circular plates using nonlinear finite element analysis

Nonlinear finite element analysis is used for the estimation of damage due to low-velocity impact loading of laminated composite circular plates. The impact loading is treated as an equivalent static loading by assuming the impactor to be spherical and the contact to obey Hertzian law. The stresses in the laminate are calculated using a 48 d.o.f. laminated composite sector element. Subsequently, the Tsai-Wu criterion is used to detect the zones of failure and the maximum stress criterion is used to identify the mode of failure. Then the material properties of the laminate are degraded in the failed regions. The stress analysis is performed again using the degraded properties of the plies. The iterative process is repeated until no more failure is detected in the laminate. The problem of a typical T300/N5208 composite [45 degrees/0 degrees/-45 degrees/90 degrees](s) circular plate being impacted by a spherical impactor is solved and the results are compared with experimental and analytical results available in the literature. The method proposed and the computer code developed can handle symmetric, as well as unsymmetric, laminates. It can be easily extended to cover the impact of composite rectangular plates, shell panels and shells.

Translational velocity oscillations of piston generated vortex rings

Experimental results are presented that show that the translational velocities of piston generated vortex rings often undergo oscillations, similar to those recently discovered for drop generated rings. An attempt has been made to minimize uncertainties by utilizing both dye and hydrogen bubbles for visualization and carefully repeating measurements on the same ring and on different realizations under the same nominal piston conditions. The results unambiguously show that under most conditions, both for laminar and turbulent rings and for rings generated from pipes and orifices, the oscillations are present. The present results, together with the earlier results on drop generated rings, give support to the view that translational velocity oscillations are probably an inherent feature of translating vortex ring fields. (C) 1995 American Institute of Physics.

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.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Damped Belief Propagation Based Near-Optimal Equalization of Severely Delay-Spread UWB MIMO-ISI Channels

In this paper, we present a belief propagation (BP) based equalizer for ultrawideband (UWB) multiple-input multiple-output (MIMO) inter-symbol interference (ISI) channels characterized by severe delay spreads. We employ a Markov random field (MRF) graphical model of the system on which we carry out message passing. The proposed BP equalizer is shown to perform increasingly closer to optimal performance for increasing number of multipath components (MPC) at a much lesser complexity than that of the optimum equalizer. The proposed equalizer performs close to within 0.25 dB of SISO AWGN performance at 10-3 bit error rate on a severely delay-spread MIMO-ISI channel with 20 equal-energy MPCs. We point out that, although MIMO/UWB systems are characterized by fully/densely connected graphical models, the following two proposed features are instrumental in achieving near-optimal performance for large number of MPCs at low complexities: i) use of pairwise compatibility functions in densely connected MRFs, and ii) use of damping of messages.

Improved large-MIMO detection based on damped belief propagation

In this paper, we consider the application of belief propagation (BP) to achieve near-optimal signal detection in large multiple-input multiple-output (MIMO) systems at low complexities. Large-MIMO architectures based on spatial multiplexing (V-BLAST) as well as non-orthogonal space-time block codes(STBC) from cyclic division algebra (CDA) are considered. We adopt graphical models based on Markov random fields (MRF) and factor graphs (FG). In the MRF based approach, we use pairwise compatibility functions although the graphical models of MIMO systems are fully/densely connected. In the FG approach, we employ a Gaussian approximation (GA) of the multi-antenna interference, which significantly reduces the complexity while achieving very good performance for large dimensions. We show that i) both MRF and FG based BP approaches exhibit large-system behavior, where increasingly closer to optimal performance is achieved with increasing number of dimensions, and ii) damping of messages/beliefs significantly improves the bit error performance.

A study on determining the theoretical dispersion curve for Rayleigh wave propagation

A method has been presented to establish the theoretical dispersion curve for performing the inverse analysis for the Rayleigh wave propagation. The proposed formulation is similar to the one available in literature, and is based on the finite difference formulation of the governing partial differential equations of motion. The method is framed in such a way that it ultimately leads to an Eigen value problem for which the solution can be obtained quite easily with respect to unknown frequency. The maximum absolute value of the vertical displacement at the ground surface is formed as the basis for deciding the governing mode of propagation. With the proposed technique, the numerical solutions were generated for a variety of problems, comprising of a number of different layers, associated with both ground and pavements. The results are found to be generally satisfactory. (C) 2011 Elsevier Ltd. All rights reserved.

Energy released during arc crack propagation

The problem of circular arc cracks in a homogeneous medium is revisited. An unusual but simple method to calculate the energy change due to arc crack propagation along a circle is illustrated based on the earlier work of Sih and Liebowitz (1968). The limiting case of crack of angle 27pi is shown to correspond with the problem of a circular hole in a large plate under remote loading.

Incremental-angle and angular velocity estimation using a star sensor

A method is described for estimating the incremental angle and angular velocity of a spacecraft using integrated rate parameters with the help of a star sensor alone. The chief advantage of this method is that the measured stars need not be identified, whereas the identification of the stars is necessary in earlier methods. This proposed estimation can be carried out with all of the available measurements by a simple linear Kalman filter, albeit with a time-varying sensitivity matrix. The residuals of estimated angular velocity by the proposed spacecraft incremental-angle and angular velocity estimation method are as accurate as the earlier methods. This method also enables the spacecraft attitude to be reconstructed for mapping the stars into an imaginary unit sphere in the body reference frame, which will preserve the true angular separation of the stars. This will pave the way for identification of the stars using any angular separation or triangle matching techniques applied to even a narrow field of view sensor that is made to sweep the sky. A numerical simulation for inertial as well as Earth pointing spacecraft is carried out to establish the results.

MHD flow over a moving plate in a rotating fluid with magnetic field, Hall currents and free stream velocity

The non-similar boundary layer flow of a viscous incompressible electrically conducting fluid over a moving surface in a rotating fluid, in the presence of a magnetic field, Hall currents and the free stream velocity has been studied. The parabolic partial differential equations governing the flow are solved numerically using an implicit finite-difference scheme. The Coriolis force induces overshoot in the velocity profile of the primary flow and the magnetic field reduces/removes the velocity overshoot. The local skin friction coefficient for the primary flow increases with the magnetic field, but the skin friction coefficient for the secondary flow reduces it. Also the local skin friction coefficients for the primary and secondary flows are reduced due to the Hall currents. The effects of the magnetic field, Hall currents and the wall velocity, on the skin friction coefficients for the primary and secondary flows increase with the Coriolis force. The wall velocity strongly affects the flow field. When the wall velocity is equal to the free stream velocity, the skin friction coefficients for the primary and secondary flows vanish, but this does not imply separation. (C) 2002 Published by Elsevier Science Ltd.

A nonlinear spectral finite element model for analysis of wave propagation in solid with internal friction and dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.

Spherical piston problem in water

In this paper, we study the propagation of a shock wave in water, produced by the expansion of a spherical piston with a finite initial radius. The piston path in the x, t plane is a hyperbola. We have considered the following two cases: (i) the piston accelerates from a zero initial velocity and attains a finite velocity asymptotically as t tends to infinity, and (ii) the piston decelerates, starting from a finite initial velocity. Since an analytic approach to this problem is extremely difficult, we have employed the artificial viscosity method of von Neumann & Richtmyer after examining its applicability in water. For the accelerating piston case, we have studied the effect of different initial radii of the piston, different initial curvatures of the piston path in the x, t plane and the different asymptotic speeds of the piston. The decelerating case exhibits the interesting phenomenon of the formation of a cavity in water when the deceleration of the piston is sufficiently high. We have also studied the motion of the cavity boundary up to 550 cycles.