Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Theory of laminar flame propagation with non-normal diffusion.

A study has been made of the problem of steady, one-dimensional, laminar flame propagation in premixed gases, with the Lewis number differing from (and equal to) unity. Analytical solutions, using the method of matched asymptotic expansions, have been obtained for large activation energies. Numerical solutions have been obtained for a wide range of the reduced activation temperature parameter (n {geometrically equal to} E/RTb), and the Lewis number δ. The studies reveal that the flame speed eigenvalue is linear in Lewis number for first order and quadratic in Lewis number for second order reactions. For a quick determination of flame speeds, with reasonable accuracy, a simple rule, expressing the flame speed eigenvalue as a function of the Lewis number and the centroid of the reaction rate function, is proposed. Comparisons have been made with some of the earlier works, for both first and second order reactions.

On the study of a third-order mechanical oscillator

In this paper, the study of a third-order mechanical oscillator is presented by demonstrating its equivalence to the well-known R.C. multivibrator with two additional reactive elements. The conditions for the oscillator's possession of periodic solutions are presented. It is also shown that under certain conditions, the study of the given third-order autonomous system can be reduced to the study of an equivalent second-order, non-autonomous system.

Non-Linear vibration of an elastic string

Equations proposed in previous work on the non-linear motion of a string show a basic disagreement, which is here traced to an assumption about the longitudinal displacement u. It is shown that it is neither necessary nor justifiable to assume that u is zero; and also that the velocity of propagation of u disturbances in a string is different from that in an infinite medium, although this difference is usually negligible. After formulating the exact equations of motion for the string, a systematic procedure is described for obtaining approximations to these equations to any order, making only the assumption that the strain in the material of the string is small. The lowest order equations in this scheme are non-linear, and are used to describe the response of a string near resonance. Finally, it is shown that in the absence of damping, planar motion of a string is always unstable at sufficiently high amplitudes, the critical amplitude falling to zero at the natural frequency and its subharmonics. The effect of slight damping on this instability is also discussed.

Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.

Simple kinetic sensor to structural transitions

Driven nonequilibrium structural phase transformation has been probed using time-varying resistance fluctuations or noise. We demonstrate that the non-Gaussian component (NGC) of noise obtained by evaluating the higher-order statistics of fluctuations, serves as a simple kinetic detector of these phase transitions. Using the Martensite transformation in free-standing wires of nickel-titanium binary alloys as a prototype, we observe clear deviations from the Gaussian background in the transformation zone, indicative of the long-range correlations in the system as the phase transforms. The viability of non-Gaussian statistics as a robust probe to structural phase transition was also confirmed by comparing the results from differential scanning calorimetry measurements. We further studied the response of the NGC to the modifications in the microstructure on repeated thermal cycling, as well as the variations in the temperature-drive rate, and explained the results using established simplistic models based on the different competing time scales. Our experiments (i) suggest an alternative method to estimate the transformation temperature scales with high accuracy and (ii) establish a connection between the material-specific evolution of microstructure to the statistics of its linear response. Since the method depends on an in-built long-range correlation during transformation, it could be portable to other structural transitions, as well as to materials of different physical origin and size.

Design of first- and second-order sliding mode observers for induction motors using a stator-flux model

We extend current research in the area of 'sensorless' control of induction motors by presenting two observers based on first- and second-order sliding mode control theories for the simultaneous estimation of flux and speed. We base the observers on the stator-flux model of the motor instead of the usual rotor-flux model mainly because of the uncertain rotor resistance that plays a significant role in the latter. By designing the observers as if they are sliding mode controllers, we lend the properties of parameter insensitive closed-loop dynamics and finite time convergence to the stator flux and speed estimation schemes. We also present simulation and experimental results to validate the operation of the observers.

Charge-order driven multiferroic and magneto-dielectric properties of rare earth manganates

Charge-order driven magnetic ferroelectricity is shown to occur in several rare earth manganates of the general formula, Ln(1-x)A(x)MnO(3) (In = rare earth, A = alkaline earth). Charge-ordered manganates exhibit dielectric constant anomalies around the charge-ordering or the antiferromagnetic transition temperature. Magnetic fields have a marked effect on the dielectric properties of these compounds, indicating the presence of coupling between the magnetic and electrical order parameters. Magneto-dielectric properties are retained in small particles of the manganates. The observation of magneto-ferroelectricity in these manganates is in accordance with theoretical predictions.

Charge-order driven multiferroic and magneto-dielectric properties of rare earth manganates

