Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Micro-fracture instabilities in granular solids

We study the aspect of unstable behavior (like strain localization bands) in elastic solids as a consequence of micro-fracturing. A two-scale approach of computational homogenization is considered. The macroscopic behavior is investigated by finite element computations on a unit cell. At the micro-level, we consider a granular structure with elastic grains. The inter-granular boundaries are modeled with cohesive laws, friction and unilateral contact. We show that decohesion between grains gives rise to macro-instabilities, indicated by the loss of ellipticity, typical for deformation localization bands. The relation between the microscopic softening on inter-granular boundaries and the onset of macro-instabilities is studied through numerical examples. The influence of the cohesive law and friction parameters is analyzed. For periodic distributions of granular structures, we prove the loss of periodicity by failure and the corresponding size dependence effect in the homogenized response. We present numerical examples of bifurcation of solutions for granular cell structures and of particular solutions specific to elementary volumes with periodic cell distribution. Size dependence appears in the unstable regime and is strongly influenced by cohesion and friction parameters.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

Effects of onshore and offshore topography on sea breeze circulation: An observational study at eastern Attica,Greece

Field measurements of wind, air temperature and humidity were taken at the eastern part of the Attika district in June 1991, to examine the topographic influences exerted on the local sea breezes. These influences are due either to the elongated Evia island, faced by the northern half of Attica coastline some tens of kilometers offshore, or to the coast-parallel range of Hymettos mountain, rising steeply 12 km onshore. The instrumentation consisted mainly of three tethered meteorological balloons released at characteristic sires (i.e., the coast, a location between shoreline and mountain foot and the mountain top) and three autographic ground-based anemometers operating at selected locations. Data from the ground-based and upper air stations of the Hellenic National Meteorological Service, as well as the diurnal weather maps were also obtained and analyzed. Observations were made under different synoptic wind and the latter was found to determine remarkably the significance of the topographic effects. A preliminary two-dimensional numerical approach was also made concerning the sea breeze capability to reach the Hymettos mountain top in the case of a weak opposing geostrophic flow.     

The “Zero-order Model” revisited: Paul Schindler's influence on the development of trace metal scavenging models

In 1975 Paul Schindler produced the first oceanic trace metal scavenging model to explicitly include the role of surface chemistry as a control on trace metal water column residence times. The eighteen years that have elapsed since the publication of Schindler's seminal paper have seen the development of a variety of oceanic scavenging models; yet, the fundamental insight of his Zero-order Model remains the benchmark. This paper describes the role of Paul Schindler's work on surface chemistry in providing a framework for the current generation of trace element scavenging models.     

首都减灾圈的思考与建议

首都减灾圈,系首都减轻自然灾害预测防治圈。1991年12月20—21日在北京召开了《首都圈自然灾害及其减灾对策研讨会》。本文根据这次会议所提供的材料,在从整体上实现减灾的思想指导下,就首都减灾圈的“成灾背景”、“首都减灾圈的组成”、“首都减灾圈的灾害预测与防治状况”、“首都减灾圈的灾害关联性分析”、“首都减灾圈的减灾实效预估”和“首都减灾圈的减灾对策与实施”等六个带有共同性的问题,进行了讨论。供制定首都圈减灾方案时参考。     

Estimation of runoff distribution within river systems: Application to the Lena Basin (Siberia)

A channel account approach is proposed to estimate longitudinal changes in runoff along large river systems. This new method provides a quantitative basis for describing the fluvial transport of suspended particulate material and dissolved substances. This method includes an evaluation of basic elements of water balance in separate sections of the river network and subsequent correction of channel accounting equations for the entire system using a maximum likelihood principle. To calculate water discharges of tributaries that have no hydrological information, structural analysis of river network is performed. This approach provides less error in comparison with traditional methods of estimating lateral inflow. The method is used to trace water discharge with increasing distance along the Lena river basin and to evaluate the contribution of geologically and lithologically uneven sub-basins in water discharge formation during a summer low water period.     

Site Effects Study in the Serravalle Village (NE Italy) by Means of Spectral Ratios Analyses Relatively to Multiple Reference Sites

The area of Serravalle, sited in the northern part of the town of Vittorio Veneto (TV), NE Italy, has been the target of a seismic microzonation campaign. 10 seismic stations have been deployed for a 7 months period to record in continuous mode. Three stations were installed on bedrock outcrops and seven on sedimentary sites with variable cover thickness. Spectral analyses have been performed on the collected data-set using the Generalized Inversion Technique (GIT, e.g. Andrews, 1986). In particular, spectral ratios have been calculated for each station relatively to the average of the three reference, bedrock sites. The spectral ratios provide quantitative estimates of the seismic motion amplifications which occur in each of the monitored sites. Two sites show high values of amplification, 5 times larger than signal amplitude at the reference sites, in correspondence of well discernible peak frequencies of 5 Hz. Results for the other stations show smaller amounts of site amplification spreading over a broad range of frequencies. Sites where the highest amplifications were recorded all lie on the left bank of the Meschio River and in areas farther away from its outlet into the plain correlating with the presence of thick layers of Quaternary deposits.     

Temperature dependence of large-scale water-retention curves: a case study

A local-scale model for temperature-dependence of water-retention curves may be applicable to large scales. Consideration of this temperature dependence is important for modeling unsaturated flow and transport in the subsurface in numerous cases. Although significant progress has been made in understanding and modeling this temperature effect, almost all the previous studies have been limited to small scales (on the order of several centimeters). Numerical experiments were used to investigate the possibility of extending a local-scale model for the temperature-dependence of water retention curves to large scales (on the order of meters). Temperature effects on large-scale hydraulic properties are of interest in many practical applications. Numerical experiment results indicate that the local-scale model can indeed be applicable to large-scale problems for special porous media with high air entry values. A typical porous medium of this kind is the porous tuff matrix in the unsaturated zone of Yucca Mountain, Nevada, the proposed geologic disposal site for national high-level nuclear wastes. Whether this finding can approximately hold for general cases needs to be investigated in future studies.