Click on the advert above to visit the company web site

Product category: Simulation, modelling and validation software
News Release from: Vector Fields | Subject: Opera
Edited by the Engineeringtalk Editorial Team on 16 September 2005

Software simulates fusion reactor on a
standard PC

Electromagnetic simulation is playing a major role in the design of a Tokamak fusion reactor by calculating eddy currents in the diverter that controls plasma particle exhaust and power removal.

Electromagnetic simulation is playing a major role in the conceptual design of a Tokamak fusion reactor by calculating eddy currents in the diverter, a device that controls plasma particle exhaust and power removal in the machine A preconceptual design study for the Fusion Ignition Research Experiment (FIRE) is underway to assess near term opportunities for advancing the scientific understanding of self-heated fusion plasmas

One of the key issues is to design plasma facing components (PFCs) that are capable of withstanding what is known as a disruption where plasma pressure and toroidal current are quenched on a fast time scale.

The eddy currents produced during these transients are sufficiently high that the resulting stresses are capable of damaging plasma facing components.

Engineers at Sandia National Laboratories, Albuquerque, New Mexico, used PC-based magnetic simulation software to predict the eddy currents and resulting Lorentz forces caused by disruptions.

These calculated forces were then input to stress analysis programs used to optimise the mechanical design of diverter components.

For many decades, scientists have tried to harness the power of nuclear fusion, the process that occurs when the nuclei of hydrogen and other light elements collide and unite, releasing enormous amounts of energy.

The sun is heated by the fusion of hydrogen in its core and the energy radiates as light at the surface.

In this fusion reaction, hydrogen isotopes combine to form hot helium nuclei.

The fuel needed for fusion is widely distributed eliminating the possibility of supply cartels.

Fusion reactors contain only enough fuel for about 1 second of Operation making a meltdown impossible.

Another advantage of fusion is that it produces only a small fraction of the highly radioactive waste that is produced by today's nuclear plants.

Furthermore, the nuclear isotopes of hydrogen produced by fusion stay toxic for only tens of years rather than thousands of years for the fuel of traditional fission nuclear reactors.

One of the most challenging problems with potential fusion reactors is how the plasma is confined.

The hot nuclei in the plasma have to be held together at high densities long enough to allow enough of them to fuse to produce more energy than was used in heating the plasma.

Material containers cannot confine such hot particles because of energy loss to the walls.

Instead, scientists use a magnetic field to confine the hot particles.

Magnetic fields exert force only on charged or magnetised objects.

Room temperature gas atoms have as many electrons as protons, and thus are neutral.

When the gas is heated to a very high temperature, the electrons have so much energy that they are no longer bound to the nuclei.

The resulting gas is called a plasma because the energetic electrons and nuclei move freely.

Magnetic fields can bend the path of hot, fast-moving nuclei and electrons into loops.

Instead of hitting the wall, the particles in a plasma can be confined to travel in circles within a chamber.

If the particles in the plasma never collided with each other, then they would be trapped in loops indefinitely by the magnetic field.

However, collisions between the plasma particles do occur, causing the plasma to leak out of the confining magnetic fields.

Over many years, fusion scientists have been perfecting the magnetic containers that hold the hot plasma while they apply the energy needed to sustain the fusion reaction.

The most successful design, the Tokamak, uses a toroidal or doughnut-shaped chamber surrounded by strong electromagnets.

The main components in a Tokamak are the magnet system, external heating and drive system and the diverter system.

The magnet system is used to generate the magnetic field used to confine the plasma.

The external heating system, which can consist of ohmic, RF and neutral beam heating, raises the plasma temperature to a level needed to initiate fusion reactions.

The diverter is used to control the plasma edge properties and remove ionised impurities and "helium ash" that has already given up its energy to the plasma.

The diverter, the first solid surface that makes contact with the plasma, must be designed to withstand not only normal Operations but also events known as disruptions.

These disruptions generate hundreds of thousands of amps of eddy currents within milliseconds that in turn create high Lorentz forces.

Such forces determine the size and strength of the support plates, the connections and the actively cooled, tungsten-brush plasma-facing armour.

Initially, a simple L/R analysis was performed in order to estimate the eddy currents induced in the FIRE diverter structures for the case of a plasma disruption.

