Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

High-Speed Laser Processing of Thin Single and Multi-Layer Films

Theoretical models are developed for the continuous-wave and pulsed laser incision and cut of thin single and multi-layer films. A one-dimensional steady-state model establishes the theoretical foundations of the problem by combining a power-balance integral with heat flow in the direction of laser motion. In this approach, classical modelling methods for laser processing are extended by introducing multi-layer optical absorption and thermal properties. The calculation domain is consequently divided in correspondence with the progressive removal of individual layers. A second, time-domain numerical model for the short-pulse laser ablation of metals accounts for changes in optical and thermal properties during a single laser pulse. With sufficient fluence, the target surface is heated towards its critical temperature and homogeneous boiling or "phase explosion" takes place. Improvements are seen over previous works with the more accurate calculation of optical absorption and shielding of the incident beam by the ablation products. A third, general time-domain numerical laser processing model combines ablation depth and energy absorption data from the short-pulse model with two-dimensional heat flow in an arbitrary multi-layer structure. Layer removal is the result of both progressive short-pulse ablation and classical vaporisation due to long-term heating of the sample. At low velocity, pulsed laser exposure of multi-layer films comprising aluminium-plastic and aluminium-paper are found to be characterised by short-pulse ablation of the metallic layer and vaporisation or degradation of the others due to thermal conduction from the former. At high velocity, all layers of the two films are ultimately removed by vaporisation or degradation as the average beam power is increased to achieve a complete cut. The transition velocity between the two characteristic removal types is shown to be a function of the pulse repetition rate. An experimental investigation validates the simulation results and provides new laser processing data for some typical packaging materials.

Aortic rupture in high-speed skiing crashes

BACKGROUND: New equipment and techniques in winter sports, such as carving skis and snowboards, have brought up new trauma patterns into the spectrum of leisure trauma. The injuries resemble high-energy trauma known from road crashes. The aim of the present study was to assess the incidence of acute traumatic descending aortic rupture in recreational skiing-crashes. MATERIAL: Between January 1995 and December 2004, 22 patients were admitted to our hospital for aortic rupture. Four patients had skiing crashes (18.2%). Mean age was 31 years, all patients were male. In two cases, aortic rupture was associated with fractures of the upper and lower extremities. One patient additionally had a cerebral contusion with an initial Glasgow Coma Scale score of 13. In two patients, isolated aortic rupture was diagnosed. RESULTS: Two patients were treated by graft interposition, and one by endograft. One patient arrived under mechanical resuscitation without blood pressure. He died at admission. He had been observed for 5 hours in another hospital, complaining of severe intrascapular back pain, before transport to our trauma unit for unknown bleeding. In the other three cases, treatment was successful. CONCLUSION: Rescue services and paramedics should be aware of this new type of injury. Acute aortic rupture has to be considered as possible injury in high velocity skiing crashes.

High fidelity simulations of failure in fiber-reinforced composites

