Least Cost Vehicle Charging in a Smart Neighborhood Considering Uncertainty and Battery Degradation

Abstract

The electricity landscape is constantly evolving, with intermittent and distributed electricity supply causing increased variability and uncertainty. The growth in electric vehicles, and electrification on the demand side, further intensifies this issue. Managing the increasing volatility and uncertainty is of critical importance to secure and minimize costs for the energy supply. Smart neighborhoods offer a promising solution to locally manage the supply and demand of energy, which can ultimately lead to cost savings while addressing intermittency features. This study assesses the impact of different electric vehicle charging strategies on smart grid energy costs, specifically accounting for battery degradation due to cycle depths, state of charge, and uncertainties in charging demand and electricity prices. Employing a comprehensive evaluation framework, the research assesses the impacts of different charging strategies on operational costs and battery degradation. Multi-stage stochastic programming is applied to account for uncertainties in electricity prices and electric vehicle charging demand. The findings demonstrate that smart charging can significantly reduce expected energy costs, achieving a 10% cost decrease and reducing battery degradation by up to 30%. We observe that the additional cost reductions from allowing Vehicle-to-Grid supply compared to smart charging are small. Using the additional flexibility aggravates degradation, which reduces the total cost benefits. This means that most benefits are obtainable just by optimized the timing of the charging itself.

1. Introduction

Electricity supply and demand management are becoming increasingly complex. This complexity arises from the increasing contribution of decentralized renewable energy sources (RESs), like solar and wind power, feeding into local distribution grids [1]. Furthermore, the increase in electricity consumption types, such as electric vehicles, often resulting in high demand peaks, further adds to these challenges [2]. As this trend advances, the integration of distributed RESs through optimal energy management within smart buildings and neighborhoods is gaining heightened importance [3,4]. The feasibility of intermittent RESs such as solar photovoltaic (PV) as standalone power sources is limited due to the inherent fluctuations in electricity generation [5]. A straightforward solution to address RES intermittency is by using storage units. Co-generation units and (legacy) connections to electricity and heat grids may help balance generation peaks and ease the integration of RES. Dispatchable energy units such as storage units for heat and power, and combined heat and power (CHP) are especially important in disconnected energy communities [6]. The increased flexibility benefits all actors in the energy system at large; however, it comes at the expense of significant investment costs. The optimal use of storage capacity allows investing in smaller, cheaper units. Moreover, considering battery degradation allows trading off arbitrage opportunities with battery lifetime.

The increasing number of electric vehicles (EVs), most of which are idle over $90\%$ of the time, may have a role as a flexible source in terms of shiftable demand and storable supply. Managing charging schedules can be beneficial not only for the charging costs of the EVs but also the operational cost of the entire system. Within this landscape, two distinct paradigms come into focus: passive, uncontrolled charging, and active, flexible charging (e.g., [7]). In uncontrolled strategies, vehicles are charged upon plugging in, charging at their maximum capacity until the battery is fully charged. Conversely, active strategies align with electricity price signals and broader system costs, potentially encompassing battery degradation considerations. Taking active charging a step further is the concept of Vehicle-to-Grid (V2G), which represents the most flexible utilization of EV batteries. V2G allows for discharges too, albeit at the cost of added stress and a reduction in battery lifetime [7].

Battery lifetime consists of two distinct components: calendar life and cycle life, both of which contribute to capacity reduction through chemical and mechanical stresses [8]. For batteries that are used frequently, cycle life is often the key factor affecting lifespan [9]. In contrast, for EV batteries, which may only be charged a few times per week, calendar aging is usually the dominant factor. Avoiding battery degradation in techno-economic models of smart grids or energy systems is commonly achieved by imposing a cost per kWh throughput and establishing (often conservative) boundaries for state of charge (SOC) [10,11]. However, these models often overlook the non-linear increase in aging caused by cycle depth (CD) and other aging factors. This oversight typically results in relatively deep charge–discharge cycles and high resting SOC levels, leading to substantially higher battery degradation than the model predicts. Alternatively, excessively penalizing smaller cycles or setting overly tight SOC boundaries can hinder leveraging arbitrage opportunities due to inter-temporal price differences. Several studies focus specifically on modeling battery degradation. For instance, refs. [9,12] analyze the impact of degradation on electric vehicle batteries and its associated cost implications. Similarly, ref. [13] investigates the effects of battery degradation in electric vehicles and stationary batteries on the investment costs of distribution networks for renewable energy-based distributed generation, energy storage systems, and electric vehicle charging stations. This study only includes the impacts of depth of discharge on the degradation cost. However, many other studies examining the interplay between EVs and smart neighborhoods oversimplify the degradation costs, often reducing them to a constant value [14,15,16,17].

Moreover, the challenge of (foreseeable) intermittency, as well as uncertainty arising from sources such as RES, market prices, and demand, impacts least-cost operation [14,15,18,19]. Uncertainty in electricity prices and demand impacts least-cost operation and is a key challenge in energy systems. For instance, virtual power plants have been proposed to aggregate small-scale battery energy storage systems for frequency regulation, demonstrating how managing uncertainties in market prices and storage degradation is crucial for their operational success [20]. Despite the significant impact of uncertainty, most researchers deploy deterministic models (cf., [21,22,23]). An exception is [24], who propose a two-stage stochastic model aimed at optimizing the operation of a RES-based micro grid (MG) equipped with an EV parking lot. However, in energy systems, uncertain events typically unfold across multiple temporal instances, making stochastic programs with multiple (more than two) stages more appropriate [25,26]. At present, the number of multi-stage models for operational scheduling in energy systems is relatively limited. Most of the smart grid literature focuses on deterministic or two-stage stochastic models to deal with the optimal operational problems of smart energy neighborhoods [16,17,27]. However, ref. [28] employs a multi-stage stochastic model to enhance energy-efficient scheduling of a CHP unit, PV, heat boiler, and energy storage systems within a building, addressing both electricity and thermal demands. Ref. [29] applies stochastic dual dynamic programming and compares a two-stage and a multi-stage model to minimize the expected total cost of an IEEE 118-bus system with integrated RES. Additionally, many papers overlook or simplify the degradation costs associated with stationary storage and EV batteries. With this paper, we aim to contribute to the literature by developing a multi-stage stochastic model of a local energy system with a detailed battery degradation representation.

This paper investigates how electric vehicle charging strategies interact with behavioral and price uncertainties in a building complex equipped with various energy generation and storage options. We aim to minimize the expected operational costs, including battery degradation costs. To achieve this, we adopt the improved battery degradation model as developed in [30]. In [30], a battery and a spot market for electricity are the two considered elements. The multi-stage operational model in the current paper considers a system with many more elements, and it provides a comprehensive framework for accommodating uncertainty in EV owners’ behavior (charging demand) and electricity prices.

This study is conducted to assess the impact of different charging strategies for EVs on the operational energy costs while taking into account battery degradation. The method is applied to a neighborhood within an urban setting, representing a building complex in Pfreimd, Germany. It consists of four buildings with a total of 80 apartments and approximately 300 residents. The complex has its own energy resources and connections to the main energy grids for power, heat, and natural gas. The buildings are outfitted with PV modules. One of the buildings features a small-scale CHP, a TES, and a stationary battery. External electricity and district heating grids both supply energy to this community, and surplus electricity can be fed back into the grid; however, heat cannot. In this analysis, we account for battery degradation in both the stationary battery and the EVs batteries with the representations developed in [30]. We also consider uncertainty in electricity prices and demand. In this study, we consider nickel manganese cobalt (NMC) batteries only due to their widespread usage in electric vehicles [31]. Due to the lacking data, other battery types, such as lithium iron phosphate batteries, as promising as they may be in light of sustainability and ethical considerations [32], are not in the scope of the current paper.

In summary, the key contributions of this paper are as follows:

1. The interaction between electric vehicle charging strategies, behavioral uncertainty, and price fluctuations in a building complex equipped with various energy generation and storage options is analyzed.

2. A comprehensive framework is provided for considering both electric vehicle and stationary battery degradation, based on calendar aging, cycle depth, and SOC levels, in the context of minimizing operational costs.

3. A multi-stage stochastic operational model is developed to minimize the expected operational costs, including battery degradation costs, by accounting for uncertainty in EV charging demand and electricity prices.

4. Three charging strategies—passive charging, smart charging, and V2G charging—are considered to investigate the impact on degradation costs and evaluate the differences in operational energy costs for each strategy.

5. We show that smart charging achieves a significant reduction in energy cost and battery degradation, and that most benefits are obtainable just by the optimized timing of the charging itself; Vehicle-to-Grid delivery adds relatively little value.

The remainder of the paper is structured as follows. In Section 2, we explain the scenario tree structure used to represent uncertainty and present the mathematical formulation of the proposed model, including the cost function, energy balances, and operational constraints of generation units as well as storage facilities. Section 2.4 provides an electric battery degradation model used in the optimization, detailing three degradation mechanisms. Following this, Section 3 discusses the numerical results and findings of the conducted case study. Lastly, Section 4 summarizes our findings and identifies future research directions.

