The problem of uneven ventilation in pig houses must be solved to effectively improve the winter environmental quality of pig nurseries in cold regions. In this study, the airflow field and airflow uneven coefficients of pig nurseries with duct ventilation were simulated and calculated using computational fluid dynamics, and compared with pig nurseries with different duct diameters, inlet and outlet air velocities and air supply angles. The average relative error between the simulated and measured values was
. Comparison of simulation results and airflow uneven coefficients showed that the airflow uneven coefficients of the fences were reduced, and the airflow field was uniformly distributed with a duct diameter of 0.3 m , inlet and outlet air velocities of 1.5 and
, respectively, and an air supply angle of
. These improvements resulted in a more homogeneous ventilation, which led to more uniform ventilation and contributed to discharging dirty air outdoors. Then the test pig nursery was modified based on the simulation and analysis results. Difference analyses were carried out between the control and the test pig nurseries. Comparative analyses showed the differences between the test data and the monitoring data were smaller, and the duct ventilation was more uniform, which was suitable for the healthy growth of piglets.
Received: 19/07/2024
Accepted: 10/10/2024
DOI
10.23967/j.rimni.2024.10.056321
Pig nursery
airflow field
numerical simulation
duct ventilation
Northeastern China is distinguished by its extensive land resource base and ample provision of feed resources, which are conducive to the flourishing growth of the swine farming sector [1]. However, because of the cold and long winter in this region, most pig nurseries have little or no ventilation in winter to avoid cold stress on pigs and reduce energy consumption [2-4]. If only thermal insulation and reduced ventilation are considered, harmful indoor gases (e.g., ammonia, hydrogen sulfide, and carbon dioxide) and excessive moisture and dust interfere with the healthy growth of pigs and largely
harm the health of breeders [5,6]. Therefore, innovative design of ventilation systems to guarantee the environmental quality of pig nurseries in cold regions is crucial for healthy breeding of pigs.
The engineering equipment and technology used in pig-raising facilities are always improving with the increasing demands from the global pig-raising industry. A single ventilation mode (e.g., longitudinal, horizontal, or vertical mode) can no longer meet the needs of healthy pig production. Hence, the basic theory of ventilation and the innovation of ventilation modes need to be studied further [7]. Mostafa et al. analytically tested the effects of two ventilation modes on the indoor environment of a large pig nursery house in Guangxi in winter [8]. The environmental quality of the nursery house with the underground duct air intake mode was better compared with the ceiling air intake mode [8]. In southern China, the cold air temperature in winter is generally not lower than
, so the two modes are more suitable for farm buildings, and the outside fresh air can directly enter the buildings. However, in Northeast China, due to the low winter temperature, tunnel ventilation will bring in abundant cold air and rapidly decrease the temperature inside the barn, which will induce cold stress on piglets and may cause ice formation in the faecal tract, so that faeces cannot be discharged normally. Thalfeldt et al. simulated and studied the ducted exhaust system of a nursery pig house based on the computational fluid dynamics (CFD) theory, and found this system effectively improved the environmental quality of the pig house [9]. Compared with the above ventilation modes, evenlyarranged air vents are more suitable for both cold and hot areas [10]. At the same time, the ducted ventilation model introducing fresh air into a sealed space in the form of supply air facilitates air heating, filtering, cooling, and other treatments. Additional ducts can be installed and ventilation openings be evenly distributed along the duct walls to ensure an even distribution of airflow within the space [11]. Moreover, exhaust air ventilation removes polluted air from the enclosed environment quickly and expels it outdoors. However, areas with poor airflow circulation or variations in airflow velocity are prone to dead zones or cross drafts [12]. The combined supply and exhaust air system can effectively integrate the advantages of both supply air and exhaust air ventilation, and thus improves ventilation efficiency and facilitates precise environment control in weaned pig nurseries [13]. This system combines the benefits of both supply air and exhaust air, allows for efficient air exchange and meets the specific requirements for maintaining the optimal environment for pigs. The duct ventilating mode better ensures consistent ventilation in a pig nursery and helps accurately control ventilation
.
The design of ventilation structures requires much field test data, but the pig nursery is restricted by a complicated environment and only a few monitoring points. At the same time, the environmental conditions cannot be stable and controllable for long periods, which leads to large margins of errors. Application and research with the CFD technology can accurately simulate the distribution laws of temperature and airflow fields in livestock and poultry houses, and theoretically underlie the optimization of ventilation structures [16-18]. Five two-way turbulence models were utilized to model a mechanically ventilated pig nursery and predict the airflow distribution in the pigsty. The choice of the turbulence model depended on the specific situation [19]. CFD-based simulation of airflow in pregnant sow barns under centralized ventilation revealed that the temperature field uniformity increased when the air outlet was located in the middle of the enclosure and that changes in the ventilation structure caused changes in the airflow field distribution in the pig nursery [20]. The effects of air inlet on airflow and temperature in a summer laying hen house were investigated using CFD, and the uniformity of airflow and temperature distributions increased with the location of air inlet and the distance from the hen house [21]. After the conditions of air inlet and outlet were optimized through CFD simulation, the environmental quality of the house was improved and the cold stress of pigs in winter was effectively alleviated [22]. CFD was used to study the influence of wind speed
on temperature fields in pigsties and the influence of pigsty structural changes on airflow fields, and the results theoretically underlay the optimization of an environmental control scheme [23]. Based on the above research methods, CFD was used to simulate the airflow field in a pig nursery under duct ventilation, and the simulation results were analyzed, which provide a reference basis for optimal ventilation system design and on-site renovation.
In this study, the factors affecting uniform air supply were analyzed, and a CFD model of a pig nursery under the duct ventilation mode was built on basis of field measurement data. The model was used to simulate the airflow field under different duct diameters, inlet and outlet air velocities, and air inlet angles, and analyze the influence of the ventilation structure on the airflow field. Then the main parameters of the ventilation structure were selected, and reference data were provided for ventilation structure optimization and on-site renovation. This study has a certain reference value for the construction of large-scale pig nurseries in cold areas to improve the pig nursery rearing environment and pig breeding quality.
The test-simulated pig nursery was located in Qiqihar, Heilongjiang Province, China (
N,
). In this region, the average maximum and minimum temperatures are
and
respectively in summer, and
and
respectively in winter. Fig. 1a shows the status of the pig nursery. The pig nursery was 17 m long, 10 m wide, and 3 m from the roof to the shed. There were 12 fences in two columns. Each fence was 3.5 m long and 2.5 m wide and housed 10 to 12 piglets. The fence floor was 0.5 m away from the cement floor of the pig nursery and was composed of the slatted floor and the heated cement floor. The gap width of the floor for the slatted floor was 15 mm . Fig. 1 b shows the three-dimensional effect diagram of the renovated pig nursery.
The neighboring pig nursery was selected as the control, which had the same structure and layout as the test pig nursery before the renovation. On-site measurements were conducted from December 2019 to January 2020. Sensor monitoring nodes were arranged in the same way as shown in Fig. 1. The monitoring data were extracted from the test and control pig nurseries at a
interval. Then groups of data were sampled continuously for a single monitoring data, and then the average value was recorded. To ensure the validity of data comparison and analysis results, the processes of heating, feeding, and manure cleaning during the monitoring period were kept the same in the two nurseries. The ventilation time of the two nurseries for duct ventilation and ventilation opening of the door and window on one side of the corridor is shown in Table 1. Real-time temperature, humidity,
, and
in the existing ventilation mode of the pig nursery were monitored. The monitoring equipment used at the site is listed in Table 2.
A total of sensor nodes were placed inside the pig nurseries at a height corresponding to the breathing level of the weaned pigs ( 0.2 m above the slatted floor). The environmental data obtained at this height can represent the quality of exhaled air by weaned pigs.
Figure 1: The pig nursery in the field experiment; (a) floor plan; (b) three-dimensional effect diagram of the renovated pig nursery (two lines of ducts for both air inlet and air outlet; six air inlets or outlets in each row)
Table 1: Opening time of duct ventilation and door/window ventilation
Ventilation method and
duration |
Duration of door
opening in two nurseries/(min) |
Duration of window
opening in two nurseries/(min) |
Duration of ducted
air exchange in the test pig nursery/(min) |
1 | 5 | - | |
10 | 20 | - | |
10 | 120 | - | |
- | - | 4 | |
5 | 30 | - | |
- | - | 6 | |
5 | 20 | - |
Table 2: Opening time of duct ventilation and door/window ventilation
Name | Type | Company | Range | Precision |
Humidity and
temperature measuring instrument |
Testo 625 | Testo Group, Black
Forest, Germany |
Temperature:
Humidity: |
|
Pump gas detector | HYE 2000 | Taizhou HanYao
Electronic Technology Co., Taizhou, China |
||
Fluid flow in nursery in the pig nursery obeys the laws of conservation of mass, momentum and energy and the law of composition. The governing equations are as follows [24]:
Mass conservation:
Energy conservation:
Momentum conservation:
where
is the fluid density,
, and
are the velocity vectors in the direction of
, and
, respectively,
is the heat transfer coefficient,
is the thermodynamic temperature,
is the specific heat capacity,
is the internal heat source,
is the pressure fluid microelement,
is the volume concentration of component
is the diffusion coefficient of component
.
The pig nursery has a symmetrical structure. To improve computational efficiency and save computer simulation time, a half-scale physical model of the nursery house was created and the rest of the house was symmetrical. For simplification of modeling, the effects of fences, slatted floor and piglets on the airflow in the house were ignored (Fig. 2). The air inlet was 0.9 m above the concrete floor, and the middle aisle was the air outlet and was 0.4 m from the ground.
Figure 2: Model of pig nursery
The air in the pig nursery was an incompressible ideal gas, and the fences, ducts, and walls were all non-slip wall boundaries. The walls on both sides were set as symmetric boundary conditions without friction or energy flow. The outlet and inlet of the air duct were arranged as the velocity inlet and outlet, respectively. The simulation was conducted under a steady-state. Bjerg et al. showed that CFD simulations of livestock buildings using a standard
turbulence model resulted in the smallest absolute and relative errors between simulation and monitoring results [25]. The standard
turbulence model was selected. Standard wall functions were used for simulations near the wall. The control equation was discretized using the finite volume method. The SIMPLEC algorithm with good convergence was used to solve the pressure velocity coupling equation [26]. The initial boundary conditions were listed in Table 3 [27]. The solution was assumed to converge when the sum of residuals for all cells in the computational domain was less than
for continuity and energy, and
for other variables such as
-,
-,
-velocity, k, and
.
Table 3: Boundary conditions for the CFD simulation model
Items | Boundary conditions |
Inlet | Velocity inlet |
Outlet | Velocity outlet |
Floor | Non-slip wall, |
North wall | Non-slip wall, |
South wall | Non-slip wall, |
Roof | Non-slip wall, |
The main factors affecting the ventilation are duct diameter
, air inlet angle
, air inlet velocity
, and air inlet velocity
[28,29]. The values of
, and
were changed to simulate the airflow field in the pig nursery while boundary conditions, initial conditions, and other relevant parameters were the same.
In general, the piglets were 0.4 m high, so the cross-section of the airflow field was selected as the backplane of the piglets. The suitable air velocity for a pig nursery is
, and uniform ventilation must be ensured. Airflow uniformity in pig nurseries was assessed using the airflow uneven coefficient [30]. With the fence as the unit, 10 points from each unit of F1-F6 were taken to read the air velocity evenly as the air dispersed, and the
of each case was calculated. Smaller
implies the ventilation in the pig nursery is more uniform. The airflow uneven coefficient of the cross-section at height
in the pig nursery
, and the average air velocity of the cross-section at height
are calculated as follows:
Failed to parse (unknown function "\mkern"): {\mathit{J}}_{\mathit{h}}=\frac{\sqrt{\frac{{\sum }_{\mathit{i}=1}^{\mathit{n}}\mkern1mu {\left( {\mathit{v}}_{\mathit{i}}-{\mathit{v}}_{\mathit{h}}\right) }^{2}}{\mathit{n}}}}{{\mathit{v}}_{\mathit{h}}} |
Failed to parse (unknown function "\mkern"): {\mathit{v}}_{\mathit{h}}=\frac{1}{\mathit{n}}{\sum }_{\mathit{i}=1}^{\mathit{n}}\mkern1mu {\mathit{v}}_{\mathit{i}} |
where
is the air velocity of monitoring point
is the number of monitoring points.
The efficiency of calculations is affected by the number of grids and can be improved by choosing the appropriate grid number. The results for certain characteristic physical quantities obtained with a particular grid size stabilize when the variability of the results due to an increase in grid number is acceptable. Unstructured grids with high flexibility are often used to simulate airflow [31].
Three different grids were used for independence tests, and the changing trends of wind speed were similar among the different grids, except that the wind speeds were different. The average wind speed of monitoring points 1-4 in F1 was taken as the reference. The maximum relative errors of coarse, normal and fine grids from the simulation results were
and
, respectively (Table 4). The wind speed fluctuation amplitude of a normal grid or a fine grid was very small, and did not largely impact the results. At the same time, computer performance and time cost were comprehensively considered during simulation and analysis with normal grid number.
Table 4: Grid independence validation
Number of grids | Average wind speed | Relative error (%) | |
Grid 1 (Coarse) | 538,527 | 0.1053 | 10.4 |
Grid 2 (Normal) | 0.1125 | 4.25 | |
Grid 3 (Fine) | 0.1128 | 3.96 |
Based on the simulation of airflow field, the 0.4 m high section in the living area of piglets was chosen. Four test points were selected for each fence to monitor air velocity, which was compared to the simulation data at the same location. The data curves are shown in Fig. 3. As the south fence was close to the window, external air entered the pig nursery through the gaps in the windows. As a result, the monitored air velocity and the simulated air velocity from monitoring points 21 to 24 were largely different. For the other 20 monitoring points, the relative errors calculated using the relative error formula were between 0 and 0.2 , and the average error was
. These results indicate the simulation results agree well with the monitoring data, and the simulations are effective.
When the air velocity was constant, the cross-sectional area of the duct was the main factor affecting the ventilation rate in unit time. At the same time, the duct diameter directly affected both the area and construction cost of the pig nursery. After comprehensive consideration, the ducts with diameters
and
were selected for simulation. The ducts were perpendicular to the horizontal plane, and the inlet and outlet air velocities were set as
and
, respectively.
Figure 3: Comparison of simulated and monitored air velocities
Fig. 4 reveals the changing rules of airflow field corresponding to different duct diameters. The air velocities of middle fences F3 and F4 were higher than those of other fences. The reason for this phenomenon is that the static pressure of the fluid inside the ventilation duct gradually increases along the flow direction. In contrast, the dynamic pressure decreases near the fence, resulting in higher outlet speeds at the rear end compared to the front end [32,33]. When the duct diameter increased from 0.2 to 0.4 m , the airflow field distribution was more uniform, but the maximum local air velocity reached 0.39
, which exceeded the ventilation requirement for piglets. At the same time, many weak ventilation areas existed in the pig nursery, and the air velocity on the side of slatted floor was generally low.
(a)
Figure 4: (Continued)
Figure 4: Airflow field at different duct diameters. (a)
The variation curves of airflow uneven coefficient
for each fence calculated from Eqs. (5) and (6) are shown in Fig. 5.
Figure 5: Effects of duct diameters on airflow uneven coefficient
The uneven coefficient
of fences F3 and F4 was higher (Fig. 5) owing to the higher local air velocities of fences F3 and F4. As the duct diameter rose from 0.2 to 0.4 m , the difference in air velocity
for each fence increased, resulting in an increment in
. When the duct diameter was
was generally smaller. Given the airflow field distribution, airflow uniformity, the space occupied by the duct, and the initial cost, we selected the 0.3 m duct as the inlet and outlet ventilation duct.
Negative pressure ventilation was adopted in the pig nursery, and the duct diameter was set at 0.3 m . The situations
, and
were simulated, and the air inlet and outlet velocities (
) were set as 0.5 and
and
and
, respectively. The simulation results at
and
are shown in Fig. 4 b . Other simulation results are shown in Fig. 6.
At
and
, many windless areas appeared, and the air velocity on the side of the slatted floor was too low (Fig. 6a). The air inlet and outlet velocities as well as the local air velocities on the concrete floor increased. At
and
, the local air velocity maximized to
and the air velocity on the side of the slatted floor also increased (Fig. 6b), but the weak ventilation on the side of the slatted floor was not fixed.
Figure 6: Airflow fields at different air inlet and outlet velocities. (a)
and
and
The variation curves of airflow uneven coefficient
for each fence calculated from Eqs. (5) and (6) are shown in Fig. 7.
Figure 7: Effects of air velocities on the airflow uneven coefficient
When the air velocity increased, the airflow uneven coefficient
decreased. However, at
and
, the local air velocity difference for fence F 4 was large, resulting in an increase of
in this area. Given the airflow field distribution and airflow uniformity analysis results, the ventilation requirements of piglets were better met at
and
.
According to the analysis results in Sections 3.1 and 3.2, there were problems of excessive local air velocity in fences F3 and F4 and weak ventilation on the side of the slatted floor, regardless of the changes to duct diameter and air inlet and outlet velocities. The air inlet angle also changed the airflow field distribution. To improve this problem, the air inlet angle was changed. The air inlet angles at
and
were simulated when the duct diameter was 0.3 m , the air inlet and outlet velocities were
and
, the vertical ground angle was
, and the inlet of the duct was inclined to the side of the slatted floor. The simulation results are shown in Fig. 8.
(a)
Figure 8: (Continued)
Figure 8: Airflow field at different air inlet angles. (a)
Comparison of the left side of the colorimetric map showed as the ventilation angle decreased, the area of maximum air velocity in the airflow field moved from the cement floor to the slatted floor, and the maximum air velocity for local ventilation decreased from 0.42 to
. When the ventilation angle was
, many weakly ventilated areas appeared. At a ventilation angle of
, the ventilation was more uniform, and a few areas with weak ventilation were formed.
The variation curves of airflow uneven coefficient
for each fence calculated as per Eqs. (5) and (6) are shown in Fig. 9.
Figure 9: Effects of air inlet angles on the airflow uneven coefficient
When the air inlet angle was
, the
of the other fences was smaller, except for F3 and F4 in which
was large because of the large local air velocity difference of the slatted floor. Nevertheless, the overall air velocity met the requirements of the piglets, and the air velocity on the side of the slatted floor increased, which was more conducive to the discharge of dirty gas. Given the reduction of the weak ventilation areas, an air inlet angle of
is more suitable for raising piglets.
Based on the numerical simulation results from the airflow field of the pig nursery under various conditions, the experimental pig nursery was redesigned by considering the variation in airflow uneven coefficient curves. The redesigned pig nursery is shown in Fig. 10. Each fence was provided with an air inlet right above the concrete floor, and an air outlet in the middle aisle.
Figure 10: The redesigned nursery
The main parameters are as follows: the diameters of both the inlet and outlet ducts were 0.3 m ; the air velocities of the fans in the inlet and outlet ducts were 1.5 and
, respectively; the air inlet was equipped with a hood to adjust the ventilation angle to
.
The environmental conditions in the pig nursery from 6:00 on 19 January 2020 to 6:00 the next day were monitored. The outdoor temperature maximized to
during the day, and minimized to
at night. The monitored parameters included temperature, relative humidity, and
and
concentrations. Data were sampled every 15 min , and 10 sets of data for each monitoring session were continuously sampled. The average values were then recorded. Some of the data are shown in Table 2.
The significance of differences between the monitoring data at each node was calculated to further analyze whether the airflow in the pig nursery was uniform under the duct ventilation mode. The variability of temperature, humidity, and
and
concentrations among the nodes was tested on SPSS. The results are shown in Table 5.
Because node A was closer to the door and window and was more affected by the outdoor air, the data at node A were significantly different from the data at other nodes in both the control and test pig nurseries. The main reason for the consistency in the difference of
concentrations at different nodes in the two nurseries was that the amount and moment of fecal discharge were more random, resulting in large changes in
concentrations at a certain time point.
Table 5: Monitoring data and difference analysis
A | B | C | D | E | F | ||
Control pig nursery | Temperature | ||||||
Humidity (%) | |||||||
) | |||||||
) | |||||||
Test pig nursery | Temperature ( ) | ||||||
Humidity (%) | |||||||
Note: Differences between different letters in the same row are significant (
) and differences between the same mothers are not significant (
).
In the control pig nursery, the data of neighboring nodes were more significantly different. For instance, the differences in temperature and humidity between nodes B and C were not significant (
), but the differences in
and
concentrations were significant (
). In contrast, the data of neighboring nodes in the test pig nursery were less significantly different, such as the insignificant differences in temperature, humidity, and
concentration between nodes B and C (
). This result indicates the test barn has better ventilation uniformity and is more suitable for piglet growth.
The duct-ventilated system of a pig nursery in a cold region was studied. CFD was used to simulate the airflow field of the pig nursery with different duct diameters, air velocities, and air inlet angles. The ventilation state in each case was evaluated and analyzed using the uneven coefficient of airflow field. The accuracy of the simulation results was verified through field measurement.
(1) The maximum relative error between the simulated and measured data is
, and the average relative error is
. The simulation results well agree with the actual environment, and thus can be used as a basis for optimizing the ventilation system design. Based on the CFD simulation and data analysis, the duct diameter was determined to be 0.3 m , the air inlet and outlet velocities were 1.5 and
, respectively, and the ventilation angle was
.
(2) The ventilation system in the test pig nursery was modified, and the monitoring data in the control and test nurseries were compared. The data at neighboring nodes of the control pig nursery were significantly more different. The test pig nursery had significantly less difference in the data of neighboring nodes. Hence, the modified test pig nursery had better ventilation uniformity and was more suitable for the healthy growth of piglets.
The pig nursery ventilation system was optimally designed, but the control strategy still needs to be improved. To better control the pig nursery environment in the future, the coupling mechanism of the environment inside the pig house can be investigated, and the environmental control strategy related to the environment outside the pig nursery can be established. The seasonality of the pig nursery environment can also be analyzed to achieve autonomous and precise control of the pig house environment in different seasons.
Acknowledgement: The authors would like to thank the anonymous reviewers for their constructive comments, which helped improve the quality of this paper.
Funding Statement: This work was supported by Collaborative Innovation Achievement Project of "Double First-Class" Disciplines in Heilongjiang Province, grant number LJGXCG2024-P25, Key Projects of Qiqihar City Scientific and Technological Plan, grant number ZDGG-202202, the Basic Research Fund for State-Owned Universities in Heilongjiang Province, grant number 145309319, General Research Project on Higher Education Teaching Reform in Heilongjiang Province, grant number SJGY20220410, and Educational Science Research Project of Qiqihar University, grant number GJQTYB202212.
Author Contributions: Zhidong Wu: Writing review, project administration; Kaixiang Xu: Writing original draft; Yanwei Chen: Formal analysis; Yonglan Liu: Validation; Meiqi Liu: Conceptualization. All authors reviewed the results and approved the final version of the manuscript.
Availability of Data and Materials: The experimental data are original and not convenient for public disclosure. If there is a reasonable request, please contact the authors.
Ethics Approval: Not applicable.
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
Published on 13/12/24
Licence: CC BY-NC-SA license
Are you one of the authors of this document?