Monte Carlo Neutronics and Thermal-hydraulics coupling applied to Fast Reactors

Un escenario habitualmente considerado para el uso sostenible y prolongado de la energía nuclear contempla un parque de reactores rápidos refrigerados por metales líquidos (LMFR) dedicados al reciclado de Pu y la transmutación de actínidos minoritarios (MA). Otra opción es combinar dichos reactores con algunos sistemas subcríticos asistidos por acelerador (ADS), exclusivamente destinados a la eliminación de MA. El diseño y licenciamiento de estos reactores innovadores requiere herramientas computacionales prácticas y precisas, que incorporen el conocimiento obtenido en la investigación experimental de nuevas configuraciones de reactores, materiales y sistemas. A pesar de que se han construido y operado un cierto número de reactores rápidos a nivel mundial, la experiencia operacional es todavía reducida y no todos los transitorios se han podido entender completamente. Por tanto, los análisis de seguridad de nuevos LMFR están basados fundamentalmente en métodos deterministas, al contrario que las aproximaciones modernas para reactores de agua ligera (LWR), que se benefician también de los métodos probabilistas. La aproximación más usada en los estudios de seguridad de LMFR es utilizar una variedad de códigos, desarrollados a base de distintas teorías, en busca de soluciones integrales para los transitorios e incluyendo incertidumbres. En este marco, los nuevos códigos para cálculos de mejor estimación ("best estimate") que no incluyen aproximaciones conservadoras, son de una importancia primordial para analizar estacionarios y transitorios en reactores rápidos. Esta tesis se centra en el desarrollo de un código acoplado para realizar análisis realistas en reactores rápidos críticos aplicando el método de Monte Carlo. Hoy en día, dado el mayor potencial de recursos computacionales, los códigos de transporte neutrónico por Monte Carlo se pueden usar de manera práctica para realizar cálculos detallados de núcleos completos, incluso de elevada heterogeneidad material. Además, los códigos de Monte Carlo se toman normalmente como referencia para los códigos deterministas de difusión en multigrupos en aplicaciones con reactores rápidos, porque usan secciones eficaces punto a punto, un modelo geométrico exacto y tienen en cuenta intrínsecamente la dependencia angular de flujo. En esta tesis se presenta una metodología de acoplamiento entre el conocido código MCNP, que calcula la generación de potencia en el reactor, y el código de termohidráulica de subcanal COBRA-IV, que obtiene las distribuciones de temperatura y densidad en el sistema. COBRA-IV es un código apropiado para aplicaciones en reactores rápidos ya que ha sido validado con resultados experimentales en haces de barras con sodio, incluyendo las correlaciones más apropiadas para metales líquidos. En una primera fase de la tesis, ambos códigos se han acoplado en estado estacionario utilizando un método iterativo con intercambio de archivos externos. El principal problema en el acoplamiento neutrónico y termohidráulico en estacionario con códigos de Monte Carlo es la manipulación de las secciones eficaces para tener en cuenta el ensanchamiento Doppler cuando la temperatura del combustible aumenta. Entre todas las opciones disponibles, en esta tesis se ha escogido la aproximación de pseudo materiales, y se ha comprobado que proporciona resultados aceptables en su aplicación con reactores rápidos. Por otro lado, los cambios geométricos originados por grandes gradientes de temperatura en el núcleo de reactores rápidos resultan importantes para la neutrónica como consecuencia del elevado recorrido libre medio del neutrón en estos sistemas. Por tanto, se ha desarrollado un módulo adicional que simula la geometría del reactor en caliente y permite estimar la reactividad debido a la expansión del núcleo en un transitorio. éste módulo calcula automáticamente la longitud del combustible, el radio de la vaina, la separación de los elementos de combustible y el radio de la placa soporte en función de la temperatura. éste efecto es muy relevante en transitorios sin inserción de bancos de parada. También relacionado con los cambios geométricos, se ha implementado una herramienta que, automatiza el movimiento de las barras de control en busca d la criticidad del reactor, o bien calcula el valor de inserción axial las barras de control. Una segunda fase en la plataforma de cálculo que se ha desarrollado es la simulació dinámica. Puesto que MCNP sólo realiza cálculos estacionarios para sistemas críticos o supercríticos, la solución más directa que se propone sin modificar el código fuente de MCNP es usar la aproximación de factorización de flujo, que resuelve por separado la forma del flujo y la amplitud. En este caso se han estudiado en profundidad dos aproximaciones: adiabática y quasiestática. El método adiabático usa un esquema de acoplamiento que alterna en el tiempo los cálculos neutrónicos y termohidráulicos. MCNP calcula el modo fundamental de la distribución de neutrones y la reactividad al final de cada paso de tiempo, y COBRA-IV calcula las propiedades térmicas en el punto intermedio de los pasos de tiempo. La evolución de la amplitud de flujo se calcula resolviendo las ecuaciones de cinética puntual. Este método calcula la reactividad estática en cada paso de tiempo que, en general, difiere de la reactividad dinámica que se obtendría con la distribución de flujo exacta y dependiente de tiempo. No obstante, para entornos no excesivamente alejados de la criticidad ambas reactividades son similares y el método conduce a resultados prácticos aceptables. Siguiendo esta línea, se ha desarrollado después un método mejorado para intentar tener en cuenta el efecto de la fuente de neutrones retardados en la evolución de la forma del flujo durante el transitorio. El esquema consiste en realizar un cálculo cuasiestacionario por cada paso de tiempo con MCNP. La simulación cuasiestacionaria se basa EN la aproximación de fuente constante de neutrones retardados, y consiste en dar un determinado peso o importancia a cada ciclo computacial del cálculo de criticidad con MCNP para la estimación del flujo final. Ambos métodos se han verificado tomando como referencia los resultados del código de difusión COBAYA3 frente a un ejercicio común y suficientemente significativo. Finalmente, con objeto de demostrar la posibilidad de uso práctico del código, se ha simulado un transitorio en el concepto de reactor crítico en fase de diseño MYRRHA/FASTEF, de 100 MW de potencia térmica y refrigerado por plomo-bismuto. ABSTRACT Long term sustainable nuclear energy scenarios envisage a fleet of Liquid Metal Fast Reactors (LMFR) for the Pu recycling and minor actinides (MAs) transmutation or combined with some accelerator driven systems (ADS) just for MAs elimination. Design and licensing of these innovative reactor concepts require accurate computational tools, implementing the knowledge obtained in experimental research for new reactor configurations, materials and associated systems. Although a number of fast reactor systems have already been built, the operational experience is still reduced, especially for lead reactors, and not all the transients are fully understood. The safety analysis approach for LMFR is therefore based only on deterministic methods, different from modern approach for Light Water Reactors (LWR) which also benefit from probabilistic methods. Usually, the approach adopted in LMFR safety assessments is to employ a variety of codes, somewhat different for the each other, to analyze transients looking for a comprehensive solution and including uncertainties. In this frame, new best estimate simulation codes are of prime importance in order to analyze fast reactors steady state and transients. This thesis is focused on the development of a coupled code system for best estimate analysis in fast critical reactor. Currently due to the increase in the computational resources, Monte Carlo methods for neutrons transport can be used for detailed full core calculations. Furthermore, Monte Carlo codes are usually taken as reference for deterministic diffusion multigroups codes in fast reactors applications because they employ point-wise cross sections in an exact geometry model and intrinsically account for directional dependence of the ux. The coupling methodology presented here uses MCNP to calculate the power deposition within the reactor. The subchannel code COBRA-IV calculates the temperature and density distribution within the reactor. COBRA-IV is suitable for fast reactors applications because it has been validated against experimental results in sodium rod bundles. The proper correlations for liquid metal applications have been added to the thermal-hydraulics program. Both codes are coupled at steady state using an iterative method and external files exchange. The main issue in the Monte Carlo/thermal-hydraulics steady state coupling is the cross section handling to take into account Doppler broadening when temperature rises. Among every available options, the pseudo materials approach has been chosen in this thesis. This approach obtains reasonable results in fast reactor applications. Furthermore, geometrical changes caused by large temperature gradients in the core, are of major importance in fast reactor due to the large neutron mean free path. An additional module has therefore been included in order to simulate the reactor geometry in hot state or to estimate the reactivity due to core expansion in a transient. The module automatically calculates the fuel length, cladding radius, fuel assembly pitch and diagrid radius with the temperature. This effect will be crucial in some unprotected transients. Also related to geometrical changes, an automatic control rod movement feature has been implemented in order to achieve a just critical reactor or to calculate control rod worth. A step forward in the coupling platform is the dynamic simulation. Since MCNP performs only steady state calculations for critical systems, the more straight forward option without modifying MCNP source code, is to use the flux factorization approach solving separately the flux shape and amplitude. In this thesis two options have been studied to tackle time dependent neutronic simulations using a Monte Carlo code: adiabatic and quasistatic methods. The adiabatic methods uses a staggered time coupling scheme for the time advance of neutronics and the thermal-hydraulics calculations. MCNP computes the fundamental mode of the neutron flux distribution and the reactivity at the end of each time step and COBRA-IV the thermal properties at half of the the time steps. To calculate the flux amplitude evolution a solver of the point kinetics equations is used. This method calculates the static reactivity in each time step that in general is different from the dynamic reactivity calculated with the exact flux distribution. Nevertheless, for close to critical situations, both reactivities are similar and the method leads to acceptable practical results. In this line, an improved method as an attempt to take into account the effect of delayed neutron source in the transient flux shape evolutions is developed. The scheme performs a quasistationary calculation per time step with MCNP. This quasistationary simulations is based con the constant delayed source approach, taking into account the importance of each criticality cycle in the final flux estimation. Both adiabatic and quasistatic methods have been verified against the diffusion code COBAYA3, using a theoretical kinetic exercise. Finally, a transient in a critical 100 MWth lead-bismuth-eutectic reactor concept is analyzed using the adiabatic method as an application example in a real system.

