Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Unsteady natural convection from a horizontal annulus filled with a porous medium

The unsteady natural convection flow from a horizontal cylindrical annulus filled with a non-Darcy porous medium has been studied. The unsteadiness in the problem arises due to the impulsive change in the wall temperature of the outer cylinder. The Navier–Stokes equations along with the energy equation governing the unsteady natural convection flow have been solved by the finite-volume method. The effect of time variation on the heat transfer is more pronounced only in a small time interval immediately after the start of the impulsive motion and the steady state is reached after certain time. The results show that the annulus completely filled with a porous medium has the best insulating effectiveness. Convection in the horizontal annulus is confined mostly at top and bottom regions. Hence, only these regions should be insulated. In case of annulus partially filled with a porous material, insulating the region near the outer cylinder is more effective than insulating the region near the inner cylinder. The effect of Darcy number on the heat transfer is more pronounced than that of the Grashof number.

Applications of dimensional reduction to electroweak and QCD matter

When heated to high temperatures, the behavior of matter changes dramatically. The standard model fields go through phase transitions, where the strongly interacting quarks and gluons are liberated from their confinement to hadrons, and the Higgs field condensate melts, restoring the electroweak symmetry. The theoretical framework for describing matter at these extreme conditions is thermal field theory, combining relativistic field theory and quantum statistical mechanics. For static observables the physics is simplified at very high temperatures, and an effective three-dimensional theory can be used instead of the full four-dimensional one via a method called dimensional reduction. In this thesis dimensional reduction is applied to two distinct problems, the pressure of electroweak theory and the screening masses of mesonic operators in quantum chromodynamics (QCD). The introductory part contains a brief review of finite-temperature field theory, dimensional reduction and the central results, while the details of the computations are contained in the original research papers. The electroweak pressure is shown to converge well to a value slightly below the ideal gas result, whereas the pressure of the full standard model is dominated by the QCD pressure with worse convergence properties. For the mesonic screening masses a small positive perturbative correction is found, and the interpretation of dimensional reduction on the fermionic sector is discussed.

Growth Mechanism of Cadmium Sulfide Nanocrystals

We investigate the growth kinetics of CdS nanocrystals in the quantum confinement regime using time-resolved small-angle X-ray scattering. In contrast to earlier reports for similar systems, we establish that the growth kinetics in this case follows the Lifshitz-Slyozov-Wagner theory, for not only growth of the average diameter of the nanocrystals but also the time dependence of the size distribution and the temperature dependence of the rate constant. This is the first rigorous example of the coarsening process in the quantum confinement (< 5 nm)regime. Ab initio studies for the reaction pathways provide a microscopic understanding of this finding.

Comparative study of some constitutive equations characterising non-Newtonian fluids

This paper compares, in a general way, the predictions of the constitutive equations given by Rivlin and Ericksen, Oldroyd, and Walters. Whether we consider the rotational problems in cylindrical co-ordinates or in spherical polar co-ordinates, the effect of the non-Newtonicity on the secondary flows is collected in a single parameterα which can be explicitly expressed in terms of the non-Newtonian parameters that occur in each of the above-mentioned constitutive equations. Thus, for a given value ofα, all the three fluids will have identical secondary flows. It is only through the study of appropriate normal stresses that a Rivlin-Ericksen fluid can be distinguished from the other two fluids which are indistinguishable as long as this non-Newtonian parameter has the same value.

Two-stream instability in plasma

The problem of two-stream instability in plasma is studied by specifying the importance of initial magnetic field associated with the motion of the charged particles and the boundary effects. In Part I the accurate initial steady state is studied when the streams of electrons and ions move with different uniform speeds in plasmas with plane and cylindrical geometry. In Part II, in order to show the effects of finiteness and inhomogeneity of the system, small transverse plasma oscillations are studied in the case of plane plasmas. The role of plasma-sheath oscillations at the boundaries is found to be very important in driving the instabilities associated with the electromagnetic modes. The numerical estimates of the growth rates of the instability are given for the specific case of the physical data in discharge tubes.

Free convection heat transfer in vertical annuli

Free convection heat transfer in vertical concentric, cylindrical annuli is investigated analytically and experimentally. The approximate double boundary layer model used by Emery and Chu for the case of vertical parallel plates is extended to the present case in obtaining heat transfer correlations in laminar free convection. Different correlations for the inner cylinder depending on the radius to the length ratio of the inner cylinder and the Rayleigh number, were used in the derivation of correlations for the annuli. The results for the case of short cylinders inside tubes are in agreement (within about 10 per cent) with the existing correlations. For other cases, namely long cylinders in annuli and wires in annuli, experiments conducted show the agreement of the analysis with experiments.

Alignment of the ALICE Inner Tracking System with cosmic-ray tracks

ALICE (A Large Ion Collider Experiment) is the LHC (Large Hadron Collider) experiment devoted to investigating the strongly interacting matter created in nucleus-nucleus collisions at the LHC energies. The ALICE ITS, Inner Tracking System, consists of six cylindrical layers of silicon detectors with three different technologies; in the outward direction: two layers of pixel detectors, two layers each of drift, and strip detectors. The number of parameters to be determined in the spatial alignment of the 2198 sensor modules of the ITS is about 13,000. The target alignment precision is well below 10 micron in some cases (pixels). The sources of alignment information include survey measurements, and the reconstructed tracks from cosmic rays and from proton-proton collisions. The main track-based alignment method uses the Millepede global approach. An iterative local method was developed and used as well. We present the results obtained for the ITS alignment using about 10^5 charged tracks from cosmic rays that have been collected during summer 2008, with the ALICE solenoidal magnet switched off.

