High Fidelity Models for Near-Earth Object Dynamics

Motivado por los últimos hallazgos realizados gracias a los recientes avances tecnológicos y misiones espaciales, el estudio de los asteroides ha despertado el interés de la comunidad científica. Tal es así que las misiones a asteroides han proliferado en los últimos años (Hayabusa, Dawn, OSIRIX-REx, ARM, AIMS-DART, ...) incentivadas por su enorme interés científico. Los asteroides son constituyentes fundamentales en la evolución del Sistema Solar, son además grandes concentraciones de valiosos recursos naturales, y también pueden considerarse como objectivos estratégicos para la futura exploración espacial. Desde hace tiempo se viene especulando con la posibilidad de capturar objetos próximos a la Tierra (NEOs en su acrónimo anglosajón) y acercarlos a nuestro planeta, permitiendo así un acceso asequible a los mismos para estudiarlos in-situ, explotar sus recursos u otras finalidades. Por otro lado, las asteroides se consideran con frecuencia como posibles peligros de magnitud planetaria, ya que impactos de estos objetos con la Tierra suceden constantemente, y un asteroide suficientemente grande podría desencadenar eventos catastróficos. Pese a la gravedad de tales acontecimientos, lo cierto es que son ciertamente difíciles de predecir. De hecho, los ricos aspectos dinámicos de los asteroides, su modelado complejo y las incertidumbres observaciones hacen que predecir su posición futura con la precisión necesaria sea todo un reto. Este hecho se hace más relevante cuando los asteroides sufren encuentros próximos con la Tierra, y más aún cuando estos son recurrentes. En tales situaciones en las cuales fuera necesario tomar medidas para mitigar este tipo de riesgos, saber estimar con precisión sus trayectorias y probabilidades de colisión es de una importancia vital. Por ello, se necesitan herramientas avanzadas para modelar su dinámica y predecir sus órbitas con precisión, y son también necesarios nuevos conceptos tecnológicos para manipular sus órbitas llegado el caso. El objetivo de esta Tesis es proporcionar nuevos métodos, técnicas y soluciones para abordar estos retos. Las contribuciones de esta Tesis se engloban en dos áreas: una dedicada a la propagación numérica de asteroides, y otra a conceptos de deflexión y captura de asteroides. Por lo tanto, la primera parte de este documento presenta novedosos avances de apliación a la propagación dinámica de alta precisión de NEOs empleando métodos de regularización y perturbaciones, con especial énfasis en el método DROMO, mientras que la segunda parte expone ideas innovadoras para la captura de asteroides y comenta el uso del “ion beam shepherd” (IBS) como tecnología para deflectarlos. Abstract Driven by the latest discoveries enabled by recent technological advances and space missions, the study of asteroids has awakened the interest of the scientific community. In fact, asteroid missions have become very popular in the recent years (Hayabusa, Dawn, OSIRIX-REx, ARM, AIMS-DART, ...) motivated by their outstanding scientific interest. Asteroids are fundamental constituents in the evolution of the Solar System, can be seen as vast concentrations of valuable natural resources, and are also considered as strategic targets for the future of space exploration. For long it has been hypothesized with the possibility of capturing small near-Earth asteroids and delivering them to the vicinity of the Earth in order to allow an affordable access to them for in-situ science, resource utilization and other purposes. On the other side of the balance, asteroids are often seen as potential planetary hazards, since impacts with the Earth happen all the time, and eventually an asteroid large enough could trigger catastrophic events. In spite of the severity of such occurrences, they are also utterly hard to predict. In fact, the rich dynamical aspects of asteroids, their complex modeling and observational uncertainties make exceptionally challenging to predict their future position accurately enough. This becomes particularly relevant when asteroids exhibit close encounters with the Earth, and more so when these happen recurrently. In such situations, where mitigation measures may need to be taken, it is of paramount importance to be able to accurately estimate their trajectories and collision probabilities. As a consequence, advanced tools are needed to model their dynamics and accurately predict their orbits, as well as new technological concepts to manipulate their orbits if necessary. The goal of this Thesis is to provide new methods, techniques and solutions to address these challenges. The contributions of this Thesis fall into two areas: one devoted to the numerical propagation of asteroids, and another to asteroid deflection and capture concepts. Hence, the first part of the dissertation presents novel advances applicable to the high accuracy dynamical propagation of near-Earth asteroids using regularization and perturbations techniques, with a special emphasis in the DROMO method, whereas the second part exposes pioneering ideas for asteroid retrieval missions and discusses the use of an “ion beam shepherd” (IBS) for asteroid deflection purposes.

