# Temperature Field-Wind Velocity Field Optimum Control of Greenhouse Environment Based on CFD Model

Authors: Yongbo Li; Guoxiang Sun; Xiaochan Wang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The computational fluid dynamics technology is applied as the environmental control model, which can include the greenhouse space. Basic environmental factors are set to be the control objects, the field information is achieved via the division of layers by height, and numerical characteristics of each layer are used to describe the field information. Under the natural ventilation condition, real-time requirements, energy consumption, and distribution difference are selected as index functions. The optimization algorithm of adaptive simulated annealing is used to obtain optimal control outputs. A comparison with full-open ventilation shows that the whole index can be reduced at 44.21% and found that a certain mutual exclusiveness exists between the temperature and velocity field in the optimal course. All the results indicate that the application of CFD model has great advantages to improve the control accuracy of greenhouse.


---

## Body

## 1. Introduction

In the greenhouse practice, the difference of growth and quality for the crop is an ordinary phenomenon, which could lead to a reduction of overall crop production [1]. A simulation based on Tomgro model also clearly shows that crops generally display the fast growth rate, high yield, and excellent quality in the central area; whereas crops in the peripheral area have the relatively slow growth rate, low yield, and worsening quality [2]. This unequilibrium has a lot of reasons, but the control system must be involved. Because it is generally assumed that the internal air in greenhouse is thoroughly mixed and homogeneous during the design of control system [3]. Based on this assumption, all models are applied in the design process.In the early stage, the design of one-factor control system is generally based on the transfer function model [4]; and later, the design of the multifactor control system is developed in terms of the state space model [5]. Recently, the intelligent model has been more commonly applied [6]. Although these models are different, their designs are consistently based on the aforementioned assumptions. As a result, during the closed-loop feedback process with environmental factors, the “point” form is commonly used (sampling via a single sensor), and thus one feedback value is used to represent environmental conditions of the entire greenhouse space. However crops in the greenhouse are spatially distributed, and environmental factors at different locations could be manifested by different numerical values, the point forms to describe crops’ requirements are inaccurate. So the field forms should be chosen to fit the real situations; in order to further improve the control precision of greenhouse environment, the selection of control model to explicitly reflect the “field” form is badly needed.In recent years, the CFD technology has been fully developed and applied in the field of greenhouse environment modeling. The earliest exploration of the applicability and effectiveness of CFD technology for the simulation of greenhouse environment was studied by the literature [7]. During this period, greenhouse models were generally used in two-dimensional forms due to limitations of the computer hardware and software [8, 9]. Later, with the dramatic increase of computational speed and improvement of algorithm efficiency, the greenhouse CFD model was extended to three-dimensional form [10].Meanwhile, a variety of scenarios, such as greenhouse structures [11–13], geographical features [14, 15], environmental factors [16–18], and the existence/absence of crops [19, 20], were included in the construction of greenhouse CFD model, and the results generally displayed that utilizing CFD technology to establish greenhouse environment model is effective. In theory, CFD model is based on basic principles such as the conservation of energy, mass, and momentum. During the establishment of the model, the nonlinearity and coupling features of greenhouse have been fully taken into account; hence, the description precision of the model is high. Another advantage lies in the fact that, through mesh generation, the factor state of any points within the greenhouse space could be expressed theoretically. Therefore, the integration of CFD in greenhouse environment model could meet the requirement of “field” form, which could help planters to control the overall greenhouse state, not to limit it to one point for part of crops.Currently, the use of CFD as the greenhouse environment model is commonly applied for the short-term forecasting and indirect guidance of environment control [11, 13, 14, 16]. In the literature [21–25], the design of experiment (DOE) and genetic algorithm (GA) to control greenhouse temperature field under CFD model were applied. However, these studies are merely limited to a single factor, which could not provide a complete description of the state of greenhouse environment. Hence, this paper selects two factors (temperature and wind velocity) as the greenhouse control object under the CFD greenhouse model. Different from the common design, the wind velocity was chosen as a control variable. The reason is that the wind velocity is a very important factor if we want to control the distributions of humidity and CO2 concentration. This point has been testified by our experiments and former research literatures. Therefore, the control for the distributions of humidity or CO2 concentration can be changed accordingly to control the wind velocity.The paper’s structure is as follows. Firstly, the natural ventilation is chosen as the control pattern to set up the whole system. The field information which includes the maximum value, average value, minimum value, and standard deviation is achieved by digitizing the height of greenhouse space. Secondly, with a careful consideration, a multiple index system for greenhouse control is obtained, which contained the energy confinement, real-time control, temperature, wind velocity field, and so forth. Finally, adaptive simulated annealing (ASA) algorithm is applied as the optimization algorithm to solve the optimal control problem and get the natural ventilation strategy of using the CFD greenhouse model for the temperature field-wind velocity field (TF-WVF) under the condition of multiconstraints and multitargets.


## 2. Greenhouse Model

### 2.1. Physical Model

The greenhouse is located in Nanjing Agricultural University (E 118°46′, E 32°03′). It is 2-ridge Venlo structure with aluminium alloy material, and size is 8 m (width) × 16 m (length); gutter and ridge height are 4.2 m and 5.2 m, respectively; the total area is 128 m2. The greenhouse is covered with 4 mm float glass whose light transmission is over 88%. Environment control devices include window opening, external/internal sun shading, spraying, and drip irrigation. The natural ventilation consists of 6 windows that open to the north and 6 that open in the opposite direction and an inlet which is located in west front; all of windows have a synchronous moving mode. Geometric structure of greenhouse is shown in Figure 1.Figure 1
Greenhouse physics structure.