A low temperature Weissenberg camera

Details of the design and operation of a Weissenberg camera suitable for x-ray investigations between -180°c and 200°c are presented. The camera employs a novel arrangement of spur and bevel gears to couple the goniometer spindle to the worm rod which controls the range of oscillation. The entire drive system and the goniometer assembly are mounted on a support which permits the insertion or removal of a cylindrical cassette from the gear-box side without disturbing the cooling assembly and the layer-line screen. The cassette can also be inserted from the opposite side. The specimen can be cooled either directly by a stream of liquid air or by the cold gas from its evaporation. Condensation of moisture at low temperatures is prevented by heating the layer-line tubes internally.

A coaxial capacitive voltage divider for oscilloscopic recording of nanosecond high voltage pulses

A coaxial capacitance voltage divider with a ratio of 110 and a rise time much less than 2.5 ns was developed for use with a transmission line pulse generator capable of producing 100 kV rectangular pulses of 2 mu s duration. The low voltage arm of the divider is a 3 cm long tube of titania (TiO2) turned out from a cylindrical compact. The compact was made by first pressing titania powder using a suitable binder and then sintering at controlled temperatures. The tube was slipped over the terminating end of the pulse-forming cable to form the divider with the cable capacitance.

Confined concrete—Its application in R.C. beams and frames

Stress-strain characteristics of concrete confined in steel binders have been determined. A new factor “confinement index” has been introduced for a quantitative measure of the confinement and using these results a “stress-block” has been developed. Tests have been made on simply supported reinforced concrete beams with spiral binder confinement and analysed on the basis of the proposed stress-block. Tests have also been made oon reinforced concrete portal frames and continuous beams with spiral binder confinement at sections of possible plastic hinge formation. An analysis of these tests indicates that a full redistribution of moments has taken place at ultimate.

Helical FDK Algorithms With No Backprojection Weight

In this paper, we present two new filtered backprojection (FBP) type algorithms for cylindrical detector helical cone-beam geometry with no position dependent backprojection weight. The algorithms are extension of the recent exact Hilbert filtering based 2D divergent beam reconstruction with no backprojection weight to the FDK type algorithm for reconstruction in 3D helical trajectory cone-beam tomography. The two algorithms named HFDK-W1 and HFDK-W2 result in better image quality, noise uniformity, lower noise and reduced cone-beam artifacts.

Analysis of solar radiation pressure induced coupled liberations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Structure and dynamics of benzene in one-dimensional channels

Molecular dynamics investigation of benzene in one-dimensional channel systems A1PO(4)-5, VPI-5, and carbon nanotube is reported. The results suggest that, in all the three host systems, the plane of benzene is almost perpendicular to the channel axis when the molecule is near the center of the channel and the plane of benzene is parallel to the channel axis when the molecule is near the wall of the channel. The density distribution of benzene as a function of channel length, z and the radial distance, r, from the channel axis is also different in the three host structures. Anisotropy in translational diffusion coefficient, calculated in body-fixed frame of benzene, suggests that benzene prefers to move with its plane parallel to the direction of motion in A1PO(4)-5 and VPI-5 whereas in carbon nanotube the motion occurs predominantly with the plane of the benzene perpendicular to the direction of motion.;Anisotropy associated with the rotational motion is seen to alter significantly in confinement as compared to liquid benzene. In A1PO(4)-5, the rotational anisotropy is reversed as compared to liquid benzene thereby suggesting that anisotropy arising out of molecular geometry can be reduced. Reorientational correlation times for C-6 and C-2 axes Of benzene are reported. Apart from the inertial decay of reorientational correlation function due to free, rotation, two other distinct regimes of decay are observed in narrower channels (AIPO(4)-5 and carbon nanotube): (i) an initial fast decay (0.5-2 ps) and (ii) a slower decay (>2 ps) of reorientational correlation function where C-6 decays slower than C-2 Similar to what is observed in liquid benzene. In the initial fast decay, it is seen that the decay for C-6 is faster than C-2 which is in contrast to what is observed in liquid benzene or for benzene confined in VPI-5.

Flow Generated By Pitched Blade Turbines Ii: Simulation Using Κ-Ε Model

Experimental data on average velocity and turbulence intensity generated by pitched blade downflow turbines (PTD) were presented in Part I of this paper. Part II presents the results of the simulation of flow generated by PTD The standard κ-ε model along with the boundary conditions developed in the Part 1 have been employed to predict the flow generated by PTD in cylindrical baffled vessel. This part describes the new software FIAT (Flow In Agitated Tanks) for the prediction of three dimensional flow in stirred tanks. The basis of this software has been described adequately. The influence of grid size, impeller boundary conditions and values of model parameters on the predicted flow have been analysed. The model predictions successfully reproduce the three dimensionality and the other essential characteristics of the flow. The model can be used to improve the overall understanding about the relative distribution of turbulence by PTD in the agitated tank