JuMP.@NLconstraint(mod, real_flow[1] == G[br_idx,n1,n1]*v[1]^2 + G[br_idx,n1,n2]*v[1]*v[2]*cos(δ) + B[br_idx,n1,n2]*v[1]*v[2]*sin(δ));
JuMP.@NLconstraint(mod, real_flow[2] == G[br_idx,n2,n2]*v[2]^2 + G[br_idx,n2,n1]*v[1]*v[2]*cos(δ) + B[br_idx,n2,n1]*v[1]*v[2]*sin(-δ));
JuMP.@NLconstraint(mod, imag_flow[1] == -B[br_idx,n1,n1]*v[1]^2 - B[br_idx,n1,n2]*v[1]*v[2]*cos(δ) + G[br_idx,n1,n2]*v[1]*v[2]*sin(δ));
JuMP.@NLconstraint(mod, imag_flow[2] == -B[br_idx,n2,n2]*v[2]^2 - B[br_idx,n2,n1]*v[1]*v[2]*cos(δ) + G[br_idx,n2,n1]*v[1]*v[2]*sin(-δ));
JuMP.@NLconstraint(mod, cur_mag_sqr == (g^2 + b^2)*(v[1]^2/tm^2 + v[2]^2 - 2*(tr*v[1]*v[2]*cos(δ) + ti*v[1]*v[2]*sin(δ))/tm^2)
- ym_sh_sqr*(v[1]^2/tm^2) + 2*(g_fr*real_flow[1] - b_fr*imag_flow[1]));
JuMP.@NLconstraint(mod, real_flow[1]^2 + imag_flow[1]^2 ≤ v[1]^2/tm^2*cur_mag_sqr);
JuMP.@NLconstraint(mod, real_flow[1]^2 + imag_flow[1]^2 ≤ branch["rate_a"]^2);
JuMP.@NLconstraint(mod, real_flow[2]^2 + imag_flow[2]^2 ≤ branch["rate_a"]^2);
It seems me that some nonlinear constraints in my code is not supported by SCIP.jl
The nonlinear constraints are like the following:
Please see the following error message: