Electric Load Center Distribution Manager[LINK]
The electric load center distribution manager (object name: ElectricLoadCenter:Distribution) operates on-site generators specified in the simulation and reports the amount of generated and purchased electricity. Electric load centers really serve as a “load” to the generators and “supply” power to the rest of the building. The internal meters used by EnergyPlus for reporting do all of the demand tracking. For more details on the individual inputs required see the Input Output Reference document.
The electric load center manager sums all of the building and system electric loads and provides operating schemes for the generators. The electric load center objects are operated in the order they are defined in the input data file (IDF), and generators are dispatched sequentially in the order they are specified within each load center object. The electricity produced from photovoltaic arrays is handled somewhat separately and is always “used” to reduced the demand that the generators will try to meet for that time step. What is not provided by the on-site generation equipment, and electric storage units if specified, is met by (purchasing) off-site electricity.
The available operating schemes are “Baseload”, “DemandLimit”, “TrackElectrical,” “TrackSchedule,” “TrackMeter,” “FollowThermal” and “FollowThermalLimitElectrical.” These operating schemes affect how loads are dispatched to the generators, in effect telling the generators whether or not to run and requesting power levels.
The BASELOAD scheme operates the generators at their rated (requested) electric power output when the generator is scheduled ON (ref. ElectricLoadCenter:Generators in the Input Output Reference). The Baseload scheme requests all generators scheduled ON (available) to operate, even if the amount of electric power generated exceeds the total facility electric power demand.
The DEMAND LIMIT scheme limits the amount of purchased electricity from the utility to the amount specified in the input. The Demand Limit scheme tries to have the generators meet all of the demand above the purchased electric limit defined by the user.
The TRACK ELECTRICAL scheme tries to have the generators meet all of the electrical demand for the building.
The TRACK METER scheme tries to have the generators meet all the electrical demand from a meter chosen by the user rather than the usual meter for the entire facility. The meter can be a custom meter so that generators are tied to only certain loads in the building.
The TRACK SCHEDULE scheme tries to have the generators meet all of the electrical demand determined by a user-defined schedule.
The FOLLOW THERMAL and FOLLOW THERMAL LIMIT ELECTRICAL schemes run the generators to meet thermal demand. The thermal demand is determined from the plant modeling and depends on the flow requested by other components on the demand side of the plant loop, the loop temperatures, and the loop temperature setpoint. The electric load center distribution manager converts the thermal load to an electrical load using a nominal ratio of the thermal to electrical power production for each generator. For these schemes, the generator needs to be connected to the supply side of a plant loop and serve components that use hot water on the demand side of the plant loop. The thermal load request is obtained from the plant data structure (structure location in code is PlantLoop%LoopSide%Branch%Comp%MyLoad). The distribution manager converts the thermal load, qthermal , to an electrical load using:
qElect=qthermalThermElectRatio
where,
ThermElectRatio is a nominal, constant, user-defined value for the ratio of thermal production to electrical production for a cogenerator. This ratio is used for supervisory control and dispatch of the electric power request to the generator; however, the cogenerator model may determine that actual performance varies from this nominal value at different times in the simulation when operating conditions differ from those used for the nominal ratio.
For all operating schemes except BASELOAD, a total electric load reduction target (or thermal load converted to electrical equivalent) is established for the load center based on the specific operating scheme. The load center then requests that its generators operate, one-by-one in the order specified, until the target is met or exceeded. Generators that are not scheduled as ‘available’ for the simulation time step are not called to operate. The requested power demand to be met by each generator is the smaller of the nominal ‘rated’ electric power output (as specified in the ElectricLoadCenter:Generators object) or the remaining total electric load reduction target for the load center. After each electric generator is requested to operate, the actual electric power delivered by the generator, which may be greater than or less than the requested amount due to inputs specified in the generator performance model (e.g., Generator:CombustionTurbine, Generator:MicroTurbine, etc.), is used to update the remaining total electric power target for the other generators associated with this load center.
Most of the operating schemes will sequentially load the available electric load centers and generators. EnergyPlus can accept multiple “ElectricLoadCenter:Distribution” objects with different operating schemes. Because of this, there are two levels of reporting, one for the whole building and a second for each load center. The whole-building results are managed with the internal meters for the entire model. The individual load-center results are summed for those generators connected to a particular load center. The total electricity purchased is reported both in power and energy units. This value is positive when the amount of energy is purchased from the utility. This value can be negative when the total electricity produced is greater than the facility electrical needs. The excess will either be available for storage or to sell back to the electric utility company.
The order of input objects (ElectricLoadCenter:Distribution) in the input file is significant and used to structure how generators are dispatched with the first load centers and generators managed before the later ones. Therefore, load centers listed earlier in the file effectively have a higher priority.
Electric load centers can have one of five different configurations. Load centers can get fairly complicated and include power conditioning and storage. Separate inverter models are used to condition DC power from photovoltaics into AC power for the building and utility. Load centers serving PV need to be specified with a direct current buss. The other generators may have inverters inside the devices but these are already embedded in the generator models. The load center can also manage electrical storage (ref. Electrical Storage below).
The most basic configuration is selected with the keyword “Alternating Current” for the Electrical Buss Type, shown in the following diagram.
The AlternatingCurrent load centers have AC generators with no storage and behave in the following way. All electric demand not met by the sum of the electrical power produced by the available generators will be met by purchased electricity. If a generator is needed in the simulation for a small load and the load is less than the generator’s minimum part load ratio, the generator will operate at the minimum part load ratio and the excess will either reduce demand or the excess energy will be exported back to the electric utility company. The purchased electrical demand limit is the user input for the demand limit above which the generators will try and meet the entire electrical load on the building. It is possible to prescribe a set of ElectricLoadCenter:Distribution objects with inconsistent or conflicting operating schemes, so users need to be careful.
A configuration with AC generators with on-site electrical storage is selected with the keyword “AlternatingCurrentWithStorage” and is shown in the following diagram.
The AlternatingCurrentWithStorage load centers attempt to augment the generator electricity production so that the power requests are met. Storage control logic is discussed below under Electrical Storage.
The basic configuration for photovoltaic generators is selected using the “DirectCurrentWithInverter” keyword and is shown in the following diagram.
The DirectCurrentWithInverter load centers collect DC power from various PV arrays, run the DC power through an inverter and produce AC power. The PV arrays produce DC power based on the availability of sunshine and do not respond to load requests made by the electric load center. The AC output from the inverter is what is recorded as electricity production.
If the PV-based load center is equipped with DC electrical storage that is connected before the inverter, then the buss type should be “DirectCurrentWithInverterDCStorage” and is shown in the following diagram.
The DirectCurrentWithInverterDCStorage load centers charge or draw DC power to meet the requested electrical load.
If the PV-based load center is equipped with AC electrical storage that is connected after the inverter, then the buss type should be “DirectCurrentWithInverterACStorage” and is shown in the following diagram.
The DirectCurrentWithInverterACStorage load centers charge or draw AC power to meet the requested electrical load.
Electric Load Center Generators[LINK]
The electric load center generators (object name: ElectricLoadCenter:Generators) provide a set of scheduled electric generators for electric power generation. Here is where the user lists what generators and PVs are available at any given time. For more details on the individual inputs required see the EnergyPlus Input Output Reference.
EnergyPlus includes three models for converting Direct Current (DC) electrical power into Alternating Current (AC) electrical power. The DC power into the inverter, PDC−in , is converted to AC power out, PAC−out , of the inverter using:
PAC−out=PDC−in⋅εinverter
The inverter efficiency is determined using one of the three models. For the“Simple” inveter model, efficiency is constant and input by the user. For the “Look Up Table” model, the efficiency is calculated using linear interpolation. For the “Function of Power” model, the efficiency is calculating using a single-variable curve object. For both the Look Up Table and Function of Power models, the power production is normalized by PDC−in .
The thermal losses are calculated from the difference between PDC−in and PAC−out .
Electrical Storage[LINK]
EnergyPlus includes two models for storing electrical energy: a simple model that is not intended to represent any specific type of storage technology and a battery model that represents the kinetic battery model originally developed by Manwell and McGowan.
The simple model might be called “constrained bucket with energy losses.” The “bucket” holds a quantity of Joules of electrical energy, refered to as the state of charge. There are losses and limits to storing and drawing power but otherwise the bucket just holds electricity. The user sets constraints on the rates of charging, Pstor−charge−max , and drawing, Pstor−draw−max . The user defines efficiencies for charging, εcharge , and drawing, εdraw . The user defines an initial state of charge and a maximum state of charge.
The storage control algorithms determine a value for the charging power, Pstor−charge , or the drawing power, Pstor−draw . The basic storage control algorithm is to compare the requested generator electricity loads to the current available supply and make up the difference with storage. If extra power is generated, then store it. If there is a shortage, then attempt to draw from storage to meet the remaining electricity request. The load center dispatchs a requested electric load for each generator, runs each generator, and then stores the actual power. This power dispatch can be a function of many different things depending on the operating scheme. The sum of the generator load requests, Pload−request , is then compared to the sum of the generator production, Pgen−supply
IF(Pload−request<Pgen−supply)Pstor−charge=Pgen−supply−Pload−request
IF(Pload−request>Pgen−supply)Pstor−draw=Pload−request−Pgen−supply
The limits Pstor−charge−max and Pstor−draw−max are applied.
If charging, the new state of charge, Qt+Δtstor , is determined using:
Qt+Δtstor=Qtstor+Pstor−charge⋅εcharge⋅Δt
If drawing, the new state of charge is:
Qt+Δtstor=Qtstor−Pstor−draw⋅Δtεdraw
Where Δt is the length of the system time step in seconds.
The storage device has an availability schedule. If it is not available then no power can be drawn or stored.
The gross electric power drawn and stored includes losses in the form of heat. These thermal losses are calculated from the user-specified efficiencies for charging and drawing and gross electric power stored and drawn. The thermal (heat) losses are included in a zone heat balance if the user specifies a thermal zone. A user-defined radiative split is used to divide thermal losses into radiation and convection portions. If no zone is specified, then the thermal losses are simply disregarded (e.g., rejected to outdoors and do not impact the zone air heat balance).
Electrical Storage - Kinetic Battery Model[LINK]
The Kinetic Battery Model (KiBaM) (object: ElectricLoadCenter:Storage:Battery) was originally developed by Manwell and McGowan (1993) for use in time series performance models of hybrid energy systems. The model is called kinetic because it is based on a chemical kinetics process to simulate the battery charging and discharging behavior. The model, with different improvements and modifications, has been incorporated into the software Hybrid2 and HOMER as the electrical storage module of hybrid and distributed power systems. In 2005, KiBaM was implemented as a stand-alone application in support of the European Union Benchmarking research project (Bindner et al. 2005).
The Kinetic Battery Model assumes that the battery charge is distributed over two tanks: an available-charge tank and a bound-charge tank. The tank for available charges can supply electrons directly to the load, whereas the tank for chemically bound charges can only supply electrons to the available-charge tank. At any time, the total charge q in the battery is the sum of the available charge ( q1 ) and bound charge ( q2 ). That is:
q=q1+q2
Based on the governing equations on the change of charge in both tanks (Manwell and McGowan 1993), the battery capacity can be related to a constant charge/discharge current ( I ) as the following equation:
qmax(I)=qmaxk⋅c⋅t1−e−kt+c(kt−1+e−kt)
where,
qmax(I) : Maximum capacity (Ah) at charge or discharge current I
qmax : Maximum capacity (Ah) at infinitesimal current
t : Charge or discharge time (hr), defined by t=qmax(I)I
k : Constant coefficient (hr-1)
c : Parameter indicating the ratio of available charge capacity to total capacity
Assuming that a constant current is used in any time step for charging and discharging, the available charge ( q1 ) and bound charge ( q2 ) at any time step are given by:
q1=q1,0e−kΔt+(q0kc−I)(1−e−kΔt)k−Ic(kΔt−1+e−kΔt)k
q2=q2,0e−kΔt+q0(1−c)(1−e−kΔt)−I(1−c)(kΔt−1+e−kΔt)k
where,
q1,0 : Available charge at the beginning of time step (Ah)
q2,0 : Bound charge at the beginning of time step (Ah)
q0 : Total charge at the beginning of time step (Ah), q0=q1,0+q2,0
Δt : Length of time step (hr)
KiBaM views the battery as a voltage source in series with an electric resistance (Figure 338). The internal resistance is assumed to be constant and the open circuit voltage varies with current and state of charge.
The battery’s open circuit voltage is modeled in the same form for charging and discharging, but with different coefficients. The open circuit voltage in charging ( Ec ) and in discharging ( Ed ) can be respectively expressed as:
Ec=E0,d+AcXc+CcXcDc−Xc
Ed=E0,c+AdXd+CdXdDd−Xd
where,
E0,c: Open circuit voltage for a fully charged battery
E0,d: Open circuit voltage for a fully discharged battery
Ac, Cc, Dc: Constant parameters for charging
Ad, Cd, Dd: Constant parameters for discharging
Xc, Xd: Normalized maximum capacity at a given charging or discharging current, calculated as:
X={q0/qmax(I)(charging)(qmax−q0)/qmax(I)(discharging)
It needs to be noted that the performance curve (Curve:RectangularHyperbola2) used in the model input covers the 2nd and the 3rd item of the open circuit voltage equation. Due to the reformatting of performance curve, the voltage function regression coefficients can map to the curve coefficients as follows:
C1=−C; C2=−D; C3=A
With open circuit voltage, the battery terminal voltage (V) can be calculated as:
V=E−IR
where, R is the battery internal resistance in Ohms; the current is positive for discharging and negative for charging.
Given desired power in/out of the battery, the desired charge or discharge current can be calculated from the basic power equation: P=VI . In this calculation, iteration is needed to ensure the electric current has converged and the battery operation satisfies all specified technical constraints such as maximum discharge current and charge rate limit.
KiBaM assumes that battery life is a primary function of charge/discharge cycles. One cycle is defined as the process of starting from a certain state of charge (SOC), the battery is discharged to a lower SOC and then recharged back to the starting SOC. It is regarded that the magnitude of cycle plays more important than the average of SOC during the cycle. This means that in terms of the impact on battery life, the cycle from 90% to 70% and then recharge back to 90% of SOC is equivalent to another cycle from 50% to 30% and then recharge back to 50% of SOC. Battery life in terms of the number of cycles is predicted as a function of the cycle range measured by the fractional depth of discharge. A double exponential equation is used to capture the damage to batteries due to cycling. The equation takes the following form where the coefficients need to be derived from battery test data via curve fitting.
CF=C1+C2eC3R+C4eC5R
where,
CF: Cycles to failure
C1 - C5: Regression coefficients
R: Cycle range in terms of fractional SOC
Following Hybrid2, the rainflow counting method (Downing and Socie 1982) is used to count battery cycles within a state of charge time series. Based on the number of cycles for each fractional SOC range, the battery damage is estimated as:
D=∑i=1(Ni1CF,i)
where,
D: Fractional battery damage. For example, a value of 0.5 at the end of simulation means that half of the battery life is used up after the length of the simulation period.
CF,i: Number of cycles to failure for the i-th cycle range
Ni: Total number of cycles over the simulation with the i-th cycle range
It needs to be noted that the temperature effects on battery performance and battery self-discharge are not supported in the current model.
Bindner H., Cronin T., Lundsager P., Manwell J.F., Abdulwahid U., and Baring-Gould I. 2005. Lifetime Modeling of Lead Acid Batteries. Riso National Laboratory, Roskilde, Denmark.
Downing S. D. and Socie D. F. 1982. Simple rainflow counting algorithms, International Journal of Fatigue, 1982.
Manwell J. F. and McGowan J. G. 1993. A lead acid battery storage model for hybrid energy systems, Solar Energy 50(5): 399- 405.
Transformers (object name: ElectricLoadCenter:Transformer) are an integral part of the electric distribution system. They have two broad applications closely related to building energy simulation. First, transformers are used to lower the voltage of electricity from utility primary circuits to customer secondary circuits, and in this case they are called distribution transformers. Second, transformers are used to output the surplus power from onsite generators to the electricity grid.
Distribution transformers reduce the voltage on utility distribution lines (34.5 kV or less) to a lower secondary voltage (600 V or less) suitable for customer equipment. Distribution transformers are usually categorized according to the medium used for cooling and insulation (liquid or air), the voltage class that they serve (low or medium), and the number of phases (single phase or three phase).
Liquid-immersed transformers rely on oil or other fire resistant liquid around the coils for cooling. In contrast, dry type transformers rely only on the natural convection of air for insulation and cooling. Medium-voltage transformers step from utility line voltage down to a lower secondary voltage, depending on the application. The secondary voltage from a medium-voltage transformer is usually at 277 V for single phase and 480 V for three phase. This secondary voltage can be directly used as 480 V three-phase power for large motors or as 277 V single-phase power for some fluorescent lighting. However, for most industrial and commercial facilities, low-voltage transformers are needed to reduce the above voltages further to 208/120 V. Common 120 V loads are wall plugs and incandescent lighting.
Most liquid-immersed transformers are owned by utilities and they are of the medium-voltage type. Virtually all dry type transformers are owned by commercial and industrial customers (Barnes et al. 1996). Of the dry type transformers, those of the medium-voltage type are generally special-order items while those of the low-voltage type are commodity items. The efficiency requirement of distribution transformers is covered by the NEMA (National Electrical Manufactures Association) Standard TP 1. ASHRAE 90.1-2010 will cite the NEMA Standard TP 1 to stipulate the efficiency requirement for the low-voltage dry type distribution transformers.
There are two main types of energy losses in transformers: no load loss and load loss. The no load loss comes primarily from the switching of the magnetic fields in the core material. Hence, it is also called the core loss. The no load (core) loss is roughly constant and exists continuously in the core material as long as the transformer is energized. The load loss comes from the electrical resistance in the windings when there is a load on the transformer. Hence, the load loss is also called the winding loss. The load (winding) loss is proportional to the load squared with a small temperature correction.
Given the no load loss (NL) and the load loss (LL) at rated load and conditions, the total energy losses in a transformer at time t is calculated as:
TL(t)=NL+LL(t)=NL+LL∗P(t)2∗fT(t)
where,
TL(t) Total energy loss at time t (W)
LL(t) Load loss at time t (W)
P(t) Per unit load at time t
fT(t) Temperature correction factor for the load loss at time t
The per unit load at time t is calculated as:
P(t)=Load(t)SB
where,
Load(t) Transformer load at time t (W)
SB Transformer nameplate rating (VA)
The temperature correction factor at time t is calculated as (NEMA 2002):
fT(t)=Ldc∗R(t)Rref+Leddy∗RrefR(t)
where,
Ldc Per unit load loss due to electrical resistance
Leddy Per unit load loss due to eddy currents
R(t) Winding electrical resistance at time t
Rref Winding electrical resistance at the full load reference conditions
The ratio of winding electrical resistance is calculated as:
R(t)Rref=F+Twinding(t)+Tamb(t)F+Twinding,ref+Tamb,ref
where,
F Thermal coefficient of resistance for the winding material (=225 for aluminum and 234.5 for copper)
Twinding,ref Winding temperature rise at the full load reference conditions (°C)
Twinding(t) Winding temperature rise at time t (°C)
Tamb,ref Ambient temperature at the reference condition (=20 °C)
Tamb(t) Ambient temperature at time t (°C)
The Ambient temperature Tamb(t) is equal to the zone temperature if a thermal zone is specified in the input; otherwise, it is assumed equal to 20 °C. The winding temperature rise at time t is calculated as (Barnes et al. 1997):
Twinding(t)=P(t)1.6∗Twinding.ref
Based on the derived total energy losses in a transformer, the transformer efficiency at time t can be calculated according to the following equation:
η(t)=Load(t)Load(t)+NL+LL(t)
The above procedure describes how to calculate the total transformer energy losses based on the no load loss and load loss at rated conditions. The transformer model also supports the case when the nominal transformer efficiency is given. In this case, the user needs to provide the nameplate efficiency and the corresponding per unit load, the maximum efficiency and the corresponding per unit load, and the reference conductor temperature at which the nameplate efficiency is measured. Given these information, both no load loss and load loss at rated conditions can be derived as below.
The nameplate efficiency can be expressed as:
ηnp=SB∗PnpSB∗Pnp+NL+LL∗P2np∗fT,np
where,
ηnp Nameplate efficiency
SB Nameplate rating (VA)
Pnp Per unit load at which the nameplate efficiency is measured
fT,np Applied temperature correction factor for the nameplate efficiency
Maximum efficiency generally occurs when the load loss is equal to the no-load loss. Because the no-load loss does not vary with the load on the transformer, the following relationship can be established:
NL=LLmax,η=LL∗P2max,η∗fT,max−η
where,
Pmax,η Per unit load at which the maximum efficiency is obtained
fT,max−η Applied temperature correction factor for the maximum efficiency
Transformers typically have close per unit loads for the nameplate efficiency and the maximum efficiency. Therefore, it is reasonable to assume that the applied temperature correction factors are equal at those two efficiencies. This implies that:
LLnpLLmax,η=LL∗Pnp2∗fT,npLL∗Pmax,η2∗fT,max−η=(PnpPmax,η)2
Rearranging Equation and combining it with Equation leads to:
LLnp=LLmax,η∗(PnpPmax,η)2=NL∗(PnpPmax,η)2
Combining Equations and , we can obtain the no load loss as:
NL=SB∗Pnp∗(1−ηnp)ηnp∗[1+(PrefPmax,η)2]
Substitute NL into Equation , we can calculate the load loss at rated conditions as:
LL=NLfT,max−η∗P2max,η=NLfT,np∗P2max,η
Since both no load and load losses at rated conditions are known, the total energy losses in a transformer at time t can then be calculated according to Equation
Barnes, PR., JW. Van Dyke, BW. McConnell, and S. Das. 1996. Determination Analysis of Energy Conservation Standards for Distribution Transformer, ORNL-6847. Oak Ridge National Laboratory, Oak Ridge, TN.
Barnes, PR., S. Das, BW. McConnell, and JW. Van Dyke. 1997. Supplement to the “Determination Analysis” (ORNL-6847) and Analysis of the NEMA Efficiency Standard for Distribution Transformer, ORNL-6925. Oak Ridge National Laboratory, Oak Ridge, TN.
NEMA. 2002. NEMA Standards Publication TP 1-2002: Guide for Determining Energy Efficiency for Distribution Transformers. National Electrical Manufactures Association, Rosslyn, VA.
Electric Load Center Distribution Manager[LINK]
Overview[LINK]
The electric load center distribution manager (object name: ElectricLoadCenter:Distribution) operates on-site generators specified in the simulation and reports the amount of generated and purchased electricity. Electric load centers really serve as a “load” to the generators and “supply” power to the rest of the building. The internal meters used by EnergyPlus for reporting do all of the demand tracking. For more details on the individual inputs required see the Input Output Reference document.
The electric load center manager sums all of the building and system electric loads and provides operating schemes for the generators. The electric load center objects are operated in the order they are defined in the input data file (IDF), and generators are dispatched sequentially in the order they are specified within each load center object. The electricity produced from photovoltaic arrays is handled somewhat separately and is always “used” to reduced the demand that the generators will try to meet for that time step. What is not provided by the on-site generation equipment, and electric storage units if specified, is met by (purchasing) off-site electricity.
The available operating schemes are “Baseload”, “DemandLimit”, “TrackElectrical,” “TrackSchedule,” “TrackMeter,” “FollowThermal” and “FollowThermalLimitElectrical.” These operating schemes affect how loads are dispatched to the generators, in effect telling the generators whether or not to run and requesting power levels.
The BASELOAD scheme operates the generators at their rated (requested) electric power output when the generator is scheduled ON (ref. ElectricLoadCenter:Generators in the Input Output Reference). The Baseload scheme requests all generators scheduled ON (available) to operate, even if the amount of electric power generated exceeds the total facility electric power demand.
The DEMAND LIMIT scheme limits the amount of purchased electricity from the utility to the amount specified in the input. The Demand Limit scheme tries to have the generators meet all of the demand above the purchased electric limit defined by the user.
The TRACK ELECTRICAL scheme tries to have the generators meet all of the electrical demand for the building.
The TRACK METER scheme tries to have the generators meet all the electrical demand from a meter chosen by the user rather than the usual meter for the entire facility. The meter can be a custom meter so that generators are tied to only certain loads in the building.
The TRACK SCHEDULE scheme tries to have the generators meet all of the electrical demand determined by a user-defined schedule.
The FOLLOW THERMAL and FOLLOW THERMAL LIMIT ELECTRICAL schemes run the generators to meet thermal demand. The thermal demand is determined from the plant modeling and depends on the flow requested by other components on the demand side of the plant loop, the loop temperatures, and the loop temperature setpoint. The electric load center distribution manager converts the thermal load to an electrical load using a nominal ratio of the thermal to electrical power production for each generator. For these schemes, the generator needs to be connected to the supply side of a plant loop and serve components that use hot water on the demand side of the plant loop. The thermal load request is obtained from the plant data structure (structure location in code is PlantLoop%LoopSide%Branch%Comp%MyLoad). The distribution manager converts the thermal load, qthermal , to an electrical load using:
qElect=qthermalThermElectRatio
where,
ThermElectRatio is a nominal, constant, user-defined value for the ratio of thermal production to electrical production for a cogenerator. This ratio is used for supervisory control and dispatch of the electric power request to the generator; however, the cogenerator model may determine that actual performance varies from this nominal value at different times in the simulation when operating conditions differ from those used for the nominal ratio.
For all operating schemes except BASELOAD, a total electric load reduction target (or thermal load converted to electrical equivalent) is established for the load center based on the specific operating scheme. The load center then requests that its generators operate, one-by-one in the order specified, until the target is met or exceeded. Generators that are not scheduled as ‘available’ for the simulation time step are not called to operate. The requested power demand to be met by each generator is the smaller of the nominal ‘rated’ electric power output (as specified in the ElectricLoadCenter:Generators object) or the remaining total electric load reduction target for the load center. After each electric generator is requested to operate, the actual electric power delivered by the generator, which may be greater than or less than the requested amount due to inputs specified in the generator performance model (e.g., Generator:CombustionTurbine, Generator:MicroTurbine, etc.), is used to update the remaining total electric power target for the other generators associated with this load center.
Most of the operating schemes will sequentially load the available electric load centers and generators. EnergyPlus can accept multiple “ElectricLoadCenter:Distribution” objects with different operating schemes. Because of this, there are two levels of reporting, one for the whole building and a second for each load center. The whole-building results are managed with the internal meters for the entire model. The individual load-center results are summed for those generators connected to a particular load center. The total electricity purchased is reported both in power and energy units. This value is positive when the amount of energy is purchased from the utility. This value can be negative when the total electricity produced is greater than the facility electrical needs. The excess will either be available for storage or to sell back to the electric utility company.
The order of input objects (ElectricLoadCenter:Distribution) in the input file is significant and used to structure how generators are dispatched with the first load centers and generators managed before the later ones. Therefore, load centers listed earlier in the file effectively have a higher priority.
Electric load centers can have one of five different configurations. Load centers can get fairly complicated and include power conditioning and storage. Separate inverter models are used to condition DC power from photovoltaics into AC power for the building and utility. Load centers serving PV need to be specified with a direct current buss. The other generators may have inverters inside the devices but these are already embedded in the generator models. The load center can also manage electrical storage (ref. Electrical Storage below).
The most basic configuration is selected with the keyword “Alternating Current” for the Electrical Buss Type, shown in the following diagram.
Basic Alternating Current Schematic
The AlternatingCurrent load centers have AC generators with no storage and behave in the following way. All electric demand not met by the sum of the electrical power produced by the available generators will be met by purchased electricity. If a generator is needed in the simulation for a small load and the load is less than the generator’s minimum part load ratio, the generator will operate at the minimum part load ratio and the excess will either reduce demand or the excess energy will be exported back to the electric utility company. The purchased electrical demand limit is the user input for the demand limit above which the generators will try and meet the entire electrical load on the building. It is possible to prescribe a set of ElectricLoadCenter:Distribution objects with inconsistent or conflicting operating schemes, so users need to be careful.
A configuration with AC generators with on-site electrical storage is selected with the keyword “AlternatingCurrentWithStorage” and is shown in the following diagram.
AC Generators with On-site Electrical Storage Schematic
The AlternatingCurrentWithStorage load centers attempt to augment the generator electricity production so that the power requests are met. Storage control logic is discussed below under Electrical Storage.
The basic configuration for photovoltaic generators is selected using the “DirectCurrentWithInverter” keyword and is shown in the following diagram.
Direct Current With Inverter Photovoltaic Generators Schematic
The DirectCurrentWithInverter load centers collect DC power from various PV arrays, run the DC power through an inverter and produce AC power. The PV arrays produce DC power based on the availability of sunshine and do not respond to load requests made by the electric load center. The AC output from the inverter is what is recorded as electricity production.
If the PV-based load center is equipped with DC electrical storage that is connected before the inverter, then the buss type should be “DirectCurrentWithInverterDCStorage” and is shown in the following diagram.
PV based Load Center with DC Electrical Storage Schematic
The DirectCurrentWithInverterDCStorage load centers charge or draw DC power to meet the requested electrical load.
If the PV-based load center is equipped with AC electrical storage that is connected after the inverter, then the buss type should be “DirectCurrentWithInverterACStorage” and is shown in the following diagram.
PV Based Load Center with AC Electrical Storage Schematic
The DirectCurrentWithInverterACStorage load centers charge or draw AC power to meet the requested electrical load.
Electric Load Center Generators[LINK]
The electric load center generators (object name: ElectricLoadCenter:Generators) provide a set of scheduled electric generators for electric power generation. Here is where the user lists what generators and PVs are available at any given time. For more details on the individual inputs required see the EnergyPlus Input Output Reference.
Inverters[LINK]
EnergyPlus includes three models for converting Direct Current (DC) electrical power into Alternating Current (AC) electrical power. The DC power into the inverter, PDC−in , is converted to AC power out, PAC−out , of the inverter using:
PAC−out=PDC−in⋅εinverter
The inverter efficiency is determined using one of the three models. For the“Simple” inveter model, efficiency is constant and input by the user. For the “Look Up Table” model, the efficiency is calculated using linear interpolation. For the “Function of Power” model, the efficiency is calculating using a single-variable curve object. For both the Look Up Table and Function of Power models, the power production is normalized by PDC−in .
The thermal losses are calculated from the difference between PDC−in and PAC−out .
Electrical Storage[LINK]
EnergyPlus includes two models for storing electrical energy: a simple model that is not intended to represent any specific type of storage technology and a battery model that represents the kinetic battery model originally developed by Manwell and McGowan.
The simple model might be called “constrained bucket with energy losses.” The “bucket” holds a quantity of Joules of electrical energy, refered to as the state of charge. There are losses and limits to storing and drawing power but otherwise the bucket just holds electricity. The user sets constraints on the rates of charging, Pstor−charge−max , and drawing, Pstor−draw−max . The user defines efficiencies for charging, εcharge , and drawing, εdraw . The user defines an initial state of charge and a maximum state of charge.
The storage control algorithms determine a value for the charging power, Pstor−charge , or the drawing power, Pstor−draw . The basic storage control algorithm is to compare the requested generator electricity loads to the current available supply and make up the difference with storage. If extra power is generated, then store it. If there is a shortage, then attempt to draw from storage to meet the remaining electricity request. The load center dispatchs a requested electric load for each generator, runs each generator, and then stores the actual power. This power dispatch can be a function of many different things depending on the operating scheme. The sum of the generator load requests, Pload−request , is then compared to the sum of the generator production, Pgen−supply
IF(Pload−request<Pgen−supply)Pstor−charge=Pgen−supply−Pload−request
IF(Pload−request>Pgen−supply)Pstor−draw=Pload−request−Pgen−supply
The limits Pstor−charge−max and Pstor−draw−max are applied.
If charging, the new state of charge, Qt+Δtstor , is determined using:
Qt+Δtstor=Qtstor+Pstor−charge⋅εcharge⋅Δt
If drawing, the new state of charge is:
Qt+Δtstor=Qtstor−Pstor−draw⋅Δtεdraw
Where Δt is the length of the system time step in seconds.
The storage device has an availability schedule. If it is not available then no power can be drawn or stored.
The gross electric power drawn and stored includes losses in the form of heat. These thermal losses are calculated from the user-specified efficiencies for charging and drawing and gross electric power stored and drawn. The thermal (heat) losses are included in a zone heat balance if the user specifies a thermal zone. A user-defined radiative split is used to divide thermal losses into radiation and convection portions. If no zone is specified, then the thermal losses are simply disregarded (e.g., rejected to outdoors and do not impact the zone air heat balance).
Electrical Storage - Kinetic Battery Model[LINK]
The Kinetic Battery Model (KiBaM) (object: ElectricLoadCenter:Storage:Battery) was originally developed by Manwell and McGowan (1993) for use in time series performance models of hybrid energy systems. The model is called kinetic because it is based on a chemical kinetics process to simulate the battery charging and discharging behavior. The model, with different improvements and modifications, has been incorporated into the software Hybrid2 and HOMER as the electrical storage module of hybrid and distributed power systems. In 2005, KiBaM was implemented as a stand-alone application in support of the European Union Benchmarking research project (Bindner et al. 2005).
The Kinetic Battery Model assumes that the battery charge is distributed over two tanks: an available-charge tank and a bound-charge tank. The tank for available charges can supply electrons directly to the load, whereas the tank for chemically bound charges can only supply electrons to the available-charge tank. At any time, the total charge q in the battery is the sum of the available charge ( q1 ) and bound charge ( q2 ). That is:
q=q1+q2
Based on the governing equations on the change of charge in both tanks (Manwell and McGowan 1993), the battery capacity can be related to a constant charge/discharge current ( I ) as the following equation:
qmax(I)=qmaxk⋅c⋅t1−e−kt+c(kt−1+e−kt)
where,
qmax(I) : Maximum capacity (Ah) at charge or discharge current I
qmax : Maximum capacity (Ah) at infinitesimal current
t : Charge or discharge time (hr), defined by t=qmax(I)I
k : Constant coefficient (hr-1)
c : Parameter indicating the ratio of available charge capacity to total capacity
Assuming that a constant current is used in any time step for charging and discharging, the available charge ( q1 ) and bound charge ( q2 ) at any time step are given by:
q1=q1,0e−kΔt+(q0kc−I)(1−e−kΔt)k−Ic(kΔt−1+e−kΔt)k
q2=q2,0e−kΔt+q0(1−c)(1−e−kΔt)−I(1−c)(kΔt−1+e−kΔt)k
where,
q1,0 : Available charge at the beginning of time step (Ah)
q2,0 : Bound charge at the beginning of time step (Ah)
q0 : Total charge at the beginning of time step (Ah), q0=q1,0+q2,0
Δt : Length of time step (hr)
KiBaM views the battery as a voltage source in series with an electric resistance (Figure 338). The internal resistance is assumed to be constant and the open circuit voltage varies with current and state of charge.
Electrical equivalent model for KiBaM
The battery’s open circuit voltage is modeled in the same form for charging and discharging, but with different coefficients. The open circuit voltage in charging ( Ec ) and in discharging ( Ed ) can be respectively expressed as:
Ec=E0,d+AcXc+CcXcDc−Xc
Ed=E0,c+AdXd+CdXdDd−Xd
where,
E0,c: Open circuit voltage for a fully charged battery
E0,d: Open circuit voltage for a fully discharged battery
Ac, Cc, Dc: Constant parameters for charging
Ad, Cd, Dd: Constant parameters for discharging
Xc, Xd: Normalized maximum capacity at a given charging or discharging current, calculated as:
X={q0/qmax(I)(charging)(qmax−q0)/qmax(I)(discharging)
It needs to be noted that the performance curve (Curve:RectangularHyperbola2) used in the model input covers the 2nd and the 3rd item of the open circuit voltage equation. Due to the reformatting of performance curve, the voltage function regression coefficients can map to the curve coefficients as follows:
C1=−C; C2=−D; C3=A
With open circuit voltage, the battery terminal voltage (V) can be calculated as:
V=E−IR
where, R is the battery internal resistance in Ohms; the current is positive for discharging and negative for charging.
Given desired power in/out of the battery, the desired charge or discharge current can be calculated from the basic power equation: P=VI . In this calculation, iteration is needed to ensure the electric current has converged and the battery operation satisfies all specified technical constraints such as maximum discharge current and charge rate limit.
KiBaM assumes that battery life is a primary function of charge/discharge cycles. One cycle is defined as the process of starting from a certain state of charge (SOC), the battery is discharged to a lower SOC and then recharged back to the starting SOC. It is regarded that the magnitude of cycle plays more important than the average of SOC during the cycle. This means that in terms of the impact on battery life, the cycle from 90% to 70% and then recharge back to 90% of SOC is equivalent to another cycle from 50% to 30% and then recharge back to 50% of SOC. Battery life in terms of the number of cycles is predicted as a function of the cycle range measured by the fractional depth of discharge. A double exponential equation is used to capture the damage to batteries due to cycling. The equation takes the following form where the coefficients need to be derived from battery test data via curve fitting.
CF=C1+C2eC3R+C4eC5R
where,
CF: Cycles to failure
C1 - C5: Regression coefficients
R: Cycle range in terms of fractional SOC
Following Hybrid2, the rainflow counting method (Downing and Socie 1982) is used to count battery cycles within a state of charge time series. Based on the number of cycles for each fractional SOC range, the battery damage is estimated as:
D=∑i=1(Ni1CF,i)
where,
D: Fractional battery damage. For example, a value of 0.5 at the end of simulation means that half of the battery life is used up after the length of the simulation period.
CF,i: Number of cycles to failure for the i-th cycle range
Ni: Total number of cycles over the simulation with the i-th cycle range
It needs to be noted that the temperature effects on battery performance and battery self-discharge are not supported in the current model.
References[LINK]
Bindner H., Cronin T., Lundsager P., Manwell J.F., Abdulwahid U., and Baring-Gould I. 2005. Lifetime Modeling of Lead Acid Batteries. Riso National Laboratory, Roskilde, Denmark.
Downing S. D. and Socie D. F. 1982. Simple rainflow counting algorithms, International Journal of Fatigue, 1982.
Manwell J. F. and McGowan J. G. 1993. A lead acid battery storage model for hybrid energy systems, Solar Energy 50(5): 399- 405.
Electric Load Center Transformers[LINK]
Transformers (object name: ElectricLoadCenter:Transformer) are an integral part of the electric distribution system. They have two broad applications closely related to building energy simulation. First, transformers are used to lower the voltage of electricity from utility primary circuits to customer secondary circuits, and in this case they are called distribution transformers. Second, transformers are used to output the surplus power from onsite generators to the electricity grid.
Distribution transformers reduce the voltage on utility distribution lines (34.5 kV or less) to a lower secondary voltage (600 V or less) suitable for customer equipment. Distribution transformers are usually categorized according to the medium used for cooling and insulation (liquid or air), the voltage class that they serve (low or medium), and the number of phases (single phase or three phase).
Liquid-immersed transformers rely on oil or other fire resistant liquid around the coils for cooling. In contrast, dry type transformers rely only on the natural convection of air for insulation and cooling. Medium-voltage transformers step from utility line voltage down to a lower secondary voltage, depending on the application. The secondary voltage from a medium-voltage transformer is usually at 277 V for single phase and 480 V for three phase. This secondary voltage can be directly used as 480 V three-phase power for large motors or as 277 V single-phase power for some fluorescent lighting. However, for most industrial and commercial facilities, low-voltage transformers are needed to reduce the above voltages further to 208/120 V. Common 120 V loads are wall plugs and incandescent lighting.
Most liquid-immersed transformers are owned by utilities and they are of the medium-voltage type. Virtually all dry type transformers are owned by commercial and industrial customers (Barnes et al. 1996). Of the dry type transformers, those of the medium-voltage type are generally special-order items while those of the low-voltage type are commodity items. The efficiency requirement of distribution transformers is covered by the NEMA (National Electrical Manufactures Association) Standard TP 1. ASHRAE 90.1-2010 will cite the NEMA Standard TP 1 to stipulate the efficiency requirement for the low-voltage dry type distribution transformers.
There are two main types of energy losses in transformers: no load loss and load loss. The no load loss comes primarily from the switching of the magnetic fields in the core material. Hence, it is also called the core loss. The no load (core) loss is roughly constant and exists continuously in the core material as long as the transformer is energized. The load loss comes from the electrical resistance in the windings when there is a load on the transformer. Hence, the load loss is also called the winding loss. The load (winding) loss is proportional to the load squared with a small temperature correction.
Given the no load loss (NL) and the load loss (LL) at rated load and conditions, the total energy losses in a transformer at time t is calculated as:
TL(t)=NL+LL(t)=NL+LL∗P(t)2∗fT(t)
where,
TL(t) Total energy loss at time t (W)
LL(t) Load loss at time t (W)
P(t) Per unit load at time t
fT(t) Temperature correction factor for the load loss at time t
The per unit load at time t is calculated as:
P(t)=Load(t)SB
where,
Load(t) Transformer load at time t (W)
SB Transformer nameplate rating (VA)
The temperature correction factor at time t is calculated as (NEMA 2002):
fT(t)=Ldc∗R(t)Rref+Leddy∗RrefR(t)
where,
Ldc Per unit load loss due to electrical resistance
Leddy Per unit load loss due to eddy currents
R(t) Winding electrical resistance at time t
Rref Winding electrical resistance at the full load reference conditions
The ratio of winding electrical resistance is calculated as:
R(t)Rref=F+Twinding(t)+Tamb(t)F+Twinding,ref+Tamb,ref
where,
F Thermal coefficient of resistance for the winding material (=225 for aluminum and 234.5 for copper)
Twinding,ref Winding temperature rise at the full load reference conditions (°C)
Twinding(t) Winding temperature rise at time t (°C)
Tamb,ref Ambient temperature at the reference condition (=20 °C)
Tamb(t) Ambient temperature at time t (°C)
The Ambient temperature Tamb(t) is equal to the zone temperature if a thermal zone is specified in the input; otherwise, it is assumed equal to 20 °C. The winding temperature rise at time t is calculated as (Barnes et al. 1997):
Twinding(t)=P(t)1.6∗Twinding.ref
Based on the derived total energy losses in a transformer, the transformer efficiency at time t can be calculated according to the following equation:
η(t)=Load(t)Load(t)+NL+LL(t)
The above procedure describes how to calculate the total transformer energy losses based on the no load loss and load loss at rated conditions. The transformer model also supports the case when the nominal transformer efficiency is given. In this case, the user needs to provide the nameplate efficiency and the corresponding per unit load, the maximum efficiency and the corresponding per unit load, and the reference conductor temperature at which the nameplate efficiency is measured. Given these information, both no load loss and load loss at rated conditions can be derived as below.
The nameplate efficiency can be expressed as:
ηnp=SB∗PnpSB∗Pnp+NL+LL∗P2np∗fT,np
where,
ηnp Nameplate efficiency
SB Nameplate rating (VA)
Pnp Per unit load at which the nameplate efficiency is measured
fT,np Applied temperature correction factor for the nameplate efficiency
Maximum efficiency generally occurs when the load loss is equal to the no-load loss. Because the no-load loss does not vary with the load on the transformer, the following relationship can be established:
NL=LLmax,η=LL∗P2max,η∗fT,max−η
where,
Pmax,η Per unit load at which the maximum efficiency is obtained
fT,max−η Applied temperature correction factor for the maximum efficiency
Transformers typically have close per unit loads for the nameplate efficiency and the maximum efficiency. Therefore, it is reasonable to assume that the applied temperature correction factors are equal at those two efficiencies. This implies that:
LLnpLLmax,η=LL∗Pnp2∗fT,npLL∗Pmax,η2∗fT,max−η=(PnpPmax,η)2
Rearranging Equation and combining it with Equation leads to:
LLnp=LLmax,η∗(PnpPmax,η)2=NL∗(PnpPmax,η)2
Combining Equations and , we can obtain the no load loss as:
NL=SB∗Pnp∗(1−ηnp)ηnp∗[1+(PrefPmax,η)2]
Substitute NL into Equation , we can calculate the load loss at rated conditions as:
LL=NLfT,max−η∗P2max,η=NLfT,np∗P2max,η
Since both no load and load losses at rated conditions are known, the total energy losses in a transformer at time t can then be calculated according to Equation
References:[LINK]
Barnes, PR., JW. Van Dyke, BW. McConnell, and S. Das. 1996. Determination Analysis of Energy Conservation Standards for Distribution Transformer, ORNL-6847. Oak Ridge National Laboratory, Oak Ridge, TN.
Barnes, PR., S. Das, BW. McConnell, and JW. Van Dyke. 1997. Supplement to the “Determination Analysis” (ORNL-6847) and Analysis of the NEMA Efficiency Standard for Distribution Transformer, ORNL-6925. Oak Ridge National Laboratory, Oak Ridge, TN.
NEMA. 2002. NEMA Standards Publication TP 1-2002: Guide for Determining Energy Efficiency for Distribution Transformers. National Electrical Manufactures Association, Rosslyn, VA.
Documentation content copyright © 1996-2015 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.