Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Rarefied gas dynamics: fundamentals, simulations and micro flows

IP方法模拟硬盘读写头气膜润滑问题

本论文中采用IP方法模拟硬盘读写头的气膜润滑问题，并选择耗时但精准的DSMC方法作为该问题中IP方法的检验标准，IP方法计算得到的读写头表面压力分布及读写头所受净升力为读写头的设计提供可靠的参考数据。 首先，将读写头复杂表面简化为平面并假定尺寸为微米（远小于真实的毫米），微米尺寸时可以得到DSMC的模拟结果用于比较。IP法对该问题的模拟结果与DSMC的模拟结果完全吻合。 其次，选定一个真实的读写头表面几何，但仍然假定尺寸为微米。IP法在该问题的应用过程中遇到四点困难，我们提出或采取了相应的解决办法。在复杂表面几何的读写头气膜润滑问题中，IP法的模拟结果与DSMC的模拟结果也相符一致。 最后完成真实几何、真实尺寸（毫米）的读写头气膜润滑问题的模拟，这也是本论文的最终研究目标。IP法模拟该问题时需要采用大网格，为此我们构造二维检验模型验证了大网格的合理性。 真实几何读写头气膜润滑问题的模拟结果表明，在流场的局部区域压力沿着高度方向变化显著，这意味着读写头问题的传统计算方法（即概括性雷诺方程）在该局部区域不成立，故用它计算得到的结果值得认真检验，但文献中尚未给出过这类检验。这正是该论文的研究背景和采用新方法（IP方法）的原因。 本论文研究的问题属于过渡领域、低速流动问题，故本论文还介绍了该领域另外一种粒子模拟方法（LVDSMC方法）以及我们对该方法提出的一些改进。 关键词：硬盘读写头，气膜润滑，概括性雷诺方程，IP方法，DSMC方法，过渡领域，LVDSMC方法

EU Criminal Law and Enlargement

Since creation of the European Communities the number of Member States has gradually increased from the original six to current twenty-eight. Enlargement has become an EU’s flagship external policy, demonstrating the EU’s ability to shape its neighbourhood and to serve as a catalyst of deep and multilayered reforms. The consecutive seven enlargement rounds went in parallel with widespread internal developments, culminating with the creation of the European Union and, most recently, entry into force of the Treaty of Lisbon. As this volume demonstrates, EU criminal law has evolved considerably from its early days under the legal framework laid down by the Treaty of Maastricht to its current post-Lisbon shape. On 1 December 2014, that is with expiry of a five year transitional regime for the jurisdiction of the Court of Justice, Police and Judicial Co-operation in Criminal Matters became a fully fledged EU policy, governed largely by the same modus operandi as other areas of EU competence and with compulsory jurisdiction of the Court of Justice. As EU criminal law developed internally, so did its external dimension, including the role it plays in the enlargement policy. In case of the latter the expiry of the same transitional period has brought to an end a rather anomalous situation whereby the European Union had more enforcement tools before and after accession vis-à-vis its future/new Member States than it could employ against the old ones. This bifurcation, quite rightly, triggered a lot of discussions about double standards used by the European Union in its pre-accession policy. This is exacerbated by the fact that some of those standards are neither defined in EU law, nor pursued vis-à-vis the existing EU’s Member States. The aim of this chapter is to demonstrate that evolution with particular emphasis on the role of EU Criminal Law in the policy currently employed by the European Union vis-à-vis candidate and potential candidate countries of the Western Balkans and to Turkey. Arguably, together with political conditionality, it has become one of the pillars of the enlargement process and, as the examples of accession negotiations with Montenegro and Serbia prove, its role is likely to increase as rapprochement of other candidates and potential candidates progresses to the next stages.

El régimen de impuesto a la renta y la devolución del IVA a las empresas públicas bajo la perspectiva de la Ley Orgánica de Empresas Públicas (Tema Central)