### 2.2. Greenhouse CFD Model

For setting up the CFD model in Figure1, continuity equation, momentum equation, energy equation, three-dimensional turbulence equation, radiation equation, and solar load model are embedded in greenhouse model. During the CFD calculation process, it is assumed that the temperature is a well-distribution even through glass layers of the outer covering; translucent materials are used for the construction of the roof and sides area; the concrete structure is used for the ground; the absorption rate of radiation energy remains constant; and the heat is transferred via the form of fixed absorption [21].


### 2.3. Description of Calculation Software and Hardware

The mesh and fluid calculations were made by using the software: Gambit 2.4.6 and Fluent 6.3.26, respectively. One velocity inlet, twelve skylight pressure outlets, five glass walls, and one concrete wall (ground) are set up. The computer mode is Lenovo Ideapad Z370, with CPU as Intel (R) Celeron (TM) I3-2350M 4 × 2.30 GHz, memory as DDR3 1333 MHz 4 GB, and hard disk as ST500LM012 with the total of 500 G storage. The operating system is 32-bit ultimate version of Windows 7, and four-thread parallel computing is adopted during the process of Fluent iterative calculation. In the experiments, the initial conditions and boundary conditions in Fluent and Gambit were set as Table1.Table 1
Initial conditions and boundary conditions.
Parameters name (unit)
Value
Air
Density (kg·m−3)
1.225
Heat coefficient (W·(m·K)−1)
0.0242
Thermal expansion coefficient (K−1)
3.356
e
-
3
Specific heat capacity (kJ·(kg·K)−1)
1006.43
Glass
Density (kg·m−3)
2220
Heat coefficient (W·(m·K)−1)
1.15
Specific heat capacity (kJ·(kg·K)−1)
830
Thickness (m)
0.004
Boundary conditions
Inlet velocity (m·s−1)
2.0
Inlet temperature (K)
300
Outdoor temperature (K)
305
Glass wall temperature (K)
300
Ground wall temperature (K)
290


## 2.1. Physical Model

The greenhouse is located in Nanjing Agricultural University (E 118°46′, E 32°03′). It is 2-ridge Venlo structure with aluminium alloy material, and size is 8 m (width) × 16 m (length); gutter and ridge height are 4.2 m and 5.2 m, respectively; the total area is 128 m2. The greenhouse is covered with 4 mm float glass whose light transmission is over 88%. Environment control devices include window opening, external/internal sun shading, spraying, and drip irrigation. The natural ventilation consists of 6 windows that open to the north and 6 that open in the opposite direction and an inlet which is located in west front; all of windows have a synchronous moving mode. Geometric structure of greenhouse is shown in Figure 1.Figure 1
Greenhouse physics structure.


## 2.2. Greenhouse CFD Model

For setting up the CFD model in Figure1, continuity equation, momentum equation, energy equation, three-dimensional turbulence equation, radiation equation, and solar load model are embedded in greenhouse model. During the CFD calculation process, it is assumed that the temperature is a well-distribution even through glass layers of the outer covering; translucent materials are used for the construction of the roof and sides area; the concrete structure is used for the ground; the absorption rate of radiation energy remains constant; and the heat is transferred via the form of fixed absorption [21].


## 2.3. Description of Calculation Software and Hardware

The mesh and fluid calculations were made by using the software: Gambit 2.4.6 and Fluent 6.3.26, respectively. One velocity inlet, twelve skylight pressure outlets, five glass walls, and one concrete wall (ground) are set up. The computer mode is Lenovo Ideapad Z370, with CPU as Intel (R) Celeron (TM) I3-2350M 4 × 2.30 GHz, memory as DDR3 1333 MHz 4 GB, and hard disk as ST500LM012 with the total of 500 G storage. The operating system is 32-bit ultimate version of Windows 7, and four-thread parallel computing is adopted during the process of Fluent iterative calculation. In the experiments, the initial conditions and boundary conditions in Fluent and Gambit were set as Table1.Table 1
Initial conditions and boundary conditions.
Parameters name (unit)
Value
Air
Density (kg·m−3)
1.225
Heat coefficient (W·(m·K)−1)
0.0242
Thermal expansion coefficient (K−1)
3.356
e
-
3
Specific heat capacity (kJ·(kg·K)−1)
1006.43
Glass
Density (kg·m−3)
2220
Heat coefficient (W·(m·K)−1)
1.15
Specific heat capacity (kJ·(kg·K)−1)
830
Thickness (m)
0.004
Boundary conditions
Inlet velocity (m·s−1)
2.0
Inlet temperature (K)
300
Outdoor temperature (K)
305
Glass wall temperature (K)
300
Ground wall temperature (K)
290


## 3. Control Objects in Field Form

Under the “point” feedback form, the control objective is a constant, whereas under the “field” feedback form, the control object needs to be expanded that it should thoroughly reflect the spatial distribution of greenhouse environmental factors. So we will decide the control object from two aspects: constraint and index.


### 3.1. Control Constraints

#### 3.1.1. Time Constraint of CFD Computing

