Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Radiative Transfer in Inhomogeneous Atmospheres. III

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.     

On cross-correlating weak lensing surveys

The present generation of weak lensing surveys will be superseded by surveys run from space with much better sky coverage and high level of signal-to-noise ratio, such as the Supernova/Acceleration Probe ( SNAP ). However, removal of any systematics or noise will remain a major cause of concern for any weak lensing survey. One of the best ways of spotting any undetected source of systematic noise is to compare surveys that probe the same part of the sky. In this paper we study various measures that are useful in cross-correlating weak lensing surveys with diverse survey strategies. Using two different statistics – the shear components and the aperture mass – we construct a class of estimators which encode such cross-correlations. These techniques will also be useful in studies where the entire source population from a specific survey can be divided into various redshift bins to study cross-correlations among them. We perform a detailed study of the angular size dependence and redshift dependence of these observables and of their sensitivity to the background cosmology. We find that one-point and two-point statistics provide complementary tools which allow one to constrain cosmological parameters and to obtain a simple estimate of the noise of the survey.     

Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz and Klu’zniak (2001). In a first paper (P'etri, 2005a, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. After a discussion on the magnitude of this deformation applied to neutron stars, we show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to {three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric resonance}. In a second part, we focus on the linear response of a thin accretion disk in the 2D limit. {Waves are launched at the aforementioned resonance positions and propagate in some permitted regions inside the disk, according to the dispersion relation obtained by a WKB analysis}. In a last part, these results are confirmed and extended via non linear hydrodynamical numerical simulations performed with a pseudo-spectral code solving Euler's equations in a 2D cylindrical coordinate frame. {We found that for a weak potential perturbation, the Lindblad resonance is the only effective mechanism producing a significant density fluctuation}. In a last step, we replaced the Newtonian potential by the so called logarithmically modified pseudo-Newtonian potential in order to take into account some general-relativistic effects like the innermost stable circular orbit (ISCO). The latter potential is better suited to describe the close vicinity of a neutron star or a black hole. However, from a qualitative point of view, the resonance conditions remain the same. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Discrimination of sedimentary facies by association analysis

The ability of association analysis to discriminate sedimentary facies was tested on Purdy's modal analyses of modern sediments of the Great Bahama Bank. Purdy's data set has served in the past as a standard reference for evaluating various multivariate classification algorithms. In order to adapt Purdy's data to association analysis, the percent abundance of the 12 constituents was converted to binary form by dichotomizing each variable on its mean value. The results obtained by association analysis are virtually identical to those obtained by Purdy and other authors. The same four main sedimentary facies were discriminated; 86% of the samples were identically classified (97% when misclassified borderline cases are counted as matches); the total partition variance of the classification is only negligibly greater (4%); and the grouping of the variables yielded the same four groups. The rank order of the three division-attributes responsible for the sample classification is fines, oolites, and corals. Association analysis has been employed by other authors to differentiate meaningful facies groups in studies of ancient reef carbonates, modern reef sediments, and heavy minerals in stream sediments. In all these studies, the results were found to be compatible with those obtained by using the continuous quantitative measurements, indicating that qualitative binary data may often be sufficient for the purpose of facies discrimination in many branches of geology and that association analysis is an effective method for this purpose.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.