La devolución del IVA consiste en un régimen tributario especial que se aplica en función de consideraciones específicas y en atención a la calidad que ostentan los sujetos pasivos que se benefician del mismo. Para que esta devolución proceda –en el caso de empresas públicas–, es indispensable que estas se encuentren constituidas de conformidad con las disposiciones de la Ley orgánica de Empresas Públicas o, en su defecto, hayan cumplido con las disposiciones específicas del régimen de transición contemplado para que las empresas públicas o estatales existen tes, las sociedades anónimas en las que el Estado es accionista, las empresas subsidiarias y las empresas de las Fuerzas armadas y de la Policía Nacional se conviertan en nuevas empresas públicas. En este contexto, el presente trabajo pretende analizar, de conformidad con el contenido de la Ley Orgánica de Empresas Públicas, el momento a partir del cual se genera el derecho a la devolución del Impuesto al valor agregado, para las empresas públicas, así como cuestionar ciertas inconsistencias que se desprenden de esta Ley.

Numerical simulations of microphysical processes in pyro-convective clouds

Hochreichende Konvektion über Waldbränden ist eine der intensivsten Formen von atmosphärischer Konvektion. Die extreme Wolkendynamik mit hohen vertikalen Windgeschwindigkeiten (bis 20 m/s) bereits an der Wolkenbasis, hohen Wasserdampfübersättigungen (bis 1%) und die durch das Feuer hohen Anzahlkonzentration von Aerosolpartikeln (bis 100000 cm^-3) bilden einen besonderen Rahmen für Aerosol-Wolken Wechselwirkungen.Ein entscheidender Schritt in der mikrophysikalischen Entwicklung einer konvektiven Wolke ist die Aktivierung von Aerosolpartikeln zu Wolkentropfen. Dieser Aktivierungsprozess bestimmt die anfängliche Anzahl und Größe der Wolkentropfen und kann daher die Entwicklung einer konvektiven Wolke und deren Niederschlagsbildung beeinflussen. Die wichtigsten Faktoren, welche die anfängliche Anzahl und Größe der Wolkentropfen bestimmen, sind die Größe und Hygroskopizität der an der Wolkenbasis verfügbaren Aerosolpartikel sowie die vertikale Windgeschwindigkeit. Um den Einfluss dieser Faktoren unter pyro-konvektiven Bedingungen zu untersuchen, wurden numerische Simulationen mit Hilfe eines Wolkenpaketmodells mit detaillierter spektraler Beschreibung der Wolkenmikrophysik durchgeführt. Diese Ergebnisse können in drei unterschiedliche Bereiche abhängig vom Verhältnis zwischen vertikaler Windgeschwindigkeit und Aerosolanzahlkonzentration (w/NCN) eingeteilt werden: (1) ein durch die Aerosolkonzentration limitierter Bereich (hohes w/NCN), (2) ein durch die vertikale Windgeschwindigkeit limitierter Bereich (niedriges w/NCN) und (3) ein Übergangsbereich (mittleres w/NCN). Die Ergebnisse zeigen, dass die Variabilität der anfänglichen Anzahlkonzentration der Wolkentropfen in (pyro-) konvektiven Wolken hauptsächlich durch die Variabilität der vertikalen Windgeschwindigkeit und der Aerosolkonzentration bestimmt wird. rnUm die mikrophysikalischen Prozesse innerhalb der rauchigen Aufwindregion einer pyrokonvektiven Wolke mit einer detaillierten spektralen Mikrophysik zu untersuchen, wurde das Paketmodel entlang einer Trajektorie innerhalb der Aufwindregion initialisiert. Diese Trajektore wurde durch dreidimensionale Simulationen eines pyro-konvektiven Ereignisses durch das Model ATHAM berechnet. Es zeigt sich, dass die Anzahlkonzentration der Wolkentropfen mit steigender Aerosolkonzentration ansteigt. Auf der anderen Seite verringert sich die Größe der Wolkentropfen mit steigender Aerosolkonzentration. Die Reduzierung der Verbreiterung des Tropfenspektrums stimmt mit den Ergebnissen aus Messungen überein und unterstützt das Konzept der Unterdrückung von Niederschlag in stark verschmutzen Wolken.Mit Hilfe des Models ATHAM wurden die dynamischen und mikrophysikalischen Prozesse von pyro-konvektiven Wolken, aufbauend auf einer realistischen Parametrisierung der Aktivierung von Aerosolpartikeln durch die Ergebnisse der Aktivierungsstudie, mit zwei- und dreidimensionalen Simulationen untersucht. Ein modernes zweimomenten mikrophysikalisches Schema wurde in ATHAM implementiert, um den Einfluss der Anzahlkonzentration von Aerosolpartikeln auf die Entwicklung von idealisierten pyro-konvektiven Wolken in US Standardamtosphären für die mittleren Breiten und den Tropen zu untersuchen. Die Ergebnisse zeigen, dass die Anzahlkonzentration der Aerosolpartikel die Bildung von Regen beeinflusst. Für geringe Aerosolkonzentrationen findet die rasche Regenbildung hauptsächlich durch warme mikrophysikalische Prozesse statt. Für höhere Aerosolkonzentrationen ist die Eisphase wichtiger für die Bildung von Regen. Dies führt zu einem verspäteten Einsetzen von Niederschlag für verunreinigtere Atmosphären. Außerdem wird gezeigt, dass die Zusammensetzung der Eisnukleationspartikel (IN) einen starken Einfluss auf die dynamische und mikrophysikalische Struktur solcher Wolken hat. Bei sehr effizienten IN bildet sich Regen früher. Die Untersuchung zum Einfluss des atmosphärischen Hintergrundprofils zeigt eine geringe Auswirkung der Meteorologie auf die Sensitivität der pyro-konvektiven Wolken auf diernAerosolkonzentration. Zum Abschluss wird gezeigt, dass die durch das Feuer emittierte Hitze einen deutlichen Einfluss auf die Entwicklung und die Wolkenobergrenze von pyro-konvektive Wolken hat. Zusammenfassend kann gesagt werden, dass in dieser Dissertation die Mikrophysik von pyrokonvektiven Wolken mit Hilfe von idealisierten Simulation eines Wolkenpaketmodell mit detaillierte spektraler Mikrophysik und eines 3D Modells mit einem zweimomenten Schema im Detail untersucht wurde. Es wird gezeigt, dass die extremen Bedingungen im Bezug auf die vertikale Windgeschwindigkeiten und Aerosolkonzentrationen einen deutlichen Einfluss auf die Entwicklung von pyro-konvektiven Wolken haben.