2. Model Formulation

In this section, we present the mathematical framework developed to evaluate the impact of EV charging strategies within a complex building environment as discussed in Section 1. Consider B to represent the set of buildings within the complex, with b indexing the individual buildings. Each building b has distinct heat and electricity demands. To supply this demand, the model considers two types of grids: local, L, and main, M. The model includes three energy carriers within set E: electricity, heat, and gas denoted as $e=el,ht,gs$, respectively. Gas is used exclusively as an input for the CHP. In the following subsection, we begin by introducing the scenario tree structure, a key component for addressing uncertainty.

2.1. Structure of Scenario Trees

To account for stochastic scenarios, we utilize a tree structure. A scenario tree, which features interconnected nodes h, is similar to a decision tree but differs in its branching. A scenario tree branches due to uncertain events, not based on distinct decisions, cf., Figure 1.

All scenarios start from the root node (the first hour in the planning horizon, node 1 in the example figure, Figure 1), with each scenario leading to a unique terminal node ($\widehat{h}$, a possible realization of the last hour in the planning horizon, cf., one of the nodes 20–23 in the example figure). In a scenario tree, predecessor and successor nodes determine the possible paths from the root node to all terminal nodes. Each node in the scenario tree has one unique predecessor but may have multiple successors (e.g., node 6 has predecessor node 4 and successors 8 and 9; we use $\tilde{h}$ to denote predecessors, e.g., in Figure 1, $\tilde{6}$ is node 4). A ‘stage’ in the tree represents a specific time period and may include several parallel nodes in different scenarios (e.g., nodes 4 and 5 are both in stage three). We allocate probabilities ${\alpha }_h$ to each node and ensure their sum equals 1 in each stage to enable calculating the expected values. This structure systematically represents the outcomes of uncertain events, enabling decisions to anticipate different future outcomes, each weighted by their assigned probabilities. In any scenario tree node, outcomes and decisions in its predecessor nodes are known and can be accounted for. In Figure 1), decisions in all nodes 1–7 consider multiple future scenarios. However, decisions in nodes with number 8 and higher do not face uncertainty anymore.

2.2. Objective Function

The objective function is to minimize the expected energy costs as shown in Equation (1):

(1)$\begin{array}{l}\min {\displaystyle \sum _{h\in H}}{\alpha }_h·\left[{\displaystyle \sum _{b\in B,e\in E}}(p_{e,h}·G_{b,e,M,h}\right.-v_{e,M}·(p^{RES}·I_{b,e,M,h}^{RES}+p^{CHP}·I_{b,e,M,h}^{CHP}))\\ \left.+{\displaystyle \sum _{b\in B}}(p_{gs}·CH{P}_{b,gs,h}^{IN}+k^{CHP}·CH{P}_{b,el,h}^{OUT})+{\displaystyle \sum _{b\in B,w\in W}}{K}_{b,w,el,h}\right]\\ +{\displaystyle \sum _{b\in B,w\in W,\widehat{h}\in H}}{\alpha }_h·(SO{C}_{b,w,el,0}-SO{C}_{b,w,el,\widehat{h}})·k^{EOH}\end{array}$

The objective function spans three lines. The first term on the first line represents the cost of purchasing energy from the main grid ($G_{b,e,M,h}$) at the price of $p_{e,h}$ for each building b and energy carrier e at hour h. The second term on the first line calculates the feed-in remuneration from injecting energy into the main grids, accounting for transmission losses ($v_{e,M}$) and with varying tariffs depending on the energy source. It includes remuneration for the RES at a rate of $p^{RES}$ and for CHP generation at a rate of $p^{CHP}$, with associated injections represented by $I_{b,e,M,h}$. On the second line, the first term represents the costs of feed gas for the CHP at price $p_{gs}$ and the operational costs $k^{CHP}$ for each $kW{h}_{el}$ produced. Additionally, the second terms $K_{b,w,h}$ represents the costs associated with battery degradation type w for each building b and hour h. The term on the third line is a penalty for any deviation in SOC levels at the end of the scheduling period from their initial values. In active charging schemes, deviations between the end-of-horizon and starting SOC levels are allowed but penalized by a value $k^{EOH}$ (positive deviations are credited by the same value per kWh as negative deviations are penalized with). Values for parameters can be found in Appendix A. The value we have used for $k^{EOH}=p_{el,0}$, to connect the value of the electricity in the battery at the end of the day to its value at the start of the day. See [33] for further details. This allows flexibility in the response to significant electricity price variations at similar hours on different days. Each value is weighted according to its probability, denoted as ${\alpha }_h$. Since our planning horizon is short, we do not apply discounting.

2.3. Constraints

This section presents the various types of constraints.

2.3.1. Energy Balances

Equations (2) and (3) detail the energy balances for the buildings and local grids, respectively. These balances include electricity and heat, covering the supply from RES and CHP units. They also account for interactions with local and main grids, as well as storage. We factor in grid exchange losses during energy injection. There is no balance for gas; its usage is completely determined by the CHP output as discussed below:

(2)$\begin{array}{cc}& {\displaystyle \sum _{t\in T}}{g}_{b,t,e,h}^{RES}+CH{P}_{b,e,h}^{OUT}+{\displaystyle \sum _{n\in N}}{G}_{b,e,n,h}=\hfill \\ & de{m}_{b,e,h}+C{U}_{b,e,h}+{\displaystyle \frac{I_{b,e,M,h}}{v_{e,M}}}+I_{b,e,L,h}\hfill \\ & -{\displaystyle \sum _{w\in W,j\in J_e}}{D}_{b,w,e,j,h}+{\displaystyle \sum _{w\in W,j\in J_e}}{C}_{b,w,e,j,h}\hfill & & \forall b\in B,e=\{el,ht\},h\in H\hfill \end{array}$

(3)$\begin{array}{cccc}& {\displaystyle \sum _{b\in B}}{G}_{b,e,L,h}=v_{e,L}·{\displaystyle \sum _b}{I}_{b,e,L,h}\hfill & & \forall e=el,ht,h\in H\hfill \end{array}$

The upper limit for grid injections from $RES$ is set by the RES production itself as shown in Equation (4). The total injections into the main grid are the sum of those from CHP and RES, detailed in Equation (5). Any excess heat is either stored in the thermal energy storage (TES) or curtailed; it cannot be injected into the main grid (cf., Equation (6)):

(4)$\begin{array}{cccc}& I_{b,e,M,h}^{RES}\le v_{e,M}·{\displaystyle \sum _{t\in T}}{g}_{b,t,e,h}^{RES}\hfill & & \forall b\in B,e=\{el,ht\},h\in H\hfill \end{array}$

(5)$\begin{array}{cccc}& I_{b,e,M,h}=I_{b,e,M,h}^{CHP}+I_{b,e,M,h}^{RES}\hfill & & \forall b\in B,e=\{el,ht\},h\in H\hfill \end{array}$

(6)$\begin{array}{cccc}& I_{b,ht,M,h}=0\hfill & & \forall b\in B,h\in H\hfill \end{array}$

2.3.2. CHP Unit

To model CHP units and battery storage, further down in the paper, piecewise linearization is used. Line segments correspond to the operating ranges. We apply an $SOS2$ approach which means that we use binary variables to select the active operating range, while weighting variables determine the actual operating level on that segment. The line segmentation is defined by breakpoints, including an off state, a minimum operating level, and some other operating levels. Equation (7) ensures that only one line segment can be active at a time. Equation (8) ensures that the breakpoint representing the off state can only be active alone. Equation (9) guarantees that these weights can only be positive for the selected active segment. Finally, Equation (10) ensures that the sum of the weights $X_{u,h}$ determining the operating level equals 1. The energy input and output for CHP units are then defined by Equations (11) and (12):

(7)$\begin{array}{cccc}& B_{u,h}^{CHP}+B_{u+2,h}^{CHP}\le 1\hfill & & \forall u\in U,h\in H\hfill \end{array}$

(8)$\begin{array}{cccc}& B_{1,h}^{CHP}+B_{u,h}^{CHP}\le 1\hfill & & \forall u\in U|u>1,h\in H\hfill \end{array}$

(9)$\begin{array}{cccc}& X_{u,h}\le B_{u,h}^{CHP}\hfill & & \forall u\in U,h\in H\hfill \end{array}$

(10)$\begin{array}{cccc}& {\displaystyle \sum _{u\in U}}{X}_{u,h}=1\hfill & & \forall h\in H\hfill \end{array}$

(11)$\begin{array}{cccc}& CH{P}_{b,h}^{IN}={\displaystyle \sum _{u\in U}}{o}_{u,gs}^{CHP}·X_{u,h}\hfill & & \forall b\in B,h\in H\hfill \end{array}$

(12)$\begin{array}{cccc}& CH{P}_{b,e,h}^{OUT}={\displaystyle \sum _{u\in U}}{o}_{u,e}^{CHP}·X_{u,h}\hfill & & \forall b\in B,e=\{el,ht\},h\in H\hfill \end{array}$

