Object Kinetic Monte Carlo calculations of irradiated Fe-Cr dilute alloys: The effect of the interaction radius between substitutional Cr and self-interstitial Fe

ObjectKineticMonteCarlo models allow for the study of the evolution of the damage created by irradiation to time scales that are comparable to those achieved experimentally. Therefore, the essential ObjectKineticMonteCarlo parameters can be validated through comparison with experiments. However, this validation is not trivial since a large number of parameters is necessary, including migration energies of point defects and their clusters, binding energies of point defects in clusters, as well as the interactionradii. This is particularly cumbersome when describing an alloy, such as the Fe–Cr system, which is of interest for fusion energy applications. In this work we describe an ObjectKineticMonteCarlo model for Fe–Cr alloys in the dilute limit. The parameters used in the model come either from density functional theory calculations or from empirical interatomic potentials. This model is used to reproduce isochronal resistivity recovery experiments of electron irradiateddiluteFe–Cr alloys performed by Abe and Kuramoto. The comparison between the calculated results and the experiments reveal that an important parameter is the capture radius between substitutionalCr and self-interstitialFe atoms. A parametric study is presented on the effect of the capture radius on the simulated recovery curves.

MMonCa: A flexible and powerful new Kinetic Monte Carlo simulator

Kinetic Monte Carlo (KMC) is a widely used technique to simulate the evolution of radiation damage inside solids. Despite de fact that this technique was developed several decades ago, there is not an established and easy to access simulating tool for researchers interested in this field, unlike in the case of molecular dynamics or density functional theory calculations. In fact, scientists must develop their own tools or use unmaintained ones in order to perform these types of simulations. To fulfil this need, we have developed MMonCa, the Modular Monte Carlo simulator. MMonCa has been developed using professional C++ programming techniques and has been built on top of an interpreted language to allow having a powerful yet flexible, robust but customizable and easy to access modern simulator. Both non lattice and Lattice KMC modules have been developed. We will present in this conference, for the first time, the MMonCa simulator. Along with other (more detailed) contributions in this meeting, the versatility of MMonCa to study a number of problems in different materials (particularly, Fe and W) subject to a wide range of conditions will be shown. Regarding KMC simulations, we have studied neutron-generated cascade evolution in Fe (as a model material). Starting with a Frenkel pair distribution we have followed the defect evolution up to 450 K. Comparison with previous simulations and experiments shows excellent agreement. Furthermore, we have studied a more complex system (He-irradiated W:C) using a previous parametrization [1]. He-irradiation at 4 K followed by isochronal annealing steps up to 500 K has been simulated with MMonCa. The He energy was 400 eV or 3 keV. In the first case, no damage is associated to the He implantation, whereas in the second one, a significant Frenkel pair concentration (evolving into complex clusters) is associated to the He ions. We have been able to explain He desorption both in the absence and in the presence of Frenkel pairs and we have also applied MMonCa to high He doses and fluxes at elevated temperatures. He migration and trapping dominate the kinetics of He desorption. These processes will be discussed and compared to experimental results. [1] C.S. Becquart et al. J. Nucl. Mater. 403 (2010) 75

Bare-tether sheath and current: comparison of asymptotic theory and kinetic simulations in stationary plasma

Analytical expressions for current to a cylindrical Langmuir probe at rest in unmagnetized plasma are compared with results from both steady-state Vlasov and particle-in-cell simulations. Probe bias potentials that are much greater than plasma temperature (assumed equal for ions and electrons), as of interest for bare conductive tethers, are considered. At a very high bias, both the electric potential and the attracted-species density exhibit complex radial profiles; in particular, the density exhibits a minimum well within the plasma sheath and a maximum closer to the probe. Excellent agreement is found between analytical and numerical results for values of the probe radiusR close to the maximum radius Rmax for orbital-motion-limited (OML) collection at a particular bias in the following number of profile features: the values and positions of density minimum and maximum, position of sheath boundary, and value of a radius characterizing the no-space-charge behavior of a potential near the high-bias probe. Good agreement between the theory and simulations is also found for parametric laws jointly covering the following three characteristic R ranges: sheath radius versus probe radius and bias for Rmax; density minimum versus probe bias for Rmax; and (weakly bias-dependent) current drop below the OML value versus the probe radius for R > Rmax.

Optical Multistability and Oscillations in Hybrid Optical Bistable Systems with Short Delay Times

Optical instabilities in the output light from a bistable optical device (BOD) with a delayed feedback was predicted by Ikeda [1]. Gibbs et al. [2] gave the first experimental verification of this type of instabilities. From that time several groups have studied the instabilities of the BOD for different relations between the delay time tR and the time constant ح of the system. In a previous paper [3] an empirical and analytical study of instabilities in hybrid BOD was reported by us. The employed set up is shown in Fig. 1.

MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion

In this work, we introduce the Object Kinetic Monte Carlo (OKMC) simulator MMonCa and simulate the defect evolution in three different materials. We start by explaining the theory of OKMC and showing some details of how such theory is implemented by creating generic structures and algorithms in the objects that we want to simulate. Then we successfully reproduce simulated results for defect evolution in iron, silicon and tungsten using our simulator and compare with available experimental data and similar simulations. The comparisons validate MMonCa showing that it is powerful and flexible enough to be customized and used to study the damage evolution of defects in a wide range of solid materials.

First order mem-circuits: modeling, nonlinear oscillations and bifurcations

This paper presents a theoretical framework intended to accommodate circuit devices described by characteristics involving more than two fundamental variables. This framework is motivated by the recent appearance of a variety of so-called mem-devices in circuit theory, and makes it possible to model the coexistence of memory effects of different nature in a single device. With a compact formalism, this setting accounts for classical devices and also for circuit elements which do not admit a two-variable description. Fully nonlinear characteristics are allowed for all devices, driving the analysis beyond the framework of Chua and Di Ventra We classify these fully nonlinear circuit elements in terms of the variables involved in their constitutive relations and the notions of the differential- and the state-order of a device. We extend the notion of a topologically degenerate configuration to this broader context, and characterize the differential-algebraic index of nodal models of such circuits. Additionally, we explore certain dynamical features of mem-circuits involving manifolds of non-isolated equilibria. Related bifurcation phenomena are explored for a family of nonlinear oscillators based on mem-devices.

Temperature-depnedent disruption of circadian oscillations in woody plants-A mathematical approach

Poster

Vortex pinning vs superconducting wire network: origin of periodic oscillations induced by applied magnetic fields in superconducting films with arrays of nanomagnets

Hybrid magnetic arrays embedded in superconducting films are ideal systems to study the competition between different physical (such as the coherence length) and structural length scales such as are available in artificially produced structures. This interplay leads to oscillation in many magnetically dependent superconducting properties such as the critical currents, resistivity and magnetization. These effects are generally analyzed using two distinct models based on vortex pinning or wire network. In this work, we show that for magnetic dot arrays, as opposed to antidot (i.e. holes) arrays, vortex pinning is the main mechanism for field induced oscillations in resistance R(H), critical current Ic(H), magnetization M(H) and ac-susceptibility χ ac(H) in a broad temperature range. Due to the coherence length divergence at Tc, a crossover to wire network behaviour is experimentally found. While pinning occurs in a wide temperature range up to Tc, wire network behaviour is only present in a very narrow temperature window close to Tc. In this temperature interval, contributions from both mechanisms are operational but can be experimentally distinguished.

Numerical study of ultrasound induced non-linear shape and size bubble oscillations in viscoelastic media

The theoretical study of forced bubble oscillations is motivated by the importance of cavitation bubbles and oscillating encapsulated microbubbles (i.e. contrast agents) in medical sciences. In more details,theoretical studies on bubble dynamics addressing the sound-bubble interaction phenomenon provide the basis for understanding the dynamics of contrast agent microbubbles used in medical diagnosis and of non-linearly oscillating cavitation bubbles in the case of high-intensity ultrasound therapy. Moreover, the inclusion of viscoelasticity is of vital importance for an accurate theoretical analysis since most biological tissues and fluids exhibit non-Newtonian behavior.

Electron transport and azimuthal oscillations in hall thrusters

The Hall Effect Thruster (HET) is a type of satellite electric propulsion device initially developed in the 1960’s independently by USA and the former USSR. The development continued in the shadow during the 1970’s in the Soviet Union to reach a mature status from the technological point of view in the 1980’s. In the 1990’s the advanced state of this Russian technology became known in western countries, which rapidly restarted the analysis and development of modern Hall thrusters. Currently, there are several companies in USA, Russia and Europe manufacturing Hall thrusters for operational use. The main applications of these thrusters are low-thrust propulsion of interplanetary probes, orbital raising of satellites and stationkeeping of geostationary satellites. However, despite the well proven in-flight experience, the physics of the Hall Thruster are not completely understood yet. Over the last two decades large efforts have been dedicated to the understanding of the physics of Hall Effect thrusters. However, the so-called anomalous diffusion, short name for an excessive electron conductivity along the thruster, is not yet fully understood as it cannot be explained with classical collisional theories. One commonly accepted explanation is the existence of azimuthal oscillations with correlated plasma density and electric field fluctuations. In fact, there is experimental evidence of the presence of an azimuthal oscillation in the low frequency range (a few kHz). This oscillation, usually called spoke, was first detected empirically by Janes and Lowder in the 1960s. More recently several experiments have shown the existence of this type of oscillation in various modern Hall thrusters. Given the frequency range, it is likely that the ionization is the cause of the spoke oscillation, like for the breathing mode oscillation. In the high frequency range (a few MHz), electron-drift azimuthal oscillations have been detected in recent experiments, in line with the oscillations measured by Esipchuk and Tilinin in the 1970’s. Even though these low and high frequency azimuthal oscillations have been known for quite some time already, the physics behind them are not yet clear and their possible relation with the anomalous diffusion process remains an unknown. This work aims at analysing from a theoretical point of view and via computer simulations the possible relation between the azimuthal oscillations and the anomalous electron transport in HET. In order to achieve this main objective, two approaches are considered: local linear stability analyses and global linear stability analyses. The use of local linear stability analyses shall allow identifying the dominant terms in the promotion of the oscillations. However, these analyses do not take into account properly the axial variation of the plasma properties along the thruster. On the other hand, global linear stability analyses do account for these axial variations and shall allow determining how the azimuthal oscillations are promoted and their possible relation with the electron transport.

