A parallel multistep predictor-corrector algorithm for solving ordinary differential equations

In this paper we give a generalized predictor-corrector algorithm for solving ordinary differential equations with specified initial values. The method uses multiple correction steps which can be carried out in parallel with a prediction step. The proposed method gives a larger stability interval compared to the existing parallel predictor-corrector methods. A method has been suggested to implement the algorithm in multiple processor systems with efficient utilization of all the processors.

Analytical modeling of quantum threshold voltage for triple gate MOSFET

In this work a physically based analytical quantum threshold voltage model for the triple gate long channel metal oxide semiconductor field effect transistor is developed The proposed model is based on the analytical solution of two-dimensional Poisson and two-dimensional Schrodinger equation Proposed model is extended for short channel devices by including semi-empirical correction The impact of effective mass variation with film thicknesses is also discussed using the proposed model All models are fully validated against the professional numerical device simulator for a wide range of device geometries (C) 2010 Elsevier Ltd All rights reserved

Effect of wall thickness on the end corrections of the extended inlet and outlet of a double-tuned expansion chamber

Simple expansion chambers, the simplest of the muffler configurations, have very limited practical application due to the presence of periodic troughs in the transmission loss spectrum which drastically lower the overall transmission loss of the muffler. Tuned extended inlet and outlet can be designed to nullify three-fourths of these troughs, making use of the plane wave theory. These cancellations would not occur unless one altered the geometric lengths for the extended tube in order to incorporate the effect of evanescent higher-order modes (multidimensional effect) through end corrections or lumped inertance approximation at the area discontinuities or junctions. End corrections of the extended inlet and outlet have been studied by several researchers. However the effect of wall thickness of the inlet/outlet duct on end correction has not been studied explicitly. This has significant effect on the tuning of an extended inlet/outlet expansion chamber. It is investigated here experimentally as well as numerically (through use of 3-D FEM software) for stationary medium. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

Optimal explicit guidance of multistage launch vehicle along three-dimensional trajectory

Details of an efficient optimal closed-loop guidance algorithm for a three-dimensional launch are presented with simulation results. Two types of orbital injections, with either true anomaly or argument of perigee being free at injection, are considered. The resulting steering-angle profile under the assumption of uniform gravity lies in a canted plane which transforms a three-dimensional problem into an equivalent two-dimensional one. Effects of thrust are estimated using a series in a recursive way. Encke's method is used to predict the trajectory during powered flight and then to compute the changes due to actual gravity using two gravity-related vectors. Guidance parameters are evaluated using the linear differential correction method. Optimality of the algorithm is tested against a standard ground-based trajectory optimization package. The performance of the algorithm is tested for accuracy, robustness, and efficiency for a sun-synchronous mission involving guidance for a multistage vehicle that requires large pitch and yaw maneuver. To demonstrate applicability of the algorithm to a range of missions, injection into a geostationary transfer orbit is also considered. The performance of the present algorithm is found to be much better than others.

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)

Chemical potential of oxygen for iron-rutile-ilmenite and iron-ilmenite-ulvospinel equilibria

The chemical potential of oxygen corresponding to the iron-rutile-ilmenite (IRI) and iron-ilmenite-ulvospinel (IIU) equilibria has been measured employing solid-state galvanic cells,$$Pt, Fe + TiO_2 + FeTiO_3 //(Y_2 O_3 ) ZrO_2 //Fe + FeO, Pt$$ and $${\text{Pt, Fe + FeTiO}}_{\text{3}} {\text{ + Fe}}_{\text{2}} {\text{TiO}}_{\text{4}} {\text{//(Y}}_{\text{2}} {\text{0}}_{\text{3}} {\text{) ZrO}}_{\text{2}} {\text{//Fe + FeO, Pt}}$$ in the temperature range of 875 to 1275 K and 900 to 1373 K, respectively. The cells are written such that the right-hand electrodes are positive. The electromotive force (emf) of both the cells was found to be reversible and to vary linearly with temperature over the entire range of measurement. The chemical potential of oxygen for IRI equilibrium is represented by Δμo2(IRI) = -550,724 - 29.445T + 20.374T InT(±210) J mol−1 (875 <-T<- 1184 K) = -620,260 + 369.593T - 27.716T lnT(±210) J mol−1 (1184 <-T<- 1275 K) and that for IIU equilibrium by Δμo2(IIU) = -501,800 - 49.035T + 20.374T lnT(±210) J mol−1 (900 <-T<- 1184 K) = -571,336 + 350.003T− 27.716T lnT(=−210) J mol-1 (1184 <-T<- 1373 K) The standard Gibbs energy changes for IRI and IIU equilibria have been deduced from the measured oxygen potentials. Since ilmenite contains small amounts of Ti³+ ions, a correction for the activity of FeTiO3 has been incorporated by assuming ideal mixing on each cation sublattice in the FeTiO3-Ti2O3 system. Similarly, the ulvospinel contains some Fe³+ ions and a correction for the activity of Fe2TiO4 has been included by modeling the Fe2TiO4-Fe3O4 system. The third-law analysis of the results obtained for IRI equilibrium gives ΔH 298 0 = -575 (±1.0) kJ mol-1 and for IIU equilibrium yields ΔH 298 0 = -523.7 (±0.7) kJ mol−1}. The present results suggest that Fe2+ and Ti4+ cations mix almost ideally on the octahedral site of spinel lattice in Fe2TiO4, giving rise to a configurational contribution of 2R In 2 (11.5256 J mol-1 K-1) to the entropy of Fe2TiO4.

