Create structure and evolution for build currents #18
Conversation
|
@lstagner The goal here is to keep QED very lightweight. We should build the inductance matrix, resistivity vector, and plasma-coil mutual/mutual-derivative vectors externally (from an IMAS dd or TokSys model), then just pass it to a I've made the After this, we need to couple it to the diffusion inside the plasma, which amounts to combining the matrix equations defined in |
|
As of this commit, the build evolution gives the same result as @lstagner's DifferentialEquations.jl implementation. |
|
7a2ecca is untested, but should evolve the coupled plasma-build system. |
- Move plasma boundary condition in coupled problem to after adding non-inductive source - Renormalize build part of A matrix to improve condition number - Debug outputs plasma current with coil currents - Add some examples
|
As of the last commit, the plasma-as-coil boundary condition produces the same results for the plasma and coil currents as a system where an extra coil is added with the same L, R, and V_noninductive as the plasma. If L, R, and V_ni are correct, then this should work. |
bclyons12
changed the title
This will allow for the evolution of currents in active & passive conducting structures (hereafter all referred to as coils) with implicit coupling to the changing plasma current and explicit coupling to the changing plasma shape (through provided, fixed time-derivatives of the plasma-coil mutuals).