From GridLAB-D Wiki
Jump to: navigation, search

WECC/LMTF San Diego Meeting (4 March 2009) Comments

  1. High priority
  2. Low priority
    • Temperature units in the load composition table still show the TMY units in the labels. This is confusing. Try hiding everything in [].
    • Latitude unit is N. It should be degrees (N is positive).
    • Check units of solar flux Wh/m^2 vs. Btu/sf.h. Probably ok, because of hourly conversion, but check anyway.
    • Summer Phoenix 65% Motor D is right per Richard Bravo.
    • How to use the ZIP fraction?
    • Confirm that composition can be scaled by energy by using the “scalar” cell on “composition” tab.
    • AMI data are available. Explain how that might be incorporated.
    • Rule of association: engineering judgment. Current data are based on experience. Need to be reviewed by MVWG. Emphasize importance of reviewing these in instructions.
    • Commercial buildings: sf is set to the same sf for one category.
    • Remarks on commercial tab need to be revised. Some are not limited.
    • PF is for ZIP only on the feeder tab. Change the label to be “ZIP PF”.
    • Design conditions are used as extreme conditions on sensitivity tab. Summer – 1 degree down; winter – 1 degree up.
    • The table layout is important as formulas rely on the layout. Changing the layout is dangerous.
    • Need to distribute to MVWG for review and test.
    • Step-by-step instructions to get the spreadsheet from the internet. Website is SourceForge project GridLab-D, search for "load composition".


A/C unit testing results
Current P (W) Q(VAR) Stall threshold
80°F 100°F 115°F 80°F 100°F 115°F 80°F 100°F 115°F 80°F 100°F 115°F
EPRI_Coleman_2.5_R410_S 9 11 13 1.85 2.46 2.87 0.65 0.68 0.71 0.55 0.61 0.66
EPRI_Lennox_3_R22_S 10 12 13 2.22 2.64 2.95 0.76 0.79 0.82 0.55 0.57 0.60
SCE_Carrier_3_R22_R 14 16 17 3.10 3.50 3.80 1.30 1.30 1.30 0.64 0.65 0.68
SCE_Rheem_3_R22_S 14 17 20 3.20 3.90 4.70 1.20 1.25 1.35 0.58 0.65 0.70
SCE_Carrier_3_R22_S 12 15 17 2.75 3.35 4.05 0.85 0.85 1.00 0.55 0.60 0.65
SCE_Carrier_3_R410_S 12 15 18 2.80 3.50 4.30 0.65 0.70 0.85 0.60 0.65 0.73
EPRI_Goodman_3.5_R22_S 14 16 19 2.90 3.53 4.15 1.30 1.36 1.44 0.55 0.65 0.66
EPRI_Carrier_3.5_R410_S 17 20 21 3.66 4.41 4.70 1.59 1.66 1.66 0.55 0.55 0.64
EPRI_Coleman_3.5_R410_S 18 20 23 3.95 4.38 5.07 1.78 1.84 1.93 0.60 0.62 0.66
EPRI_Lennox_4_R22_S 12 16 18 2.86 3.60 4.24 0.29 0.40 0.53 0.60 0.65
EPRI_Trane_4_R22_S 13 15 17 2.89 3.56 3.73 0.54 0.73 0.79 0.57 0.63 0.64
EPRI_Armstrong_4_R22_R 18 19 19 4.20 4.40 4.49 0.55 0.58 0.66 0.62 0.65 0.66
EPRI_Carrier_4_R410_S 19 21 24 4.07 4.58 5.10 2.12 2.09 2.12 0.56 0.60 0.64
SCE_Goodman_4_R22_R 17 19 21 3.80 4.25 4.65 1.70 1.75 1.75 0.60 0.65 0.68
SCE_Day & Night_4_R22_R 21 23 25 4.40 5.00 5.40 2.15 2.25 2.30 0.60 0.70 0.70
SCE_Carrier_4_R22_S 16 19 22 3.60 4.40 5.10 1.05 1.10 1.15 0.55 0.60 0.65
SCE_Coleman_4_R410_S 17 22 26 4.20 5.00 6.05 1.10 1.25 1.40 0.55 0.60 0.65
SCE_Rheem_4_R22_S 19 22 25 3.90 4.70 5.45 2.50 2.55 2.60 0.58 0.63 0.68
SCE_Rheem OC_4_R22_S 22 26 32 4.50 5.60 6.85 2.60 2.65 2.75 0.60 0.70 0.75
SCE_Carrier_5_R410_S 21 26 31 4.75 6.00 7.10 1.60 1.75 1.95 0.63 0.68
EPRI_Goodman_5_R22_S 16 20 25 3.61 4.62 5.61 1.17 1.27 1.42 0.57 0.60 0.60
EPRI_Carrier_5_R410_R 8 9 25 1.80 2.07 5.70 -0.04 0.00 1.14 0.55 0.57 0.60

