Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

model to keep q>1 in stationary transport #785

Open
wants to merge 5 commits into
base: master
Choose a base branch
from
Open

model to keep q>1 in stationary transport #785

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
faae825
initial setup for raising q>1 for stationary transport actor
jmcclena Dec 12, 2024
ffc8ee9
add back in refreezes
jmcclena Dec 12, 2024
c708605
clean cruft
jmcclena Dec 12, 2024
b320584
Merge branch 'master' into steady_sawtooth
jmcclena Dec 12, 2024
5c5bd96
small updates to qed jardin sawteeth
orso82 Dec 16, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion src/actors/current/current_actor.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ Base.@kwdef mutable struct FUSEparameters__ActorCurrent{T<:Real} <: ParametersAc
_parent::WeakRef = WeakRef(nothing)
_name::Symbol = :not_set
_time::Float64 = NaN
model::Switch{Symbol} = Switch{Symbol}([:SteadyStateCurrent, :QED, :none], "-", "Current actor to run"; default=:SteadyStateCurrent)
model::Switch{Symbol} = Switch{Symbol}([:SteadyStateCurrent, :QED, :none], "-", "Current actor to run"; default=:QED)
allow_floating_plasma_current::Entry{Bool} = Entry{Bool}("-", "Zero loop voltage if non-inductive fraction exceeds 100% of the target Ip"; default=true)
#== data flow parameters ==#
ip_from::Switch{Symbol} = switch_get_from(:ip)
Expand Down
21 changes: 18 additions & 3 deletions src/actors/current/qed_actor.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,6 +14,7 @@ Base.@kwdef mutable struct FUSEparameters__ActorQED{T<:Real} <: ParametersActor{
#== data flow parameters ==#
ip_from::Switch{Symbol} = switch_get_from(:ip)
vloop_from::Switch{Symbol} = switch_get_from(:vloop)
qmin_desired::Entry{Float64} = Entry{Float64}("-", "Desired minimum magnitude of q-profile"; default=1.0)
end

mutable struct ActorQED{D,P} <: SingleAbstractActor{D,P}
Expand Down Expand Up @@ -120,8 +121,20 @@ function _step(actor::ActorQED)
Ip = nothing
Vedge = IMAS.get_from(dd, Val{:vloop}, par.vloop_from)
end
B0 = eqt.global_quantities.vacuum_toroidal_field.b0

actor.QO = QED.steady_state(actor.QO, η_imas(dd.core_profiles.profiles_1d[]); Vedge, Ip)
j_total = QED.JB(actor.QO; ρ=cp1d.grid.rho_tor_norm) ./ B0

# using Jardin's model for stationary sawteeth to raise q>1
# requires running QED twice (the first time to evaluate q, the second to relax using updated q profile)
if par.qmin_desired > 0
qval = 1.0 ./ abs.(actor.QO.ι.(cp1d.grid.rho_tor_norm))
i_qdes = findlast(qval .< par.qmin_desired)

actor.QO = QED.steady_state(actor.QO, η_imas(dd.core_profiles.profiles_1d[]; i_qdes); Vedge, Ip)
j_total = QED.JB(actor.QO; ρ=cp1d.grid.rho_tor_norm) ./ B0
end
end

return actor
Expand Down Expand Up @@ -191,8 +204,10 @@ function qed_init_from_imas(eqt::IMAS.equilibrium__time_slice, cp1d::IMAS.core_p
return QED.initialize(rho_tor, B0, gm1, f, dvolume_drho_tor, q, j_tor, gm9; ρ_j_non_inductive, ρ_grid)
end

function η_imas(cp1d::IMAS.core_profiles__profiles_1d; use_log::Bool=true)
rho = cp1d.grid.rho_tor_norm
function η_imas(cp1d::IMAS.core_profiles__profiles_1d; use_log::Bool=true, i_qdes::Int=0)
η = 1.0 ./ cp1d.conductivity_parallel
return QED.η_FE(rho, η; use_log)
if i_qdes != 0
η[1:i_qdes] .= η[i_qdes]
end
return QED.η_FE(cp1d.grid.rho_tor_norm, η; use_log)
end
Loading