Los ensayos virtuales de materiales compuestos han aparecido como un nuevo concepto dentro de la industria aeroespacial, y disponen de un vasto potencial para reducir los enormes costes de certificación y desarrollo asociados con las tediosas campañas experimentales, que incluyen un gran número de paneles, subcomponentes y componentes. El objetivo de los ensayos virtuales es sustituir algunos ensayos por simulaciones computacionales con alta fidelidad. Esta tesis es una contribución a la aproximación multiescala desarrollada en el Instituto IMDEA Materiales para predecir el comportamiento mecánico de un laminado de material compuesto dadas las propiedades de la lámina y la intercara. La mecánica de daño continuo (CDM) formula el daño intralaminar a nivel constitutivo de material. El modelo de daño intralaminar se combina con elementos cohesivos para representar daño interlaminar. Se desarrolló e implementó un modelo de daño continuo, y se aplicó a configuraciones simples de ensayos en laminados: impactos de baja y alta velocidad, ensayos de tracción, tests a cortadura. El análisis del método y la correlación con experimentos sugiere que los métodos son razonablemente adecuados para los test de impacto, pero insuficientes para el resto de ensayos. Para superar estas limitaciones de CDM, se ha mejorado la aproximación discreta de elementos finitos enriqueciendo la cinemática para incluir discontinuidades embebidas: el método extendido de los elementos finitos (X-FEM). Se adaptó X-FEM para un esquema explícito de integración temporal. El método es capaz de representar cualitativamente los mecanismos de fallo detallados en laminados. Sin embargo, los resultados muestran inconsistencias en la formulación que producen resultados cuantitativos erróneos. Por último, se ha revisado el método tradicional de X-FEM, y se ha desarrollado un nuevo método para superar sus limitaciones: el método cohesivo X-FEM estable. Las propiedades del nuevo método se estudiaron en detalle, y se concluyó que el método es robusto para implementación en códigos explícitos dinámicos escalables, resultando una nueva herramienta útil para la simulación de daño en composites. Virtual testing of composite materials has emerged as a new concept within the aerospace industry. It presents a very large potential to reduce the large certification costs and the long development times associated with the experimental campaigns, involving the testing of a large number of panels, sub-components and components. The aim of virtual testing is to replace some experimental tests by high-fidelity numerical simulations. This work is a contribution to the multiscale approach developed in Institute IMDEA Materials to predict the mechanical behavior of a composite laminate from the properties of the ply and the interply. Continuum Damage Mechanics (CDM) formulates intraply damage at the the material constitutive level. Intraply CDM is combined with cohesive elements to model interply damage. A CDM model was developed, implemented, and applied to simple mechanical tests of laminates: low and high velocity impact, tension of coupons, and shear deformation. The analysis of the results and the comparison with experiments indicated that the performance was reasonably good for the impact tests, but insuficient in the other cases. To overcome the limitations of CDM, the kinematics of the discrete finite element approximation was enhanced to include mesh embedded discontinuities, the eXtended Finite Element Method (X-FEM). The X-FEM was adapted to an explicit time integration scheme and was able to reproduce qualitatively the physical failure mechanisms in a composite laminate. However, the results revealed an inconsistency in the formulation that leads to erroneous quantitative results. Finally, the traditional X-FEM was reviewed, and a new method was developed to overcome its limitations, the stable cohesive X-FEM. The properties of the new method were studied in detail, and it was demonstrated that the new method was robust and can be implemented in a explicit finite element formulation, providing a new tool for damage simulation in composite materials.

A deformation model of flexible, high area-to-mass ratio debris under perturbations and validation method