Comments from Anatoliy

Note
This comment was receive by email and posted for reference. It refers to an Excel spreadsheet that cannot be uploaded. Contact David Chassin to obtain a copy. --Dchassin 17:43, 23 March 2009 (UTC)

Please consider a practical method to accommodate user knowledge about load class proportions for a given feeder/substation/area. This method is compatible with your “from bottom up” approach. The needed additions are illustrated in the two attached spreadsheets which supposed to be added to your Excel files.

The idea is in considering typical customer numbers for load categories in the “composition” sheet (column B) as variables. To find a number of customers (or scaling multipliers to initially assumed numbers), we need to calculate monthly kWh’s for the initial mixture of typical customers and correct the number of customers to satisfy the proportions given by the user for the entire feeder/substation/area. The attached Sheet 5 includes two load shape tables I obtained using your tool for the same set of conditions and differed only by the day of a week (weekday and weekend). Sheet 5 calculates kWh’s for 24 hours of the given day by summarizing every hour kWs. Doing that for weekdays and weekends of the given month and multiplying totals by say 22 and 9 and dividing the sum by 31, we get the needed average monthly energy for the given mixture of typical customers.

Sheet 5 includes also fractions of average daily energy obtained for 6 transmission buses for the same month as the load shapes (in general for unlimited number of buses).

Multipliers for adjusting customer numbers are calculated to satisfy the following:

<math>E_{ti}/E_{t\Sigma} \cdot M_i = E_{bi}/E_{bS}\,\!</math>

where <math>E_{tS}</math> - daily average energy for original mixture of typical customers <math>E_{tS}</math> - daily average energy for class <math>i</math> of typical customers <math>E_{bi}</math> - billing data based fraction of daily average energy for class <math>i</math> bus load <math>E_{bS} =1</math>, <math>M_i</math> - multiplier to reduce/increase number of all customer categories inside class <math>i</math> loads. <math>M_i</math> values are calculated in cells L37:O42 using formula:

<math>M_i = E_{bi} \cdot E_{tS} / E_{ti} \,\!</math>

Sheet 5 calculates multipliers for customer numbers assuming that the original set of numbers correctly reflects customer mixture inside load classes (we would not need to make this assumption having load shapes for each customer category and corresponding billing data based fractions for those categories).

Multipliers from Sheet 5 are used in Sheet 4 to scale numbers of customers. Besides this scaling, Sheet 4 is a duplicate of the “composition” spreadsheet.

Cell P28 in Sheet 5 is used for selection of a bus to process. Having say 3 in P28 puts the corresponding multipliers (with “3” in column K) to buffer 1 for using in Sheet 4. The results from Sheet 4 are delivered to the output model composition table via buffer 2. All the results can be obtained in one Excel execution with automatic increase of numbers in P28 (not done).

Another comment from Anatoliy

Note
This comment was received by email and posted for reference. --Dchassin 17:43, 23 March 2009 (UTC)

One my comment in the previous e-mail was about defining component P fractions relatively to total real power. The formulas for fractions should also take into account the exact meaning of CMPLDW inputs for zip load components. The coefficients P1c, P2c, Q1c and Q2c are defined in CMPLDW as fractions of only ZIP real and reactive powers (not total powers). Btw, do we really want to use LCM for deriving zip real and reactive load compositions? How sure are we LCM provides this composition for nominal voltages? I would fill better using default values for P1c, P2c, Q1c, Q2c and PFs.

