Analysis of Airflow Uniformity in Pig Nurseries Using Duct Ventilation in Northeast China_prova_word

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-8e0f0d51e5c8e8e534a6496110217c84de6ede65-c.jpg|600px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Analysis of Airflow Uniformity in Pig Nurseries Using Duct Ventilation in Northeast China'''</big></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Zhidong Wu¹,2,3,*, Kaixiang Xu¹, Yanwei Chen¹, Yonglan Liu1 <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1}</math>
|}
</div> and Meiqi Liu1</div>

1 School of Mechanical and Electrical Engineering, Qiqihar University, Qiqihar, 161006, China<br/>2 The Engineering Technology Research Center for Precision Manufacturing Equipment and Industrial Perception of Heilongjiang Province, Qiqihar, 161006, China<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{3}</math>
|}
Institute of Economic Crops, Heilongjiang Academy of Agricultural Sciences, Harbin, 150086, China<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-bbd9e5c95b2f3d860a0301bdb798dd03aefaab76-c.jpg|600px]] </div>

==INFORMATION==

Keywords:<br/>Pig nursery<br/>airflow field<br/>numerical simulation<br/>duct ventilation<br/>DOI: 10.23967/j.rimni.2024.10.056321

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Analysis of Airflow Uniformity in Pig Nurseries Using Duct Ventilation in Northeast China'''</big></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Zhidong Wu <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1,2,3,}</math>
|}
</div>, Kaixiang Xu <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1}</math>
|}
</div>, Yanwei Chen <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1}</math>
|}
</div>, Yonglan Liu <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1}</math>
|}
</div> and Meiqi Liu</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{1}</math>
|}
</div> School of Mechanical and Electrical Engineering, Qiqihar University, Qiqihar, 161006, China</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{2}</math>
|}
</div> The Engineering Technology Research Center for Precision Manufacturing Equipment and Industrial Perception of Heilongjiang Province, Qiqihar, 161006, China</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\, }^{3}</math>
|}
</div> Institute of Economic Crops, Heilongjiang Academy of Agricultural Sciences, Harbin, 150086, China</div>
-->

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Abstract'''</big></div>

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>12\%</math>
|}
. 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>2.0\, m/s</math>
|}
, respectively, and an air supply angle of
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
. 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.

==OPEN ACCESS==

Received: 19/07/2024<br/>Accepted: 10/10/2024<br/>DOI<br/>10.23967/j.rimni.2024.10.056321

==Keywords:==

Pig nursery<br/>airflow field<br/>numerical simulation<br/>duct ventilation

=1 Introduction=

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{0}^{\circ }C</math>
|}
, 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>[14,15]</math>
|}
.

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<br/>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.

=2 Materials and Methods=

==2.1 Test Pig Nursery==

The test-simulated pig nursery was located in Qiqihar, Heilongjiang Province, China (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{47}^{\circ }{44}^{\mathit{{'}}}</math>
|}
N,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{124}^{\circ }{04}^{\mathit{{'}}}E</math>
|}
). In this region, the average maximum and minimum temperatures are
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{27.3}^{\circ }C</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{18}^{\circ }C</math>
|}
respectively in summer, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-{10.7}^{\circ }C</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-{21.7}^{\circ }C</math>
|}
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.

==2.2 Model Parameter Measurement==

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>15-min</math>
|}
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,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{CO}_{2}</math>
|}
, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
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.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-5fa70be81a37804fc13db5846f90522715868f12-c.jpg|600px]] </div>

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

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Ventilation method and

duration
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Duration of door

opening in two

nurseries/(min)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Duration of window

opening in two

nurseries/(min)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Duration of ducted

air exchange in the

test pig nursery/(min)
|-
| <math>07:35</math>
| 1
| 5
| -
|-
| <math>10:00</math>
| 10
| 20
| -
|-
| <math>12:00</math>
| 10
| 120
| -
|-
| <math>12:35</math>
| -
| -
| 4
|-
| <math>16:00</math>
| 5
| 30
| -
|-
| <math>17:54</math>
| -
| -
| 6
|-
|  style="border-bottom: 1pt solid black;"|<math>22:40</math>
|  style="border-bottom: 1pt solid black;"|5
|  style="border-bottom: 1pt solid black;"|20
|  style="border-bottom: 1pt solid black;"|-
|}


Table 2: Opening time of duct ventilation and door/window ventilation

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Name
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Type
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Company
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Range
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Precision
|-
| Humidity and

temperature measuring

instrument
| Testo 625
| Testo Group, Black

Forest, Germany
| Temperature:

<math>-{10}^{\circ }C-{60}^{\circ }C</math>

Humidity: <math>0\%\sim 100\%</math>
| <math>\pm {0.5}^{\circ }C</math>
|-
| Pump gas detector
|  rowspan='1'|HYE 2000
| Taizhou HanYao

Electronic Technology

Co., Taizhou, China
| <math>{CO}_{2}:0-5398mg/{m}^{3}</math>
| 
|-
| <math>{NH}_{3}:0-69mg/{m}^{3}</math>
| 
| 
| 
|}

==2.3 Governing Equations==

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:<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\frac{\mathit{\partial }\mathit{\rho }\mathit{u}}{\mathit{\partial }\mathit{x}}+</math><math>\frac{\mathit{\partial }\mathit{\rho }\mathit{v}}{\mathit{\partial }\mathit{y}}+</math><math>\frac{\mathit{\rho }\mathit{w}}{\mathit{\partial }\mathit{z}}=0</math>
|}
<br/>Energy conservation:<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\frac{\mathit{\partial }(\mathit{\rho }\mathit{T})}{\mathit{\partial }\mathit{t}}+</math><math>\frac{\mathit{\partial }(\mathit{\rho }\mathit{u}\mathit{T})}{\mathit{\partial }\mathit{x}}+</math><math>\frac{\mathit{\partial }(\mathit{\rho }\mathit{v}\mathit{T})}{\mathit{\partial }\mathit{y}}+</math><math>\frac{\mathit{\partial }(\mathit{\rho }\mathit{w}\mathit{T})}{\mathit{\partial }\mathit{z}}</math>
|}
<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>=\frac{\mathit{\partial }}{\mathit{\partial }\mathit{x}}\left( \frac{\mathit{k}}{{\mathit{c}}_{\mathit{p}}}\frac{\mathit{\partial }\mathit{T}}{\mathit{\partial }\mathit{x}}\right) +</math><math>\frac{\mathit{\partial }}{\mathit{\partial }\mathit{y}}\left( \frac{\mathit{k}}{{\mathit{c}}_{\mathit{p}}}\frac{\mathit{\partial }\mathit{T}}{\mathit{\partial }\mathit{y}}\right) +</math><math>\frac{\mathit{\partial }}{\mathit{\partial }\mathit{z}}\left( \frac{\mathit{k}}{{\mathit{c}}_{\mathit{p}}}\frac{\mathit{\partial }\mathit{T}}{\mathit{\partial }\mathit{z}}\right) +</math><math>{\mathit{S}}_{\mathit{T}}</math>
|}


Momentum conservation:
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\left\{ \begin{matrix}\frac{\mathit{\partial }(\mathit{\rho }\mathit{u})}{\mathit{\partial }\mathit{t}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{u}\mathit{u})}{\mathit{\partial }\mathit{x}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{u}\mathit{v})}{\mathit{\partial }\mathit{y}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{u}\mathit{w})}{\mathit{\partial }\mathit{z}}\\=\frac{\mathit{\partial }}{\mathit{\partial }\mathit{x}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{u}}{\mathit{\partial }\mathit{x}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{y}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{u}}{\mathit{\partial }\mathit{y}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{z}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{u}}{\mathit{\partial }\mathit{z}}\right) -\frac{\mathit{\partial }\mathit{p}}{\mathit{\partial }\mathit{x}}\\\frac{\mathit{\partial }(\mathit{\rho }\mathit{v})}{\mathit{\partial }\mathit{t}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{v}\mathit{u})}{\mathit{\partial }\mathit{x}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{v}\mathit{v})}{\mathit{\partial }\mathit{y}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{v}\mathit{\omega })}{\mathit{\partial }\mathit{z}}\\=\frac{\mathit{\partial }}{\mathit{\partial }\mathit{x}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{v}}{\mathit{\partial }\mathit{x}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{y}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{v}}{\mathit{\partial }\mathit{y}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{z}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{v}}{\mathit{\partial }\mathit{z}}\right) -\frac{\mathit{\partial }\mathit{p}}{\mathit{\partial }\mathit{y}}\\\frac{\mathit{\partial }(\mathit{\rho }\mathit{w})}{\mathit{\partial }\mathit{t}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{w}\mathit{u})}{\mathit{\partial }\mathit{x}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{w}\mathit{v})}{\mathit{\partial }\mathit{y}}+\frac{\mathit{\partial }(\mathit{\rho }\mathit{w}\mathit{w})}{\mathit{\partial }\mathit{z}}\\=\frac{\mathit{\partial }}{\mathit{\partial }\mathit{x}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{w}}{\mathit{\partial }\mathit{x}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{y}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{w}}{\mathit{\partial }\mathit{y}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{z}}\left( \mathit{\mu }\frac{\mathit{\partial }\mathit{w}}{\mathit{\partial }\mathit{z}}\right) -\frac{\mathit{\partial }\mathit{p}}{\mathit{\partial }\mathit{z}}-\mathit{\rho }\mathit{g}\end{matrix}\right.</math> 
|}

==Conservation of composition:==
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\begin{matrix}&\frac{\mathit{\partial }\left( \mathit{\rho }{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{t}}+\frac{\mathit{\partial }\left( \mathit{\rho }\mathit{u}{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{x}}+\frac{\mathit{\partial }\left( \mathit{\rho }\mathit{v}{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{y}}+\frac{\mathit{\partial }\left( \mathit{\rho }\mathit{w}{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{z}}\\&\mathit{\, }=\frac{\mathit{\partial }}{\mathit{\partial }\mathit{x}}\left( {\mathit{D}}_{\mathit{s}}\frac{\mathit{\partial }\left( \mathit{\rho }{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{x}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{y}}\left( {\mathit{D}}_{\mathit{s}}\frac{\mathit{\partial }\left( \mathit{\rho }{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{y}}\right) +\frac{\mathit{\partial }}{\mathit{\partial }\mathit{z}}\left( {\mathit{D}}_{\mathit{s}}\frac{\mathit{\partial }\left( \mathit{\rho }{\mathit{c}}_{\mathit{s}}\right) }{\mathit{\partial }\mathit{z}}\right) \end{matrix}</math>
|}


where
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{\rho }</math>
|}
is the fluid density,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>kg/{m}^{3};\mathit{u},\mathit{v}</math>
|}
, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{w}</math>
|}
are the velocity vectors in the direction of
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{x},\mathit{y}</math>
|}
, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{z}</math>
|}
, respectively,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>m/s;\mathit{k}</math>
|}
is the heat transfer coefficient,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>W/(m\cdot K),\mathit{T}</math>
|}
is the thermodynamic temperature,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>K;{\mathit{C}}_{\mathit{P}}</math>
|}
is the specific heat capacity,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>J/(kg\cdot K),{\mathit{S}}_{\mathit{T}}</math>
|}
is the internal heat source,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>W;\mathit{P}</math>
|}
is the pressure fluid microelement,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>Pa;{\mathit{C}}_{s}</math>
|}
is the volume concentration of component
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{s},\, kg/kg;{\mathit{D}}_{\mathit{s}}</math>
|}
is the diffusion coefficient of component
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{s},{\, m}^{2}/s</math>
|}
.

==2.4 Pig Nursery Model==

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.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-fd9e48b80a43a7b3fb37313136c58a0a09720f6f-c.jpg|600px]] </div>

Figure 2: Model of pig nursery

==2.5 Boundary Conditions and Numerical Calculation==

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>k-\mathit{\epsilon }</math>
|}
turbulence model resulted in the smallest absolute and relative errors between simulation and monitoring results [25]. The standard
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>k-\mathit{\epsilon }</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>1\times {10}^{-6}</math>
|}
for continuity and energy, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>1\times</math> 
|}
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{10}^{-3}</math>
|}
for other variables such as
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{X}</math>
|}
-,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{Y}</math>
|}
-,
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{Z}</math>
|}
-velocity, k, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{\epsilon }</math>
|}
.

Table 3: Boundary conditions for the CFD simulation model

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Items
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Boundary conditions
|-
| Inlet
| Velocity inlet <math display="inline">\left( 0.5\, m/s,{20}^{\circ }C\right)</math> 
|-
| Outlet
| Velocity outlet <math display="inline">(1\, m/s)</math>
|-
| Floor
| Non-slip wall, <math display="inline">{12}^{\circ }C</math>
|-
| North wall
| Non-slip wall, <math display="inline">{18}^{\circ }C</math>
|-
| South wall
| Non-slip wall, <math display="inline">{15}^{\circ }C</math>
|-
|  style="border-bottom: 1pt solid black;"|Roof
|  style="border-bottom: 1pt solid black;"|Non-slip wall, <math display="inline">{16}^{\circ }C</math>
|}

==2.6 Evaluation Indicators==

The main factors affecting the ventilation are duct diameter
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>(\mathit{D})</math>
|}
, air inlet angle
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>(\mathit{A})</math>
|}
, air inlet velocity
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\left( {\mathit{v}}_{in~}\right)</math> 
|}
, and air inlet velocity
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\left( {\mathit{v}}_{\mathit{o}\mathit{u}\mathit{t}}\right)</math> 
|}
[28,29]. The values of
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{D},\mathit{A}</math>
|}
, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{v}</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>0.2-0.3\, m/s</math>
|}
, 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{\mathit{h}}</math>
|}
of each case was calculated. Smaller
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{\mathit{h}}</math>
|}
implies the ventilation in the pig nursery is more uniform. The airflow uneven coefficient of the cross-section at height
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{h}</math>
|}
in the pig nursery
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\left( {\mathit{J}}_{\mathit{h}}\right)</math> 
|}
, and the average air velocity of the cross-section at height
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{h}\left( {\mathit{v}}_{\mathit{h}},\, m/s\right)</math> 
|}
are calculated as follows:<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\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}}}</math>
|}
<br/>
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{h}}=\frac{1}{\mathit{n}}{\sum }_{\mathit{i}=1}^{\mathit{n}}\mkern1mu {\mathit{v}}_{\mathit{i}}</math>
|}
<br/>where
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{i}}</math>
|}
is the air velocity of monitoring point
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{i},\, m/s;\mathit{n}</math>
|}
is the number of monitoring points.

=3 Results and Discussion=

==3.1 Model Validation==

===3.1.1 Grid Independence Validation===

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>10.4\%,4.25\%</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>3.96\%</math>
|}
, 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

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Number of grids
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Average wind speed <math display="inline">\left( m\cdot {s}^{-1}\right)</math> 
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;"|Relative error (%)
|-
| Grid 1 (Coarse)
| 538,527
| 0.1053
| 10.4
|-
| Grid 2 (Normal)
| <math>1,463,425</math>
| 0.1125
| 4.25
|-
|  style="border-bottom: 1pt solid black;"|Grid 3 (Fine)
|  style="border-bottom: 1pt solid black;"|<math>2,325,956</math>
|  style="border-bottom: 1pt solid black;"|0.1128
|  style="border-bottom: 1pt solid black;"|3.96
|}

===3.1.2 CFD Model Validation===

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>11.79\%</math>
|}
. These results indicate the simulation results agree well with the monitoring data, and the simulations are effective.

==3.2 Influence of Duct Diameter on Simulation Results==

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{D}}_{1}=0.2\, m,{\mathit{D}}_{2}=0.3\, m</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{D}}_{3}=0.4\, m</math>
|}
were selected for simulation. The ducts were perpendicular to the horizontal plane, and the inlet and outlet air velocities were set as
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{i}\mathit{n}}=1.0\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-1.5\, m/s</math>
|}
, respectively.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-e0d5a79be9c526c949ccdd967cb3242344d4ec1b-c.jpg|600px]] </div>

Figure 3: Comparison of simulated and monitored air velocities<br/>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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>m/s</math>
|}
, 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.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-43042191218514b21ee9161e472f60d8faa1d663-c.jpg|600px]] </div>

<br/>(a)

Figure 4: (Continued)<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-e6f5818836789bf2b3825d84812c09ce87bc9f5a-c.jpg|600px]] </div>

Figure 4: Airflow field at different duct diameters. (a)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{D}}_{1}=0.2\, m</math>
|}
; (b)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{D}}_{2}=0.3\, m</math>
|}
; (c)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{D}}_{1}=0.4\, m</math>
|}
<br/>The variation curves of airflow uneven coefficient
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
for each fence calculated from Eqs. (5) and (6) are shown in Fig. 5.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-a319816ce99f76376e1e621c8376d846cf5e3f22-c.jpg|600px]] </div>

Figure 5: Effects of duct diameters on airflow uneven coefficient

The uneven coefficient
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
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

==SCIPEDIA==

for each fence increased, resulting in an increment in
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
. When the duct diameter was
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>0.3\, m,{\mathit{J}}_{0.4}</math>
|}
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.

==3.3 Influence of Wind Velocity on Simulation Results==

Negative pressure ventilation was adopted in the pig nursery, and the duct diameter was set at 0.3 m . The situations
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{V}1,\mathit{V}2</math>
|}
, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{V}3</math>
|}
were simulated, and the air inlet and outlet velocities (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{i}\mathit{n}},{\mathit{v}}_{out~}</math>
|}
) were set as 0.5 and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-1.0\, m/s;1.0</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-1.5\, m/s;1.5</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-2.0\, m/s</math>
|}
, respectively. The simulation results at
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=1.0\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-1.5\, m/s</math>
|}
are shown in Fig. 4 b . Other simulation results are shown in Fig. 6.

At
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=0.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-1.0\, m/s</math>
|}
, 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=1.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-2.0\, m/s</math>
|}
, the local air velocity maximized to
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>0.4\, m/s</math>
|}
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.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-bb5bd33ccf2be09b2d2e06da3dcb3feb45a138d9-c.jpg|600px]] </div>

Figure 6: Airflow fields at different air inlet and outlet velocities. (a)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=0.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-1.0\, m/s</math>
|}
; (b)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=1.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-2.0\, m/s</math>
|}


The variation curves of airflow uneven coefficient
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
for each fence calculated from Eqs. (5) and (6) are shown in Fig. 7.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-829c44aec321dae60f725ad792d838200ac0eda4-c.jpg|600px]] </div>

Figure 7: Effects of air velocities on the airflow uneven coefficient<br/>When the air velocity increased, the airflow uneven coefficient
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
decreased. However, at
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{i}\mathit{n}}=</math>
|}
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>1.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-2.0\, m/s</math>
|}
, the local air velocity difference for fence F 4 was large, resulting in an increase of
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
in this area. Given the airflow field distribution and airflow uniformity analysis results, the ventilation requirements of piglets were better met at
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{\mathit{i}\mathit{n}}=1.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-2.0\, m/s</math>
|}
.

==3.4 Influence of Ventilation Inlet Angle on Simulation Results==

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{60}^{\circ },{45}^{\circ }</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{30}^{\circ }</math>
|}
were simulated when the duct diameter was 0.3 m , the air inlet and outlet velocities were
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{in~}=1.5\, m/s</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{v}}_{out~}=-2.0\, m/s</math>
|}
, the vertical ground angle was
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{90}^{\circ }</math>
|}
, and the inlet of the duct was inclined to the side of the slatted floor. The simulation results are shown in Fig. 8.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-ce483029d1a13827692d1206b9ecb17876b767cd-c.jpg|600px]] </div>

<br/>(a)

Figure 8: (Continued)<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-8461d4da6e25c9e96433b906f7dc6120db1f733c-c.jpg|600px]] </div>

Figure 8: Airflow field at different air inlet angles. (a)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{A}}_{1}={60}^{\circ }</math>
|}
; (b)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{A}}_{1}={45}^{\circ }</math>
|}
; (c)
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{A}}_{1}={30}^{\circ }</math>
|}
<br/>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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>0.21\, m/s</math>
|}
. When the ventilation angle was
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{30}^{\circ }</math>
|}
, many weakly ventilated areas appeared. At a ventilation angle of
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
, the ventilation was more uniform, and a few areas with weak ventilation were formed.

The variation curves of airflow uneven coefficient
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
for each fence calculated as per Eqs. (5) and (6) are shown in Fig. 9.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-d2eafddaca4ff8a2354abb975e087ea806a82fb5-c.jpg|600px]] </div>

Figure 9: Effects of air inlet angles on the airflow uneven coefficient

When the air inlet angle was
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
, the
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
of the other fences was smaller, except for F3 and F4 in which
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{\mathit{J}}_{0.4}</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
is more suitable for raising piglets.

==3.5 Structural Optimization and Testing==

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.<br/>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Sanchez Pinedo_937639649-image-27eab0abec456cd1f7926e3a7ac28fc604794353-c.jpg|600px]] </div>

Figure 10: The redesigned nursery<br/>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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>2.0\, m/s</math>
|}
, respectively; the air inlet was equipped with a hood to adjust the ventilation angle to
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
.

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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-{15}^{\circ }C</math>
|}
during the day, and minimized to
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>-{28}^{\circ }C</math>
|}
at night. The monitored parameters included temperature, relative humidity, and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{CO}_{2}</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{CO}_{2}</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
concentrations at a certain time point.

Table 5: Monitoring data and difference analysis

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|A
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|B
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|C
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|D
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|E
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|F
|-
|  rowspan='1' style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Control pig nursery
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Temperature <math display="inline">\left( {\, }^{\circ }C\right)</math> 
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>24.00\pm {0.66}^{c~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>24.91\pm {0.48}^{ac~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>24.96\pm {0.51}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>24.91\pm {0.46}^{ac}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>24.60\pm {0.48}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.68\pm {0.53}^{d~}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Humidity (%)
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>78.19\pm {2.47}^{c~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>88.46\pm {0.85}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>85.53\pm {1.36}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>90.24\pm {0.80}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>90.35\pm {0.66}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>88.32\pm {0.88}^{\mathit{b}}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math display="inline">{CO}_{2}\left( mg/{m}^{3}\right.</math>  )
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1660.04\pm \\&{124.82}^{d}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1746.71\pm \\&{111.27}^{c}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1784.72\pm \\&{128.58}^{b}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1817.06\pm \\&{134.12}^{ab}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1834.63\pm \\&{139.06}^{a}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1754.94\pm \\&{119.27}^{b}\end{matrix}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math display="inline">{NH}_{3}\left( mg/{m}^{3}\right.</math>  )
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.42\pm {0.51}^{C}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.59\pm {0.48}^{b}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.75\pm {0.44}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.74\pm {0.49}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.79\pm {0.42}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.72\pm {0.50}^{ab}</math>
|-
|  rowspan='1' style="border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Test pig nursery
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Temperature ( <math display="inline">{\, }^{\circ }C</math> )
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>22.91\pm {0.80}^{c}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.47\pm {0.48}^{b}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.51\pm {0.51}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.75\pm {0.45}^{a~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.91\pm {0.51}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>23.41\pm {0.51}^{b~}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|Humidity (%)
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>53.21\pm {3.05}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>64.56\pm {2.21}^{a~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>64.62\pm {2.45}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>64.74\pm {2.41}^{a~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>64.87\pm {2.45}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>64.56\pm {2.42}^{a~}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>{CO}_{2}\left( mg/{m}^{3}\right)</math> 
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>941.93\pm {36.31}^{b~}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1018.34\pm \\&{41.11}^{a}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1024.73\pm \\&{41.08}^{a}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1023.08\pm \\&{39.05}^{a}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1023.53\pm \\&{40.24}^{a}\end{matrix}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>\begin{matrix}&1019.15\pm \\&{40.21}^{a}\end{matrix}</math>
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>{NH}_{3}\left( mg/{m}^{3}\right)</math> 
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>1.10\pm {0.28}^{C}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.10\pm {0.46}^{b}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.22\pm {0.45}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.21\pm {0.41}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.27\pm {0.39}^{a}</math>
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<math>2.16\pm {0.40}^{ab}</math>
|}


Note: Differences between different letters in the same row are significant (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{p}<0.05</math>
|}
) and differences between the same mothers are not significant (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{p}>0.05</math>
|}
).

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 (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{p}>0.05</math>
|}
), but the differences in
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{CO}_{2}</math>
|}
and
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{NH}_{3}</math>
|}
concentrations were significant (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{p}<0.05</math>
|}
). 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{CO}_{2}</math>
|}
concentration between nodes B and C (
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>\mathit{p}>0.05</math>
|}
). This result indicates the test barn has better ventilation uniformity and is more suitable for piglet growth.

=4 Conclusions=

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.<br/>(1) The maximum relative error between the simulated and measured data is
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>20\%</math>
|}
, and the average relative error is
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>11.79\%</math>
|}
. 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
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>2.0\, m/s</math>
|}
, respectively, and the ventilation angle was
{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| <math>{45}^{\circ }</math>
|}
.<br/>(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.<br/>Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

==References==

:1. Wang B, Xiao H. Regional analysis on structure and competitiveness of China's animal husbandry-Based on shift-share analysis space model. Chin J Agric Resour Reg Plan. 2021;42(5):142-8 (In Chinese).

:2. Lesser T, Braun C, Wolfram F, Gottschall R. A special double lumen tube for use in pigs is suitable for different lung ventilation conditions. Res Vet Sci. 2020;133:111-6. doi:10.1016/j.rvsc.2020.09.007.

:3. Costa A, Borgonovo F, Leroy T, Berckmans D, Guarino M. Dust concentration variation in relation to animal activity in a pig barn. Biosyst Eng. 2009;104(1):118-24.

:4. Zong C, Li H, Zhang G. Ammonia and greenhouse gas emissions from fattening pig house with two types of partial pit ventilation systems. Agri Ecosyst Environ. 2015;208:94-105. doi:10.1016/j.agee.2015.04.031.

:5. Ma H, Xie Y, Li A, Zhang T, Liu Y, Luo X. A review on the effect of light-thermal-humidity environment in sow houses on sow reproduction and welfare. Reprod Domest Anim. 2023;8(8):1023-45. doi:10.1111/rda. 14400 .

:6. Dai X, Ni JQ, Pan Q, Wang K, Heber AJ. Monitoring of temperature, humidity and air quality inside pig weaner house in Eastern China. Trans Chin Soc Agric Eng. 2016;47(7):315-22.

:7. Wang D, Huang H, Zhang H, Liu J, Dong X. Analysis of research status and development on engineering technology of swine farming facilities. Trans Chin Soc Agric. 2018;49(11):1-14 (In Chinese). doi:10.6041/j.issn. 1000-1298.2018.11.001.

:8. Mostafa E, Hoelscher R, Diekmann B, Ghaly AE, Buescher W. Evaluation of two indoor air pollution abatement techniques in forced-ventilation fattening pig barns. Atmos Pollut Res. 2017;8(3):428-38. doi:10.1016/j.apr.2016.11.003.

:9. Thalfeldt M, Kurnitski J, Latôšov E. Exhaust air heat pump connection schemes and balanced heat recovery ventilation effect on district heat energy use and return temperature. Appl Ther Eng. 2018;128:402-14. doi:10.1016/j.applthermaleng.2017.09.033.

:10. Li X, Ye Z, Li B, Feng X, Zhu S, Shen P. Influence of different ventilation system on environment of nursery piggery in winter. Trans Chin Soc Agricult Mach. 2020;51:317-25 (In Chinese). doi:10.6041/j.issn.1000-1298.2020.03.036.

:11. Kwon KS, Lee IB, Ha T. Identification of key factors for dust generation in a nursery pig house and evaluation of dust repipeion efficiency using a CFD technique. Biosyst Eng. 2016;151(3):28-52. doi:10.1016/j.biosystemseng.2016.08.020.

:12. Xu X, Ge B, Tao D, Han J, Wang L. 3D temperature field of high-temperature gas cooling reactor cooling medium drive motor and ventilation structure improvement. IET Electr Pow Appl. 2018;12(7):1020-6. doi:10.1049/iet-epa.2017.0661.

:13. Shi L, et al. Experimental study on the application effect of targeted ventilation inmulti-floor farrowing sow house. J China Agricult Univ. 2023;28(2):217-26 (In Chinese). doi:10.11841/j.issn.1007-4333.2023.02.19.

:14. Zong C, Feng Y, Zhang G, Hansen MJ. Effects of different air inlets on indoor air quality and ammonia emission from two experimental fattening pig rooms with partial pit ventilation system-Summer condition. Biosyst Eng. 2014;122:163-73.

:15. Hu Z, Yang Q, Tao Y, Shi L, Tu J, Wang Y. A review of ventilation and cooling systems for large-scale pig farms. Sustain Cities Soc. 2023;89:104372. doi:10.1016/j.scs.2022.104372.

:16. Kim RW, Kim JG, Lee IB, Yeo UH, Lee SY. Development of a VR simulator for educating CFD-computed internal environment of piglet house. Biosyst Eng. 2019;188(1):243-64. doi:10.1016/j.biosystemseng.2019.10.024.

:17. Wang K, Pan Q, Li K. Computational fluid dynamics simulation of the hygrothermal conditions in a weaner house in Eastern China. Trans ASABE. 2017;60(1):195-205. doi:10.13031/trans.11655.

:18. Qi F, Zhao X, Shi Z, Rong L, Zhang G, Li H. Applicability evaluation of innovative simplified methods of slatted floor in pig houses-A CFD study. Comput Electron Agric. 2024;216(4):108532. doi:10.1016/j.compag.2023.108532.

:19. Li H, Rong L, Zhang G. "Reliability of turbulence models and mesh types for CFD simulations of a mechanically ventilated pig house containing animals. Biosyst Eng. 2017;161:37-52. doi:10.1016/j.biosystemseng.2017.06.012.

:20. Wei X, Li B, Lu H, Lü E, Guo J, Jiang Y, et al. Numerical simulation of airflow distribution in a pregnant sow piggery with centralized ventilation. Appl Sci. 2022;12(22):11556.

:21. Cheng Q, Wu W, Li H, Zhang G, Li B. CFD study of the influence of laying hen geometry, distribution and weight on airflow resistance. Comput Electron Agric. 2018;144:181-9.

:22. Yeo UH, Lee IB, Kim RW, Lee SY, Kim JG. Computational fluid dynamics evaluation of pig house ventilation systems for improving the internal rearing environment. Biosyst Eng. 2019;186:259-78.

:23. Mossad RR. Numerical predictions of air temperature and velocity distribution to assist in the design of natural ventilation piggery buildings. Aust J Multi-Discip Eng. 2011;8(2):181-7. doi:10.1080/14488388.2011.11464837.

:24. Guo J, Wei X, Du X, Ren J, Lü E. Numerical simulation of liquid nitrogen injection in a container with controlled atmosphere. Biosyst Eng. 2019;187(1):53-68. doi:10.1016/j.biosystemseng.2019.08.015.

:25. Bjerg B, Cascone G, Lee IB, Bartzanas T, Norton T, Hong SW, et al. Modelling of ammonia emissions from naturally ventilated livestock buildings. Part 3: CFD modelling. Biosyst Eng. 2013;116(3):259-75. doi:10.1016/j.biosystemseng.2013.06.012.

:26. Qin J, Pan H, Rahman MM, Tian X, Zhu Z. Introducing compressibility with SIMPLE algorithm. Math Comput Simulat. 2021;180(1):328-53. doi:10.1016/j.matcom.2020.09.010.

:27. Wu Z, Fang J, Ba W, Wu S, Li H, Cai Y, et al. Design and evaluation of combined ventilation system with air supply and exhaust ducts for a pig nursery heating piggery in winter. Trans Chin Soc Agric Eng. 2021;37(10):152-8 (In Chinese).

:28. Liu H, Fang H, Cheng R, Yang Q. Simulation and optimization of the air flow and temperature in plant factory with artificial light based on CFD. J China Agric Univ. 2018;23(5):108-16 (In Chinese).

:29. Wang X, Chen Z, Wang M, Liu J. Numerical simulation of heat supply for heat recovery ventilation system of piggery in winter. Trans Chin Soc Agric Eng. 2011;27(12):227-33+438 (In Chinese).

:30. Deng S, Shi Z, Li B, Zhao S, Ding T, Zheng W. CFD simulation of airflow distribution in low profile cross ventilated dairy cattle barn. Trans Chin Soc Agric Eng. 2014;30(6):139-46 (In Chinese).

:31. Wang H, Liu J, Wu Z, Feng G, Shen Z, Wang M. Airflow characteristics of attachment ventilation in a nursery pig house under heating mode. Biosyst Eng. 2022;224(9):346-60. doi:10.1016/j.biosystemseng.2022.10.018.

:32. Gu Y, Mu J, Zheng S, Gang Z, Jing R, Wang C. Effect of jet hole arrangement on drag reduction characteristics of jet surface. Trans Chin Soc Agric. 2014;45(10):340-6 (In Chinese).

:33. Xing D, Yan C, Wang C, Sun L. Effects of aspect ratio of rectangular channel on characteristics of singlephase laminar flow. Chin J Theoretical Appl Mech. 2013;45(3):331-6 (In Chinese).

:1. *Correspondence: Zhidong Wu ([mailto:wzd139446@163.com wzd139446@163.com]). This is an article distributed under the terms of the Creative Commons BY-NC-SA license