Análisis paramétrico determinista y aleatorio de problemas dinámicos en estructuras aeroespaciales

Los fenómenos dinámicos pueden poner en peligro la integridad de estructuras aeroespaciales y los ingenieros han desarrollado diferentes estrategias para analizarlos. Uno de los grandes problemas que se plantean en la ingeniería es cómo atacar un problema dinámico estructural. En la presente tesis se plantean distintos fenómenos dinámicos y se proponen métodos para estimar o simular sus comportamientos mediante un análisis paramétrico determinista y aleatorio del problema. Se han propuesto desde problemas sencillos con pocos grados de libertad que sirven para analizar las diferentes estrategias y herramientas a utilizar, hasta fenómenos muy dinámicos que contienen comportamientos no lineales, daños y fallos. Los primeros ejemplos de investigación planteados cubren una amplia gama de los fenómenos dinámicos, como el análisis de vibraciones de elementos másicos, incluyendo impactos y contactos, y el análisis de una viga con carga armónica aplicada a la que también se le añaden parámetros aleatorios que pueden responder a un desconocimiento o incertidumbre de los mismos. Durante el desarrollo de la tesis se introducen conceptos y se aplican distintos métodos, como el método de elementos finitos (FEM) en el que se analiza su resolución tanto por esquemas implícitos como explícitos, y métodos de análisis paramétricos y estadísticos mediante la técnica de Monte Carlo. Más adelante, una vez ya planteadas las herramientas y estrategias de análisis, se estudian fenómenos más complejos, como el impacto a baja velocidad en materiales compuestos, en el que se busca evaluar la resistencia residual y, por lo tanto, la tolerancia al daño de la estructura. Se trata de un suceso que puede producirse por la caída de herramienta, granizo o restos en la pista de aterrizaje. Otro de los fenómenos analizados también se da en un aeropuerto y se trata de la colisión con un dispositivo frangible, el cual tiene que romperse bajo ciertas cargas y, sin embargo, soportar otras. Finalmente, se aplica toda la metodología planteada en simular y analizar un posible incidente en vuelo, el fenómeno de la pérdida de pala de un turbohélice. Se trata de un suceso muy particular en el que la estructura tiene que soportar unas cargas complejas y excepcionales con las que la aeronave debe ser capaz de completar con éxito el vuelo. El análisis incluye comportamientos no lineales, daños, y varios tipos de fallos, y en el que se trata de identificar los parámetros clave en la secuencia del fallo. El suceso se analiza mediante análisis estructurales deterministas más habituales y también mediante otras técnicas como el método de Monte Carlo con el que se logran estudiar distintas incertidumbres en los parámetros con variables aleatorias. Se estudian, entre otros, el tamaño de pala perdida, la velocidad y el momento en el que se produce la rotura, y la rigidez y resistencia de los apoyos del motor. Se tiene en cuenta incluso el amortiguamiento estructural del sistema. Las distintas estrategias de análisis permiten obtener unos resultados valiosos e interesantes que han sido objeto de distintas publicaciones. ABSTRACT Dynamic phenomena can endanger the integrity of aerospace structures and, consequently, engineers have developed different strategies to analyze them. One of the major engineering problems is how to deal with the structural dynamics. In this thesis, different dynamic phenomena are introduced and several methods are proposed to estimate or simulate their behaviors. The analysis is considered through parametric, deterministic and statistical methods. The suggested issues are from simple problems with few degrees of freedom, in order to develop different strategies and tools to solve them, to very dynamic phenomena containing nonlinear behaviors failures, damages. The first examples cover a wide variety of dynamic phenomena such as vibration analysis of mass elements, including impacts and contacts, and beam analysis with harmonic load applied, in which random parameters are included. These parameters can represent the unawareness or uncertainty of certain variables. During the development of the thesis several concepts are introduced and different methods are applied, such as the finite element method (FEM), which is solved through implicit and explicit schemes, and parametrical and statistical methods using the Monte Carlo analysis technique. Next, once the tools and strategies of analysis are set out more complex phenomena are studied. This is the case of a low-speed impact in composite materials, the residual strength of the structure is evaluated, and therefore, its damage tolerance. This incident may occur from a tool dropped, hail or debris throw on the runway. At an airport may also occur, and it is also analyzed, a collision between an airplane and a frangible device. The devise must brake under these loads, however, it must withstand others. Finally, all the considered methodology is applied to simulate and analyze a flight incident, the blade loss phenomenon of a turboprop. In this particular event the structure must support complex and exceptional loads and the aircraft must be able to successfully complete the flight. Nonlinear behavior, damage, and different types of failures are included in the analysis, in which the key parameters in the failure sequence are identified. The incident is analyzed by deterministic structural analysis and also by other techniques such as Monte Carlo method, in which it is possible to include different parametric uncertainties through random variables. Some of the evaluated parameters are, among others, the blade loss size, propeller rotational frequency, speed and angular position where the blade is lost, and the stiffness and strength of the engine mounts. The study does also research on the structural damping of the system. The different strategies of analysis obtain valuable and interesting results that have been already published.