2.3.3. Energy Storage Devices

The stored energy level for both electrical and thermal storages at the start of a time step equals the level of the previous time step adjusted for loss-corrected charges or discharges as indicated by Equation (13). In this context, $\tilde{h}$ denotes the immediate predecessor node of h in the scenario tree. For arriving EVs, the storage level from the previous time-step indicates the SOC upon arrival. To assure a fully charged EV at departure, we impose $S_{b,EV,el,j,h^{DEP}}=e_{b,EV,el}^{RAT}$. For passive charging, we fix appropriate values for charging levels $C_{b,EV,el,j,h}$. Naturally, for smart charging and V2G, $C_{b,EV,el,j,h}$ is determined by the model as part of the cost minimization. While the discharge $D_{b,EV,el,j,h}=0$ for both passive and smart, for V2G $D_{b,EV,el,j,h}$, it is a model-determined variable.

In the following, we extend the battery degradation formulation developed in [30], incorporating additional indices to account for a more granular problem setting as well as uncertainty. For non-linear CD-based degradation, we follow the approach to apply piecewise linearization on virtual electric battery segments $J_e$, with equal sizes and segment-specific costs. We do not consider the degradation of TES. To obtain a more succinct presentation, the equations below apply for heat storage with one segment only: $|J_{ht}|=1$.

Equation (14) enforces the storage capacity limits (at the segment level):

(13)$\begin{array}{cccc}& \begin{array}{c}\hfill S_{b,w,e,j,h}=S_{b,w,e,j,\tilde{h}}+v_{b,w,e}^c·C_{b,w,e,j,\tilde{h}}\\ \hfill -{\displaystyle \frac{D_{b,w,e,j,\tilde{h}}}{v_{b,w,e}^d}}\end{array}\hfill & & \forall b\in B,w\in W,e=\{el,ht\},j\in J_e,h\in H\hfill \end{array}$

(14)$\begin{array}{cccc}& S_{b,w,e,j,h}\le {\displaystyle \frac{e_{b,w,e}^{RAT}}{|J_e|}}\hfill & & \forall b\in B,w\in W,e=\{el,ht\},j\in J_e,h\in H\hfill \end{array}$

A storage device can only be in one mode at a time—charging, discharging, or steady state—which is enforced by Equation (15). Equations (16) and (17) restrict the charging and discharging rates to their maximum values. Here, parameter $a_{b,w,h}$ defines the presence of batteries, particularly for EVs that may or may not be connected. For stationary storages, $a_{b,w,h}=1$, always. Lastly, Equation (18) determines if the battery is in a steady state:

(15)$\begin{array}{cccc}& B_{b,w,e,h}^C+B_{b,w,e,h}^D+B_{b,w,e,h}^S=1\hfill & & \forall b\in B,w\in W,e=\{el,ht\},h\in H\hfill \end{array}$

(16)$\begin{array}{cccc}& {\displaystyle \sum _{j\in J_e}}{C}_{b,w,e,j,h}\le c_{b,w,e}^{MAX}·a_{b,w,h}·B_{b,w,e,h}^C\hfill & & \forall b\in B,w\in W,e=\{el,ht\},h\in H\hfill \end{array}$

(17)$\begin{array}{cccc}& {\displaystyle \sum _{j\in J_e}}{D}_{b,w,e,j,h}\le d_{b,w,e}^{MAX}·a_{b,w,h}·B_{b,w,e,h}^D\hfill & & \forall b\in B,w\in W,e=\{el,ht\},h\in H\hfill \end{array}$

(18)$\begin{array}{cccc}& 1-{\displaystyle \sum _{j\in J_e}}(D_{b,w,e,j,h}+C_{b,w,e,j,h})\le B_{b,w,e,h}^S\hfill & & \forall b\in B,w\in W,e=\{el,ht\},h\in H\hfill \end{array}$

2.4. Electric Battery Degradation

Variables $K_{b,w,e,h}$ in Equation (19) aggregate costs for cycle-depth degradation, SOC-cycle degradation, and calendar degradation. The component-wise relative capacity losses (whose computations are presented below) are added and multiplied by battery replacement unit cost (r) and the rated battery capacity ($e^{RAT}$). To be explicit, the subscript ‘$el$’ is used in Equations (19)–(37) to denote that only degradation is considered for the electricity battery:

(19)$\begin{array}{cccc}& K_{b,w,el,h}=r·e_{b,w,el}^{RAT}·(Q_{b,w,el,h}^{CD}+Q_{b,w,el,h}^{SOC}+Q_{b,w,el,h}^{CAL})\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

2.4.1. Cycle Depth Degradation

CD aging accounts for the degradation of a battery due to repeated discharging cycles, with the degradation rate increasing non-linearly with deeper discharges. A cycle counting mechanism, such as the rainflow algorithm, is employed to identify and measure these cycles, ensuring consecutive discharges within the same cycle are aggregated appropriately [34]. The CD degradation model used in this study is based on a piecewise linear stress function based on material stress properties as shown in Equation (20) [30,35]. This function captures the relationship between capacity loss per cycle and the depth of discharge, where a is the maximum capacity loss per cycle and m is the fatigue strength exponent:

(20)$\begin{array}{c}\hfill Q^{CD}\left(CD\right)=a^{CD}·C{D}^{\frac{1}{m}}\end{array}$

In this paper, the CD-based degradation after each discharging period is computed using a piecewise linear approximation as defined in Equation (21). In this equation, the first term determines the amount discharged from each battery segment during the time step, normalized by the segment’s rated energy. The second term calculates the relative degradation for the segment using the stress function $Q^{CD}(.)$, effectively summing the degradation contributions across all segments:

(21)$\begin{array}{cccc}& Q_{b,w,el,h}^{CD}={\displaystyle \sum _{j\in J_{el}}}\left({\displaystyle \frac{D_{b,w,el,j,h}}{\frac{e_{b,w,el}^{RAT}}{|J_{el}|}}}\left[Q^{CD}\left({\displaystyle \frac{j}{|J_{el}|}}\right)-Q^{CD}\left({\displaystyle \frac{j-1}{|J_{el}|}}\right)\right]\right)\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

By combining a robust cycle counting mechanism with the piecewise linear degradation model, this approach enables the accurate computation of capacity loss for each segment, providing a detailed and implementable framework for battery aging analysis. For more details, refer to [30].

2.4.2. Average Cycle SOC Degradation

SOC-based aging requires auxiliary values that are used to impose a penalty at the end of each discharge cycle. Equation (22) computes the relative SOC compared to the battery capacity. This model and the corresponding constraints are adapted from [30]:

(22)$\begin{array}{cccc}& SO{C}_{b,w,el,h}={\displaystyle \frac{{\sum }_{j\in J_{el}}{S}_{b,w,el,j,h}}{e_{b,w,el}^{RAT}}}\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

Equations (23)–(25) and (26)–(28) determine the start and end points of specific discharging cycles, respectively:

(23)$\begin{array}{cccc}& B_{b,w,el,h}^{ST}\ge B_{b,w,el,h}^D-B_{b,w,el,\tilde{h}}^D\hfill & & \forall b\in B,w\in W,h,\tilde{h}\in H\hfill \end{array}$

(24)$\begin{array}{cccc}& B_{b,w,el,h}^{ST}\le B_{b,w,el,h}^D\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

(25)$\begin{array}{cccc}& B_{b,w,el,h}^{ST}\le 1-B_{b,w,el,\tilde{h}}^D\hfill & & \forall b\in B,w\in W,h,\tilde{h}\in H\hfill \end{array}$

(26)$\begin{array}{cccc}& B_{b,w,el,h}^{END}\ge B_{b,w,el,\tilde{h}}^D-B_{b,w,el,h}^D\hfill & & \forall b\in B,w\in W,h,\tilde{h}\in H\hfill \end{array}$

(27)$\begin{array}{cccc}& B_{b,w,el,h}^{END}\le B_{b,w,el,\tilde{h}}^D\hfill & & \forall b\in B,w\in W,h,\tilde{h}\in H\hfill \end{array}$

(28)$\begin{array}{cccc}& B_{b,w,el,h}^{END}\le 1-B_{b,w,el,h}^D\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

The binary variables $B_{b,w,el,h}^{ST}$ and $B_{b,w,el,h}^{END}$ represent the start and end of a discharge cycle, respectively. Next, Equations (29)–(32) compute and penalize the average deviation from 50% SOC within each discharge cycle:

(29)$\begin{array}{cccc}& \begin{array}{c}\hfill A_{b,w,el,h}^{ST}=SO{C}_{b,w,el,h}·B_{b,w,el,h}^{ST}\\ \hfill +(1-B_{b,w,el,h}^{ST}-B_{b,w,el,\tilde{h}}^{END})·A_{b,w,el,\tilde{h}}^{ST}\end{array}\hfill & & \forall b\in B,w\in W,h,\tilde{h}\in H\hfill \end{array}$