Caracterización de las Propiedades Dinámicas de la Seda de Araña

El término Biomimética se ha hecho común en los medios científicos, se refiere al trabajo de diversos científicos (ingenieros, químicos, físicos, biólogos, etc.) que tratan de copiar los procesos biológicos y aplicarlos en distintas áreas tecnológicas y científicas. En este campo científico, uno de los productos naturales que llama más la atención es la telaraña. Numerosos científicos en todo el mundo tratan de copiar las propiedades de la seda que produce la araña, y lo más interesante es que hasta intentan reproducir el método que usan las arañas para fabricar la seda De la bibliografía consultada, se desprende la importancia de la seda en la vida de las arañas, pues toda actividad que realizan tiene que ver de alguna manera con este elemento. Uno de estos elementos es la tela de araña orbicular, que representa el objeto principal para la supervivencia de la araña y su especie. Las investigaciones realizadas en esta línea nos proporcionan información sobre sus magníficas propiedades mecánicas como de resistencia, elasticidad y tenacidad del hilo de seguridad (seda MA) segregada por una araña de la especie Argiope Argentata. El enfoque de la presente tesis se realiza desde una perspectiva analítica-experimental, tomando a la tela de araña como una clase especial de sistemas pretensados, llamados Tensegrity Structures. Se desarrolla un modelo conceptual que describe en forma aproximada el comportamiento dinámico de una estructura hecha de seda MA. Haciendo uso de las técnicas experimentales de vibraciones libres se realizan los ensayos experimentales. La evaluación de los resultados analíticos y experimentales reflejan claramente que la función principal de la tela de araña es la de convertir energía cinética en energía de deformación y primordialmente en energía de disipación, el cual se efectúa gracias a las propiedades viscoelásticas de la seda. La araña en forma instintiva recurre a la ayuda del aire (como elemento disipador) para el buen funcionamiento de la tela de araña al momento de la captura de las presas, disipándose el 99% de la energía total en los tres primeros ciclos de oscilación de la tela de araña luego del impacto de la presa. ABSTRACT The term Biomimetic has become a very common word in the scientific world to describe the reproduction of the biological processes and its application in the different technogycal and scientific areas. One of the most notably natural product of this field is the Spiderweb. In the present days, many scientists of the world are active working in the reproduction of the proprieties of the Spiderweb. Most interesting even more, is the attempt to reproduce the process of the production of the Spiderweb by the spider. Most of the bibliography references deals whit the importance of the Spiderweb silk in the life of the spiders and the orbicular Spiderweb represents the spider survival and of the species. The research conducted in this field provide information about the excellent mechanical proprieties such us strength, elasticity and tenacity of the safety fiber (silk MA Drag-line) segregated by a spider of the Argiope Argentata species. The present work is oriented to an analytical and experimental study considering the Spiderweb as a special class of pre-stressed systems called Tensegrity (tensional integrity) structures. A conceptual model maked up by a cord und a point mass has been developed. This model approximates the dynamics performance of the structure made of the silk MA. The evaluation of the analytical and experimental results clear described that the main function of the Spiderweb is the transformation of the kinetic energy in deformation energy, and mainly in dissipation energy thank to the viscoelastic proprieties of the Spiderweb. With the help of the Spiderweb, the spider instinctively resorts to the help of the surroundings air as a dissipation element. This permits to dissipative the 99% of the total energy during the three first oscillations cycles of the web after the impact of the victim.

Global instability analysis and control of three-dimensional long open cavity flow