Estudo da cristalinidade de filmes finos de nitreto de índio e simulado pelo pacote de programas Wien2k

Neste trabalho, foi utilizado o método de deposição assistida por feixe de íons (IBAD na sigla em inglês) para produção de filmes finos de nitreto de índio em substratos de silício (111) e Safira-C. Variando as condições de deposição e utlilizando a técnica de difração de raios-X, investigou-se com o intuito de obter os parâmetros que resultam em filmes finos com melhor grau de cristalinidade. Os filmes produzidos a 380C apresentaram alta cristalinidade, superior àqueles a 250C. Temperaturas muito superiores a 380C não ocasionam a formação de filme cristalino de InN, como foi observado ao utilizar a temperatura de 480C; o mesmo se observa ao utilizar temperatura ambiente. Na temperatura considerada adequada ,de 380C, obteve-se que a utilização de Ra, ou seja, a razão de fluxo de partículas entre o nitrogênio e índio, em torno de 2,3 permite obter um melhor grau de cristalinização, o qual decresce conforme se diverge desse valor. A comparação entre difratogramas de amostras produzidas com e sem a evaporação prévia de titânio, o qual é possível observar um deslocamento dos picos do InN, indicam que o efeito Gettering permite a redução de impurezas no filme, principalmente de oxigênio. Utilizou-se a técnica de Retroespalhamento de Rutherford para obtenção da composição dos elementos e o perfil de profundidade. Notou-se uma forte mistura dos elementos do substrato de silício e safira com o nitreto de índio mesmo próximos a superfície. A presença indesejável de impurezas, principalmente o oxigênio, durante a deposição de filmes finos é praticamente inevitável. Desta forma, cálculos ab initio baseados na Teoria do Funcional da Densidade (DFT) foram realizados para investigar defeitos isolados e complexos de oxigênio no nitreto de índio e a sua influência nas propriedades óticas. Considerou-se diferentes concentrações de oxigênio (x=2,76, 8,32, 11,11 e 22,22%) aplicando-se o método PBEsolGGA e TB-mBJ para o tratamento da energia e potencial de troca e correlação. Obteve-se que é energeticamente favorável o oxigênio existir principalmente como defeito carregado e isolado. Os resultados utilizando a aproximação de TB-mBJ indicam um estreitamento do bandgap conforme a concentração de oxigênio aumenta. Entretanto, a alta contribuição do efeito de Moss-Burstein resulta num efetivo alargamento do band gap, gerando valores de band gap ótico maiores que no do bulk de nitreto de índio.

Systematic investigation of projectile fragmentation using beams of unstable B and C isotopes

