Modeling and Simulation

Brambilla Spectrum

LH Coupling Code

Originally written by M. Brambilla[1] in Fortran, this code is well benchmarked by experiment and is used for modeling of lower hybrid wave coupling efficiency and spectra. This code was modified to use C-Mod's phase settings and produce an improved output file. Results from the code are compared directly to detailed experimental measurements.
DQL Wave Damping

CQL3D

A very complete Fokker-Planck solver for waves in plasmas written in Fortran by R.W. Harvey and M.G. McCoy[2]. It also has a coupled x-ray module for comparison to HXR camera results [3]. Requires Genray input.
DKE HXR Chords

DKE

Matlab code written by Y. Peysson and J. Decker[4] to solve the electron drift kinetic equation. In general a simpler model than CQL3D, which is also easier and faster to run. It also has a great HXR emission package that is useful for interpreting HXR camera measurements. Can use input rays from Genray, or a simple model for the damping may be input by the user.
Maxwellian Damping

Genray

Fortran code written by Smirnov[5] that solves the ray equations for plasma waves and computes their trajectories, using one of a number of dispersion relations. The ray trajectories and Maxwellian damping results are needed as input for CQL3D. It also includes a module for calculating the ECE emission.
MCNP Photon Spectra

MCNP

The standard Monte Carlo neutron and photon transport code[5] in the field of nuclear engineering. This code was used to determine the background energy spectra and flux of neutrons and photons at the location of the HXR camera. It was also used to study shield thickness requirements.
Profiles

TRANSP

TRANSP is a time dependent global plasma transport code which can be used to either analyze transport coefficients for a given experimental disharge or used to predict plasma performance. TRANSP is maintained at PPPL by D. McCune and includes a number of modules such as NCLASS, TORIC, and LSC. Advanced tokamak scenarios for C-Mod were studied and web archived.
[1] M. Brambilla Nucl. Fusion 16, 1976 p47
[2] R.W. Harvey and M.G. McCoy. IAEA TCM 1992
[3] R.W. Harvey. CompX Internal Report, CompX-2000-2, February 2001
[4] J. Decker and Y. Peysson. DKE: A Fast Numerical Solver for the 3-d Relativistic Bounce-Averaged Electron Drift Kinetic Equation. Technical Report EUR-CEA-FC-1736, EURATOM CEA-Cadarache, December 2004
[4] A.P. Smirnov and R.W. Harvey. CompX Internal Report, CompX-2000-01, Ver. 2, March 2003
[5] J. Briemeister et. al. MCNP- A General Monte Carlo N-Particle Transport Code, Version 4C, la-13709-m edition, April 2000