(30)$\begin{array}{cccc}& A_{b,w,el,h}^{CYC}\ge {\displaystyle \frac{SO{C}_{b,w,el,h}+A_{b,w,el,h}^{ST}}{2}}-{\displaystyle \frac{1}{2}}\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

(31)$\begin{array}{cccc}& A_{b,w,el,h}^{CYC}\ge {\displaystyle \frac{1}{2}}-{\displaystyle \frac{SO{C}_{b,w,el,h}+A_{b,w,el,h}^{ST}}{2}}\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

(32)$\begin{array}{cccc}& Q_{b,w,el,h}^{SOC}=f·A_{b,w,el,h}^{CYC}·B_{b,w,el,h}^{END}\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

$A_{b,w,el,h}^{ST}$ represents the SOC level at the start of a discharge cycle. This is calculated by considering the current SOC and the binary variables indicating the start and end of the cycle. $A_{b,w,el,h}^{CYC}$ measures the average deviation from an (optimal) 50% SOC level during each cycle. Lastly, Equation (32) uses these deviations to calculate the SOC-based degradation cost $Q_{b,w,el,h}^{SOC}$, highlighting the economic implications of SOC fluctuations on battery lifetime.

2.4.3. Calendar Degradation

Calendar aging depends on factors such as time, ambient temperature, and the state of charge, with prolonged periods at high SOC levels accelerating degradation. In this work, the calendar degradation model builds on the approach presented in [30], incorporating an empirically derived degradation rate for lithium-ion batteries as described in [36]:

(33)$\begin{array}{cccc}& B_{b,w,el,i,h}^{CAL}+B_{b,w,el,m,h}^{CAL}\le 1\hfill & & \forall b\in B,w\in W,i,m\in I|m>i+1,h\in H\hfill \end{array}$

(34)$\begin{array}{cccc}& Z_{b,w,el,i,h}\le B_{b,w,el,i,h}^{CAL}\hfill & & \forall b\in B,w\in W,i\in I,h\in H\hfill \end{array}$

(35)$\begin{array}{cccc}& {\displaystyle \sum _{i\in I}}{Z}_{b,w,el,i,h}=1\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

(36)$\begin{array}{cccc}& {\displaystyle \sum _{i\in I}}(o_{b,w,i}^{CAL}·Z_{b,w,el,i,h})=SO{C}_{b,w,el,h}\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

(37)$\begin{array}{cccc}& Q_{b,w,el,h}^{CAL}={\displaystyle \sum _{i\in I}}(y_{b,w,i}·Z_{b,w,el,i,h})\hfill & & \forall b\in B,w\in W,h\in H\hfill \end{array}$

Equation (33) ensures that only one line segment is active at a time. Equation (34) ensures that only the weights $Z_{b,w,el,i,h}$ for the active line segment can be positive. Equation (35) guarantees that the weights sum to 1. Equation (36) determines for an SOC level the weight factors on the relevant line segment. Lastly, the calendar degradation value is calculated using Equation (37).

This concludes the presentation of the model formulation. The next section provides the input data description for the executed numerical experiments.

3. Numerical Results and Discussion

First, we present a set of deterministic cases to provide a clear understanding of the operational dynamics of the system when accounting for battery degradation and a specific charging scheme. Next, we consider uncertainty in grid electricity prices and demand through varying the SOC of EV batteries when they arrive at the building complex to be recharged. All cases are solved using Pyomo with Gurobi version 9.1.2 [37]. The solution time for the deterministic cases is between one and two minutes each on a regular desktop computer with the following specifications: Intel(R) Core(TM) i5-6200U CPU @ 2.30 GHz 2.40 GHz Processor with 8 GB RAM.

3.1. Input Data

The building complex is modeled based on an existing complex in Pfreimd, Germany. Apart from the local grid between the buildings, the buildings are connected to the main grids for electricity, gas, and district heating. We introduce a stationary electric battery and EV charging stations to the setup. We consider an hourly time resolution. For the analysis, we select a day late in spring, Monday, 3 June 2019, when the complex has a peak electricity demand of 68 kW and a peak heat demand of 43 kW. The PV peak generation is 73 kW (for most other parts of the year, the solar panels generate less electricity and the heat demand is high; in such cases, the CHP runs the entire day, which makes the operational dynamics uninteresting, while, in contrast, in the hottest summer months, the heat demand is negligible and solar generation high, which results in very low to no CHP operation). Appendix A.1 details the hourly input data of electricity price, electricity and heat demand, and PV generation of the building complex.

Table 1 presents the CHP output and efficiency for four different operation modes, which form the breakpoints of the piecewise linear segments. Table 2 displays the storage data. We include six EV with battery dimensions of a Nissan Leaf and replacement cost r at 150 EUR/kWh. All degradation data are for NMC batteries. The values for CD parameters a, m, and f are sourced from [30,38]. The calendar degradation parameters are based on [39] with four segments reported in Table 3.

Hourly electricity prices are German spot prices for Monday, 3 June 2019; we add the network fee ($7.39$ ¢/kWh [40]) and value-added tax ($19\%$). The gas price is constant at $6.145$ ¢/kWh and the heat price is 8 ¢/kWh. The CHP has operational costs of $1.7$ ¢/${\mathrm{kWh}}_{el}$. The PV power can be fed into the main grid at a sales price of $6.59$ ¢/kWh independent of time [41]. Electricity generated by CHP can be fed in at the prevailing spot price.

We consider three EV groups with two cars each, and the respective arrival times of 16:00, 17:00, and 20:00, and the respective SOC levels upon arrival of $50\%$, $60\%$, and $40\%$.

3.2. Deterministic Case

We consider EV charging strategies passive, smart, and V2G.

3.2.1. Passive Charging

Figure 2a displays own electricity generation, battery discharges, and grid imports as positive values, while demand, battery charges, and grid exports are indicated by negative values. At the start of the day, from 0:00 to 3:00, the CHP runs at full capacity, then slightly reduces its output until 5:00 when solar generation begins, and continues to operate at a lower capacity until 9:00. As such, the CHP covers most of the electricity and the entire heat demand during the night, while also charging the TES. By assumption, at both the start and the end of the day, the TES is at exactly $30\%$ of its capacity. Charging by the CHP brings the TES to $64\%$ capacity (see Figure 2b), enabling it to meet the heat demand between 09:00 and 15:00. No CHP generated electricity is sold to the grid at any time.

After 9:00, PV generation continues to rise, and the CHP shuts down. Surplus PV power is used for two purposes. First, some is sold to the main grid. Here we note that due to the fixed feed-in remuneration per kWh, the specific selling hour does not matter. Second, batteries are charged. The feed-in tariff is lower than the lowest hourly purchasing price from the grid, making storage for own use attractive. The stationary battery is charged gradually over the day, with a substantial charge to reach the maximum SOC right before the arrival of the first EVs. With grid prices declining gradually between 7:00 and 16:00, the charging schedule is optimized to minimize both costs and calendar aging by keeping the battery at a low SOC during most of the day.

Figure 2c shows the SOC profiles of all battery groups. By assumption, the start and end SOC levels of the stationary battery are $30\%$. For EVs, an SOC of zero indicates that they are not present. The initial SOC levels of EVs at 0:00 equals their respective levels at 24:00.

During the entire evening, there is a high residential and EV demand. With passive charging, as soon as an EV connects, it starts charging at maximum capacity. Additionally, PV generation declines rapidly after 18:00, leading to a spike in grid imports around 20:00 when electricity prices are at their highest for the day. Therefore, while the EV batteries are charging, the stationary battery discharges during these hours to avoid import costs from the main grid. Car groups EV1 and EV2 charge until they are full at 20:00. They both start the next day at $100\%$, while EV3 starts the day at $90\%$ and charges for one more hour to reach $100\%$. Due to the passive charging strategy, the charging cannot be postponed until later hours with lower prices.

EV demand is about $12\%$ of the total electricity demand; however, the cost increase to meet this additional load is almost $30\%$ (increasing from EUR66.86 to EUR86.45). Unfortunately, the passive charging strategy imposes that EVs are charged to a large extent when electricity prices are highest.

3.2.2. Smart Charging

Here, we allow the model to determine when to charge the EVs, given the same charges arrival SOC levels as above, and also impose a full charge upon departure. All other parameter values stay the same.

Figure 3a shows a significant change in the charging schedule compared to passive charging. The earliest arriving EV1 charges for an hour from 16:00 to 17:00 right after arrival. The other EVs do not charge in the evening and start the day with their arrival SOC from the day before. All EVs charge during the lowest price period between 02:00 and 05:00 (At the main grid or system scale, this can be considered to be filling the valley). Locally, this shift has two benefits. The EVs charge at the lowest electricity price, and as the SOC is lower for much of the night, it also reduces calendar aging significantly. Table 4 shows that the annualized battery degradation is reduced by $29\%$, from $1.17\%$ to $0.83\%$.

3.2.3. V2G Charging

