On the Selection of Optimum Savitzky-Golay Filters

Savitzky-Golay (S-G) filters are finite impulse response lowpass filters obtained while smoothing data using a local least-squares (LS) polynomial approximation. Savitzky and Golay proved in their hallmark paper that local LS fitting of polynomials and their evaluation at the mid-point of the approximation interval is equivalent to filtering with a fixed impulse response. The problem that we address here is, ``how to choose a pointwise minimum mean squared error (MMSE) S-G filter length or order for smoothing, while preserving the temporal structure of a time-varying signal.'' We solve the bias-variance tradeoff involved in the MMSE optimization using Stein's unbiased risk estimator (SURE). We observe that the 3-dB cutoff frequency of the SURE-optimal S-G filter is higher where the signal varies fast locally, and vice versa, essentially enabling us to suitably trade off the bias and variance, thereby resulting in near-MMSE performance. At low signal-to-noise ratios (SNRs), it is seen that the adaptive filter length algorithm performance improves by incorporating a regularization term in the SURE objective function. We consider the algorithm performance on real-world electrocardiogram (ECG) signals. The results exhibit considerable SNR improvement. Noise performance analysis shows that the proposed algorithms are comparable, and in some cases, better than some standard denoising techniques available in the literature.

Time dependent superhydrophobicity of drag reducing surfaces

Air can be trapped on the crevices of specially textured hydrophobic surfaces immersed in water. This heterogenous state of wetting in which the water is in contact with both the solid surface and the entrapped air is not stable. Diffusion of air into the surrounding water leads to gradual reduction in the size and numbers of the air bubbles. The sustainability of the entrapped air on such surfaces is important for many underwater applications in which the surfaces have to remain submersed for longer time periods. In this paper we explore the suitability of different classes of surface textures towards the drag reduction application by evaluating the time required for the disappearance of the air bubbles under hydrostatic conditions. Different repetitive textures consisting of holes, pillars and ridges of different sizes have been generated in silicon, aluminium and brass by isotropic etching, wire EDM and chemical etching respectively. These surfaces were rendered hydrophobic with self-assembled layer of fluorooctyl trichlorosilane for silicon and aluminium surfaces and 1-dodecanethiol for brass surfaces. Using total internal reflection the air bubbles are visualized with the help of a microscope and time lapse photography. Irrespective of the texture, both the size and the number of air pockets were found to decrease with time gradually and eventually disappear. In an attempt to reverse the diffusion we explore the possibility of using electrolysis to generate gases at the textured surfaces. The gas bubbles are nucleated everywhere on the surface and as they grow they coalesce with each other and get pinned at the texture edges.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Transient Momentum and Enthalpy Transfer in Packed Beds at High Reynolds Numbers

An experimental study for transient temperature response and pressure drop in a randomly packed bed at high Reynolds numbers is presented.The packed bed is used as a compact heat exchanger along with a solid-propellant gas generator, to generate room-temperature gases for use in control actuation, air bottle pressurization, etc. Packed beds of lengths 200 and 300 mm were characterized for packing-sphere-based Reynolds numbers ranging from 0.8 x 10(4) to 8.5 x 10(4).The solid packing used in the bed consisted of phi 9.5 mm steel spheres. The bed-to-particle diameter ratio was with the average packed-bed porosity around 0.43. The inlet flow temperature was unsteady and a mesh of spheres was used at either end to eliminate flow entrance and exit effects. Gas temperature and pressure were measured at the entry, exit,and at three axial locations along centerline in the packed beds. The solid packing temperature was measured at three axial locations in the packed bed. A correlation based on the ratio of pressure drop and inlet-flow momentum (Euler number) exhibited an asymptotically decreasing trend with increasing Reynolds number. Axial conduction across the packed bed was found to he negligible in the investigated Reynolds number range. The enthalpy absorption rate to solid packing from hot gases is plotted as a function of a nondimensional time constant for different Reynolds numbers. A longer packed bed had high enthalpy absorption rate at Reynolds number similar to 10(4), which decreased at Reynolds number similar to 10(5). The enthalpy absorption plots can be used for estimating enthalpy drop across packed bed with different material, but for a geometrically similar packing.

Algebraic and geometric intersection numbers for free groups

We show that the algebraic intersection number of Scott and Swarup for splittings of free groups Coincides With the geometric intersection number for the sphere complex of the connected sum of copies of S-2 x S-1. (C) 2009 Elsevier B.V. All rights reserved.

Generalized Burgers equations and Euler-Painleve transcendents. II