Since the CFD model is a nonlinear system, its calculation process is time-consuming. In order to enhance the feasibility in practice, constraints should contain model parameters, the selection of equations, the choice of appropriate algorithms, the determination of the convergence level, and so forth. Here the calculation time of CFD model should meet the requirements of (1) as follows:
(1)
T
CFD
<
T
s
,
where T
CFD represents the required computing time of CFD model during the optimization process and T
s represents the sampling period of the control system. Additionally, (1) requires sufficient calculation effectiveness of CFD model in Section 2.2 so that the manipulation time under the realistic setting could be met, and the optimization requirement of control algorithm for subsequent steps could be provided.


#### 3.1.2. Constraint of TF-WVF

In order to achieve the extended form of “point” feedback to “field” form, this paper uses height as the reference and chooses the extreme values (maximum and minimum value) and average value obtained from different horizontal planes as numerical characteristics (NC) to describe the TF-WVF information database based on CFD model. Expressions of TF and WVF are given as follows.(
1
) NC of TF:
(2)
max
⁡
T
=
max
⁡
T
i
,
min
⁡
T
=
min
⁡
T
i
,
ave
T
=
ave
T
i
,
where max
⁡
T, min
⁡
T, and ave
T represent the maximum temperature, minimum temperature, and average temperature, respectively; these values correspond to the height of h of the horizontal plane, i represents the mesh point number after mesh generation, and T
i represents the temperature feedback value of point i.(
2
) NC of WVF:
(3)
max
⁡
V
=
max
⁡
V
→
i
,
min
⁡
V
=
min
⁡
V
→
i
,
ave
V
=
ave
V
→
i
,
where max
⁡
V, min
⁡
V, and ave
V represent the maximum velocity, minimum velocity, and average velocity, respectively, which corresponds to the horizontal plane with a height of h and V
→
i represents the wind velocity vector of point i after mesh generation of that plane.It can be seen that, different from the temperature field, the wind velocity field is in a vector format, and thus its amplitude value is calculated in accordance with the following:(4)
V
→
i
=
V
i
x
2
+
V
i
y
2
+
V
i
z
2
,
where V
i
x, V
i
y, and V
i
z represent the velocity components of point i along three directions.Although the range of heighth
∈
[
0
y
max
⁡
], which is continuous (y
max
⁡ represents the maximum height of the greenhouse), growers’ major concern is about some key heights in practice, such as ground layer, intermediate layer, leaf layer, and canopy layer. Considering our experimental object is the tomato in greenhouse. In this paper, we divided the horizontal planes into six different levels including 0.2 m, 0.4 m, 0.8 m, 1.0 m, 1.2 m, and 2.0 m.


### 3.2. Control Indexes Set

During the process of greenhouse production, many aspects should be involved in control system. This paper explicitly considers three aspects, including the control cost index, distribution difference index, and overall index; the details of indexes are given as follows.


#### 3.2.1. Control Cost Index

The control cost index refers to the amount of energy consumption of related equipment under natural ventilation conditions. In this paper, both the greenhouse entrance and the skylight are controlled by motors which are the same mode. Therefore, the control cost index could be described as follows:(5)
J
u
=
∑
k
=
1
3
λ
k
u
k
,
where J
u represents the overall cost index. k
=
1,2
,
3 represents the control output of entrance, north skylight, and south skylight, respectively. λ
k represents respective proportion of motor power with ∑
k
=
1
3
λ
k
=
1. The considered three kinds of motors belong to the same mode with λ
1
=
λ
2
=
λ
3
=
1
/
3. u
k
∈
0
100
% represents respective openness.Note. The experimental greenhouse is small and the cost of electric motor is less. However, in the practical greenhouse, especially greenhouse-group, the cost of running electric motor is a very important factor for planters’ income.


#### 3.2.2. Distribution Difference Index

The distribution difference index refers to the difference in the distribution of environmental factors under a certain level of heighth; smaller index values indicate that the wind velocity or temperature is evenly distributed at the plane; thus the variation in the growth status of crops could be reduced. Therefore, if the control system can regulate the devices to decrease this index which satisfied the demand of all crops in greenhouse, the greenhouse yield will be enlarged and the precision of control system will be dramatically enhanced. This paper uses the mean standard deviation to describe the distribution difference level of TF-WVF.(
1
) TF distribution difference:
(6)
J
T
=
1
N
∑
i
=
1
N
T
i
-
ave
T
2
,
where J
T represents the distribution difference of temperature and N represents the number of nodes at the height of h after being meshed by Gambit.Considering layers are limited, stratified (6) can get
(7)
J
T
=
∑
m
=
1
6
a
m
J
T
m
,
where a
m represents the weight coefficient at layer m with ∑
m
=
1
6
a
m
=
1 and J
T
m represents the distribution difference index at layer m.(
2
) WVF distribution difference:
(8)
J
V
=
1
N
∑
i
=
1
N
V
→
i
-
ave
V
2
,
where J
V represents WVF distribution difference and the calculation of V
→
i is referred to (4). After similar stratification and description process, we have
(9)
J
V
=
∑
m
=
1
6
b
m
J
V
m
,
where b
m represents the weight coefficient of layer m with ∑
m
=
1
6
b
m
=
1 and J
V
m represents distribution difference index at layer m.


#### 3.2.3. Overall Index