Charge-order driven magnetic ferroelectricity is shown to occur in several rare earth manganates of the general formula, Ln(1-x)A(x)MnO(3) (In = rare earth, A = alkaline earth). Charge-ordered manganates exhibit dielectric constant anomalies around the charge-ordering or the antiferromagnetic transition temperature. Magnetic fields have a marked effect on the dielectric properties of these compounds, indicating the presence of coupling between the magnetic and electrical order parameters. Magneto-dielectric properties are retained in small particles of the manganates. The observation of magneto-ferroelectricity in these manganates is in accordance with theoretical predictions.

Systematic bias in MODIS dust aerosol retrieval at Kanpur (AERONET), Indo-Gangetic Basin - art.no.640816

We have compared the spectral aerosol optical depth (AOD) and aerosol fine mode fraction (AFMF) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) with those of Aerosol Robotic Network (AERONET) at Kanpur (26.45N, 80.35E), northern India for the pre-monsoon season (March to June, 2001-2005). We found that MODIS systematically overestimates AOD during pre-monsoon season (known to be influenced by dust transport from north-west of India). The errors in AOD were correlated with the MODIS top-of-atmosphere apparent surface reflectance in 2.1 mu m channel (rho*(2.1)). MODIS aerosol algorithm uses p*(2.1) to derive the surface reflectance in visible channels (rho(0.47), rho(0.66)) using an empirical mid IR-visible correlation (rho(0.47) = rho(2.1)/4, rho(0.66) = rho(2.1)/2). The large uncertainty in estimating surface reflectance in visible channels (Delta rho(0.66)+/- 0.04, Delta rho(0.47)+/- 0.02) at higher values of p*(2.1) (p*(2.1) > 0.18) leads to higher aerosol contribution in the total reflected radiance at top-of atmosphere to compensate for the reduced surface reflectance in visible channels and thus leads to overestimation of AOD. This was also reflected in the very low values of AFMF during pre-monsoon whose accuracy depends on the aerosol path radiance in 0.47 and 0.66 mu m channels and aerosol models. The errors in AOD were also high in the scattering angle range 110 degrees-140 degrees, where the effect of dust non-spherity on its optical properties is significant. The direct measurements of spectral surface reflectance are required over the Indo-Gangetic basin in order to validate the mid IR-visible relationship. MODIS aerosol models should also be modified to incorporate the effect of non-spherity of dust aerosols.

Robust fixed order lateral H-2 controller for micro air vehicle

This paper presents a robust fixed order H-2 controller design using Strengthened discrete optimal projection equations, which approximate the first order necessary optimality condition. The novelty of this work is the application of the robust H-2 controller to a micro aerial vehicle named Sarika2 developed in house. The controller is designed in discrete domain for the lateral dynamics of Sarika2 in the presence of low frequency atmospheric turbulence (gust) and high frequency sensor noise. The design specification includes simultaneous stabilization, disturbance rejection and noise attenuation over the entire flight envelope of the vehicle. The resulting controller performance is comprehensively analyzed by means of simulation.

Extending lagrange interpolation to develop a 4-node quadrilateral element

Today finite element method is a well established tool in engineering analysis and design. Though there axe many two and three dimensional finite elements available, it is rare that a single element performs satisfactorily in majority of practical problems. The present work deals with the development of 4-node quadrilateral element using extended Lagrange interpolation functions. The classical univariate Lagrange interpolation is well developed for 1-D and is used for obtaining shape functions. We propose a new approach to extend the Lagrange interpolation to several variables. When variables axe more than one the method also gives the set of feasible bubble functions. We use the two to generate shape function for the 4-node arbitrary quadrilateral. It will require the incorporation of the condition of rigid body motion, constant strain and Navier equation by imposing necessary constraints. The procedure obviates the need for isoparametric transformation since interpolation functions are generated for arbitrary quadrilateral shapes. While generating the element stiffness matrix, integration can be carried out to the accuracy desired by dividing the quadrilateral into triangles. To validate the performance of the element which we call EXLQUAD4, we conduct several pathological tests available in the literature. EXLQUAD4 predicts both stresses and displacements accurately at every point in the element in all the constant stress fields. In tests involving higher order stress fields the element is assured to converge in the limit of discretisation. A method thus becomes available to generate shape functions directly for arbitrary quadrilateral. The method is applicable also for hexahedra. The approach should find use for development of finite elements for use with other field equations also.

Second order cone programming approaches for handling missing and uncertain data

We propose a novel second order cone programming formulation for designing robust classifiers which can handle uncertainty in observations. Similar formulations are also derived for designing regression functions which are robust to uncertainties in the regression setting. The proposed formulations are independent of the underlying distribution, requiring only the existence of second order moments. These formulations are then specialized to the case of missing values in observations for both classification and regression problems. Experiments show that the proposed formulations outperform imputation.