Renewable Energy Systems

Module: Sustainable Tech | Difficulty: Beginner

Overview

Explore the fundamentals of renewable energy including solar, wind, and hydroelectric power systems with performance modeling.

Learning Objectives

• Understand the core scientific principles underlying this sustainable technology.
• Apply quantitative methods to analyze system performance and efficiency.
• Implement Python-based models for simulation and optimization.
• Evaluate environmental impact using sustainability metrics.
• Design solutions that integrate renewable energy and resource conservation.

System Architecture

Solar Energy Fundamentals

Solar photovoltaic systems convert sunlight directly into electricity using semiconductor materials.

where is the panel efficiency, is panel area in , and is global horizontal irradiance in .

import numpy as np

def solar_power_output(efficiency, area, irradiance):
    return efficiency * area * irradiance

irradiance_profile = np.random.normal(500, 150, 8760)
irradiance_profile = np.clip(irradiance_profile, 0, 1200)

eta = 0.20
area = 40.0
hourly_power = solar_power_output(eta, area, irradiance_profile)
annual_energy_kwh = np.sum(hourly_power) / 1000
print(f"Annual energy: {annual_energy_kwh:.0f} kWh")

Wind Energy Conversion

The Betz limit restricts maximum extractable power to :

def wind_power(rho, area, velocity, cp=0.45):
    betz_limit = 16 / 27
    cp = min(cp, betz_limit)
    return 0.5 * rho * area * velocity**3 * cp

air_density = 1.225
rotor_radius = 40
swept_area = np.pi * rotor_radius**2
wind_speeds = np.random.weibull(2.0, 8760) * 8

power_curve = [wind_power(air_density, swept_area, v) for v in wind_speeds]
annual_wind_energy = sum(power_curve) / 1000
print(f"Wind annual output: {annual_wind_energy:.0f} kWh")

Hydroelectric Power

Large-scale hydro achieves efficiencies above .

Renewable Mix Optimization

The objective function minimizes total operating cost subject to power balance constraints:

Capacity constraints enforce generation limits for each source:

Ramp rate constraints limit how quickly generation can change:

import numpy as np

def renewable_dispatch(solar_cf, wind_cf, demand, solar_cap, wind_cap):
    """Simple merit-order dispatch for renewable mix."""
    solar_gen = solar_cf * solar_cap
    wind_gen = wind_cf * wind_cap
    total_renewable = solar_gen + wind_gen
    met_by_renewable = np.minimum(total_renewable, demand)
    curtailment = np.maximum(total_renewable - demand, 0)
    deficit = demand - met_by_renewable
    penetration = np.mean(met_by_renewable) / np.mean(demand)
    return penetration, np.sum(curtailment), np.sum(deficit)

pen, curt, defc = renewable_dispatch(
    solar_cf=0.25, wind_cf=0.35, demand=np.full(8760, 5000),
    solar_cap=8000, wind_cap=6000
)
print(f"Renewable penetration: {pen:.1%}")
print(f"Annual curtailment: {curt/1000:.0f} MWh")

Levelized Cost of Energy

LCOE provides a metric for comparing generation technologies:

where is total cost in year , is energy production, is the discount rate, and is project lifetime.

Grid Integration Challenges

Renewable integration introduces variability. The duck curve illustrates the mismatch between solar generation peaks and evening demand peaks.

Net load ramping requirements can exceed in regions with high solar penetration. Flexible resources including battery storage, demand response, and fast-ramping gas plants are essential for maintaining grid stability.

import numpy as np

def duck_curve_analysis(solar_penetration):
    """Model the duck curve with increasing solar penetration."""
    hours = np.arange(24)
    base_demand = 15 + 5 * np.sin(2 * np.pi * (hours - 6) / 24)
    solar_profile = np.maximum(0, np.sin(2 * np.pi * (hours - 6) / 24))
    solar_gen = solar_penetration * solar_profile
    net_load = base_demand - solar_gen
    ramp = np.diff(net_load)
    max_ramp_up = np.max(ramp)
    min_net_load = np.min(net_load)
    return max_ramp_up, min_net_load

for sp in [0.1, 0.3, 0.5]:
    ramp, minimum = duck_curve_analysis(sp)
    print(f"Solar {sp:.0%}: max ramp {ramp:.2f} GW/hr, min net load {minimum:.1f} GW")

Hands-On Exercise

import numpy as np

def optimize_mix(solar_cap, wind_cap, hydro_cap, demand):
    solar_gen = np.random.uniform(0, solar_cap, 8760)
    wind_gen = np.random.weibull(2.0, 8760) * wind_cap * 0.3
    hydro_gen = np.full(8760, hydro_cap * 0.8)

    total_gen = solar_gen + wind_gen + hydro_gen
    curtailment = np.maximum(total_gen - demand, 0)
    deficit = np.maximum(demand - total_gen, 0)

    capacity_factor = np.mean(total_gen) / np.max(total_gen)
    self_sufficiency = 1 - np.sum(deficit) / np.sum(demand)

    return capacity_factor, self_sufficiency, np.sum(curtailment)

result = optimize_mix(5000, 3000, 2000, np.full(8760, 4000))
print(f"Capacity factor: {result[0]:.2%}")
print(f"Self-sufficiency: {result[1]:.2%}")

Key Takeaways

1. Solar irradiance and wind speed variability directly affect renewable energy output.
2. The Betz limit of 59.3% constrains maximum wind energy extraction.
3. Hydropower achieves the highest conversion efficiencies among renewables.
4. Hybrid systems combining multiple sources improve reliability and reduce curtailment.
5. Grid integration requires sophisticated forecasting and flexible dispatch strategies.

Review Questions

1. What factors determine the capacity factor of a solar installation?
2. How does the cubic relationship between wind speed and power affect turbine siting?
3. Explain the duck curve and its implications for grid operators.
4. Compare the LCOE across solar, wind, and hydro sources.
5. Design a control strategy for a hybrid renewable microgrid.

References and Further Reading

• Jacobson, M. Z. et al. (2017). 100% Clean and Renewable Wind, Water, and Sunlight.
• IRENA (2023). Renewable Power Generation Costs in 2022.
• Smil, V. (2017). Energy Transitions: Global and National Perspectives.
• Masters, G. M. (2013). Renewable and Efficient Electric Power Systems.