(Table 1) Stratigraphic datums from ODP Site 130-806

A 2 m.y. oxygen isotope record of Globigerinoides sacculifer from the Ontong Java Plateau, based on cores from Ocean Drilling Program Leg 130, is dated by matching variations to an orbital template. The procedure allows us to present the most complete Quaternary record available for the western equatorial Pacific. The template-generating algorithm describes a balance between growth and melting of ice. Following basic Milankovitch theory, ice growth is taken as constant, while melting is taken to depend on summer insolation, current ice mass, and average past ice mass. Template settings must be changed once, between 1 and 1.2 Ma, to reflect a major shift in climate. Template fits are strikingly good over much of the record and can be used to detect and fill gaps from core breaks and other disturbances. One result of template dating is an exact age for the Brunhes-Matuyama magnetic reversal boundary, at 790+/-5 ka, as well as several other precise dates (900 ka for the middle Pleistocene climate shift; 1070, 1240, and 1450 ka for isotope stages 31, 37, and 47, respectively). Sedimentation rates fluctuate between 18 and 28 m/m.y., a ca. 400 ka cycle being the most prominent. Major anomalies arise within the transitional regime (1.2 to 1 Ma). The origin of the cycles is unknown; we propose productivity variations in the western equatorial Pacific.

Código híbrido avanzado de motores de plasma de efecto Hall

