The Adomian method applied to some extraordinary differential equations

Several ''extraordinary'' differential equations are considered for their solutions via the decomposition method of Adomian. Verifications are made with the solutions obtained by other methods.

Thomas-Fermi Method For Particles Obeying Generalized Exclusion Statistics

We use the Thomas-Fermi method to examine the thermodynamics of particles obeying Haldane exclusion statistics. Specifically, we study Calogero-Sutherland particles placed in a given external potential in one dimension. For the case of a simple harmonic potential (constant density of states), we obtain the exact one-particle spatial density and a {\it closed} form for the equation of state at finite temperature, which are both new results. We then solve the problem of particles in a $x^{2/3} ~$ potential (linear density of states) and show that Bose-Einstein condensation does not occur for any statistics other than bosons.

Graphical method to determine the parameters of the two-parameter fracture model for concrete

To evaluate the parameters in the two-parameter fracture model, i.e. the critical stress intensity factor and critical crack tip opening displacement for the fracture of plain concrete in Mode 1 for the given test configuration and geometry, considerable computational effort is necessary. A simple graphical method has been proposed using normalized fracture parameters for the three-point bend (3PB) notched specimen and the double-edged notched (DEN) specimen. A similar graphical method is proposed to compute the maximum load carrying capacity of a specimen, using the critical fracture parameters both for 3PB and DEN configurations.

Improving the convergence rate of the inverse iteration method

Solution of generalized eigenproblem, K phi = lambda M phi, by the classical inverse iteration method exhibits slow convergence for some eigenproblems. In this paper, a modified inverse iteration algorithm is presented for improving the convergence rate. At every iteration, an optimal linear combination of the latest and the preceding iteration vectors is used as the input vector for the next iteration. The effectiveness of the proposed algorithm is demonstrated for three typical eigenproblems, i.e. eigenproblems with distinct, close and repeated eigenvalues. The algorithm yields 29, 96 and 23% savings in computational time, respectively, for these problems. The algorithm is simple and easy to implement, and this renders the algorithm even more attractive.

Block coherence: a method for measuring the interdependence between two blocks of neurobiological time series

Multisensor recordings are becoming commonplace. When studying functional connectivity between different brain areas using such recordings, one defines regions of interest, and each region of interest is often characterized by a set (block) of time series. Presently, for two such regions, the interdependence is typically computed by estimating the ordinary coherence for each pair of individual time series and then summing or averaging the results over all such pairs of channels (one from block 1 and other from block 2). The aim of this paper is to generalize the concept of coherence so that it can be computed for two blocks of non-overlapping time series. This quantity, called block coherence, is first shown mathematically to have properties similar to that of ordinary coherence, and then applied to analyze local field potential recordings from a monkey performing a visuomotor task. It is found that an increase in block coherence between the channels from V4 region and the channels from prefrontal region in beta band leads to a decrease in response time.

Real-time Computer Simulation of Three Dimensional Elastostatics using the Finite Point Method

Real-time simulation of deformable solids is essential for some applications such as biological organ simulations for surgical simulators. In this work, deformable solids are approximated to be linear elastic, and an easy and straight forward numerical technique, the Finite Point Method (FPM), is used to model three dimensional linear elastostatics. Graphics Processing Unit (GPU) is used to accelerate computations. Results show that the Finite Point Method, together with GPU, can compute three dimensional linear elastostatic responses of solids at rates suitable for real-time graphics, for solids represented by reasonable number of points.

Preparation of the layered double hydroxide (LDH) LiAl2(OH)(7)center dot 2H(2)O, by gel to crystallite conversion and a hydrothermal method, and its conversion to lithium aluminates

