Abstract
Disruption and vertical instability are simulated by use of the multi-layers method. The tokamak system is not closed in the energy since it has an energy source to contorol the coil current. Then plasma equilibrium is found by minimizing the energy functional, which is free energy including the energy flow from the power supply. Time evolution of the system is presumed the continuous equilibrium. Because this model, which includes no cold plasma surrounding the main plasma, can take longer computational time step than MHD model, it enables us to design a PF coil system incorporating the transient phenomenon of plasma such as disruption, and to study the efficient operation and the reliable control methods of next-generation tokamak-type reactors.