Do we consider electronic load 100% real? If not, where is PFel needed for CMPLDW?

It looks like I figured out how to count the bus-by-bus load composition for unlimited number of loads in one Excel execution. It takes 20-30 seconds to get adjusted load model compositions for 1,400 loads.

I also made some corrections in additional sheets to count the situations when assumed number of customers in one or more load classes is zero but it is not zero in the known composition. If such situation occurs, there is no multiplier to the customer numbers which could provide a desirable proportion. To exclude such situations, I corrected formulas for customer numbers by always adding 0.1 customer (could be 0.01 or even less) for each load category in the "composition" sheet. This provides non-zero load shapes and energies for all categories and always derives the finite value scaling multipliers. Unfortunately, I need to make the described formula corrections in the main "composition" sheet. I would prefer to do this in a copy if this sheet, but coping does not work for your "Update Loadshapes" button. Is there simple way to correct your macros and get this button in a copy sheet?

Comment from Richard Bravo

Note
This comment was received by email and posted for reference. --Dchassin 17:46, 23 March 2009 (UTC)

The fraction of Motor D is somewhat similar for all regions. We believe that the current spreadsheet assumes that every house in each area has an air-conditioner. This assumption is not true, particularly for coastal and northern areas. I agree that new or remodeled homes will have A/C but older might not.

For Portland, while newer residential construction is very likely to have an air-conditioner, most of older construction does not.

  1. There is a value of adding a parameter “residential AC penetration” to the calculations. I think the same thing happens here in SCE.
  2. How can we get information on residential AC penetration by region? We might have to do this by area instead by region. In our case we have low penetration of A/C near the coast (~20%) where the temperature does not increase above 85F. We have high penetration of A/C inland (40% to 65%) where the summer temperatures goes from 90 to 110F. The case of Portland might be on about 20% but if you go east this numbers might rise drastically.
A less important question
Is there a value of adding diversity in the thermostat settings? Not everyone will turn their AC on during a normal (80-85F) summer day. What about a coincidental factor? We know that some air conditioners are "shut down" but others will be cycling (momentarily off) by thermostats.
Fixed in V1.6.6
See #Air-conditioner penetration factor

SEER temperature correction

[--Dchassin 20:46, 11 May 2009 (UTC)]

The approach to including outdoor air temperature effects in the cooling load calculation is as follows.

The following inputs are required

  • <math>SEER</math> is the (Seasonal Energy Efficiency Rating) value for the average equipment in use in Btu/W.h.
  • <math>T_{design}</math> is the cooling design temperature used for sizing the equipment in °F.
  • <math>T_{set}</math> is the thermostat cooling setpoint in °F.
  • <math>T_{out}</math> is the outdoor temperature in °F.

The maximum EER (energy efficiency ratio in Btu/W.h) is typically used to calculate the SEER value using the assumption <math>SEER = 0.875 EER_{max}</math>. The EER is determined as follows:

<math>

EER_{max} = EER_{coeff} \times 3.412 \frac{T_{set}+459.7}{{T_{design}}-T_{set}} \qquad (1)

</math>

where the <math>EER_{coeff}</math> is typically around <math>0.1</math>. This is the relatively invariant ratio of the Carnot efficiency to the actual efficiency and can be used to compute <math>EER_{coeff}</math>:

<math>

EER_{coeff} = \frac{SEER}{0.875} \frac{T_{design}-T_{set}}{3.412 (T_{set}+459.7) } \qquad (2)

</math>

The theoretical EER, <math>EER_{theoretical}</math> can thus be computed for any outdoor temperature

<math>

EER_{theoretical} = EER_{coeff} \times 3.412 \frac{T_{set}+459.7}{{T_{out}}-T_{set}} \qquad (3)

</math>

The test results obtained do not show this large variation in EER and <math>EER_{theoretical}</math> must be adjusted using the following formula to obtain the <math>EER_{model}</math>

<math>

EER_{model} = EER_{theoretical} \left ( 1 - E_{factor}(T_{design}-T_{out}) \right ) \qquad (4)

</math>

