Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] absg_v1_0.gz(8 Kbytes)
Manuscript Title: HEATER: a 2D laser propagation subroutine for underdense plasmas.
Authors: J.N. McMullin, C.E. Capjack, C.R. James
Program title: HEATER
Catalogue identifier: ABSG_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 23(1981)31
Programming language: Fortran.
Computer: AMDAHL 470/V7.
RAM: 37K words
Word size: 32
Keywords: Laser physics, Plasma physics, Laser propagation, Laser heating, Refraction, Diffraction, Absorption, Cubic spline, Inertial confinement.
Classification: 15, 19.7.

Nature of problem:
The physical problem is to describe the propagation of a laser beam in a cylindrically symmetric, underdense plasma given the electron density, electron temperature, ionization number and characteristics of the incoming beam. The solution is to be calculated in a subroutine which can be used with 2D simulation codes such as CASTOR 2.

Solution method:
The subroutine defines its own 2D mesh based on the mesh of the main simulation code and uses a cubic spline technique to solve the propagation equation of the laser.

Cylindrical symmetry is assumed for the plasma and the beam, and the incoming beam is assumed to have a Gaussian cross-section. The approximations used to derive the propagation equation break down when the density is close to the critical density. Anomalous effects are not included.

Running time:
If the internal mesh is calculated, one call to the subroutine requires approximately 0.7 s on the Amdahl 470/V7 (approximately equivalent to the IBM 3033). Approximately 0.6 s are required if the mesh has been calculated in a previous call. These times will increase if the internal mesh is made finer in an attempt to gain more accuracy but do not depend on the coarseness of the external mesh of the calling program.