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] abus_v1_0.gz(13 Kbytes)
Manuscript Title: The computation of steady state arcs in nozzle flow.
Authors: M.T.C. Fang, S.K. Chan, R.D. Wright
Program title: DCANF
Catalogue identifier: ABUS_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 13(1977)363
Programming language: Fortran.
Computer: CDC 7600.
Operating system: SCOPE 2.1.
RAM: 12K words
Word size: 60
Keywords: Plasma physics, Arc, Electric, Gas flow, Thermal plasma, Nozzle, Discharge.
Classification: 19.5.

Nature of problem:
DCANF finds the supercritical solution of a high pressure arc burning in a converging-diverging nozzle where the primary standard shape factors are assumed constant and the external flow one-dimensional.

Solution method:
The three first-order ordinary differential equations, [Q]dx/dzeta= b, describing the d.c. arcs in nozzle flow are solved for a given solution parameter, K assuming constant standard primary shape factors. To obtain the supercritical solution, it is necessary to ensure that |Q| = 0 and the compatibility conditions are satisfied simultaneously at the critical point, zeta c. This is ensured by guessing the arcing conditions at the critical point (i.e. the velocity magnitude factor gc, and Beta c, the ratio of the arc displacement area to the nozzle area). The governing equations are integrated forward from the upstream electrode and backward from the critical point to the socalled patching point. The sum of the squares of the differences of the dependent variables at the patching point is then minimised with respect gc and Beta c using Powell's method.

Running time:
With reasonable guesses for gc and Beta c, a solution can usually be found within 1.5s on CDC7600.