La propulsión eléctrica constituye hoy una tecnología muy competitiva y de gran proyección de futuro. Dentro de los diversos motores de plasma existentes, el motor de efecto Hall ha adquirido una gran madurez y constituye un medio de propulsión idóneo para un rango amplio de misiones. En la presente Tesis se estudian los motores Hall con geometría convencional y paredes dieléctricas. La compleja interacción entre los múltiples fenómenos físicos presentes hace que sea difícil la simulación del plasma en estos motores. Los modelos híbridos son los que representan un mejor compromiso entre precisión y tiempo de cálculo. Se basan en utilizar un modelo fluido para los electrones y algoritmos de dinámica de partículas PIC (Particle-In- Cell) para los iones y los neutros. Permiten hacer uso de la hipótesis de cuasineutralidad del plasma, a cambio de resolver separadamente las capas límite (o vainas) que se forman en torno a las paredes de la cámara. Partiendo de un código híbrido existente, llamado HPHall-2, el objetivo de la Tesis doctoral ha sido el desarrollo de un código híbrido avanzado que mejorara la simulación de la descarga de plasma en un motor de efecto Hall. Las actualizaciones y mejoras realizadas en las diferentes partes que componen el código comprenden tanto aspectos teóricos como numéricos. Fruto de la extensa revisión de la algoritmia del código HPHall-2 se han conseguido reducir los errores de precisión un orden de magnitud, y se ha incrementado notablemente su consistencia y robustez, permitiendo la simulación del motor en un amplio rango de condiciones. Algunos aspectos relevantes a destacar en el subcódigo de partículas son: la implementación de un nuevo algoritmo de pesado que permite determinar de forma más precisa el flujo de las magnitudes del plasma; la implementación de un nuevo algoritmo de control de población, que permite tener suficiente número de partículas cerca de las paredes de la cámara, donde los gradientes son mayores y las condiciones de cálculo son más críticas; las mejoras en los balances de masa y energía; y un mejor cálculo del campo eléctrico en una malla no uniforme. Merece especial atención el cumplimiento de la condición de Bohm en el borde de vaina, que en los códigos híbridos representa una condición de contorno necesaria para obtener una solución consistente con el modelo de interacción plasma-pared, y que en HPHall-2 aún no se había resuelto satisfactoriamente. En esta Tesis se ha implementado el criterio cinético de Bohm para una población de iones con diferentes cargas eléctricas y una gran dispersión de velocidades. En el código, el cumplimiento de la condición cinética de Bohm se consigue por medio de un algoritmo que introduce una fina capa de aceleración nocolisional adyacente a la vaina y mide adecuadamente el flujo de partículas en el espacio y en el tiempo. Las mejoras realizadas en el subcódigo de electrones incrementan la capacidad de simulación del código, especialmente en la región aguas abajo del motor, donde se simula la neutralización del chorro del plasma por medio de un modelo de cátodo volumétrico. Sin abordar el estudio detallado de la turbulencia del plasma, se implementan modelos sencillos de ajuste de la difusión anómala de Bohm, que permiten reproducir los valores experimentales del potencial y la temperatura del plasma, así como la corriente de descarga del motor. En cuanto a los aspectos teóricos, se hace especial énfasis en la interacción plasma-pared y en la dinámica de los electrones secundarios libres en el interior del plasma, cuestiones que representan hoy en día problemas abiertos en la simulación de los motores Hall. Los nuevos modelos desarrollados buscan una imagen más fiel a la realidad. Así, se implementa el modelo de vaina de termalización parcial, que considera una función de distribución no-Maxwelliana para los electrones primarios y contabiliza unas pérdidas energéticas más cercanas a la realidad. Respecto a los electrones secundarios, se realiza un estudio cinético simplificado para evaluar su grado de confinamiento en el plasma, y mediante un modelo fluido en el límite no-colisional, se determinan las densidades y energías de los electrones secundarios libres, así como su posible efecto en la ionización. El resultado obtenido muestra que los electrones secundarios se pierden en las paredes rápidamente, por lo que su efecto en el plasma es despreciable, no así en las vainas, donde determinan el salto de potencial. Por último, el trabajo teórico y de simulación numérica se complementa con el trabajo experimental realizado en el Pnnceton Plasma Physics Laboratory, en el que se analiza el interesante transitorio inicial que experimenta el motor en el proceso de arranque. Del estudio se extrae que la presencia de gases residuales adheridos a las paredes juegan un papel relevante, y se recomienda, en general, la purga completa del motor antes del modo normal de operación. El resultado final de la investigación muestra que el código híbrido desarrollado representa una buena herramienta de simulación de un motor Hall. Reproduce adecuadamente la física del motor, proporcionando resultados similares a los experimentales, y demuestra ser un buen laboratorio numérico para estudiar el plasma en el interior del motor. Abstract Electric propulsion is today a very competitive technology and has a great projection into the future. Among the various existing plasma thrusters, the Hall effect thruster has acquired a considerable maturity and constitutes an ideal means of propulsion for a wide range of missions. In the present Thesis only Hall thrusters with conventional geometry and dielectric walls are studied. The complex interaction between multiple physical phenomena makes difficult the plasma simulation in these engines. Hybrid models are those representing a better compromise between precision and computational cost. They use a fluid model for electrons and Particle-In-Cell (PIC) algorithms for ions and neutrals. The hypothesis of plasma quasineutrality is invoked, which requires to solve separately the sheaths formed around the chamber walls. On the basis of an existing hybrid code, called HPHall-2, the aim of this doctoral Thesis is to develop an advanced hybrid code that better simulates the plasma discharge in a Hall effect thruster. Updates and improvements of the code include both theoretical and numerical issues. The extensive revision of the algorithms has succeeded in reducing the accuracy errors in one order of magnitude, and the consistency and robustness of the code have been notably increased, allowing the simulation of the thruster in a wide range of conditions. The most relevant achievements related to the particle subcode are: the implementation of a new weighing algorithm that determines more accurately the plasma flux magnitudes; the implementation of a new algorithm to control the particle population, assuring enough number of particles near the chamber walls, where there are strong gradients and the conditions to perform good computations are more critical; improvements in the mass and energy balances; and a new algorithm to compute the electric field in a non-uniform mesh. It deserves special attention the fulfilment of the Bohm condition at the edge of the sheath, which represents a boundary condition necessary to match consistently the hybrid code solution with the plasma-wall interaction, and remained as a question unsatisfactory solved in the HPHall-2 code. In this Thesis, the kinetic Bohm criterion has been implemented for an ion particle population with different electric charges and a large dispersion in their velocities. In the code, the fulfilment of the kinetic Bohm condition is accomplished by an algorithm that introduces a thin non-collisional layer next to the sheaths, producing the ion acceleration, and measures properly the flux of particles in time and space. The improvements made in the electron subcode increase the code simulation capabilities, specially in the region downstream of the thruster, where the neutralization of the plasma jet is simulated using a volumetric cathode model. Without addressing the detailed study of the plasma turbulence, simple models for a parametric adjustment of the anomalous Bohm difussion are implemented in the code. They allow to reproduce the experimental values of the plasma potential and the electron temperature, as well as the discharge current of the thruster. Regarding the theoretical issues, special emphasis has been made in the plasma-wall interaction of the thruster and in the dynamics of free secondary electrons within the plasma, questions that still remain unsolved in the simulation of Hall thrusters. The new developed models look for results closer to reality, such as the partial thermalization sheath model, that assumes a non-Maxwellian distribution functions for primary electrons, and better computes the energy losses at the walls. The evaluation of secondary electrons confinement within the chamber is addressed by a simplified kinetic study; and using a collisionless fluid model, the densities and energies of free secondary electrons are computed, as well as their effect on the plasma ionization. Simulations show that secondary electrons are quickly lost at walls, with a negligible effect in the bulk of the plasma, but they determine the potential fall at sheaths. Finally, numerical simulation and theoretical work is complemented by the experimental work carried out at the Princeton Plasma Physics Laboratory, devoted to analyze the interesting transitional regime experienced by the thruster in the startup process. It is concluded that the gas impurities adhered to the thruster walls play a relevant role in the transitional regime and, as a general recomendation, a complete purge of the thruster before starting its normal mode of operation it is suggested. The final result of the research conducted in this Thesis shows that the developed code represents a good tool for the simulation of Hall thrusters. The code reproduces properly the physics of the thruster, with results similar to the experimental ones, and represents a good numerical laboratory to study the plasma inside the thruster.

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.

