OJAppS  Vol.6 No.7 , July 2016
Magnetohydrodynamic Calculations of Toroidal Fusion Reactor to Ensure Stable Control
The development of magnetic configurations to confine the stability fluid plasmas for fusion energy is a challenge that is a mixture of basic fusion engineering and invention. In order to keep the fusion reactions in the plasma to be continuing in the fusion reactors, the speed of tritium breeding (TBR) should be kept above a certain value. At the Apex fusion reactor, a fast flowing thin liquid wall has replaced the solid first wall concept of the traditional reactors. Behind the fast flowing thin liquid wall, a slower and thicker second liquid wall (coat) is present. Monte Carlo Random method (MCRS) is the general name for the solution of experimental and statistical problems with a random approach. This method is dependent upon the theory of probability. In the present work, Mhd impacts are investigated quite unimportant for Flibe salt solutions. In this study, the fissile fuel production calculations are done for a neutron wall load of 10 MW/m2 fissile fuel production rates of 238U(n, γ)239Pu and 232Th(n,γ)233U increases almost linearly with increased heavy metal content.

Received 31 May 2016; accepted 11 July 2016; published 14 July 2016

1. Introduction

Apex fusion reactor has an advantage of high neutron wall load capability due to its high beta. Controlled nuclear fusion can be a very important energy for future. At the Apex fusion reactor, a fast flowing thin liquid wall has replaced the solid first wall concept of the traditional reactors. In this study, Flibe salt solutions are used in the liquid wall of the reactor. The results obtained for both salt solutions are compared. At the construction of the liquid wall, it is possible to use liquid metals but it might be necessary to use electric isolators to manage magnetic-hydrodynamic drafts at the liquid metal flows, and as this will affect the plasma stability which maintains vital importance for fusion reactors, it is an important disadvantage for the usage of liquid metals. The application of the method to a problem is based on the idea of an approximate calculation of the parameter that is aimed to be calculated, using the random numbers of the problem by considering the results of this simulation [1] - [3] . A random number generator is a critical component in modern cryptographic communication systems, statistical simulation systems and scientific area incorporating Monte Carlo methods [4] . The general look of the linear equation that will form the random number series is as follows: xn+1 = (axn + b) mode m. Here mode means modular arithmetic procedure. The first element of the series xo is called the seed [5] - [7] . Constants a and b are chosen so that the elements of the series will be distributed anarchic and orderly. The same research company has used the following famous Random algorithm. (a = 69,069, b = 0): xn+1 = (69,069xn) mod 231 − 1. xn+1 = (16,807xn) mod 231. However, with appropriate algorithms, a computer may be forced to generate Random numbers chosen from an orderly scatter and takes real numbers between (0 - 1). Especially in nuclear field, using the nuclear transport codes, by the aid of MCNP (Monte Carlo n-particle transport code), physical and dynamical events are resulted by repeated utilization of the random numbers.

2. Method

2.1. Control of Plasma Position

Plasma is a collection of charged particles which is macroscopically neutral over a volume that is small compared to its dimensions. Nuclear fusion research has been most actively pursued by the use of hot plasma. In principle, the motion of each particle can be determined from Newton’s second Law and the electrostatic force that each particle exerts on all other particles. The magnetic plasma method must be heated and confined in order to utilize the nuclear fusion energy released from the plasma with a high temperature. The plasma consists of charged particles. Charged particles cannot cross magnetic field lines, but rotate around the lines. There is the outline of the idea of magnetic field lines in which to confine the plasma. In Figure 1 as shown, elliptic torus of fusion design cross-section, the particle drifts in the presence of a purely torodial magnetic field. The B magnetic gradient and curvature drifts produce a charge separation, which, in turn, produces an electric field with an associated outward drift for both ions and electrons and loss of confinement. In the geometry, a torus is a surface of motion generated by revolving a circle in 3D space.

