High power beam dump project for the accelerator prototype LIPAc: Cooling design and analysis

En el campo de la fusión nuclear y desarrollándose en paralelo a ITER (International Thermonuclear Experimental Reactor), el proyecto IFMIF (International Fusion Material Irradiation Facility) se enmarca dentro de las actividades complementarias encaminadas a solucionar las barreras tecnológicas que aún plantea la fusión. En concreto IFMIF es una instalación de irradiación cuya misión es caracterizar materiales resistentes a condiciones extremas como las esperadas en los futuros reactores de fusión como DEMO (DEMOnstration power plant). Consiste de dos aceleradores de deuterones que proporcionan un haz de 125 mA y 40 MeV cada uno, que al colisionar con un blanco de litio producen un flujo neutrónico intenso (1017 neutrones/s) con un espectro similar al de los neutrones de fusión [1], [2]. Dicho flujo neutrónico es empleado para irradiar los diferentes materiales candidatos a ser empleados en reactores de fusión, y las muestras son posteriormente examinadas en la llamada instalación de post-irradiación. Como primer paso en tan ambicioso proyecto, una fase de validación y diseño llamada IFMIFEVEDA (Engineering Validation and Engineering Design Activities) se encuentra actualmente en desarrollo. Una de las actividades contempladas en esta fase es la construcción y operación de una acelarador prototipo llamado LIPAc (Linear IFMIF Prototype Accelerator). Se trata de un acelerador de deuterones de alta intensidad idéntico a la parte de baja energía de los aceleradores de IFMIF. Los componentes del LIPAc, que será instalado en Japón, son suministrados por diferentes países europeos. El acelerador proporcionará un haz continuo de deuterones de 9 MeV con una potencia de 1.125 MW que tras ser caracterizado con diversos instrumentos deberá pararse de forma segura. Para ello se requiere un sistema denominado bloque de parada (Beam Dump en inglés) que absorba la energía del haz y la transfiera a un sumidero de calor. España tiene el compromiso de suministrar este componente y CIEMAT (Centro de Investigaciones Energéticas Medioambientales y Tecnológicas) es responsable de dicha tarea. La pieza central del bloque de parada, donde se para el haz de iones, es un cono de cobre con un ángulo de 3.5o, 2.5 m de longitud y 5 mm de espesor. Dicha pieza está refrigerada por agua que fluye en su superficie externa por el canal que se forma entre el cono de cobre y otra pieza concéntrica con éste. Este es el marco en que se desarrolla la presente tesis, cuyo objeto es el diseño del sistema de refrigeración del bloque de parada del LIPAc. El diseño se ha realizado utilizando un modelo simplificado unidimensional. Se han obtenido los parámetros del agua (presión, caudal, pérdida de carga) y la geometría requerida en el canal de refrigeración (anchura, rugosidad) para garantizar la correcta refrigeración del bloque de parada. Se ha comprobado que el diseño permite variaciones del haz respecto a la situación nominal siendo el flujo crítico calorífico al menos 2 veces superior al nominal. Se han realizado asimismo simulaciones fluidodinámicas 3D con ANSYS-CFX en aquellas zonas del canal de refrigeración que lo requieren. El bloque de parada se activará como consecuencia de la interacción del haz de partículas lo que impide cualquier cambio o reparación una vez comenzada la operación del acelerador. Por ello el diseño ha de ser muy robusto y todas las hipótesis utilizadas en la realización de éste deben ser cuidadosamente comprobadas. Gran parte del esfuerzo de la tesis se centra en la estimación del coeficiente de transferencia de calor que es determinante en los resultados obtenidos, y que se emplea además como condición de contorno en los cálculos mecánicos. Para ello por un lado se han buscado correlaciones cuyo rango de aplicabilidad sea adecuado para las condiciones del bloque de parada (canal anular, diferencias de temperatura agua-pared de decenas de grados). En un segundo paso se han comparado los coeficientes de película obtenidos a partir de la correlación seleccionada (Petukhov-Gnielinski) con los que se deducen de simulaciones fluidodinámicas, obteniendo resultados satisfactorios. Por último se ha realizado una validación experimental utilizando un prototipo y un circuito hidráulico que proporciona un flujo de agua con los parámetros requeridos en el bloque de parada. Tras varios intentos y mejoras en el experimento se han obtenido los coeficientes de película para distintos caudales y potencias de calentamiento. Teniendo en cuenta la incertidumbre de las medidas, los valores experimentales concuerdan razonablemente bien (en el rango de 15%) con los deducidos de las correlaciones. Por motivos radiológicos es necesario controlar la calidad del agua de refrigeración y minimizar la corrosión del cobre. Tras un estudio bibliográfico se identificaron los parámetros del agua más adecuados (conductividad, pH y concentración de oxígeno disuelto). Como parte de la tesis se ha realizado asimismo un estudio de la corrosión del circuito de refrigeración del bloque de parada con el doble fin de determinar si puede poner en riesgo la integridad del componente, y de obtener una estimación de la velocidad de corrosión para dimensionar el sistema de purificación del agua. Se ha utilizado el código TRACT (TRansport and ACTivation code) adaptándalo al caso del bloque de parada, para lo cual se trabajó con el responsable (Panos Karditsas) del código en Culham (UKAEA). Los resultados confirman que la corrosión del cobre en las condiciones seleccionadas no supone un problema. La Tesis se encuentra estructurada de la siguiente manera: En el primer capítulo se realiza una introducción de los proyectos IFMIF y LIPAc dentro de los cuales se enmarca esta Tesis. Además se describe el bloque de parada, siendo el diseño del sistema de rerigeración de éste el principal objetivo de la Tesis. En el segundo y tercer capítulo se realiza un resumen de la base teórica así como de las diferentes herramientas empleadas en el diseño del sistema de refrigeración. El capítulo cuarto presenta los resultados del relativos al sistema de refrigeración. Tanto los obtenidos del estudio unidimensional, como los obtenidos de las simulaciones fluidodinámicas 3D mediante el empleo del código ANSYS-CFX. En el quinto capítulo se presentan los resultados referentes al análisis de corrosión del circuito de refrigeración del bloque de parada. El capítulo seis se centra en la descripción del montaje experimental para la obtención de los valores de pérdida de carga y coeficiente de transferencia del calor. Asimismo se presentan los resultados obtenidos en dichos experimentos. Finalmente encontramos un capítulo de apéndices en el que se describen una serie de experimentos llevados a cabo como pasos intermedios en la obtención del resultado experimental del coeficiente de película. También se presenta el código informático empleado para el análisis unidimensional del sistema de refrigeración del bloque de parada llamado CHICA (Cooling and Heating Interaction and Corrosion Analysis). ABSTRACT In the nuclear fusion field running in parallel to ITER (International Thermonuclear Experimental Reactor) as one of the complementary activities headed towards solving the technological barriers, IFMIF (International Fusion Material Irradiation Facility) project aims to provide an irradiation facility to qualify advanced materials resistant to extreme conditions like the ones expected in future fusion reactors like DEMO (DEMOnstration Power Plant). IFMIF consists of two constant wave deuteron accelerators delivering a 125 mA and 40 MeV beam each that will collide on a lithium target producing an intense neutron fluence (1017 neutrons/s) with a similar spectra to that of fusion neutrons [1], [2]. This neutron flux is employed to irradiate the different material candidates to be employed in the future fusion reactors, and the samples examined after irradiation at the so called post-irradiative facilities. As a first step in such an ambitious project, an engineering validation and engineering design activity phase called IFMIF-EVEDA (Engineering Validation and Engineering Design Activities) is presently going on. One of the activities consists on the construction and operation of an accelerator prototype named LIPAc (Linear IFMIF Prototype Accelerator). It is a high intensity deuteron accelerator identical to the low energy part of the IFMIF accelerators. The LIPAc components, which will be installed in Japan, are delivered by different european countries. The accelerator supplies a 9 MeV constant wave beam of deuterons with a power of 1.125 MW, which after being characterized by different instruments has to be stopped safely. For such task a beam dump to absorb the beam energy and take it to a heat sink is needed. Spain has the compromise of delivering such device and CIEMAT (Centro de Investigaciones Energéticas Medioambientales y Tecnológicas) is responsible for such task. The central piece of the beam dump, where the ion beam is stopped, is a copper cone with an angle of 3.5o, 2.5 m long and 5 mm width. This part is cooled by water flowing on its external surface through the channel formed between the copper cone and a concentric piece with the latter. The thesis is developed in this realm, and its objective is designing the LIPAc beam dump cooling system. The design has been performed employing a simplified one dimensional model. The water parameters (pressure, flow, pressure loss) and the required annular channel geometry (width, rugoisty) have been obtained guaranteeing the correct cooling of the beam dump. It has been checked that the cooling design allows variations of the the beam with respect to the nominal position, being the CHF (Critical Heat Flux) at least twice times higher than the nominal deposited heat flux. 3D fluid dynamic simulations employing ANSYS-CFX code in the beam dump cooling channel sections which require a more thorough study have also been performed. The beam dump will activateasaconsequenceofthe deuteron beam interaction, making impossible any change or maintenance task once the accelerator operation has started. Hence the design has to be very robust and all the hypotheses employed in the design mustbecarefully checked. Most of the work in the thesis is concentrated in estimating the heat transfer coefficient which is decisive in the obtained results, and is also employed as boundary condition in the mechanical analysis. For such task, correlations which applicability range is the adequate for the beam dump conditions (annular channel, water-surface temperature differences of tens of degrees) have been compiled. In a second step the heat transfer coefficients obtained from the selected correlation (Petukhov- Gnielinski) have been compared with the ones deduced from the 3D fluid dynamic simulations, obtaining satisfactory results. Finally an experimental validation has been performed employing a prototype and a hydraulic circuit that supplies a flow with the requested parameters in the beam dump. After several tries and improvements in the experiment, the heat transfer coefficients for different flows and heating powers have been obtained. Considering the uncertainty in the measurements the experimental values agree reasonably well (in the order of 15%) with the ones obtained from the correlations. Due to radiological reasons the quality of the cooling water must be controlled, hence minimizing the copper corrosion. After performing a bibligraphic study the most adequate water parameters were identified (conductivity, pH and dissolved oxygen concentration). As part of this thesis a corrosion study of the beam dump cooling circuit has been performed with the double aim of determining if corrosion can pose a risk for the copper beam dump , and obtaining an estimation of the corrosion velocitytodimension the water purification system. TRACT code(TRansport and ACTivation) has been employed for such study adapting the code for the beam dump case. For such study a collaboration with the code responsible (Panos Karditsas) at Culham (UKAEA) was established. The work developed in this thesis has supposed the publication of three articles in JCR journals (”Journal of Nuclear Materials” y ”Fusion Engineering and Design”), as well as presentations in more than four conferences and relevant meetings.

