6.4 Renewable Generation Integration: Solar Photovoltaics (PV) & Wind Power Systems

6.4 Renewable Generation Integration: Solar Photovoltaics (PV) & Wind Power Systems

Executive Overview: Renewable energy systems—primarily solar photovoltaics (PV) and wind power generation—comprise the fastest-growing generation sector connected to transmission and distribution grids. On the NCEES PE Electrical and Computer: Power examination, renewable energy questions focus on photovoltaic semiconductor physics, I-V and P-V curve parameters, temperature and irradiance derating calculations, NEC Article 690 maximum DC string voltage sizing, grid-interconnection standards (IEEE 1547 ride-through and anti-islanding), Betz law wind aerodynamics, and the structural differences among Type 1, 2, 3 (DFIG), and 4 (Full-Converter) wind turbine generators.

---

1. Solar Photovoltaic (PV) Cell Physics & Characteristic Curves

A solar photovoltaic cell is a large-area semiconductor p-n junction diode. When incident photons with energy greater than the semiconductor bandgap ($h\nu > E_g$, where $E_g \approx 1.12\text{ eV}$ for silicon) strike the cell, they generate electron-hole pairs. The built-in electric field of the depletion region sweeps electrons toward the n-side and holes toward the p-side, establishing a photo-generated current ($I_{ph}$).

                 SINGLE-DIODE PV CELL EQUIVALENT CIRCUIT
                     I_ph (Photo Current)
                 +--------->----------+
                 |                    |
               ( | )                +---+ 
              I_ph(^)         I_d   |   | D (Ideal Diode)     R_s (Series Resistor)
                 |          ------> +---+                  +--[ZZZZ]--+-----> + I
                 |                    |                    |          |
                 |                  +---+                +---+        |
                 |                  |   | R_sh (Shunt)   |   |      + | 
                 |                  |   | Resistor       |   |     (~) V
                 |                  +---+                +---+      - | 
                 |                    |                    |          |
                 +--------------------+--------------------+----------+-----> -

Single-Diode Master Current Equation

$I = I_{ph} - I_0 \left[ \exp\left(\frac{q(V + I R_s)}{n k T}\right) - 1 \right] - \frac{V + I R_s}{R_{sh}}$

Where:

• $I_{ph}$ = Light-generated photo-current (proportional to solar irradiance $G$) $[\text{A}]$
• $I_0$ = Diode reverse saturation dark current $[\text{A}]$
• $q = 1.602 \times 10^{-19}\text{ C}$ (electron charge)
• $k = 1.381 \times 10^{-23}\text{ J/K}$ (Boltzmann constant)
• $T$ = Absolute cell temperature in Kelvin ($T [\text{K}] = T [^\circ\text{C}] + 273.15$)
• $n$ = Diode ideality factor ($1.0 \le n \le 1.5$)
• $R_s$ = Internal series resistance (conductor and contact resistance) $[\Omega]$
• $R_{sh}$ = Internal shunt resistance (p-n junction edge leakage) $[\Omega]$

                    PV MODULE I-V AND P-V CHARACTERISTIC CURVES
     Current (A)                                            Power (W)
        ^
    I_sc+-------------------\                                   ^
        |                    \                                  |          P_mp (Max Power Point)
        |                     \  MPP (V_mp, I_mp)               |               /\ 
    I_mp+......................X                                |              /  \ 
        |                      |                                |             /    \ 
        |                      |                                |            /      \ 
        |       I-V Curve      |        P-V Curve               |           /        \ 
        |                      |                                |          /          \ 
        +----------------------+--------------------+->         +---------+------------+-----> Voltage (V)
        0                     V_mp                 V_oc         0        V_mp         V_oc

Key I-V Curve Parameters & Fill Factor

