An FETI-preconditioned conjuerate gradient method for large-scale stochastic finite element problems

In the spectral stochastic finite element method for analyzing an uncertain system. the uncertainty is represented by a set of random variables, and a quantity of Interest such as the system response is considered as a function of these random variables Consequently, the underlying Galerkin projection yields a block system of deterministic equations where the blocks are sparse but coupled. The solution of this algebraic system of equations becomes rapidly challenging when the size of the physical system and/or the level of uncertainty is increased This paper addresses this challenge by presenting a preconditioned conjugate gradient method for such block systems where the preconditioning step is based on the dual-primal finite element tearing and interconnecting method equipped with a Krylov subspace reusage technique for accelerating the iterative solution of systems with multiple and repeated right-hand sides. Preliminary performance results on a Linux Cluster suggest that the proposed Solution method is numerically scalable and demonstrate its potential for making the uncertainty quantification Of realistic systems tractable.

Entire Solutions of Hydrodynamical Equations with Exponential Dissipation

We consider a modification of the three-dimensional Navier-Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as e(vertical bar k vertical bar/kd) at high wavenumbers vertical bar k vertical bar. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than e-(C(k/kd) ln(vertical bar k vertical bar/kd)) for any C < 1/(2 ln 2). The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with C = C-* = 1/ ln 2. The same behavior with a universal constant C-* is conjectured for the Navier-Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier-Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.

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.

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.

Block Convolution Using Discrete Trigonometric Transforms and Discrete Fourier Transform

We derive expressions for convolution multiplication properties of discrete cosine transform II (DCT II) starting from equivalent discrete Fourier transform (DFT) representations. Using these expressions, a method for implementing linear filtering through block convolution in the DCT II domain is presented. For the case of nonsymmetric impulse response, additional discrete sine transform II (DST II) is required for implementing the filter in DCT II domain, where as for a symmetric impulse response, the additional transform is not required. Comparison with recently proposed circular convolution technique in DCT II domain shows that the proposed new method is computationally more efficient.

Breath Testing by Fourier Transform Infrared Spectroscopy for Solvent Intoxication Diagnostics

Technical or contaminated ethanol products are sometimes ingested either accidentally or on purpose. Typical misused products are black-market liquor and automotive products, e.g., windshield washer fluids. In addition to less toxic solvents, these liquids may contain the deadly methanol. Symptoms of even lethal solvent poisoning are often non-specific at the early stage. The present series of studies was carried out to develop a method for solvent intoxication breath diagnostics to speed up the diagnosis procedure conventionally based on blood tests. Especially in the case of methanol ingestion, the analysis method should be sufficiently sensitive and accurate to determine the presence of even small amounts of methanol from the mixture of ethanol and other less-toxic components. In addition to the studies on the FT-IR method, the Dräger 7110 evidential breath analyzer was examined to determine its ability to reveal a coexisting toxic solvent. An industrial Fourier transform infrared analyzer was modified for breath testing. The sample cell fittings were widened and the cell size reduced in order to get an alveolar sample directly from a single exhalation. The performance and the feasibility of the Gasmet FT-IR analyzer were tested in clinical settings and in the laboratory. Actual human breath screening studies were carried out with healthy volunteers, inebriated homeless men, emergency room patients and methanol-intoxicated patients. A number of the breath analysis results were compared to blood test results in order to approximate the blood-breath relationship. In the laboratory experiments, the analytical performance of the Gasmet FT-IR analyzer and Dräger 7110 evidential breath analyzer was evaluated by means of artificial samples resembling exhaled breath. The investigations demonstrated that a successful breath ethanol analysis by Dräger 7110 evidential breath analyzer could exclude any significant methanol intoxication. In contrast, the device did not detect very high levels of acetone, 1-propanol and 2-propanol in simulated breath. The Dräger 7110 evidential breath ethanol analyzer was not equipped to recognize the interfering component. According to the studies the Gasmet FT-IR analyzer was adequately sensitive, selective and accurate for solvent intoxication diagnostics. In addition to diagnostics, the fast breath solvent analysis proved feasible for controlling the ethanol and methanol concentration during haemodialysis treatment. Because of the simplicity of the sampling and analysis procedure, non-laboratory personnel, such as police officers or social workers, could also operate the analyzer for screening purposes.

Synthesis, spectral characterization, in-vitro microbiological evaluation and cytotoxic activities of novel macrocyclic bis hydrazone