Here, we allow the model to decide when EVs charge but also allow discharge (we impose a minimum SOC at 0:00 and 24:00 to $30\%$ for all batteries). Overall, the hourly electricity balance is very similar to that of smart charging. The battery charging hours are all the same. The major difference is in the evening, when EV1 discharges from 19:00 to 21:00. Thereby, it reduces the grid imports during the period of the day with the highest prices. Additionally, EV3 discharges in the last hour of the day (not visible in the figure) and recharges using cheap PV between 06:00 and 09:00 in the morning. Compared to smart charging, the reduction in total system costs from V2G is very small, while EV battery degradation increases by $23\%$ due to the discharges. The operational dynamics of our case study with deterministic data are logical, as shown in Figure 4, and can help to appreciate the results obtained from the case study with uncertainty in the prices and EVs charge demand in the next section.

3.3. Stochastic Case

We explore the effects of two sources of uncertainty: EV arrival SOCs and electricity price. We consider three moments during the day where uncertainty is revealed: two for prices, and one for charge demand. At each price branching point in the scenario tree, we consider upward and downward deviations of $10\%$ compared to the average value in the medium scenario. We assign a probability of $40\%$ to the medium scenario, whereas both the low- and high-outcome scenarios have a $30\%$ probability. Price uncertainty during the middle part of the day would not have much impact, as the EVs are not connected, and surplus PV generation is sold at a fixed feed-in tariff. We choose to let price branches occur at 05:00 in the morning, and 20:00 in the evening, based on battery activity and price levels in the deterministic cases. The price deviations are additive, resulting in five different price levels after 20:00, varying from $-20\%$ to $+20\%$ compared to average values (cf., Figure 5). The probabilities corresponding to these scenarios are documented in Appendix A.2.

We assume that the branching on demand uncertainty, defined by the EVs arrival SOC, happens at 16:00, despite EV2 arriving one hour later. EV1 and EV2 have expected arrival SOC levels of $50\%$ and $60\%$, respectively, that occur with a probability of $80\%$. With $10\%$ probability, EV1–EV2 can arrive both with 8 kWh higher or lower combined SOC level (the volume uncertainty is therefore $4\times 4=16$ kWh, which is $10\%$ of the total battery capacity). We ignore the SOC uncertainty for EV3, as we have seen in the deterministic cases that its late departure time in the morning allows charging in the morning hours using cheap power, and limit the number of scenarios to reduce the needed time to find model solutions. To further reduce the needed solution time, in the following analysis, we deactivate the functionality to compute $Q^{SOC}$ (Equations (23)–(32)), as this affects the solution time very much, but it is the degradation component with the least impact on the total degradation. We consider three different SOC-level and five price scenarios, which combine the results in 15 scenarios in total. We define scenario identification codes PXN. The first letter represents the first price branch and the second letter is the arrival SOC branch, both with three possibilities: low (L), medium (M), or high (H). The third symbol is the number for the second price branch, where 1 is the lowest and 5 the highest (cf., Figure 5). Thus, “LL1” is the scenario: low first price, low arrival SOC, and low second price. Also, “MM3” is the expected value scenario. The scenario tree is illustrated in Figure 6, providing an overview of the 15 scenarios covering a 24 h period. To reduce the length of the analysis, we illustrate the system behavior by contrasting the results of the five selected scenarios only. We emphasize that all 15 scenarios are solved simultaneously at once, which means that decisions early in the day are made anticipating all possible scenarios later in the day.

Again, we evaluate the three different charging strategies but now considering uncertainty in prices and charge demand. The solution time for each stochastic case is several hours.

3.3.1. Stochastic—Passive Charging

EVs have to charge at maximum power as soon as they arrive, until they are fully charged. Consequently, different charging patterns constitute a shorter or longer charge period and are due only to different arrival SOCs (these results are not shown in a figure). Not only the EVs battery charge but also the stationary battery charge and discharge vary dependent on the scenario. Firstly, the first price branch at 05:00 leads to three SOC paths. In all scenarios when prices between 5:00 and 20:00 are lower than average, “LXN”, the stationary battery reaches a maximum SOC of $70\%$ during the day, $80\%$ in all “MXN” scenarios, and $90\%$ in all “HXN” scenarios. The effect of uncertainty is that in most cases, the battery charges significantly more than in the deterministic case. In all scenarios, the stationary battery discharges start at 19:00, when the period starts with the highest electricity prices.

3.3.2. Stochastic—Smart Charging

Figure 7 shows the stationary battery SOC levels for smart charging for the hours between 13:00 and 23:00, with the three maximum levels being the same as for passive charging, with the same time of day to start the stationary battery discharges, at 19:00. EV2 and EV3 charge in the morning, while EV1 charges at the same amount of PV power from 16:00 to 17:00 in all scenarios (see Figure 8). Between 19:00 and 21:00, the stationary battery discharges (until its minimum level of $30\%$) to supply residential demand. Thus, for smart charging, the arrival SOC and the second price level do not affect battery operation compared to passive charging.

3.3.3. Stochastic—Vehicle to Grid

Under V2G charging, the stationary battery profile closely resembles the smart charging profile but with some distinctions arising from EV discharges. As depicted in Figure 9, which shows EV1 SOC profiles between 16:00 and 23:00, in scenarios with above average prices, EV1 begins discharging at 19:00, while in scenarios with continuing low prices in the evening, no discharges occur. The stored electricity is saved for the next day. EV2 engage exclusively in discharging in the “HH5” scenario with consistently high prices, as illustrated in Figure 10.

We observe that the discharge schedules prioritize the utilization of the stationary battery and are driven more by price differences than by battery degradation. The EVs discharge can be beneficial if they arrive early enough to charge with cheap PV power. However, in most situations, PV, CHP, and the stationary battery are cheaper supply options to meet the demand. However, if EV discharges are used, they reduce the degradation of the stationary battery.

3.3.4. Cost Comparison

Figure 11 presents the system cost for the three charging schemes. It shows that CHP costs are the same in all scenarios and form the main part of the total cost. In contrast, grid exchange costs vary significantly. Naturally, net import costs from the main grid are lowest in the low-price scenario “LL1” and highest in the high-price scenario “HH5”. Especially in scenarios with high electricity prices early on (“HXN”), battery degradation costs are relatively high due to high battery discharges to mitigate some of the high price consequences. Battery degradation is up to $5\%$ of the total costs in some scenarios, which is substantial, and highlights the importance of trading off degradation costs and the grid electricity price variations. Considering that CHP costs are such a large component, the relative share of battery degradation cost may be even higher in a setup relying completely on RES. Interestingly, in the middle ground scenario “MM3”, the V2G charging scheme reduces expected overall costs by $2.1\%$ compared to smart charging, much more than in the deterministic case. The added flexibility has higher value when uncertainty is considered. Furthermore, due to the rather modest price variations over the day, the stationary battery cycles only once, which limits the impact of degradation and the interplay with EV charging and discharging. For all charging schemes, calendar degradation is the primary EV battery degradation factor in this analysis, while for the stationary battery, CD degradation is the biggest contributor. Higher price volatility and uncertainty would possibly lead to higher degradation costs, while at the same time reducing total energy costs.

In the middle ground scenario “MM3”, the V2G scheme reduces overall costs by $2.1\%$ compared to smart charging, much more than in the deterministic case.

4. Conclusions

We have developed a novel optimization model that can easily be parameterized to analyze and optimize smart-energy building complexes and neighborhoods. We have examined the impact of different charging strategies for EV batteries under demand and price uncertainties while considering a detailed model for battery degradation. This analysis focuses on how charging strategies affect the overall operational energy cost within a smart urban neighborhood. Our results show an about $10\%$ cost decrease from smart charging compared to passive charging. This reduction is mainly achieved by shifting the charging time frame from the expensive evening period to the early morning hours with lower prices. This delayed charging also reduces EV battery degradation, with reductions of up to $29\%$. When considering V2G and allowing EV discharges, total system costs hardly reduce further compared to smart charging, as the modest direct energy cost decrease is partly undone by higher degradation costs. However, uncertainty in electricity prices and EV arrival SOCs allows a more significant cost improvement of V2G compared to smart charging. Its increased flexibility makes V2G $2.1\%$ cheaper than smart charging.

Battery degradation is up to $5\%$ of total costs in some scenarios, which is substantial, and highlights the importance of accurately trading off degradation costs and grid electricity price variations. Accounting for degradation affects the timing of charges and discharges, not only considering price variations, but weighing them against degradation costs, thereby inducing a slight preference for charging later and (dis) charging less aggressively.

The smart management of a building complex reduces demand in high-price periods and increases demand in lower-price periods. This provides value to the energy system and society at large, increasing the utilization of cheaper generation options, and allowing the higher-value utilization of generated electricity by other sources.

The flexibility of V2G is likely to become increasingly beneficial in a future with much more renewable electricity supply and potentially more volatile prices. As suggested by [42], benefiting from demand shifting should not rely on manual actions and be automated as much as possible. Artificial intelligence will play a key-role in customizing charging patterns for specific users and situations through the real-time monitoring of battery degradation [43,44] Apps and bidirectional charging equipment should empower consumers to easily or automatically set, accept or override suggested charge–discharge patterns, resulting in lower energy bills and contributing to higher system efficiency.

