ShuntCompensatorReactiveFlowEquationTerm.java
/**
* Copyright (c) 2019, RTE (http://www.rte-france.com)
* This Source Code Form is subject to the terms of the Mozilla Public
* License, v. 2.0. If a copy of the MPL was not distributed with this
* file, You can obtain one at http://mozilla.org/MPL/2.0/.
* SPDX-License-Identifier: MPL-2.0
*/
package com.powsybl.openloadflow.ac.equations;
import com.powsybl.math.matrix.DenseMatrix;
import com.powsybl.openloadflow.equations.Variable;
import com.powsybl.openloadflow.equations.VariableSet;
import com.powsybl.openloadflow.network.LfBus;
import com.powsybl.openloadflow.network.LfShunt;
import java.util.List;
import java.util.Objects;
/**
* @author Geoffroy Jamgotchian {@literal <geoffroy.jamgotchian at rte-france.com>}
*/
public class ShuntCompensatorReactiveFlowEquationTerm extends AbstractShuntCompensatorEquationTerm {
private Variable<AcVariableType> bVar;
private final List<Variable<AcVariableType>> variables;
public ShuntCompensatorReactiveFlowEquationTerm(LfShunt shunt, LfBus bus, VariableSet<AcVariableType> variableSet, boolean deriveB) {
super(shunt, bus, variableSet);
if (deriveB) {
bVar = variableSet.getVariable(shunt.getNum(), AcVariableType.SHUNT_B);
variables = List.of(vVar, bVar);
} else {
variables = List.of(vVar);
}
}
@Override
public List<Variable<AcVariableType>> getVariables() {
return variables;
}
private double b() {
return bVar != null ? sv.get(bVar.getRow()) : element.getB();
}
public static double q(double v, double b) {
return -b * v * v;
}
public static double dqdv(double v, double b) {
return -2 * b * v;
}
public static double dqdb(double v) {
return -v * v;
}
@Override
public double eval() {
return q(v(), b());
}
@Override
public double der(Variable<AcVariableType> variable) {
Objects.requireNonNull(variable);
if (variable.equals(vVar)) {
return dqdv(v(), b());
} else if (variable.equals(bVar)) {
return dqdb(v());
} else {
throw new IllegalStateException("Unknown variable: " + variable);
}
}
public static double calculateSensi(double v, double b, double dv, double db) {
return dqdv(v, b) * dv + dqdb(v) * db;
}
@Override
public double calculateSensi(DenseMatrix dx, int column) {
double dv = dx.get(vVar.getRow(), column);
double db = bVar != null ? dx.get(bVar.getRow(), column) : 0;
double v = v();
double b = b();
return calculateSensi(v, b, dv, db);
}
@Override
protected String getName() {
return "ac_q_shunt";
}
}