Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

Effect of pouring temperature on cooling slope casting of semi-solid Al-Si-Mg alloy

Present trend of semi-solid processing is directed towards rheocasting route which allows manufacturing of near-net-shape cast components directly from the prepared semi-solid slurry. Generation of globular equi-axed grains during solidification of rheocast components, compared to the columnar dendritic structure of conventional casting routes, facilitates the manufacturing of components with improved mechanical properties and structural integrity. In the present investigation, a cooling slope has been designed and indigenously fabricated to produce semi solid slurry of Al-Si-Mg (A356) alloy and successively cast in a metallic mould. The scope of the present work discusses about development of a numerical model to simulate the liquid metal flow through cooling slope using Eulerian two-phase flow approach and to investigate the effect of pouring temperature on cooling slope semi-solid slurry generation process. The two phases considered in the present model are liquid metal and air. Solid fraction evolution of the solidifying melt is tracked at different locations of the cooling slope, following Schiel's equation. The continuity equation, momentum equation and energy equation are solved considering thin wall boundary condition approach. During solidification of the liquid metal, a modified temperature recovery scheme has been employed taking care of the latent heat release and change of fraction of liquid. The results obtained from simulations are compared with experimental findings and good agreement has been found.

Characterization of groundwater chemistry under the influence of lithologic and anthropogenic factors along a climatic gradient in Upper Cauvery basin, South India

Hydrogeological and climatic effect on chemical behavior of groundwater along a climatic gradient is studied along a river basin. `Semi-arid' (500-800 mm of mean annual rainfall), `sub-humid' (800-1,200 mm/year) and `humid' (1,200-1,500 mm/year) are the climatic zones chosen along the granito-gneissic plains of Kabini basin in South India for the present analysis. Data on groundwater chemistry is initially checked for its quality using NICB ratio (<+/- 5 %), EC versus TZ+ (similar to 0.85 correlation), EC versus TDS and EC versus TH analysis. Groundwater in the three climatic zones is `hard' to `very hard' in terms of Ca-Mg hardness. Polluted wells are identified (> 40 % of pollution) and eliminated for the characterization. Piper's diagram with mean concentrations indicates the evolution of CaNaHCO3 (semi-arid) from CaHCO3 (humid zone) along the climatic gradient. Carbonates dominate other anions and strong acids exceeded weak acids in the region. Mule Hole SEW, an experimental watershed in sub-humid zone, is characterized initially using hydrogeochemistry and is observed to be a replica of entire sub-humid zone (with 25 wells). Extension of the studies for the entire basin (120 wells) showed a chemical gradient along the climatic gradient with sub-humid zone bridging semi-arid and humid zones. Ca/Na molar ratio varies by more than 100 times from semi-arid to humid zones. Semi-arid zone is more silicaceous than sub-humid while humid zone is more carbonaceous (Ca/Cl similar to 14). Along the climatic gradient, groundwater is undersaturated (humid), saturated (sub-humid) and slightly supersaturated (semi-arid) with calcite and dolomite. Concentration-depth profiles are in support of the geological stratification i.e., not approximate to 18 m of saprolite and similar to 25 m of fracture rock with parent gneiss beneath. All the wells are classified into four groups based on groundwater fluctuations and further into `deep' and `shallow' based on the depth to groundwater. Higher the fluctuations, larger is its impact on groundwater chemistry. Actual seasonal patterns are identified using `recharge-discharge' concept based on rainfall intensity instead of traditional monsoon-non-monsoon concept. Non-pumped wells have low Na/Cl and Ca/Cl ratios in recharge period than in discharge period (Dilution). Few other wells, which are subjected to pumping, still exhibit dilution chemistry though water level fluctuations are high due to annual recharge. Other wells which do not receive sufficient rainfall and are constantly pumped showed high concentrations in recharge period rather than in discharge period (Anti-dilution). In summary, recharge-discharge concept demarcates the pumped wells from natural deep wells thus, characterizing the basin.

