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.

Characteristics of the turbulent/nonturbulent interface in boundary layers, jets and shear-free turbulence

The characteristics of turbulent/nonturbulent interfaces (TNTI) from boundary layers, jets and shear-free turbulence are compared using direct numerical simulations. The TNTI location is detected by assessing the volume of turbulent flow as function of the vorticity magnitude and is shown to be equivalent to other procedures using a scalar field. Vorticity maps show that the boundary layer contains a larger range of scales at the interface than in jets and shear-free turbulence where the change in vorticity characteristics across the TNTI is much more dramatic. The intermittency parameter shows that the extent of the intermittency region for jets and boundary layers is similar and is much bigger than in shear-free turbulence, and can be used to compute the vorticity threshold defining the TNTI location. The statistics of the vorticity jump across the TNTI exhibit the imprint of a large range of scales, from the Kolmogorov micro-scale to scales much bigger than the Taylor scale. Finally, it is shown that contrary to the classical view, the low-vorticity spots inside the jet are statistically similar to isotropic turbulence, suggesting that engulfing pockets simply do not exist in jets

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

Low-resolution spectroscopy and spectral energy distributions of selected sources towards σ Orionis

Aims. We study in detail nine sources in the direction of the young σ Orionis cluster, which is considered to be a unique site for studying stellar and substellar formation. The nine sources were selected because of their peculiar properties, such as extremely-red infrared colours or excessively strong Hα emission for their blue optical colours. Methods. We acquired high-quality, low-resolution spectroscopy (R ∼ 500) of the nine targets with ALFOSC at the Nordic Optical Telescope. We also re-analysed [24]-band photometry from MIPS/Spitzer and compiled the highest quality photometric dataset available at the ViJHK_s passbands and the four IRAC/Spitzer channels, for constructing accurate spectral energy distributions between 0.55 and 24 μm. Results. The nine targets were classified into: one Herbig Ae/Be star with a scattering edge-on disc; two G-type stars; one X-ray flaring, early-M, young star with chromospheric Hα emission; one very low-mass, accreting, young spectroscopic binary; two young objects at the brown-dwarf boundary with the characteristics of classical T Tauri stars; and two emission-line galaxies, one undergoing star formation, and another whose spectral energy distribution is dominated by an active galactic nucleus. We also discovered three infrared sources associated with overdensities in a cold cloud of the cluster centre. Conclusions. Low-resolution spectroscopy and spectral energy distributions are a vital tool for measuring the physical properties and evolution of young stars and candidates in the σ Orionis cluster.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.

Tropically-driven climate shifts in southwestern Europe during MIS 19, a low eccentricity interglacial

The relative roles of high- versus low-latitude forcing of millennial-scale climate variability are still not well understood. Here we present terrestrial–marine climate profiles from the southwestern Iberian margin, a region particularly affected by precession, that show millennial climate oscillations related to a nonlinear response to the Earth's precession cycle during Marine Isotope Stage (MIS) 19. MIS 19 has been considered the best analogue to our present interglacial from an astronomical point of view due to the reduced eccentricity centred at 785 ka. In our records, seven millennial-scale forest contractions punctuated MIS 19 superimposed to two orbitally-driven Mediterranean forest expansions. In contrast to our present interglacial, we evidence for the first time low latitude-driven 5000-yr cycles of drying and cooling in the western Mediterranean region, along with warmth in the subtropical gyre related to the fourth harmonic of precession. These cycles indicate repeated intensification of North Atlantic meridional moisture transport that along with decrease in boreal summer insolation triggered ice growth and may have contributed to the glacial inception, at ∼774 ka. The freshwater fluxes during MIS 19ab amplified the cooling events in the North Atlantic promoting further cooling and leading to MIS 18 glaciation. The discrepancy between the dominant cyclicity observed during MIS 1, 2500-yr, and that of MIS 19, 5000-yr, challenges the similar duration of the Holocene and MIS 19c interglacials under natural boundary conditions.

Benthic foraminiferal abundances at the Cretaceous/Paleogene boundary of ODP Hole 117-1049C at Blake Nose, Northwest Atlantic (Appendix A)

