Home
last modified time | relevance | path

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

/petsc/src/ts/tutorials/power_grid/stability_9bus/
H A Dex9busdmnetwork.c33 PetscScalar Xdp; /* d-axis transient reactance */ member
129 static const PetscScalar Xdp[3] = {0.0608, 0.1198, 0.1813}; /* d-axis transient reactance */ in read_data() local
213 gen[i * 3 + j].Xdp = Xdp[j]; in read_data()
429 Eqp = Vq + gen->Rs * Iq + gen->Xdp * Id; in SetInitialGuess()
440 Efd = Eqp + (gen->Xd - gen->Xdp) * Id; in SetInitialGuess()
605 PetscScalar Xd, Xdp, Td0p, Xq, Xqp, Tq0p, TM, D, M, Rs; /* Generator parameters */ in FormIFunction() local
620 Xdp = gen->Xdp; in FormIFunction()
634 farr[offsetgen] = (Eqp + (Xd - Xdp) * Id - Efd) / Td0p + xdotarr[offsetgen]; 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()
[all …]
H A Dex9busadj.c36 const PetscScalar Xdp[3] = {0.0608, 0.1198, 0.1813}; /* d-axis transient reactance */ variable
161 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in SetInitialGuess()
179 Efd = Eqp + (Xd[i] - Xdp[i]) * Id; in SetInitialGuess()
253 fgen[idx] = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[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()
268 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
478 val[1] = (Xd[i] - Xdp[i]) / Td0p[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()
[all …]
H A Dex9busoptfd.c33 const PetscScalar Xdp[3] = {0.0608, 0.1198, 0.1813}; /* d-axis transient reactance */ variable
163 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in SetInitialGuess()
181 Efd = Eqp + (Xd[i] - Xdp[i]) * Id; in SetInitialGuess()
255 fgen[idx] = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[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()
269 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
479 val[1] = (Xd[i] - Xdp[i]) / Td0p[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()
[all …]
H A Dex9busopt.c43 const PetscScalar Xdp[3] = {0.0608, 0.1198, 0.1813}; /* d-axis transient reactance */ variable
200 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in SetInitialGuess()
218 Efd = Eqp + (Xd[i] - Xdp[i]) * Id; in SetInitialGuess()
288 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in InitialGuess()
304 Efd = Eqp + (Xd[i] - Xdp[i]) * Id; in InitialGuess()
400 fgen[idx] = (Eqp + (Xd[i] - Xdp[i]) * Id - Efd) / Td0p[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()
414 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
622 val[1] = (Xd[i] - Xdp[i]) / Td0p[i]; in ResidualJacobian()
[all …]
H A Dex9bus.c63 const PetscScalar Xdp[3] = {0.0608, 0.1198, 0.1813}; /* d-axis transient reactance */ variable
387 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in SetInitialGuess()
405 Efd = Eqp + (Xd[i] - Xdp[i]) * Id; in SetInitialGuess()
481 fgen[idx] = (-Eqp - (Xd[i] - Xdp[i]) * Id + Efd) / Td0p[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()
495 Zdq_inv[2] = -Xdp[i] / det; in ResidualFunction()
737 val[1] = -(Xd[i] - Xdp[i]) / Td0p[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()
[all …]