NEA-1140 MEDUSA-1B.

1. NAME OR DESIGNATION OF PROGRAM:  MEDUSA-1B.

3. DESCRIPTION OF PROGRAM OR FUNCTION

The MEDUSA-1B code performs implosion and thermonuclear burn calculations of an ion beam driven  ICP 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, 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, alpha-particle and neutron spectra as a function of     time

4. METHOD OF SOLUTION

MEDUSA-IB code is based on the one-dimensional  Lagrangian hydrodynamic code MEDUSA for implosion and thermo- nuclear burn calculations. The collision probability method is used  to solve neutron transport. The subroutines for this process have be taken from the code MEDUSA-PIJ. The variable Eddington method was adopted to solve the multigroup radiation transport. The particl tracking method is used to solve both ion beam and alpha-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.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

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

6. TYPICAL RUNNING TIME

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

NEA-DB ran the test case included in this package on an IBM 3083 computer in 63 seconds of CPU time.

7. UNUSUAL FEATURES OF THE PROGRAM

MEDUSA-IB is an upgraded version of MEDUSA-LIB for any ion beam species. It also includes multigroup  X-ray, alpha-particle and neutron transport processes etc.

8. RELATED AND AUXILIARY PROGRAMS:  MEDUSA, MEDUSA-LIB.

10. REFERENCES

- J.P. Christiansen et al.
  "MEDUSA, A One-Dimensional Laser Fusion Code" Computer Physics
  Communication, 7 271 (1974).
- "UTLIF(2)" Chapter 4,
  Nuclear Engineering Research Laboratory, University of Tokyo,
  UTNL-R-0150.
- M. Uchida et al.
Parametric Study of Light Ion Beam Fusion Target Uusing MEDUSA-LIB    and MEDUSA-IB, INS International Symposium on Heavy Ion
  Accelerators and their Applications to ICF, Jan. 23-27, 1984,
  Institute of Nuclear Study, University of Tokyo, Tokyo.
- H. Takano and Y. Ishiguro
  "MEDUSA-PIJ: A Code for One-Dimensional Laser Fusion Analysis
  Taking Account of Neutron Heating Effect",
  JAERI-M-8186 (1979).
- P.M. Cambell and J.J. Kubis
  A Variable Eddington Method for Radiation Transport in Dense
  Fusion Plasmas, KMS Fusion, Inc.,
  Report KMSF-U458 (Jan. 1976).
- Moses, Nucl. Sci. Eng. 64 49 (1977).

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

11. MACHINE REQUIREMENTS

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

To run the test case on an IBM 3083, 1036K bytes of main storage were required.

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  HITAC VOS 3 (OS like MVS on IBM machines).

MVS/XA (IBM 3038).

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

The compilers should be able to handle NAMELIST format, because it is used in the input scheme. A subroutine CLOCK is called to sample CPU time. It should be replaced om other systems by an appropriate subroutine

15. NAME AND ESTABLISHMENT OF AUTHORS

     Mitsuhiko UCHIDA, Yoshiaki OKA, Shigehiro AN
     Nuclear Engineering Research Laboratory
     Faculaty of Engineering
     Tokai-mura, Ibaraki
     JAPAN

File name	File description	Records
NEA1140_01.001	Information file	87
NEA1140_01.002	MEDUSA FORTRAN source	14652
NEA1140_01.003	JCL used at Data Bank	31
NEA1140_01.004	Ionization potential data set	27
NEA1140_01.005	Neutron elastic scattering cross section	20
NEA1140_01.006	Sample input data of MEDUSA	113
NEA1140_01.007	Sample output list	4645