Background: Models describing nuclear fragmentation and fragmentation fission deliver important input for planning nuclear physics experiments and future radioactive ion beam facilities. These models are usually benchmarked against data from stable beam experiments. In the future, two-step fragmentation reactions with exotic nuclei as stepping stones are a promising tool for reaching the most neutron-rich nuclei, creating a need for models to describe also these reactions. Purpose: We want to extend the presently available data on fragmentation reactions towards the light exotic region on the nuclear chart. Furthermore, we want to improve the understanding of projectile fragmentation especially for unstable isotopes. Method: We have measured projectile fragments from (10,12-18C) and B10-15 isotopes colliding with a carbon target. These measurements were all performed within one experiment, which gives rise to a very consistent data set. We compare our data to model calculations. Results: One-proton removal cross sections with different final neutron numbers (1 pxn) for relativistic C-10,C-12-18 and B10-15 isotopes impinging on a carbon target. Comparing model calculations to the data, we find that the EPAX code is not able to describe the data satisfactorily. Using ABRABLA07 on the other hand, we find that the average excitation energy per abraded nucleon needs to be decreased from 27 MeV to 8.1 MeV. With that decrease ABRABLA07 describes the data surprisingly well. Conclusions: Extending the available data towards light unstable nuclei with a consistent set of new data has allowed a systematic investigation of the role of the excitation energy induced in projectile fragmentation. Most striking is the apparent mass dependence of the average excitation energy per abraded nucleon. Nevertheless, this parameter, which has been related to final-state interactions, requires further study.

Rare earth and trace elements in igneous and metamorphic minerals of oceanic gabbros from the Mid-Atlantic Ridge

SIMS analyses have been carried out on clinopyroxenes, plagioclases and amphiboles of six gabbroic samples from Holes 921-924 of the Ocean Drilling Program Leg 153 sited in the MARK area of the Mid-Atlantic Ridge at the ridge-transform intersection, to investigate the rare earth, trace and volatile element distribution in the lower ocean crust during igneous crystallization and higher grade metamorphic conditions. The metagabbros underwent granulite to subgreenschist facies conditions through three main tectono-metamorphic phases: (1) ductile regime (750 < T < 1000 °C and P = 0.3 GPa); (2) transitional regime (600 < T < 700 °C and P = 0.2 GPa); (3) brittle regime (350 < T < 600 °C and P < 0.2 GPa). Igneous clinopyroxenes show Cl-chondrite normalized patterns depleted in LREE, and nearly flat for HREE. The rare earth and trace element distributions in igneous clinopyroxenes and plagioclases indicate that these minerals act as REE reservoirs, and comprise the main contribution to the overall rock content. The abundances in igneous minerals reflect the degree of fractionation of the parent liquids. In metamorphic clinopyroxenes recrystallized in anhydrous assemblages, the REE and trace elements patterns mimic those of the primary ones. Conversely, clinopyroxerie re-equilibrated in amphibolebearing assemblages shows a significant increase in REE, Ti, Zr, Y and V, a negative Eu anomaly, and slight decreases in Sr and Ba. An overall increase of REE and some trace elements is evident in hydrous assemblages, with preferential partitioning in the amphibole. It shows high Ti (18196-22844 ppm), LREE depleted patterns and LaN/SmN = 0.10-0.33, LaN/YbN = 0.10-0.30. Amphiboles from granoblastic assemblages show homogeneous patterns with no or a positive anomaly for TiN and negative anomalies for SrN and ZrN. Volatiles in amphibole are low, with Cl/F < 1; H2O% is significantly lower than the stoichiometric ratio (1.33-1.53%). The composition of the clinopyroxene and amphibole recrystallized in low-strain domains records evidence of incomplete re-equilibration, and element diffusion and partitioning is in part controlled by the textural site. The possible origins of the fluids involved in the metamorphic recrystallization are discussed: (1) remobilization from igneous amphibole; (2) exsolution from evolved melts; (3) introduction of seawater-derived fluids modified in rock-dominated systems; (4) injection of highly evolved hydrous melts during the metamorphic process.

Chemistry and boron isotope ratios of Izu arc front volcanic rocks

