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.

Intramembrane charge movements and excitation– contraction coupling expressed by two-domain fragments of the Ca2+ channel

To investigate the molecular basis of the voltage sensor that triggers excitation–contraction (EC) coupling, the four-domain pore subunit of the dihydropyridine receptor (DHPR) was cut in the cytoplasmic linker between domains II and III. cDNAs for the I-II domain (α1S 1–670) and the III-IV domain (α1S 701-1873) were expressed in dysgenic α1S-null myotubes. Coexpression of the two fragments resulted in complete recovery of DHPR intramembrane charge movement and voltage-evoked Ca2+ transients. When fragments were expressed separately, EC coupling was not recovered. However, charge movement was detected in the I-II domain expressed alone. Compared with I-II and III-IV together, the charge movement in the I-II domain accounted for about half of the total charge (Qmax = 3 ± 0.23 vs. 5.4 ± 0.76 fC/pF, respectively), and the half-activation potential for charge movement was significantly more negative (V1/2 = 0.2 ± 3.5 vs. 22 ± 3.4 mV, respectively). Thus, interactions between the four internal domains of the pore subunit in the assembled DHPR profoundly affect the voltage dependence of intramembrane charge movement. We also tested a two-domain I-II construct of the neuronal α1A Ca2+ channel. The neuronal I-II domain recovered charge movements like those of the skeletal I-II domain but could not assist the skeletal III-IV domain in the recovery of EC coupling. The results demonstrate that a functional voltage sensor capable of triggering EC coupling in skeletal myotubes can be recovered by the expression of complementary fragments of the DHPR pore subunit. Furthermore, the intrinsic voltage-sensing properties of the α1A I-II domain suggest that this hemi-Ca2+ channel could be relevant to neuronal function.

Indirect coupling of phosphate release to de novo tension generation during muscle contraction.

A key question in muscle contraction is how tension generation is coupled to the chemistry of the actomyosin ATPase. Biochemical and mechanochemical experiments link tension generation to a change in structure associated with phosphate release. Length-jump and temperature-jump experiments, on the other hand, implicate phase 2slow, a significantly faster, markedly strain-sensitive kinetic process in tension generation. We use a laser temperature jump to probe the kinetics and mechanism of tension generation in skinned rabbit psoas fibers--an appropriate method since both phosphate release and phase 2slow are readily perturbed by temperature. Kinetics characteristic of the structural change associated with phosphate release are observed only when phosphate is added to fibers. When present, it causes a reduction in fiber tension; otherwise, no force is generated when it is perturbed. We therefore exclude this step from tension generation. The kinetics of de novo tension generation by the temperature-jump equivalent of phase 2slow appear unaffected by phosphate binding. We therefore propose that phosphate release is indirectly coupled to de novo tension generation via a steady-state flux through an irreversible step. We conclude that tension generation occurs in the absence of chemical change as the result of an entropy-driven transition between strongly bound crossbridges in the actomyosin-ADP state. The mechanism resembles the operation of a clock, with phosphate release providing the energy to tension the spring, and the irreversible step functions as the escapement mechanism, which is followed in turn by tension generation as the movement of the hands.

Sprint patterns in rugby union players during competition

The purpose of this study was to characterize sprint patterns of rugby union players during competition. Velocity profiles (60 m) of 28 rugby players were initially established in testing from standing, walking, jogging, and striding starts. During competition, the individual sprinting patterns of 17 rugby players were determined from video by using the individual velocity profiles. Forwards commenced sprints from a standing start most frequently (41%), whereas backs sprinted from standing (29%), walking (29%),jogging (29%), and occasionally striding (13%) starts. Forwards and backs achieved speeds in excess of 90% maximal velocity (Vmax) on 5 +/- 4 and 9 +/- 4 occasions (similar to 50% of the sprints performed), respectively, during competition. The higher frequency of sprinting for the backs compared with the forwards highlights the importance of speed training for this positional group. The similar relative distribution of velocities achieved during competition for forwards and backs suggests both positional groups should train acceleration and Vmax qualities. The backs should have a higher total volume of sprint training. Sprinting efforts should be performed from a variety of starting speeds to mimic the movement patterns of competition.

Force sensing for measuring human body movement

The research developed in this thesis explores the sensing and inference of human movement in a dynamic way, as opposed to conventional measurement systems, that are only concerned with discrete evaluations of stimuli in sequential time. Typically, conventional approaches are used to infer the dynamic movement of the body; such as vision and motion tracking devices, with either a human diagnosis or complex image processing algorithm to classify the movement. This research is therefore the first of its kind to attempt and provide a movement classifying algorithm through the use of minimal sensing points, with the application for this novel system, to classify human movement during a golf swing. There are two main categories of force sensing. Firstly, array-type systems consisting of many sensing elements, and are the most commonly researched and commercially available. Secondly, reduced force sensing element systems (RFSES) also known as distributive systems have only been recently exploited in the academic world. The fundamental difference between these systems is that array systems handle the data captured from each sensor as unique outputs and suffer the effects of resolution. The effect of resolution, is the error in the load position measurement between sensing elements, as the output is quantized in terms of position. This can be compared to a reduced sensor element system that maximises that data received through the coupling of data from a distribution of sensing points to describe the output in discrete time. Also this can be extended to a coupling of transients in the time domain to describe an activity or dynamic movement. It is the RFSES that is to be examined and exploited in the commercial sector due to its advantages over array-based approaches such as reduced design, computational complexity and cost.

