One-dimensional transport code for field-reversed configurations

A radial transport code for Field-Reversed Configurations (FRCs) is described. Assuming quasi-neutrality and neglecting inertia and viscous force, the evolution of particle densities, temperatures, and magnetic field through a series of equilibrium states is simulated. To solve the equations, transformations of all dependent variables and the one independent variable are carried out. The processes of interest can be decoupled into two distinctive sets of equations that describe adiabatic and nonadiabatic processes. These sets of equations are then solved by two alternating steps adiabatic and nonadiabatic