Coupling EnergyPlus with the Building Controls Virtual Test Bed[LINK]
Algorithm for data exchange[LINK]
The process in which at least two simulators solve initialvalue differential equations that are coupled to each other is called cosimulation. Various algorithms are possible for the data exchange. In the BCVTB, data are exchanged between its client programs, including EnergyPlus, using a fixed synchronization time step. There is no iteration between the clients. In the cosimulation literature, this coupling scheme is referred to as quasidynamic coupling, loose coupling or pingpong coupling (Hensen 1999, Zhai and Chen 2005).
The algorithm for exchanging data is as follows: Suppose we have a system with two clients, with client 1 being EnergyPlus and client 2 being, for example, the Simulink program from Mathworks. Suppose each client solves an initialvalue ordinary differential equation that is coupled to the differential equation of the other client. Let N∈N denote the number of time steps and let k∈{1,...,N} denote the time steps. We will use the subscripts 1 and 2 to denote the state variable and the function that computes the next state variable of the simulator 1 and 2, respectively.
The simulator 1 computes, for k∈{1,...,N−1} the sequence
x1(k+1) = f1(x1(k), x2(k))
and, similarly, the simulator 2 computes the sequence
x2(k+1) = f2(x2(k), x1(k))
with initial conditions x1(0) = x1,0 and x2(0) = x2,0.
To advance from time k to k+1, each simulator uses its own time integration algorithm. At the end of the time step, the simulator 1 sends the new state x1(k+1) to the BCVTB and it receives the state x2(k+1) from the BCVTB. The same procedure is done with the simulator 2. The BCVTB synchronizes the data in such a way that it does not matter which of the two simulators is called first.
In comparison to numerical methods of differential equations, this scheme is identical to an explicit Euler integration, which is an integration algorithm that computes for an ordinary differential equation with specified initial values,
dx/dt = h(x),
x(0) = x0,
on the time interval t ∈ [0, 1], the following sequence:
Step 0: 
Initialize counter k=0 and number of steps N∈N. 

Set initial state x(k) = x and set time step Dt = 1/N. 
Step 1: 
Compute new state x(k+1) = x(k) + h(x(k)) Dt. 

Replace k by k+1. 
Step 2: 
If k=N stop, else go to Step 1. 
In the situation where the differential equation is solved using cosimulation, the above algorithm becomes
Step 0: 
Initialize counter k=0 and number of steps N∈N . 

Set initial state x (k) = x and x (k) = x . Set the time step Dt = 1/N. 
Step 1: 
Compute new states x (k+1) = x (k) + f (x (k), x (k)) Dt, and x (k+1) = x (k) + f (x (k), x (k)) Dt. 

Replace k by k+1. 
Step 2: 
If k=N stop, else go to Step 1. 
This algorithm is implemented in the BCVTB. Note that there is no iteration between the two simulators.
Coupling EnergyPlus with the Building Controls Virtual Test Bed[LINK]
Algorithm for data exchange[LINK]
The process in which at least two simulators solve initialvalue differential equations that are coupled to each other is called cosimulation. Various algorithms are possible for the data exchange. In the BCVTB, data are exchanged between its client programs, including EnergyPlus, using a fixed synchronization time step. There is no iteration between the clients. In the cosimulation literature, this coupling scheme is referred to as quasidynamic coupling, loose coupling or pingpong coupling (Hensen 1999, Zhai and Chen 2005).
The algorithm for exchanging data is as follows: Suppose we have a system with two clients, with client 1 being EnergyPlus and client 2 being, for example, the Simulink program from Mathworks. Suppose each client solves an initialvalue ordinary differential equation that is coupled to the differential equation of the other client. Let N∈N denote the number of time steps and let k∈{1,...,N} denote the time steps. We will use the subscripts 1 and 2 to denote the state variable and the function that computes the next state variable of the simulator 1 and 2, respectively.
The simulator 1 computes, for k∈{1,...,N−1} the sequence
x1(k+1) = f1(x1(k), x2(k))
and, similarly, the simulator 2 computes the sequence
x2(k+1) = f2(x2(k), x1(k))
with initial conditions x1(0) = x1,0 and x2(0) = x2,0.
To advance from time k to k+1, each simulator uses its own time integration algorithm. At the end of the time step, the simulator 1 sends the new state x1(k+1) to the BCVTB and it receives the state x2(k+1) from the BCVTB. The same procedure is done with the simulator 2. The BCVTB synchronizes the data in such a way that it does not matter which of the two simulators is called first.
In comparison to numerical methods of differential equations, this scheme is identical to an explicit Euler integration, which is an integration algorithm that computes for an ordinary differential equation with specified initial values,
dx/dt = h(x),
x(0) = x0,
on the time interval t ∈ [0, 1], the following sequence:
In the situation where the differential equation is solved using cosimulation, the above algorithm becomes
This algorithm is implemented in the BCVTB. Note that there is no iteration between the two simulators.
Documentation content copyright © 19962016 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.