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.

Solution of MDPS using simulation-based value iteration

This article proposes a three-timescale simulation based algorithm for solution of infinite horizon Markov Decision Processes (MDPs). We assume a finite state space and discounted cost criterion and adopt the value iteration approach. An approximation of the Dynamic Programming operator T is applied to the value function iterates. This 'approximate' operator is implemented using three timescales, the slowest of which updates the value function iterates. On the middle timescale we perform a gradient search over the feasible action set of each state using Simultaneous Perturbation Stochastic Approximation (SPSA) gradient estimates, thus finding the minimizing action in T. On the fastest timescale, the 'critic' estimates, over which the gradient search is performed, are obtained. A sketch of convergence explaining the dynamics of the algorithm using associated ODEs is also presented. Numerical experiments on rate based flow control on a bottleneck node using a continuous-time queueing model are performed using the proposed algorithm. The results obtained are verified against classical value iteration where the feasible set is suitably discretized. Over such a discretized setting, a variant of the algorithm of [12] is compared and the proposed algorithm is found to converge faster.

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.

Effect of boundary condition description on convergence of solution in a boundary-value problem

It is well known that the numerical accuracy of a series solution to a boundary-value problem by the direct method depends on the technique of approximate satisfaction of the boundary conditions and on the stage of truncation of the series. On the other hand, it does not appear to be generally recognized that, when the boundary conditions can be described in alternative equivalent forms, the convergence of the solution is significantly affected by the actual form in which they are stated. The importance of the last aspect is studied for three different techniques of computing the deflections of simply supported regular polygonal plates under uniform pressure. It is also shown that it is sometimes possible to modify the technique of analysis to make the accuracy independent of the description of the boundary conditions.

Double time window targeting technique: Real-time DMRG dynamics in the Pariser-Parr-Pople model

We present a generalized adaptive time-dependent density matrix renormalization-group (DMRG) scheme, called the double time window targeting (DTWT) technique, which gives accurate results with nominal computational resources, within reasonable computational time. This procedure originates from the amalgamation of the features of pace keeping DMRG algorithm, first proposed by Luo et al. [Phys. Rev. Lett. 91, 049701 (2003)] and the time-step targeting algorithm by Feiguin and White [Phys. Rev. B 72, 020404 (2005)]. Using the DTWT technique, we study the phenomena of spin-charge separation in conjugated polymers (materials for molecular electronics an spintronics), which have long-range electron-electron interactions and belong to the class of strongly correlated low-dimensional many-body systems. The issue of real-time dynamics within the Pariser-Parr-Pople (PPP) model which includes long-range electron correlations has not been addressed in the literature so far. The present study on PPP chains has revealed that, (i) long-range electron correlations enable both the charge and spin degree of freedom of the electron, to propagate faster in the PPP model compared to Hubbard model, (ii) for standard parameters of the PPP model as applied to conjugated polymers, the charge velocity is almost twice that of the spin velocity, and (iii) the simplistic interpretation of long-range correlations by merely renormalizing the U value of the Hubbard model fails to explain the dynamics of doped holes/electrons in the PPP model.

Role of aluminum-magnesium alloys in humid aging of composite solid propellants

The humid aging of composite propellants containing a terpolymer of polybutadiene, acrylic acid, and acrylonitrile (PBAN) as a binder has been studied as a function of aging temperature, relative humidity, and aging time. Three composite types - AP-PBAN, AP-Al-PBAN, and AP-(Al-Mg) alloy- PBAN - have been studied. The burning rates of all three propellant types were unaffected by aging. The calorimetric values of composites containing aluminum-magnesium alloy decreased on aging, and the lattice parameter of the alloy decreased to a value close to that of aluminum. Water absorption in all of the samples increased with increases in the temperature, relative humidity, and aging time. The compression strength of the nonmetalized and aluminized samples decreased on aging, whereas that of the composites containing the alloy increased. The latter effect has been traced to reaction of residual carboxyl groups on the polymer chains with magnesium, leading to cross-linking. The reaction between the -COOH groups and magnesium has been proved using infrared spectroscopy. (Author)

Crystal structure of L-lysyl-L-glutamic acid dihydrate

L-Lysyl-L-glutamic acid dihydrate, C11N3O5H21·2H2O, crystallizes in the monoclinic space group P21 with a = 12.474(2), b = 5.020(1), c = 13.157(2) Å, β= 114.69(1)° and Z = 2. The crystal structure was solved by direct methods and refined to an R value of 0.037 using full matrix least-squares method. The molecule exists as a double zwitterion with both the amino and carboxyl groups ionised. The peptide has a folded conformation with its Lys residue trans and Glu residue gauche−gauche+. The side chains of the Lys and Glu residues correspond to all trans and folded (g−g−g−) conformations respectively. The terminal carboxyl group forms hydrogen bonds with the ξ-amino group of the lysine side chain. The head-to-tail interaction often seen in peptide crystals is absent in the present structure. In the extended crystal structure water molecules form channels along the b direction and are enclosed within helically arranged hydrogen bonds formed by the lysine side chain and the peptide backbone.

