RSICC CODE PACKAGE CCC-505

1. NAME AND TITLE

MEDUSA-IB: One-Dimensional Lagrangian Code for Plasma Hydrodynamic Analysis of a Fusion Pellet Driven by Ion Beams.

The plasma hydrodynamic routine of MEDUSA-IB is based on MEDUSA code for laser fusion by J. P. Christiansen et al. A preliminary version, MEDUSA-LIB was developed for the analysis of a hollow single shell target driven by proton beams. MEDUSA-IB is an upgraded version of MEDUSA-LIB for any ion beam species. It also includes multigroup X-ray, -particle and neutron transport processes etc.

2. CONTRIBUTOR

University of Tokyo, Tokyo, Japan.

3. CODING LANGUAGE AND COMPUTER

Fortran 77; HITAC-M-200H (C00505HM20000).

4. NATURE OF PROBLEM SOLVED

The MEDUSA-IB code performs implosion and thermonuclear burn calculations of an ion beam driven ICF target, based on one-dimensional plasma hydrodynamics and transport theory. It can calculate the following values in spherical geometry through the progress of implosion and fuel burnup of a multi-layered target. (1) Hydrodynamic velocities, density, ion, electron and radiation temperature, radiation energy density, Rs and burn rate of target as a function of coordinates and time, (2) Fusion gain as a function of time, (3) Ionization degree, (4) Temperature dependent ion beam energy deposition, (5) Radiation, -particle and neutron spectra as a function of time.

5. METHOD OF SOLUTION

MEDUSA-IB code is based on the one-dimensional Lagrangian hydrodynamic code MEDUSA for implosion and thermonuclear burn calculations. Collision probability method is used to solve neutron transport. The subroutines for this process has been taken from the MEDUSA-PIJ (CCC-349) code. The variable Eddington method is adopted to solve the multigroup radiation transport. The particle tracking method is used to solve both ion beam and -particle transport. These transport processes can be considered only in spherical geometry. Implosion and thermonuclear burn calculations are carried out by performing these processes with controlled time steps.

6. RESTRICTIONS OR LIMITATIONS

The maximum numbers of both X-ray and -particle energy groups for the transport calculation are 20. The maximum neutron energy groups is 22. Total spatial mesh size is below 150. These limitations can be modified easily.

7. TYPICAL RUNNING TIME

The sample problem included in the code package requires about 10 minutes by HITAC-M-280H with NOIAP (No Integrated Array Processor). Level of optimization should be OPT=2 or OPT=1.

8. COMPUTER HARDWARE REQUIREMENTS

MEDUSA-IB is designed and operable in HITAC M - series computers, such as HITAC-M-150H, M-200H and M-280H. Core memory storage requirements for the sample problem are below 1.3 MB. Maximum auxiliary storage requirements are 8 storage devices in addition to the standard input and output devices.

9. COMPUTER SOFTWARE REQUIREMENTS

The code is written in FORTRAN 77. The compilers must accept NAMELIST format, because it is used in input scheme. A system dependent subroutine CLOCK is called to sample CPU time.

10. REFERENCES

M. Uchida, Y. Oka, and S. An, "MEDUSA-IB: One Dimensional Implosion and Burnup Calculation Code for Ion Beam Driven Inertial Confinement Fusion Target," UTNL-R-0168 (October 1984).

M. Uchida, et al., "Simulation of Multigroup X-ray Alpha-Particle and Neutron Transport in Ion Beam Driven ICF Target," J. Nucl. Sci. Technol, 22(9), 683-696 (1985).

11. CONTENTS OF CODE PACKAGE

Included in the package are the referenced documents and one DOS diskette.

12. DATE OF ABSTRACT

December 1986.

KEYWORDS: CTR; HYDRONAMICS; BURNUP