where

  • <math>E_{factor}</math> is estimated by minimizing the weighted RMS error between the observed power used in the tests and the model's prediction at various cooling design temperatures, using the weights for the equipment sizes, as shown in Table 1. The results for the model calibaration are shown in Table 2.
Note
If the outdoor temperature is below 83°F then the SEER values are for 83 °F are used.
Table 1. Residential air-conditioner size weighting factors
Size (tons) Weight
2.5 10%
3.0 20%
3.5 30%
4.0 20%
5.0 20%
Table 2: Model calibration results
Tdesign (°F) Efactor
90 0.05914
95 0.04613
100 0.03782
105 0.03197
110 0.02764
115 0.02433
Figure 1: Diversified Single-family AC load curves

A quadratic fit of this data yields the following relation

<math>

E_{factor} = 4.6\times10^{-5} T_{design}^2 - 1.08\times10^{-2} T_{design} + 0.656 \qquad (5)

</math>

with <math>R^2 = 0.9969</math>.

Figure 1 shows the diversified AC load curves (in kWh/h) as a function of outdoor temperature (in °F) for various cooling design conditions.

Caveat
This approach assumes that the thermostat setpoint is maintained. If the duty cycle is 100%, that assumption is not valid and the efficiency may be better than that predicted based on the change in outdoor temperature alone.

Air-conditioner penetration factor

[--Dchassin 19:21, 12 May 2009 (UTC)]

The air-conditioner penetration factor specifies what fraction of homes have operating air-conditioning units. This number is based on the cooling design temperature established for the model, as follows:

<math>

r_{AC} = \begin{cases}

 0.1 : T_{design} \le 80 \\
 0.1 + 0.9\frac{T_{design}-80}{20} : 80 < T_{design} < 100 \\
 1.0 : T_{design} \ge 100 \\

\end{cases}

</math>

The air-conditioning penetration factor is used to discount the diversified cooling load as follows:

<math>

D_{cooling} = r_{AC}\ D_{cooling}\ P_{cooling}\ A/1000

</math>

where

  • <math>D_{cooling}</math> is the cooling duty cycle
  • <math>P_{cooling}</math> is the cooling power density
  • <math>A</math> is the floor area
Fixed
Done in version 1.6.6 [--Dchassin 15:59, 12 June 2009 (UTC)]

Retail process load density is too high

[--Dchassin 22:08, 12 May 2009 (UTC)]

This problem seems to stem from a very low square footage for building with process loads in the CEUS data. As a result, the normalization process give it too much weight. The solution is to make the square footage representative of reality, rather than the sampling of the CEUS data.

Fixed
At this time, I recommend using 50,000 sf instead of 1,000 square, but this number is just a guess. This will need to be researched more carefully before we commit to any value. Done in version 1.6.6 [--Dchassin 15:59, 12 June 2009 (UTC)].

Comments from Dmitry

These comments were received from Dmitry Kosterev (BPA) on Fri Jun 5, 2006.

Rules of Association

Please see [below] suggested "rules of association". Your guess may be as good (or better) than mine.

Single-family Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Cooking 0.25 0.75
Hotwater 1.00
Lighting 0.20 0.80
Plugs 1.00
Washing 0.25 0.75
Heating 0.25 0.75
Cooling 1.00
Refrigeration 1.00


Multi-family Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Cooking 0.25 0.75
Hotwater 1.00
Lighting 0.20 0.80
Plugs 1.00
Washing 0.25 0.75
Heating 0.25 0.75
Cooling 1.00
Refrigeration 1.00


Small office Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Heating 0.50 0.40 0.10 0.015
Cooling 0.75 0.25
Vent 0.30 0.70
WaterHeat 1.00 0.15
Cooking 0.20 0.20 0.60
Refrig 0.40 0.60
ExtLight 1.00 -1.20 0.20
IntLight 1.00 -1.20 0.20
OfficeEquip 1.00
Misc 1.00
Process 0.50 0.50
Motors 0.50 0.50
AirComp 1.00


Large office Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Heating 0.75 0.15 0.10 0.015
Cooling 1.0
Vent 0.30 0.70
WaterHeat 1.00 0.15
Cooking 0.20 0.20 0.60
Refrig 0.40 0.60
ExtLight 1.00 -1.20 0.20
IntLight 1.00 -1.20 0.20
OfficeEquip 1.00
Misc 1.00
Process 0.50 0.50
Motors 0.50 0.50
AirComp 1.00