When the magnetic field changes slowly, the magnetic moment is conserve. Therefore, if B is increased, mv2/2 = μB is also increased and the particles are heated. This kind of heating is called adiabatic heating. This solution shows that the guiding center motion consists of usual ExB drift and the new drift along E. Plasmas can be described as magnetohyrodynamic two fluids of ions and electrons with mass densities ρmi, ρme, charge density ρ, current density j, flow velocities Vi, Ve and pressures Pi, Pe.

The number of ions dNi in a small phase space volume,

Figure 1. Elliptic torus model of apex fusion design.

Magnetohydrodynamic (MHD) equations are the following;

2.2. The Utilization of the MCSR for Reactor Design

In the Monte Carlo Method, random numbers are generated for each element and for each random number, an examination is conducted in order to determine if the system obtained or not. There are practical manual methods for the generation of the Random numbers as well as a number of mathematical methods. According to the status of the analyzed system, there are practical methods like tossing dices, heads or tails or bingo, as will be understood from the name Monte Carlo. MCNP code includes Monte Carlo simulation and physics and mathematics subjects that have nuclear characteristics including certain subjects. During particle transport, neutrons, photons and electrons are used. MCNP in fact is a code that has been generalized by the Monte Carlo group in los-alamos labs for theoretic physics [5] [6] [8] [9] . MCNP possesses a general block that contains approximately 40,000 lines of Fortran and 1000 lines C source coder that includes interpretations and a general block that applies the program Monte Carlo Simulation is used in the simulation of some physical and mathematical systems. Generally, it is used for analyzing the effect of the random variables on the system and operates mostly with an algorithm that generates random numbers for the variables. The variables for which Monte Carlo simulation is to be conducted may be generated in a way that will suit a specific statistical distribution [4] [8] - [11] . The opposite of complementary probability distribution function (CDF) is used. The random numbers that are uniformly generated by the computer at a certain interval (0 − 1) may easily be transformed to any kind of distribution using the opposite of the

The design of 3d torus of fusion reactor the plasma and the blanket lithium around the plasma and the structural element that contains U and Th are shown in Figure 2.

Figure 2. 3D torus geometry of the apex reactor for flibe liquid flow [12] [13] .

The neutrons that are radiated as a result of the fusion reactions that occur at the plasma are captured at the blanket by uranium, thorium and the production of fissile fuel is provided. It is not possible to combust U238 in most reactors because the fusion impact incision is very low and the activation energy is too high. Th232 also is not an appropriate material for combustion [12] [14] [15] . Monte Carlo Random Simulation (MCRS), models the probable systems based on randoms. It is possible to calculate the integrals with the Monte Carlo Method through phase transitions [5] [6] [8] - [10] .

f(x) function at the [a, b] interval is considered : din order to show the size of the integral, quadratic total demonstration of such an integral shall be,

. According to the Average value theorem, the integral of a

f(x) function at [a, b] interval is equal to the area of the rectangle which is drawn at the average value of the function as seen on Figure 2; then average value is multiplied by (b − a) and the integral is calculated.

If the arithmetic average of N number of dots that are chosen randomly at the, [a, b] interval of the function is calculated, it will be: So, reaches the Monte Carlo integration formula:, the dots have not been arranged with equal intervals, they are chosen

completely randomly. As number N increases, the numeric calculation gets closer to the value of the I integral. But this approach is not as fast as expected. This slow approach is a general characteristic of the random methods. Here, according to the result that is provided without proof, the margin of error decreases slowly as [6] [11] [16] .

3. Results and Discussion

The apex fusion reactor used in the study has designed by using MCNP-4b (Endf/b-iv-vi) computer code and the Monte-Carlo technique, as 3-D torus [8] - [10] [16] . The dimensions for the apex reactor were taken from the aries-rs reactor design which has made in the framework of apex studies. A cross-sectional view of apex designed by using MCNP-4b computer code is shown in Figure 4. The inner region consists of plasma and vacuum. The first liquid wall is enclosed by the liquid blanket and then ferritic steel; these are further enclosed by a shield a stainless steel and ferritic steel zone. In this model, the radius of torus is 552 cm and minor radius starting from inner surface of first wall is 143 cm . By preserving a self sufficient tritium breeding ratio (TBR > 1.12) for a mole fraction up to 12% of UF4, the modified APEX reactor can produce up to ∼3000 kg of 233 U/year or ∼5000 kg of 239 Pu/year, assuming the same baseline fusion power production of 4000 MWth, as in the original APEX concept 50% UF4 in the coolant, the total energy output will increase to 6560 MWth or 9440 MWth, respectively. For a plant operation period of 30 full power years, the atomic displacement and helium production rates remain well below the presumable limits [17] - [19] .