Semi-analytical finite element analysis of a strip footing on an elastic reinforced soil

A plane strain elastic interaction analysis of a strip footing resting on a reinforced soil bed has been made by using a combined analytical and finite element method (FEM). In this approach the stiffness matrix for the footing has been obtained using the FEM, For the reinforced soil bed (halfplane) the stiffness matrix has been obtained using an analytical solution. For the latter, the reinforced zone has been idealised as (i) an equivalent orthotropic infinite strip (composite approach) and (ii) a multilayered system (discrete approach). In the analysis, the interface between the strip footing and reinforced halfplane has been assumed as (i) frictionless and (ii) fully bonded. The contact pressure distribution and the settlement reduction have been given for different depths of footing and scheme of reinforcement in soil. The load-deformation behaviour of the reinforced soil obtained using the above modelling has been compared with some available analytical and model test results. The equivalent orthotropic approach proposed in this paper is easy to program and is shown to predict the reinforcing effects reasonably well.

Rheology of particle-loaded semi-dilute polymer solutions

We present measurements of the rheology of suspensions of rigid spheres in a semi-dilute polymer solution from experiments of steady and oscillatory shear. For a given value of the shear rate gamma, addition of particles enhances the viscosity and the first normal stress difference but decreases the magnitude of the second normal stress difference. The viscosity eta exhibits a power law variation in gamma for a range of gamma that grows with phi. The first normal stress N-1 is positive and its value grows with phi; it exhibits a clear power law variation for the entire range of gamma that was studied. The second normal stress difference N-2 is negative for the pure polymer solution and much smaller in magnitude than N-1; on addition of particles, its magnitude further decreases, and it appears to change sign at large phi. The behavior of N-1 and N-2 is at odds with the findings of recent studies on particle-loaded dilute polymer solutions and polymer melts. The small-amplitude oscillatory shear experiments show the linear viscoelastic properties, G(') and G('), increasing with phi at a given value of the angular frequency omega. The dynamic viscosity of the suspension differs substantially from its steady shear viscosity, and the difference increases as gamma, omega -> 0.

An upper bound for Cubicity in terms of Boxicity

An axis-parallel b-dimensional box is a Cartesian product R-1 x R-2 x ... x R-b where each R-i (for 1 <= i <= b) is a closed interval of the form [a(i), b(i)] on the real line. The boxicity of any graph G, box(G) is the minimum positive integer b such that G can be represented as the intersection graph of axis-parallel b-dimensional boxes. A b-dimensional cube is a Cartesian product R-1 x R-2 x ... x R-b, where each R-i (for 1 <= i <= b) is a closed interval of the form [a(i), a(i) + 1] on the real line. When the boxes are restricted to be axis-parallel cubes in b-dimension, the minimum dimension b required to represent the graph is called the cubicity of the graph (denoted by cub(G)). In this paper we prove that cub(G) <= inverted right perpendicularlog(2) ninverted left perpendicular box(G), where n is the number of vertices in the graph. We also show that this upper bound is tight.Some immediate consequences of the above result are listed below: 1. Planar graphs have cubicity at most 3inverted right perpendicularlog(2) ninvereted left perpendicular.2. Outer planar graphs have cubicity at most 2inverted right perpendicularlog(2) ninverted left perpendicular.3. Any graph of treewidth tw has cubicity at most (tw + 2) inverted right perpendicularlog(2) ninverted left perpendicular. Thus, chordal graphs have cubicity at most (omega + 1) inverted right erpendicularlog(2) ninverted left perpendicular and circular arc graphs have cubicity at most (2 omega + 1)inverted right perpendicularlog(2) ninverted left perpendicular, where omega is the clique number.

Incompressible micropolar fluid flow over a semi-infinite plate