Retail Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Heating 0.75 0.15 0.10 0.015
Cooling 1.0
Vent 0.30 0.70
WaterHeat 1.00 0.15
Cooking 0.20 0.20 0.60
Refrig 0.40 0.60
ExtLight 1.00 -1.20 0.20
IntLight 1.00 -1.20 0.20
OfficeEquip 1.00
Misc 1.00
Process 0.50 0.50
Motors 0.50 0.50
AirComp 1.00


Lodging Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Heating 0.25 0.25 0.5 0.01
Cooling 0.40 0.60
Vent 0.30 0.70
WaterHeat 1.00 0.15
Cooking 0.20 0.20 0.60
Refrig 0.40 0.60
ExtLight 1.00 -1.20 0.20
IntLight 1.00 -1.20 0.20
OfficeEquip 1.00
Misc 1.00
Process 0.50 0.50
Motors 0.50 0.50
AirComp 1.00


All other commercial Rules of Association
Electronic Motor-A Motor-B Motor-C Motor-D ZIP Ip ZIP Iq ZIP Pp (P) ZIP Pq (Q) ZIP Zp (G) ZIP Zq (B)
Heating 0.75 0.25 0.10 0.015
Cooling 0.40 0.60
Vent 0.30 0.70
WaterHeat 1.00 0.15
Cooking 0.20 0.20 0.60
Refrig 0.40 0.60
ExtLight 1.00 -1.20 0.20
IntLight 1.00 -1.20 0.20
OfficeEquip 1.00
Misc 1.00
Process 0.50 0.50
Motors 0.50 0.50
AirComp 1.00

AC Load Calculation

There seem to be an issue with AC load calculation. In load shape calculation, the kW amount of Motor D (corresponding to AC load) is nearly constant for the day.

Fixed
The problem is caused by how the temperature is determined based on the study type. Peak load studies used the peak temperature regardless of the time of day. This was fixed in version 1.6.7 by having the temperature (and all other conditions) be taken for the time of day given. [--Dchassin 15:57, 12 June 2009 (UTC)]

Enduse Electrification by City

Is it possible to have a support file (or a sheet) with default entries for each city:

End-use electrical defaults
Enduse electrification Fraction of customers that use electric…
Resistive heating fraction 80%
Heatpump fraction 20%
Electric hotwater fraction 80%
Electric cooking fraction 80%
Electric drying fraction 50%
Air-conditioning fraction 100%

The default data is to be used in "Feeder" tab calculations.

The reason is that Phoenix may have 100% houses with air-conditioners, while Portland say 50% and Seattle 25%.

Done
This was fixed in version 1.6.6 [--Dchassin 15:57, 12 June 2009 (UTC)]

Rob Pratt comments

Single family

  1. Building height default to 8 feet
  2. Roof insulation should be composition less than R60, suggest R40
  3. Window should be about R2.5 (double w/low-E)
  4. Window/wall ratio is high (probably less than 10%)
  5. Suggest adding door wall ratio with door Rvalue (2 doors minimum)
  6. ACH should not be below 0.5 and no usually above 1.0 on average
  7. heating and cooling capacity needs to include internal gains
  8. check solarexposure calcs for fully diversified orientation
  9. building UA should include doors
  10. need to modernize ELCAP gains
  11. add UA heat gain to latent load numerator
  12. in rules of association
    1. what about CFLs
    2. why ZIP in washing? isn't washing 100% motor-D?
    3. there should be 1.3 fridges per house and 0.3 freezers per house
    4. ELCAP shapes should be normalized before being includes; the magnitudes don't mean much anyway.
  13. daily energy use is missing drying and refrigeration/freezer (table 5)
  14. the load calc equation is not shown
  15. needs to explain heating temperature COP discount

Multi family

  1. can't tell how many apts exclude dryer/washer or second fridge/freezer
  2. electrification is generally higher in apts

Commercial

  1. shouldn't assume that 100% of interior lighting is not T8 (electronic ballast)