An analytical approximation of the inductance of a rectangular plate was used to estimate induced eddy currents for the outer and inner diverter.

The force on the diverter modules would then be 1.9MN for the outer diverter and 2.8MN for the inner.

Dennis Youchison, a senior staff member at Sandia was given the assignment of a more detailed follow-on analysis.

"These back of the envelope calculations were useful for determining the rough magnitude of the forces that were involved", Youchison said.

"But because of the assumptions that were involved in producing them and particularly the fact that they did not take the geometry of the diverter into account, we realised we needed more accurate predictions".

"Fusion researchers have developed electromagnetic simulation codes but those codes evolved over many years with less than intuitive user interfaces making it very difficult to model complex geometries such as the diverter and quickly obtain results needed for the design analysis".

"Another problem with these codes is that they run on expensive, high-end mainframes with restricted access".

"For these reasons", Youchison said, "I decided to run the simulation on Opera, a personal-computer based electromagnetic simulation program from Vector Fields that I have used on several previous projects".

"This program has a modern Windows user interface that makes it easy to model complicated geometries and also generate intuitive graphical output".

"Just as important, Opera can analyse complicated problems on a Pentium PC and that provides ease of access and reduces support costs".

Youchison input the complicated model geometry using the Opera preprocessor in a matter of a few days.

Taking advantage of the symmetrical design of the Tokamak, he modelled a 1:16 section of the torus including the structure of the reactor and the electromagnetic coils.

He modelled the plasma as a compilation of 668 different current filaments.

The currents in each filament change over time to simulate how the plasma responds during a disruption event.

Output from a plasma simulation code, TSC developed at Princeton Plasma Physics Laboratory, was used to generate the time response functions used to drive the filaments.

The finite element mesh contains 92,000 elements and solves the time response equations for each of the filaments.

It took 46 hours to complete a 300ms simulation on a 1.7GHz PC.

About 90% of the total processing time was required to compute the superimposed fields from the 710 coils including the plasma filaments.

The analysis provided graphical output consisting of vector and contour plots that showed the magnetic fields, current densities and force densities superimposed on the surfaces of the model.

Several disruption scenarios were analysed.

The most severe case, a vertical disruption event, resulted in a maximum total force to an outer diverter plate of 86.7kN produced by 131kA of current intercepting a 12T magnetic field.

Youchison pointed out that this is much reduced from the initial scoping calculations, because the toroidally continuous diverter was segmented into 16 insulated plates.

It was obvious from the initial Opera runs that this reduced the size of the eddy current loops and hence the resultant Lorentz force on any one plate.

The baffle and inner diverter were likewise segmented.

Each inner diverter segment experiences a maximum total force of 14.3kN and each baffle segment endures a 22.4kN force.

However, the sum of forces on each segment over the entire outer diverter is 1.4MN.

Insulation that isolated the diverter from the passive plates and the inner vacuum vessel also forced the eddy currents to flow along less deleterious paths.

The toroidally continuous passive plates carry most of the eddy currents produced during the disruption, about 3.1MA, but as these currents are parallel to the toroidal magnetic field they only see the contribution of the much lower poloidal field to the Lorentz forces.

In contrast, due to their location, the diverter PFCs must contend with a large toroidal field component and higher forces.

Youchison used Opera to calculate the Lorentz force densities and exported these values in tabular format into Mathcad where he performed four-point Gauss-Legendre quadrature analysis to determine the total forces.

The total forces were then exported to a mechanical finite element analysis software package used by the Boeing engineers performing the mechanical design of the diverter.

"The ability to accurately estimate eddy currents is making it possible to design the diverter based on accurate load conditions", Youchison said.

"This is just one of many design tasks that must be completed in order to carry out the design of any burning plasma experiment such as FIRE that may eventually explore issues that are critical to the development of commercial fusion reactors".

• Vector Fields: contact details and other news
• Email this article to a colleague
• Register for the free Engineeringtalk email newsletter
• Engineeringtalk Home Page

Search the Pro-Talk network of sites

Put your news on Engineeringtalk  |   Advertise  |   Get the free Engineeringtalk newsletter  |   Engineeringtalk Home
About Engineeringtalk  |   Copyright © 2000-2008 Pro-Talk Ltd, UK