Computationally efficient algorithms for the two-dimensional Kolmogorov-Smirnov test

Goodness-of-fit statistics measure the compatibility of random samples against some theoretical or reference probability distribution function. The classical one-dimensional Kolmogorov-Smirnov test is a non-parametric statistic for comparing two empirical distributions which defines the largest absolute difference between the two cumulative distribution functions as a measure of disagreement. Adapting this test to more than one dimension is a challenge because there are 2^d-1 independent ways of ordering a cumulative distribution function in d dimensions. We discuss Peacock's version of the Kolmogorov-Smirnov test for two-dimensional data sets which computes the differences between cumulative distribution functions in 4n^2 quadrants. We also examine Fasano and Franceschini's variation of Peacock's test, Cooke's algorithm for Peacock's test, and ROOT's version of the two-dimensional Kolmogorov-Smirnov test. We establish a lower-bound limit on the work for computing Peacock's test of Omega(n^2.lg(n)), introducing optimal algorithms for both this and Fasano and Franceschini's test, and show that Cooke's algorithm is not a faithful implementation of Peacock's test. We also discuss and evaluate parallel algorithms for Peacock's test

A novel method for predicting the response variability of friction-damped gas turbine blades

Predicting the response of gas turbine blades with underplatform friction dampers is challenging due to the combination of frictional nonlinearity and system uncertainty: a traditional Monte Carlo approach to predicting response distributions requires a large number of nonlinear simulations which is computationally expensive. This paper presents a new approach based on the principle of Maximum Entropy that provides an estimate of the response distribution that is approximately two orders of magnitude faster than Monte Carlo Harmonic Balance Method simulations. The premise is to include the concept of `computational uncertainty': incorporating lack of knowledge of the solution as part of the uncertainty, on the basis that there are diminishing returns in computing precise solutions to an uncertain system. To achieve this, the method uses a describing function approximation of the friction-damped part of the system; chooses an ignorance prior probability density function for the complex value of the describing function based on Coulombs friction law; updates the distribution using an estimate of the mean solution, the admissible domain of solutions, and the principle of Maximum Entropy; then carries out a linear Monte Carlo simulation to estimate the response distribution. The approach is validated by comparison with HBM simulations and experimental tests, using an idealised academic system consisting of a periodic array of beams (with controllable uncertainty) coupled by single-point friction dampers. Comparisons with two- and eight-blade systems show generally good agreement. Predicting the response statistics of the maximum blade amplitude reveals specific well-understood circumstances when the method is less effective. Predictions of the overall blade response statistics agree with Monte Carlo HBM extremely well across a wide range of excitation amplitudes and uncertainty levels. Critically, experimental comparisons reveal the care that is needed in accurately characterising uncertainty in order to obtain agreement of response percentiles. The new method allowed fast iteration of uncertainty parameters and correlations to achieve good agreement, which would not have been possible using traditional methods.Mitsubishi Heavy Industrie

Predicting response bounds for friction-damped gas turbine blades with uncertain friction coupling

Friction dampers are often used to reduce high amplitude vibration within gas turbines: they are a robust solution that are able to withstand extreme operating environments. Although the turbine blades are manufactured to tight tolerances, there can be significant variability in the overall response of the assembly. Uncertainties associated with the frictional contact properties are a major factor contributing to this variability. This paper applies a recently developed method for predicting response bounds to friction-damped gas turbines when the characteristics of the friction dampers are unknown, including uncertainty regarding the underlying functional form of the friction law. The approach taken is to represent the frictional contact using a describing function, and formulate an optimisation problem to seek upper and lower bounds on a chosen response metric, such as displacement amplitude. Constraints are chosen that describe known properties of the frictional nonlinearity, without needing to specify a particular constitutive law. The method was validated by comparison with numerical and experimental results from an idealised test system. The experimental test rig consisted of an array of eight beams coupled by pin-contact friction dampers. A modal description of this test rig formed the basis of a numerical model, which uses the Harmonic Balance Method (HBM) for nonlinear simulations. A set of Monte Carlo tests was carried out numerically and experimentally for both a two-beam sub-assembly as well as for the full eight-beam assembly. Comparisons with numerical results showed excellent agreement providing confident verification of the implementation, and comparisons with experimental results revealed that the bounds became less conservative as the system complexity increased. Overall the results are promising: upper and lower response bounds for an array of friction-damped systems can be computed at similar cost to a single HBM simulation, giving reliable bounds that are valid for both parametric and model uncertainties associated with the friction couplings.Mitsubishi Heavy Industrie

The enhancement of ferromagnetism in uniaxially stressed diluted magnetic semiconductors

We predict a new mechanism of enhancement of ferromagnetic phase transition temperature $T_c$ in uniaxially stressed diluted magnetic semiconductors (DMS) of p-type. Our prediction is based on comparative studies of both Heisenberg (inherent to undistorted DMS with cubic lattice) and Ising (which can be applied to strongly enough stressed DMS) models in a random field approximation permitting to take into account the spatial inhomogeneity of spin-spin interaction. Our calculations of phase diagrams show that area of parameters for existence of DMS-ferromagnetism in Ising model is much larger than that in Heisenberg model.Comment: Accepted for publication in Phys. Rev.

Fractalkine: A Cellular Link Between Adipose Tissue Inflammation and Vascular Pathologies

It is hard to imagine, given the wealth of new datareported over the recent past, that adipose tissue atone time was primarily considered as a passive res-ervoir for energy deposition and storage. However, research beginning in the early 1990s on the role of tumor necrosis factor (TNF)-a ushered in a new era of inves-tigation, and since that time, there has been an incredibly rapid and substantive increase in our understanding of underlying physiologic systems and molecular pathways linking obesity, inflammation, and insulin action (1,2). Spe-cifically, our understanding of the link between obesity and carbohydrate metabolism has been significantly enhanced with the elucidation of key regulators of energy balance and cellular insulin signaling that are complex and highly in-tegrated (3–6). We now readily accept adipose tissue as a key endocrine organ regulating processes throughout the body with its significant number of adipocyte secre-tions. What now appears to be emerging is the elucidatio

Sharp lines in the absorption edge of EuTe and Pb0.1_{0.1}Eu0.9_{0.9}Te in high magnetic fields

The optical absorption spectra in the region of the \fd transition energies of epitaxial layers of of EuTe and \PbEuTe, grown by molecular beam epitaxy, were studied using circularly polarized light, in the Faraday configuration. Under \sigmam polarization a sharp symmetric absorption line (full width at half-maximum 0.041 eV) emerges at the low energy side of the band-edge absorption, for magnetic fields intensities greater than 6 T. The absorption line shows a huge red shift (35 meV/T) with increasing magnetic fields. The peak position of the absorption line as a function of magnetic field is dominated by the {\em d-f} exchange interaction of the excited electron and the \Euion spins in the lattice. The {\em d-f} exchange interaction energy was estimated to be $J_{df}S=0.15\pm 0.01$ eV. In \PbEuTe the same absorption line is detected, but it is broader, due to alloy disorder, indicating that the excitation is localized within a finite radius. From a comparison of the absorption spectra in EuTe and \PbEuTe the characteristic radius of the excitation is estimated to be $\sim 10$\AA.Comment: Journal of Physics: Condensed Matter (2004, at press

The effect of maximum normal impact load, absorbed energy and contact impulse on the impact craters volume/depth of DLC coating under repetitive impacts

Recently, the requirements for measuring dynamic responses have become severe and varied in many industrial and research applications such as material testing, model analysis and crash testing1). Surface degradation often occurs due to this dynamic response. This phenomenon also appears in the DLC coatings material. In this present work, a self-developed horizontal impact tester can provides this type of response, which gives an effect to the impact craters volume/depth of DLC coating

Raman Spectroscopy Study of Impacted DLC Coatings

Raman scattering is an excellent tool to characterize the structure of carbon atoms in diamond-like carbon (DLC). The study of DLC coatings in the sliding conditions has been conducted for a decade using Raman spectroscopy analysis. However, there is still insufficient information about how the structure of DLC coatings changes during repetitive impact. In this paper changes in the structure under 90o repetitive impact at a large number of cycles are presented