Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

Manuscript Title: GTNEUT: A code for the calculation of neutral particle transport in plasmas based on the Transmission and Escape Probability method.
Authors: J. Mandrekas.
Program title: GTNEUT
Catalogue identifier: ADTX
Journal reference: Comput. Phys. Commun. 161(2004)36
Programming language: Fortran 77.
Computer: PCs and Unix workstations.
Operating system: Solaris 8, 9, HP-UX 11i, Linux Red Hat v8.0, Windows NT/2000/XP.
RAM: 6.5 Mbytes
Keywords: Neutral transport in plasmas, Escape probability methods.
PACS: 52.25.Ya, 52.65-y, 52.55.Fa.
Classification: 19.11.

Nature of problem:
This code calculates the transport of neutral particles in thermonuclear plasmas in two-dimensional geometric configurations.

Solution method:
The code is based on the Transmission and Escape Probability (TEP) methodology [1], which is part of the family of integral transport methods for neutral particles and neutrons. The resulting linear system of equations is solved by standard direct linear system solvers (sparse and non-sparse versions are included).

The current version of the code can handle only one species of atomic neutrals.

Unusual features:
The program requires linking with the publicly available LAPACK linear algebra library which is usually included with the Fortran compilers of many UNIX vendors or can be obtained from NETLIB (www.netlib.org). To use the optional sparse matrix solver, the UMFPACK library is required which can be obtained from www.cise.ufl.edu/research/sparse/umfpack

Running time:
It depends on the size of the problem and the computing platform. For example, it takes 15.6 seconds of user time to run the second test problem on a SUN Ultra 10 workstation, using the sparse linear matrix solver.

[1] W. M. Stacey and J. Mandrekas, Nucl. Fusion 35(1994)1385.