A new type of space debris was recently discovered by Schildknecht in near -geosynchronous orbit (GEO). These objects were later identified as exhibiting properties associated with High Area-to-Mass ratio (HAMR) objects. According to their brightness magnitudes (light curve), high rotation rates and composition properties (albedo, amount of specular and diffuse reflection, colour, etc), it is thought that these objects are multilayer insulation (MLI). Observations have shown that this debris type is very sensitive to environmental disturbances, particularly solar radiation pressure, due to the fact that their shapes are easily deformed leading to changes in the Area-to-Mass ratio (AMR) over time. This thesis proposes a simple effective flexible model of the thin, deformable membrane with two different methods. Firstly, this debris is modelled with Finite Element Analysis (FEA) by using Bernoulli-Euler theory called “Bernoulli model”. The Bernoulli model is constructed with beam elements consisting 2 nodes and each node has six degrees of freedom (DoF). The mass of membrane is distributed in beam elements. Secondly, the debris based on multibody dynamics theory call “Multibody model” is modelled as a series of lump masses, connected through flexible joints, representing the flexibility of the membrane itself. The mass of the membrane, albeit low, is taken into account with lump masses in the joints. The dynamic equations for the masses, including the constraints defined by the connecting rigid rod, are derived using fundamental Newtonian mechanics. The physical properties of both flexible models required by the models (membrane density, reflectivity, composition, etc.), are assumed to be those of multilayer insulation. Both flexible membrane models are then propagated together with classical orbital and attitude equations of motion near GEO region to predict the orbital evolution under the perturbations of solar radiation pressure, Earth’s gravity field, luni-solar gravitational fields and self-shadowing effect. These results are then compared to two rigid body models (cannonball and flat rigid plate). In this investigation, when comparing with a rigid model, the evolutions of orbital elements of the flexible models indicate the difference of inclination and secular eccentricity evolutions, rapid irregular attitude motion and unstable cross-section area due to a deformation over time. Then, the Monte Carlo simulations by varying initial attitude dynamics and deformed angle are investigated and compared with rigid models over 100 days. As the results of the simulations, the different initial conditions provide unique orbital motions, which is significantly different in term of orbital motions of both rigid models. Furthermore, this thesis presents a methodology to determine the material dynamic properties of thin membranes and validates the deformation of the multibody model with real MLI materials. Experiments are performed in a high vacuum chamber (10-4 mbar) replicating space environment. A thin membrane is hinged at one end but free at the other. The free motion experiment, the first experiment, is a free vibration test to determine the damping coefficient and natural frequency of the thin membrane. In this test, the membrane is allowed to fall freely in the chamber with the motion tracked and captured through high velocity video frames. A Kalman filter technique is implemented in the tracking algorithm to reduce noise and increase the tracking accuracy of the oscillating motion. The forced motion experiment, the last test, is performed to determine the deformation characteristics of the object. A high power spotlight (500-2000W) is used to illuminate the MLI and the displacements are measured by means of a high resolution laser sensor. Finite Element Analysis (FEA) and multibody dynamics of the experimental setups are used for the validation of the flexible model by comparing with the experimental results of displacements and natural frequencies.

Bond Strength of Multicomponent White Cast Iron Coatings Applied by HVOF Thermal Spray Process

Multicomponent white cast iron is a new alloy that belongs to system Fe-C-Cr-W-Mo-V, and because of its excellent wear resistance it is used in the manufacture of hot rolling mills rolls. To date, this alloy has been processed by casting, powder metallurgy, and spray forming. The high-velocity oxyfuel process is now also considered for the manufacture of components with this alloy. The effects of substrate, preheating temperature, and coating thickness on bond strength of coatings have been determined. Substrates of AISI 1020 steel and of cast iron with preheating of 150 A degrees C and at room temperature were used to apply coatings with 200 and 400 mu m nominal thickness. The bond strength of coatings was measured with the pull-off test method and the failure mode by scanning electron microscopic analysis. Coatings with thickness of 200 mu m and applied on substrates of AISI 1020 steel with preheating presented bond strength of 87 +/- A 4 MPa.

The qWR star HD 45166 - II. Fundamental stellar parameters and evidence of a latitude-dependent wind

Context. The enigmatic object HD 45166 is a qWR star in a binary system with an orbital period of 1.596 day, and presents a rich emission-line spectrum in addition to absorption lines from the companion star (B7 V). As the system inclination is very small (i = 0.77 degrees +/- 0.09 degrees), HD 45166 is an ideal laboratory for wind-structure studies. Aims. The goal of the present paper is to determine the fundamental stellar and wind parameters of the qWR star. Methods. A radiative transfer model for the wind and photosphere of the qWR star was calculated using the non-LTE code CMFGEN. The wind asymmetry was also analyzed using a recently-developed version of CMFGEN to compute the emerging spectrum in two-dimensional geometry. The temporal-variance spectrum (TVS) was calculated to study the line-profile variations. Results. Abundances and stellar and wind parameters of the qWR star were obtained. The qWR star has an effective temperature of T(eff) = 50 000 +/- 2000 K, a luminosity of log(L/L(circle dot)) = 3.75 +/- 0.08, and a corresponding photospheric radius of R(phot) = 1.00 R(circle dot). The star is helium-rich (N(H)/N(He) = 2.0), while the CNO abundances are anomalous when compared either to solar values, to planetary nebulae, or to WR stars. The mass-loss rate is. M = 2.2 x 10(-7) M(circle dot) yr(-1), and the wind terminal velocity is v(infinity) = 425 km s(-1). The comparison between the observed line profiles and models computed under different latitude-dependent wind densities strongly suggests the presence of an oblate wind density enhancement, with a density contrast of at least 8: 1 from equator to pole. If a high velocity polar wind is present (similar to 1200 km s(-1)), the minimum density contrast is reduced to 4:1. Conclusions. The wind parameters determined are unusual when compared to O-type stars or to typical WR stars. While for WR stars v(infinity)/v(esc) > 1.5, in the case of HD 45166 it is much smaller (v(infinity)/v(esc) = 0.32). In addition, the efficiency of momentum transfer is eta = 0.74, which is at least 4 times smaller than in a typical WR. We find evidence for the presence of a wind compression zone, since the equatorial wind density is significantly higher than the polar wind. The TVS supports the presence of such a latitude-dependent wind and a variable absorption/scattering gas near the equator.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Proper motions of thermally emitting isolated neutron stars measured with Chandra