The height of torus starting from center of first wall is 250 cm. The radius and thicknesses in one dimension are shown in detail in Figure 3 [10] [11] [16] [20] .

3.1. Neutron Tracking

Neutron tracking refers to the process of simulating the movement of a single neutron through the different material regions of the geometry. The shortest length of path that the neutron makes between the points of interaction is called a track. The track can also be cut short by a material boundary surface located between the two points. The entire set of tracks made from the initial emission to the final absorption or escape is called a neutron history.

The number of neutron histories is hence equal to the total number of emitted source neutrons. The basis of the simulation process is the sampling of the free path length between two collision points. Figure 4 shows the estimation of Monte Carlo random errors of number of neutron histories.

Figure 3. Inside-out cell and surface cross section of Apex fusion reactor model designed in MCNP-4b [10] [16] .

Figure 4. Estimation of Monte Carlo errors of number of neutron histories.

3.2. Tritium Breeding and Fissile Productions

A commercial fusion reactor must have a tritium breeding ratio of (TBR) > 1.05 to be self sustaining. Figure 5 shows the variation of total fissile material production speed for Flibe due to the ratio of UF4 and ThF4. As expected, the TBR decreases with increased heavy metal content. Tritium self sufficiency has been maintained in the range of molten salt mixtures (Flibe with up to 10% - 50% UF4) are shown results in Table 1. These molten salt mixtures would make the highest level of fissile breeding possible. Uranium has a higher neutron multiplication than thorium. On the other hand, the excess tritium gain for external fusion reactors is reduced significantly with the introduction of heavy metal into the coolant. When the fast and slow neutrons interaction with the lithium blanket, Li and T are produced as shown [15] [21] [22] ;

Table 1. Modeling of the apex for nuclear problems by simulating.

Figure 5. Total fissile material production speed for Flibe due to the ratio of UF4 and ThF4.

where, M is the energy multiplication factor, is fissile fuel production ratio, is ratio of total fission and υ is total number of fusion neutron. In Figure 5 as shown fissile fuel breeding are done fissile fuel breeding rates of 238U(n, ɣ)239Pu and 232 Th (n, ɣ)233U increase almost linearly with increased heavy metal content.

4. Conclusion

In fusion and hybrid reactors, it is necessary to provide the Tbr > 1.05 terms in order to maintain the continuity of the operation of the reactors. At the situation when there are no additions made to the blanket, the value Tbr = 1.20 for flibe is observed. The most important feature of hybrid reactor is that the first wall surrounding the plasma is liquid. For the D-T fusion driver, the energy multiplication factor, M, is defined in Figure 5, as the ratio of the total energy releases in the blanket to the incident fusion neutron energy [12] - [15] . For uranium, it also has been calculated that UF4 added hybrid reactor with UF4 ratios between 10% - 50% working with a wall load of 10 MW/m2 also works about twice as powerful as the original APEX Fusion reactor. Besides, the fissile fuel production provided at the reactor is also of great importance.

Cite this paper
Hançerlioğullari, A. , Kurnaz, A. and Madee, Y. (2016) Magnetohydrodynamic Calculations of Toroidal Fusion Reactor to Ensure Stable Control. Open Journal of Applied Sciences, 6, 402-408. doi: 10.4236/ojapps.2016.67041.
[1]   Hançerliogullari, A. (2006) Monte Carlo Simulation Method and MCNP Code System. Kastamonu Education Journal, 14, 545-556.

[2]   Sobol Ilya, M. (1994) A Primer for the Monte Carlo Method. CRC Press, London.