In this paper, two involved environmental factors are the wind velocity and temperature. Meanwhile, the control indexes include real-time, cost, and distribution difference. Therefore, the control system is multiconstraint and multi-index. In order to take all these characteristics into account, the following overall index is proposed:(10)
J
=
α
J
T
+
β
J
V
+
γ
J
u
,
where α, β, and γ are weight coefficient of three indexes. Here, we set α
=
β
=
0.4 and γ
=
0.2.


## 3.1. Control Constraints

### 3.1.1. Time Constraint of CFD Computing

Since the CFD model is a nonlinear system, its calculation process is time-consuming. In order to enhance the feasibility in practice, constraints should contain model parameters, the selection of equations, the choice of appropriate algorithms, the determination of the convergence level, and so forth. Here the calculation time of CFD model should meet the requirements of (1) as follows:
(1)
T
CFD
<
T
s
,
where T
CFD represents the required computing time of CFD model during the optimization process and T
s represents the sampling period of the control system. Additionally, (1) requires sufficient calculation effectiveness of CFD model in Section 2.2 so that the manipulation time under the realistic setting could be met, and the optimization requirement of control algorithm for subsequent steps could be provided.


### 3.1.2. Constraint of TF-WVF

In order to achieve the extended form of “point” feedback to “field” form, this paper uses height as the reference and chooses the extreme values (maximum and minimum value) and average value obtained from different horizontal planes as numerical characteristics (NC) to describe the TF-WVF information database based on CFD model. Expressions of TF and WVF are given as follows.(
1
) NC of TF:
(2)
max
⁡
T
=
max
⁡
T
i
,
min
⁡
T
=
min
⁡
T
i
,
ave
T
=
ave
T
i
,
where max
⁡
T, min
⁡
T, and ave
T represent the maximum temperature, minimum temperature, and average temperature, respectively; these values correspond to the height of h of the horizontal plane, i represents the mesh point number after mesh generation, and T
i represents the temperature feedback value of point i.(
2
) NC of WVF:
(3)
max
⁡
V
=
max
⁡
V
→
i
,
min
⁡
V
=
min
⁡
V
→
i
,
ave
V
=
ave
V
→
i
,
where max
⁡
V, min
⁡
V, and ave
V represent the maximum velocity, minimum velocity, and average velocity, respectively, which corresponds to the horizontal plane with a height of h and V
→
i represents the wind velocity vector of point i after mesh generation of that plane.It can be seen that, different from the temperature field, the wind velocity field is in a vector format, and thus its amplitude value is calculated in accordance with the following:(4)
V
→
i
=
V
i
x
2
+
V
i
y
2
+
V
i
z
2
,
where V
i
x, V
i
y, and V
i
z represent the velocity components of point i along three directions.Although the range of heighth
∈
[
0
y
max
⁡
], which is continuous (y
max
⁡ represents the maximum height of the greenhouse), growers’ major concern is about some key heights in practice, such as ground layer, intermediate layer, leaf layer, and canopy layer. Considering our experimental object is the tomato in greenhouse. In this paper, we divided the horizontal planes into six different levels including 0.2 m, 0.4 m, 0.8 m, 1.0 m, 1.2 m, and 2.0 m.


## 3.1.1. Time Constraint of CFD Computing

Since the CFD model is a nonlinear system, its calculation process is time-consuming. In order to enhance the feasibility in practice, constraints should contain model parameters, the selection of equations, the choice of appropriate algorithms, the determination of the convergence level, and so forth. Here the calculation time of CFD model should meet the requirements of (1) as follows:
(1)
T
CFD
<
T
s
,
where T
CFD represents the required computing time of CFD model during the optimization process and T
s represents the sampling period of the control system. Additionally, (1) requires sufficient calculation effectiveness of CFD model in Section 2.2 so that the manipulation time under the realistic setting could be met, and the optimization requirement of control algorithm for subsequent steps could be provided.


## 3.1.2. Constraint of TF-WVF

In order to achieve the extended form of “point” feedback to “field” form, this paper uses height as the reference and chooses the extreme values (maximum and minimum value) and average value obtained from different horizontal planes as numerical characteristics (NC) to describe the TF-WVF information database based on CFD model. Expressions of TF and WVF are given as follows.(
1
) NC of TF:
(2)
max
⁡
T
=
max
⁡
T
i
,
min
⁡
T
=
min
⁡
T
i
,
ave
T
=
ave
T
i
,
where max
⁡
T, min
⁡
T, and ave
T represent the maximum temperature, minimum temperature, and average temperature, respectively; these values correspond to the height of h of the horizontal plane, i represents the mesh point number after mesh generation, and T
i represents the temperature feedback value of point i.(
2
) NC of WVF:
(3)
max
⁡
V
=
max
⁡
V
→
i
,
min
⁡
V
=
min
⁡
V
→
i
,
ave
V
=
ave
V
→
i
,
where max
⁡
V, min
⁡
V, and ave
V represent the maximum velocity, minimum velocity, and average velocity, respectively, which corresponds to the horizontal plane with a height of h and V
→
i represents the wind velocity vector of point i after mesh generation of that plane.It can be seen that, different from the temperature field, the wind velocity field is in a vector format, and thus its amplitude value is calculated in accordance with the following:(4)
V
→
i
=
V
i
x
2
+
V
i
y
2
+
V
i
z
2
,
where V
i
x, V
i
y, and V
i
z represent the velocity components of point i along three directions.Although the range of heighth
∈
[
0
y
max
⁡
], which is continuous (y
max
⁡ represents the maximum height of the greenhouse), growers’ major concern is about some key heights in practice, such as ground layer, intermediate layer, leaf layer, and canopy layer. Considering our experimental object is the tomato in greenhouse. In this paper, we divided the horizontal planes into six different levels including 0.2 m, 0.4 m, 0.8 m, 1.0 m, 1.2 m, and 2.0 m.


## 3.2. Control Indexes Set

During the process of greenhouse production, many aspects should be involved in control system. This paper explicitly considers three aspects, including the control cost index, distribution difference index, and overall index; the details of indexes are given as follows.


### 3.2.1. Control Cost Index

The control cost index refers to the amount of energy consumption of related equipment under natural ventilation conditions. In this paper, both the greenhouse entrance and the skylight are controlled by motors which are the same mode. Therefore, the control cost index could be described as follows:(5)
J
u
=
∑
k
=
1
3
λ
k
u
k
,
where J
u represents the overall cost index. k
=
1,2
,
3 represents the control output of entrance, north skylight, and south skylight, respectively. λ
k represents respective proportion of motor power with ∑
k
=
1
3
λ
k
=
1. The considered three kinds of motors belong to the same mode with λ
1
=
λ
2
=
λ
3
=
1
/
3. u
k
∈
0
100
% represents respective openness.Note. The experimental greenhouse is small and the cost of electric motor is less. However, in the practical greenhouse, especially greenhouse-group, the cost of running electric motor is a very important factor for planters’ income.


### 3.2.2. Distribution Difference Index

The distribution difference index refers to the difference in the distribution of environmental factors under a certain level of heighth; smaller index values indicate that the wind velocity or temperature is evenly distributed at the plane; thus the variation in the growth status of crops could be reduced. Therefore, if the control system can regulate the devices to decrease this index which satisfied the demand of all crops in greenhouse, the greenhouse yield will be enlarged and the precision of control system will be dramatically enhanced. This paper uses the mean standard deviation to describe the distribution difference level of TF-WVF.(
1
) TF distribution difference:
(6)
J
T
=
1
N
∑
i
=
1
N
T
i
-
ave
T
2
,
where J
T represents the distribution difference of temperature and N represents the number of nodes at the height of h after being meshed by Gambit.Considering layers are limited, stratified (6) can get
(7)
J
T
=
∑
m
=
1
6
a
m
J
T
m
,
where a
m represents the weight coefficient at layer m with ∑
m
=
1
6
a
m
=
1 and J
T
m represents the distribution difference index at layer m.(
2
) WVF distribution difference:
(8)
J
V
=
1
N
∑
i
=
1
N
V
→
i
-
ave
V
2
,
where J
V represents WVF distribution difference and the calculation of V
→
i is referred to (4). After similar stratification and description process, we have
(9)
J
V
=
∑
m
=
1
6
b
m
J
V
m
,
where b
m represents the weight coefficient of layer m with ∑
m
=
1
6
b
m
=
1 and J
V
m represents distribution difference index at layer m.


### 3.2.3. Overall Index

In this paper, two involved environmental factors are the wind velocity and temperature. Meanwhile, the control indexes include real-time, cost, and distribution difference. Therefore, the control system is multiconstraint and multi-index. In order to take all these characteristics into account, the following overall index is proposed:(10)
J
=
α
J
T
+
β
J
V
+
γ
J
u
,
where α, β, and γ are weight coefficient of three indexes. Here, we set α
=
β
=
0.4 and γ
=
0.2.


## 3.2.1. Control Cost Index

The control cost index refers to the amount of energy consumption of related equipment under natural ventilation conditions. In this paper, both the greenhouse entrance and the skylight are controlled by motors which are the same mode. Therefore, the control cost index could be described as follows:(5)
J
u
=
∑
k
=
1
3
λ
k
u
k
,
where J
u represents the overall cost index. k
=
1,2
,
3 represents the control output of entrance, north skylight, and south skylight, respectively. λ
k represents respective proportion of motor power with ∑
k
=
1
3
λ
k
=
1. The considered three kinds of motors belong to the same mode with λ
1
=
λ
2
=
λ
3
=
1
/
3. u
k
∈
0
100
% represents respective openness.Note. The experimental greenhouse is small and the cost of electric motor is less. However, in the practical greenhouse, especially greenhouse-group, the cost of running electric motor is a very important factor for planters’ income.


## 3.2.2. Distribution Difference Index

The distribution difference index refers to the difference in the distribution of environmental factors under a certain level of heighth; smaller index values indicate that the wind velocity or temperature is evenly distributed at the plane; thus the variation in the growth status of crops could be reduced. Therefore, if the control system can regulate the devices to decrease this index which satisfied the demand of all crops in greenhouse, the greenhouse yield will be enlarged and the precision of control system will be dramatically enhanced. This paper uses the mean standard deviation to describe the distribution difference level of TF-WVF.(
1
) TF distribution difference:
(6)
J
T
=
1
N
∑
i
=
1
N
T
i
-
ave
T
2
,
where J
T represents the distribution difference of temperature and N represents the number of nodes at the height of h after being meshed by Gambit.Considering layers are limited, stratified (6) can get
(7)
J
T
=
∑
m
=
1
6
a
m
J
T
m
,
where a
m represents the weight coefficient at layer m with ∑
m
=
1
6
a
m
=
1 and J
T
m represents the distribution difference index at layer m.(
2
) WVF distribution difference:
(8)
J
V
=
1
N
∑
i
=
1
N
V
→
i
-
ave
V
2
,
where J
V represents WVF distribution difference and the calculation of V
→
i is referred to (4). After similar stratification and description process, we have
(9)
J
V
=
∑
m
=
1
6
b
m
J
V
m
,
where b
m represents the weight coefficient of layer m with ∑
m
=
1
6
b
m
=
1 and J
V
m represents distribution difference index at layer m.