Micropolar fluid flow over a semi-infinite flat plate has been described by using the parabolic co-ordinates and the method of series truncation in order to study the flow for low to large Reynolds numbers. These co-ordinates permit to study the flow regime at the leading edge. Numerical results have been presented for different Reynolds numbers. Results show a reduction in skin friction.

Unsteady flow of a dilute suspension in a semi-infinite contracting or expanding pipe

In the present paper an exact similar solution of the Navier-Stokes equation for unsteady flow of a dilute suspension in a semi-infinite contracting or expanding circular pipe is presented. The effects of the Schmidt number (Sc), Reynolds number (|ε|), the volume fraction (α) and the relaxation time (τ) of the particulate phase on the flow characteristics are examined. The presence of the solid particles has been observed to influence the flow behaviour significantly. These solutions are valid down to the state of a completely collapsed pipe, since the nonlinearity is retained fully. The results may help understanding the flow near the heart and certain forced contractions or expansions of valved veins.

Semi-similar solution of an unsteady compressible three-dimensional stagnation point boundary layer flow with massive blowing

A semi-similar solution of an unsteady laminar compressible three-dimensional stagnation point boundary layer flow with massive blowing has been obtained when the free stream velocity varies arbitrarily with time. The resulting partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme with a quasi-linearization technique in the nodal point region and an implicit finite-difference scheme with a parametric differentiation technique in the saddle point region. The results have been obtained for two particular unsteady free stream velocity distributions: (i) an accelerating stream and (ii) a fluctuating stream. Results show that the skin-friction and heat-transfer parameters respond significantly to the time dependent arbitrary free stream velocity. Velocity and enthalpy profiles approach their free stream values faster as time increases. There is a reverse flow in the y-wise velocity profile, and overshoot in the x-wise velocity and enthalpy profiles in the saddle point region, which increase as injection and wall temperature increase. Location of the dividing streamline increases as injection increases, but as the wall temperature and time increase, it decreases.

Vertical uplift capacity of two interfering horizontal anchors in sand using an upper bound limit analysis

The vertical uplift resistance of two interfering rigid rough strip anchors embedded horizontally in sand at shallow depths has been examined. The analysis is performed by using an upper bound theorem o limit analysis in combination with finite elements and linear programming. It is specified that both the anchors are loaded to failure simultaneously at the same magnitude of the failure load. For different clear spacing (S) between the anchors, the magnitude of the efficiency factor (xi(gamma)) is determined. On account of interference, the magnitude of xi(gamma) is found to reduce continuously with a decrease in the spacing between the anchors. The results from the numerical analysis were found to compare reasonably well with the available theoretical data from the literature.

Heat transfer with solid liquid phase change initiating at a point in a semi infinite mold

Short-time analytical solutions of temperature and moving boundary in two-dimensional two-phase freezing due to a cold spot are presented in this paper. The melt occupies a semi-infinite region. Although the method of solution is valid for various other types of boundary conditions, the results in this paper are given only for the prescribed flux boundary conditions which could be space and time dependent. The freezing front propagations along the interior of the melt region exhibit well known behaviours but the propagations along the surface are of new type. The freezing front always depends on material parameters. Several interesting results can be obtained as particular cases of the general results.

Risk Minimizing Option Pricing in a Semi-Markov Modulated Market

We address risk minimizing option pricing in a semi-Markov modulated market where the floating interest rate depends on a finite state semi-Markov process. The growth rate and the volatility of the stock also depend on the semi-Markov process. Using the Föllmer–Schweizer decomposition we find the locally risk minimizing price for European options and the corresponding hedging strategy. We develop suitable numerical methods for computing option prices.

Portfolio Optimization in a Semi-Markov Modulated Market

We address a portfolio optimization problem in a semi-Markov modulated market. We study both the terminal expected utility optimization on finite time horizon and the risk-sensitive portfolio optimization on finite and infinite time horizon. We obtain optimal portfolios in relevant cases. A numerical procedure is also developed to compute the optimal expected terminal utility for finite horizon problem.