modelNaturalVentilationRate_1

Air exchange rate from the greenhouse air to the outside air function of wind and temperature by De Jong (1990)

Information

The natural ventilation rate due to roof ventilation is described by De Jong (1990):

f_vent = 1/2*fr_window*sqrt(f_wind^2+f_temp^2) [m^3*s^(-1)*m^(-2)]

where

  • fr_window (m^(-2)) is the number of windows per square meter of greenhouse floor,
  • f_wind (m^3*s^(-1)*window^(-1)) is the wind driven ventilation,
  • f_temp (m^3*s^(-1)*window^(-1)) is the temperature driven ventilation.


The contribution of the temperature driven vetilation in the total ventilation is small but can be important during nighttime and in winter. It is described by:

f_temp = C_f*l/3*sqrt(abs(g*beta*Delta_T))*H^1.5 [m^3*s^(-1)*window^(-1)]

where

  • C_f (-) is the discharge coefficient caused by friction in the window opening,
  • l (m) is the length of the window
  • () is the gravitational acceleration,
  • beta (K^(-1)) is the thermal expansion coefficient,
  • Delta_T (K) is the temperature difference of the greenhouse air and outside air,
  • H (m) is the height of the front opening of the window (H = h*(sin(psi) - sin(psi-theta)) with h the window width, psi the roof slope and theta the opening angle of the window).

The wind speed driven ventilation is computed differently for vents in the windward side and the leeside side. De Jong used 'Window Functions' to relate the air exhange with the wind speed and the opening of a window. The air exchange is described by:

f_wind = (2.29e-2*(1-exp(-theta_l/21.1)) + 1.2e-3*theta_w*exp(theta_w/211))*A_0*u_wind [m^3*s^(-1)*window^(-1)]

where

  • theta_l (rad) is the opening angle of the leeside ventilators,
  • theta_w (rad) is the opening angle of the windward ventilators,
  • A_0 (m^2) is the surface of a window (A_0 = l*h),
  • () is the wind speed.

Parameters

TypeNameDefaultDescription
Modelica.Units.SI.Anglephi25/180*Modelica.Constants.piinclination of the roof (0 if horizontal, 25 for typical cover)
Realfr_window0.078number of windows per m2 greenhouse
Modelica.Units.SI.Lengthllength of the window
Modelica.Units.SI.Lengthhwidth of the window
BooleanthermalScreenfalsepresence of a thermal screen in the greenhouse

Components

TypeNameDefaultDescription
Modelica.Units.SI.Velocityu0Wind speed
Modelica.Units.SI.Angletheta_l0window opening at the leeside side
Modelica.Units.SI.Angletheta_w0window opening at the windward side
Modelica.Units.SI.TemperatureDifferencedT10
Realf_vent
Realf_vent_wind
Realf_vent_temp
RealC_fconstant of discharge of energy caused by friction in the window opening
Realbetathermal expansion coefficient
Realf_vent_topAir exchange rate at the top air compartment in presence of a thermal screen
Realf_vent_airAir exchange rate at the main air compartment