It was proposed earlier [P. L. Sachdev, K. R. C. Nair, and V. G. Tikekar, J. Math. Phys. 27, 1506 (1986)] that the Euler Painlevé equation yy[script `]+ay[script ']2+ f(x)yy[script ']+g(x) y2+by[script ']+c=0 represents the generalized Burgers equations (GBE's) in the same manner as Painlevé equations do the KdV type. The GBE was treated with a damping term in some detail. In this paper another GBE ut+uaux+Ju/2t =(gd/2)uxx (the nonplanar Burgers equation) is considered. It is found that its self-similar form is again governed by the Euler Painlevé equation. The ranges of the parameter alpha for which solutions of the connection problem to the self-similar equation exist are obtained numerically and confirmed via some integral relations derived from the ODE's. Special exact analytic solutions for the nonplanar Burgers equation are also obtained. These generalize the well-known single hump solutions for the Burgers equation to other geometries J=1,2; the nonlinear convection term, however, is not quadratic in these cases. This study fortifies the conjecture regarding the importance of the Euler Painlevé equation with respect to GBE's. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Studies on friction and formation of transfer layer when Al-4Mg alloy pins slid at various numbers of cycles on steel plates of different surface texture

In the present investigation, various kinds of textures, namely, unidirectional, 8-ground, and random were attained on the die surfaces. Roughness of the textures was varied using different grits of emery papers or polishing powders. Then pins made of Al-4Mg alloys were slid against steel plates at various numbers of cycles, namely 1, 2, 6, 10 and 20 under both dry and lubricated conditions using an inclined pin-on-plate sliding tester. The morphologies of the worn surfaces of the pins and the formation of transfer layer on the counter surfaces were observed using a scanning electron microscope. Surface roughness parameters of the plate were measured using an optical profilometer. It was observed that the coefficient of friction and formation of transfer layer during the first few cycles depend on the die surface textures under both dry and lubricated conditions. It was also observed that under lubricated condition, the coefficient of friction decreases with number of cycles for all kinds of textures. However, under dry condition, it ecreases for unidirectional and 8-ground surfaces while for random surfaces it increases with number of cycles

Generalized Burgers equations and Euler-Painlev transcendents. I

Initial-value problems for the generalized Burgers equation (GBE) ut+u betaux+lambdaualpha =(delta/2)uxx are discussed for the single hump type of initial data both continuous and discontinuous. The numerical solution is carried to the self-similar ``intermediate asymptotic'' regime when the solution is given analytically by the self-similar form. The nonlinear (transformed) ordinary differential equations (ODE's) describing the self-similar form are generalizations of a class discussed by Euler and Painlevé and quoted by Kamke. These ODE's are new, and it is postulated that they characterize GBE's in the same manner as the Painlev equations categorize the Kortweg-de Vries (KdV) type. A connection problem for some related ODE's satisfying proper asymptotic conditions at x=±[infinity], is solved. The range of amplitude parameter is found for which the solution of the connection problem exists. The other solutions of the above GBE, which display several interesting features such as peaking, breaking, and a long shelf on the left for negative values of the damping coefficient lambda, are also discussed. The results are compared with those holding for the modified KdV equation with damping. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Hydration numbers of some metal acetates, monochloroacetates and trichloroacetates in solution from ultrasonic velocity and conpressibility measurements

Ultrasonic velocities in aqueous solutions of some metal acetates, monochloroacelates and trichloroacetates, and the respective acids have been measured at 1 MHz frequency using the pulse technique. The ultrsonic velocity, adiabatic compressibility and apperent molal compressibility were measured as a function of concentration. The apparent molal compressibility values at infinite dilution were calculated and used to determine the hydration numbers.

Mass transfer with surface reaction in arrays of cylinders at low reynolds and high peclet numbers

Creeping flow hydrodynamics combined with diffusion boundary layer equation are solved in conjunction with free-surface cell model to obtain a solution of the problem of convective transfer with surface reaction for flow parallel to an array of cylindrical pellets at high Peclet numbers and under fast and intermediate kinetics regimes. Expressions are derived for surface concentration, boundary layer thickness, mass flux and Sherwood number in terms of Damkoehler number, Peclet number and void fraction of the array. The theoretical results are evaluated numerically.

Two-way converter logic for positive and negative binary numbers

This paper describes a hardware implementation of a two-way converter logic by which conversion between numbers from positive to negative binary representation is possible. Index terms: (i) Negative radix, (ii) Positive radix, (iii) Two-way conversion.

The Formation, Numbers And Radio Output Of Giant Radio Galaxies

The numbers and mean radio luminosities of giant radio galaxies (GRGs) have been calculated for redshifts up to z = 0.6, assuming a sensitivity limit of 1 Jy at 1 GHz for the observations. The estimates are obtained with a model for the beam propagation, first through the hot gaseaous halo around the parent galaxy, and thereafter, through the even hotter but less dense intergalactic medium. The model is able to accurately reproduce the observed numbers and mean radio luminosities of GRGs at redshifts of less than 0.1, and it predicts that a somewhat larger number of GRGs should be found at redshifts of greater than 0.1.

PVU and wave-particle splitting schemes for Euler equations of gas dynamics

A new way of flux-splitting, termed as the wave-particle splitting is presented for developing upwind methods for solving Euler equations of gas dynamics. Based on this splitting, two new upwind methods termed as Acoustic Flux Vector Splitting (AFVS) and Acoustic Flux Difference Splitting (AFDS) methods are developed. A new Boltzmann scheme, which closely resembles the wave-particle splitting, is developed using the kinetic theory of gases. This method, termed as Peculiar Velocity based Upwind (PVU) method, uses the concept of peculiar velocity for upwinding. A special feature of all these methods that the unidirectional and multidirectional parts of the flux vector are treated separately. Extensive computations done using these schemes demonstrate the soundness of the ideas.