Stable isotope record of Cibicidoides wuellerstorfi from Pleistocene sediments of DSDP Site 86-502 in the Caribbean Sea

Variations in the stable isotopic composition of benthic foraminifera from Deep-Sea Drilling Project (DSDP) site 502B in the Caribbean Sea are used to reconstruct Atlantic intermediate water circulation variability over the last 1.2 m.y. Comparison of this record with other North Atlantic benthic isotope records indicates that Atlantic intermediate water circulation was relatively enhanced during glacial maxima when North Atlantic deep water (NADW) production was reduced. However, a simple, compensatory relationship between intermediate and deepwater circulation is not apparent. Geochemical models have shown that such changes in ocean circulation can affect atmospheric CO2 levels by changing vertical nutrient and alkalinity profiles. The Delta delta13C difference between Caribbean site 502B and deep equatorial Pacific site 677 is highly coherent and in phase with ice volume. Like the delta18O record, there is an increase in amplitude (40%) and a large increase in 100 kyr power after 0.7 Ma. The 1.2? Delta delta13C amplitude scales to 70 ppm V in atmospheric CO2 using Boyle's (1986) box model result. The implied increase in CO2 amplitude after 0.7 Ma may suggest a positive feedback role in effecting the higher-amplitude climatic fluctuations which characterize the last 0.7 m.y.

