<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Use of Branch and Bound Algorithms for Greenhouse Climate Control</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>George Dimokas</string-name>
          <email>gedimokas@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Constantinos Kittas</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Department of Agriculture Crop Production and Rural Environment, University of Thessaly</institution>
          ,
          <country country="GR">Greece</country>
        </aff>
      </contrib-group>
      <fpage>136</fpage>
      <lpage>145</lpage>
      <abstract>
        <p>Optimization of greenhouse climate management during winter period is an issue that intensely preoccupied researchers over the last two decades as it is directly linked to energy saving, products quality, and reduction of chemical inputs. Goal of this project was the use of Branch and Bound algorithms in order to obtain the optimize greenhouse climate control. For this purpose a biophysical simulator were used and experiments were carried out in the farm of the University of Thessaly in the region of Volos (Greece), during the autumn and winter periods of 2005 and 2007. Branch and Bound algorithms used for two different climate scenarios and the results showed the difference between the classical greenhouse climate control and the control according B &amp; B optimization technique. Finally results showed the contribution of optimization technique to increase tomato production and to reduce energy consumption.</p>
      </abstract>
      <kwd-group>
        <kwd>Branch and Bound Algorithms</kwd>
        <kwd>Climate Control</kwd>
        <kwd>Greenhouse</kwd>
        <kwd>tomato production</kwd>
        <kwd>optimization</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1 Introduction</title>
      <p>
        Most of the greenhouse climate control problems show large number of possible
solutions, which are therefore entered the need for finding the best "route" or optimal
solution
        <xref ref-type="bibr" rid="ref4 ref5">(Dimokas 2009)</xref>
        . The use of advanced optimization models, may contribute
to the variation of the classical greenhouse climate management during winter period
which usually consists in the management of heating by specifying two desired
thermostat setting one temperature, for the night and one for day period, based
mainly on the producer experience
        <xref ref-type="bibr" rid="ref9">(Tap et al. 1993)</xref>
        . The desired values of the
temperature setting (set- points) are depending on the type and the age of crop.
Nevertheless, experimental work showed that the growth and development of many
vegetable and horticultural species appears to respond more to an average daily
temperature than on accurate temperature evolution during the day
        <xref ref-type="bibr" rid="ref12 ref6 ref8">(Heuvelink 1989,
Vogelezang et al. 2000)</xref>
        .
      </p>
      <p>
        The heating of the greenhouse by the sum of the temperatures (Integrated
Temperature Control) during the day has already been applied in several floriculture
        <xref ref-type="bibr" rid="ref12 ref8">(Rijsdijk and Vogelezang 2000)</xref>
        and horticultural species
        <xref ref-type="bibr" rid="ref3">(De Konning 1988)</xref>
        ,
allowing an energy saving of about 10 - 20%
        <xref ref-type="bibr" rid="ref1">(Bailey and Seginer 1989)</xref>
        .
      </p>
      <p>
        Many experimental projects are referred to optimization of greenhouse climate.
The optimal control is one of the processes studied further
        <xref ref-type="bibr" rid="ref10 ref11">(van Henten 1994,
Tchamitchian and Tantau 1996)</xref>
        . According to this method, a model function of the
optimization system is used, consisting of differential first degree equations and an
algebraic criterion that measures the result quality of this operation.
      </p>
      <p>
        The principle of Bellman
        <xref ref-type="bibr" rid="ref2">(Bellman 1957)</xref>
        , allows to solve the problem with
dynamic programming, while the principle of Pontryagin
        <xref ref-type="bibr" rid="ref7">(Pontryagin et al. 1962)</xref>
        uses the Lagrange multipliers to convert the power system model and the algebraic
criterion to one function for minimization. These experimental works have not yet
reached to commercial systems. Goal of this project was the use of Branch and
Bound algorithms in order to obtain the optimize greenhouse climate control for
energy saving.
2
      </p>
      <p>Material &amp; Methods</p>
      <p>Branch and bound (BB or B&amp;B) is an algorithm design paradigm for discrete and
combinatorial optimization problems, as well as general real valued problems. A
branch-and-bound algorithm consists of a systematic enumeration of candidate
solutions by means of state space search: the set of candidate solutions is thought of
as forming a rooted tree (Figure 1) with the full set at the root. The algorithm
explores branches of this tree, which represent subsets of the solution set. Before
enumerating the candidate solutions of a branch, the branch is checked against upper
and lower estimated bounds on the optimal solution, and is discarded if it cannot
produce a better solution than the best one found so far by the algorithm.
S
(Α)</p>
      <p>
        The method of branching and bounding (Branch &amp; Bound) has found application
