MHD simulations on an unstructured mesh

Conference · Tue Dec 31 04:00:00 EST 1996

OSTI ID:489525

Strauss, H R; Park, W

We describe work on a full MHD code using an unstructured mesh. MH3D++ is an extension of the PPPL MH3D resistive full MHD code. MH3D++ replaces the structured mesh and finite difference / fourier discretization of MH3D with an unstructured mesh and finite element / fourier discretization. Low level routines which perform differential operations, solution of PDEs such as Poisson`s equation, and graphics, are encapsulated in C++ objects to isolate the finite element operations from the higher level code. The high level code is the same, whether it is run in structured or unstructured mesh versions. This allows the unstructured mesh version to be benchmarked against the structured mesh version. As a preliminary example, disruptions in DIIID reverse shear equilibria are studied numerically with the MH3D++ code. Numerical equilibria were first produced starting with an EQDSK file containing equilibrium data of a DIII-D L-mode negative central shear discharge. Using these equilibria, the linearized equations are time advanced to get the toroidal mode number n = 1 linear growth rate and eigenmode, which is resistively unstable. The equilibrium and linear mode are used to initialize 3D nonlinear runs. An example shows poloidal slices of 3D pressure surfaces: initially, on the left, and at an intermediate time, on the right.

• DOE Contract Number: FG02-86ER53223; AC02-76CH03073
• OSTI ID: 489525
• Report Number(s): CONF-960354--
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
99 GENERAL AND MISCELLANEOUS
DOUBLET-3 DEVICE
M CODES
MESH GENERATION
MHD EQUILIBRIUM
PLASMA DISRUPTION
PLASMA SIMULATION