A macrocyclic hydrazone Schiff base was synthesized by reacting 1,4-dicarbonyl phenyl dihydrazide with 2,6-diformyl-4-methyl phenol and a series of metal complexes with this new Schiff base were synthesized by reaction with Co(II), Ni(II) and Cu(II) metal salts. The Schiff base and its complexes have been characterized by elemental analyses, IR, H-1 NMR, UV-vis, FAB mass, ESR spectra, fluorescence, thermal, magnetic and molar conductance data. The analytical data reveal that the Co(II), Ni(II) and Cu(II) complexes possess 2:1 metal-ligand ratios. All the complexes are non-electrolytes in DMF and DMSO due to their low molar conductance values. Infrared spectral data suggest that the hydrazone Schiff base behaves as a hexadentate ligand with NON NON donor sequence towards the metal ions. The ESR spectral data shows that the metal-ligand bond has considerable covalent character. The electrochemical behavior of the copper(II) complex was investigated by cyclic voltammetry. The Schiff base and its complexes have also been screened for their antibacterial (Escherichia coli, Staphylococcus aureus, Shigella dysentery, Micrococcus, Bacillus subtilis, Bacillus cereus and Pseudomonas aeruginosa) and antifungal activities (Aspergillus niger, Penicillium and Candida albicans) by MIC method. The brine shrimp bioassay was also carried out to study their in-vitro cytotoxic properties. (C) 2009 Elsevier Masson SAS. All rights reserved.

Synthesis of BaSO4 nanoparticles by precipitation method using sodium hexa metaphosphate as a stabilizer

We describe the synthesis and structure of Barium sulfate nanoparticles by precipitation method in the presence of water soluble inorganic stabilizing agent, sodium hexametaphosphate, (NaPO3)(6). The structural parameters were refined by the Rietveld refinement method using powder X-ray diffraction data. Barium sulfate nanoparticles were crystallized in the orthorhombic structure with space group Pbnm (No. 62) having the lattice parameters a = 7.215(1) (angstrom), b = 8.949(1) (angstrom) and c = 5.501 (1) (angstrom) respectively. Transmission electron microscopy study reveals that the nanoparticles are size range, 30-50 nm. Fourier transform infrared spectra showed distinct absorption due to the SO42- moiety at 1115 and 1084 cm(-1) indicating formation of barium sulfate nanoparticles free from the phosphate group from the stabilizer used in the synthesis. (C) 2009 Elsevier Ltd. All rights reserved.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Automatic Road Extraction using High Resolution Satellite Image Based on Texture Progressive Analysis and Normalized Cut Method

In this paper the approach for automatic road extraction for an urban region using structural, spectral and geometric characteristics of roads has been presented. Roads have been extracted based on two levels: Pre-processing and road extraction methods. Initially, the image is pre-processed to improve the tolerance by reducing the clutter (that mostly represents the buildings, parking lots, vegetation regions and other open spaces). The road segments are then extracted using Texture Progressive Analysis (TPA) and Normalized cut algorithm. The TPA technique uses binary segmentation based on three levels of texture statistical evaluation to extract road segments where as, Normalizedcut method for road extraction is a graph based method that generates optimal partition of road segments. The performance evaluation (quality measures) for road extraction using TPA and normalized cut method is compared. Thus the experimental result show that normalized cut method is efficient in extracting road segments in urban region from high resolution satellite image.

Spectral unmixing with estimated adaptive endmember index using Extended Support Vector Machine

The most difficult operation in the flood inundation mapping using optical flood images is to separate fully inundated areas from the ‘wet’ areas where trees and houses are partly covered by water. This can be referred as a typical problem the presence of mixed pixels in the images. A number of automatic information extraction image classification algorithms have been developed over the years for flood mapping using optical remote sensing images. Most classification algorithms generally, help in selecting a pixel in a particular class label with the greatest likelihood. However, these hard classification methods often fail to generate a reliable flood inundation mapping because the presence of mixed pixels in the images. To solve the mixed pixel problem advanced image processing techniques are adopted and Linear Spectral unmixing method is one of the most popular soft classification technique used for mixed pixel analysis. The good performance of linear spectral unmixing depends on two important issues, those are, the method of selecting endmembers and the method to model the endmembers for unmixing. This paper presents an improvement in the adaptive selection of endmember subset for each pixel in spectral unmixing method for reliable flood mapping. Using a fixed set of endmembers for spectral unmixing all pixels in an entire image might cause over estimation of the endmember spectra residing in a mixed pixel and hence cause reducing the performance level of spectral unmixing. Compared to this, application of estimated adaptive subset of endmembers for each pixel can decrease the residual error in unmixing results and provide a reliable output. In this current paper, it has also been proved that this proposed method can improve the accuracy of conventional linear unmixing methods and also easy to apply. Three different linear spectral unmixing methods were applied to test the improvement in unmixing results. Experiments were conducted in three different sets of Landsat-5 TM images of three different flood events in Australia to examine the method on different flooding conditions and achieved satisfactory outcomes in flood mapping.