• Open-Circuit Voltage ($V_{oc}$): Voltage across terminals when $I = 0\text{ A}$: $V_{oc} \approx \frac{n k T}{q} \ln\left(\frac{I_{ph}}{I_0} + 1\right)$
• Short-Circuit Current ($I_{sc}$): Current flowing when terminals are shorted ($V = 0\text{ V}$): $I_{sc} \approx I_{ph}$.
• Maximum Power Point (MPP: $V_{mp}, I_{mp}, P_{mp}$): The unique operating point on the I-V curve where the product of voltage and current is maximized: $P_{mp} = V_{mp} \cdot I_{mp}$.
• Fill Factor ($FF$): A figure of merit measuring cell squareness and junction quality ($0.70 \le FF \le 0.85$ for silicon): $\text{FF} = \frac{P_{mp}}{V_{oc} \cdot I_{sc}} = \frac{V_{mp} \cdot I_{mp}}{V_{oc} \cdot I_{sc}}$
• Module Conversion Efficiency ($\eta$): $\eta = \frac{P_{mp}}{G \cdot A_{\text{module}}} = \frac{V_{oc} \cdot I_{sc} \cdot \text{FF}}{G \cdot A_{\text{module}}}$ Where $G$ is solar irradiance (Standard Test Conditions: $G_{STC} = 1,000\text{ W/m}^2$) and $A_{\text{module}}$ is surface area in $\text{m}^2$.

---

2. Environmental Effects: Temperature and Irradiance Coefficients

Photovoltaic performance varies dynamically with ambient temperature and solar irradiance.

+---------------------------------------------------------------------------------------------------+
|                         ENVIRONMENTAL PARAMETER SENSITIVITY MATRIX                                |
+---------------------------------------------------------------------------------------------------+
| Environmental Variable | Impact on Short-Circuit Current (I_sc) | Impact on Open-Circuit Voltage (V_oc) |
| :---                   | :---                                   | :---                                  |
| **Higher Irradiance**  | **Scales Linearly:** $I_{sc} \propto G$ | **Increases Logarithmically:** Minor  |
| **Higher Temperature** | **Increases Slightly:** $+0.03\%/^\circ\text{C}$| **Decreases Dramatically:** $-0.28\%$ to $-0.35\%/^\circ\text{C}$|
+---------------------------------------------------------------------------------------------------+

Cell Operating Temperature Calculation

Because dark PV cells absorb solar radiation, operating cell temperature ($T_{cell}$) exceeds ambient temperature ($T_{amb}$):

$T_{cell} = T_{amb} + \left(\frac{\text{NOCT} - 20^\circ\text{C}}{800\text{ W/m}^2}\right) \times G \quad [^\circ\text{C}]$

Where $\text{NOCT}$ is the Nominal Operating Cell Temperature (typically $45^\circ\text{C} \pm 2^\circ\text{C}$, measured at $T_{amb} = 20^\circ\text{C}$, $G = 800\text{ W/m}^2$, and wind speed $= 1\text{ m/s}$).

Temperature Derating Formulations

• Open-Circuit Voltage at Operating Temperature: $V_{oc}(T) = V_{oc,STC} \times \left[1 + \frac{\beta_{Voc}}{100} \times (T_{cell} - 25^\circ\text{C})\right] \quad [\text{Volts}]$
• Short-Circuit Current at Operating Temperature & Irradiance: $I_{sc}(G, T) = I_{sc,STC} \times \left(\frac{G}{1000}\right) \times \left[1 + \frac{\alpha_{Isc}}{100} \times (T_{cell} - 25^\circ\text{C})\right] \quad [\text{Amperes}]$
• Maximum Power Output at Operating Temperature: $P_{mp}(T) = P_{mp,STC} \times \left[1 + \frac{\gamma_{Pmp}}{100} \times (T_{cell} - 25^\circ\text{C})\right] \quad [\text{Watts}]$ (Typical power temperature coefficient: $\gamma_{Pmp} = -0.35%$ to $-0.45%/^\circ\text{C}$).

---

3. Maximum DC Voltage & String Sizing per NEC Article 690