in solving various and important optimization problems, eg, in integer programming,
nonlinear problems, programming problems, plant sitting problems
        <xref ref-type="bibr" rid="ref4 ref5">(Dimokas 2009)</xref>
        .
S2
S21
S22
      </p>
      <p>S
S1
S3
S1
S3</p>
      <p>S4</p>
      <p>S31
S32</p>
      <p>S1</p>
      <p>S2</p>
      <p>S3</p>
      <p>S4 (Β)
S21</p>
      <p>S22 S31</p>
      <p>S32
(C)
2.1</p>
      <p>Optimization Method for Greenhouse Climate Control</p>
      <p>Optimization of the biophysical simulator with the use of Branch &amp; Bound
method uses a space selection strategy in order to be investigated in accordance with
the algorithms (1), (2) below. The algorithm (1), is responsible for setting the (Vh)
ventilation (system controller) inside the greenhouse
where pb1, pc1, pRΤ, the values of the variables that have to be optimized in order to
give the best solution. At the same time the algorithm (1), uses the biophysical
simulator and more specifically the results obtained for the relative humidity (RH)
and the air temperature inside the greenhouse (Ti), during the process of
optimization.</p>
      <p>The algorithm (2) below is responsible for determining the operation of the (Vt)
heating system (system controller) within the greenhouse.</p>
      <p>Vt = 0.1⋅ (1 + (</p>
      <p>RH − pRT
( Abs ⋅ (RH − pRT )) + 0.5
)) +
e ⋅</p>
      <p>1
( pb1−Ti )</p>
      <p>pc1
e ⋅</p>
      <p>1
( pb2p−c2DTi )
(1)
(2)
(3)
(4)
where pb2, pc2, pRΤ, the values of the variables that have to be optimized in order to
give the best solution. At the same time the algorithm (2), uses the biophysical
simulator and more specifically the results obtained for the relative humidity (RH)
and the air temperature difference inside and outside the greenhouse (DTi), during
the process of optimization.</p>
      <p>The problem to be solved is to minimize the objective function J to a range of
possible solutions, S:</p>
      <p>min J = ∑ (Fcon) − ∑ (D.W .F )
where</p>
      <p>∑ (Fcon) is the total energy gives the heating system inside the greenhouse,
while respectively ∑(D.W.F) the resulting dry weight of mature fruit.</p>
      <p>The limiting function used to optimize the biophysical simulator using Branch &amp;
Bound method shown below:</p>
      <p>min J (one) &gt; max J (other )</p>
      <p>Function (4) reject some subset of possible solutions by further exploring inside,
to find the possible solution when the above condition is true, that the optimal
minimal solution that offers the particular subset is greater than the largest value of
another subset.
2.2 Climate and Biological Measurements</p>
      <p>
        The measurements that used to optimize the model were data for tomato
production and development simultaneously with the greenhouse climatic data,
during the autumn and winter periods of 2005 and 2007. For the experimental
periods were used also calculated values by the modified TOMGRO
        <xref ref-type="bibr" rid="ref4 ref5">(Dimokas et. al.,
2009)</xref>
        , for: i) plant development and the number of leaves, fruits, flowers, ii) biomass
and fruit production. The aim was to identify differences between the optimization
method, and the climate control during the experimental period.
      </p>
    </sec>
    <sec id="sec-2">
      <title>3 Results</title>
      <p>This section presents results of treatment followed during the experimental
measurements and optimum proposed in accordance with the branch &amp; bound
method. The results are giving detailed greenhouse climatic conditions and all
features related to the development and production of tomato plants. From these were
selected and are presented in the following sections the results concerning: a. number
of plants node, b. shoot dry weight, c. leaves dry weight, d. whole plant dry weight,
e. simulation curve of windows opening, f. air temperature inside the greenhouse, g.
temperature of the greenhouse cover. Values that are used as input variables for the
number of nodes, dry weight of leaves, stems and fruits are the same that used of the
modified TOMGRO.
3.1 Results of the first simulated period</p>
      <p>Initially a variation was observed in greenhouse air - cover temperature and
presented in Figures 2 and 3. The optimum calculated values fall short of measured
values and this is due to the diversification arising from the way of windows
opening. Control of when and how much, windows are opening through the use of
the algorithm (1) are indicated in Figure 5. Time variation of window opening rate
leads to both reduce the temperature and reduce the humidity inside emissions (data
not shown).</p>
      <p>Figures 4 (a-d) are showing the variation of calculated values using the modified