A layered double hydroxide (LDH) with chemical composition LiAl2(OH)(7) . 2H(2)O was prepared via a wet chemical route of gel to crystallite (G-C) conversion at 80 degrees C involving the reaction of hydrated alumina gel, Al2O3.yH(2)O (80 < y < 120) with LiOH (Li2O/Al2O3 greater than or equal to 0.5) in presence of hydrophilic solvents such as ethanol under refluxing conditions. The hydrothermal synthesis was carried out using the same reactants by heating to less than or equal to 140 degrees C in a Teflon-lined autoclave under autogenerated pressure (less than or equal to 20 MPa). Transmission electron microscopy showed needle-shaped aggregates of size 0.04-0.1 mu m for the gel to crystallite conversion product, whereas the hydrothermal products consisted of individual lamellar crystallites of size 0.2-0.5 mu m with hexagonal morphology. The LDH prepared through the gel to crystallite conversion could be converted into LiAl(OH)(4) . H2O or LiAl(OH)(3)NO3 . H2O by imbibition of LiOH or LiNO3, respectively, under hydrothermal conditions. Thermal decomposition of LDH above 1400 degrees C gave rise to LiAl5O8 accompanied by the evaporation of Li2O. LiAl(OH)(4) . H2O and LiAl(OH)(3)NO3 . H2O decomposed in the temperature range 400-1000 degrees C to alpha- or beta-LiAlO2. The compositional dependence of the product, the intermediate phases formed during the heat treatment and the possible reactions involved are described in detail.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Comparison of electrical resistivity by geophysical method and neutron probe logging for soil moisture monitoring in a forested watershed

