Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.

Stimulated interdiffusion and optical recording in Sb/As2S3 nanomultilayers

The role of the compositional modulation at nano-scale dimensions (similar to 2-10 nm) in the enhancement of optical recording parameters in nanomultilayers, which contain Sb as active, optical absorbing and diffusing layers and As2S3 as barrier (matrix) layers was investigated. Comparison was made with single homogeneous layers made of ternary (As2S3)(x)Sb1-x glasses and co-deposited from Sb and As2S3. It was shown that essential increase of the recording efficiency, sensitivity of the bleaching process, broadening of its spectral range occurs due to the stimulated interdiffusion of adjacent components in Sb/As2S3 nanomultilayers with optimized Sb layer thickness.

Sound Velocity Anomaly at the Mott Transition: Application to Organic Conductors and V2O3

Close to the Mott transition, lattice degrees of freedom react to the softening of electron degrees of freedom. This results in a change of lattice spacing, a diverging compressibility, and a critical anomaly of the sound velocity. These effects are investigated within a simple model, in the framework of dynamical mean-field theory. The results compare favorably to recent experiments on the layered organic-conductor kappa-(BEDT-TTF)(2)Cu[N(CN)(2)]Cl. We predict that effects of a similar magnitude are expected for V2O3, despite the much larger value of the elastic modulus of this material.

Crumpling and second sound in lyotropic lamellar phases

We study relative concentration fluctuations in two component lamellar smectic liquid crystals consisting of surfactant layers of width w separated by a background fluid and show that these fluctuations are dominated by crumpling fluctuations of the surfactant layers when (w/l)2

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Continuous Recording of Retardation and Intensity of Echoes from the Ionosphere

Pulse retardation method of Breit and Tuve has been modified to record continuously the equivalent height as well as the intensity of reflections from the ionosphere. Synchronized pulses are transmitted, and the received ground pulse and the reflected pulses, after amplification and suitable distortion, are applied to the focusing cylinder of a cathode ray tube the horizontal deflecting plates of which are connected to a synchronized linear time base circuit. The pattern on the screen is composed of a bright straight line corresponding to the time base with dark gaps corresponding to the received pulses. The distance between the initial points of the gaps represents retardation while the widths of the gaps correspond to the intensity of the pulses. The pattern is photographed on a vertically moving film. One of the first few records taken at Bangalore on 4 megacycles is reproduced. It shows, among other things, that the less retarded component of magneto-ionic splitting from the F layer is present most of the time. Whenever the longer retardation component does occur, it has stronger intensity than the former. Towards the late evening hours, just before disappearing, when the F layer rises and exhibits magnetoionic splitting, the intensity of the less retarded component is extremely low compared with the other component.

Changing resonator geometry to boost sound power decouples size and song frequency in a small insect

Despite their small size, some insects, such as crickets, can produce high amplitude mating songs by rubbing their wings together. By exploiting structural resonance for sound radiation, crickets broadcast species-specific songs at a sharply tuned frequency. Such songs enhance the range of signal transmission, contain information about the signaler's quality, and allow mate choice. The production of pure tones requires elaborate structural mechanisms that control and sustain resonance at the species-specific frequency. Tree crickets differ sharply from this scheme. Although they use a resonant system to produce sound, tree crickets can produce high amplitude songs at different frequencies, varying by as much as an octave. Based on an investigation of the driving mechanism and the resonant system, using laser Doppler vibrometry and finite element modeling, we show that it is the distinctive geometry of the crickets' forewings (the resonant system) that is responsible for their capacity to vary frequency. The long, enlarged wings enable the production of high amplitude songs; however, as a mechanical consequence of the high aspect ratio, the resonant structures have multiple resonant modes that are similar in frequency. The drive produced by the singing apparatus cannot, therefore, be locked to a single frequency, and different resonant modes can easily be engaged, allowing individual males to vary the carrier frequency of their songs. Such flexibility in sound production, decoupling body size and song frequency, has important implications for conventional views of mate choice, and offers inspiration for the design of miniature, multifrequency, resonant acoustic radiators.

Monaural sound source separation using covariance profile of partials

This paper addresses the problem of separation of pitched sounds in monaural recordings. We present a novel feature for the estimation of parameters of overlapping harmonics which considers the covariance of partials of pitched sounds. Sound templates are formed from the monophonic parts of the mixture recording. A match for every note is found among these templates on the basis of covariance profile of their harmonics. The matching template for the note provides the second order characteristics for the overlapped harmonics of the note. The algorithm is tested on the RWC music database instrument sounds. The results clearly show that the covariance characteristics can be used to reconstruct overlapping harmonics effectively.

An 8-to-1 bit 1-MS/s SAR ADC With VGA and Integrated Data Compression for Neural Recording

Low power consumption per channel and data rate minimization are two key challenges which need to be addressed in future generations of neural recording systems (NRS). Power consumption can be reduced by avoiding unnecessary processing whereas data rate is greatly decreased by sending spike time-stamps along with spike features as opposed to raw digitized data. Dynamic range in NRS can vary with time due to change in electrode-neuron distance or background noise, which demands adaptability. An analog-to-digital converter (ADC) is one of the most important blocks in a NRS. This paper presents an 8-bit SAR ADC in 0.13-mu m CMOS technology along with input and reference buffer. A novel energy efficient digital-to-analog converter switching scheme is proposed, which consumes 37% less energy than the present state-of-the-art. The use of a ping-pong input sampling scheme is emphasized for multichannel input to alleviate the bandwidth requirement of the input buffer. To reduce the data rate, the A/D process is only enabled through the in-built background noise rejection logic to ensure that the noise is not processed. The ADC resolution can be adjusted from 8 to 1 bit in 1-bit step based on the input dynamic range. The ADC consumes 8.8 mu W from 1 V supply at 1 MS/s speed. It achieves effective number of bits of 7.7 bits and FoM of 42.3 fJ/conversion-step.

A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

Timing Recovery Algorithms and Architectures for 2-D Magnetic Recording Systems

We investigate the problem of timing recovery for 2-D magnetic recording (TDMR) channels. We develop a timing error model for TDMR channel considering the phase and frequency offsets with noise. We propose a 2-D data-aided phase-locked loop (PLL) architecture for tracking variations in the position and movement of the read head in the down-track and cross-track directions and analyze the convergence of the algorithm under non-separable timing errors. We further develop a 2-D interpolation-based timing recovery scheme that works in conjunction with the 2-D PLL. We quantify the efficiency of our proposed algorithms by simulations over a 2-D magnetic recording channel with timing errors.

A model supersonic buried-nozzle jet: instability and acoustic wave scattering and the far-field sound

We consider sound source mechanisms involving the acoustic and instability modes of dual-stream isothermal supersonic jets with the inner nozzle buried within an outer shroud-like nozzle. A particular focus is scattering into radiating sound waves at the shroud lip. For such jets, several families of acoustically coupled instability waves exist, beyond the regular vortical Kelvin-Helmholtz mode, with different shapes and propagation characteristics, which can therefore affect the character of the radiated sound. In our model, the coaxial shear layers are vortex sheets while the incident acoustic disturbances are the propagating shroud modes. The Wiener-Hopf method is used to compute their scattering at the sharp shroud edge to obtain the far-field radiation. The resulting far-field directivity quantifies the acoustic efficiency of different mechanisms, which is particularly important in the upstream direction, where the results show that the scattered sound is more intense than that radiated directly by the shear-layer modes.