Study of third-order non-linear systems- Variation of parameters approach

In this paper, the transient response of a third-order non-linear system is obtained by first reducing the given third-order equation to three first-order equations by applying the method of variation of parameters. On the assumption that the variations of amplitude and phase are small, the functions are expanded in ultraspherical polynomials. The expansion is restricted to the constant term. The resulting equations are solved to obtain the response of the given third-order system. A numerical example is considered to illustrate the method. The results show that the agreement between the approximate and digital solution is good thus vindicating the approximation.

A critical analysis of the shear difference method

The shear difference method which is commonly used for the separation of normal stresses using photoelastic techniques depends on the step-by-step integration of one of the differential equations of equilibrium. It is assumed that the isoclinic and the isochromatic parameters measured by the conventional methods pertain to the state of stress at the midpoint of the light path. In practice, a slice thin enough for the above assumption to be true and at the same time thick enough to give differences in the shear-stress values over the thickness is necessary. The paper discusses the errors introduced in the isoclinic and isochromatic values by the conventional methods neglecting the variation of stresses along the light path. It is shown that while the error introduced in the measurement of the isochromatic parameter may not be serious, the error caused in the isoclinic measurement may lead to serious errors. Since the shear-difference method involves step-by-step integration the error introduced will be of a cumulative nature.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

The line transect method for estimating densities of large mammals in a tropical deciduous forest: An evaluation of models and field experiments

We have evaluated techniques of estimating animal density through direct counts using line transects during 1988-92 in the tropical deciduous forests of Mudumalai Sanctuary in southern India for four species of large herbivorous mammals, namely, chital (Axis axis), sambar (Cervus unicolor), Asian elephant (Elephas maximus) and gaur (Bos gauras). Density estimates derived from the Fourier Series and the Half-Normal models consistently had the lowest coefficient of variation. These two models also generated similar mean density estimates. For the Fourier Series estimator, appropriate cut-off widths for analysing line transect data for the four species are suggested. Grouping data into various distance classes did not produce any appreciable differences in estimates of mean density or their variances, although model fit is generally better when data are placed in fewer groups. The sampling effort needed to achieve a desired precision (coefficient of variation) in the density estimate is derived. A sampling effort of 800 km of transects returned a 10% coefficient of variation on estimate for chital; for the other species a higher effort was needed to achieve this level of precision. There was no statistically significant relationship between detectability of a group and the size of the group for any species. Density estimates along roads were generally significantly different from those in the interior af the forest, indicating that road-side counts may not be appropriate for most species.

Crystallization kinetics of SrBi2B2O7 glasses by non-isothermal methods

Transparent glasses of SrBi2B2O7 (SBBO) were fabricated via the conventional melt-quenching technique. The amorphous and the glassy nature of the as-quenched samples were, respectively, confirmed by X-ray powder diffraction (XRD) and differential scanning calorimetry (DSC). The glass transition (T (g)) and the crystallization parameters [crystallization activation energy (E (cr)) and Avrami exponent (n)] were evaluated under non-isothermal conditions using DSC. There was a close agreement between the activation energies for the crystallization process determined by Augis and Bennet and Kissinger methods. The variation of local activation energy [E (c)(x)] that was determined by Ozawa method, decreased with the fraction of crystallization (x). The Avrami exponent (n(x)) increased with the increase in fraction of crystallization (x) suggesting that there was a change over in the crystallization process from the surface to the bulk.

Beyond fractional derivatives: local approximation of other convolution integrals

Dynamic systems involving convolution integrals with decaying kernels, of which fractionally damped systems form a special case, are non-local in time and hence infinite dimensional. Straightforward numerical solution of such systems up to time t needs O(t(2)) computations owing to the repeated evaluation of integrals over intervals that grow like t. Finite-dimensional and local approximations are thus desirable. We present here an approximation method which first rewrites the evolution equation as a coupled in finite-dimensional system with no convolution, and then uses Galerkin approximation with finite elements to obtain linear, finite-dimensional, constant coefficient approximations for the convolution. This paper is a broad generalization, based on a new insight, of our prior work with fractional order derivatives (Singh & Chatterjee 2006 Nonlinear Dyn. 45, 183-206). In particular, the decaying kernels we can address are now generalized to the Laplace transforms of known functions; of these, the power law kernel of fractional order differentiation is a special case. The approximation can be refined easily. The local nature of the approximation allows numerical solution up to time t with O(t) computations. Examples with several different kernels show excellent performance. A key feature of our approach is that the dynamic system in which the convolution integral appears is itself approximated using another system, as distinct from numerically approximating just the solution for the given initial values; this allows non-standard uses of the approximation, e. g. in stability analyses.