Future work may consider the degradation of other distributed energy resources such as a CHP and TES. Moreover, employing decomposition methods could shorten solution times and enable the examination of extended scheduling periods.

Author Contributions

Conceptualization, C.S., R.E.-B. and P.A.; methodology, C.P., C.S., R.E.-B. and P.A.; software, C.S.; validation, C.S., R.E.-B. and P.A.; formal analysis, C.S. and P.A.; investigation, C.S.; data curation, C.P. and C.S.; writing—original draft preparation, C.S.; writing—review and editing, R.E.-B. and P.A.; visualization, C.S. and P.A.; supervision, R.E.-B. and P.A. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The data presented in this study is available in the paper itself and the data document Schade and Egging-Bratseth (2023).

Conflicts of Interest

The authors declare no conflicts of interest. The funder had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:

CHP	Combined Heat and Power
CD	Cycle Depth
EV	Electric Vehicle
MG	Micro Grid
PV	Photovoltaic
RES	Renewable Energy Sources
SOC	State of Charge
TES	Thermal Energy Storage
V2G	Vehicle-to-Grid
NMC	Nickel Manganese Cobalt

Nomenclature

Index set	Description
$b\in B$	Buildings index
$e\in E$	Energy carriers electricity, heat, and gas: $e=\{\mathrm{el},\mathrm{ht},\mathrm{gs}\}$
$h\in H$	Stochastic scenario tree nodes
$\tilde{h}\in H$	Predecessor scenario tree node of h
$\widehat{h}\in H$	Terminal nodes of scenarios
$h^{ARR}\in H$	EV arrival time scenario node index
$h^{DEP}\in H$	EV departure time scenario node index
$i,m\in I$	Calendar aging breakpoints
$j\in J_e$	(Virtual) segments in storages of type e
$n\in N$	Grid types (main and local grids): $n=\{\mathrm{M},\mathrm{L}\}$
$t\in T$	Renewable generation technologies
$u\in U$	CHP line segment breakpoints
$w\in W$	Storage types, stationary or EVs: $w=\{\mathrm{ST},\mathrm{EV}1,\mathrm{EV}2,\mathrm{EV}3\}$
Variable	Description
$A_{b,w,e,h}^{ST}$	SOC level at the start of the charging cycle.
$A_{b,w,e,h}^{CYC}$	Deviation to 50% of SOC.
$B_{b,w,e,h}^C$	Binary, battery charge mode.
$B_{b,w,e,i,h}^{CAL}$	Binary, SOS2 battery calendar aging.
$B_{u,h}^{CHP}$	Binary, SOS2 CHP operation.
$B_{b,w,e,h}^D$	Binary, battery discharge mode.
$B_{b,w,e,h}^{END}$	Binary, end of a discharge cycle.
$B_{b,w,e,h}^S$	Binary, battery steady-state.
$B_{b,w,e,h}^{ST}$	Binary, start of a discharge cycle.
$C_{b,w,e,j,h}$	Battery charge.
$CH{P}_{b,e,h}^{IN}$	CHP energy input (where $e=gs$).
$CH{P}_{b,e,h}^{OUT}$	CHP energy output (where $e=\mathrm{el},\mathrm{ht}$).
$C{U}_{b,e,h}$	Generation curtailment.
$D_{b,w,e,j,h}$	Battery discharge.
$G_{b,e,n,h}$	Grid extractions.
$I_{b,e,n,h}$	Grid injections.
$I_{b,e,M,h}^{CHP}$	Main grid injections, generated by CHP.
$I_{b,e,M,h}^{RES}$	Main grid injections, generated by RES.
$K_{b,w,e,h}$	Collected cost of electric battery degradation (where $e=\mathrm{el}$).
$Q_{b,w,e,h}^{CAL}$	Calendar aging-based battery capacity loss (where $e=\mathrm{el}$) [%].
$Q_{b,w,e,h}^{CD}$	Cycle depth-based battery capacity loss (where $e=\mathrm{el}$) [%].
$Q_{b,w,e,h}^{SOC}$	Average cycle SOC-based battery capacity loss (where $e=\mathrm{el}$) [%].
$S_{b,w,e,j,h}$	Storage level.
$SO{C}_{b,w,e,h}$	State of charge [%] (where $e=\mathrm{el},\mathrm{ht}$).
$X_{u,h}$	SOS2 segment breakpoint weight for CHP operation.
$Z_{b,w,e,i,h}$	SOS2 segment breakpoint weight for battery aging (where $e=\mathrm{el}$).
Parameter	Description
$a_{b,w,h}$	EV connected to charger.
${\alpha }_h$	Scenario probability.
$c_{b,w,e}^{MAX}$	Maximum charge rate (where $e=\mathrm{el},\mathrm{ht}$).
$d_{b,w,e}^{MAX}$	Maximum discharge rate (where $e=\mathrm{el},\mathrm{ht}$).
$de{m}_{b,e,h}$	Residential energy demand (where $e=\mathrm{el},\mathrm{ht}$).
$e_{b,w,e}^{RAT}$	Rated energy of the storage (where $e=\mathrm{el},\mathrm{ht}$).
f	Fitting parameter for average SOC stress of storage.
$g_{b,t,e,h}^{RES}$	Renewable energy production (where $t=\mathrm{PV},e=\mathrm{el}$).
$k^{CHP}$	Operational costs of CHP.
$k^{EOH}$	Deviation penalty of end-of-horizon SOC-level’s storage.
$p_{e,h}$	Grid import price (where $e=\mathrm{el},\mathrm{ht}$) [EUR/kWh].
$p^{CHP}$	Feed-in tariff CHP power [EUR/kWh].
$p_{gs}$	Gas price [EUR/kWh].
$p^{RES}$	Feed-in tariff RES power [EUR/kWh].
r	Replacement cost of the battery [EUR/kWh].
$o_{b,w,e,i,h}^{CAL}$	Calendar aging breakpoint (where $e=\mathrm{el}$).
$o_{u,gs}^{CHP}$	Breakpoint value CHP gas intake.
$o_{u,e}^{CHP}$	Breakpoint value CHP energy generation (where $e=\mathrm{el},\mathrm{ht}$).
$v_{b,w,e}^c$	Storage charging efficiency (where $e=\mathrm{el},\mathrm{ht}$).
$v_{b,w,e}^d$	Storage discharging efficiency (where $e=\mathrm{el},\mathrm{ht}$).
$v_{e,n}$	Grid injection efficiency (where $e=\mathrm{el},\mathrm{ht},n=\mathrm{M},\mathrm{L}$).
$y_{b,w,e,i}$	Breakpoint value calendar aging (where $e=\mathrm{el}$).

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Scenario tree example illustration.

Figure 2. Passive charging results—deterministic. (a) Hourly mass balance electricity. (b) Hourly mass balance heat. (c) Battery and EV SOC profiles.

Figure 3. Smart charging results—deterministic. (a) Hourly mass balance electricity. (b) Battery and EV SOC profiles.

Figure 4. V2G charging results—deterministic. (a) Hourly mass balance electricity. (b) Hourly mass balance heat. (c) Battery and EV SOC profiles.

Figure 5. Price scenarios in the stochastic analysis. The colored dots represent the price paths in different scenarios.

Figure 6. Scenario tree showing 15 scenarios based on price and EVs’ SOC uncertainties over 24 h.

Figure 7. SOC profiles, stationary battery, smart charging, and stochastic. Note that line LXN shows the profile for both LL1 and LH3. The same for line HXN with regard to HL3 and HH5.

Figure 8. SOC profiles, EV1, smart charging, and stochastic.Note that the EV is disconnected from 6:00–16:00, after which it can arrive with one of three different SOC-levels. Additionally, lines LL1 and HL3 overlap, as do lines LH3 and HH5.

Figure 9. SOC profiles, EV1, V2G charging—stochastic.

Figure 10. EV2, V2G charging—stochastic.

Figure 11. System cost breakdown. Each column group: left, passive; middle, smart; right, V2G. Net grid import costs main grid (-feedin) is the difference between power purchase cost and remuneration for solar power feed-in. Vertical axis truncated for readability.

Table 1

CHP segment breakpoint data.

Operating Level	0%	50%	75%	100%
Gas consumption [kW]	0	37.1	48.1	61.1
Electric output [kW]	0	10.0	15.0	20.0
Heat output [kW]	0	28.1	34.2	42.2

The technical CHP specifications document is available at https://www.ecpower.eu/files/ec-power/customer/EN/Downloads_EN/Information_Material/EC_POWER_EN_technical_data_XRGI20.pdf, accessed on 29 April 2022.

Table 2

Storage characteristics.