Under NEC Article 690.7, the maximum photovoltaic system DC voltage ($V_{max}$) must be calculated based on the lowest expected ambient temperature ($T_{min}$) at the installation site. Because $V_{oc}$ rises as temperature falls, cold winter mornings present peak voltage stress that can destroy inverter input stages and violate dielectric ratings ($600\text{ V}$ residential, $1,000\text{ V}$ commercial, $1,500\text{ V}$ utility-scale).

+---------------------------------------------------------------------------------------------------+
|                       NEC 690.7 DC STRING SIZING GOVERNING EQUATIONS                              |
+---------------------------------------------------------------------------------------------------+
| Maximum Single-Module Voltage (Cold):                                                             |
|       V_oc,max = V_oc,STC * [ 1 + (β_Voc / 100) * (T_min - 25°C) ]                                |
|                                                                                                   |
| Maximum Series Modules per String:                                                                |
|       N_series,max = floor( V_dc,max / V_oc,max )                                                 |
|                                                                                                   |
| Minimum Series Modules per String (Hot Summer MPPT Lower Limit):                                 |
|       N_series,min = ceil( V_mppt,min / V_mp,hot )                                                |
+---------------------------------------------------------------------------------------------------+

---

4. Grid-Tied Inverter Architectures & IEEE 1547 Interconnection Standards

+---------------------------------------------------------------------------------------------------+
|                         SOLAR INVERTER TOPOLOGY CLASSIFICATIONS                                   |
+---------------------------------------------------------------------------------------------------+
| 1. Central Inverters (1.0 - 5.0 MW): Massive single inverter for utility-scale solar farms. High  |
|    efficiency (98.5%+), lowest capital cost ($/W), but requires DC combiner boxes and presents a  |
|    single point of failure for large array blocks.                                                |
| 2. String Inverters (25 - 125 kW): Distributed commercial/industrial inverters managing 2 to 12   |
|    strings each. Multiple independent Maximum Power Point Trackers (MPPT) optimize string energy  |
|    yield under partial shading and uneven orientations.                                           |
| 3. Micro-Inverters (300 - 500 W): Module-Level Power Electronics (MLPE) mounted directly to each  |
|    individual PV module, converting DC to AC at the panel. Completely eliminates high-voltage DC  |
|    wiring and natively complies with NEC 690.12 rapid shutdown requirements.                      |
+---------------------------------------------------------------------------------------------------+

IEEE 1547-2018 Interconnection & Interoperability Requirements

IEEE 1547-2018 governs the interconnection of Distributed Energy Resources (DER) with electric power systems:

+---------------------------------------------------------------------------------------------------+
|                             IEEE 1547 MANDATORY GRID FUNCTIONS                                   |
+---------------------------------------------------------------------------------------------------+
| 1. Voltage Ride-Through (LVRT / HVRT): DER inverters must remain connected and ride through      |
|    specified transmission/distribution voltage sags and swells without tripping.                 |
| 2. Frequency Ride-Through (LFRT / HFRT): DER must remain operational across specified grid       |
|    frequency deviations (e.g., 58.5 Hz - 61.2 Hz continuous operation).                           |
| 3. Frequency-Watt (Droop) Control: Inverters automatically curtail active power P during grid    |
|    over-frequency events (f > 60.036 Hz) to stabilize system frequency.                           |
| 4. Volt-VAR Control: Inverters dynamically absorb or inject reactive power Q as a function of    |
|    local bus voltage V to assist in voltage regulation.                                           |
| 5. Anti-Islanding Protection: In the event of an upstream utility feeder trip, the DER must      |
|    detect loss of mains and cease to energize the grid within 2.0 SECONDS to prevent creating    |
|    an energized electrical island dangerous to utility line crews.                               |
+---------------------------------------------------------------------------------------------------+

---

5. Wind Power Aerodynamics & The Betz Limit

Wind turbines extract kinetic energy from moving air masses. The theoretical total power available in an unperturbed wind stream of cross-sectional area $A$ is:

$P_{wind} = \frac{1}{2} \dot{m} v^2 = \frac{1}{2} (\rho A v) v^2 = \frac{1}{2} \rho A v^3 \quad [\text{Watts}]$

Where:

• $\rho$ = Air density (standard sea-level density $\rho = 1.225\text{ kg/m}^3$ at $15^\circ\text{C}$)
• $A = \pi R^2 = \frac{\pi D^2}{4}$ = Rotor swept area $[\text{m}^2]$
• $v$ = Upstream free wind velocity $[\text{m/s}]$

The Betz Limit (Albert Betz, 1919)

An ideal wind turbine cannot extract $100%$ of wind kinetic energy because air must retain velocity to exit the downstream side of the rotor. Applying mass and momentum conservation across an actuator disk yields the maximum theoretical aerodynamic power coefficient ($C_{p,max}$):

$C_{p,max} = \frac{16}{27} \approx 0.5926 \implies 59.26\%$

$\text{Actual Extracted Mechanical Power: } P_{mech} = \frac{1}{2} \rho A v^3 C_p(\lambda, \beta) \quad [\text{Watts}]$

Where $C_p(\lambda, \beta)$ is the turbine aerodynamic power coefficient (typically $0.40 - 0.50$ for modern 3-bladed turbines), $\lambda = \frac{\omega_{rotor} R}{v}$ is the Tip-Speed Ratio (TSR), and $\beta$ is the blade pitch angle.

                    WIND TURBINE OPERATIONAL POWER CURVE
     Power Output (kW)
        ^
  Rated |                                      /-------------------- Region 3: Constant Rated Power
  Power |                                     /                      (Pitch Control Feathering)
        |                                    / 
        |                                   /  Region 2: Variable Speed MPPT
        |                                  /   (Cp maximized)
        |                                 / 
        |                                / 
        +-------------------------------+--------------------+-----> Wind Speed (m/s)
        0                             v_cut-in            v_rated   v_cut-out (25 m/s, Brakes Engage)
                                      (3-4 m/s)          (11-14 m/s)

---

6. Wind Turbine Generator Topologies (Types 1 through 4)

Utility-scale wind turbines utilize four standardized generator topologies classified by IEEE and IEC standards:

+---------------------------------------------------------------------------------------------------+
|                         WIND TURBINE GENERATOR CLASSIFICATION MATRIX                              |
+---------------------------------------------------------------------------------------------------+
| Type | Generator Description           | Power Converter Rating    | Speed Variation Range        |
| :--- | :---                            | :---                      | :---                         |
| **1**| Squirrel Cage Induction (SCIG)  | None (Direct Grid Tie)    | Fixed Speed (\pm 1% slip)   |
| **2**| Wound Rotor Induction (WRIG)    | Dynamic Rotor Resistor    | Limited Variable (0-10%)     |
| **3**| Doubly-Fed Induction (DFIG)     | Partial Scale (~30%)      | Variable Speed (±30%)        |
| **4**| Synchronous / PMSG (Full Conv)  | Full Scale (100%)         | Full Variable (0-100%)       |
+---------------------------------------------------------------------------------------------------+

              TYPE 3: DOUBLY-FED INDUCTION GENERATOR (DFIG)
                     Stator (70% Power) -> Direct to Grid
               +--------------------------------------------------+ (Grid Bus)
               |                                                  |
         [Turbine Rotor] === [Gearbox] === [DFIG Generator]       |
                                                | (Rotor 30%)     |
                                                v                 |
                                         +-------------+          |
                                         | Rotor Conv  |          |
                                         +-------------+          |
                                                | DC Link         |
                                         +-------------+          |
                                         | Grid Conv   | --------+
                                         +-------------+

Topology Details:

1. Type 1 (Fixed-Speed Induction Generator): Standard squirrel-cage induction machine coupled via a multi-stage gearbox directly to the $60\text{ Hz}$ grid. Draws significant lagging reactive power; requires dedicated shunt capacitor banks and provides zero Low-Voltage Ride-Through (LVRT) capability.
2. Type 2 (Variable-Slip Induction Generator): Wound-rotor induction machine with an internal power electronic chopper controlling external rotor resistance (OptiSlip). Dissipates excess rotor power as heat in resistors to achieve narrow speed flexibility ($10%$).
3. Type 3 (Doubly-Fed Induction Generator / DFIG): The stator connects directly to the $60\text{ Hz}$ grid ($70%$ power flow), while the wound rotor is connected through slip rings to a partial-scale back-to-back AC-DC-AC converter ($30%$ rating). Enables full four-quadrant active and reactive power control ($P-Q$), sub-synchronous and super-synchronous operation, and high efficiency at lower converter cost. Requires an active crowbar resistor to protect the rotor converter during grid fault sags.
4. Type 4 (Full-Converter Wind Turbine): Employs a synchronous generator (often a multi-pole Permanent Magnet Synchronous Generator / PMSG eliminating the gearbox, or a standard geared synchronous machine) connected to the grid through a 100% full-scale back-to-back AC-DC-AC power electronic converter. The generator is completely decoupled from grid frequency, providing full variable-speed operation, maximum aerodynamic efficiency, STATCOM-grade reactive power delivery at zero wind, and unmatched Low-Voltage Ride-Through (LVRT) compliance.

---

7. Comprehensive Step-by-Step Worked Mathematical Example

Problem Statement

An electrical engineering firm is designing a commercial rooftop solar PV array and an adjacent utility-scale wind turbine installation.

Part A: Solar PV Array String Sizing (NEC 690.7)

• Module STC Parameters ($1,000\text{ W/m}^2$, $25^\circ\text{C}$):
  • $P_{mp} = 400\text{ W}$, $V_{oc} = 49.5\text{ V}$, $V_{mp} = 41.2\text{ V}$, $I_{sc} = 10.35\text{ A}$, $I_{mp} = 9.71\text{ A}$
  • Temperature coefficient of open-circuit voltage: $\beta_{Voc} = -0.28%/^\circ\text{C}$
  • $\text{NOCT} = 45.0^\circ\text{C}$
• Environmental & Inverter Site Parameters:
  • Lowest expected winter ambient temperature: $T_{min} = -15.0^\circ\text{C}$
  • Highest summer ambient temperature: $T_{max,amb} = 38.0^\circ\text{C}$ with peak irradiance $G = 1,000\text{ W/m}^2$
  • Inverter maximum DC input voltage: $V_{dc,max} = 1,000\text{ V}$
  • Inverter MPPT operating window: $V_{mppt,min} = 550\text{ V}$ to $V_{mppt,max} = 850\text{ V}$

Part B: Wind Turbine Power Generation

• Turbine Parameters:
  • Rotor diameter: $D = 120.0\text{ m}$ (radius $R = 60.0\text{ m}$)
  • Air density: $\rho = 1.225\text{ kg/m}^3$
  • Operating wind speed: $v = 11.0\text{ m/s}$
  • Aerodynamic power coefficient: $C_p = 0.45$
  • Drivetrain gearbox and generator electrical efficiency: $\eta_{total} = 93.5%$

Calculate:

1. The maximum open-circuit voltage per PV module at $T_{min} = -15^\circ\text{C}$ and the maximum number of series modules ($N_{series,max}$) allowed per string.
2. The operating cell temperature ($T_{cell,hot}$) on a peak summer day ($38^\circ\text{C}$) and the resulting string voltage ($V_{string,mp,hot}$) to verify it stays within the inverter's MPPT window ($550 - 850\text{ V}$).
3. The total electrical power in Megawatts ($P_{elec}$) injected into the grid by the wind turbine at $v = 11.0\text{ m/s}$.

=========================================================================================
CALCULATION WORKFLOW & DETAILED STEP-BY-STEP SOLUTION:
=========================================================================================

PART A: SOLAR PV STRING SIZING (NEC 690.7)

Step 1: Compute Maximum Module Open-Circuit Voltage at Coldest Temperature
  Delta_T_cold = T_min - 25°C = -15.0°C - 25.0°C = -40.0°C
  Voltage temperature multiplier:
    Multiplier = 1 + [ (beta_Voc / 100) * Delta_T_cold ]
               = 1 + [ (-0.0028) * (-40.0) ]
               = 1 + 0.1120 = 1.1120

  Maximum module open-circuit voltage:
    V_oc,max = V_oc,STC * Multiplier
             = 49.50 V * 1.1120
             = 55.044 V