## 3.2.3. Overall Index

In this paper, two involved environmental factors are the wind velocity and temperature. Meanwhile, the control indexes include real-time, cost, and distribution difference. Therefore, the control system is multiconstraint and multi-index. In order to take all these characteristics into account, the following overall index is proposed:(10)
J
=
α
J
T
+
β
J
V
+
γ
J
u
,
where α, β, and γ are weight coefficient of three indexes. Here, we set α
=
β
=
0.4 and γ
=
0.2.


## 4. The Optimal Control

Constraints and indexes obtained under field form conditions; the natural ventilation optimization control system based on CFD model could be described as shown in Figure2.Figure 2
Principle of control system.Figure2 displays a closed-loop control system. The CFD model and the control system are linked via the extraction of field information which is based on Section 3. The optimization control system applies ASA algorithm to solve the following optimization problem with a wide stochastic optimal course as follows:
(11)
min
⁡
J
=
min
⁡
α
J
T
+
β
J
V
+
γ
J
u
s
.
t
.
h
i
T
CFD
<
T
S
hhhhh
h
ave
T
∈
Ω
T
hhhhh
h
ave
V
∈
Ω
V
,
where Ω
T
=
295
±
2
K represents the average temperature requirement of the corresponding layers and Ω
V
∈
0.3
±
0.1
m
/
s represents the average wind velocity. All the parameters set in this section are only to testify the feasibility of the CFD technology.


## 5. Results and Discussions

The natural ventilation is an ordinary regulation method for the greenhouse control, and in practice, it is often in full-action (full open or close) or inaccurate controlled by the planter. So in this paper, the continuous and accurate action mode will be studied. To make a comparison, the field information and related indexes under full open natural ventilation (FONV) conditions are researched firstly. Next, under the same indexes the optimal control system will be analyzed completely. By comparison and analysis, the effective information would be collected for the guidance of system design in practice.


### 5.1. FONV

The meaning of FONV isu
1
=
u
2
=
u
3
=
100
%. Under this condition, the previous mentioned greenhouse environment information which includes NCs of wind velocity-temperature factors and related indexes could be calculated. Some results are given in Table 2.Table 2
NCs of FONV.
h (m)
0.2
0.4
0.8
1.0
1.2
2.0
T (K)
Min
294.35
294.43
294.54
294.59
294.64
294.70
Ave
294.75
294.89
294.94
294.96
294.98
295.01
Max
296.20
295.95
295.73
295.71
295.67
295.80
V (m/s)
Min
0.014299
0.018991
0.031150
0.032261
0.036340
0.056740
Ave
0.832410
0.817559
0.803519
0.799759
0.796762
0.756358
Max
3.042729
3.131846
3.047643
3.035799
3.035799
3.017646
SDT (K)
0.18608
0.15697
0.13566
0.12649
0.12165
0.11381
SDV (m/s)
0.76187
0.75538
0.74742
0.74337
0.74227
0.71905Table2 represents the numerical characteristics and related indexes of the six layers. It can be seen from this table that for either temperature or wind velocity, the corresponding average value closely relates with the height. That is to say, the average value of temperature increases with the increasing height, whereas the average value of wind velocity decreases with the increasing height, indicating that among greenhouse environmental factors, the temperature and wind velocity are negatively correlated. Meanwhile, there is a negative correlation between the temperature and wind velocity distribution too. From Table 2, it is clear that J
T
=
0.1357, J
V
=
0.7461, and J
u
=
1.0. For Table 2, the overall index under the optimal scenario is 0.5528. All of these results under full ventilation conditions will be used as a reference to compare.For a visual expression, we choose the middle layer (height is 1.0 m) to display the cloud map which contains the temperature and wind velocity in Figure3.Cloud map of FONV at 1.0 m.
(a)
Temperature
(b)
Velocity ofx direction
(c)
Velocity ofz directionFigure3(a) shows the temperature cloud map with h
=
1.0
m. It can be seen from this figure that, for most part of the layer, the average temperature is around 295 K with some local variation. Figures 3(b) and 3(c) show the wind velocity contours along x and z directions, respectively. Clearly, the velocity along x direction is slower than that along the z direction; the reason could be explained by the direction of the greenhouse entrance. Meanwhile, two contours display significant turbulence phenomena, which are quite similar to Figure 3(a), indicating that temperature and wind velocity factors are coupled strongly.


### 5.2. Optimal Scenario

