Computer Programs
NEA-1140 MEDUSA-1B.
last modified: 11-OCT-1989 | catalog | categories | new | search |

NEA-1140 MEDUSA-1B.

MEDUSA-1B, 1-D Plasma Hydrodynamic Analysis of Fusion Pellet Ion Beams

top ]
1. NAME OR DESIGNATION OF PROGRAM:  MEDUSA-1B.
top ]
2. COMPUTERS
To submit a request, click below on the link of the version you wish to order. Rules for end-users are available here.
Program name Package id Status Status date
MEDUSA-1B NEA-1140/01 Tested 11-OCT-1989

Machines used:

Package ID Orig. computer Test computer
NEA-1140/01 HITAC M-200 H IBM 3083
top ]
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
top ]
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.
top ]
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.
top ]
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-1140/01
NEA-DB ran the test case included in this package on an IBM 3083 computer in 63 seconds of CPU time.
top ]
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.
top ]
8. RELATED AND AUXILIARY PROGRAMS:  MEDUSA, MEDUSA-LIB.
top ]
9. STATUS
Package ID Status date Status
NEA-1140/01 11-OCT-1989 Tested at NEADB
top ]
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).
NEA-1140/01, included references:
- 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)
top ]
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.
NEA-1140/01
To run the test case on an IBM 3083, 1036K bytes of main storage were required.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NEA-1140/01 FORTRAN-77
top ]
13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  HITAC VOS 3 (OS like MVS on IBM machines).
NEA-1140/01
MVS/XA (IBM 3038).
top ]
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
top ]
15. NAME AND ESTABLISHMENT OF AUTHORS

     Mitsuhiko UCHIDA, Yoshiaki OKA, Shigehiro AN
     Nuclear Engineering Research Laboratory
     Faculaty of Engineering
     Tokai-mura, Ibaraki
     JAPAN
top ]
16. MATERIAL AVAILABLE
NEA-1140/01
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
top ]
17. CATEGORIES
  • X. Magnetic Fusion Research

Keywords: fusion reactions, fusion reactors, hydrodynamics, ion beams, lasers, plasma, thermonuclear reactions, two-dimensional.