[3]   Rod Blais, J.A. and Zhang, Z. (2011) Exploring Pseudo and Chaotic Random Monte Carlo Simulations. Computer & Geosciences, 37, 928-934. http://dx.doi.org/10.1016/j.cageo.2011.01.009

[4]   Behnia, S., Akhavan, A., Akhshani, A. and Samsudin, A. (2011) A Novel Dynamic Model of Pseudo Random Number Generator. Journal of Computational and Applied Mathematics, 235, 3455-3463.

[5]   Binder, K. (1997) Applications of Monte Carlo Methods to Statistical Physics. Reports on Progress in Physics, 60, 487-559

[6]   Carter Jr., E.F. (1999) Random Walks, Markov Chains and the Monte Carlo Method, Monte Carlo Integration. Taygeta Scientific Inc. https://www.taygeta.com/rwalks/node5.html

[7]   Pevey, R.E. (1998) Monte Carlo Course Notes. Department of Nuclear Engineering, the University of Tennessee, Knoxville. Spring Semester. http://web.utk.edu/~rpevey/

[8]   Foster, D. and Artur, G. (1982) Average Neutronic Properties of “Prompt” Fission Poducts. Los Alamos National Laboraty Report la-9168-ms, February.

[9]   Johston, R.A. (1963) General Monte Carlo Neutronics Code. Lams-2856, Los Alamos.

[10]   Briesmeister, J. (1993) RSIC Computer Code Collection MCNP 4a, Monte Carlo N-Particle Transport Code System. Los Alamos National Laboratory, New Mexico.

[11]   Hançerliogullari, A. (2011) Determining of Energy Multiplication in the Apex Hybrid Reactor by Using ThF4 and UF4 Heavy Metals Salts. International Journal of Energy Research, 36, 1375-1382.

[12]   Abdou, M.A. (1999) Exploring Novel High Power Density Concepts for Attractive Fusion System. Fusion Engineering & Design, 4, 67-145.

[13]   Moir, R.W. (1997) Liquid First Walls for Magnetic Fusion Energy Configurations. Nuclear Fusion, 37, 557. http://dx.doi.org/10.1088/0029-5515/37/4/I13

[14]   Abdou, M.A. and The Apex Team (1999) Exploring Novel High Power Density Concepts for Attractive Fusion Systems. Fusion Engineering and Design, 45, 145-167.

[15]   Abdou, M.A., Ying, A. and Morley, N. and The Apex Team (2001) On the Exploration of Innovative Concepts for Fusion Chamber Technology. Fusion Engineering and Design, 54, 181-247.

[16]   Sarer, B., Hançerliogullari, A. and Savruk, N. (2005) A Design for Apex Hybrid Reactor Model by Using Monte Carlo Method. G.U. Journal of Science, 18, 17-25.

[17]   Sahin, H.M. (2007) Monte Carlo Calculation of Radiation Damage in First Wall of an Experimental Hybrid Reactor. Annals of Nuclear Energy, 34, 861-870.

[18]   Doucet, A., Johansen, M.A. and Tadic, B.V. (2010) On Solving Integral Equations Using Markov Chain Monte Carlo Methods. Applied Mathematics and Computation, 216, 2869-2880.

[19]   Beyli, M.U. and Acir, A. (2007) Utilization of Thorium in a High Power Density Hybrid Reactor with ?nnovative Coolants. Energy Conversion and Management, 48, 576-582.

[20]   Maziar, N., Matthew, C.W. and Needs, J.R. (2003) Quantum Monte Carlo Studies of Density Functional Theory. Mathematics and Computers in Simulation, 62, 463-470.

[21]   Mahdavi, M. and Asadi, E. (2013) Estimates of Tritium Produced Ratio in the Blancet of Fusion Reactors. Open Journal of Microphysics, 3, 8-11. http://dx.doi.org/10.4236/ojm.2013.31002

[22]   Hançerliogullari, A. (2014) Thermodynamics Properties of Molten Salt Technology Assessment for New Generation Fusion Reactors. Journal of Fusion Energy, 33, 5.