Seafloor spreading magnetic anomaly isochron map compilation for the Weddell Sea and Scotia Sea

The deep sea sedimentary record is an archive of the pre-glacial to glacial development of Antarctica and changes in climate, tectonics and ocean circulation. Identification of the pre-glacial, transitional and full glacial components in the sedimentary record is necessary for ice sheet reconstruction and to build circum-Antarctic sediment thickness grids for past topography and bathymetry reconstructions, which constrain paleoclimate models. A ~3300 km long Weddell Sea to Scotia Sea transect consisting of multichannel seismic reflection data from various organisations, were used to interpret new horizons to define the initial basin-wide seismostratigraphy and to identify the pre-glacial to glacial components. We mapped seven main units of which three are in the inferred Cretaceous-Paleocene pre-glacial regime, one in the Eocene-Oligocene transitional regime and three units in the Miocene-Pleistocene full glacial climate regime. Sparse borehole data from ODP leg 113 and SHALDRIL constrain the ages of the upper three units. Compiled seafloor spreading magnetic anomalies constrain the basement ages and the hypothetical age model. In many cases, the new horizons and stratigraphy contradict the interpretations in local studies. Each seismic sedimentary unit and its associated base horizon are continuous and traceable for the entire transect length, but reflect a lateral change in age whilst representing the same deposition process. The up to 1240 m thick pre-glacial seismic units form a mound in the central Weddell Sea basin and, in conjunction with the eroded flank geometry, support the interpretation of a Cretaceous proto-Weddell Gyre. The base reflector of the transitional seismic unit, which marks the initial ice sheet advances to the outer shelf, has a lateral model age of 26.6-15.5 Ma from southeast to northwest. The Pliocene-Pleistocene glacial deposits reveals lower sedimentations rates, indicating a reduced sediment supply. Sedimentation rates for the pre-glacial, transitional and full glacial components are highest around the Antarctic Peninsula, indicating higher erosion and sediment supply of a younger basement. We interpret an Eocene East Antarctic Ice Sheet expansion, Oligocene grounding of the West Antarctic Ice Sheet and Early Miocene grounding of the Antarctic Peninsula Ice Sheet.