In this pattern, the optimization processes spend about 61 minutes to attain the control outputs (u
1
=
22.62
%, u
2
=
93.59
%, and u
3
=
10.95
%). The temperature-wind velocity cloud map at height of h
=
1.0
m is shown in Figure 4, and the NCs of the six layers are shown in Table 3.Table 3
NCs of optimal scenario.
h (m)
0.2
0.4
0.8
1.0
1.2
2.0
T (K)
Min
293.65
293.97
294.20
294.24
294.33
294.46
Ave
294.88
295.02
295.09
295.11
295.13
295.17
Max
296.33
296.53
296.21
296.13
296.09
296.52
V (m/s)
Min
0.001656
0.002452
0.007214
0.008335
0.010763
0.014362
Ave
0.270686
0.254144
0.233365
0.224360
0.215042
0.157519
Max
2.997552
3.063393
3.002424
2.955838
2.941845
3.117379
SDT (K)
0.32906
0.28879
0.24432
0.22823
0.21794
0.21176
SDV (m/s)
0.33999
0.33691
0.32172
0.31735
0.30061
0.22776Cloud map of scenario 1 at 1.0 m.
(a)
Temperature
(b)
Velocity ofx direction
(c)
Velocity ofz directionIt can be seen from the cloud map that the isotherms include 294.5 K, 295 K, and 295.5 K. Although the distribution is similarly compared to FONV, the distribution difference is reduced, and the double-eddy phenomenon could be found in the wind field that the amplitude is reduced as well. NCs of six different layers are calculated, including the average temperature, the average velocity, SDT, and SDV. Table3 has a similar regulation with Table 2.


### 5.3. Comparing and Conclusions

To make a comparison, we choose (12) to define the difference between FONV and optimal scenario as follows:(12)
Difference
=
Optimal
-
FONV
.Note. The horizontal axis is about the height of the layer, and the vertical layer is the average temperature, the average velocity, SDT, or SDV, respectively.Figure5 shows that both the average temperature and SDT are positive. Therefore, during the optimization process, the average temperature and the temperature distribution are both positive at each layer; these results show that the indexes of temperature are not better than the FONV. In other words, full ventilation mode represents the best ventilation form to reduce distribution for temperature, but only in temperature. Meanwhile, the average velocity was decreased because the u
1 is less than 100%; it led to the decrease of air exchange. The SDV is negative and indicates that the velocity distribution became better than FONV; these results can be explained by turbulent theory that the turbulence intensity has the same change direction with u. This phenomenon can also be conformed from an extreme case: when the greenhouse is closed the wind will not flow, so the ave
V
≈
0 and SDV
≈
0; this is the best way to reduce SDV but does not suit the greenhouse control. To further enrich the information differences between FONV and optimal scenario, we make Table 4 to describe it.Table 4
Compared result of indexes with FONV.
FONV
Optimal scenario
Value
Variation
J
T
0.1357
0.2457
+80.95%
J
V
0.7461
0.3154
−57.72%
J
u
1.0000
0.4196
−58.04%
J
0.5528
0.3084
−44.21%Figure 5
NCs difference between optimal scenario and FONV.In the optimal scenario, the values ofα
=
β
=
0.4 and γ
=
0.2 indicate that this mode belongs to a comprehensive optimization system, which includes the temperature, velocity, and control cost. Therefore, this scenario has a wide variety of applications in practice. Except for J
T, values of all other indexes are reduced, with the overall index reduced 44.21%, indicating that the optimization system in this paper has a great superiority for application in terms of greenhouse environment control.


## 5.1. FONV

The meaning of FONV isu
1
=
u
2
=
u
3
=
100
%. Under this condition, the previous mentioned greenhouse environment information which includes NCs of wind velocity-temperature factors and related indexes could be calculated. Some results are given in Table 2.Table 2
NCs of FONV.
h (m)
0.2
0.4
0.8
1.0
1.2
2.0
T (K)
Min
294.35
294.43
294.54
294.59
294.64
294.70
Ave
294.75
294.89
294.94
294.96
294.98
295.01
Max
296.20
295.95
295.73
295.71
295.67
295.80
V (m/s)
Min
0.014299
0.018991
0.031150
0.032261
0.036340
0.056740
Ave
0.832410
0.817559
0.803519
0.799759
0.796762
0.756358
Max
3.042729
3.131846
3.047643
3.035799
3.035799
3.017646
SDT (K)
0.18608
0.15697
0.13566
0.12649
0.12165
0.11381
SDV (m/s)
0.76187
0.75538
0.74742
0.74337
0.74227
0.71905Table2 represents the numerical characteristics and related indexes of the six layers. It can be seen from this table that for either temperature or wind velocity, the corresponding average value closely relates with the height. That is to say, the average value of temperature increases with the increasing height, whereas the average value of wind velocity decreases with the increasing height, indicating that among greenhouse environmental factors, the temperature and wind velocity are negatively correlated. Meanwhile, there is a negative correlation between the temperature and wind velocity distribution too. From Table 2, it is clear that J
T
=
0.1357, J
V
=
0.7461, and J
u
=
1.0. For Table 2, the overall index under the optimal scenario is 0.5528. All of these results under full ventilation conditions will be used as a reference to compare.For a visual expression, we choose the middle layer (height is 1.0 m) to display the cloud map which contains the temperature and wind velocity in Figure3.Cloud map of FONV at 1.0 m.
(a)
Temperature
(b)
Velocity ofx direction
(c)
Velocity ofz directionFigure3(a) shows the temperature cloud map with h
=
1.0
m. It can be seen from this figure that, for most part of the layer, the average temperature is around 295 K with some local variation. Figures 3(b) and 3(c) show the wind velocity contours along x and z directions, respectively. Clearly, the velocity along x direction is slower than that along the z direction; the reason could be explained by the direction of the greenhouse entrance. Meanwhile, two contours display significant turbulence phenomena, which are quite similar to Figure 3(a), indicating that temperature and wind velocity factors are coupled strongly.