Sediments recovered at lower bathyal ODP Site 1049 on Blake Nose (Northwestern Atlantic) offer an opportunity to study environmental changes at the Cretaceous/Paleogene (K/P) boundary relatively close to the Chicxulub impact structure on the Yucatan peninsula, Mexico. In Hole 1049C, the boundary is located at the base of a 9-cm-thick layer with abundant spherules, considered to be impact ejecta. Uppermost Maastrichtian oozes below, and lowermost Danian pelagic oozes above the spherulebed contain well-preserved bathyal benthic foraminifera. The spherule-bed itself, in contrast, contains a mixture of shallow (neritic) and deeper (bathyal) species, and specimens vary strongly in preservation. This assemblage was probably formed by reworking and down-slope transport triggered by the K/P impact. Across the spherule-bed (i.e., the K/P boundary) only ~7% of benthic foraminiferal species became extinct, similar to the low extinction rates of benthic foraminifera worldwide. Quantitative analysis of benthic foraminiferal assemblages and morphogroups in the >63-µm size fraction indicates a relatively eutrophic, stable environment during the latest Maastrichtian, interrupted by a sudden decrease in the food supply to the benthos at the K/P boundary and a decrease in diversity of the faunas, followed by a stepped recovery during the earliest Danian. The recovery was probably linked to the gradual recovery of surface-dwelling primary producers.

Abundance of calcareous nannofossils across the K/T boundary at ODP Hole 113-690C (Table 1)

A biostratigraphically continuous, but intensely bioturbated, Cretaceous/Tertiary boundary sequence was cored during Ocean Drilling Program (ODP) Leg 113 on Maud Rise (65°S) in the Weddell Sea off East Antarctica. This interval is the first recovered by ODP/DSDP in the South Atlantic sector of the Southern Ocean and offers a unique opportunity to study the nannofossil sequences leading up to and beyond the terminal Cretaceous event at a high southern latitude. The K/T boundary lies just within Chron 29R and is placed at ODP Sample 113-690C-15X-4, 41.5 cm. An iridium anomaly was independently noted at about this level as well. Upper Maestrichtian-lower Paleocene sediments consist mostly of light-colored nannofossil chalks. Dark brown sediments at the base of the Danian (Zone CPla) are characterized by an increased clay content attributed to a drop in calcareous microplankton productivity following the terminal Cretaceous event. Although delineation of the boundary is hampered by intense bioturbation, the sharp color contrast between overlying clay-rich, dark brown chalks of the Tertiary and light cream colored chalks of the Cretaceous aids in the selection of the K/T horizon. Several dark colored burrows sampled at intervals as far as 1.3 m below the boundary and within the light colored Cretaceous chalk were found to contain up to 17% Tertiary nannofossils. Calcareous nannofossils from the boundary interval were divided into three groups for quantitative study. The three groups, "Cretaceous," "Tertiary," and "Survivor," exhibit a sequential change across the boundary with the Cretaceous forms giving way to a Survivor-dominated assemblage beginning at the boundary followed shortly thereafter by the appearance of the Tertiary taxa, Cruciplacolithus and Hornibrookina. The species, H. edwardsii, comprises nearly 50% of the assemblage just above the Zone CPla/CPlb boundary, an abundance not reported elsewhere at this level. Calculation of individual species abundances reveals several additional differences between this K/T boundary interval and those studied from middle and low latitude sections. The percentage of Thoracosphaera is much lower at the boundary in this section and a small form, Prediscosphaera stoveri, is extremely abundant in Cretaceous sediments just below the boundary.

(Table 1) Conditions of vertical near-bottom profiling and the main parameters of the benthic boundary layer in several regions of the Black Sea shelf

Measurements of 14 vertical profiles of currents and hydrological parameters in the near-bottom layer with depth resolution of 0.1 m were carried out in several regions of the Black Sea shelf, at five points over the continental slope, and in three deep water regions. The upper boundary of the benthic boundary layer (BBL) was reliably determined at a point at distance from 5-7 to 35-40 m from the bottom where the gradients of density and current velocity changed. Experimental data obtained were used to determine the coefficient of bottom friction, friction velocity, coefficients of vertical diffusion of momentum and density, and vertical fluxes of temperature and salinity in the BBL.

Concentrations and isotopic compositions of low molecular weight hydrocarbons in hydrothermal fluids from ODP Field, Middle Valley, northern Juan de Fuca Ridge