Railgun

Com a presente evolução das railguns na Marinha dos Estados Unidos da América e possível instalação em seus navios num futuro muito próximo, outras marinhas se seguirão. Creio que será do interesse da Marinha Portuguesa acompanhar esta evolução tecnológica, considerando as vantagens que advém da adoção deste tipo de armamento. Neste documento são abordados os princípios básicos subjacentes ao funcionamento da railgun, com principal foco nas questões eletrodinâmicas. Pretende-se adquirir familiaridade com este novo tipo de armamento através do estudo crítico dos seus princípios de funcionamento. O princípio básico de funcionamento de uma railgun, à primeira vista, parece bastante simples, à luz da aplicação imediata da expressão da força de Lorentz sobre um condutor percorrido por corrente elétrica. No entanto, tudo se torna mais complicado no caso de uma variação rápida dos parâmetros envolvidos (regime transitório), que exige uma análise mais aprofundada do comportamento da corrente, campos elétrico e magnético, e todos os materiais envolvidos neste sistema. Este trabalho envolveu ainda a construção de duas railguns, uma primeira de dimensões mais pequenas para ganhar familiaridade com o sistema, e uma última de dimensões de laboratório na qual foram feitos vários disparos para testar diferentes tipos de material e dimensões de projétil. Em suma, é demonstrado neste documento uma análise, no domínio do tempo, da distribuição espacial do campo eletromagnético, corrente elétrica e consequente fluxo de energia, complementados por uma parte experimental.

From Absence to Presence, from Silence to Voice: Victims in Transitional Justice Since the Nuremberg Trials

As much as victims have been absent in traditional and national criminal justice for a long time, they were invisible in transitional and international criminal justice after World War II. The Nuremberg Trials were dominated by the perpetrators, and documents were mainly used instead of victim testimony. Contemporaries shared the perspective that transitional justice, both international and national procedures should channel revenge by the victims and their families into the more peaceful venues of courts and legal procedures. Since then, victims have gained an ever more important role in transitional, post-conflict and international criminal justice. Non-judicial tribunals, Truth and Reconciliation Commissions, and international criminal courts and tribunals are relying on the testimony of victims and thus provide a prominent role for victims who often take centre stage in such procedures and trials. International criminal law and the human rights regime have provided victims with several routes to make themselves heard and fight against impunity. This paper tracks the road from absence to presence, and from invisibility to the visibility of victims during the second half of the last and the beginning of the present century. It shows in which ways their presence has shaped and changed transitional and international justice, and in particular how their absence or presence is linked to amnesties.