Geophysical methods are becoming more popular nowadays in the field of hydrology due to their time and space efficiency. So an attempt has been made here to relate electrical resistivity with soil moisture content in the field. The experiments were carried out in an experimental watershed `Mulehole' in southern India, which is a forested watershed with approximately 80% red soil. Five auger holes were drilled to perform the soil moisture and electrical resistivity measurements in a toposequence having red and black soils, with sandy weathered soil at the bottom. Soil moisture was measured using neutron probe and electrical resistivity was measured using electrical logging tool. The results indicate that electrical resistivity measurements can be used to measure soil moisture content for red soils only.

Extended Muskingum method for flood routing

Routing of floods is essential to control the flood flow at the flood control station such that it is within the specified safe limit. In this paper, the applicability of the extended Muskingum method is examined for routing of floods for a case study of Hirakud reservoir, Mahanadi river basin, India. The inflows to the flood control station are of two types-one controllable which comprises of reservoir releases for power and spill and the other is uncontrollable which comprises of inflow from lower tributaries and intermediate catchment between the reservoir and the flood control station. Muskingum model is improved to incorporate multiple sources of inflows and single outflow to route the flood in the reach. Instead of time lag and prismoidal flow parameters, suitable coefficients for various types of inflows were derived using Linear Programming. Presently, the decisions about operation of gates of Hirakud dam are being taken once in 12 h during floods. However, four time intervals of 24, 18, 12 and 6 h are examined to test the sensitivity of the routing time interval on the computed flood flow at the flood control station. It is observed that mean relative error decreases with decrease in routing interval both for calibration and testing phase. It is concluded that the extended Muskingum method can be explored for similar reservoir configurations such as Hirakud reservoir with suitable modifications. (C) 2010 International Association of Hydro-environment Engineering and Research. Asia Pacific Division. Published by Elsevier By. All rights reserved.

A checking method for probabilistic seismic-hazard assessment: case studies on three cities

The conventional Cornell's source-based approach of probabilistic seismic-hazard assessment (PSHA) has been employed all around the world, whilst many studies often rely on the use of computer packages such as FRISK (McGuire FRISK-a computer program for seismic risk analysis. Open-File Report 78-1007, United States Geological Survey, Department of Interior, Washington 1978) and SEISRISK III (Bender and Perkins SEISRISK III-a computer program for seismic hazard estimation, Bulletin 1772. United States Geological Survey, Department of Interior, Washington 1987). A ``black-box'' syndrome may be resulted if the user of the software does not have another simple and robust PSHA method that can be used to make comparisons. An alternative method for PSHA, namely direct amplitude-based (DAB) approach, has been developed as a heuristic and efficient method enabling users to undertake their own sanity checks on outputs from computer packages. This paper experiments the application of the DAB approach for three cities in China, Iran, and India, respectively, and compares with documented results computed by the source-based approach. Several insights regarding the procedure of conducting PSHA have also been obtained, which could be useful for future seismic-hazard studies.

Discovery of Application Workloads from Network File Traces

An understanding of application I/O access patterns is useful in several situations. First, gaining insight into what applications are doing with their data at a semantic level helps in designing efficient storage systems. Second, it helps create benchmarks that mimic realistic application behavior closely. Third, it enables autonomic systems as the information obtained can be used to adapt the system in a closed loop.All these use cases require the ability to extract the application-level semantics of I/O operations. Methods such as modifying application code to associate I/O operations with semantic tags are intrusive. It is well known that network file system traces are an important source of information that can be obtained non-intrusively and analyzed either online or offline. These traces are a sequence of primitive file system operations and their parameters. Simple counting, statistical analysis or deterministic search techniques are inadequate for discovering application-level semantics in the general case, because of the inherent variation and noise in realistic traces.In this paper, we describe a trace analysis methodology based on Profile Hidden Markov Models. We show that the methodology has powerful discriminatory capabilities that enable it to recognize applications based on the patterns in the traces, and to mark out regions in a long trace that encapsulate sets of primitive operations that represent higher-level application actions. It is robust enough that it can work around discrepancies between training and target traces such as in length and interleaving with other operations. We demonstrate the feasibility of recognizing patterns based on a small sampling of the trace, enabling faster trace analysis. Preliminary experiments show that the method is capable of learning accurate profile models on live traces in an online setting. We present a detailed evaluation of this methodology in a UNIX environment using NFS traces of selected commonly used applications such as compilations as well as on industrial strength benchmarks such as TPC-C and Postmark, and discuss its capabilities and limitations in the context of the use cases mentioned above.

Modeling of concrete cracking-A hybrid technique of using displacement discontinuity element method and direct boundary element method

This paper presents a new approach by making use of a hybrid method of using the displacement discontinuity element method and direct boundary element method to model concrete cracking by incorporating fictitious crack model. Fracture mechanics approach is followed using the Hillerborg's fictitious crack model. A boundary element based substructure method and a hybrid technique of using displacement discontinuity element method and direct boundary element method are compared in this paper. In order to represent the process zone ahead of the crack, closing forces are assumed to act in such a way that they obey a linear normal stress-crack opening displacement law. Plain concrete beams with and without initial crack under three-point loading were analyzed by both the methods. The numerical results obtained were shown to agree well with the results from existing finite element method. The model is capable of reproducing the whole range of load-deflection response including strain-softening and snap-back behavior as illustrated in the numerical examples. (C) 2011 Elsevier Ltd. All rights reserved.

A method for the calculation of the adsorbed phase volume and pseudo-saturation pressure from adsorption isotherm data on activated carbon

We propose a new method for evaluating the adsorbed phase volume during physisorption of several gases on activated carbon specimens. We treat the adsorbed phase as another equilibrium phase which satisfies the Gibbs equation and hence assume that the law of rectilinear diameters is applicable. Since invariably the bulk gas phase densities are known along measured isotherms, the constants of the adsorbed phase volume can be regressed from the experimental data. We take the Dubinin-Astakhov isotherm as the model for verifying our hypothesis since it is one of the few equations that accounts for adsorbed phase volume changes. In addition, the pseudo-saturation pressure in the supercritical region is calculated by letting the index of the temperature term in Dubinin's equation to be temperature dependent. Based on over 50 combinations of activated carbons and adsorbates (nitrogen, oxygen, argon, carbon dioxide, hydrocarbons and halocarbon refrigerants) it is observed that the proposed changes fit experimental data quite well.