Skip to content

Derive the candidate-line Big-M from the nodal angle bound #166

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

Merged
erikfilias merged 1 commit into from
Aug 30, 2026
+5 −2
Merged

Derive the candidate-line Big-M from the nodal angle bound #166

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
7 changes: 5 additions & 2 deletions openTEPES/openTEPES_DataConfiguration.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -938,8 +938,11 @@ def Create_ESS_RES_Sets(mTEPES) -> None:
for lea in mTEPES.lea:
par['pBigMFlowBck'].loc[lea] = par['pLineNTCBck'][lea]
par['pBigMFlowFrw'].loc[lea] = par['pLineNTCFrw'][lea]
# Delta-theta across a line is bounded by twice the nodal angle bound, so the angle side of the
# Big-M follows from pMaxTheta rather than from an independent literal.
pMaxThetaValue = math.pi / 2
for lca in mTEPES.lca:
M_angle_lca = (1.0 + pMBigMEpsilon) * max(par['pLineNTCBck'][lca], math.pi * par['pSBase'] / par['pLineX'][lca])
M_angle_lca = (1.0 + pMBigMEpsilon) * max(par['pLineNTCBck'][lca], 2.0 * pMaxThetaValue * par['pSBase'] / par['pLineX'][lca])
par['pBigMFlowBck'].loc[lca] = M_angle_lca
par['pBigMFlowFrw'].loc[lca] = M_angle_lca

Expand All @@ -948,7 +951,7 @@ def Create_ESS_RES_Sets(mTEPES) -> None:
par['pBigMFlowFrw'] = par['pBigMFlowFrw'].where(par['pBigMFlowFrw'] != 0.0, 1.0)

# maximum voltage angle
par['pMaxTheta'] = par['pDemandElec']*0.0 + math.pi/2
par['pMaxTheta'] = par['pDemandElec']*0.0 + pMaxThetaValue
par['pMaxTheta'] = par['pMaxTheta'].loc[mTEPES.psn]

# this option avoids a warning in the following assignments
Expand Down
Loading