Short-term impact of soybean management on ammonia oxidizers in a Brazilian savanna under restoration as revealed by coupling different techniques

Acknowledgments The authors are very grateful to Mr. Fabiano Bielefeld Nardotto, owner of the Tabapuã dos Pireneus farm, for allowing our free movement around the farm and collection of soil samples, as well as providing information about soybean cultivation. The authors also thank Dr. Plínio de Camargo, who performed the isotopic analysis in the CENA laboratory at the University of São Paulo (USP). This work was supported by grants from the National Council of Technological and Scientific Development (CNPq), Brazilian Federal Agency for Support and Evaluation of Graduate Education (CAPES), and Foundation for Research Support of Distrito Federal (FAP-DF).

Neural Network Evidence for the Coupling of Presaccadic Visual Remapping to Predictive Eye Position Updating.

As we look around a scene, we perceive it as continuous and stable even though each saccadic eye movement changes the visual input to the retinas. How the brain achieves this perceptual stabilization is unknown, but a major hypothesis is that it relies on presaccadic remapping, a process in which neurons shift their visual sensitivity to a new location in the scene just before each saccade. This hypothesis is difficult to test in vivo because complete, selective inactivation of remapping is currently intractable. We tested it in silico with a hierarchical, sheet-based neural network model of the visual and oculomotor system. The model generated saccadic commands to move a video camera abruptly. Visual input from the camera and internal copies of the saccadic movement commands, or corollary discharge, converged at a map-level simulation of the frontal eye field (FEF), a primate brain area known to receive such inputs. FEF output was combined with eye position signals to yield a suitable coordinate frame for guiding arm movements of a robot. Our operational definition of perceptual stability was "useful stability," quantified as continuously accurate pointing to a visual object despite camera saccades. During training, the emergence of useful stability was correlated tightly with the emergence of presaccadic remapping in the FEF. Remapping depended on corollary discharge but its timing was synchronized to the updating of eye position. When coupled to predictive eye position signals, remapping served to stabilize the target representation for continuously accurate pointing. Graded inactivations of pathways in the model replicated, and helped to interpret, previous in vivo experiments. The results support the hypothesis that visual stability requires presaccadic remapping, provide explanations for the function and timing of remapping, and offer testable hypotheses for in vivo studies. We conclude that remapping allows for seamless coordinate frame transformations and quick actions despite visual afferent lags. With visual remapping in place for behavior, it may be exploited for perceptual continuity.

A comparison of the physiological demands of rugby union match-play when determined by absolute (ABS) or individual (IND) velocity bands in Global Positioning System (GPS) analysis

Since turning professional in 1995 there have been considerable advances in the research on the demands of rugby union, largely using Global Positioning System (GPS) analysis over the last 10 years. A systematic review on the use of GPS, particularly the setting of absolute (ABS) and individual (IND) velocity bands in field based, intermittent, high-intensity (HI) team sports was undertaken. From 3669 records identified, 38 studies were included for qualitative analysis. Little agreement on the definition of movement intensities within team sports was found, only three papers, all on rugby union, had used IND bands, with only one comparing ABS and IND methods. Thus, the aim of this study was to determine if there is a difference in the demands within positions when comparing ABS and IND methods for GPS analysis and if these differences are significantly different between the forward and back positional groups. A total of 214 data files were recorded from 26 players in 17 matches of the 2015/2016 Scottish BT Premiership. ABS velocity zones 1-7 were set at 1) 0-6, 2) 6.1-11, 3) 11.1-15, 4) 15.1-18, 5) 18.1-21, 6) 21.1-15 and 7) 25.1-40km.h-1 while IND zones 1-7 were 1) <20, 2) 20-40, 3) 40-50, 4) 50-70, 5) 70-80, 6) 80-95 and 7) 95-100% of player’s individually determined maximum velocity (Vmax). A 40m sprint test measured Vmax using OptaPro S4 10 Hz (catapult, Australia) GPS units to derive IND bands. The same GPS units were worn during matches. GPS outputs analysed were % distance, % time, high intensity efforts (HIEs) over 18.1 km.h-1 / 70% max velocity and repeated high intensity efforts (RHIEs) which consists of three HIEs in 21secs. General linear model (GLM) analysis identified a significant difference in the measurement of % total distance covered, between the ABS and IND methods in all zones for forwards (p<0.05) and backs (p<0.05). This difference was also significant between forwards and backs in zones 1, shown as mean difference ± standard deviation (3.7±0.7%), 6 (1.2±0.4%) and 7 (1.0±0.0%) respectively (p<0.05). Percentage time estimations were significantly different between ABS and IND analysis within forwards in zones 1 (1.7±1.7%), 2 (-2.9±1.3%), 3 (1.9±0.8%), 4 (-1.4±0.8%) and 5 (0.2±0.4%), and within backs in zones 1 (-10±1.5%), 2 (-1.2±1.1%), 3 (1.8±0.9%) and 5 (0.6±0.5%) (p<0.05). The difference between groups was significant in zones 1, 2, 4 and 5 (p<0.05). The number of HIEs was significantly different between forwards and backs in zones 6 (6±2) and 7 (3±2). RHIEs were significantly different between ABS and IND for forwards (1±2, p<0.05) although not between groups. Until more research on the differences in ABS and IND methods is carried out, then neither can be deemed a criterion method. In conclusion, there are significant differences between the ABS and IND methods of GPS analysis of the physical demands of rugby union, which must be considered when used to inform training load and recovery to improve performance and reduce injuries.