TOMGRO, and the optimum values obtained after the use of branch &amp; bound
method for the number of nodes, dry weight of stem, leaves and whole plant
correspondingly. In Figure 4 (a) observed that the climate change does not change
the number of plant nodes formed. Correspondingly there is no change at the number
of leaves and the number of produced flowers (data not shown). Variation is
observed in Figure 4 (b) showing the dry weight of the shoot. The calculated optimal
values are above those calculated by the modified TOMGRO. Reverse change is
observed in the dry weight of leaves Figure 4 (c), wherein the optimum values are
below those estimated by the modified TOMGRO.</p>
      <p>Smaller differences observed in Figure 4 (d) illustrating the total output of the
plant biomass. The differences are due to the diversification of air temperature
values, Figure 2. Reduction of air temperature leads to a small hysteresis of growth
and development of tomato plants. The reduction in biomass production when the
plants are in the initial stage of development helps to create more robust plants and
there is a regular practice for the producers that mainly use the chemical composition
of the nutrient solution and the irrigation dose to achieve it.
3.2 Results of the second simulated period
4 Discussion &amp; Conclusions</p>
      <p>The results presented above were according two different climate scenarios within
greenhouses, with simultaneous display of the changes caused in the development
and production of the tomato crop. The periods used were selected of producer’s
interest and was at the start of the growing season and the second in the medium.</p>
      <p>The study of the results for the first scenario, it is established that the reduction of
air temperature contributes to the reduction of both dry weight of leaves and whole
plant. The decrease in air temperature was due to the different treatment of the
windows operating system (system controller and B &amp; B method).
10/10/2007
11/10/2007
12/10/2007
13/10/2007</p>
      <p>14/10/2007
Days</p>
      <p>The practice followed in accordance with the use of branch &amp; bound method leads
to an opening of a greenhouse window for longer period than that followed in the
experimental procedure. However, the reduction of biomass produced when the
plants are in the initial stage of development helps to create more robust plants in the
future.</p>
      <p>Days
Days
3
3
4
4
5
5</p>
      <p>Days