Wave propagation analysis in carbon nanotube embedded composite using wavelet based spectral finite elements

In this paper, elastic wave propagation is studied in a nanocomposite reinforced with multiwall carbon nanotubes (CNTs). Analysis is performed on a representative volume element of square cross section. The frequency content of the exciting signal is at the terahertz level. Here, the composite is modeled as a higher order shear deformable beam using layerwise theory, to account for partial shear stress transfer between the CNTs and the matrix. The walls of the multiwall CNTs are considered to be connected throughout their length by distributed springs, whose stiffness is governed by the van der Waals force acting between the walls of nanotubes. The analyses in both the frequency and time domains are done using the wavelet-based spectral finite element method (WSFEM). The method uses the Daubechies wavelet basis approximation in time to reduce the governing PDE to a set of ODEs. These transformed ODEs are solved using a finite element (FE) technique by deriving an exact interpolating function in the transformed domain to obtain the exact dynamic stiffness matrix. Numerical analyses are performed to study the spectrum and dispersion relations for different matrix materials and also for different beam models. The effects of partial shear stress transfer between CNTs and matrix on the frequency response function (FRF) and the time response due to broadband impulse loading are investigated for different matrix materials. The simultaneous existence of four coupled propagating modes in a double-walled CNT-composite is also captured using modulated sinusoidal excitation.

Synthesis, characterization and photoluminescence properties of Gd2O3:Eu3+ nanophosphors prepared by solution combustion method

Gd2O3:Eu3+ (0.5-8.0 mol%) nanophosphors have been prepared by low temperature solution combustion method using metal nitrates as oxidizers and oxalyl dihydrazide (ODH) as a fuel. The phosphors are well characterized by powder X-ray diffraction (PXRD), scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FT-IR) and photoluminescence (PL) techniques. PXRD patterns of as-formed and calcined (800 degrees C, 3 h) Gd2O3 powders exhibit monoclinic phase with mean crystallite sizes ranging from 20 to 50 nm. Eu3+ doping changes the structure from monoclinic to mixed phase of monoclinic and cubic. SEM micrographs shows the products are foamy, agglomerated and fluffy in nature due to the large amount of gases liberated during combustion reaction. Upon 254 nm excitation the photoluminescence of the Gd2O3:Eu3+ particles show red emission at 611 nm corresponding to D-5(0)-> F-7(2) transition. It is observed that PL intensity increases with calcination temperature. This might be attributed to better crystallization and eliminates the defects, which serve as centers of non-radiative relaxation for nanomaterials. It is observed that the optical energy gap (E-g) is widened with increase Eu3+ content. (C) 2010 Elsevier B.V. All rights reserved.

A hybrid electrochemical-thermal method for the preparation of large ZnO nanoparticles

A simple and efficient two-step hybrid electrochemical-thermal route was developed for the synthesis of large quantity of ZnO nanoparticles using aqueous sodium bicarbonate electrolyte and sacrificial Zn anode and cathode in an undivided cell under galvanostatic mode at room temperature. The bath concentration and current density were varied from 30 to 120 mmol and 0.05 to 1.5 A/dm(2). The electrochemically generated precursor was calcined for an hour at different range of temperature from 140 to 600 A degrees C. The calcined samples were characterized by XRD, SEM/EDX, TEM, TG-DTA, FT-IR, and UV-Vis spectral methods. Rietveld refinement of X-ray data indicates that the calcined compound exhibits hexagonal (Wurtzite) structure with space group of P63mc (No. 186). The crystallite sizes were in the range of 22-75 nm based on Debye-Scherrer equation. The TEM results reveal that the particle sizes were in the order of 30-40 nm. The blue shift was noticed in UV-Vis absorption spectra, the band gaps were found to be 5.40-5.11 eV. Scanning electron micrographs suggest that all the samples were randomly oriented granular morphology.

Spectral assignments based on the 13C NMR spectra of mixed liquid crystals of opposite diamagnetic anisotropies

The proton-decoupled 13C NMR spectra of mixtures of liquid crystals with opposite diamagnetic anisotropies have been studied in the natural abundance of 13C. A new method to assign the spectral lines to specific carbons in the liquid crystalline phase has been developed. For this purpose, the assignments of lines in the isotropic media are required, and they were obtained from two-dimensional hetero-COSY experiments. From the spectra in the �critical� mixtures where both the orientations of the liquid crystal directors, with the alignments along and perpendicular to the direction of the magnetic field, �coexist,� the 13C chemical-shift anisotropies have been determined, assuming uniaxial symmetry.