Mathematical Models Driving Solar, Wind, and Storage Systems

Introduction

As global energy structures shift toward carbon neutrality, mathematical modeling becomes essential for designing, optimizing, and managing renewable energy systems. Models provide scientific bases for photovoltaic efficiency prediction, wind power output optimization, and storage capacity allocation.

1. Solar Photovoltaic System Modeling

1.1 Solar Irradiance Model

The horizontal‑plane solar irradiance can be approximated by:

I_h = I_{sc} \cdot r_n \cdot (\cos\theta_z + \text{diffuse} + \text{ground\,reflected})

where I_{sc} is the solar constant (~1367 W/m²), r_n is the Earth‑Sun distance correction factor, and \theta_z is the solar zenith angle calculated from latitude, solar declination, and hour angle. For tilted arrays, tilt angle β and azimuth γ introduce additional geometric factors for direct, diffuse, and reflected components.

1.2 PV Cell Conversion Efficiency Model

The single‑diode equivalent circuit yields the output current:

I = I_{ph} - I_0\left(\exp\frac{V+IR_s}{nV_t}-1\right) - \frac{V+IR_s}{R_{sh}}

Key parameters include photo‑current I_{ph} (proportional to irradiance), diode saturation current I_0, series resistance R_s (0.1–1 Ω), shunt resistance R_{sh} (100–1000 Ω), ideality factor n (1–1.5), and thermal voltage V_t = kT/q (~26 mV at 25 °C). Temperature effects are modeled with the NOCT method, adjusting efficiency by a temperature coefficient of about –0.0045 /°C.

2. Wind Power System Modeling

2.1 Wind Turbine Power Characteristics

The power captured by a wind turbine is: P_{wind}=\frac{1}{2}\rho A v^3 C_p where \rho is air density (~1.225 kg/m³), A the rotor swept area, v the upstream wind speed, and C_p the power coefficient (theoretical maximum 0.593, typical 0.4–0.5 after empirical fitting).

2.2 Weibull Wind Speed Distribution

Long‑term wind speed is often described by a two‑parameter Weibull distribution with probability density:

f(v)=\frac{k}{\lambda}\left(\frac{v}{\lambda}\right)^{k-1}\exp\left[-\left(\frac{v}{\lambda}\right)^k\right]

k

is the shape parameter (1.5–3) and \lambda the scale parameter (m/s). The mean wind speed and annual energy production are obtained by integrating the turbine power curve over this distribution, considering cut‑in and cut‑out speeds and the total operating hours (8760 h per year).

3. Energy Storage Capacity Optimization

3.1 Battery Storage Dynamic Model

The state of charge (SOC) evolves as:

SOC_{t+1}=SOC_t+\frac{\eta_{ch} P_{ch}\Delta t - \frac{1}{\eta_{dis}} P_{dis}\Delta t}{E_{rated}}

where \eta_{ch} and \eta_{dis} are charge/discharge efficiencies (0.95–0.98), P_{ch} / P_{dis} are power flows (kW), and E_{rated} is the rated capacity (kWh). Self‑discharge rate is typically 0.001–0.005 / day. Cycle life relates to depth of discharge (DOD) via a power‑law relationship.

3.2 Capacity Allocation Optimization

Capacity sizing is formulated as a multi‑objective problem minimizing lifecycle cost and maximizing system efficiency, subject to constraints such as:

No simultaneous charging and discharging.

Power balance between generation, load, and storage.

Device capacity limits.

4. Microgrid Energy Management System

4.1 Energy Balance and Power Flow

The instantaneous power balance for a microgrid can be expressed as: \sum P_{gen} - \sum P_{load} - P_{storage}=0 Positive storage power denotes discharge, negative denotes charge. System losses include inverter, line, and transformer losses.

4.2 Economic Optimal Dispatch

A multi‑period objective combines electricity purchase/sale prices, storage degradation cost, and carbon emission cost:

\min \sum_{t}\big( c_{buy}\,P_{buy}(t) - c_{sell}\,P_{sell}(t) + c_{deg}\,\Delta SOC(t) + c_{CO2}\,E_{CO2}(t) \big)

Constraints enforce power balance, equipment capacity, and SOC limits.

5. Prediction and Control Models

5.1 ARMA Time‑Series Forecasting

Short‑term PV and wind power can be forecast with an ARMA(p,q) model:

X_t = \mu + \sum_{i=1}^{p}\phi_i X_{t-i} + \sum_{j=1}^{q}\theta_j \epsilon_{t-j} + \epsilon_t

If the series is non‑stationary, differencing of order d yields an ARIMA(p,d,q) model. Exogenous variables (e.g., temperature, irradiance) lead to ARIMAX formulations.

5.2 Model Predictive Control (MPC)

MPC solves a finite‑horizon optimization at each control step. The standard objective is:

\min_{u}\sum_{k=0}^{N_p}\|y_{k|t}-r_k\|_Q^2 + \|\Delta u_{k|t}\|_R^2

where N_p is the prediction horizon, N_c the control horizon, y the predicted outputs, r the reference trajectory, and Q,R are positive‑definite weighting matrices. Constraints include power limits, SOC bounds, and ramp rates.

With the rise of artificial intelligence and big‑data techniques, data‑driven models such as deep learning and reinforcement learning are increasingly combined with physical models, enabling multi‑scale, multi‑physics, and uncertainty‑aware simulations that further improve renewable energy system efficiency, reliability, and economics.