The process of fluid release from the subducting slab beneath the Izu arc volcanic front (Izu VF) was examined by measuring B concentrations and B isotope ratios in the Neogene fallout tephra (ODP Site 782A). Both were measured by secondary ion mass spectrometry, in a subset of matrix glasses and glassy plagioclase-hosted melt inclusions selected from material previously analyzed for major and trace elements (glasses) and radiogenic isotopes (Sr, Nd, Pb; bulk tephra). These tephra glasses have high B abundances (~10-60 ppm) and heavy delta11B values (+4.5? to +12.0?), extending the previously reported range for Izu VF rocks (delta11B, +7.0? to +7.3?). The glasses show striking negative correlations of delta11B with large ion lithophile element (LILE)/Nb ratios. These correlations cannot be explained by mixing two separate slab fluids, originating from the subducting sediment and the subducting basaltic crust, respectively (model A). Two alternative models (models B and C) are proposed. Model B proposes that the inverse correlations are inherited from altered oceanic crust (AOC), which shows a systematic decrease of B and LILE with increasing depth (from basaltic layer 2A to layer 3), paralleled by an increase in delta11B (from ~ +1? to +10? to +24?). In this model, the contribution of sedimentary B is insignificant (<4% of B in the Izu VF rocks). Model C explains the correlation as a mixture of a low-delta11B (~ +1?) 'composite' slab fluid (a mixture of metasediment- and metabasalt-derived fluids) with a metasomatized mantle wedge containing elevated B (~1-2 ppm) and heavy delta11B (~ +14?). The mantle wedge was likely metasomatized by 11B-rich fluids beneath the outer forearc, and subsequently down dragged to arc front depths by the descending slab. Pb-B isotope systematics indicate that, at arc front depths, ~ 53% of the B in the Izu VF is derived from the wedge. This implies that the heavy delta11B values of Izu VF rocks are largely a result of fluid fractionation, and do not reflect variations in slab source provenance (i.e. subducting sediment vs. basaltic crust). Since the B content of the peridotite at the outer forearc (7-58 ppm B, mean 24 +/- 16 ppm) is much higher than beneath the arc front (~1-2 ppm B), the hydrated mantle wedge must have released a B-rich fluid on its downward path. This 'wedge flux' can explain (1) the across-arc decrease in B and delta11B (e.g. Izu, Kuriles), without requiring a progressive decrease in fluid flux from the subducting slab, and (2) the thermal structure of volcanic arcs, as reflected in the B and delta11B variations of volcanic arc rocks.

Major and trace element composition of clinopyroxenes from an olivine gabbro of ODP Hole 176-735B (Table 1)

Major and trace element profiles of clinopyroxene grains in oceanic gabbros from ODP Hole 735B have been investigated by a combined in situ analytical study with ion probe, and electron microprobe. In contrast to the homogeneous major element compositions, trace elements (REE, Y, Cr, Sr, and Zr) show continuous core to rim zoning profiles. The observed trace element systematics in clinopyroxene cannot be explained by a simple diffusive exchange between melts and gabbros along grain boundaries. A simultaneous modification of the melt composition is required to generate the zoning, although Rayleigh fractional crystallization modelling could mimic the general shape of the profiles. Simultaneous metasomatism between the cumulate crystal and the porous melt during crystal accumulation is the most likely process to explain the zoning. Deformation during solidification of the crystal mush could have caused squeezing out of the incompatible element enriched residual melts (interstitial liquid). Migration of the melt along grain boundaries might carry these melt out of the system. This process named as synkinematic differentiation or differentiation by deformation (Natland and Dick, 2001, doi:10.1016/S0377-0273(01)00211-6) may act as an important magma evolution mechanism in the oceanic crust, at least at slow-spreading ridges.

Major and trace element concentration of melt inclusions from Zhupanovsky volcano, Kamchatka

The paper presents data on naturally quenched melt inclusions in olivine (Fo 69-84) from Late Pleistocene pyroclastic rocks of Zhupanovsky volcano in the frontal zone of the Eastern Volcanic Belt of Kamchatka. The composition of the melt inclusions provides insight into the latest crystallization stages (~70% crystallization) of the parental melt (~46.4 wt % SiO2, ~2.5 wt % H2O, ~0.3 wt % S), which proceeded at decompression and started at a depth of approximately 10 km from the surface. The crystallization temperature was estimated at 1100 ± 20°C at an oxygen fugacity of deltaFMQ = 0.9-1.7. The melts evolved due to the simultaneous crystallization of olivine, plagioclase, pyroxene, chromite, and magnetite (Ol: Pl: Cpx : (Crt-Mt) ~ 13 : 54 : 24 : 4) along the tholeiite evolutionary trend and became progressively enriched in FeO, SiO2, Na2O, and K2O and depleted in MgO, CaO, and Al2O3. Melt crystallization was associated with the segregation of fluid rich in S-bearing compounds and, to a lesser extent, in H2O and Cl. The primary melt of Zhupanovsky volcano (whose composition was estimated from data on the most primitive melt inclusions) had a composition of low-Si (~45 wt % SiO2) picrobasalt (~14 wt % MgO), as is typical of parental melts in Kamchatka and other island arcs, and was different from MORB. This primary melt could be derived by ~8% melting of mantle peridotite of composition close to the MORB source, under pressures of 1.5 ± 0.2 GPa and temperatures 20-30°C lower than the solidus temperature of 'dry' peridotite (1230-1240°C). Melting was induced by the interaction of the hot peridotite with a hydrous component that was brought to the mantle from the subducted slab and was also responsible for the enrichment of the Zhupanovsky magmas in LREE, LILE, B, Cl, Th, U, and Pb. The hydrous component in the magma source of Zhupanovsky volcano was produced by the partial slab melting under water-saturated conditions at temperatures of 760-810°C and pressures of ~3.5 GPa. As the depth of the subducted slab beneath Kamchatkan volcanoes varies from 100 to 125 km, the composition of the hydrous component drastically changes from relatively low-temperature H2O-rich fluid to higher temperature H2O-bearing melt. The geothermal gradient at the surface of the slab within the depth range of 100-125 km beneath Kamchatka was estimated at 4°C/km.

