Numerical modeling of two-phase hysteresis combined with an interface condition for heterogeneous porous media

This paper presents a numerical implementation of two-phase capillary hysteresis and its combination with a capillary interface condition for the treatment of heterogeneities. The hysteresis concepts chosen in this work are first implemented in a node-centered FV discretization scheme and subsequently combined with the interface condition that predicts sharp saturation discontinuities at material interfaces, based on a pressure equilibrium concept. This approach allows for the approximation of history-dependent, and at the same time discontinuous, saturations at material interfaces. The resulting model provides a well-defined evolution of the hysteretic capillary pressure–saturation relationships at material interfaces that is independent of the grid spacing. As demonstrated with a simple 1-D example, this concept therefore offers the advantage that the solution of a two-phase flow problem involving hysteresis does not relate to the grid resolution at the material interfaces.     

A simulation‐based framework for risk assessment and probabilistic sensitivity analysis of base‐isolated structures

A versatile, simulation‐based framework for risk assessment and probabilistic sensitivity analysis of base‐isolated structures is discussed in this work. A probabilistic foundation is used to address the various sources of uncertainties, either excitation or structural, and to characterize seismic risk. This risk is given, in this stochastic setting, by some statistics of the system response over the adopted probability models and stochastic simulation is implemented for its evaluation. An efficient, sampling‐based approach is also introduced for establishing a probabilistic sensitivity analysis to identify the importance of each of the uncertain model parameters in affecting the overall risk. This framework facilitates use of complex models for the structural system and the excitation. The adopted structural model explicitly addresses nonlinear characteristics of the isolators and of any supplemental dampers, and the effect of seismic pounding of the base to the surrounding retaining walls. An efficient stochastic ground motion model is also discussed for characterizing future near‐fault ground motions and relating them to the seismic hazard for the structural site. An illustrative example is presented that emphasizes the results from the novel probabilistic sensitivity analysis and their dependence on seismic pounding occurrences and on addition of supplemental dampers. Copyright © 2011 John Wiley & Sons, Ltd.     

Parametric analysis of single pile response in laterally spreading ground

Extensive damage to pile-supported structures has been witnessed in several recent earthquakes (Chi-Chi, 1999; Kobe, 1995, etc.), as a result of liquefaction-induced lateral spreading of slightly sloping ground or free-face topographic irregularities. This paper presents a parametric analysis of the basic pile and soil parameters, as well as the pile-soil interaction mechanisms affecting the response of single piles subjected to such lateral spreading, based on numerical simulation with the nonlinear P-y method. In parallel, a set of design charts and analytical relations is established, for approximate computation of maximum pile deflections and bending moments, using a “theory guided” multi-variable statistical analysis of the numerical predictions. Three different combinations (design cases) of pile head constraints and soil conditions were considered, which are commonly encountered in practice. The overall accuracy of the proposed analytical relations is evaluated against experimental results from seven centrifuge and five large shaking table experiments.     

Analyses of the warm season rainfall climatology of the northeastern US using regional climate model simulations and radar rainfall fields

We examine the warm season (April-September) rainfall climatology of the northeastern US through analyses of high-resolution radar rainfall fields from the Hydro-NEXRAD system and regional climate model simulations using the weather research and forecasting (WRF) model. Analyses center on the 5-year period from 2003 to 2007 and the study area includes the New York-New Jersey metropolitan region covered by radar rainfall fields from the Fort Dix, NJ WSR-88D. The objective of this study is to develop and test tools for examining rainfall climatology, with a special focus on heavy rainfall. An additional emphasis is on rainfall climatology in regions of complex terrain, like the northeastern US, which is characterized by land-water boundaries, large heterogeneity in land use and cover, and mountainous terrain in the western portion of the region. We develop a 5-year record of warm season radar rainfall fields for the study region using the Hydro-NEXRAD system. We perform regional downscaling simulations for the 5-year study period using the WRF model. Radar rainfall fields are used to characterize the interannual, seasonal and diurnal variation of rainfall over the study region and to examine spatial heterogeneity of rainfall. Regional climate model simulations are characterized by a wet bias in the rainfall fields, with the largest bias in the high-elevation regions of the model domain. We show that model simulations capture broad features of the interannual, seasonal, and diurnal variation of rainfall. Model simulations do not capture spatial gradients in radar rainfall fields around the New York metropolitan region and land-water boundaries to the east. The model climatology of convective available potential energy (CAPE) is used to interpret the regional distribution of warm season rainfall and the seasonal and diurnal variability of rainfall. We use hydrologic and meteorological observations from July 2007 to examine the interactions of land surface processes and rainfall from a regional perspective.     

Precursory seismic crustal deformation in the area of southern Albanides

On the basis of growing evidence thatstrong earthquakes are preceded by a periodof accelerating seismicity of moderatemagnitude earthquakes, an attempt is madeto search for such seismicity pattern in NWAegean area. Accelerating seismic crustaldeformation has been identified in the areaof southern Albanides mountain range(border region between Greece, formerYugoslavia and Albania). Based on certainproperties of this activity and on itssimilarity with accelerating seismicdeformation observed before a strongearthquake which occurred in the sameregion on 26 May 1960 (M = 6.5), we canconclude that a similar earthquake may begenerated in the same region during thenext few years. This conclusion is inagreement with independent results whichhave been derived on the basis of the timepredictable model.     

A review of inverse methods in seismic site characterization

Journal of Seismology - Seismic site characterization attempts to quantify seismic wave behavior at a specific location based on near-surface geophysical properties, for the purpose of mitigating...     

Geo-Characterization at selected accelerometric stations in Crete (Greece) and comparison of earthquake data recordings with EC8 elastic spectra

To estimate the seismic response according to Eurocode (EC8) and almost all other national codes, site conditions have to be properly characterized so that soil amplification and the corresponding peak ground motion can be calculated. In this work, different geophysical and geotechnical methods are combined in order to define the detailed ground conditions in selected sites of the Hellenic Accelerometric Network (HAN) in Crete. For this purpose, the geological information of the sites and shear wave velocity, calculated from surface wave measurements, is used. Additionally, ground acceleration data recorded through HAN have been utilized from intermediate depth earthquakes in the broader area of South Aegean Sea. Using the recorded ground motion data and the procedure defined in EC8, the corresponding elastic response spectrum is calculated for the selected sites. The resulting information is compared to the values defined in the corresponding EC8 spectrum for the seismic zone that includes the island of Crete. The comparison shows that accurate definition of ground type through geological, geotechnical and geophysical investigations is important. However, our current comparison focuses on the distribution of values rather than the absolute values of EC8-prescribed spectra, and the results should be considered in this context.     

Supported Axisymmetric Tunnels Within Linear Viscoelastic Burgers Rocks

An exact closed form solution is derived for the mechanical behaviour of a linear viscoelastic Burgers rock around an axisymmetric tunnel, supported by a linear elastic ring. Analytical formulae are provided for the displacement of the rock/lining interface and for the pressure exerted by the rock on the lining, taking into account the stiffness and its installation time. Results calculated from these formulae do validate the corresponding numerical results of a 2D finite differences code. Further, comparison to previous existing solutions for the same viscoelastic model indicates similarities and differences. A parametric study is performed to investigate the effect of the viscoelastic constants, the stiffness and installation time of the support. The derived closed form solution is used to construct the time-dependent Supported Ground Reaction Curves of the viscoelastic rock, i.e. the time contour plots on the convergence confinement diagram. The importance of the effect of the support on the restrained rock creep and the exerted pressure on the lining, during the design life of a structure, is examined.