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. Only liaison officers are authorised to submit online requests. Rules for requesters 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.