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.

Method for measurement of the angles of 90-, 45-, 45-deg and 60-, 30-, 90-deg prisms

A simple technique is devised to measure the angles of 90-, 45-, 45-deg and 60-, 30-, 90-deg prisms without using expensive spectrometers, autocollimators, and angle gauges. The method can be extended to unpolished and opaque prisms made of materials other then glass. (C) 1997 Society of Photo-Optical instrumentation Engineers.

Method for measurement of the angles of an equilateral (60-deg) prism

A simple technique is devised io measure the angles of equilateral (60-deg) prisms, without using the expensive spectrometers, autocollimators, and angle gauges. The method can be extended to unpolished and opaque prisms made out of materials other than glass. (C) 1997 Society of Photo-Optical Instrumentation Engineers.

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.

The log delta log t method for the determination of the coefficient of consolidation

The behaviour of saturated soils undergoing consolidation is very complex, It may not follow Terzaghi's theory over the entire consolidation process, Different soils may behave in such a way as to fit into Terzaghi's theory over some specific stages of the consolidation process (percentage of consolidation), This may be one of the reasons for the difficulties faced by the existing curve-fitting procedures in obtaining the coefficient of consolidation, c(v). It has been shown that the slope of the initial linear portion of the theoretical log U-log T curve is constant over a wider range of degree of consolidation, U, when compared with the other methods in use, This initial well-defined straight line in the log U-log T plot intersects the U = 100% line at T = pi/4, which corresponds to U = 88.3%, The proposed log delta-log t method is based on this observation, which gives the value of c(v) through simple graphical construction, In the proposed method, which is more versatile, identification of the characteristic straight lines is very clear; the intersection of these lines is more precise and the method does not depend upon the initial compression for the determination of c(v).

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.

Parallel implementation of alternate quadrant interlocking factorisation method on star topology

This paper discusses the parallel implementation of the solution of a set of linear equations using the Alternative Quadrant Interlocking Factorisation Methods (AQIF), on a star topology. Both the AQIF and LU decomposition methods are mapped onto star topology on an IBM SP2 system, with MPI as the internode communicator. Performance parameters such as speedup, efficiency have been obtained through experimental and theoretical means. The studies demonstrate (i) a mismatch of 15% between the theoretical and experimental results, (ii) scalability of the AQIF algorithm, and (iii) faster executing AQIF algorithm.

Rapid method of consolidation testing

The determination of consolidation characteristics forms an important aspect in the design of foundations and other earth-retaining structures. The conventional consolidation test as originally proposed by Teaaghi takes considerable time (more than 15 days in highly compressible soils with low coefficient of consolidation) and effort. Any effort to reduce the duration of testing will be desirable from several considerations. In this paper, an attempt has been made to propose a rapid method of consolidation testing. In the proposed method, the next load increment is applied as soon as the necessary time required to identify the percent consolidation is reached and to evaluate the coefficient of consolidation by one of the popular curve-fitting procedures. The rectangular hyperbola method has been used to identify the percent consolidation reached after any load increment, and to determine the coefficient of consolidation, before making the next load increment. The time required to complete the test using the rapid consolidation method could be as low as 4-5 h compared with 1 or 2 weeks in the case of the conventional consolidation test.

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.

A differential geometric method for kinematic analysis of two- and three-degree-of-freedom rigid body motions

In this paper, we present a novel differential geometric characterization of two- and three-degree-of-freedom rigid body kinematics, using a metric defined on dual vectors. The instantaneous angular and linear velocities of a rigid body are expressed as a dual velocity vector, and dual inner product is defined on this dual vector, resulting in a positive semi-definite and symmetric dual matrix. We show that the maximum and minimum magnitude of the dual velocity vector, for a unit speed motion, can be obtained as eigenvalues of this dual matrix. Furthermore, we show that the tip of the dual velocity vector lies on a dual ellipse for a two-degree-of-freedom motion and on a dual ellipsoid for a three-degree-of-freedom motion. In this manner, the velocity distribution of a rigid body can be studied algebraically in terms of the eigenvalues of a dual matrix or geometrically with the dual ellipse and ellipsoid. The second-order properties of the two- and three-degree-of-freedom motions of a rigid body are also obtained from the derivatives of the elements of the dual matrix. This results in a definition of the geodesic motion of a rigid body. The theoretical results are illustrated with the help of a spatial 2R and a parallel three-degree-of-freedom manipulator.