Magnetic Attitude Control System for a Small Satellite. Impact on the Thermal Performance

El principal objetivo de la tesis es estudiar el acoplamiento entre los subsistemas de control de actitud y de control térmico de un pequeño satélite, con el fin de buscar la solución a los problemas relacionados con la determinación de los parámetros de diseño. Se considera la evolución de la actitud y de las temperaturas del satélite bajo la influencia de dos estrategias de orientación diferentes: 1) estabilización magnética pasiva de la orientación (PMAS, passive magnetic attitude stabilization), y 2) control de actitud magnético activo (AMAC, active magnetic attitude control). En primer lugar se presenta el modelo matemático del problema, que incluye la dinámica rotacional y el modelo térmico. En el problema térmico se considera un satélite cúbico modelizado por medio de siete nodos (seis externos y uno interno) aplicando la ecuación del balance térmico. Una vez establecido el modelo matemático del problema, se estudia la evolución que corresponde a las dos estrategias mencionadas. La estrategia PMAS se ha seleccionado por su simplicidad, fiabilidad, bajo coste, ahorrando consumo de potencia, masa coste y complejidad, comparado con otras estrategias. Se ha considerado otra estrategia de control que consigue que el satélite gire a una velocidad requerida alrededor de un eje deseado de giro, pudiendo controlar su dirección en un sistema inercial de referencia, ya que frecuentemente el subsistema térmico establece requisitos de giro alrededor de un eje del satélite orientado en una dirección perpendicular a la radiación solar incidente. En relación con el problema térmico, para estudiar la influencia de la velocidad de giro en la evolución de las temperaturas en diversos puntos del satélite, se ha empleado un modelo térmico linealizado, obtenido a partir de la formulación no lineal aplicando un método de perturbaciones. El resultado del estudio muestra que el tiempo de estabilización de la temperatura y la influencia de las cargas periódicas externas disminuye cuando aumenta la velocidad de giro. Los cambios de temperatura se reducen hasta ser muy pequeños para velocidades de rotación altas. En relación con la estrategia PMAC se ha observado que a pesar de su uso extendido entre los micro y nano satélites todavía presenta problemas que resolver. Estos problemas están relacionados con el dimensionamiento de los parámetros del sistema y la predicción del funcionamiento en órbita. Los problemas aparecen debido a la dificultad en la determinación de las características magnéticas de los cuerpos ferromagnéticos (varillas de histéresis) que se utilizan como amortiguadores de oscilaciones en los satélites. Para estudiar este problema se presenta un modelo analítico que permite estimar la eficiencia del amortiguamiento, y que se ha aplicado al estudio del comportamiento en vuelo de varios satélites, y que se ha empleado para comparar los resultados del modelo con los obtenidos en vuelo, observándose que el modelo permite explicar satisfactoriamente el comportamiento registrado. ABSTRACT The main objective of this thesis is to study the coupling between the attitude control and thermal control subsystems of a small satellite, and address the solution to some existing issues concerning the determination of their parameters. Through the thesis the attitude and temperature evolution of the satellite is studied under the influence of two independent attitude stabilization and control strategies: (1) passive magnetic attitude stabilization (PMAS), and (2) active magnetic attitude control (AMAC). In this regard the mathematical model of the problem is explained and presented. The mathematical model includes both the rotational dynamics and the thermal model. The thermal model is derived for a cubic satellite by solving the heat balance equation for 6 external and 1 internal nodes. Once established the mathematical model of the problem, the above mentioned attitude strategies were applied to the system and the temperature evolution of the 7 nodes of the satellite was studied. The PMAS technique has been selected to be studied due to its prevalent use, simplicity, reliability, and cost, as this strategy significantly saves the overall power, weight, cost, and reduces the complexity of the system compared to other attitude control strategies. In addition to that, another control law that provides the satellite with a desired spin rate along a desired axis of the satellite, whose direction can be controlled with respect to the inertial reference frame is considered, as the thermal subsystem of a satellite usually demands a spin requirement around an axis of the satellite which is positioned perpendicular to the direction of the coming solar radiation. Concerning the thermal problem, to study the influence of spin rate on temperature evolution of the satellite a linear approach of the thermal model is used, which is based on perturbation theory applied to the nonlinear differential equations of the thermal model of a spacecraft moving in a closed orbit. The results of this study showed that the temperature stabilization time and the periodic influence of the external thermal loads decreases by increasing the spin rate. However, the changes become insignificant for higher values of spin rate. Concerning the PMAS strategy, it was observed that in spite of its extended application to micro and nano satellites, still there are some issues to be solved regarding this strategy. These issues are related to the sizing of its system parameters and predicting the in-orbit performance. The problems were found to be rooted in the difficulties that exist in determining the magnetic characteristics of the ferromagnetic bodies (hysteresis rods) that are applied as damping devices on-board satellites. To address these issues an analytic model for estimating their damping efficiency is proposed and applied to several existing satellites in order to compare the results with their respective in-flight data. This model can explain the behavior showed by these satellites.