The remarkable astrometric capabilities of the Chandra Observatory offer the possibility to measure proper motions of X-ray sources with an unprecedented accuracy in this wavelength range. We recently completed a proper motion survey of three of the seven thermally emitting radio-quiet isolated neutron stars (INSs) discovered in the ROSAT all-sky survey. These INSs (RXJ0420.0-5022, RXJ0806.4-4123 and RXJ1308.6+2127) either lack an optical counterpart or have one so faint that ground based or space born optical observations push the current possibilities of the instrumentation to the limit. Pairs of ACIS observations were acquired 3 to 5 years apart to measure the displacement of the sources on the X-ray sky using as a reference the background of extragalactic or remote Galactic X-ray sources. We derive 2 sigma upper limits of 123 mas yr(-1) and 86 mas yr(-1) on the proper motion of RXJ0420.0-5022 and RXJ0806.4-4123, respectively. RXJ1308.6+2127 exhibits a very significant displacement (similar to 9 sigma) yielding mu = 220 +/- 25 mas yr(-1), the second fastest measured among all ROSAT-discovered INSs. The source is probably moving away rapidly from the Galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of similar to 740 (d/700 pc) km s(-1) might be the largest of all ROSAT INSs and its corresponding spatial velocity lies among the fastest recorded for neutron stars. RXJ1308.6+2127 is thus a middle-aged (age similar to 1 My) high velocity cooling neutron star. We investigate its possible origin in nearby OB associations or from a field OB star. In most cases, the flight time from birth place appears significantly shorter than the characteristic age derived from spin down rate. Overall, the distribution in transverse velocity of the ROSAT INSs is not statistically different from that of normal radio pulsars.

Effect of binders and surface finish on wear resistance of HVOF coatings

The high velocity oxygen fuel (HVOF) thermal spray process produces highly wear and/or corrosion resistant coatings. Tungsten carbide with a metallic binder is often used for this purpose. In this work, tungsten carbide coatings containing cobalt or nickel binder were produced by HVOF and characterised by optical and electron microscopy, hardness and a dry sand/rubber wheel abrasion test. The HVOF process produced dense coatings with low porosity levels and high hardness. The wear resistance of the specimens, which were surface treated, increased as the roughness percentage decreased. Tungsten carbide nickel based coating yielded the best wear resistance in the as sprayed condition. However, the wear rate and wear of the two coatings converged to the same values as the number of revolutions increased. Wear behaviour in the ground condition was similar, although the tungsten carbide cobalt based coating yielded better performance with increasing distance travelled during the wear test.

Influence of parameters of the HVOF thermal spray process on the properties of multicomponent white cast iron coatings

High velocity oxi-fuel (HVOF) thermal spray process has been used in order to deposit a new alloy known as multicomponent white cast iron. The coatings were characterized in terms of macrostructure, phase composition, porosity and hardness. Coating characteristics and properties were found to be dependent on the particles size range, spray distance, gases flow rate and oxygen to propane ratio. For set of parameters utilized in this job a narrow particle size range between 20 and 45 gm with a spray distance of 200 mm and oxygen to propane ratio of 4.6 are the preferred coating parameters. Coating porosity of 0.9% and hardness of 766 HV were obtained under these conditions. (c) 2007 Elsevier B.V. All rights reserved.

