Critical earthquake input power spectral density function models for engineering structures

The problem of determining optimal power spectral density models for earthquake excitation which satisfy constraints on total average power, zero crossing rate and which produce the highest response variance in a given linear system is considered. The solution to this problem is obtained using linear programming methods. The resulting solutions are shown to display a highly deterministic structure and, therefore, fail to capture the stochastic nature of the input. A modification to the definition of critical excitation is proposed which takes into account the entropy rate as a measure of uncertainty in the earthquake loads. The resulting problem is solved using calculus of variations and also within linear programming framework. Illustrative examples on specifying seismic inputs for a nuclear power plant and a tall earth dam are considered and the resulting solutions are shown to be realistic.

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.

Maximum weight independent sets in hole- and dart-free graphs

The Maximum Weight Independent Set (MWIS) problem on graphs with vertex weights asks for a set of pairwise nonadjacent vertices of maximum total weight. The complexity of the MWIS problem for hole-free graphs is unknown. In this paper, we first prove that the MWIS problem for (hole, dart, gem)-free graphs can be solved in O(n(3))-time. By using this result, we prove that the MWIS problem for (hole, dart)-free graphs can be solved in O(n(4))-time. Though the MWIS problem for (hole, dart, gem)-free graphs is used as a subroutine, we also give the best known time bound for the solvability of the MWIS problem in (hole, dart, gem)-free graphs. (C) 2012 Elsevier B.V. All rights reserved.

Microearthquake activity near the Idukki Reservoir, south India: A rare example of renewed triggered seismicity

Earthquakes triggered by artificial reservoirs have been documented for more than seven decades and the processes leading to this phenomenon are fairly well understood. Larger among such earthquakes are known to occur within a few years of reservoir impoundment and usually the activity decreases with time. A documented example of Reservoir Triggered Seismicity (RTS), the Idukki Reservoir in Kerala, south India, impounded in 1975, is an exception wherein the triggered activity has been revived in 2011, nearly 35 years after the initial burst of activity in 1977, two years after the dam was filled. The magnitude of the largest shock in the 2011 sequence exceeded that of the previously documented largest microearthquake. Presence of faults that are close to failure and vulnerable to increase in pore pressure due to reservoir loading or increased rainfall, or a combination of both seems to trigger shocks in this area. The renewed burst of earthquakes after a prolonged period of reduced activity at the Idukki Reservoir is a rare example of RTS. (C) 2012 Elsevier B.V. All rights reserved.

Assessment of hydrologic impacts of climate change in Tunga-Bhadra river basin, India with HEC-HMS and SDSM

Climate change would significantly affect many hydrologic systems, which in turn would affect the water availability, runoff, and the flow in rivers. This study evaluates the impacts of possible future climate change scenarios on the hydrology of the catchment area of the TungaBhadra River, upstream of the Tungabhadra dam. The Hydrologic Engineering Center's Hydrologic Modeling System version 3.4 (HEC-HMS 3.4) is used for the hydrological modelling of the study area. Linear-regression-based Statistical DownScaling Model version 4.2 (SDSM 4.2) is used to downscale the daily maximum and minimum temperature, and daily precipitation in the four sub-basins of the study area. The large-scale climate variables for the A2 and B2 scenarios obtained from the Hadley Centre Coupled Model version 3 are used. After model calibration and testing of the downscaling procedure, the hydrological model is run for the three future periods: 20112040, 20412070, and 20712099. The impacts of climate change on the basin hydrology are assessed by comparing the present and future streamflow and the evapotranspiration estimates. Results of the water balance study suggest increasing precipitation and runoff and decreasing actual evapotranspiration losses over the sub-basins in the study area.

The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India

We investigated area changes in glaciers covering an area of similar to 200 km(2) in the Tista basin, Sikkim, Eastern Indian Himalaya, between similar to 1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 +/- 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 +/- 0.82, 2.56 +/- 0.61 and 2.28 +/- 2.01 km(2) for the periods 1989-97, 1997-2004/05 and 2004-2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987-2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0 degrees C rise in the summer minimum temperature and a 2.0 degrees C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

A constant factor approximation algorithm for boxicity of circular arc graphs