The presence of sedimentary organic matter blanketing midocean ridge crests has a potentially strong impact on metal transport in hydrothermal vent fluids. To constrain the role of organic matter in metal mobility during hydrothermal sediment alteration, we reacted organic-rich diatomaceous ooze from Guaymas Basin, Gulf of California, and organic-poor hemipelagic mud from Middle Valley, northern Juan de Fuca Ridge, with seawater and a Na-Ca-K-Cl fluid of seawater chlorinity, at 275° to 400°C, 350 to 500 bars, and initial fluid: sediment mass ratios ranging from 1.6 to 9.8. Reaction of these fluids with both sediment types released CO2 and high concentrations of ore-forming metals (Fe, Mn, Zn, Pb) to solution. Relatively low concentrations of Cu were observed in solution and likely reflect the reducing conditions that resulted from the presence of sedimentary organic matter. Both the concentrations of CO2 and dissolved metals were lower in fluids reacted with Middle Valley sediment compared with aqueous concentrations in fluids reacted with Guaymas Basin sediment. During alteration of both sediment types, metal concentrations varied strongly as a function of temperature, increasing by up to an order of magnitude over the 75°C range of each experiment. Major element fluid chemistry and observed alteration assemblages suggest that during hydrothermal alteration of organic-lean sediment from Middle Valley a feldspar-quartz-illite mineral assemblage buffered in situ pH. In contrast, data from the experimental alteration of organic-rich Guaymas Basin sediment suggest that a calcite-plagioclase-quartz assemblage regulated in situ pH. Fluid speciation calculations suggest that in situ pH during Guaymas Basin sediment alteration was lower than during alteration of Middle Valley sediment and accounts for the substantially greater metal mobility at a given temperature and pressure during the former experiment. Comparison of our results with the results of basalt alteration experiments indicate that except for Cu, hydrothermal sediment alteration results in equal or greater concentrations of ore-forming metals at a given temperature and pressure. Accordingly, the presence of ore-forming metals in fluids currently venting from sediment-covered hydrothermal systems at concentrations substantially lower than in fluids from bare-rock systems may reflect chemical reequilibration during subsurface cooling within the sediment pile.

Selection of species and provenances for low-rainfall areas: physiological responses of Eucalyptus cloeziana and Eucalyptus argophloia to seasonal conditions in subtropical Queensland

Responses of stomatal conductance (g(s)) and net photosynthesis (A) to changes in soil water availability, photosynthetic photon flux density (Q), air temperature (1) and leaf-to-air vapour pressure deficit (D) were investigated in 4-year-old trees of a dry inland provenance of Eucalyptus argophloia Blakely, and two dry inland provenances (Coominglah and Hungry Hills) and a humid coastal provenance (Wolvi) of Eucalyptus cloeziana F. Muell. between April 2001 and April 2002 in southeast Queensland, Australia. There were minimal differences in A, g, and water relations variables among the coastal and inland provenances of E. cloeziana but large differences between E. argophloia and E. cloeziana. E. argophloia and to a lesser extent the Hungry Hills (inland) provenance of E. cloeziana maintained relatively higher pre-dawn water potential (psi(pd)) during the dry season suggesting possible access to water at depth. Simple phenomenological models of stomatal conductance as a function of Q, T and D explained 60% of variation in gs in E. cloeziana and more than 75% in E. argophloia, when seasonal effect was incorporated in the model. A Ball-Berry model for net photosynthesis explained between 70 and 80% of observed variation in A in both species. These results have implications in matching the dry and humid provenances of E. cloeziana and E. argophloia to suitable sites in subtropical environments. (C) 2004 Elsevier B.V. All rights reserved.

Multigradient field-active contour model for multilayer boundary detection of ultrasound rectal wall image