Storage Type	Battery	TES	EV (Each)
Capacity [kWh]	150	172	40
Power [kW]	90	42.2	7.0
Round trip efficiency [%]	90	90	88
Assumed operating temperature [°C]	25	–	25
Max. capacity loss [%/cycle]	0.04519	–	0.04519
Fatigue strength exponent	0.4926	–	0.4926
Fitting parameter SOC stress [%/cycle]	0.0085	–	0.0085

Table 3

Battery degradation segment breakpoint data $\left({10}^{-7}\right)$.

Breakpoint oCAL	0%	30%	60%	70%	100%
Aging val OLD	15.01	35.02	40.02	73.62	89.34
Aging value y	3.75	8.76	10.01	18.41	22.34

Table 4

Deterministic case KPI.

KPI	Passive	Smart	V2G
Objective value [EUR]	86.45	80.16	80.13
Battery throughput [kWh]	64.12	64.12	64.12
PV generation [kWh]	681.13	681.13	681.13
Electric CHP generation [kWh]	241.31	242.50	242.50
Peak demand [kW]	80.46	68.08	68.08
Peak grid import [kW]	67.64	53.81	53.81
Total grid import [kWh]	294.01	273.42	268.82
Total grid export [kWh]	161.08	136.34	136.34
Annualized stationary battery aging [%]	4.12	4.11	4.11
EV battery aging while connected [%] 1	1.17	0.83	1.02

¹ Average annualized EV battery aging while EVs are connected. We ignore degradation from discharges while driving.

Appendix A. Hourly Data and Scenario Specifics

This appendix includes the hourly data and scenario specifics.

Appendix A.1. Hourly Data and PV Generation

Table A1 lists the hourly demand and PV generation data for the building complex and the spot price of electricity for 3 June 2019.

Table A1

Hourly data for demand, PV generation, and electricity prices (including network charges and VAT).

Time	Electricity Price [EUR/kWh]	Total Electricity Demand [kW]	Total Heat Demand [kW]	PV Generation [kW]
0	0.120	26.479	23.245	00.000
1	0.116	24.282	30.809	00.000
2	0.113	23.337	39.999	00.000
3	0.112	19.982	41.164	00.000
4	0.111	19.140	43.635	00.000
5	0.117	20.185	29.156	00.001
6	0.136	26.059	27.257	08.926
7	0.141	26.390	23.605	25.874
8	0.142	32.282	24.549	40.363
9	0.138	33.548	19.357	52.565
10	0.135	32.011	16.052	62.179
11	0.133	35.344	13.691	68.833
12	0.131	37.014	14.635	72.292
13	0.130	40.673	15.108	72.585
14	0.130	36.176	14.163	69.800
15	0.130	44.011	12.275	63.963
16	0.130	49.439	11.803	55.149
17	0.139	49.209	12.747	43.695
18	0.142	46.269	17.468	30.095
19	0.148	66.893	19.356	14.801
20	0.146	68.080	20.300	00.003
21	0.145	66.630	18.412	00.000
22	0.143	53.811	15.579	00.000
23	0.137	40.221	10.859	00.000

Appendix A.2. Scenario Specifics

Table A2 provides information regarding the start and end nodes of branches for the multi-stage scenario programming model. Table A3 reports the probabilities of scenarios.

Table A2

First and last scenario tree nodes in each stage.

Stage (Hours)	First Node	Last Node
0	0	0
1	1	1
2	2	2
3	3	3
4	4	4
5	5	7
6	8	10
7	11	13
8	14	16
9	17	19
10	20	22
11	23	25
12	26	28
13	29	31
14	32	34
15	35	37
16	38	46
17	47	55
18	56	64
19	65	73
20	74	100
21	101	127
22	128	154
23	155	181

Table A3

Probability of scenarios.

P	X	N
		0.009
	0.03	0.012
		0.009
		0.072
0.3	0.24	0.096
		0.072
		0.009
	0.03	0.012
		0.009
		0.012
	0.04	0.016
		0.012
		0.096
0.4	0.32	0.128
		0.096
		0.012
	0.04	0.016
		0.012
		0.009
	0.03	0.012
		0.009
		0.072
0.3	0.24	0.096
		0.072
		0.009
	0.03	0.012
		0.009

References

1. Vadari, S. Smart Grid Redefined: Transformation of the Electric Utility; Artech House: Norwood, MA, USA, 2018.

2. Gabbar, H. Smart Energy Grid Engineering; Academic Press: Cambridge, MA, USA, 2016.