Impaired conflict monitoring in Parkinson's disease patients during an oculomotor redirect task

Fallibility is inherent in human cognition and so a system that will monitor performance is indispensable. While behavioral evidence for such a system derives from the finding that subjects slow down after trials that are likely to produce errors, the neural and behavioral characterization that enables such control is incomplete. Here, we report a specific role for dopamine/basal ganglia in response conflict by accessing deficits in performance monitoring in patients with Parkinson's disease. To characterize such a deficit, we used a modification of the oculomotor countermanding task to show that slowing down of responses that generate robust response conflict, and not post-error per se, is deficient in Parkinson's disease patients. Poor performance adjustment could be either due to impaired ability to slow RT subsequent to conflicts or due to impaired response conflict recognition. If the latter hypothesis was true, then PD subjects should show evidence of impaired error detection/correction, which was found to be the case. These results make a strong case for impaired performance monitoring in Parkinson's patients.

Nonaxisymmetric unsteady motion over a rotating disk in the presence of free convection and magnetic field

The nonaxisymmetric unsteady motion produced by a buoyancy-induced cross-flow of an electrically conducting fluid over an infinite rotating disk in a vertical plane and in the presence of an applied magnetic field normal to the disk has been studied. Both constant wall and constant heat flux conditions have been considered. It has been found that if the angular velocity of the disk and the applied magnetic field squared vary inversely as a linear function of time (i.e. as (1??t*)?1, the governing Navier-Stokes equation and the energy equation admit a locally self-similar solution. The resulting set of ordinary differential equations has been solved using a shooting method with a generalized Newton's correction procedure for guessed boundary conditions. It is observed that in a certain region near the disk the buoyancy induced cross-flow dominates the primary von Karman flow. The shear stresses induced by the cross-flow are found to be more than these of the primary flow and they increase with magnetic parameter or the parameter ? characterizing the unsteadiness. The velocity profiles in the x- and y-directions for the primary flow at any two values of the unsteady parameter ? cross each other towards the edge of the boundary layer. The heat transfer increases with the Prandtl number but reduces with the magnetic parameter.

Carbon Dioxide emissions from an Acacia plantation on peatland in Sumatra, Indonesia

Abstract. Peat surface CO2 emission, groundwater table depth and peat temperature were monitored for two years along transects in an Acacia plantation on thick tropical peat (>4 m) in Sumatra, Indonesia. A total of 2300 emission measurements were taken at 144 locations. The autotrophic root respiration component of the CO2 emission was separated from heterotrophic emissions caused by peat oxidation in three ways: (i) by comparing CO2 emissions within and beyond the tree rooting zone, (ii) by comparing CO2 emissions with and without peat trenching (i.e. cutting any roots remaining in the peat beyond the tree rooting zone), and (iii) by comparing CO2 emissions before and after Acacia tree harvesting. On average, the contribution of root respiration to daytime CO2 emission is 21 % along transects in mature tree stands. At locations 0.5 m from trees this is up to 80 % of the total emissions, but it is negligible at locations more than 1.3 m away. This means that CO2 emission measurements well away from trees are free of any root respiration contribution and thus represent only peat oxidation emission. We find daytime mean annual CO2 emission from peat oxidation alone of 94 t ha−1 yr−1 at a mean water table depth of 0.8 m, and a minimum emission value of 80 t ha−1 yr−1 after correction for the effect of diurnal temperature fluctuations, which resulted in a 14.5 % reduction of the daytime emission. There is a positive correlation between mean long-term water table depths and peat oxidation CO2 emission. However, no such relation is found for instantaneous emission/water table depth within transects and it is clear that factors other than water table depth also affect peat oxidation and total CO2 emissions. The increase in the temperature of the surface peat due to plantation development may explain over 50 % of peat oxidation emissions.

Water quality estimation by optical remote sensing in boreal lakes

In lake-rich regions, the gathering of information about water quality is challenging because only a small proportion of the lakes can be assessed each year by conventional methods. One of the techniques for improving the spatial and temporal representativeness of lake monitoring is remote sensing from satellites and aircrafts. The experimental material included detailed optical measurements in 11 lakes, air- and spaceborne remote sensing measurements with concurrent field sampling, automatic raft measurements and a national dataset of routine water quality measurements from over 1100 lakes. The analyses of the spatially high-resolution airborne remote sensing data from eutrophic and mesotrophic lakes showed that one or a few discrete water quality observations using conventional monitoring can yield a clear over- or underestimation of the overall water quality in a lake. The use of TM-type satellite instruments in addition to routine monitoring results substantially increases the number of lakes for which water quality information can be obtained. The preliminary results indicated that coloured dissolved organic matter (CDOM) can be estimated with TM-type satellite instruments, which could possibly be utilised as an aid in estimating the role of lakes in global carbon budgets. Based on the results of reflectance modelling and experimental data, MERIS satellite instrument has optimal or near-optimal channels for the estimation of turbidity, chlorophyll a and CDOM in Finnish lakes. MERIS images with a 300 m spatial resolution can provide water quality information in different parts of large and medium-sized lakes, and in filling in the gaps resulting from conventional monitoring. Algorithms that would not require simultaneous field data for algorithm training would increase the amount of remote sensing-based information available for lake monitoring. The MERIS Boreal Lakes processor, trained with the optical data and concentration ranges provided by this study, enabled turbidity estimations with good accuracy without the need for algorithm correction with field measurements, while chlorophyll a and CDOM estimations require further development of the processor. The accuracy of interpreting chlorophyll a via semi empirical algorithms can be improved by classifying lakes prior to interpretation according to their CDOM level and trophic status. Optical modelling indicated that the spectral diffuse attenuation coefficient can be estimated with reasonable accuracy from the measured water quality concentrations. This provides more detailed information on light attenuation from routine monitoring measurements than is available through the Secchi disk transparency. The results of this study improve the interpretation of lake water quality by remote sensing and encourage the use of remote sensing in lake monitoring.

Microbial identification by detection of ligation probes on DNA microarray

Microbes in natural and artificial environments as well as in the human body are a key part of the functional properties of these complex systems. The presence or absence of certain microbial taxa is a correlate of functional status like risk of disease or course of metabolic processes of a microbial community. As microbes are highly diverse and mostly notcultivable, molecular markers like gene sequences are a potential basis for detection and identification of key types. The goal of this thesis was to study molecular methods for identification of microbial DNA in order to develop a tool for analysis of environmental and clinical DNA samples. Particular emphasis was placed on specificity of detection which is a major challenge when analyzing complex microbial communities. The approach taken in this study was the application and optimization of enzymatic ligation of DNA probes coupled with microarray read-out for high-throughput microbial profiling. The results show that fungal phylotypes and human papillomavirus genotypes could be accurately identified from pools of PCR amplicons generated from purified sample DNA. Approximately 1 ng/μl of sample DNA was needed for representative PCR amplification as measured by comparisons between clone sequencing and microarray. A minimum of 0,25 amol/μl of PCR amplicons was detectable from amongst 5 ng/μl of background DNA, suggesting that the detection limit of the test comprising of ligation reaction followed by microarray read-out was approximately 0,04%. Detection from sample DNA directly was shown to be feasible with probes forming a circular molecule upon ligation followed by PCR amplification of the probe. In this approach, the minimum detectable relative amount of target genome was found to be 1% of all genomes in the sample as estimated from 454 deep sequencing results. Signal-to-noise of contact printed microarrays could be improved by using an internal microarray hybridization control oligonucleotide probe together with a computational algorithm. The algorithm was based on identification of a bias in the microarray data and correction of the bias as shown by simulated and real data. The results further suggest semiquantitative detection to be possible by ligation detection, allowing estimation of target abundance in a sample. However, in practise, comprehensive sequence information of full length rRNA genes is needed to support probe design with complex samples. This study shows that DNA microarray has the potential for an accurate microbial diagnostic platform to take advantage of increasing sequence data and to replace traditional, less efficient methods that still dominate routine testing in laboratories. The data suggests that ligation reaction based microarray assay can be optimized to a degree that allows good signal-tonoise and semiquantitative detection.

A Feature-Based Solution to Forward Problem in Electrical Capacitance Tomography of Conductive Materials

A new feature-based technique is introduced to solve the nonlinear forward problem (FP) of the electrical capacitance tomography with the target application of monitoring the metal fill profile in the lost foam casting process. The new technique is based on combining a linear solution to the FP and a correction factor (CF). The CF is estimated using an artificial neural network (ANN) trained using key features extracted from the metal distribution. The CF adjusts the linear solution of the FP to account for the nonlinear effects caused by the shielding effects of the metal. This approach shows promising results and avoids the curse of dimensionality through the use of features and not the actual metal distribution to train the ANN. The ANN is trained using nine features extracted from the metal distributions as input. The expected sensors readings are generated using ANSYS software. The performance of the ANN for the training and testing data was satisfactory, with an average root-mean-square error equal to 2.2%.

Quantification of enantiomeric excess by H-1-detected heteronuclear refocusing and homonuclear multiple quantum NMR experiments

The omega(1)-heterodecoupled-C-13-filtered proton detected NMR experiments are reported for the accurate quantification of enantiomeric excess in chiral molecules embedded in chiral liquid crystal. The differential values of both H-1-H-1 and C-13-H-1 dipolar couplings in the direct dimension and only H-1-H-1 dipolar couplings in the indirect dimension enable unraveling of overlapped enantiomeric peaks. The creation of unequal C-13-bound proton signal for each enantiomer in the INEPT block and non-uniform excitation of coherences in homonuclear multiple quantum experiments do not yield accurate quantification of enantiomeric excess. In circumventing these difficulties, a coupling dependent intensity correction factor has been invoked. (C) 2010 Elsevier B.V. All rights reserved.

Autoregressive spectral array for graphical display of EEG data

A graphical display of the frequency content of background,electroencephalogram (EEG) activity is obtained by calculating the spectral estimates using autocorrelation autoregressive method and the classical Fourier transform method, Display of spectral content of consecutive data segments is made using hidden-line suppression technique so as to get a spectral array, The autoregressive spectral array (ASA) is found to be sensitive to baseline drift, Following baseline correction the autoregressive technique is found to be superior to the Fourier method of compressed spectral array (CSA) in detecting the transitions in the frequencies of the signal. The smoothed ASA gives a better picture of transitions and changes in the background activity, The ASA can be made to adapt to specific changes of dominant frequencies while eliminating unnecessary peaks in the spectrum. The utility,of the ASA for background EEG analysis is discussed,

Xenon in sodium Y zeolite. 2. Arrhenius relation, mechanism, and barrier height distribution for cage-to-cage diffusion

Various geometrical and energetic distribution functions and other properties connected with the cage-to-cage diffusion of xenon in sodium Y zeolite have been obtained from long molecular dynamics calculations. Analysis of diffusion pathways reveals two interesting mechanisms-surface-mediated and centralized modes for cage-to-cage diffusion. The surface-mediated mode of diffusion exhibits a small positive barrier, while the centralized diffusion exhibits a negative barrier for the sorbate to diffuse across the 12-ring window. In both modes, however, the sorbate has to be activated from the adsorption site to enable it to gain mobility. The centralized diffusion additionally requires the sorbate to be free of the influence of the surface of the cage as well. The overall rate for cage-to-cage diffusion shows an Arrhenius temperature dependence with E(a) = 3 kJ/mol. It is found that the decay in the dynamical correction factor occurs on a time scale comparable to the cage residence time. The distributions of barrier heights have been calculated. Functions reflecting the distribution of the sorbate-zeolite interaction at the window and the variations of the distance between the sorbate and the centers of the parent and daughter cages are presented.