## 5.2. Optimal Scenario

In this pattern, the optimization processes spend about 61 minutes to attain the control outputs (u
1
=
22.62
%, u
2
=
93.59
%, and u
3
=
10.95
%). The temperature-wind velocity cloud map at height of h
=
1.0
m is shown in Figure 4, and the NCs of the six layers are shown in Table 3.Table 3
NCs of optimal scenario.
h (m)
0.2
0.4
0.8
1.0
1.2
2.0
T (K)
Min
293.65
293.97
294.20
294.24
294.33
294.46
Ave
294.88
295.02
295.09
295.11
295.13
295.17
Max
296.33
296.53
296.21
296.13
296.09
296.52
V (m/s)
Min
0.001656
0.002452
0.007214
0.008335
0.010763
0.014362
Ave
0.270686
0.254144
0.233365
0.224360
0.215042
0.157519
Max
2.997552
3.063393
3.002424
2.955838
2.941845
3.117379
SDT (K)
0.32906
0.28879
0.24432
0.22823
0.21794
0.21176
SDV (m/s)
0.33999
0.33691
0.32172
0.31735
0.30061
0.22776Cloud map of scenario 1 at 1.0 m.
(a)
Temperature
(b)
Velocity ofx direction
(c)
Velocity ofz directionIt can be seen from the cloud map that the isotherms include 294.5 K, 295 K, and 295.5 K. Although the distribution is similarly compared to FONV, the distribution difference is reduced, and the double-eddy phenomenon could be found in the wind field that the amplitude is reduced as well. NCs of six different layers are calculated, including the average temperature, the average velocity, SDT, and SDV. Table3 has a similar regulation with Table 2.


## 5.3. Comparing and Conclusions

To make a comparison, we choose (12) to define the difference between FONV and optimal scenario as follows:(12)
Difference
=
Optimal
-
FONV
.Note. The horizontal axis is about the height of the layer, and the vertical layer is the average temperature, the average velocity, SDT, or SDV, respectively.Figure5 shows that both the average temperature and SDT are positive. Therefore, during the optimization process, the average temperature and the temperature distribution are both positive at each layer; these results show that the indexes of temperature are not better than the FONV. In other words, full ventilation mode represents the best ventilation form to reduce distribution for temperature, but only in temperature. Meanwhile, the average velocity was decreased because the u
1 is less than 100%; it led to the decrease of air exchange. The SDV is negative and indicates that the velocity distribution became better than FONV; these results can be explained by turbulent theory that the turbulence intensity has the same change direction with u. This phenomenon can also be conformed from an extreme case: when the greenhouse is closed the wind will not flow, so the ave
V
≈
0 and SDV
≈
0; this is the best way to reduce SDV but does not suit the greenhouse control. To further enrich the information differences between FONV and optimal scenario, we make Table 4 to describe it.Table 4
Compared result of indexes with FONV.
FONV
Optimal scenario
Value
Variation
J
T
0.1357
0.2457
+80.95%
J
V
0.7461
0.3154
−57.72%
J
u
1.0000
0.4196
−58.04%
J
0.5528
0.3084
−44.21%Figure 5
NCs difference between optimal scenario and FONV.In the optimal scenario, the values ofα
=
β
=
0.4 and γ
=
0.2 indicate that this mode belongs to a comprehensive optimization system, which includes the temperature, velocity, and control cost. Therefore, this scenario has a wide variety of applications in practice. Except for J
T, values of all other indexes are reduced, with the overall index reduced 44.21%, indicating that the optimization system in this paper has a great superiority for application in terms of greenhouse environment control.


## 6. Conclusions

By theoretical design and a specific calculation example, the following conclusions have been obtained.(1)
The optimization control method of greenhouse TF-WVF based on the CFD model is effective, which could extend the conventional “point” feedback forms to “field” feedback forms. Therefore, the new method will be a good strategy to control the distribution of greenhouse environmental factors with multi-index constrained.
(2)
The calculation example shows that the distribution of overall optimization index is reduced by 44.21%. The optimal scenario indicates that it is advantageous to explore the potentials of greenhouse control system when the CFD model is applied.
(3)
Among the NCs of six layers in the optimal scenario and FONV, the average temperature increases with the height, whereas the average velocity decreases with the height. Furthermore, both the temperature and velocity distributions are increased with the height; the result indicates that the height and distribution are negatively correlated.
(4)
In comparison with the FONV, the change of absolute value is different,Δ
J
V
=
0.431, Δ
J
T
=
0.110, so we can conclude that the WVF optimal has a more potential to dig than TF.A note must be explained. Although the single CFD cycle is set less thanT
CFD, the optimum course has a long time (61 minutes). In order to practice this way to control greenhouse environment, we must pay more attention to enhance the algorithm efficiency and adopt multicore parallel computing to reduce the time cost. If the optimal course can be limited in about 15 minutes, the CFD model for greenhouse control will be applied in practice.



---

*Source: 949128-2014-11-06.xml*