3. Paterakis, N.G.; Erdinç, O.; Pappi, I.N.; Bakirtzis, A.G.; Catalão, J.P. Coordinated operation of a neighborhood of smart households comprising electric vehicles, energy storage and distributed generation. IEEE Trans. Smart Grid; 2016; 7, pp. 2736-2747. [DOI: https://dx.doi.org/10.1109/TSG.2015.2512501]

4. Aliabadi, F.E.; Agbossou, K.; Kelouwani, S.; Henao, N.; Hosseini, S.S. Coordination of smart home energy management systems in neighborhood areas: A systematic review. IEEE Access; 2021; 9, pp. 36417-36443. [DOI: https://dx.doi.org/10.1109/ACCESS.2021.3061995]

5. Lund, P.D.; Mikkola, J.; Ypyä, J. Smart energy system design for large clean power schemes in urban areas. J. Clean. Prod.; 2015; 103, pp. 437-445. [DOI: https://dx.doi.org/10.1016/j.jclepro.2014.06.005]

6. Pulazza, G.; Orozco, C.; Borghetti, A.; Tossani, F.; Napolitano, F. Procurement cost minimization of an energy community with biogas, photovoltaic and storage units. Proceedings of the 2021 IEEE Madrid PowerTech; Madrid, Spain, 28 June–2 July 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1-6. [DOI: https://dx.doi.org/10.1109/PowerTech46648.2021.9494878]

7. Amara-Ouali, Y.; Goude, Y.; Massart, P.; Poggi, J.M.; Yan, H. A review of electric vehicle load open data and models. Energies; 2021; 14, 2233. [DOI: https://dx.doi.org/10.3390/en14082233]

8. Leippi, A.; Fleschutz, M.; Murphy, M.D. A review of ev battery utilization in demand response considering battery degradation in non-residential vehicle-to-grid scenarios. Energies; 2022; 15, 3227. [DOI: https://dx.doi.org/10.3390/en15093227]

9. Ahmadian, A.; Sedghi, M.; Elkamel, A.; Fowler, M.; Golkar, M.A. Plug-in electric vehicle batteries degradation modeling for smart grid studies: Review, assessment and conceptual framework. Renew. Sustain. Energy Rev.; 2018; 81, pp. 2609-2624. [DOI: https://dx.doi.org/10.1016/j.rser.2017.06.067]

10. Yusuf, J.; Ula, S. A comprehensive optimization solution for buildings with distributed energy resources and V2G operation in smart grid applications. Proceedings of the 2020 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT); Washington, DC, USA, 17–20 February 2020; pp. 1-5. [DOI: https://dx.doi.org/10.1109/ISGT45199.2020.9087656]

11. Shakeri, M.; Shayestegan, M.; Abunima, H.; Reza, S.S.; Akhtaruzzaman, M.; Alamoud, A.; Sopian, K.; Amin, N. An intelligent system architecture in home energy management systems (HEMS) for efficient demand response in smart grid. Energy Build.; 2017; 138, pp. 154-164. [DOI: https://dx.doi.org/10.1016/j.enbuild.2016.12.026]

12. Ahmadian, A.; Sedghi, M.; Mohammadi-ivatloo, B.; Elkamel, A.; Golkar, M.A.; Fowler, M. Cost-benefit analysis of V2G implementation in distribution networks considering PEVs battery degradation. IEEE Trans. Sustain. Energy; 2017; 9, pp. 961-970. [DOI: https://dx.doi.org/10.1109/TSTE.2017.2768437]

13. Zhou, S.; Han, Y.; Mahmoud, K.; Darwish, M.M.; Lehtonen, M.; Yang, P.; Zalhaf, A.S. A novel unified planning model for distributed generation and electric vehicle charging station considering multi-uncertainties and battery degradation. Appl. Energy; 2023; 348, 121566. [DOI: https://dx.doi.org/10.1016/j.apenergy.2023.121566]

14. Ghasemi, A.; Banejad, M.; Rahimiyan, M.; Zarif, M. Investigation of the micro energy grid operation under energy price uncertainty with inclusion of electric vehicles. Sustain. Oper. Comput.; 2021; 2, pp. 12-19. [DOI: https://dx.doi.org/10.1016/j.susoc.2021.02.002]

15. Li, H.; Rezvani, A.; Hu, J.; Ohshima, K. Optimal day-ahead scheduling of microgrid with hybrid electric vehicles using MSFLA algorithm considering control strategies. Sustain. Cities Soc.; 2021; 66, 102681. [DOI: https://dx.doi.org/10.1016/j.scs.2020.102681]

16. Yan, Z.; Duan, X.; Chang, Y.; Xu, Z.; Sobhani, B. Optimal energy management in smart buildings with electric vehicles based on economic and risk aspects using developed whale optimization algorithm. J. Clean. Prod.; 2023; 415, 137710. [DOI: https://dx.doi.org/10.1016/j.jclepro.2023.137710]

17. Mansouri, S.A.; Nematbakhsh, E.; Jordehi, A.R.; Marzband, M.; Tostado-Véliz, M.; Jurado, F. An interval-based nested optimization framework for deriving flexibility from smart buildings and electric vehicle fleets in the TSO-DSO coordination. Appl. Energy; 2023; 341, 121062. [DOI: https://dx.doi.org/10.1016/j.apenergy.2023.121062]

18. Eghbali, N.; Hakimi, S.M.; Hasankhani, A.; Derakhshan, G.; Abdi, B. A scenario-based stochastic model for day-ahead energy management of a multi-carrier microgrid considering uncertainty of electric vehicles. J. Energy Storage; 2022; 52, 104843. [DOI: https://dx.doi.org/10.1016/j.est.2022.104843]

19. Alipour, M.; Mohammadi-Ivatloo, B.; Zare, K. Stochastic scheduling of renewable and CHP-based microgrids. IEEE Trans. Ind. Inform.; 2015; 11, pp. 1049-1058. [DOI: https://dx.doi.org/10.1109/TII.2015.2462296]

20. Esfahani, M.; Alizadeh, A.; Cao, B.; Kamwa, I.; Xu, M. A stochastic-robust aggregation strategy for VPP to participate in the frequency regulation market via backup batteries. J. Energy Storage; 2024; 98, 113057. [DOI: https://dx.doi.org/10.1016/j.est.2024.113057]

21. Rafique, S.; Hossain, M.J.; Nizami, M.S.H.; Irshad, U.B.; Mukhopadhyay, S.C. Energy management systems for residential buildings with electric vehicles and distributed energy resources. IEEE Access; 2021; 9, pp. 46997-47007. [DOI: https://dx.doi.org/10.1109/ACCESS.2021.3067950]

22. Nunna, H.K.; Battula, S.; Doolla, S.; Srinivasan, D. Energy management in smart distribution systems with vehicle-to-grid integrated microgrids. IEEE Trans. Smart Grid; 2016; 9, pp. 4004-4016. [DOI: https://dx.doi.org/10.1109/TSG.2016.2646779]

23. Al Muala, Z.A.; Bany Issa, M.A.; Bello Bugallo, P.M. Integrating Life Cycle Principles in Home Energy Management Systems: Optimal Load PV–Battery–Electric Vehicle Scheduling. Batteries; 2024; 10, 138. [DOI: https://dx.doi.org/10.3390/batteries10040138]

24. Guo, Q.; Nojavan, S.; Lei, S.; Liang, X. Economic-environmental analysis of renewable-based microgrid under a CVaR-based two-stage stochastic model with efficient integration of plug-in electric vehicle and demand response. Sustain. Cities Soc.; 2021; 75, 103276. [DOI: https://dx.doi.org/10.1016/j.scs.2021.103276]

25. Ding, T.; Qu, M.; Huang, C.; Wang, Z.; Du, P.; Shahidehpour, M. Multi-period active distribution network planning using multi-stage stochastic programming and nested decomposition by SDDIP. IEEE Trans. Power Syst.; 2020; 36, pp. 2281-2292. [DOI: https://dx.doi.org/10.1109/TPWRS.2020.3032830]

26. Hafiz, F.; de Queiroz, A.R.; Fajri, P.; Husain, I. Energy management and optimal storage sizing for a shared community: A multi-stage stochastic programming approach. Appl. Energy; 2019; 236, pp. 42-54. [DOI: https://dx.doi.org/10.1016/j.apenergy.2018.11.080]

27. Chai, Y.T.; Che, H.S.; Tan, C.; Tan, W.N.; Yip, S.C.; Gan, M.T. A two-stage optimization method for Vehicle to Grid coordination considering building and Electric Vehicle user expectations. Int. J. Electr. Power Energy Syst.; 2023; 148, 108984. [DOI: https://dx.doi.org/10.1016/j.ijepes.2023.108984]

28. Liu, W.H.; Alwi, S.R.W.; Hashim, H.; Muis, Z.A.; Klemeš, J.J.; Rozali, N.E.M.; Lim, J.S.; Ho, W.S. Optimal design and sizing of integrated centralized and decentralized energy systems. Energy Procedia; 2017; 105, pp. 3733-3740. [DOI: https://dx.doi.org/10.1016/j.egypro.2017.03.866]

29. Lu, R.; Ding, T.; Qin, B.; Ma, J.; Fang, X.; Dong, Z. Multi-stage stochastic programming to joint economic dispatch for energy and reserve with uncertain renewable energy. IEEE Trans. Sustain. Energy; 2019; 11, pp. 1140-1151. [DOI: https://dx.doi.org/10.1109/TSTE.2019.2918269]

30. Schade, C.; Egging-Bratseth, R. Battery degradation: Impact on economic dispatch. Energy Storage; 2024; 6, e588. [DOI: https://dx.doi.org/10.1002/est2.588]

31. White, C.; Thompson, B.; Swan, L.G. Comparative performance study of electric vehicle batteries repurposed for electricity grid energy arbitrage. Appl. Energy; 2021; 288, 116637. [DOI: https://dx.doi.org/10.1016/j.apenergy.2021.116637]

32. Gandoman, F.H.; Jaguemont, J.; Goutam, S.; Gopalakrishnan, R.; Firouz, Y.; Kalogiannis, T.; Omar, N.; Van Mierlo, J. Concept of reliability and safety assessment of lithium-ion batteries in electric vehicles: Basics, progress, and challenges. Appl. Energy; 2019; 251, 113343. [DOI: https://dx.doi.org/10.1016/j.apenergy.2019.113343]

33. Schade, C.; Egging-Bratseth, R. Battery Degradation: Impact on Economic Dispatch—Data documentation; Zenodo: Geneva, Switzerland, 2023; [DOI: https://dx.doi.org/10.5281/zenodo.10255081]

34. Loew, S.; Anand, A.; Szabo, A. Economic model predictive control of Li-ion battery cyclic aging via online rainflow-analysis. Energy Storage; 2021; 3, e228. [DOI: https://dx.doi.org/10.1002/est2.228]

35. Xu, B.; Zhao, J.; Zheng, T.; Litvinov, E.; Kirschen, D.S. Factoring the cycle aging cost of batteries participating in electricity markets. IEEE Trans. Power Syst.; 2017; 33, pp. 2248-2259. [DOI: https://dx.doi.org/10.1109/TPWRS.2017.2733339]

36. de Hoog, J.; Timmermans, J.M.; Ioan-Stroe, D.; Swierczynski, M.; Jaguemont, J.; Goutam, S.; Omar, N.; Van Mierlo, J.; Van Den Bossche, P. Combined cycling and calendar capacity fade modeling of a Nickel-Manganese-Cobalt Oxide Cell with real-life profile validation. Appl. Energy; 2017; 200, pp. 47-61. [DOI: https://dx.doi.org/10.1016/j.apenergy.2017.05.018]

37. Gurobi Optimization. Available online: www.gurobi.com/products/gurobi-optimizer (accessed on 29 November 2024).

38. Laresgoiti, I.; Käbitz, S.; Ecker, M.; Sauer, D.U. Modeling mechanical degradation in lithium ion batteries during cycling: Solid electrolyte interphase fracture. J. Power Sources; 2015; 300, pp. 112-122. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2015.09.033]

39. Keil, P.; Schuster, S.F.; Wilhelm, J.; Travi, J.; Hauser, A.; Karl, R.C.; Jossen, A. Calendar aging of lithium-ion batteries. J. Electrochem. Soc.; 2016; 163, A1872. [DOI: https://dx.doi.org/10.1149/2.0411609jes]

40. Electricity Prices by Type of User. Eurostat. 2019; Available online: https://ec.europa.eu/eurostat/databrowser/view/nrg_pc_204/default/table?lang=en (accessed on 29 November 2024).

41. Renewable Energy Sources Act—EEG. BMWI. 2019; Available online: www.gesetze-im-internet.de/eeg_2014/__21.html (accessed on 29 November 2024).

42. Hofmann, M.; Lindberg, K.B. Residential demand response and dynamic electricity contracts with hourly prices: A study of Norwegian households during the 2021/22 energy crisis. Smart Energy; 2024; 13, 100126. [DOI: https://dx.doi.org/10.1016/j.segy.2023.100126]

43. Teodorescu, R.; Sui, X.; Vilsen, S.B.; Bharadwaj, P.; Kulkarni, A.; Stroe, D.I. Smart battery technology for lifetime improvement. Batteries; 2022; 8, 169. [DOI: https://dx.doi.org/10.3390/batteries8100169]

44. de la Vega, J.; Riba, J.R.; Ortega-Redondo, J.A. Real-Time Lithium Battery Aging Prediction Based on Capacity Estimation and Deep Learning Methods. Batteries; 2023; 10, 10. [DOI: https://dx.doi.org/10.3390/batteries10010010]