Spectral analysis of seabed elevation from the Danish Wadden Sea

Bedforms both reflect and influence shallow water hydrodynamics and sediment dynamics. A correct characterization of their spatial distribution and dimensions is required for the understanding, assessment and prediction of numerous coastal processes. A method to parameterize geometrical characteristics using two-dimensional (2D) spectral analysis is presented and tested on seabed elevation data from the Knudedyb tidal inlet in the Danish Wadden Sea, where large compound bedforms are found. The bathymetric data were divided into 20x20 m areas on which a 2D spectral analysis was applied. The most energetic peak of the 2D spectrum was found and its energy, frequency and direction were calculated. A power-law was fitted to the average of slices taken through the 2D spectrum; its slope and y-intercept were calculated. Using these results the test area was morphologically classified into 4 distinct morphological regions. The most energetic peak and the slope and intercept of the power-law showed high values above the crest of the primary bedforms and scour holes, low values in areas without bedforms, and intermediate values in areas with secondary bedforms. The secondary bedform dimensions and orientations were calculated. An area of 700x700 m was used to determine the characteristics of the primary bedforms. However, they were less distinctively characterized compared to the secondary bedforms due to relatively large variations in their orientations and wavelengths. The method is thus appropriate for morphological classification of the seabed and for bedform characterization, being most efficient in areas characterized by bedforms with regular dimensions and directions.

