Computer Programs

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES

[ top ]

[ top ]

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

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.

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 ]

[ 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.

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 ]

[ top ]

[ 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).

- 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.

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 ]

NEA-1140/01

MVS/XA (IBM 3038).[ top ]

[ top ]

[ top ]

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 |

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