AVES: A high performance computer cluster array for the INTEGRAL satellite scientific data analysis

In this paper we describe a new computing system array, designed, built and now used at the Space Astrophysics and Planetary Institute (IAPS) in Rome, Italy, for the INTEGRAL Space Observatory scientific data analysis. This new system has become necessary in order to reduce the processing time of the INTEGRAL data accumulated during the more than 9 years of in-orbit operation. In order to fulfill the scientific data analysis requirements with a moderately limited investment the starting approach has been to use a ??cluster?? array of commercial quad-CPU computers, featuring the extremely large scientific and calibration data archive on line.     

Normal forms for the epicyclic approximations of the perturbed Kepler problem

We compute the normal forms for the Hamiltonian leading to the epicyclic approximations of the (perturbed) Kepler problem in the plane. The Hamiltonian setting corresponds to the dynamics in the Hill synodic system where, by means of the tidal expansion of the potential, the equations of motion take the form of perturbed harmonic oscillators in a rotating frame. In the unperturbed, purely Keplerian case, the post-epicyclic solutions produced with the normal form coincide with those obtained from the expansion of the solution of the Kepler equation. In all cases where the perturbed problem can be cast in autonomous form, the solution is easily obtained as a perturbation series. The generalization to the spatial problem and/or the non-autonomous case is straightforward.     

A simplified vulnerability model for the extensive liquefaction risk assessment of buildings

Bulletin of Earthquake Engineering - The quantitative assessment of natural risks offers a rational strategy to protect communities, undertake cost effective mitigation and plan the organic and...     

Thermal convection in ice-I shells of Titan and Enceladus

Cassini-Huygens observations have shown that Titan and Enceladus are geologically active icy satellites. Mitri and Showman [Mitri, G., Showman, A.P., 2005. Icarus 177, 447-460] and McKinnon [McKinnon, W.B., 2006. Icarus 183, 435-450] investigated the dynamics of an ice shell overlying a pure liquid-water ocean and showed that transitions from a conductive state to a convective state have major implications for the surface tectonics. We extend this analysis to the case of ice shells overlying ammonia-water oceans. We explore the thermal state of Titan and Enceladus ice-I shells, and also we investigate the consequences of the ice-I shell conductive-convective switch for the geology. We show that thermal convection can occur, under a range of conditions, in the ice-I shells of Titan and Enceladus. Because the Rayleigh number Ra scales with δ³/ηb, where δ is the thickness of the ice shell and ηb is the viscosity at the base of the ice-I shell, and because ammonia in the liquid layer (if any) strongly depresses the melting temperature of the water ice, Ra equals its critical value for two ice-I shell thicknesses: for relatively thin ice shell with warm, low-viscosity base (Onset I) and for thick ice shell with cold, high-viscosity base (Onset II). At Onset I, for a range of heat fluxes, two equilibrium states—corresponding to a thin, conductive shell and a thick, convective shell—exist for a given heat flux. Switches between these states can cause large, rapid changes in the ice-shell thickness. For Enceladus, we demonstrate that an Onset I transition can produce tectonic stress of ∼500 bars and fractures of several tens of km depth. At Onset II, in contrast, we demonstrate that zero equilibrium states exist for a range of heat fluxes. For a mean heat flux within this range, the satellite experiences oscillations in surface heat flux and satellite volume with periods of ∼50-800 Myr even when the interior heat production is constant or monotonically declining in time; these oscillations in the thermal state of the ice-I shell would cause repeated episodes of extensional and compressional tectonism.     

Adiabatic chaos in the spin–orbit problem

We provide evidences that the angular momentum of a symmetric rigid body in a spin–orbit resonance can perform large scale chaotic motions on time scales which increase polynomially with the inverse of the oblateness of the body. This kind of irregular precession appears as soon as the orbit of the center of mass is non-circular and the angular momentum of the body is far from the principal directions with minimum (maximum) moment of inertia. We also provide a quantitative explanation of these facts by using the theory of adiabatic invariants, and we provide numerical applications to the cases of the 1:1 and 1:2 spin–orbit resonances.     

Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

旋转扰动流体中壁面湍流不稳定及其拟序结构研究Ⅱ：数值模式

对第一部分的实验结果进行了讨论,指出局地最大、最小涡度的存在与所谓的正压不稳定和大气中经常观测到的滚动涡有关。速度廓线中的拐点和涡度极大值点对应不稳定发展的位置。另外,为了改进湍流的模拟,根据“准正则”近似,发展了一个三阶闭合模式。这个模型已被证明可以描述与湍流过程有关的流体动力学变量的趋势。     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

Global change and relative sea level rise at Venice: what impact in term of flooding

Relative sea level rise (RSLR) due to climate change and geodynamics represents the main threat for the survival of Venice, emerging today only 90 cm above the Northern Adriatic mean sea level (msl). The 25 cm RSLR occurred over the 20th century, consisting of about 12 cm of land subsidence and 13 cm of sea level rise, has increased the flood frequency by more than seven times with severe damages to the urban heritage. Reasonable forecasts of the RSLR expected to the century end must be investigated to assess the suitability of the Mo.S.E. project planned for the city safeguarding, i.e., the closure of the lagoon inlets by mobile barriers. Here we consider three RSLR scenarios as resulting from the past sea level rise recorded in the Northern Adriatic Sea, the IPCC mid-range A1B scenario, and the expected land subsidence. Available sea level measurements show that more than 5 decades are required to compute a meaningful eustatic trend, due to pseudo-cyclic 7–8 year long fluctuations. The period from 1890 to 2007 is characterized by an average rate of 0.12 ± 0.01 cm/year. We demonstrate that linear regression is the most suitable model to represent the eustatic process over these 117 year. Concerning subsidence, at present Venice is sinking due to natural causes at 0.05 cm/year. The RSLR is expected to range between 17 and 53 cm by 2100, and its repercussions in terms of flooding frequency are associated here to each scenario. In particular, the frequency of tides higher than 110 cm, i.e., the value above which the gates would close the lagoon to the sea, will increase from the nowadays 4 times per year to a range between 20 and 250. These projections provide a large spread of possible conditions concerning the survival of Venice, from a moderate nuisance to an intolerable aggression. Hence, complementary solutions to Mo.S.E. may well be investigated.     

Magnesite growth rates as a function of temperature and saturation state

Magnesite growth rates and step velocities have been measured systematically as a function of temperature from 80 to 105 °C and saturation state in 0.1 M NaCl solutions using hydrothermal atomic force microscopy (HAFM). The observations indicate that at these conditions magnesite precipitation is dominated by the coupling of step generation via spiral growth at screw dislocations and step advancement away from these dislocations. As these two processes occur in series the slowest of these dominates precipitation rates. At 100 °C magnesite growth rates (r) determined by HAFM are consistent with

r=k(Ω-1)²,