Progress in silicon-based quantum computation

We review progress at the Australian Centre for Quantum Computer Technology towards the fabrication and demonstration of spin qubits and charge qubits based on phosphorus donor atoms embedded in intrinsic silicon. Fabrication is being pursued via two complementary pathways: a 'top-down' approach for near-term production of few-qubit demonstration devices and a 'bottom-up' approach for large-scale qubit arrays with sub-nanometre precision. The 'top-down' approach employs a low-energy (keV) ion beam to implant the phosphorus atoms. Single-atom control during implantation is achieved by monitoring on-chip detector electrodes, integrated within the device structure. In contrast, the 'bottom-up' approach uses scanning tunnelling microscope lithography and epitaxial silicon overgrowth to construct devices at an atomic scale. In both cases, surface electrodes control the qubit using voltage pulses, and dual single-electron transistors operating near the quantum limit provide fast read-out with spurious-signal rejection.

Transmission electron microscopy sample preparation technique for sintered alloys

Powder metallurgy alloys are typically inhomogeneous with a significant amount of porosity. This complicates conventional transmission electron microscopy sample preparation. However, the use of focused ion beam milling allows site specific transmission electron microscopy samples to be prepared in a short amount of time. This paper presents a method that can be used to produce transmission electron microscopy samples from an Al-Cu-Mg PM alloy. (C) 2003 IoM Communications Ltd. Published by Maney for the Institute of Materials, Minerals and Mining.

Grain size measurements in Mg-Al high pressure die castings using electron back-scattered diffraction (EBSD)

Optical metallographic techniques for grain-size measurement give unreliable results for high pressure diecast Mg-Al alloys and electron back-scattered diffraction mapping (EBSD) provides a good tool for improving the quality of these measurements. An application of EBSD mapping to this question is described, and data for some castings are presented. Ion-beam milling was needed to prepare suitable samples, and this technique is detailed. As is well-known for high pressure die castings, the grain size distribution comprises at least two populations. The mean grain size of the fine-grained population was similar in both AZ91 and AM60 and in two casting thicknesses (2 mm and 5 mm) and, contrary to previously published reports, it did not vary with depth below the surface.

Aluminium phosphide as a eutectic grain nucleus in hypoeutectic Al-Si alloys

Aluminium phosphide (AlP) particles arc often suggested to be the nucleation site for eutectic silicon in Al-Si alloys, since both the crystal structure and lattice parameter of AlP (crystal structure: cubic K(4) over bar m; lattice parameter: 5.421 Angstrom) are close to that of silicon (cubic Fd3m, 5.431 Angstrom), and the melting point is higher than the Al-Si eutectic temperature. However, the crystallographic relationships between AlP particles and the surrounding eutectic silicon are seldom reported due to the difficulty in analysing the AlP particles, which react with water during sample preparation for polishing. in this study, the orientation relationships between AlP and Si are analysed by transmission electron microscopy using focused ion-beam milling for sample preparation to investigate the nucleation mechanism of eutectic silicon on AlP. The results show a clear and direct lattice relationship between centrally located AlP particles and the surrounding silicon in the hypoeutectic Al-Si alloy.