Home
last modified time | relevance | path

Searched refs:Xqp (Results 1 – 5 of 5) sorted by relevance

/petsc/src/ts/tutorials/power_grid/stability_9bus/
H A Dex9busdmnetwork.c35 PetscScalar Xqp; /* q-axis transient reactance */ member
131 static const PetscScalar Xqp[3] = {0.0969, 0.1969, 0.25}; /* q-axis transient reactance */ in read_data() local
215 gen[i * 3 + j].Xqp = Xqp[j]; in read_data()
426 Edp = Vd + gen->Rs * Id - gen->Xqp * Iq; in SetInitialGuess()
605 PetscScalar Xd, Xdp, Td0p, Xq, Xqp, Tq0p, TM, D, M, Rs; /* Generator parameters */ in FormIFunction() local
623 Xqp = gen->Xqp; in FormIFunction()
635 farr[offsetgen + 1] = (Edp - (Xq - Xqp) * Iq) / Tq0p + xdotarr[offsetgen + 1]; in FormIFunction()
637 …farr[offsetgen + 3] = (-TM + Edp * Id + Eqp * Iq + (Xqp - Xdp) * Id * Iq + D * (w - W_S)) / M + xd… in FormIFunction()
642 det = Rs * Rs + Xdp * Xqp; in FormIFunction()
645 Zdq_inv[1] = Xqp / det; in FormIFunction()
[all …]
H A Dex9busadj.c38 const PetscScalar Xqp[3] = {0.0969, 0.1969, 0.25}; /* q-axis transient reactance */ variable
160 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
254 fgen[idx + 1] = (Edp - (Xq[i] - Xqp[i]) * Iq) / Tq0p[i]; in ResidualFunction()
256 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / … in ResidualFunction()
264 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
267 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
488 val[1] = -(Xq[i] - Xqp[i]) / Tq0p[i]; in ResidualJacobian()
509 val[3] = (Edp + (Xqp[i] - Xdp[i]) * Iq) / M[i]; in ResidualJacobian()
510 val[4] = (Eqp + (Xqp[i] - Xdp[i]) * Id) / M[i]; in ResidualJacobian()
517 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
[all …]
H A Dex9busoptfd.c35 const PetscScalar Xqp[3] = {0.0969, 0.1969, 0.25}; /* q-axis transient reactance */ variable
162 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
256 fgen[idx + 1] = (Edp - (Xq[i] - Xqp[i]) * Iq) / Tq0p[i]; in ResidualFunction()
258 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / … in ResidualFunction()
265 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
268 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
489 val[1] = -(Xq[i] - Xqp[i]) / Tq0p[i]; in ResidualJacobian()
510 val[3] = (Edp + (Xqp[i] - Xdp[i]) * Iq) / M[i]; in ResidualJacobian()
511 val[4] = (Eqp + (Xqp[i] - Xdp[i]) * Id) / M[i]; in ResidualJacobian()
518 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
[all …]
H A Dex9bus.c65 const PetscScalar Xqp[3] = {0.0969, 0.1969, 0.25}; /* q-axis transient reactance */ variable
386 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
482 fgen[idx + 1] = (-Edp + (Xq[i] - Xqp[i]) * Iq) / Tq0p[i]; in ResidualFunction()
484 …fgen[idx + 3] = (TM[i] - Edp * Id - Eqp * Iq - (Xqp[i] - Xdp[i]) * Id * Iq - D[i] * (w - w_s)) / M… in ResidualFunction()
491 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
494 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
747 val[1] = (Xq[i] - Xqp[i]) / Tq0p[i]; in ResidualJacobian()
768 val[3] = (-Edp - (Xqp[i] - Xdp[i]) * Iq) / M[i]; in ResidualJacobian()
769 val[4] = (-Eqp - (Xqp[i] - Xdp[i]) * Id) / M[i]; in ResidualJacobian()
776 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualJacobian()
[all …]
H A Dex9busopt.c45 const PetscScalar Xqp[3] = {0.0969, 0.1969, 0.25}; /* q-axis transient reactance */ variable
199 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
287 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in InitialGuess()
401 fgen[idx + 1] = (Edp - (Xq[i] - Xqp[i]) * Iq) / Tq0p[i]; in ResidualFunction()
403 …fgen[idx + 3] = (-TM[i] + Edp * Id + Eqp * Iq + (Xqp[i] - Xdp[i]) * Id * Iq + D[i] * (w - w_s)) / … in ResidualFunction()
410 det = Rs[i] * Rs[i] + Xdp[i] * Xqp[i]; in ResidualFunction()
413 Zdq_inv[1] = Xqp[i] / det; in ResidualFunction()
632 val[1] = -(Xq[i] - Xqp[i]) / Tq0p[i]; in ResidualJacobian()
653 val[3] = (Edp + (Xqp[i] - Xdp[i]) * Iq) / M[i]; in ResidualJacobian()
654 val[4] = (Eqp + (Xqp[i] - Xdp[i]) * Id) / M[i]; in ResidualJacobian()
[all …]