Electrocardiogram, Heart Movement And Heart Rate In The Awake Gecko (hemidactylus Mabouia).

The electrocardiogram (ECG) is the simplest and most effective non-invasive method to assess the electrical activity of the heart and to obtain information on the heart rate (HR) and rhythm. Because information on the HR of very small reptiles (body mass <10 g) is still scarce in the literature, in the present work we describe a procedure for recording the ECG in non-anesthetized geckos (Hemidactylus mabouia, Moreau de Jonnès, 1818) under different conditions, namely manual restraint (MR), spontaneous tonic immobility (TI), and in the non-restrained condition (NR). In the gecko ECG, the P, QRS and T waves were clearly distinguishable. The HR was 2.83 ± 0.02 Hz under MR, which was significantly greater (p < 0.001) than the HR under the TI (1.65 ± 0.09 Hz) and NR (1.60 ± 0.10 Hz) conditions. Spontaneously beating isolated gecko hearts contracted at 0.84 ± 0.03 Hz. The in vitro beating rate was affected in a concentration-dependent fashion by adrenoceptor stimulation with noradrenaline, as well as by the muscarinic cholinergic agonist carbachol, which produced significant positive and negative chronotropic effects, respectively (p < 0.001). To our knowledge, this is the first report on the ECG morphology and HR values in geckos, particularly under TI. The methodology and instrumentation developed here are useful for non-invasive in vivo physiological and pharmacological studies in small reptiles without the need of physical restraint or anesthesia.

Linearization of Chains and Sideward Movement

360

Avaliação do desempenho anaeróbio de atletas de rugby em cadeira de rodas = : Anaerobic performance evaluation of wheelchair rugby athletes

Universidade Estadual de Campinas . Faculdade de Educação Física

Transition to complete synchronization in phase-coupled oscillators with nearest neighbor coupling

We investigate synchronization in a Kuramoto-like model with nearest neighbor coupling. Upon analyzing the behavior of individual oscillators at the onset of complete synchronization, we show that the time interval between bursts in the time dependence of the frequencies of the oscillators exhibits universal scaling and blows up at the critical coupling strength. We also bring out a key mechanism that leads to phase locking. Finally, we deduce forms for the phases and frequencies at the onset of complete synchronization.

Border figure detection using a phase oscillator network with dynamical coupling

Oscillator networks have been developed in order to perform specific tasks related to image processing. Here we analytically investigate the existence of synchronism in a pair of phase oscillators that are short-range dynamically coupled. Then, we use these analytical results to design a network able of detecting border of black-and-white figures. Each unit composing this network is a pair of such phase oscillators and is assigned to a pixel in the image. The couplings among the units forming the network are also dynamical. Border detection emerges from the network activity.

Epidermal growth factor in liposomes may enhance osteoclast recruitment during tooth movement in rats

Objective: To evaluate the effects of local administration of epidermal growth factor (EGF) located within liposomes on recruitment of osteoclasts during mechanical force in rats. Materials and Methods: An orthodontic elastic band was inserted between the left upper first and second molars, to move mesially the first molar. Rats were randomly divided into 4 groups (n = 8): EGF (2 ng/mu L) located within liposomes (group 1), liposomes only (group 2), soluble EGF (2 ng/mu L; group 3), or vehicle alone (group 4). The solutions were injected into the region of the root furcation of the left first molar after elastic band insertion. Tooth movement was measured using a plaster model of the maxilla, and the number of osteoclasts recruited at the pressure side of the first molar was histologically evaluated. Results: Intergroup analysis showed that there was no significant difference between group 2 and group 4 (P >.05) and between group 1 and group 3 (P >.05). However, group 1 and group 3 exhibited greater differences in tooth movement than group 2 and group 4 (P <.05). On the other hand, group 1 showed greater tooth movement than groups 2 and 4 with statistical significance (P <.01). The increase in the number of osteoclasts in group 1 was significantly higher than in the other groups (P <.05). Conclusion: Exogenous EGF-liposome administration has an additive effect when compared with soluble EGF on the rate of osteoclast recruitment, producing faster bone resorption and tooth movement.