.PowerSystems.Common.IronSaturation.saturationAnalytic

Information

Returns analytic saturation function.

  Input: psi0     unsaturated flux pu
  Output: psi     saturating flux as function of psi0

  psi = (c[1]/c[2])*tanh(c[2]*psi0) + c[3]*psi0

  xratio      ratio saturated/unsaturated coupling impedance
  psi_sat     saturation value of flux pu

  c[1] = 1-xratio
  c[2] = (1-xratio)/(psi_sat-xratio)
  c[3] = xratio

Interface

function saturationAnalytic
  input Real[:] psi0 "unsaturated flux pu";
  input Real[3] c;
  output Real[size(psi0, 1)] psi "n=0: saturated flux, n=1: derivative d_psi/d_psi0";
end saturationAnalytic;