The three-dimensional wall-bounded open cavity may be considered as a simplified geometry found in industrial applications such as leading gear or slotted flats on the airplane. Understanding the three-dimensional complex flow structure that surrounds this particular geometry is therefore of major industrial interest. At the light of the remarkable former investigations in this kind of flows, enough evidences suggest that the lateral walls have a great influence on the flow features and hence on their instability modes. Nevertheless, even though there is a large body of literature on cavity flows, most of them are based on the assumption that the flow is two-dimensional and spanwise-periodic. The flow over realistic open cavity should be considered. This thesis presents an investigation of three-dimensional wall-bounded open cavity with geometric ratio 6:2:1. To this aim, three-dimensional Direct Numerical Simulation (DNS) and global linear instability have been performed. Linear instability analysis reveals that the onset of the first instability in this open cavity is around Recr 1080. The three-dimensional shear layer mode with a complex structure is shown to be the most unstable mode. I t is noteworthy that the flow pattern of this high-frequency shear layer mode is similar to the observed unstable oscillations in supercritical unstable case. DNS of the cavity flow carried out at different Reynolds number from steady state until a nonlinear saturated state is obtained. The comparison of time histories of kinetic energy presents a clearly dominant energetic mode which shifts between low-frequency and highfrequency oscillation. A complete flow patterns from subcritical cases to supercritical case has been put in evidence. The flow structure at the supercritical case Re=1100 resembles typical wake-shedding instability oscillations with a lateral motion existed in the subcritical cases. Also, This flow pattern is similar to the observations in experiments. In order to validate the linear instability analysis results, the topology of the composite flow fields reconstructed by linear superposition of a three-dimensional base flow and its leading three-dimensional global eigenmodes has been studied. The instantaneous wall streamlines of those composited flows display distinguish influence region of each eigenmode. Attention has been focused on the leading high-frequency shear layer mode; the composite flow fields have been fully recognized with respect to the downstream wave shedding. The three-dimensional shear layer mode is shown to give rise to a typical wake-shedding instability with a lateral motions occurring downstream which is in good agreement with the experiment results. Moreover, the spanwise-periodic, open cavity with the same length to depth ratio has been also studied. The most unstable linear mode is different from the real three-dimensional cavity flow, because of the existence of the side walls. Structure sensitivity of the unstable global mode is analyzed in the flow control context. The adjoint-based sensitivity analysis has been employed to localized the receptivity region, where the flow is more sensible to momentum forcing and mass injection. Because of the non-normality of the linearized Navier-Stokes equations, the direct and adjoint field has a large spatial separation. The strongest sensitivity region is locate in the upstream lip of the three-dimensional cavity. This numerical finding is in agreement with experimental observations. Finally, a prototype of passive flow control strategy is applied.

The forced oscillations of submerged bodies

The increasing use of very light structures in aerospace applications are given rise to the need of taking into account the effects of the surrounding media in the motion of a structure (as for instance, in modal testing of solar panels or antennae) as it is usually performed in the motion of bodies submerged in water in marine applications. New methods are in development aiming at to determine rigid-body properties (the center of mass position and inertia properties) from the results of oscillations tests (at low frequencies during modal testing, by exciting the rigid-body modes only) by using the equations of the rigid-body dynamics. As it is shown in this paper, the effect of the surrounding media significantly modifies the oscillation dynamics in the case of light structures and therefore this effect should be taken into account in the development of the above-mentioned methods. The aim of the paper is to show that, if a central point exists for the aerodynamic forces acting on the body, the motion equations for the small amplitude rotational and translational oscillations can be expressed in a form which is a generalization of the motion equations for a body in vacuum, thus allowing to obtain a physical idea of the motion and aerodynamic effects and also significantly simplifying the calculation of the solutions and the interpretation of the results. In the formulation developed here the translational oscillations and the rotational motion around the center of mass are decoupled, as is the case for the rigid-body motion in vacuum, whereas in the classical added mass formulation the six motion equations are coupled. Also in this paper the nonsteady motion of small amplitude of a rigid body submerged in an ideal, incompressible fluid is considered in order to define the conditions for the existence of the central point in the case of a three-dimensional body. The results here presented are also of interest in marine applications.

Accurate Analysis of Subharmonic Oscillations of V^2 and V^2Ic Controls Applied to Buck Converter

V2Ic control provides very fast dynamic performance to the Buck converter both under load steps and under voltage reference steps. However, the design of this control is complex since it is prone to subharmonic oscillations and several parameters affect the stability of the system. This paper derives and validates a very accurate modeling and stability analysis of a closed-loop V2Ic control using the Floquet theory. This allows the derivation of sensitivity analysis to design a robust converter. The proposed methodology is validated on a 5-MHz Buck converter. The work is also extended to V2 control using the same methodology, showing high accuracy and robustness. The paper also demonstrates, on the V2 control, that even a low bandwidth-linear controller can affect the stability of a ripple-based control.