Comparative study on damage identification from lso-Eigen-Value-Change contours and smeared damage model

The paper proposes two methodologies for damage identification from measured natural frequencies of a contiguously damaged reinforced concrete beam, idealised with distributed damage model. The first method identifies damage from Iso-Eigen-Value-Change contours, plotted between pairs of different frequencies. The performance of the method is checked for a wide variation of damage positions and extents. The method is also extended to a discrete structure in the form of a five-storied shear building and the simplicity of the method is demonstrated. The second method is through smeared damage model, where the damage is assumed constant for different segments of the beam and the lengths and centres of these segments are the known inputs. First-order perturbation method is used to derive the relevant expressions. Both these methods are based on distributed damage models and have been checked with experimental program on simply supported reinforced concrete beams, subjected to different stages of symmetric and un-symmetric damages. The results of the experiments are encouraging and show that both the methods can be adopted together in a damage identification scenario.

Reliability models for existing structures based on dynamic state estimation and data based asymptotic extreme value analysis

The problem of time variant reliability analysis of existing structures subjected to stationary random dynamic excitations is considered. The study assumes that samples of dynamic response of the structure, under the action of external excitations, have been measured at a set of sparse points on the structure. The utilization of these measurements m in updating reliability models, postulated prior to making any measurements, is considered. This is achieved by using dynamic state estimation methods which combine results from Markov process theory and Bayes' theorem. The uncertainties present in measurements as well as in the postulated model for the structural behaviour are accounted for. The samples of external excitations are taken to emanate from known stochastic models and allowance is made for ability (or lack of it) to measure the applied excitations. The future reliability of the structure is modeled using expected structural response conditioned on all the measurements made. This expected response is shown to have a time varying mean and a random component that can be treated as being weakly stationary. For linear systems, an approximate analytical solution for the problem of reliability model updating is obtained by combining theories of discrete Kalman filter and level crossing statistics. For the case of nonlinear systems, the problem is tackled by combining particle filtering strategies with data based extreme value analysis. In all these studies, the governing stochastic differential equations are discretized using the strong forms of Ito-Taylor's discretization schemes. The possibility of using conditional simulation strategies, when applied external actions are measured, is also considered. The proposed procedures are exemplifiedmby considering the reliability analysis of a few low-dimensional dynamical systems based on synthetically generated measurement data. The performance of the procedures developed is also assessed based on a limited amount of pertinent Monte Carlo simulations. (C) 2010 Elsevier Ltd. All rights reserved.

Solid-state chemistry of high-temperature oxide superconductors:The experimental situation

High-temperature superconductivity in oxides of the type(La, Ln)2?xBax(Sr)xCuO4, Y(Ln)Ba2Cu3O7??, La3?xBa3+xCu6O14, and related systems is discussed with emphasis on aspects related to experimental solid-state chemistry. All of these oxides possess perovskite-related structures. Oxygen-excess and La-deficient La2CuO4 also exhibit superconductivity in the 20�40 K just as La2?xBax(Srx)CuO4; these oxides are orthorhombic in the superconductivity phase. The crucial role of oxygen stoichiometry in the superconductivity ofYBa2Cu3O7?? (Tc = 95 ± 5K) is examined; this oxide remains orthorhombic up to ? ? 0.6 and becomes tetragonal and nonsuperconducting beyond this value of ?. Oxygen stoichiometry in this and related oxides has to be understood in terms of structure and disorder. The structure of La3?xBa3+xCu6O14 is related to that of YBa2Cu3O7, the orthorhombic structure manifesting itself when the population of O1 oxygens (along the Cusingle bondOsingle bondCu chains) is preponderant compared to that of O5 oxygens (along thea-axis); nearly equal populations of O1 and O5 sites give rise to the tetragonal structure. A transition from a high-Tc (95 K) superconductivity regime to a low-Tc (not, vert, similar60 K) regime occurs in YBa2Cu3O7?? accompanying a change in ?. There is no evidence for Cu3+ in these nominally mixed valent copper oxides. Instead, holes are present on oxygens giving rise to O? or O2?2 species, the concentration of these species increasing with the lowering of temperature. Certain interesting aspects of the superconducting oxides such as domain or twin boundaries, Raman spectra, microwave absorption, and anomalous high-temperature resistivity drops are presented along with the important material parameters. Preparative aspects of the superconducting oxides are briefly discussed. Phase transitions seem to occur atTc as well as at not, vert, similar240 K in YBa2Cu3O7.