Days
3
3
4
4
5
5
0,6
1)tn 0,5
-a
l
(gp 0,4
s
tem 0,3
fso
tgh 0,2
i
e
ryw 0,1
D
0</p>
      <p>(a)
0</p>
      <p>1
(b)
0</p>
      <p>1
100%
)
(%90%
g
inn 80%
e
pO 70%
ow 60%
d
inW50%
feo 40%
g
tan 30%
e
rce 20%
P 10%
0%
15
16
17
18
19
20
25
26
27
28
29</p>
      <p>30
21 22 23 24</p>
      <p>Air Temperature (oC)</p>
      <p>By observing the results of the second scenario it was found that the reduction of
production costs by 19.72% led to a corresponding reduction of the biomass
production by 15.5% of shoot dry weight, 18.5% of leaves dry weight and 17.4% of
whole plant dry weight. However, the reduction caused in plants growth and
development can be balanced with an increase in air temperature inside the
greenhouse to a period prior to the reduction.
0
28/11/2007 29/11/2007 30/11/2007
1/12/2007
2/12/2007
3/12/2007</p>
      <p>4/12/2007</p>
      <p>Days</p>
      <p>Managing greenhouse climate as mentioned is a daily activity for the producers,
which despite its frequency poses many problems. Further analysis of possible
climate scenarios will help to create strategies that will shape the conditions for
reducing production costs and improving the climate inside the greenhouse units.
15,0
-1t)n14,0
oapn
l
(
se13,0
oodn
f
er12,0
b
m
u
N11,0
3,0
-1t)an2,5
lgp2,0
(se
tm1,5
fso
th1,0
eg
i
w0,5
ry
D
0,0
(a)
(b)
0
1</p>
      <p>2
0
1
2</p>
      <p>Days
Days
3
3
4
4
5
5
5,0
-1)tn4,0
a
l
p
g
(s3,0
eeva
lf2,0
o
t
h
g
ie1,0
w
ry
D
0,0
1)7,0
-ltan6,0
(gp
tn5,0
ap
l
fo4,0
t
h
ieg3,0
w
ryd2,0
lta
oT1,0
0,0
0
0
(c)
(d)
1
1
2
2</p>
      <p>Days
Days
3
3</p>
      <p>Acknowledgments. This paper is part of the 03ED526 research project, implemented
within the framework of the “Reinforcement Program of Human Research
100%
90%
80%
70%
60%
50%
40%
30%
20%
10%
0%
4
4
)
%
(
g
n
i
n
e
p
O
w
o
d
n
i
W
f
o
e
g
a
t
n
e
c
r
e
P
5
5
Manpower” (PENED) and co-financed by National and Community Funds (25%
from the Greek Ministry of Development-General Secretariat of Research and
Technology and 75% from E.U.-European Social Fund).</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Bailey</surname>
            ,
            <given-names>B.J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Seginer</surname>
            ,
            <given-names>I.</given-names>
          </string-name>
          , (
          <year>1989</year>
          )
          <article-title>Optimum control of greenhouse heating</article-title>
          .
          <source>Acta Horticulturae</source>
          <volume>245</volume>
          , p.
          <fpage>512</fpage>
          -
          <lpage>518</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Bellman</surname>
            ,
            <given-names>R.E.</given-names>
          </string-name>
          , (
          <year>1957</year>
          )
          <article-title>Dynamic programming</article-title>
          . Princeton, New Jersey, USA: Princeton University Press.
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>de Koning</surname>
            ,
            <given-names>A.N.M.</given-names>
          </string-name>
          , (
          <year>1988</year>
          )
          <article-title>The effect of different day/night temperature regimes on the growth, development and yield of glasshouse tomatoes</article-title>
          .
          <source>Journal of Horticultural Science</source>
          ,
          <volume>63</volume>
          (
          <issue>3</issue>
          ), p.
          <fpage>465</fpage>
          -
          <lpage>471</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Dimokas</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          , (
          <year>2009</year>
          )
          <article-title>Development of a biophysical simulator for tomato crop in order to optimize climate management for Mediterranean greenhouses during winter period</article-title>
          .
          <source>Ph.D. Thesis</source>
          , University of Thessaly, p.
          <fpage>1</fpage>
          -
          <lpage>332</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Dimokas</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tchamitchian</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kittas</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          , (
          <year>2009</year>
          )
          <article-title>Calibration and validation of a biological model to simulate the development and production of tomatoes in Mediterranean greenhouses during winter period</article-title>
          .
          <source>Biosystems Engineering</source>
          <volume>103</volume>
          (
          <issue>2</issue>
          ), p.
          <fpage>217</fpage>
          -
          <lpage>227</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Heuvelink</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          , (
          <year>1989</year>
          )
          <article-title>Influence of day and night temperature on the growth of young tomato plants</article-title>
          .
          <source>Scientia Horticulturae</source>
          ,
          <volume>38</volume>
          , p.
          <fpage>11</fpage>
          -
          <lpage>22</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Pontryagin</surname>
            ,
            <given-names>L.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Boltyanskii</surname>
            ,
            <given-names>V.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gamkrelidze</surname>
            ,
            <given-names>R.V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mishchenko</surname>
            ,
            <given-names>E.F.</given-names>
          </string-name>
          , (
          <year>1962</year>
          )
          <article-title>The mathematical theory of optimal processes</article-title>
          . New-York, USA: Wiley Interscience.
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Rijsdijk</surname>
            ,
            <given-names>A.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Vogelezang</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          , (
          <year>2000</year>
          )
          <article-title>Temperature integration on a 24 hour basis: a more efficient climate control strategy</article-title>
          .
          <source>Acta Horticulturae</source>
          ,
          <volume>519</volume>
          , p.
          <fpage>163</fpage>
          -
          <lpage>170</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Tap</surname>
          </string-name>
          , R.F.,
          <string-name>
            <surname>van Willigenburg</surname>
          </string-name>
          , L.G., van Straten, G., van
          <string-name>
            <surname>Henten</surname>
            ,
            <given-names>E.J.</given-names>
          </string-name>
          , (
          <year>1993</year>
          )
          <article-title>Optimal control of greenhouse climate: computation of the influence of fast and slow dynamics</article-title>
          .
          <source>12th IFAC Triennial World Congress, Sydney, Australia</source>
          , Volume
          <volume>4</volume>
          ,
          <fpage>4pp</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Tchamitchian</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tantau</surname>
            ,
            <given-names>H.J.</given-names>
          </string-name>
          , (
          <year>1996</year>
          )
          <article-title>Optimal control of the daily greenhouse climate: physical approach</article-title>
          .
          <source>13th IFAC Triennial World Congress</source>
          , San Francisco, USA, Volume B, p.
          <fpage>471</fpage>
          -
          <lpage>475</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>van Henten</surname>
            ,
            <given-names>E.J.</given-names>
          </string-name>
          , (
          <year>1994</year>
          )
          <article-title>Greenhouse climate management: An optimal control approach</article-title>
          .
          <source>Ph.D. Thesis</source>
          , Wageningen Agricultural Univericity, Wageningen, The Netherlands, p.
          <fpage>1</fpage>
          -
          <lpage>319</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Vogelezang</surname>
          </string-name>
          , J., de Hoog, J.,
          <string-name>
            <surname>Marissen</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          , (
          <year>2000</year>
          )
          <article-title>Effects of diurnal temperature strategies on carbohydrate content and flower quality of greenhouse roses</article-title>
          .
          <source>Acta Horticulturae</source>
          <volume>515</volume>
          , p.
          <fpage>111</fpage>
          -
          <lpage>118</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>