(Table 1) Calculated sea surface temperature of sediment core RC11-210

Proxy indicators of sea surface temperature and equatorial divergence based on radiolarian assemblage data, and of trade wind intensity based on eolian grain size data show similar aspects of variability during the late Pleistocene: All indicators fluctuate at higher frequencies than the 100,000-year glacial-interglacial cycle, display reduced amplitude variations since 300,000 years ago, exhibit a change in the record character at about 300,000 years ago (the mid-Brunhes climatic event), and have higher amplitude variations in sediments 300,000-850,000 years old. Time series analyses were conducted to determine the spectral character of each record (delta18O of planktonic foraminifer, sea surface temperature values, equatorial divergence indicators, and wind intensity indicators) and to quantify interrecord coherence and phase relationships. The record was divided at the 300,000-year clear change in climatic variability (nonstationarity). The delta18O-based time scale is better lower in the core so our spectral analyses concentrated on the interval from 402,000-774,000 years. The delta18O spectra show 100,000- and 41,000-year power in the younger portion, 0-300,000 years, and 100,000-, 41,000- and 23,000-year power in the older interval, all highly coherent and in phase with the SPECMAP average stacked isotope record. Unlike the isotope record the dominant period in both the eolian grain size and equatorial divergence indicators is 31,000 years. This period is also important in the sea surface temperature signal where the dominant spectral peak is 100,000 years. The 31,000-year spectral component is coherent and in phase between the eolian and divergence records, confirming the link between atmospheric and ocean surface circulation for the first time in the paleoclimate record. Since the 31,000-year power appears in independent data sets within this core and also appears in other equatorial records [J. Imbrie personal communication, 1987], we assume it to be real and representative of both a nonlinear response to orbital forcing, possibly a combination of orbital tilt and eccentricity, and some resonance phenomenon required to amplify the response at this period so that it appears as a dominant frequency component. The mid-Brunhes climatic event is an important aspect of these records, but its cause remains unknown.

The natural history of fishes and serpents, including sea-turtles, crustaceous and shell fishes, with their medicinal uses. Illustrated with cuts.

Vol. III; of his new and accurate system of natural history.

England's improvement by sea and land : to out-do the Dutch without fighting, to pay debts without moneys, to set at work all the poor of England with the growth of our own lands : to prevent unnecessary suits in law, with the benefit of a voluntary register : directions where vast quantities of timber are to be had for the building of ships, with the advantage of making the great rivers of England navigable : rules to prevent fires in London, and other great cities : with directions how the several companies of handicraftsmen in London may always have cheap bread and drink /

Vol. 2 has title: England's improvement by sea and land. The second part : Containing, I. An account of its scituation, and the growths, and manufactures thereof. II. The benefit and necessity of a voluntary-register. III. A method for improving the Royal-navy, lessening the growing power of France, and obtaining the fishery. IV. Advantageous proposals for the city of London, for the preventing of fires therein; and for lessening the great charge of the trained bands. V. The way to make New-Haven in Sussex, fit to receive ships of burthen. VI. Seasonable discourses of the tin, iron, linnen, and woollen trades; with advantageous proposals for improving them all : Illustrated with seven large copper-plates / By Andrew Yarranton, gent. London : Printed, to be sold by Tho. Parkhurst ..., 1698.

Stronger Together: The Cross-Cultural Coalition to Stop a Fossil Fuel Export Terminal in the Salish Sea