Step 2: Determine Maximum Permissible Modules per String
  N_series,max = floor( V_dc,max / V_oc,max )
               = floor( 1,000 V / 55.044 V )
               = floor( 18.167 )
               = 18 modules per string

  Verification of worst-case cold open-circuit string voltage:
    V_string,oc,cold = 18 * 55.044 V = 990.79 V <= 1,000 V (COMPLIANT)

Step 3: Verify Hot Summer Operation within MPPT Voltage Window
  Operating cell temperature at T_amb = 38°C and G = 1,000 W/m²:
    T_cell,hot = T_amb + [ (NOCT - 20°C) / 800 ] * G
               = 38.0°C + [ (45.0 - 20.0) / 800 ] * 1000
               = 38.0°C + [ 25.0 / 800 ] * 1000
               = 38.0°C + 31.25°C
               = 69.25°C

  Delta_T_hot = 69.25°C - 25.0°C = +44.25°C
  Module MPP voltage at hot cell temperature:
    V_mp,hot = V_mp,STC * [ 1 + (beta_Voc / 100) * Delta_T_hot ]
             = 41.20 V * [ 1 + (-0.0028) * 44.25 ]
             = 41.20 V * [ 1 - 0.1239 ]
             = 41.20 V * 0.8761
             = 36.095 V

  String MPP voltage under hot summer conditions (18 modules):
    V_string,mp,hot = 18 * 36.095 V = 649.71 V
    Since 550 V <= 649.71 V <= 850 V, the string operates squarely inside the
    inverter MPPT tracking window. (DESIGN VERIFIED)

-----------------------------------------------------------------------------------------
PART B: WIND TURBINE POWER GENERATION

Step 4: Compute Rotor Swept Area and Available Wind Kinetic Power
  Rotor swept area:
    A = pi * R^2 = pi * (60.0 m)^2 = 3,600 * pi = 11,309.73 m²

  Available wind power at v = 11.0 m/s:
    v^3 = (11.0)^3 = 1,331.0 m³/s³
    P_wind = 0.5 * rho * A * v^3
           = 0.5 * (1.225 kg/m³) * (11,309.73 m²) * (1,331.0 m³/s³)
           = 0.6125 * 15,053,250.63
           = 9,220,116 W = 9.2201 MW

Step 5: Compute Mechanical and Injected Grid Electrical Power
  Mechanical power extracted by rotor (Cp = 0.45):
    P_mech = P_wind * C_p
           = 9.2201 MW * 0.45
           = 4.14905 MW

  Electrical power output to grid (eta_total = 93.5%):
    P_elec = P_mech * eta_total
           = 4.14905 MW * 0.935
           = 3.87936 MW ≈ 3.88 MW
=========================================================================================

---

8. Common PE Exam Traps & Tactical Pitfalls

• Applying Cold Temperature Factor to $V_{mp}$ Instead of $V_{oc}$: NEC 690.7 explicitly mandates that maximum system voltage checks use $V_{oc}$ (open-circuit voltage). Using $V_{mp}$ severely underestimates maximum dielectric stress, creating code violations.
• Omitting Wind Velocity Cubed Relationship: Power in the wind scales with $v^3$, not $v$ or $v^2$. If wind speed increases by $20%$ ($1.20\times$), total available wind power increases by $(1.20)^3 = 1.728\times$ ($+72.8%$ increase!).
• Confusing DFIG (Type 3) Converter Rating with Generator Rating: In a Type 3 DFIG, the power electronic converter handles only rotor slip power (approximately $25% - 30%$ of machine rating), whereas in a Type 4 turbine, the converter must be rated for $100%$ full machine power.
• Anti-Islanding Clearing Time Mandate: IEEE 1547 specifies that upon complete loss of the grid reference, distributed energy resources must cease to energize the local island within a strict maximum window of 2.0 seconds.