Abrasion and corrosion resistance of new Ni-based coating deposited by HVOF thermal spray process

Coatings based on NiCrAlC intermetallic based alloy were applied on AISI 316L stainless steel substrates using a high velocity oxygen fuel torch. The influence of the spray parameters on friction and abrasive wear resistance were investigated using an instrumented rubber wheel abrasion test, able to measure the friction forces. The corrosion behaviour of the coatings were studied with electrochemical techniques and compared with the corrosion resistance of the substrate material. Specimens prepared using lower O(2)/C(3)H(8) ratios showed smaller porosity values. The abrasion wear rate of the NiCrAlC coatings was much smaller than that described in the literature for bulk as cast materials with similar composition and one order of magnitude higher than bulk cast and heat treated (aged) NiCrAlC alloy. All coatings showed higher corrosion resistance than the AISI 316L substrate in HCl (5%) aqueous solution at 40 degrees C.

Free-Surface Aeration and Momentum Exchange at a Bottom Outlet

This study aims to provide some new understanding of the air-water flow properties in high-velocity water jets discharging past an abrupt drop. Such a setup has been little studied to date despite the relevance to bottom outlets. Downstream of the step brink, the free-jet entrains air at both upper and lower air-water interfaces, as well as along the sides. An air-water shear layer develops at the lower nappe interface. At the lower nappe, the velocity redistribution was successfully modelled and the velocity field was found to be similar to that in two-dimensional wake flow. The results highlighted further two distinct flow regions. Close to the brink (Wex < 5000), the flow was dominated by momentum transfer. Further downstream (Wex > 5000), a strong competition between air bubble diffusion and momentum exchanges took place.

Surface Waves and Roughness in Self-Aerated Supercritical Flow

In high-velocity open channel flows, free-surface aeration is commonly observed. The effects of surface waves on the air-water flow properties are tested herein. The study simulates the air-water flow past a fixed-location phase-detection probe by introducing random fluctuations of the flow depth. The present model yields results that are close to experimental observations in terms of void fraction, bubble count rate and bubble/droplet chord size distributions. The results show that the surface waves have relatively little impact on the void fraction profiles, but that the bubble count rate profiles and the distributions of bubble and chord sizes are affected by the presence of surface waves.

An extragalactic HI cloud with no optical counterpart?

We report the discovery, from the H I Parkes All-Sky Survey (HIPASS), of an isolated cloud of neutral hydrogen, which we believe to be extragalactic. The H I mass of the cloud (HIPASS J1712-64) is very low, 1.7 x 10(7) M-circle dot, using an estimated distance of similar to 3.2 Mpc. Most significantly, we have found no optical companion to this object to very faint limits [mu(B) similar to 27 mag arcsec(-2)]. HIPASS J1712-64 appears to be a binary system similar to, but much less massive than, H I 1225 + 01 (the Virgo H. I cloud) and has a size of at least 15 kpc. The mean velocity dispersion measured with the Australia Telescope Compact Array (ATCA) is only 4 km s(-1) for the main component and, because of the weak or nonexistent star formation, possibly reflects the thermal line width (T < 2000 K) rather than bulk motion or turbulence. The peak column density for HIPASS J1712-64, from the combined Parkes and ATCA data, is only 3.5 x 1019 cm(-2), which is estimated to be a factor of 2 below the critical threshold for star formation. Apart from its significantly higher velocity, the properties of HIPASS J1712-64 are similar to the recently recognized class of compact high-velocity clouds. We therefore consider the evidence for a Local Group or Galactic origin, although a more plausible alternative is that HIPASS J1712-64 was ejected from the interacting Magellanic Cloud-Galaxy system at perigalacticon similar to 2 x 10(8) yr ago.