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.

Modified density matrix renormalization group algorithm for the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange: Comparison with field theory at large J(2)/J(1)

A modified density matrix renormalization group (DMRG) algorithm is applied to the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange J(1) and J(2) between first and second neighbors. The modified algorithm yields accurate results up to J(2)/J(1) approximate to 4 for the magnetic gap Delta to the lowest triplet state, the amplitude B of the bond order wave phase, the wavelength lambda of the spiral phase, and the spin correlation length xi. The J(2)/J(1) dependences of Delta, B, lambda, and xi provide multiple comparisons to field theories of the zigzag chain. The twist angle of the spiral phase and the spin structure factor yield additional comparisons between DMRG and field theory. Attention is given to the numerical accuracy required to obtain exponentially small gaps or exponentially long correlations near a quantum phase transition.

Diffusion of Be, Co and Mn impurities in compound semiconductors and glassy carbon studied by the modified radiotracer technique

The main method of modifying properties of semiconductors is to introduce small amount of impurities inside the material. This is used to control magnetic and optical properties of materials and to realize p- and n-type semiconductors out of intrinsic material in order to manufacture fundamental components such as diodes. As diffusion can be described as random mixing of material due to thermal movement of atoms, it is essential to know the diffusion behavior of the impurities in order to manufacture working components. In modified radiotracer technique diffusion is studied using radioactive isotopes of elements as tracers. The technique is called modified as atoms are deployed inside the material by ion beam implantation. With ion implantation, a distinct distribution of impurities can be deployed inside the sample surface with good con- trol over the amount of implanted atoms. As electromagnetic radiation and other nuclear decay products emitted by radioactive materials can be easily detected, only very low amount of impurities can be used. This makes it possible to study diffusion in pure materials without essentially modifying the initial properties by doping. In this thesis a modified radiotracer technique is used to study the diffusion of beryllium in GaN, ZnO, SiGe and glassy carbon. GaN, ZnO and SiGe are of great interest to the semiconductor industry and beryllium as a small and possibly rapid dopant hasn t been studied previously using the technique. Glassy carbon has been added to demonstrate the feasibility of the technique. In addition, the diffusion of magnetic impurities, Mn and Co, has been studied in GaAs and ZnO (respectively) with spintronic applications in mind.

Physicochemical studies on the gelation of hen egg yolk; separation of gelling protein components from yolk plasma

A study of the component(s) in egg yolk responsible for gelation of yolk on freezing and thawing has shown that granule-free yolk plasma, obtained by high-speed centrifugation of yolk, has the capacity to gel. As with the whole yolk, gelation of yolk plasma on freezing and thawing could be inhibited by additives such as sugars, sodium chloride, proteolytic enzymes, and phospholipase-A. Phospholipase-C, which induces gelation of whole yolk at room temperature, has a similar effect on yolk plasma. Yolk plasma has been separated into aggregating (gelling) and soluble fractions by delipidation, using formic acid. Each of these fractions consists of three or four protein components, as observed by gel filtration, ultracentrifugation, and agar electrophoresis. The proteins are glycoproteins and contain bound hexoses, hexosamine, and sialic acid. The gelation of yolk has been attributed to the interactions between protein molecules following disruption of lipid-protein bonds.

Note: Creep behaviour of high density polyethylene

The Ramberg-Osgood type constitutive law for creep suggested by Iyengar2 has been verified on high-density polyethylene. The time functions are evaluated from the experimental data of Scheweiker and Sidebottom3. It is found that the creep behaviour of the above material can be represented by the Ramberg-Osgood type law.

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.

Steady flow of a maxwell fluid in a porous cylindrical annulus with circular cross-section

When a fluid with memory is injected into any flow region some assumptions regarding the initial state of stress have to be made in order to determine the state of stress at any subsequent instant. For a Maxwell fluid, it is assumed that the fluid near the surface of injection is suddenly stressed and responds by starting flow in accordance with the mechanical model chosen. The flow of a Maxwell fluid with a single relaxation time has been determined under the above assumption in the following two cases: (i) annulus between two porous concentric circular cylinders, and (ii) space between two porous and infinitely extending parallel plates. The nature of flow in the present case is similar to that of the Reiner-Rivlin fluids obtained by Narasimhan2).

Physicochemical studies on the gelation of hen's egg yolk: Delipidation of yolk plasma by treatment with phospholipase-C and extraction with solvents