Extraction and reconstruction of rectal wall structures from an ultrasound image is helpful for surgeons in rectal clinical diagnosis and 3-D reconstruction of rectal structures from ultrasound images. The primary task is to extract the boundary of the muscular layers on the rectal wall. However, due to the low SNR from ultrasound imaging and the thin muscular layer structure of the rectum, this boundary detection task remains a challenge. An active contour model is an effective high-level model, which has been used successfully to aid the tasks of object representation and recognition in many image-processing applications. We present a novel multigradient field active contour algorithm with an extended ability for multiple-object detection, which overcomes some limitations of ordinary active contour models—"snakes." The core part in the algorithm is the proposal of multigradient vector fields, which are used to replace image forces in kinetic function for alternative constraints on the deformation of active contour, thereby partially solving the initialization limitation of active contour for rectal wall boundary detection. An adaptive expanding force is also added to the model to help the active contour go through the homogenous region in the image. The efficacy of the model is explained and tested on the boundary detection of a ring-shaped image, a synthetic image, and an ultrasound image. The experimental results show that the proposed multigradient field-active contour is feasible for multilayer boundary detection of rectal wall

Sheet flow sediment transport under waves with acceleration skewness and boundary layer streaming

The effect of acceleration skewness on sheet flow sediment transport rates (q) over bar (s) is analysed using new data which have acceleration skewness and superimposed currents but no boundary layer streaming. Sediment mobilizing forces due to drag and to acceleration (similar to pressure gradients) are weighted by cosine and sine, respectively, of the angle phi(.)(tau)phi(tau) = 0 thus corresponds to drag dominated sediment transport, (q) over bar (s)similar to vertical bar u(infinity)vertical bar u(infinity), while phi(tau) = 90 degrees corresponds to total domination by the pressure gradients, (q) over bar similar to du(infinity)/dt. Using the optimal angle, phi = 51 degrees based on that data, good agreement is subsequently found with data that have strong influence from boundary layer streaming. Good agreement is also maintained with the large body of U-tube data simulating sine waves with superimposed currents and second-order Stokes waves, all of which have zero acceleration skewness. The recommended model can be applied to irregular waves with arbitrary shape as long as the assumption negligible time lag between forcing and sediment transport rate is valid. With respect to irregular waves, the model is much easier to apply than the competing wave-by-wave models. Issues for further model developments are identified through a comprehensive data review.

Dinosaur sanctuary on the Chatham Islands, Southwest Pacific: First record of theropods from the K-T boundary Takatika Grit

Cretaceous-Tertiaty (K-T) boundary (ca. 65 Ma) sections on a Southwest Pacific island containing dinosaurs were unknown until March 2003 when theropod bones were recovered from the Takatika Grit on the remote Chatham Islands (latitude 44 degrees S, longitude 176 degrees W), along the Chatham Rise. Tectonic and palaeontologic evidence support the eastward extension of a ca. 900 km land bridge that connected the islands to what is now New Zealand prior to the K-T boundary. The Chathams terrestrial fauna inhabited coastal, temperate environments along a low-lying, narrow, crustal extension of the New Zealand subcontinent, characterised by a tectonically dynamic, volcanic landscape with eroding hills (horsts) adjacent to flood plains and deltas, all sediments accumulating in grabens. This finger-like tract was blanketed with a conifer and clubmoss (Lycopodiopsida) dominated forest. The Chatham Islands region would have, along with New Zealand, provided a dinosaur island sanctuary after separating from the Gondwana margin ca. 80 Ma. (c) 2005 Elsevier B.V. All rights reserved.

Experimental and analytical study of the effect of contact angle on liquid convective heat transfer in microchannels

In this paper we examine the effect of contact angle (or surface wettability) on the convective heat transfer coefficient in microchannels. Slip flow, where the fluid velocity at the wall is non-zero, is most likely to occur in microchannels due to its dependence on shear rate or wall shear stress. We show analytically that for a constant pressure drop, the presence of slip increases the Nusselt number. In a microchannel heat exchanger we modified the surface wettability from a contact angle of 20 degrees-120 degrees using thin film coating technology. Apparent slip flow is implied from pressure and flow rate measurements with a departure from classical laminar friction coefficients above a critical shear rate of approximately 10,000 s(-1). The magnitude of this departure is dependant on the contact angle with higher contact angles surfaces exhibiting larger pressure drop decreases. Similarly, the non-dimensional heat flux is found to decrease relative to laminar non-slip theory, and this decrease is also a function of the contact angle. Depending on the contact angle and the wall shear rate, variations in the heat transfer rate exceeding 10% can be expected. Thus the contact angle is an important consideration in the design of micro, and even more so, nano heat exchangers. (c) 2006 Elsevier Ltd. All rights reserved.