Thesis (Master's)--University of Washington, 2016-06

Georeferenced standardized apparent wave power along the coastline of Cabrera National Park for analysis of essential habitats of littoral species

Wind- induced exposure is one of the major forces shaping the geomorphology and biota in coastal areas. The effect of wave exposure on littoral biota is well known in marine environments (Ekebon et al., 2003; Burrows et al., 2008). In the Cabrera Archipelago National Park wave exposure has demostrated to have an effect on the spatial distribution of different stages of E.marginatus (Alvarez et al., 2010). Standarized average wave exposures during 2008 along the Cabrera Archipelago National park coast line were calculated to be applied in studies of littoral species distribution within the archipelago. Average wave exposure (or apparent wave power) was calculated for points located 50 m equidistant on the coastline following the EXA methodology (EXposure estimates for fragmented Archipelagos) (Ekebon et al., 2003). The average wave exposures were standardized from 1 to 100 (minimum and maximum in the area), showing coastal areas with different levels of mea wave exposure during the year. Input wind data (direction and intensity) from 2008 was registered at the Cabrera mooring located north of Cabrera Archipelago. Data were provided by IMEDEA (CSIC-UIB, TMMOS http://www.imedea.uib-csic.es/tmoos/boyas/). This cartography has been developed under the framework of the project EPIMHAR, funded by the National Park's Network (Spanish Ministry of Environment, Maritime and Rural Affairs, reference: 012/2007 ). Part of this work has been developed under the research programs funded by "Fons de Garantia Agrària i Pesquera de les Illes Balears (FOGAIBA)".

Monthly Arctic sea ice lead fraction in 0.5° x 0.5° resolution from SARAL/Altika altimeter, link to datasets in NetCDF 4 format

Sea ice leads play an essential role in ocean-ice-atmosphere exchange, in ocean circulation, geochemistry, and in ice dynamics. Their precise detection is crucial for altimetric estimations of sea ice thickness and volume. This study evaluates the performance of the SARAL/AltiKa (Satellite with ARgos and ALtiKa) altimeter to detect leads and to monitor their spatio-temporal dynamics. We show that a pulse peakiness parameter (PP) used to detect leads by Envisat RA-2 and ERS-1,-2 altimeters is not suitable because of saturation of AltiKa return echoes over the leads. The signal saturation results in loss of 6-10% of PP data over sea ice. We propose a different parameter-maximal power of waveform-and define the threshold to discriminate the leads. Our algorithm can be applied from December until May. It detects well the leads of small and medium size from 200 m to 3-4 km. So the combination of the high-resolution altimetric estimates with low-resolution thermal infra-red or radiometric lead fraction products could enhance the capability of remote sensing to monitor sea ice fracturing.

Simple heterogeneity parametrization for sea surface temperature and chlorophyll

Using satellite maps this paper offers a complex analysis of chlorophyll & SST heterogeneity in the shelf seas around the southwest of the UK. The heterogeneity scaling follows a simple power law and is consequently pa- rametrized by two parameters. It is shown that in most cases these two parameters vary only relatively little with time. The paper offers a detailed comparison of field heterogeneity between different regions. How much hetero- geneity is in each region preserved in the annual median data is also determined. The paper explicitly demon- strates how one can use these results to calculate representative measurement area for in situ networks.

Simple heterogeneity parametrization for sea surface temperature and chlorophyll

Using satellite maps this paper offers a complex analysis of chlorophyll & SST heterogeneity in the shelf seas around the southwest of the UK. The heterogeneity scaling follows a simple power law and is consequently pa- rametrized by two parameters. It is shown that in most cases these two parameters vary only relatively little with time. The paper offers a detailed comparison of field heterogeneity between different regions. How much hetero- geneity is in each region preserved in the annual median data is also determined. The paper explicitly demon- strates how one can use these results to calculate representative measurement area for in situ networks.

Designing Short Term Wave Traces to Assess Wave Power Devices

In recent years modern numerical methods have been employed in the design of Wave Energy Converters (WECs), however the high computational costs associated with their use makes it prohibitive to undertake simulations involving statistically relevant numbers of wave cycles. Experimental tests in wave tanks could also be performed more efficiently and economically if short time traces, consisting of only a few wave cycles, could be used to evaluate the hydrodynamic characteristics of a particular device or design modification. Ideally, accurate estimations of device performance could be made utilizing results obtained from investigations with a relatively small number of wave cycles. However the difficulty here is that many WECs, such as the Oscillating Wave Surge Converter (OWSC), exhibit significant non-linearity in their response. Thus it is challenging to make accurate predictions of annual energy yield for a given spectral sea state using short duration realisations of that sea. This is because the non-linear device response to particular phase couplings of sinusoidal components within those time traces might influence the estimate of mean power capture obtained. As a result it is generally accepted that the most appropriate estimate of mean power capture for a sea state be obtained over many hundreds (or thousands) of wave cycles. This ensures that the potential influence of phase locking is negligible in comparison to the predictions made. In this paper, potential methods of providing reasonable estimates of relative variations in device performance using short duration sea states are introduced. The aim of the work is to establish the shortness of sea state required to provide statistically significant estimations of the mean power capture of a particular type of Wave Energy Converter. The results show that carefully selected wave traces can be used to reliably assess variations in power output due to changes in the hydrodynamic design or wave climate. 

Panoramic Sea Happening (After Kantor): Video Installation

Panoramic Sea Happening (After Kantor) is a 7 minute durational film that reimagines part of Tadeusz Kantor's original sea happenings from 1967 in a landscape in which the sea has retreated. The conductor of Kantor’s original performance is replaced with a sound object cast adrift on a beach in Dungeness (UK). The object plays back the sound of the sea into the landscape, which was performed live and then filmed from three distinct angles. The first angle mimics the position of the conductor in Kantor’s original happening, facing outwards into the horizon of the beach and recalls the image in Kantor’s work of a human figure undertaking the absurd task of orchestrating the sound of a gigantic expanse of water. The second angle exposes the machine itself and the large cone that amplifies the sound, reinforcing the isolation of the object. The third angle reveals a decommissioned nuclear power station and sound objects used as a warning system for the power plant. Dungeness is a location where the sea has been retreating from the land, leaving traces of human activity through the disused boat winches, abandoned cabins and the decommissioned nuclear buildings. It is a place in which the footprint of the anthropocene is keenly felt. The sound object is intended to act as an anthropomorphic figure, ghosting the original conductor and offering the sound of the sea back into the landscape through a wide mouthpiece, echoing Kantor’s own load hailer in the original sequence of sea happenings. It speculates on Kantor's theory of the bio-object, which proposed a symbiotic relationship between the human and the nonhuman object in performance, as a possible instrument to access a form of geologic imagination. In this configuration, the human itself is absent, but is evoked through the objects left behind. The sound object, helpless in a red dingy, might be thought of as a co-conspirator with the viewer, enabling a looking back to the past in a landscape of an inevitable future. The work was originally commissioned by the University of Kent in collaboration with the Polish Cultural Institute for the Symposium Kantorbury Kantorbury in Canterbury (UK) to mark the 100 years since Tadeusz Kantor’s birth (15 - 19 September 2015). It should be projected and requires stereo speakers.

The three dimensional model with the affection of stratified density (Caspian Sea)

The hydro dynamical actions in big Lakes directly influence dynamic, physical and chemical affairs. The circulation's models and temperature have something to do with the movements of fluids, and analysis for circulation in Caspian sea is because of the lack of observation through which the circulations and out comings are determined. Through the studies, three dimensional simulations (Large- Scale) are planned and performed, according to Smolakiewicz and Margolin works. This is a non- hydrostatic and Boussinesq approximation is used in its formulation is used in its formulation on the basis of Lipps (1990) theorem and curve lines, the fluid is constant adiabatic and stratified, and the wind power is considered zero. The profile of speed according to previous depth and before ridge can be drawn on the basis of density available between northern and southern ridges. The circulation field is drawn from 3 cm/s to 13 cm/s on the plate z= 5 cm , the vertical changes of speed on the plate is 0.02 m/s. Vertical profile , horizontal speed in previous on, and after the ridges on are drawn on different spaces. It changes from 0.5 cm/s to 1 cm/s before ridges.

Concentrating solar power + desalination plants (CSP+D): models and performance analysis

Tese submetida à Universidade de Lisboa, Instituto Superior Técnico e aprovada em provas públicas para a obtenção do Grau de Doutor em Sistemas Sustentáveis de Energia.