The boxicity (resp. cubicity) of a graph G(V, E) is the minimum integer k such that G can be represented as the intersection graph of axis parallel boxes (resp. cubes) in R-k. Equivalently, it is the minimum number of interval graphs (resp. unit interval graphs) on the vertex set V, such that the intersection of their edge sets is E. The problem of computing boxicity (resp. cubicity) is known to be inapproximable, even for restricted graph classes like bipartite, co-bipartite and split graphs, within an O(n(1-epsilon))-factor for any epsilon > 0 in polynomial time, unless NP = ZPP. For any well known graph class of unbounded boxicity, there is no known approximation algorithm that gives n(1-epsilon)-factor approximation algorithm for computing boxicity in polynomial time, for any epsilon > 0. In this paper, we consider the problem of approximating the boxicity (cubicity) of circular arc graphs intersection graphs of arcs of a circle. Circular arc graphs are known to have unbounded boxicity, which could be as large as Omega(n). We give a (2 + 1/k) -factor (resp. (2 + log n]/k)-factor) polynomial time approximation algorithm for computing the boxicity (resp. cubicity) of any circular arc graph, where k >= 1 is the value of the optimum solution. For normal circular arc (NCA) graphs, with an NCA model given, this can be improved to an additive two approximation algorithm. The time complexity of the algorithms to approximately compute the boxicity (resp. cubicity) is O(mn + n(2)) in both these cases, and in O(mn + kn(2)) = O(n(3)) time we also get their corresponding box (resp. cube) representations, where n is the number of vertices of the graph and m is its number of edges. Our additive two approximation algorithm directly works for any proper circular arc graph, since their NCA models can be computed in polynomial time. (C) 2014 Elsevier B.V. All rights reserved.

Use of Data Driven Techniques for Short Lead Time Streamflow Forecasting in Mahanadi basin

Streamflow forecasts at daily time scale are necessary for effective management of water resources systems. Typical applications include flood control, water quality management, water supply to multiple stakeholders, hydropower and irrigation systems. Conventionally physically based conceptual models and data-driven models are used for forecasting streamflows. Conceptual models require detailed understanding of physical processes governing the system being modeled. Major constraints in developing effective conceptual models are sparse hydrometric gauge network and short historical records that limit our understanding of physical processes. On the other hand, data-driven models rely solely on previous hydrological and meteorological data without directly taking into account the underlying physical processes. Among various data driven models Auto Regressive Integrated Moving Average (ARIMA), Artificial Neural Networks (ANNs) are most widely used techniques. The present study assesses performance of ARIMA and ANNs methods in arriving at one-to seven-day ahead forecast of daily streamflows at Basantpur streamgauge site that is situated at upstream of Hirakud Dam in Mahanadi river basin, India. The ANNs considered include Feed-Forward back propagation Neural Network (FFNN) and Radial Basis Neural Network (RBNN). Daily streamflow forecasts at Basantpur site find use in management of water from Hirakud reservoir. (C) 2015 The Authors. Published by Elsevier B.V.

某水电站坝址区三维渗流参数反演

用有限元法对某水电站坝址区的渗流场参数进行了反演，结合工程地质和水文地质资料以及坝区的钻孔抽注水试验，反演得到了各覆盖层的渗透系数，并计算了整个坝址区的渗流场，为大坝建成后运营期的渗流场计算提供了基础。

Earthquake Stability Analysis about Fine Grinded Tailings Dams

Based on the broken characteristics about earthquake to tailings dams, the earthquake stability analysis methods for tailings dams are introduced. Taking fine tailings dam in Longdu Tailings Pool as an example, the stability of the dam with various situations while earthquake with seven magnitude takes place there. The results can be used by Longdu Mine for tailings pool safety management.

Influences of the Fish-Mouth Project and the Groins on the Flow and Sediment Ratio of the Yangtze River Waterway

A depth-integrated two-dimensional numerical model of current, salinity and sediment transport was proposed and calibrated by the observation data in the Yangtze River Estuary. It was then applied to investigate the flow and sediment ratio of the navigati

The MBA-LBA I period in the Kourion region : new evidences from Erimi- Laonin tou Porakou (Lemesos, Cyprus)

The evidence from Erimi-Laonin tou Porakou within the Kouris valley settlement system can offer further new data to outline a picture of the development trends from MBA to LBA I periods in the Kourion region. As to this area, such a “transitional” phase could be particularly relevant to better understand the development process of the mature LBA II–III urban centers (in Alassa as well as the most recent discovered in Erimi-Pitharka). In 2008–2009 a focus investigation (an intensive survey, geoperspections and excavations) within the Kouris Valley Project has been addressed to the site of Erimi-Laonin tou Porakou, which lies on one of the highest hilltops facing southward the Kouris dam. The excavations evidenced a double circuit wall which surrounds the top mound and the lower terrace. The top mound (Area A) has been widely occupied by a complex workshop, directly linked with a storage area. The natural limestone bedrock has been worked out to arrange a system of carved deep basins connected each other by a series of flow channels. Southward from the workshop, a large storage area has been cleared, with a relevant assemblage of big pithoi and storage devices. The first lower terrace area (Area B) has been possibly occupied by a domestic quarter; a cemetery area (Area E) has been also evidenced outside the external circuit wall, where a series of five rock-cut tombs (Tombs 228– 232) with single chamber and small incoming dromoi has been excavated.