A method for the delipidation of egg yolk plasma using phospholipase-C, n-heptane, and 1-butanol has been described. An aggregating protein fraction and a soluble protein fraction were separated by the action of phospholipase-C. The aggregating protein fraction freed of most of the lipids by treatment with n-heptane and 1-butanol was shown to be the apolipoproteins of yolk plasma, whereas the soluble proteins were identified as the livetins. Carbohydrate and the N-terminal amino acid analysis of these protein fractions are reported. A comparison of these protein fractions with the corresponding fractions obtained by formic acid delipidation of yolk plasma has been made. The gelation of yolk plasma by the action of phospholipase-C has been interpreted as an aggregation of lipoproteins caused by ionic interactions. The role of lecithin in maintaining the structural integrity of lipoproteins has been discussed.

Stability of a cylindrical plasma in the presence of a monotonic field-II

The stability of an incompressible inviscid, perfectly conducting cylindrical plasma against azimuthal disturbances in the presence of a monotonic decreasing magnetic field having a constant pitch is discussed by using energy principle. The results obtained by this principle are compared for m = 1 mode (which is a dangerous mode in which there is a lateral shift of the entire column) with that obtained by normal mode analysis. It is found that m = 1 mode is always unstable. Further, an axial line current, external axial field and the surface tension tend to stabilise m ≠ modes.

Pd and Pt ions as highly active sites for the water-gas shift reaction over combustion synthesized zirconia and zirconia-modified ceria

Noble metal substituted ionic catalysts were synthesized by solution combustion technique. The compounds were characterized by X-ray diffraction, FT-Raman spectroscopy, and X-ray photoelectron spectroscopy. Zirconia supported compounds crystallized in tetragonal phase. The solid solutions of ceria with zirconia crystallized in fluorite structure. The noble metals were substituted in ionic form.The water-gas shift reaction was carried out over the catalysts.Negligible conversions were observed with unsubstituted compounds. The substitution of a noble metal ion was found to enhance the reaction rate. Equilibrium conversion was obtained below 250 degrees C in the presence of Pt ion substituted compounds. The formation of Bronsted acid-Bronsted base pairs was proposed to explain the activity of zirconia catalysts. The effect of oxide ion vacancies on the reactions over substituted ceria-zirconia solid solutions was established. (c)2010 Elsevier B.V. All rights reserved.

Harmonic and anharmonic oscillations of a cylindrical plasma in the presence of a uniform axial magnetic field and a steady axial current

In the present note we have studied the harmonic and anharmonic oscillations of cylindrical plasma using Lagrangian formalism. In order to study the harmonic oscillations, the equations are linearized and the resulting equation for the displacement has been numerically solved. For situations present in thermonuclear reactors, the presence of axial magnetic field is found necessary to make the periods of oscillation to become comparable with the time required for the thermonuclear reactions to set in. A detailed analysis of the anharmonic oscillations reveals that the significant interaction is between the first and the second mode. The fundamental period of anharmonic oscillation is more than the corresponding period of harmonic oscillations by 9·2%. Graphs have been drawn for the amplitudes of relative variations in density and magnetic field and of the time-varying part of anharmonic oscillation.

Top quark mass measurement in the lepton plus jets channel using a modified matrix element method

We report a measurement of the top quark mass, m_t, obtained from ppbar collisions at sqrt(s) = 1.96 TeV at the Fermilab Tevatron using the CDF II detector. We analyze a sample corresponding to an integrated luminosity of 1.9 fb^-1. We select events with an electron or muon, large missing transverse energy, and exactly four high-energy jets in the central region of the detector, at least one of which is tagged as coming from a b quark. We calculate a signal likelihood using a matrix element integration method, with effective propagators to take into account assumptions on event kinematics. Our event likelihood is a function of m_t and a parameter JES that determines /in situ/ the calibration of the jet energies. We use a neural network discriminant to distinguish signal from background events. We also apply a cut on the peak value of each event likelihood curve to reduce the contribution of background and badly reconstructed events. Using the 318 events that pass all selection criteria, we find m_t = 172.7 +/- 1.8 (stat. + JES) +/- 1.2 (syst.) GeV/c^2.

Steady laminar flow with suction or injection and heat transfer through porous pipe

Following the method due to Bhatnagar (P. L.) [Jour. Ind. Inst. Sic., 1968, 1, 50, 1], we have discussed in this paper the problem of suction and injection and that of heat transfer for a viscous, incompressible fluid through a porous pipe of uniform circular cross-section, the wall of the pipe being maintained at constant temperature. The method utilises some important properties of differential equations and some transformations that enable the solution of the two-point boundary value and eigenvalue problems without using trial and error method. In fact, each integration provides us with a solution for a suction parameter and a Reynolds number without imposing the conditions of smallness on them. Investigations on non-Newtonian fluids and on other bounding geometries will be published elsewhere.