Merge pull request #44 from cse-sim/pv-power-electronics

Add ability to model PV Power Electronics

Author: nealkruis

CMakeLists.txt

@@ -14,8 +14,7 @@ set_property(GLOBAL PROPERTY USE_FOLDERS ON)
 
 add_subdirectory(vendor)
 add_subdirectory(src)
-add_subdirectory(bin)
 
 include(ctest)
 
-add_subdirectory(test)
+add_subdirectory(test)

bin/CMakeLists.txt

@@ -91,8 +91,6 @@ The type of array describes mounting and tracking options. Roof mounted arrays h
             OneAxisTracking,
             TwoAxisTracking
 
-<!-- Hide              OneAxisBacktracking \ -->
-
 **pvTilt=*float***
 
 The tilt of the photovoltaic array from horizontal.  Values outside the range 0 to 360 are first normalized to that range. For one-axis tracking, defines the tilt of the rotation axis. Not used for two-axis tracking arrays.  Should be omitted if pvVertices is given.
@@ -135,7 +133,24 @@ Photovoltaic array azimuth (0 = north, 90 = east, etc.). If a value outside the
                                           values
   ------------------------------------------------------------------
 
-  **pvMounting=*choice***
+**pvSIF=*float***
+
+  Shading Impact Factor (SIF) of the array used to represent the disproportionate impact on array output of partially shaded modules at the sub-array level. This impact is applied to the effective beam irradiance on the array:
+
+  \[I_{poa,beam,eff} = \max\left(I_{poa,beam}\cdot\left(1-SIF\cdot f_{sh}\right),0\right)\]
+
+  where $f_{sh}$ is the fraction of the array that is shaded.
+
+  Default value is 1.2, which is representative of PV systems with sub-array microinverters or DC power optimizers. For systems without sub-array power electronics, values are closer to 2.0.
+
+<%= member_table(
+  legal_range: "*x* $\\geq$ 1.0",
+  default: "1.2",
+  required: "No",
+  variability: "constant") %>
+
+
+**pvMounting=*choice***
 
   Specified mounting location of this PVARRAY.  pvMounting=Site indicates the array position is not altered by building rotation via [TOP bldgAzm](#top-general-data-items), while PVARRAYs with pvMounting=Building are assumed to rotate with the building.
 
@@ -200,17 +215,19 @@ Fraction of total DC energy lost. The total loss from a system is aggregated fro
   **Loss Type**             **Default Assumption**
   ------------------------- ----------------------
   Soiling                   0.02
-  Shading                   0.03
+  *Shading*                 *0 (handled explicitly)*
   Snow                      0
-  Mismatch                  0.02
+  *Mismatch*                *0 (shading mismatch handled explicitly [see pvSIF])*
   Wiring                    0.02
   Connections               0.005
   Light-induced degradation 0.015
   Nameplate rating          0.01
-  Age                       0
+  *Age*                     *0.05 (estimated 0.5% degradation over 20 years)*
   Availability              0.03
   **Total**                 **0.14**
 
+*Italic* lines indicate differences from PVWatts assumptions.
+
   ------------------------------------------------------------
   **Units** **Legal** **Default** **Required** **Variability**
             **Range**

bin/LICENSE.txt

@@ -181,7 +181,6 @@ set(headers
   PGPAK.H
   PP.H
   PSYCHRO.H
-  PVCalc.h
   RCCN.H
   RMKERR.H
   SLPAK.H

bin/ssc.dll

@@ -2026,8 +2026,8 @@ CULT( "pvGrndRefl",  DAT,   PVARRAY_GRNDREFL, 0,      0, VHRLY,  TYFL,  0,
 CULT( "pvGCR",       DAT,   PVARRAY_GCR,     0,       0, VEOI,   TYFL,  0,      0.4f,   N, N),
 CULT( "pvDCtoACRatio", DAT, PVARRAY_DCACRAT, 0,       0, VEOI,   TYFL,  0,      1.1f,   N, N),
 CULT( "pvInverterEff", DAT, PVARRAY_INVEFF,  0,       0, VEOI,   TYFL,  0,      0.96f,  N, N),
+CULT( "pvSIF",       DAT,   PVARRAY_SIF,     0,       0, VEOI,   TYFL,  0,      1.2f,   N, N),
 CULT( "pvSysLosses", DAT,   PVARRAY_SYSLOSS, 0,       0, VHRLY,  TYFL,  0,      0.14f,  N, N),
-CULT( "pvUsePVWatts", DAT,  PVARRAY_USEPVWATTSDLL, 0, 0, VEOI,   TYCH,  0,      C_NOYESCH_NO, N, N),
 CULT( "pvMounting",  DAT,   PVG( MOUNTING),  0,       0, VEOI,   TYCH,  0,      C_MOUNTCH_BLDG, N, N),
 CULT( "pvVertices",  DAT,   PVG( VRTINP),    ARRAY,   0, VEOI,   TYFL,  0,      0.f,    v DIM_POLYGONXYZ,  N),
 CULT( "endPVARRAY",  ENDER, 0,               0,       0, 0,      0,     0,      0.f,    N,   N),

doc/src/records/pvarray.md

@@ -419,7 +419,6 @@ DUCTSEGP -- "class DUCTSEG*"	 4	    none
 			FXDOR "FixedOpenRack"
 			FXDRF "FixedRoofMount"
 			1AXT "OneAxisTracking"
-			1AXBT "OneAxisBacktracking"
 			2AXT "TwoAxisTracking"
 		  }
 

src/CMakeLists.txt

@@ -3753,7 +3753,6 @@ RECORD DHWUSE "DHWUse" *RAT	// input / runtime DHW single draw
 RECORD PVARRAY "PVArray" *RAT		// input / runtime photovoltaics array
   *prefix pv_
   *excon
-  *exdes			// define destructor in pvwatts.cpp to free SSC data container
   *ovrcopy
   *declare "virtual PVARRAY& CopyFrom( record* src, int copyName=1, int dupPtrs=0);"
   *declare "void FixUp();"	// virtual fixup after basAnc reAl
@@ -3781,7 +3780,6 @@ RECORD PVARRAY "PVArray" *RAT		// input / runtime photovoltaics array
   *i TI pv_elecMtri;		// meter for system electricity production
   *f ENDUSECH pv_endUse		// end use of energy. defataults to "PV"
   *i FLOAT_GEZ pv_dcCap		// system capacity/size (DC nameplate), kW
-  *i NOYESCH pv_usePVWattsDLL  // use PVWatts DLL instead of CSE calculations
   *i PVMODCH pv_moduleType	// type of module (Standard, Premium, ThinFilm)
   *i FLOAT pv_tempCoeff		// temperature coefficient, 1/F
   *i FLOAT_GE1 pv_covRefrInd  // refraction index for coating applied to cover
@@ -3791,14 +3789,20 @@ RECORD PVARRAY "PVArray" *RAT		// input / runtime photovoltaics array
   *h FRAC pv_grndRefl		// ground reflectance
   *i FRAC pv_gcr			// ground coverage ratio (what fraction of the ground is covered by the array). 1.0 implies no spacing.
   *i FLOAT_GZ pv_dcacRat	// DC to AC ratio
+  *h FLOAT_GE1 pv_sif		// Shading Impact Factor
   *i FRAC pv_invEff			// inverter efficiency at rated power
   *h FRAC pv_sysLoss		// system losses
   *h *e FLOAT pv_tCell		// cell temperature, F
   *h *e FLOAT pv_aoi		// angle of incidence (radians)
   *h *e ANGLE pv_panelTilt	// tilt of pv panel (different from array tilt for tracking systems), radians
   *h *e ANGLE pv_panelAzm	// azimuth of pv panel (different from array tilt for tracking systems), radians
-  *h *e FLOAT pv_poa		// plane of array incidence, Btu/h-ft2
-  *h *e FLOAT pv_poaT		// transmitted plane of array incidence, Btu/h-ft2
+  *h *e FLOAT pv_poa		// plane of array incidence (before shading), Btu/h-ft2
+  *h *e FLOAT pv_poaBeam	// plane of array beam incidence (before shading), Btu/h-ft2
+  *h *e FLOAT pv_radIBeam	// beam radiation incident on array, Btu/h-ft2
+  *h *e FLOAT pv_radIBeamEff	// effective beam radiation incident on array (accounts for shading impact factor), Btu/h-ft2
+  *h *e FLOAT pv_radI		// total radiation incident on array, Btu/h-ft2
+  *h *e FLOAT pv_radIEff	// effective total radiation incident on array (accounts for shading impact factor), Btu/h-ft2
+  *h *e FLOAT pv_radTrans	// transmitted radiation (after accounting for shading impact), Btu/h-ft2
   *h *e FLOAT pv_dcOut		// DC power output, Btu
   *h *e FLOAT pv_acOut		// AC power output, Btu
 
@@ -3807,8 +3811,8 @@ RECORD PVARRAY "PVArray" *RAT		// input / runtime photovoltaics array
   *r FLOAT pv_convRatio		// ratio of back convection to front convection
   *r FLOAT pv_tGrndRatio	// ratio of ground-cell temperature diff. to air-cell temperature diff.
   *r *hide FLOAT_GEZ pv_thermCap  // thermal capacitance of the array, J/kg-m2
-  *h *e FLOAT pv_poaPv		// previous timestep plane of array incidence, Btu/h-ft2
-  *h *e FLOAT pv_tCellPv	// previous timestep cell temperature, F
+  *h *e FLOAT pv_radILs		// last step (curently hour) total radiation incident on array, Btu/h-ft2
+  *h *e FLOAT pv_tCellLs	// last step (curently hour) cell temperature, F
   *r *hide SI pv_modMap		// map to PVWatts module type input
   *r *hide SI pv_arrMap		// map to PVWatts array type input
 

src/CNCULT.CPP

@@ -84,6 +84,9 @@
 
 #include "cse.h"		// decls for this file
 
+#include <penumbra/penumbra.h>	// penumbraInit penumbraTerminate for GPU calculations
+
+
 // configuration defined on compiler command line (or inferred, see cnglob.h)
 #include "cseface.h"
 
@@ -644,6 +647,7 @@ LOCAL INT cse1( INT argc, const char* argv[])
 	}
 	while (argci < argc);		// repeat til all input args used
 
+	Pumbra::penumbraTerminate();	// Clean up GPU calculation memory
 	return errlvl;			// 0 ok, nz error, 255 if ^C, 1-95.
 }			// cse1run
 //------------------------------------------------------------------------