Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

On the orbits and masses of the satellites of the Pluto-Charon system

Two small satellites of Pluto, S/2005 P1 (hereafter P1) and S/2005 P2 (hereafter P2), have recently been discovered outside the orbit of Charon, and their orbits are nearly circular and nearly coplanar with that of Charon. Because the mass ratio of Charon-Pluto is ∼0.1, the orbits of P2 and P1 are significantly non-Keplerian even if P2 and P1 have negligible masses. We present an analytic theory, with P2 and P1 treated as test particles, which shows that the motion can be represented by the superposition of the circular motion of a guiding center, the forced oscillations due to the non-axisymmetric components of the potential rotating at the mean motion of Pluto-Charon, the epicyclic motion, and the vertical motion. The analytic theory shows that the azimuthal periods of P2 and P1 are shorter than the Keplerian orbital periods, and this deviation from Kepler's third law is already detected in the unperturbed Keplerian fit of Buie and coworkers. In this analytic theory, the periapse and ascending node of each of the small satellites precess at nearly equal rates in opposite directions. From direct numerical orbit integrations, we show the increasing influence of the proximity of P2 and P1 to the 3:2 mean-motion commensurability on their orbital motion as their masses increase within the ranges allowed by the albedo uncertainties. If the geometric albedos of P2 and P1 are high and of order of that of Charon, the masses of P2 and P1 are sufficiently low that their orbits are well described by the analytic theory. The variation in the orbital radius of P2 due to the forced oscillations is comparable in magnitude to that due to the best-fit Keplerian eccentricity, and there is at present no evidence that P2 has any significant epicyclic eccentricity. However, the orbit of P1 has a significant epicyclic eccentricity, and the prograde precession of its longitude of periapse with a period of 5300 days should be easily detectable. If the albedos of P2 and P1 are as low as that of comets, the large inferred masses induce significant short-term variations in the epicyclic eccentricities and/or periapse longitudes on the 400-500-day timescales due to the proximity to the 3:2 commensurability. In fact, for the maximum inferred masses, P2 and P1 may be in the 3:2 mean-motion resonance, with the resonance variable involving the periapse longitude of P1 librating. Observations that sample the orbits of P2 and P1 well on the 400-500-day timescales should provide strong constraints on the masses of P2 and P1 in the near future.     

OSIRIS: AO-assisted integral-field spectroscopy at the Keck Observatory

OSIRIS (OH-Suppressing Infra-Red Integral-field Spectrograph) is a new facility instrument for the Keck Observatory. After seeing first light in February 2005, OSIRIS is currently undergoing commissioning. OSIRIS provides the capability of performing three-dimensional spectroscopy in the near-infrared z, J, H, and K bands at the resolution limit of the Keck II telescope, which is equipped with adaptive optics and a laser guide star. The science case for OSIRIS is summarized, and the instrument and associated data reduction software are described.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

有关边界层问题中一类奇异积分方程新的存在性结果

对边界层理论新结果中出现的一类奇异积分方程w（t）=∫1 t（1-s）（λ＋λs＋s）/w（s）ds＋（1-t）∫t 0s/w（s）ds,t∈（0,1）进行讨论，并得出了上述方程在λ∈（-1/2,0）上正解存在性的新结果。     

Morphological redshift estimates for galaxy clusters in a Sunyaev–Zel'dovich effect survey

We develop a new method to estimate the redshift of galaxy clusters through resolved images of the Sunyaev–Zel'dovich effect (SZE). Our method is based on morphological observables which can be measured by actual and future SZE experiments. We test the method with a set of high-resolution hydrodynamical simulations of galaxy clusters at different redshifts. Our method combines the observables in a principal component analysis. After calibrating the method with an independent redshift estimation for some of the clusters, we show – using a Bayesian approach – how the method can give an estimate of the redshift of the galaxy clusters. Although the error bars given by the morphological redshift estimation are large, it should be useful for future SZE surveys where thousands of clusters are expected to be detected; a first preselection of the high-redshift candidates could be done using our proposed morphological redshift estimator. Although not considered in this work, our method should also be useful to give an estimate of the redshift of clusters in X-ray and optical surveys.     

Size-Flux Relation in Solar Active Regions

We present a study of the relationship between integral area and corresponding total magnetic flux for solar active regions. It is shown that some of these relationships are satisfied to simple power laws. Fractal examination showed that some of these power laws can not be justified inside the simple models of stationary magnetic flux tube aggregation. All magnetic fluxes and corresponding areas were calculated using the data measured with the Solar Magnetic Field Telescope of the Huairou Solar Observing Station in Beijing.     

Wavelet Analysis of Space Solar Telescope Images

The scientific satellite SST (Space Solar Telescope) is an important research project strongly supported by the Chinese Academy of Sciences. Every day, SST acquires 50 GB of data (after processing) but only 10GB can be transmitted to the ground because of limited time of satellite passage and limited channel volume. Therefore, the data must be compressed before transmission. Wavelets analysis is a new technique developed over the last 10 years, with great potential of application. We start with a brief introduction to the essential principles of wavelet analysis, and then describe the main idea of embedded zerotree wavelet coding, used for compressing the SST images. The results show that this coding is adequate for the job.