Experimental FeatureThe Metastability model is an experimental feature currently available only to the PlantPredict team for internal testing.
Summary
The Metastability model accounts for performance variations caused by reversible changes in material properties based on recent operating history. It computes a power scale factor from lagged moving averages of cell temperature and POA irradiance. The scale factor is split into voltage and current scalars that adjust the maximum power point output from the single-diode model. The model uses 12 coefficients (C0 through C11).
| Name | Symbol | Units | Description |
|---|
| Metastability Coefficients | C0,C1,...,C11 | various | Model coefficients |
| Metastability Voltage Scale Factor | fV,meta | — | Voltage allocation factor |
| Max Power Voltage | Vmp | V | Voltage at maximum power point |
| Max Power Current | Imp | A | Current at maximum power point |
| Cell Temperature | Tc | °C | Cell temperature at each timestep |
| Effective POA Irradiance | Geff | kW/m² | Total effective POA irradiance (front + rear, after IAM, spectral, soiling, and shading corrections) at each timestep |
Outputs
| Name | Symbol | Units | Description |
|---|
| Adjusted Voltage | Vmp,adj | V | Voltage after metastability adjustment |
| Adjusted Current | Imp,adj | A | Current after metastability adjustment |
| Adjusted Power | Pmp,adj | W | Power after metastability adjustment |
Detailed Description
The metastability correction proceeds in three steps: (1) compute lagged moving averages of cell temperature and POA irradiance that characterize recent operating history, (2) evaluate a polynomial scale factor from those averages, and (3) split the scale factor into voltage and current adjustments applied to the maximum power point.
Lagged Moving Averages
The model uses two look-back windows. The 7-day (168-hour) window captures the thermal and irradiance history that drives slow defect-state transitions, while the 6-hour window captures shorter-term irradiance conditions.
7-day lagged average (168 hours): Only timesteps where Geff≥0.1 kW/m² are included, filtering out nighttime and very low-light periods. The averages are divided by the count n of qualifying timesteps:
Tc,7d=n1i=t−168h∑tTc,iwhere Geff,i≥0.1
G7d=n1i=t−168h∑tGeff,iwhere Geff,i≥0.1
6-hour lagged average: All timesteps are included, and the average is divided by the total number of timesteps in the window:
G6h=n1i=t−6h∑tGeff,i
Scale Factor
The lagged averages feed a polynomial with 12 empirical coefficients (C0 through C11). The polynomial includes linear, quadratic, and cubic terms in Tc,7d and G7d, plus linear and quadratic terms in G6h:
fmeta=C0+C1Tc,7d+C3(Tc,7d−C2)2+C8(Tc,7d−C2)3+C4G7d
+C6(G7d−C5)2+C7(G7d−C5)3+C9G6h+C11(G6h−C10)2
Voltage and Current Adjustment
The scale factor fmeta is split into a voltage scalar and a current scalar using the user-defined voltage allocation factor fV,meta, which controls what fraction of the power change is attributed to voltage versus current:
Vscalar=1−(1−fmeta)⋅fV,meta
Iscalar=Vscalarfmeta
These scalars are then applied to the maximum power point voltage and current from the single-diode model:
Vmp,adj=Vmp×Vscalar
Imp,adj=Imp×Iscalar
Pmp,adj=Vmp,adj×Imp,adj