Nonlinear analysis of a simple model of temperature evolution in a satellite

We analyze a simple model of the heat transfer to and from a small satellite orbiting round a solar system planet. Our approach considers the satellite isothermal, with external heat input from the environment and from internal energy dissipation, and output to the environment as black-body radiation. The resulting nonlinear ordinary differential equation for the satellite’s temperature is analyzed by qualitative, perturbation and numerical methods, which prove that the temperature approaches a periodic pattern (attracting limit cycle). This approach can occur in two ways, according to the values of the parameters: (i) a slow decay towards the limit cycle over a time longer than the period, or (ii) a fast decay towards the limit cycle over a time shorter than the period. In the first case, an exactly soluble average equation is valid. We discuss the consequences of our model for the thermal stability of satellites.

In vitro evolution of horse heart myoglobin to increase peroxidase activity

Random mutagenesis and screening for enzymatic activity has been used to engineer horse heart myoglobin to enhance its intrinsic peroxidase activity. A chemically synthesized gene encoding horse heart myoglobin was subjected to successive cycles of PCR random mutagenesis. The mutated myoglobin gene was expressed in Escherichia coli LE392, and the variants were screened for peroxidase activity with a plate assay. Four cycles of mutagenesis and screening produced a series of single, double, triple, and quadruple variants with enhanced peroxidase activity. Steady-state kinetics analysis demonstrated that the quadruple variant T39I/K45D/F46L/I107F exhibits peroxidase activity significantly greater than that of the wild-type protein with k1 (for H2O2 oxidation of metmyoglobin) of 1.34 × 104 M−1 s−1 (≈25-fold that of wild-type myoglobin) and k3 [for reducing the substrate (2, 2′-azino-di-(3-ethyl)benzthiazoline-6-sulfonic acid] of 1.4 × 106 M−1 s−1 (1.6-fold that of wild-type myoglobin). Thermal stability of these variants as measured with circular dichroism spectroscopy demonstrated that the Tm of the quadruple variant is decreased only slightly compared with wild-type (74.1°C vs. 76.5°C). The rate constants for binding of dioxygen exhibited by the quadruple variant are identical to the those observed for wild-type myoglobin (kon, 22.2 × 10−6 M−1 s−1 vs. 22.3 × 10−6 M−1 s−1; koff, 24.3 s−1 vs. 24.2 s−1; KO2, 0.91 × 10−6 M−1 vs. 0.92 × 10−6 M−1). The affinity of the quadruple variant for CO is increased slightly (kon, 0.90 × 10−6 M−1s−1 vs. 0.51 × 10−6 M−1s−1; koff, 5.08 s−1 vs. 3.51 s−1; KCO, 1.77 × 10−7 M−1 vs. 1.45 × 10−7 M−1). All four substitutions are in the heme pocket and within 5 Å of the heme group.

Evolution of genome–phenome diversity under environmental stress

The genomic era revolutionized evolutionary biology. The enigma of genotypic-phenotypic diversity and biodiversity evolution of genes, genomes, phenomes, and biomes, reviewed here, was central in the research program of the Institute of Evolution, University of Haifa, since 1975. We explored the following questions. (i) How much of the genomic and phenomic diversity in nature is adaptive and processed by natural selection? (ii) What is the origin and evolution of adaptation and speciation processes under spatiotemporal variables and stressful macrogeographic and microgeographic environments? We advanced ecological genetics into ecological genomics and analyzed globally ecological, demographic, and life history variables in 1,200 diverse species across life, thousands of populations, and tens of thousands of individuals tested mostly for allozyme and partly for DNA diversity. Likewise, we tested thermal, chemical, climatic, and biotic stresses in several model organisms. Recently, we introduced genetic maps and quantitative trait loci to elucidate the genetic basis of adaptation and speciation. The genome–phenome holistic model was deciphered by the global regressive, progressive, and convergent evolution of subterranean mammals. Our results indicate abundant genotypic and phenotypic diversity in nature. The organization and evolution of molecular and organismal diversity in nature at global, regional, and local scales are nonrandom and structured; display regularities across life; and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection, including diversifying, balancing, cyclical, and purifying selective regimes, interacting with, but ultimately overriding, the effects of mutation, migration, and stochasticity.

Semivolatile and Volatile Compound Evolution during Pyrolysis and Combustion of Colombian Coffee Husk

Thermal characterization of coffee husk (Coffea arabica) from Colombian coffee has been studied. Different products, mostly volatile and semivolatile compounds, were analyzed, paying special attention to 16 polycyclic aromatic hydrocarbons (PAHs) classified by the United States Environmental Protection Agency (U.S. EPA) as priority pollutants, frequently used for checking toxicity in environmental samples. A fixed amount of raw material was exposed to different excess air ratios (λ = 0–2.33) and nominal temperature of 1123 K in a horizontal quartz reactor. The results show that coffee husk is a promising biomass for energetic exploitation with reduced formation of PAHs in a low air excess ratio. This implies reduction of carcinogenic potential in the limited presence of oxygen, demonstrated by calculating the carcinogenic potential (KE) for each experimental condition. Most volatile and semivolatile compounds followed different trends, with the oxygen presence prevailing their decomposition with increasing the air excess ratio.

A falsification of the thermal specialization paradigm: compensation for elevated temperatures in Antarctic fishes

Specialization to a particular environment is one of the main factors used to explain species distributions. Antarctic fishes are often cited as a classic example to illustrate the specialization process and are regarded as the archetypal stenotherms. Here we show that the Antarctic fish Pagothenia borchgrevinki has retained the capacity to compensate for chronic temperature change. By displaying astounding plasticity in cardiovascular response and metabolic control, the fishes maintained locomotory performance at elevated temperatures. Our falsification of the specialization paradigm indicates that the effect of climate change on species distribution and extinction may be overestimated by current models of global warming.

Adaptive responses of wild mungbean (Vigna radiata ssp sublobata) to photo-thermal environment. II. Growth, biomass, and seed yield

The leaf growth, dry matter production, and seed yield of 11 wild mungbean ( Vigna radiata ssp. sublobata) accessions of diverse geographic origin were observed under natural and artificial photoperiod temperature conditions, to determine the extent to which genotypic differences could be attributed to adaptive responses to photo-thermal environment. Environments included serial sowings in the field in SE Queensland, complemented by artificial photoperiod extension and controlled-environment growth rooms. Photo-thermal environment influenced leaf growth, total dry matter production ( TDM), and seed yield directly, through effects of ( mainly cool) temperature on growth, and indirectly, through effects on phenology. In terms of direct effects, leaf production, leaf expansion, and leaf area were all sensitive to temperature, with implied base temperatures higher than usually observed in cultivated mungbean ( V. radiata ssp. radiata). Genotypic sensitivity to temperature varied systematically with accession provenance and appeared to be of adaptive significance. In terms of the indirect effects of photo-thermal environment, genotypic and environmental effects on TDM were positively related to changes in total growth duration, and harvest index was negatively related to the period from sowing to flowering, similar to cultivated mungbean. However, seed yield was positively related to the duration of reproductive growth, reflecting the indeterminate growth habit of the wild accessions. As a consequence, the wild accessions are more responsive to favourable environments than typically observed in cultivated mungbean, which is determinate in habit. It is suggested that the introduction of the indeterminate trait into mungbean from the wild subspecies would increase the responsiveness of mungbean to favourable environments, analogous to that of black gram ( V. mungo). Although the wild subspecies appeared more sensitive to cool temperature than cultivated mungbean, it may provide a source of tolerance to the warmer temperatures experienced during the wet season in the tropics.

Coadaptation: A unifying principle in evolutionary thermal biology

Over the last 50 yr, thermal biology has shifted from a largely physiological science to a more integrated science of behavior, physiology, ecology, and evolution. Today, the mechanisms that underlie responses to environmental temperature are being scrutinized at levels ranging from genes to organisms. From these investigations, a theory of thermal adaptation has emerged that describes the evolution of thermoregulation, thermal sensitivity, and thermal acclimation. We review and integrate current models to form a conceptual model of coadaptation. We argue that major advances will require a quantitative theory of coadaptation that predicts which strategies should evolve in specific thermal environments. Simply combining current models, however, is insufficient to understand the responses of organisms to thermal heterogeneity; a theory of coadaptation must also consider the biotic interactions that influence the net benefits of behavioral and physiological strategies. Such a theory will be challenging to develop because each organism's perception of and response to thermal heterogeneity depends on its size, mobility, and life span. Despite the challenges facing thermal biologists, we have never been more pressed to explain the diversity of strategies that organisms use to cope with thermal heterogeneity and to predict the consequences of thermal change for the diversity of communities.

Effect of thermal acclimation on female resistance to forced matings in the eastern mosquitofish

All copulations in the eastern mosquitofish, Gambusia holbrooki, are coercive-and-achieved by force. Female G. holbrooki never appear to cooperate with males, but vigorously resist matings at all times. We examined the role of females within a sexually coercive mating system by investigating the ability of female G. holbrooki to resist forced copulations after acclimation to 16 degrees C and 32 degrees C for 4-5 weeks. We also examined burst swimming performance of female G. holbrooki after acclimation, as this trait is likely to underlie a female's ability to resist forced matings. We predicted that if female G. holbrooki indiscriminately resist matings from all males, acclimation would enhance female resistance at their acclimation temperature. However, we found that it did not. We also predicted that if females are able to influence the outcome of mating interactions, acclimation to an optimal thermal environment may induce females to reduce resistance. In support of this prediction, females acclimated at 32 degrees C were able to modify their resistance behaviour between exposure to 16 degrees C and 32 degrees C. The rate of copulations experienced by 32 inverted perpendicular C-acclimated females was 2.5 times greater at 32 degrees C than at 16 degrees C. In addition, acclimation at 32 degrees C significantly enhanced burst swimming performance at 32 degrees C but no effect of acclimation was observed at 16 degrees C. Our results suggest that female G. holbrooki are able to play a greater role in determining the outcome of sexual coercive mating interactions than previously thought. (c) 2006 The Association for the Shidy of Animal Behavioor. Published by Elsevier Ltd. All rights reserved.

Short- and long-term consequences of thermal variation in the larval environment of anurans

To survive adverse or unpredictable conditions in the ontogenetic environment, many organisms retain a level of phenotypic plasticity that allows them to meet the challenges of rapidly changing conditions. Larval anurans are widely known for their ability to modify behaviour, morphology and physiological processes during development, making them an ideal model system for studies of environmental effects on phenotypic traits. Although temperature is one of the most important factors influencing the growth, development and metamorphic condition of larval anurans, many studies have failed to include ecologically relevant thermal fluctuations among their treatments. We compared the growth and age at metamorphosis of striped marsh frogs Limnodynastes peronii raised in a diurnally fluctuating thermal regime and a stable regime of the same mean temperature. We then assessed the long-term effects of the larval environment on the morphology and performance of post-metamorphic frogs. Larval L. peronii from the fluctuating treatment were significantly longer throughout development and metamorphosed about 5 days earlier. Frogs from the fluctuating group metamorphosed at a smaller mass and in poorer condition compared with the stable group, and had proportionally shorter legs. Frogs from the fluctuating group showed greater jumping performance at metamorphosis and less degradation in performance during a 10-week dormancy. Treatment differences in performance could not be explained by whole-animal morphological variation, suggesting improved contractile properties of the muscles in the fluctuating group.

The effect of Cu addition and milling contaminations on the microstructure evolution of ball milled Al-Pb alloy during sintering

Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu (x = 0-7.0) bulk alloys were prepared by sintering the mechanically alloyed powders at various temperatures. The microstructure changes of the as consolidated powders in the course of sintering were analyzed by differential scanning calorimetry, scanning electron microscopy, X-ray diffraction analysis and transmission electron microscopy. It has been found that, with respect to the Al-10 wt.%Pb-x wt.%Cu alloy, CuAl2 and Cu9Al4 phases formed in the milling process, and the amount of CuAl2 phase increased while the Cu9Al4 phase disappeared gradually in the sintering process. In both Al-10 wt.%Pb and Al-10 wt.%Pb-x wt.%Cu alloys, the sintering process results in the coarsening of Pb phase and the growth rate of Pb phase fulfills the Lifshitz-Slyozov-Wagner equation even though the size of the Pb phase was in nanometer range. The Pb particle exhibits cuboctahedral morphology and has a cubic to cubic orientation relationship with the Al matrix. The addition of Cu strongly depressed the growth rate of Pb. Contamination induced by milling has apparent influence on the microstructure of the sintered alloys. Al7Cu2Fe and aluminium oxide phases were identified in the sintered alloys. The cuboctahedral morphology of Pb particles was broken up by the presence of the oxide phase. (c) 2006 Elsevier B.V. All rights reserved.

Factors affecting hot spot development in bulk coal and associated gas evolution

There are many physical factors that can affect the self-heating rate of coal. The presence of seam gas has often been referred to as inhibiting coal self-heating due to the limited access of oxidation sites created by the presence of the gas adsorbed on the coal pores. Similarly, the presence of bed moisture in the coal also acts as an inhibitor of oxidation by blocking access of air into the pores. Gas drainage of a coal seam prior to mining removes both gas and moisture from the seam. Bulk coal self-heating tests in a two-metre column on both gassy, as-mined and gas-drained, dried high volatile bituminous coal show that removal of gas and moisture from the coal accelerates the rate of self-heating to thermal runaway from 8.5 days to 4.25 days, from a start temperature of 30°C, with an airflow of 0.25 L/min. The corresponding gas evolution pattern for each of these situations is different. Therefore, it is necessary to take this change in coal condition into consideration when developing a spontaneous combustion management plan.

Bioactive glass sol-gel foam scaffolds:evolution of nanoporosity during processing and in situ monitoring of apatite layer formation using small- and wide-angle X-ray scattering

Recent work has highlighted the potential of sol-gel-derived calcium silicate glasses for the regeneration or replacement of damaged bone tissue. The work presented herein provides new insight into the processing of bioactive calcia-silica sol-gel foams, and the reaction mechanisms associated with them when immersed in vitro in a simulated body fluid (SBF). Small-angle X-ray scattering and wide-angle X-ray scattering (diffraction) have been used to study the stabilization of these foams via heat treatment, with analogous in situ time-resolved data being gathered for a foam immersed in SBF. During thermal processing, pore sizes have been identified in the range of 16.5-62.0 nm and are only present once foams have been heated to 400 degrees C and above. Calcium nitrate crystallites were present until foams were heated to 600 degrees C; the crystallite size varied from 75 to 145 nm and increased in size with heat treatment up to 300 degrees C, then decreased in size down to 95 rim at 400 degrees C. The in situ time-resolved data show that the average pore diameter decreases as a function of immersion time in SBF, as calcium phosphates grow on the glass surfaces. Over the same time, Bragg peaks indicative of tricalcium phosphate were evident after only 1-h immersion time, and later, hydroxycarbonate apatite was also seen. The hydroxycarbonate apatite appears to have preferred orientation in the (h,k,0) direction.