Load Flow Analyzer
Perform positive-sequence load flow or unbalanced load flow and initialize models containing load flow blocks
Since R2021a
Open the Load Flow Analyzer App
MATLAB® command prompt: Enter
powerLoadFlowpowergui Block Parameters dialog box: On the Tools tab, click Load Flow Analyzer.
To perform a load flow analysis and initialize your model so that it starts in steady state:
Define the model buses using Load Flow Bus blocks.
Specify the load flow parameters of all blocks that have load flow parameters. These blocks are referred to as load flow blocks.
Solve the load flow and interactively modify the load flow parameters until a satisfactory solution is obtained.
Save the load flow parameters and machine initial conditions in the model.
Examples
Positive-Sequence Load Flow
Open the Initializing a 5-Bus Network with the Load Flow Tool of Powergui example to open a model containing five Load Flow Bus blocks and six load flow blocks.
The Load Flow Bus blocks are shown in orange and the load flow blocks are shown in yellow.
The Load Flow Bus blocks specify the bus base voltages (nominal phase-to-phase rms voltage). They also specify the voltage at PV buses or the voltage and angle of the swing buses. Once the load flow is solved, the Load Flow Bus block displays the bus positive-sequence voltage magnitude and phase angle as block annotations.
The bus type (PV, PQ, or swing) is determined by the load flow blocks connected to the bus. If you have several load flow blocks with different types (specified in the Generator type parameter or in the Load type parameter) connected to the same bus, the Load Flow tool determines the resulting bus type (swing, PQ, or PV).
In the Initializing a 5-Bus Network with the Load Flow Tool of Powergui example, the bus types are determined as follows:
| Bus | Load Flow Blocks | Resulting Bus Type |
|---|---|---|
| 120 kV Three-Phase Source | swing Specify the voltage and angle in the B120 Load Flow Bus block. |
| 13.8 kV 150 MVA Synchronous Machine 3 MW 2 Mvar RLC Load | PV Specify the voltage in the B13.8 Load Flow Bus block. |
| 10 MW, 3 Mvar Dynamic Load | PQ |
| No load flow blocks | PQ |
| Asynchronous generator 9 MW | PQ Constant Z load is included in the Ybus admittance matrix. |
Some restrictions apply when you connect several source blocks and synchronous machines at the same bus:
Two swing generators cannot be connected in parallel.
A swing generator cannot be connected in parallel with a PV ideal voltage source.
When a swing voltage source with RL impedance is connected to a PV generator, the swing bus is automatically moved to the ideal voltage source connection node, behind the RL source impedance.
Only one PV generator with finite Q limits can be connected at a generation bus. However, you may have other PQ generators and loads connected on the same bus.
For more information on how to use the Load Flow Bus block in your model, see the Load Flow Bus page.
Using the Load Flow Tool to Perform Load Flow Analysis
Once you have entered the load flow parameters in the Load Flow Bus blocks and in the various load flow blocks, open the Load Flow Analyzer by clicking the Load Flow Analyzer button of the powergui block. The tool displays a summary of the load flow data of the model. The table below shows the data found in the Initializing a 5-Bus Network with the Load Flow Tool of Powergui example.
Note that the table contains seven lines, but there are only six load flow blocks in the model. This is because the bus B25_2 is not connected to any load flow block. Line 5 is added in the table for that particular bus, so that you can see all buses listed together with their bus voltages. This bus will be considered in the load flow analysis as a PQ bus with zero P and Q.
The Block name column identifies the block type. The Block type column displays the bus type of the load flow blocks. The next four columns give the bus identification label, the bus base voltage, the reference voltage (in pu of base voltage), and the voltage angle of the load flow bus where the block is connected. The following columns are the P and Q values specified in the Load Flow tab of the blocks.
The last four columns display the current load flow solution. Because the load flow has not been performed, the columns display zero values.
The load flow parameters in the Preferences tab of the powergui block are used to build the Ybus network admittance matrix and to solve the load flow. The base power is used to specify units of the normalized Ybus matrix in pu/Pbase and bus base voltages. The Initializing a 5-Bus Network with the Load Flow Tool of Powergui model contains five buses; consequently, the Ybus matrix will be a 5-by-5 complex matrix evaluated at the frequency specified by the Frequency (Hz) parameter.
The load flow algorithm uses an iterative solution based on the Newton-Raphson method. The Max iterations parameter defines the maximum number of iterations. The load flow algorithm will iterate until the P and Q mismatch at each bus is lower than the PQ tolerance parameter (in pu/Pbase). The power mismatch is defined as the difference between the net power injected into the bus by generators and PQ loads and the power transmitted on all links leaving that bus.
To avoid a badly conditioned Ybus matrix, you should select a Base power parameter value in the range of nominal powers and loads connected to the network. For a transmission network with voltages ranging from 120 kV to 765 kV, a 100 MVA base is usually selected. For a distribution network or for a small plant consisting of generators, motors, and loads that have a nominal power in the range of hundreds of kilowatts, a 1 MVA power base is better adapted.
To solve the load flow, click the Compute button. The load flow solution is displayed in the last five columns of the table.
To display the load flow report showing power flowing at each bus, click the Report button. Save this report in a file by specifying the file name at the prompt.
The report displays the summary of active and reactive powers, including the total PQ sharing between generators (SM- and Vsrc-type blocks), PQ loads (PQ-type RLC loads and DYN loads), shunt constant Z loads (Z-type RLC loads and magnetizing branches of transformers), and asynchronous machine loads (ASM):
The Load Flow converged in 2 iterations !
SUMMARY for subnetwork No 1
Total generation : P= 5.61 MW Q= 25.51 Mvar
Total PQ load : P= 13.00 MW Q= 5.00 Mvar
Total Zshunt load : P= 0.68 MW Q= -0.51 Mvar
Total ASM load : P= -8.90 MW Q= 4.38 Mvar
Total losses : P= 0.83 MW Q= 16.64 MvarThe Total losses line represents the difference between generation
and loads (PQ type + Z type +ASM) and represents series losses. After this summary, a
voltage and power report is presented for each bus:
1 : B120 V= 1.020 pu/120kV 0.00 deg ; Swing bus
Generation : P= -114.39 MW Q= 62.76 Mvar
PQ_load : P= 0.00 MW Q= 0.00 Mvar
Z_shunt : P= 0.25 MW Q= 0.23 Mvar
--> B13.8 : P= -116.47 MW Q= 53.89 Mvar
--> B25_1 : P= 1.84 MW Q= 8.63 Mvar
2 : B13.8 V= 0.980 pu/13.8kV -23.81 deg
Generation : P= 120.00 MW Q= -37.25 Mvar
PQ_load : P= 3.00 MW Q= 2.00 Mvar
Z_shunt : P= 0.17 MW Q= 0.17 Mvar
--> B120 : P= 116.83 MW Q= -39.42 Mvar
3 : B25_1 V= 0.998 pu/25kV -30.22 deg
Generation : P= 0.00 MW Q= 0.00 Mvar
PQ_load : P= 10.00 MW Q= 3.00 Mvar
Z_shunt : P= 0.25 MW Q= 0.21 Mvar
--> B120 : P= -1.83 MW Q= -8.44 Mvar
--> B25_2 : P= -8.41 MW Q= 5.23 Mvar
4 : B25_2 V= 0.967 pu/25kV -20.85 deg
Generation : P= 0.00 MW Q= 0.00 Mvar
PQ_load : P= -0.00 MW Q= -0.00 Mvar
Z_shunt : P= 0.01 MW Q= -0.03 Mvar
--> B25_1 : P= 8.87 MW Q= -3.67 Mvar
--> B575 : P= -8.88 MW Q= 3.70 Mvar
5 : B575 V= 0.953 pu/0.575kV -18.51 deg
Generation : P= 0.00 MW Q= 0.00 Mvar
PQ_load : P= -0.00 MW Q= -0.00 Mvar
Z_shunt : P= 0.01 MW Q= -1.09 Mvar
--> ASM : P= -8.90 MW Q= 4.38 Mvar
--> B25_2 : P= 8.89 MW Q= -3.29 MvarFor every bus, the bus voltage and angle are listed on the first line. The next three lines give the PQ generated at the bus (all SM and voltage sources), the PQ absorbed by the PQ type loads, and the PQ absorbed by the Z-type loads.
The last lines, preceded by an arrow (-->), list the PQ
transmitted to neighbor buses connected through lines, series impedances, and transformers,
as well power absorbed by ASM.
Apply the Load Flow Solution to Your Model
When performing a load flow analysis, you may need to iterate on P, Q, and V values until you find satisfactory voltages at all buses. This may require, for example, changing generated power, load powers, or reactive shunt compensation.
To change the load flow setup, you need to edit the parameters of the load flow blocks and the Load Flow Bus blocks. Then click the Update button to refresh the load flow data displayed by the table in the Load Flow Analyzer. The previous load flow solution is then deleted from the table. Click the Compute button to obtain a new load flow solution that corresponds to the changes you made.
Once you have obtained a satisfactory load flow, update the model initial conditions according to the load flow solution. Click the Apply to Model button to initialize the machine blocks of the model, and as the initial conditions of regulators connected to the machines.
Open the Three-Phase Parallel RLC Load block connected at the B13.8 bus.
Because the Load type specified in the Load Flow
tab is constant PQ, the nominal voltage of this block has been changed to the corresponding
bus voltage of 0.98 pu. The Nominal phase-to-phase voltage Vn (Vrms)
parameter is set to (13800)*0.98.
Open the Three-Phase Dynamic Load block connected at the B25_1bus. The
Initial positive-sequence voltage Vo [Mag(pu) Phase (deg.) is set to
[0.998241 -30.2228].
Note that the voltage magnitudes and angles obtained at each bus are written as block annotations under the Load Flow Bus blocks.
Open the Scope block and start the simulation.
The Three-Phase Fault block applies a six-cycle fault at the B120 bus.
Observe the waveforms of the SM active power, SM and ASM speeds, and PQ of DYN load, and notice that simulation starts in steady state.
Example of Unbalanced Load Flow
Open the IEEE 13 Node Test Feeder to open a model containing 12 Load Flow Bus blocks and 13 load flow blocks. This model is a benchmark network taken from the Radial Distribution Test Feeder in the Distribution System Analysis Subcommittee Report by the Power Engineering Society in pages 908–912 written in 2001.
The original benchmark system contains 13 nodes. However, because the IEEE 13 Node Test Feeder model does not include the regulating transformer, it contains only 12 nodes.
The Load Flow Bus blocks are shown in orange and the Load Flow blocks are shown in yellow.
The Load Flow Bus blocks specify the bus base voltages (nominal phase-to-ground rms voltage). They specify the voltage at PV buses or the voltage and angle of the swing buses. Once the load flow is solved, the Load Flow Bus block displays the bus voltage magnitude and phase angle as block annotations.
Note
By default, the block annotations are set in the Block Annotation tab of the Load Flow Bus block properties to display the phase A magnitude (<VLF> parameter) and the phase A angle (<angleLF> parameter). To display phase B magnitude and angle, specify <VLFb> and <angleLFb>, respectively. To display phase C magnitude and angle, specify <VLFc> and <angleLFc>, respectively.
You can also delete some block annotations. In the IEEE 13 Node Test Feeder example, only the bus identification is displayed (<ID> parameter).
The bus type (PV, PQ, or swing) is determined by the load flow blocks connected to the bus. If you have several load flow blocks with different types (specified in the Generator type parameter or in the Load type parameter) connected to the same bus, the Load Flow tool determines the resulting bus type (swing, PQ, or PV). The table shows how the bus types are determined for some of the model buses of the IEEE 13 Node Test Feeder example.
| Bus | Load Flow Blocks | Resulting Bus Type |
|---|---|---|
| 4160 V swing | 632_a=swing V=1.0210 pu -2.49 deg. Voltages and angles are specified in the ‘632’ Load Flow Bus block. |
| No load flow block | PQ |
| 634 Yg PQ load block | PQ |
| 646_Z load block | PQ Constant Z loads are included in the Ybus admittance matrix. |
| 675 Yg PQ load 675 Yg Z
load | PQ Constant Z loads are included in the Ybus admittance matrix. |
Some restrictions apply when you have several source blocks and synchronous machines connected to the same load flow bus:
You cannot connect two swing generators in parallel.
You cannot connect a swing generator in parallel with a PV ideal voltage source
You can connect only one PV generator with finite Q limits at a generation bus. However, you can have other PQ generators and loads connected on the same bus.
For more information on how to use the Load Flow Bus block in your model, see Load Flow Bus block.
Open the Load Flow Analyzer by clicking the Load Flow Analyzer button in the powergui block. The tool displays a list of the individual single-phase buses (one bus per phase) found in the IEEE 13 Node Test Feeder model. In the Load Flow Analyzer, the load flow has not yet been performed, so the V_LF (pu) and Vangle_LF (deg) columns display zero values.
The load flow parameters in the Preferences tab of the powergui block are used to build the Ybus network admittance matrix and to solve the load flow. The base power is used to specify units of the normalized Ybus matrix in pu/Pbase and bus base voltages. The IEEE 13 Node Test Feeder model contains 29 single phase buses; consequently, the Ybus matrix is a 29-by-29 complex matrix evaluated at the frequency specified by the Frequency (Hz) parameter.
The load flow algorithm uses an iterative solution based on the Newton-Raphson method. The Max iterations parameter defines the maximum number of iterations. The load flow algorithm iterates until the P and Q mismatch at each bus is lower than the PQ tolerance parameter (in pu/Pbase). The power mismatch is defined as the difference between the net power injected into the bus by generators and PQ loads and the power transmitted on all links leaving that bus.
To avoid a badly conditioned Ybus matrix, select the Base power parameter value in the range of nominal powers and loads connected to the network. For a transmission network with voltages ranging from 120 kV to 765 kV, a 100 MVA base is usually selected. For a distribution network with loads that have a nominal power in the range of tens to hundreds of kVA, a 100 kVA to 1 MVA power base is better adapted.
To solve the load flow, click Compute. The bus voltages and angles appear in the V_LF (pu) and Vangle_LF (deg) columns of the table.
To display a load flow report that shows the power flow at each bus, click Report. Save this report in a file by specifying the file name at the prompt.
The report displays the summary of active and reactive powers, including the total PQ sharing between generators (SM- and Vsrc-type blocks), PQ loads (PQ-type RLC loads, dynamic loads, and asynchronous machine loads), and shunt constant Z loads (Z-type RLC loads and magnetizing branches of transformers):
SUMMARY for subnetwork No 1
Total generation : P= 3518.74 kW Q= 1540.14 kvar
Total PQ load : P= 3101.90 kW Q= 1880.42 kvar
Total Zshunt load : P= 363.47 kW Q= -479.42 kvar
Total losses : P= 53.36 kW Q= 139.14 kvar
The Total losses line represents the difference between generation
and loads (PQ type + Zshunt type) and represents series losses. After this summary, a
voltage and power report appears for each bus. For each phase of every bus, the bus
voltage and angle are listed on the first line. The next three lines give the PQ generated
at the bus (all SM and voltage sources), the PQ absorbed by the PQ type loads, and the PQ
absorbed by the Z-type loads. The last lines, preceded by an arrow (–>), list the PQ
power transmitted on all links leaving that bus.
The last column gives the positive-sequence bus voltage V1 (magnitude and angle, for three-phase buses only) and the sum of PQ powers for all phases (PQ generated by sources, PQ absorbed by loads, and PQ transmitted through transformers, lines, and series impedances). For example, you can verify that the total PQ load absorbed at bus 634 (P = 400 kW Q = 290 kvar) corresponds to the sum of active and reactive powers specified for phases A, B, and C in the load block.
When performing a load flow analysis, you might need to try different P, Q, and V values until you find satisfactory voltages at all buses. This can require, for example, changing generated power, load powers, or reactive shunt compensation.
To change the load flow setup, edit the parameters of the load flow blocks and of the Load Flow Bus blocks. Then click Update to refresh the load flow data displayed by the table. Click Compute to get a new load flow solution that corresponds to the changes you made.
Once you have a satisfactory load flow, update the model initial conditions according to the load flow solution. Click Apply to Model to initialize the PQ-type load blocks, the source block internal voltages, the machine blocks, and the initial conditions of associated regulators.
Open the Three-Phase Series RLC Load block connected at bus 632.
Because the Load type specified in the Load Flow
tab is constant PQ, the vector of Nominal phase-to-neutral voltages [Va Vb Vc]
(Vrms) of this block has been changed to the corresponding bus voltages
[1.021 1.042 1.0174]*2401.78 Vrms. Open the Three-Phase
Source block connected at bus 632. The Line-to-neutral voltages [Va Vb
Vc] (Vrms) parameter is also set to [1.021 1.042
1.0174]*2401.78 Vrms.
Open the Load Flow Results subsystem and start the simulation.
Observe voltage magnitudes and PQ powers on the Display blocks. These values correspond to values displayed in the load flow report.
Related Examples
Parameters
Version History
Introduced in R2021a
You can also select a web site from the following list
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Asia Pacific
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)