A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

Moment methods and "restricted variation" in kinetic theory

It is shown that there is a strict one-to-one correspondence between results obtained by the use of "restricted" variational principles and those obtained by a moment method of the Mott-Smith type for shock structure.

Effect of calcium deficiency on the mechanical properties of hydroxyapatite crystals

The deterioration of the mechanical properties of bone with age is related to several factors including the structure, organization and chemistry of the constituent phases; however, the relative contribution of each of these factors is not well understood. In this study, we have investigated the effect of chemistry (calcium deficiency) on the mechanical properties of single crystals of hydroxyapatite. Single crystals of stoichiometric crystals grown by the flux method and calcium-deficient platelet crystals grown using wet chemical methods were used as model systems. Using nanoindentation, we show that calcium deficiency leads to an 80% reduction in the hardness and elastic modulus and at least a 75% reduction in toughness in plate-shaped hydroxyapatite crystals. Measurement of local mechanical properties using nanoindentation and nanoscale chemistry through elemental mapping in a transmission electron microscope points to a direct correlation between the observed spatial variation in composition and the large scatter in the measured hardness and modulus values. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

On Application Of The Method Of Kinetic Invariants To The Description Of Dissolution Accompanied By A Chemical Reaction

The dissolution, accompanied by chemical reaction, of monodisperse solid particles has been analysed. The resulting model, which accounts for the variation of mass transfer coefficient with the size of the dissolving particles, yields an approximate analytical form of a kinetic function. Rigorous numerical and approximate analytical solutions have been obtained for the governing system of nonlinear ordinary differential equations. The transient nature of the dissolution process as well as the accuracy of the analytical solution is brought out by the rigorous numerical solution. The analytical solution is fairly accurate for the major part of the range of operational times encountered in practice.

Evaluation of ductile-brittle transition temperature (DBTT) of aluminide bond coats by micro-tensile test method

A method has been suggested to accurately determine the DBTT of diffusion aluminide bond coats. Micro-tensile testing of free-standing coating samples has been carried out. The DBTT was determined based on the variation of plastic strain-to-fracture with temperature. The positive features of this method over the previously reported techniques are highlighted. (C) 2010 Elsevier B.V. All rights reserved.

Local variability and base sequence effects in DNA crystal structures.

The importance and usefulness of local doublet parameters in understanding sequence dependent effects has been described for A- and B-DNA oligonucleotide crystal structures. Each of the two sets of local parameters described by us in the NUPARM algorithm, namely the local doublet parameters, calculated with reference to the mean z-axis, and the local helical parameters, calculated with reference to the local helix axis, is sufficient to describe the oligonucleotide structures, with the local helical parameters giving a slightly magnified picture of the variations in the structures. The values of local doublet parameters calculated by NUPARM algorithm are similar to those calculated by NEWHELIX90 program, only if the oligonucleotide fragment is not too distorted. The mean values obtained using all the available data for B-DNA crystals are not significantly different from those obtained when a limited data set is used, consisting only of structures with a data resolution of better than 2.4 A and without any bound drug molecule. Thus the variation observed in the oligonucleotide crystals appears to be independent of the quality of their crystallinity. No strong correlation is seen between any pair of local doublet parameters but the local helical parameters are interrelated by geometric relationships. An interesting feature that emerges from this analysis is that the local rise along the z-axis is highly correlated with the difference in the buckle values of the two basepairs in the doublet, as suggested earlier for the dodecamer structures (Bansal and Bhattacharyya, in Structure & Methods: DNA & RNA, Vol. 3 (Eds., R.H. Sarma and M.H. Sarma), pp. 139-153 (1990)). In fact the local rise values become almost constant for both A- and B-forms, if a correction is applied for the buckling of the basepairs. In B-DNA the AA, AT, TA and GA basepair sequences generally have a smaller local rise (3.25 A) compared to the other sequences (3.4 A) and this seems to be an intrinsic feature of basepair stacking interaction and not related to any other local doublet parameter. The roll angles in B-DNA oligonucleotides have small values (less than +/- 8 degrees), while mean local twist varies from 24 degrees to 45 degrees. The CA/TG doublet sequences show two types of preferred geometries, one with positive roll, small positive slide and reduced twist and another with negative roll, large positive slide and increased twist.(ABSTRACT TRUNCATED AT 400 WORDS)

A discrete state perspective of manipulator workspaces Vorstellungen von diskreten zuständen von manipulator-arbeitsräumen

The finite resolution of joint drives or sensors imparts a discrete nature to the joints of a manipulator. Because of this an arbitrary point in the workspace cannot be reached without error even in ideal mechanical environment. This paper investigates the effect of this discrete nature of the joints on the accuracy of performance of a manipulator and develops a method to select the joint states to reach a point with least error. It is shown that the configuration leading to least error cannot, in general, be found from configuration space, especially when there is large variation in the link lengths or joint resolutions or both. The anomaly becomes severe when the gross motion of the end-effector approaches the local resolution of the workspace. The paper also shows how to distinguish two workspaces which may be identical so far as the boundary points are concerned, taking the joint resolutions into account. Finally, the concepts have been extended to define continuous space global and local performance indices for general multi degree of freedom manipulators.

A novel method for diameter estimation of small opaque objects using Fraunhofer diffraction

A novel optical method is proposed and demonstrated, for real-time dimension estimation of thin opaque cylindrical objects. The methodology relies on free-space Fraunhofer diffraction principle. The central region, of such tailored diffraction pattern obtained under suitable choice of illumination conditions, comprises of a pair of `equal intensity maxima', whose separation remains constant and independent of the diameter of the diffracting object. An analysis of `the intensity distribution in this region' reveals the following. At a point symmetrically located between the said maxima, the light intensity varies characteristically with diameter of the diffracting object, exhibiting a relatively stronger intensity modulation under spherical wave illumination than under a plane wave illumination. The analysis reveals further, that the said intensity variation with diameter is controllable by the illumination conditions. Exploiting these `hitherto unexplored' features, the present communication reports for the first time, a reliable method of estimating diameter of thin opaque cylindrical objects in real-time, with nanometer resolution from single point intensity measurement. Based on the proposed methodology, results of few simulation and experimental investigations carried-out on metallic wires with diameters spanning the range of 5 to 50 mu m, are presented. The results show that proposed method is well-suited for high resolution on-line monitoring of ultrathin wire diameters, extensively used in micro-mechanics and semiconductor industries, where the conventional diffraction-based methods fail to produce accurate results.

Anomalous X-ray scattering study of local structures in the superionic conducting class (CuI)(0.3)(Cu2O)(0.35)(MoO3)(0.35)

The anomalous X-ray scattering (AXS) method using Cu and Mo K absorption edges has been employed for obtaining the local structural information of superionic conducting glass having the composition (CuI)(0.3)(Cu2O)(0.35)(MoO3)(0.35). The possible atomic arrangements in near-neighbor region of this glass were estimated by coupling the results with the least-squares analysis so as to reproduce two differential intensity profiles for Cu and Mo as well as the ordinary scattering profile. The coordination number of oxygen around Mo is found to be 6.1 at the distance of 0.187 nm. This implies that the MoO6 octahedral unit is a more probable structural entity in the glass rather than MoO4 tetrahedra which has been proposed based on infrared spectroscopy. The pre-peak shoulder observed at about 10 nm(-1) may be attributed to density fluctuation originating from the MoO6 octahedral units connected with the corner sharing linkage, in which the correlation length is about 0.8 nm. The value of the coordination number of I- around Cu+ is estimated as 4.3 at 0.261 nm, suggesting an arrangement similar to that in molten CuI.