Impact Analysis of Environmental Conditions on Odour Dispersion Emitted from Pig House with Complex Terrain Using CFD

: With the rise in livestock production, the odour concentration emitted from livestock facilities has signiﬁcantly increased. For this reason, odour complaints have become a major issue. The dispersion of the odour into the atmosphere is a ﬀ ected by various factors (wind speed, wind direction, atmospheric stability, terrain condition, ventilation type, and so on). Thus, a thorough analysis on the factor inﬂuencing odour dispersion is necessary to establish regulations and policies for odour management. Therefore, this research aimed to evaluate odour dispersion generated from a pig house with complex terrain using computational ﬂuid dynamics (CFD) and to statistically determine the key factor for odour dispersion. By comparing CFD-computed results with ﬁeld-measured data, an appropriate grid size, time step, and turbulence model of the CFD model were determined. Considering various factors, case studies were performed using the validated CFD model. The CFD-computed results showed that odour dispersion distance at the level of 1 OU m − 3 ranged from 129.7 to 1488.1 m. The prediction of odour dispersion distance varies with the change of factor being analysed. Finally, the statistical analysis showed that the most inﬂuential factor that a ﬀ ected odour dispersion distance was the wind speed with a t-value: − 9.057.


Introduction
Pig production in the Republic of Korea accounts for 40% of the total livestock industry (USD 1.8 billion) in the country [1]. The number of pig farms was continuously reduced to 50.3% from 12,227 in 2005 to having only 6160 pig farms for the year 2019, but, the scale of pig rearing was significantly increased from 8,961,505 heads in 2005 to 11,317,000 heads in 2019 [2]. This indicates that pig farms have progressively been enlarged to meet the rise in meat demand. The increase in production intensified manure generation, causing an inevitable rise in the number of complaints related to odour emission. The number of complaints related to livestock odour upsurged to 27% from 4302 in 2005 to 6112 in 2017 [3].
To cope with this problem, the government has established policies to regulate and minimize the effect of odour in the vicinity of pig houses. This includes but not limited to the odour prevention law related to the allowable value of odour emission from livestock facilities and the establishment of setback distance. Research on odour management using physical and chemical methods has also been conducted. Past research has focused on studies to determine the optimum practice to minimize odour nuisance through field experiments [4][5][6]. Most studies mainly evaluated the efficiency of additives (chemical and biological methods) or abatement facilities (biofilters, liquid manure recharging systems, and so on) to reduce livestock odours [7][8][9][10][11][12][13][14][15]. These field experiment studies usually measured and evaluated single odorous gas such as ammonia (NH 3 ), hydrogen sulphide (H 2 S), methane (CH 4 ), and carbon dioxide (CO 2 ). However, an evaluation of complex odour is also necessary because the odour concentration of individual single gas and the complex odour emitted from animal facilities are different. Furthermore, the odour causing a nuisance to the nearby community tends to be caused by a combination of two or more substances rather than a single substance.
Sensory evaluation (dynamic olfactometry method) is used to estimate the quantitative amount of complex odour at a specific location and time. Although this method can obtain the actual value in the field, multiple panels are required for the measurement. Despite the established guidelines for the selection and training of the sensory panel members, the different odour perception ability of each panel can contribute to the errors in measuring odour. Furthermore, the sensory evaluation method has several disadvantages such as difficulties in managing the samples and being time-consuming considering multipoint measurements [16,17]. For these reasons, the evaluation of odour dispersion outside livestock buildings using the sensory evaluation method was not intensively considered and was used only for evaluating the odour concentration present at a specific location. The influence of the distance between the source and the receptor must be considered to control the odour from the outside. The impact of the weather conditions on odour dispersion should also be taken into account, especially weather conditions such as wind speed, wind direction, and atmospheric stability, which have a great influence on dispersion. The field experiment study of Yeo et al. [18] measured complex odour using sensory evaluation and measured weather conditions. However, the conduct of odour field measurements required intensive labour and excessive cost, as it considered the effect of various meteorological conditions in the odour emission. Furthermore, due to the difficulties of controlling the environmental conditions, the results of the field measurement were only valid in a specific weather condition during the measurement period. To overcome these, numerical dispersion models that mathematically calculate the dispersion of pollutants into the atmosphere and change of pollutant concentration were developed and used in various studies [19,20]. However, like field experiments, numerical dispersion models also have limitations such that they have difficulty in imitating the actual dispersion trend when not properly modelled [21].
One of the widely used numerical analysis methods for predicting the dispersion of the odour is the Gaussian dispersion model [22][23][24][25][26][27][28][29][30]. This dispersion model assumes that the horizontal and vertical odour concentration in the plume has a normal distribution (Gaussian distribution). Owing to the required short calculation time of Gaussian dispersion models and simple calculation procedure, these dispersion models have been frequently used. Specifically, the Gaussian dispersion models were found to be effective in research studies where the main objective was to know the distance and concentration ranges of pollutants blown to the atmosphere by chimneys [31] and typically used in point source emissions [32]. However, when the mesh resolution for the Gaussian dispersion model was set to a very large grid, the odour concentration prediction of the model resulted in an overestimation of odour flow. In addition, in the Gaussian dispersion model, it was difficult to analyse the odour dispersion in detail by considering the local wind environment generated from the complex terrain Agronomy 2020, 10, 1828 3 of 27 and layout of livestock facilities. Other commonly used odour dispersion models were thoroughly discussed in the paper of Conti et al., 2020 [33].
In comparison, the numerical analysis of atmospheric odour dispersion using computational fluid dynamics (CFD) allowed more realistic simulation and can compensate for the limitations of the Gaussian dispersion model even though it required detailed geometry design and more computational time. The importance of CFD has been validated in many fields including industrial applications [34][35][36], environmental applications [37][38][39], and physiological applications [40,41]. Specifically, the application of CFD in the field of agriculture such as in odour transport [29,[42][43][44][45][46][47] and dust emission [48] was recently conducted. However, these studies used flat or very small areas to simulate dispersion phenomena with CFD.
To accurately predict odour dispersion, it is necessary to consider the topographical effects in the target areas where the local wind environment is affected by complex terrain such as the mountainous area. Hence, recent research tried to simulated airflow analysis and the change of the odour concentration emitted from livestock buildings, taking into account the complex terrain in the computational domain [20,49]. Specifically, Hong et al. [20] developed and validated an initial CFD model with complex terrain to simulate the odour dispersion of the pig house, considering weather conditions such as atmospheric stability and wind direction. However, the paper assumed the odour source was area source and the odour that was emitted at the surface of the pig farm was considered as the odour source. Furthermore, these research studies mainly focused on analysing the trend of odour dispersion outside the pig house, and there was a lack of research on the analysis of key factors affecting the dispersion of odour in the atmosphere. For instance, Lin et al. [29] concluded that wind speed and wind direction have an impact on the dispersion of odour in the atmosphere. The same statement was made by Chastain [30], implying that the direction of prevailing wind and the presence of wind breaks must be considered in the model. The ventilation inside the pig farm also contributed to the amount of odour emitted from the livestock facilities. These factors were agreed by Capelli et al. [28], which emphasized in the manuscript the importance of meteorological information needed for dispersion modelling such as wind speed, wind direction, and atmospheric stability. In addition, Chastain [30] includes the influence of the topographical condition of the target farm. Papers such as those published by Ubeda et al. [50], Hong et al. [51], and Conti et al. [33] support the other claims of the previously mentioned papers relating to the factors that affected the dispersion of odour. From the literature, wind speed, wind direction, atmospheric stability, and ventilation type were used as dependent factors in our research.
Therefore, the objective of this study was to evaluate the key factors (ventilation type, atmospheric stability, wind speed, and wind direction) that affect the dispersion of livestock odour into the atmospheric using CFD. To do this, a three-dimensional CFD model of an experimental pig house with complex terrain characteristics was designed, validated, and used to predict odour dispersion. To complement the limitation of the study conducted by Hong et al. [20], the structural characteristic of the pig house model was considered in the calculation domain. Specifically, the odorous pollutants were set to be emitted from the exhaust fans of the pig houses. The CFD model was validated through the odour concentration measured at various sampling points from field experiments [18]. Then, various cases were also conducted using the conducted CFD model according to the target factor considered. Finally, a statistical analysis was performed to determine which of the factors is the main factor that influenced odour dispersion. Figure 1 shows the flow chart for determining the key factor that influenced odour dispersion distance. The study was divided into two main contents: a field experiment [18] and a CFD simulation. For the field experiment, an area with a single odour source was selected to ensure that odours were generated from experimental pig farms. A field experiment for measuring the external environmental conditions such as wind speed, wind direction, air temperature, atmospheric stability, and odour Agronomy 2020, 10, 1828 4 of 27 concentrations (inside and outside the pig house) was conducted first. The result of the field experiment was used in this research to establish the boundary conditions for the design of the three-dimensional CFD model and was also used for the validation of the designed model. The validated CFD model was used to evaluate the dispersion range according to various factors (wind direction, wind speed, atmospheric stability, ventilation types) affecting odour dispersion. The influence of each factor was evaluated using statistical analysis, wherein the factors affecting odour dispersion were set as independent variables and odour dispersion distance was set as the dependent variable.

Materials and Methods
Agronomy 2020, 10, x 4 of 27 direction, wind speed, atmospheric stability, ventilation types) affecting odour dispersion. The influence of each factor was evaluated using statistical analysis, wherein the factors affecting odour dispersion were set as independent variables and odour dispersion distance was set as the dependent variable.

Farm Description
In the selection of the target experimental farm, it was very important to exclude the interference of livestock odour that was generated from surrounding farms. A commercial pig farm to meet this requirement that was located at Bouen-gun, Chungcheongbukdo, South Korea (127°37′15″ E, 36°28′56″ N) was selected as an experimental farm. The pig farm was located in a valley formed in a north-south orientation. The experimental pig farm consisted of three pig houses, in which, in this study, were referred to as pig house 1 (P1), pig house 2 (P2), and pig house 3 (P3), as shown in Figure 2a. The experimental pig farm raised a total of 850 pigs, with lower than standard rearing density levels of 0.8 m 2 per finishing pig [52,53]. The farm adopted an all-in and all-out rearing (AIAO) method. P1 has dimensions of 28.0 m in length, 10.0 m in width, and 3.8 m in height. In contrast, P2 and P3, with a perpendicular orientation to P1, have 28.0 × 13.6 × 4.5 m and 33.0 × 10.0 × 4.3 m in length, width, and height, respectively. All pig houses were equipped with a combined mechanical ventilation and natural ventilation system, consisting of winch curtains and roof-chimney exhaust fans (diameter (D): 0.6 m, 168 cubic meters per minute (CMM)). Roof-chimney exhaust fans were installed at equal intervals in the longitudinal direction of P1 (two fans) and P3 (three fans). On the other hand, considering the relatively wide width of P2, a total of four fans were installed in two rows in the longitudinal direction on the roof of P2. The dimension of the winch curtain for P1 and P2 was 1.0 × 27.0 m, while the winch curtain for P3 was 1.0 × 32.0 m. The natural air entered the pig houses through winch curtains that were installed on both sides of each house. The location of the slurry pit was below the pigs and the manure was collected, stored, and removed at the end of every rearing period. The inside and outside conditions of the pig house are shown in Figure 2b,c.

Farm Description
In the selection of the target experimental farm, it was very important to exclude the interference of livestock odour that was generated from surrounding farms. A commercial pig farm to meet this requirement that was located at Bouen-gun, Chungcheongbukdo, South Korea (127 • 37 15" E, 36 • 28 56" N) was selected as an experimental farm. The pig farm was located in a valley formed in a north-south orientation. The experimental pig farm consisted of three pig houses, in which, in this study, were referred to as pig house 1 (P1), pig house 2 (P2), and pig house 3 (P3), as shown in Figure 2a. The experimental pig farm raised a total of 850 pigs, with lower than standard rearing density levels of 0.8 m 2 per finishing pig [52,53]. The farm adopted an all-in and all-out rearing (AIAO) method. P1 has dimensions of 28.0 m in length, 10.0 m in width, and 3.8 m in height. In contrast, P2 and P3, with a perpendicular orientation to P1, have 28.0 m × 13.6 m × 4.5 m and 33.0 m × 10.0 m × 4.3 m in length, width, and height, respectively. All pig houses were equipped with a combined mechanical ventilation and natural ventilation system, consisting of winch curtains and roof-chimney exhaust fans (diameter (D): 0.6 m, 168 cubic meters per minute (CMM)). Roof-chimney exhaust fans were installed at equal intervals in the longitudinal direction of P1 (two fans) and P3 (three fans). On the other hand, considering the relatively wide width of P2, a total of four fans were installed in two rows Agronomy 2020, 10, 1828 5 of 27 in the longitudinal direction on the roof of P2. The dimension of the winch curtain for P1 and P2 was 1.0 m × 27.0 m, while the winch curtain for P3 was 1.0 m × 32.0 m. The natural air entered the pig houses through winch curtains that were installed on both sides of each house. The location of the slurry pit was below the pigs and the manure was collected, stored, and removed at the end of every rearing period. The inside and outside conditions of the pig house are shown in Figure 2b,c.

Korean Standard Pig House
Previously constructed pig houses in the Republic of Korea had been arbitrarily designed depending on the experience of the consultant and allotted budget of the farm owners. This resulted in various structural designs and ventilation systems' designs of pig houses, causing difficulties in the implementation of technologies to easily maintain a rearing environment, working hygiene, and odour management. Therefore, the Korean government has provided a standard pig house design to be used for constructing new and for renovating old pig houses. This standard design can reduce the economic burden of farm owners at the design stage of the pig house. Moreover, the availability of these standard designs was enormous because government legislation on unauthorized livestock buildings was enacted, and plans were being made for crackdowns and sanctions on the construction of unauthorized livestock buildings. Thus, to represent the approved design of pig houses and common internal and external conditions of a pig house, the CFD model for the case study used the 2016 standard pig house design for the finishing pigs [54].
The schematic diagram of the standard pig house is shown in Figure 3a,b. The length of the Korean standard pig house, including the number and type of exhaust fans, can be modified depending on the farm area, the number of pigs grown, and the farm owner's preference. Therefore, a fully mechanically ventilated pig house with a similar area as the pig farm in the experimental farm was considered in the computational domain for numerical analysis. The orientation of three standard pig houses followed the same position as the actual experimental pig house. An inlet (1.5 × 0.6 m), which acts as the entrance for fresh air, was also provided on both sides of each structure. Each pig house has 16 sidewall exhaust fans (D: 0.7 m) and 12 roof-chimney exhaust fans (D: 0.5 m). The designed structure of the pig house is shown in Figure 3c.

Korean Standard Pig House
Previously constructed pig houses in the Republic of Korea had been arbitrarily designed depending on the experience of the consultant and allotted budget of the farm owners. This resulted in various structural designs and ventilation systems' designs of pig houses, causing difficulties in the implementation of technologies to easily maintain a rearing environment, working hygiene, and odour management. Therefore, the Korean government has provided a standard pig house design to be used for constructing new and for renovating old pig houses. This standard design can reduce the economic burden of farm owners at the design stage of the pig house. Moreover, the availability of these standard designs was enormous because government legislation on unauthorized livestock buildings was enacted, and plans were being made for crackdowns and sanctions on the construction of unauthorized livestock buildings. Thus, to represent the approved design of pig houses and common internal and external conditions of a pig house, the CFD model for the case study used the 2016 standard pig house design for the finishing pigs [54].
The schematic diagram of the standard pig house is shown in Figure 3a,b. The length of the Korean standard pig house, including the number and type of exhaust fans, can be modified depending on the farm area, the number of pigs grown, and the farm owner's preference. Therefore, a fully mechanically ventilated pig house with a similar area as the pig farm in the experimental farm was considered in the computational domain for numerical analysis. The orientation of three standard pig houses followed the same position as the actual experimental pig house. An inlet (1.5 × 0.6 m), which acts as the entrance for fresh air, was also provided on both sides of each structure. Each pig house has

Design Tool for Numerical Analysis
To accurately evaluate the odour dispersion on a complex terrain like a mountain area, a detailed topographical design should be properly designed before incorporating it into the computational domain [20]. In this study, the complex terrain was developed using commercial software such as ArcGIS, AutoCAD, SketchUp, and Rhino3D. A topographical map of the experimental area, which consisted of contours, elevations, and land use information on the scale of 1:5000 and was obtained from the National Geographic Information Institute (NGII), was used. The ArcMap tool, a software under ArcGIS, was used to edit the digital map. However, when the maps were laid together, borders between the maps were slightly misaligned. Thus, AutoCAD was used to manually modify the contours. It was also used to design the buffer area outside the contour area. The SketchUp tool was used to generate the triangular irregular network (TIN) from the contour lines. Rhino3D, on the other hand, is a modelling tool that can provide great flexibility and precision in handling curves and surfaces. This tool was used to convert the ground from the formed TIN triangulation to solid modelling. The completed ground surface was incorporated in the geometry of the CFD model to analyse the odour dispersion with complex terrain. CFD is a numerical analysis method for predicting phenomena such as fluid flow, heat and mass transfer, and chemical reactions. The governing fluid flow Navier-Stokes equations, which are a set of nonlinear partial differential equations, were formulated in a discretized form through the finite volume method (FVM). Numerical analyses for fluid and energy flow are based on mass, momentum, and energy conservation laws. The conservation equations (Equations (1)-(3)) for each of these physical variables are as follows [55]: Mass conservation equation: Energy conservation equation:

Design Tool for Numerical Analysis
To accurately evaluate the odour dispersion on a complex terrain like a mountain area, a detailed topographical design should be properly designed before incorporating it into the computational domain [20]. In this study, the complex terrain was developed using commercial software such as ArcGIS, AutoCAD, SketchUp, and Rhino3D. A topographical map of the experimental area, which consisted of contours, elevations, and land use information on the scale of 1:5000 and was obtained from the National Geographic Information Institute (NGII), was used. The ArcMap tool, a software under ArcGIS, was used to edit the digital map. However, when the maps were laid together, borders between the maps were slightly misaligned. Thus, AutoCAD was used to manually modify the contours. It was also used to design the buffer area outside the contour area. The SketchUp tool was used to generate the triangular irregular network (TIN) from the contour lines. Rhino3D, on the other hand, is a modelling tool that can provide great flexibility and precision in handling curves and surfaces. This tool was used to convert the ground from the formed TIN triangulation to solid modelling. The completed ground surface was incorporated in the geometry of the CFD model to analyse the odour dispersion with complex terrain. CFD is a numerical analysis method for predicting phenomena such as fluid flow, heat and mass transfer, and chemical reactions. The governing fluid flow Navier-Stokes equations, which are a set of nonlinear partial differential equations, were formulated in a discretized form through the finite volume method (FVM). Numerical analyses for fluid and energy flow are based on mass, momentum, and energy conservation laws. The conservation equations (Equations (1)-(3)) for each of these physical variables are as follows [55]: Mass conservation equation: Momentum conservation equation: where ρ is the density (kg m −3 ), → υ is the velocity (m s −1 ), S m is a mass source term based on the chemical reaction (kg m −2 ), E is the total energy (kg m 2 s −2 kg −1 ), P s is the atmospheric pressure (kg m −1 s −2 ) ,

Field Experiment for CFD Model Design and Validation
A data field measurement conducted by Yeo et al. [18] from August 2018 to September 2019 was used to define the initial boundary conditions and validate the designed three-dimensional CFD model. An odour assessment result carried out during four sampling periods was used. The external weather environment that affects odour dispersion and quantitative odour concentration was measured from the field experiment. Shown in Figure 4a,b are the actual setups for collecting odorous air samples inside and outside the pig house facilities. To observe the external weather environment such as wind direction and wind speed, a portable weather station (WatchDog 2700, Spectrum Technologies, Inc., Aurora, IL, USA) was placed at 1.5 m from the roof of P1 with a measuring interval of 1 min (Figure 4c). For on-site atmospheric stability, the vertical distribution of temperature was measured through the hobo sensor (UX100-003, Onset Computer Corp., Bourne, MA, USA) twined on the balloons using a light wire (Figure 4d). The measurement location of the atmospheric stability was set up to a place so that it was free from any obstacles around the experimental pig house. The atmospheric stability was measured at 1-min intervals. Odour generated from pig houses was also measured to determine the odour concentration dispersed at a certain distance at specific atmospheric stability. The measured odour concentration was used to validate the CFD model. Odour samples at different sampling locations were collected simultaneously using a portable air sampler (Odortech, Seoul, Korea). The actual setup of field instruments is illustrated in Figure 4. The sampling locations for the field measurement of odour concentration were inside the pig houses and outside the pig houses. All measuring devices were located 1.5 m high from the ground surface. For inside the pig house, sampling instruments were installed near the location of the outlet. The sampling location outside the pig houses was established at 40.0 to 50.0 m intervals from the farm's boundary line. Sampling locations inside the pig houses were designated as PS1, PS2, and PS3 depending on the pig house, while in the case of sampling points outside, S1, S2, S3, and S4 were used ( Figure 5).
Following the field experiment, a sensory evaluation was conducted to determine the concentration of odour obtained from various sampling points. The test method followed the air dilution sensory guidelines of the odour process test (ES 09301.a) [53,54]. The collected air samples were transported to the laboratory for sensory evaluation within 24 h and were kept at room temperature (15 to 25 • C) to maintain the quality of samples obtained from the field. The odourless air was made from the clean air generator (Anytech Inc., Suwon, Korea). This device contained silica gel and activated carbon to eliminate moisture and odour. The diluted air samples were made by mixing odorous and odourless air in consideration of the dilution rate, by injecting the odorous air sample into Tedlar bags with clean air samples. The diluted air samples were distributed to the panels through a 5-p distributor (Anytech Inc., Suwon, Korea). When conducting the evaluation, the same panels were employed during the experiment period to eliminate the uncertainty in the odour concentration caused by employing different panels. Each set comprised three air samples, where one sample contained odour and two with only odourless air. The sensory evaluation followed the "triangular-forced-choice" process, where panellists sniff all three samples and must select which among the three samples contain the diluted odour. The odour concentration based on the evaluation results performed by five panellists, excluding the maximum and minimum values to prevent from being overestimated or underestimated, was calculated using the geometric mean [53].
Agronomy 2020, 10, x 8 of 27 one sample contained odour and two with only odourless air. The sensory evaluation followed the "triangular-forced-choice" process, where panellists sniff all three samples and must select which among the three samples contain the diluted odour. The odour concentration based on the evaluation results performed by five panellists, excluding the maximum and minimum values to prevent from being overestimated or underestimated, was calculated using the geometric mean [53].

Design, Validation, and Case Study Using the CFD Model
A three-dimensional CFD model for the simulation of livestock odour dispersion is shown in Figure 6; it has an area of 5.0 km in diameter, including 2 km of the buffer area, and 2.0 km in height. Because the odour dispersion problem is mainly generated in the lower part of the atmospheric layer, and the height of the grid located near the ground surfaces large, it is difficult to simulate the rapid wind speed change near the ground. Therefore, the domain was discretized by a non-structural grid strategy. This means that the size and density of grids used in the domain were not uniform. Purposely, the grid near the pig houses and the ground was designed denser and became coarser one sample contained odour and two with only odourless air. The sensory evaluation followed the "triangular-forced-choice" process, where panellists sniff all three samples and must select which among the three samples contain the diluted odour. The odour concentration based on the evaluation results performed by five panellists, excluding the maximum and minimum values to prevent from being overestimated or underestimated, was calculated using the geometric mean [53].

Design, Validation, and Case Study Using the CFD Model
A three-dimensional CFD model for the simulation of livestock odour dispersion is shown in Figure 6; it has an area of 5.0 km in diameter, including 2 km of the buffer area, and 2.0 km in height. Because the odour dispersion problem is mainly generated in the lower part of the atmospheric layer, and the height of the grid located near the ground surfaces large, it is difficult to simulate the rapid wind speed change near the ground. Therefore, the domain was discretized by a non-structural grid strategy. This means that the size and density of grids used in the domain were not uniform. Purposely, the grid near the pig houses and the ground was designed denser and became coarser

Design, Validation, and Case Study Using the CFD Model
A three-dimensional CFD model for the simulation of livestock odour dispersion is shown in Figure 6; it has an area of 5.0 km in diameter, including 2 km of the buffer area, and 2.0 km in height. Because the odour dispersion problem is mainly generated in the lower part of the atmospheric layer, and the height of the grid located near the ground surfaces large, it is difficult to simulate the rapid wind speed change near the ground. Therefore, the domain was discretized by a non-structural grid strategy. This means that the size and density of grids used in the domain were not uniform. Purposely, the grid near the pig houses and the ground was designed denser and became coarser with the increase in altitude. This kind of grid generation approach ensures a high degree of grid precision in the area close to the odour sources [45]. To precisely compute odour concentration near the surface, the CFD model was designed by using the inflation function. The first layer height used 0.1 m at the ground surface, and 25 layers with a growth rate of 1.1. To validate the CFD model, a standard design of a pig house model was placed in the centre of the computational domain, as shown in Figure 6. A grid size of 0.2 m around the roof-chimney exhaust fan of the pig house was used to design a grid where the odour was mainly emitted. Other than the roof-chimney exhaust fan, 0.5 m of grid size was used around the wall of the pig house. Additionally, to calculate more accurate airflow and odour distribution around the pig house, the grid was designed with a size of 2.0 m in the 200.0 m radius around the pig house. After 200.0 m, the CFD model was designed with a growth rate of 1.1 to reduce the total number of grids, while gradually increasing the grid size ( Figure 7). In the same way, the CFD model (the 2016 standard pig house design) for the case study applied all the grid design methods used in the design of the validation model. precision in the area close to the odour sources [45]. To precisely compute odour concentration near the surface, the CFD model was designed by using the inflation function. The first layer height used 0.1 m at the ground surface, and 25 layers with a growth rate of 1.1. To validate the CFD model, a standard design of a pig house model was placed in the centre of the computational domain, as shown in Figure 6. A grid size of 0.2 m around the roof-chimney exhaust fan of the pig house was used to design a grid where the odour was mainly emitted. Other than the roof-chimney exhaust fan, 0.5 m of grid size was used around the wall of the pig house. Additionally, to calculate more accurate airflow and odour distribution around the pig house, the grid was designed with a size of 2.0 m in the 200.0 m radius around the pig house. After 200.0 m, the CFD model was designed with a growth rate of 1.1 to reduce the total number of grids, while gradually increasing the grid size ( Figure 7). In the same way, the CFD model (the 2016 standard pig house design) for the case study applied all the grid design methods used in the design of the validation model.  The proper design for the grid size is an important factor that can influence not only the computational time but also the accuracy and reliability of the results. The smaller the grid size, the more detailed and accurate the numerical value can be calculated. However, the computation time increases as the number of grids increases, which may not be efficient. Therefore, it is necessary to conduct the grid independence test, which selects a reasonable grid size by designing grids of various precision in the area close to the odour sources [45]. To precisely compute odour concentration near the surface, the CFD model was designed by using the inflation function. The first layer height used 0.1 m at the ground surface, and 25 layers with a growth rate of 1.1. To validate the CFD model, a standard design of a pig house model was placed in the centre of the computational domain, as shown in Figure 6. A grid size of 0.2 m around the roof-chimney exhaust fan of the pig house was used to design a grid where the odour was mainly emitted. Other than the roof-chimney exhaust fan, 0.5 m of grid size was used around the wall of the pig house. Additionally, to calculate more accurate airflow and odour distribution around the pig house, the grid was designed with a size of 2.0 m in the 200.0 m radius around the pig house. After 200.0 m, the CFD model was designed with a growth rate of 1.1 to reduce the total number of grids, while gradually increasing the grid size ( Figure 7). In the same way, the CFD model (the 2016 standard pig house design) for the case study applied all the grid design methods used in the design of the validation model.  The proper design for the grid size is an important factor that can influence not only the computational time but also the accuracy and reliability of the results. The smaller the grid size, the more detailed and accurate the numerical value can be calculated. However, the computation time increases as the number of grids increases, which may not be efficient. Therefore, it is necessary to conduct the grid independence test, which selects a reasonable grid size by designing grids of various The proper design for the grid size is an important factor that can influence not only the computational time but also the accuracy and reliability of the results. The smaller the grid size, the more detailed and accurate the numerical value can be calculated. However, the computation time increases as the number of grids increases, which may not be efficient. Therefore, it is necessary to conduct the grid independence test, which selects a reasonable grid size by designing grids of various sizes in the calculation domain and comparing them with field-measured data. This takes into account the computational characteristics that the grid results no longer change significantly as the grid size becomes smaller and denser. In this study, the grid sizes from the outside of the pig house to the buffer space in the calculation domain were fixed to 4.0, 10.0, 20.0, and 40.0 m, except for the inlet, outlet, and ground surface, which requires the design of a small grid in the flow analysis, and was compared with field experiments.
Similarly, the selection of an appropriate time step to be used in the validation is an important factor that influenced the accuracy and economics of the simulation procedure. In general, the smaller the size of the calculation time step, the better the ability to simulate certain phenomena. In this study, the time steps of 1, 5, 10, and 20 s were considered for a calculation time of 10 min.
In analysing the flow of air, the choice of turbulence model is important because the turbulence characteristics may vary according to the turbulence model selected in the computational domain. Turbulent flow is a type of flow that undergoes irregular fluctuation or mixing. Turbulent flow has a smaller dispersion compared to laminar flow with a low Reynolds number, convection acts largely in a certain space, and pressure and velocity change rapidly over time and space. Therefore, tests on turbulence models, which are widely applied in the field of flow analysis in agricultural facilities, were conducted. In the same way as the grid independence test, an appropriate turbulence model to be able to improve the accuracy of the CFD-computed results was determined by comparing the field-measured data with the CFD-computed results with R 2 (R-square) and RMSE (root mean square error). The tested turbulence models include RANS (Reynolds-averaged Navier-Stokes) models (standard k-epsilon model, renormalized group (RNG) k-epsilon model, realizable k-epsilon model) and the Large Eddy Simulation (LES) model. The RANS model is time-averaged equations of motion for fluid flow. RANS models would weaken and be replaced by large eddy simulation (LES) as a consequence of rapidly increasing computational power. Certainly, greater reliance on the unaveraged Navier-Stokes equations means more accurate results. The k-epsilon (k-ε) turbulence model for the RANS model is the most common model used in CFD to simulate mean flow characteristics for turbulent flow conditions. Among the k-epsilon turbulence models, there are differences in the turbulent viscosity and transport equation for turbulent dissipation. Therefore, the effort of selecting a suitable turbulence model from the k-epsilon and LES models for turbulence analysis was made. Statistical indices, such as R 2 and RSME, were used to determine the grid sizes and turbulent models were appropriate to be used in this study.
The stability or instability of the atmosphere is a condition that considers the reaction of a volume of an air parcel to vertical displacement. The stability of the atmosphere determines the probability of convective activity, cloud type, atmospheric turbulence, and degree of vertical mixing (for pollutants, etc.). The classification of atmospheric stability includes stable, neutral, and unstable conditions. Atmospheric stability is related to the wind profile expressed as Equation (4). Specifically, the wind profile was based on a log wind profile which is used to describe the vertical wind characteristics.
where k is the von Karman constant (0.42); z is the altitude above ground; u* is the friction velocity; u is the velocity magnitude; Φ m is the Monin-Obukhov universal function for non-dimensional momentum.
To express the Monin-Obukhov universal function for the stable atmospheric stability in this study, the equations (Equations (5)- (7)) were adopted from the empirical equation created by Dryer and Hicks [56] and Hicks [57] as cited by Hong et al. [20]. The Φ m was expressed using the following equations: : unstable (7) where L is the Obukhov length and z is the height from the surface.
The assumption in the Monin-Obukhov similarity theory was adopted that when the calculated value of the Monin-Obukhov universal function has a value of infinity and almost equal to 1, the atmospheric stability is neutral. Similarly, stable air happens when the surface absorbs heat and the Obukhov length (L) is positive. Conversely, the unstable air condition occurs when the Obukhov length (L) is negative. Therefore, by integration, the following equations used to describe the vertical wind profile were expressed as Equations (8)-(10): where U(z) is the magnitude of the mean horizontal velocity at the height z (z > z 0 ); u* is the friction velocity; k is the von Karman constant; z is the height from the surface; z 0 is the roughness length of the surface; L is the Obukhov length; x is the dimensionless variable defined by the equation used to describe the Monin-Obukhov universal function. Several qualitative and quantitative measures have been conducted by many researchers to create an air temperature profile. In this study, the dry adiabatic lapse rate for microclimate conditions cited by Bache and Johnstone [58] was adopted to present the temperature turbulence or the change in temperature according to height. For a neutral atmosphere, the temperature cools down by 1 • C for every 100 m increase in height (−0.01 • C m −1 ). Meanwhile, a lapse rate of 0.015 and −0.02 • C −1 was used for stable and unstable conditions, respectively. Table 1 shows the initial boundary conditions of the model for CFD simulation. The speed of the exhaust fan of the CFD model for validation was applied by measuring the wind speed of the fan at nine points when the exhaust fan ran at maximum (Figure 8). To measure the air flow rate of a circular exhaust fan, a traverse method was applied at the duct surface using hot wire anemometers to measure the discharged speed. The measurement points (9 points) were selected based on ASHRAE (American Society of Heating, Refrigerating and Air-Conditioning Engineers) standards [59].
While the wind speed of the exhaust fan for the case study was analysed by applying the ventilation flow rate provided by the manufacturer to the CFD model, in the case of the Korean standard pig house used for the case study, ventilation rates of the pig house were calculated from the ventilation requirement per pig in the summer season presented in the Midwest Plan Service (MWPS) for finishing pigs during the summer season. The summer season was preferred in this study, since odour dispersion was very dynamic during this season due to high air temperature and active microbial activities. The wind speeds of 6.8 and 11.6 m s −1 were applied to the sidewall and roof-chimney exhaust fan, respectively. Similarly, in the case of odour concentration generated in each pig house, the validation model applied the same odour value measured in the field. For the case study, the odour emissions from each pig house were calculated using emission factors established by Belgium [60]. The generation of odour for pigs in summer was defined as 32.7 OU animal −1 s −1 .
Considering the number of pigs per unit area, the total odour unit was 40,188.3 OU s −1 . The total odour concentration needs to be converted into a mass fraction value using the Equation (11) where Y s is the odour mass fraction; P a is the atmospheric pressure of 101,325 Pa at sea level; T the air temperature in K; M the molecular weight of dry air or 0.028966 kg s −1 ; R is the universal gas constant of 8.31432 J (mol·K) −1 ; O cg the odour source concentration in OU m −3 ; m c the mass of odour compound (n-butanol) required to produce one odour unit.
where U(z) is the magnitude of the mean horizontal velocity at the height z (z > z 0 ); u* is the friction velocity; k is the von Karman constant; z is the height from the surface; z is the roughness length of the surface; L is the Obukhov length; x is the dimensionless variable defined by the equation used to describe the Monin-Obukhov universal function. Several qualitative and quantitative measures have been conducted by many researchers to create an air temperature profile. In this study, the dry adiabatic lapse rate for microclimate conditions cited by Bache and Johnstone [58] was adopted to present the temperature turbulence or the change in temperature according to height. For a neutral atmosphere, the temperature cools down by 1 °C for every 100 m increase in height (−0.01 °C m −1 ). Meanwhile, a lapse rate of 0.015 and −0.02 °C −1 was used for stable and unstable conditions, respectively. Table 1 shows the initial boundary conditions of the model for CFD simulation. The speed of the exhaust fan of the CFD model for validation was applied by measuring the wind speed of the fan at nine points when the exhaust fan ran at maximum (Figure 8). To measure the air flow rate of a circular exhaust fan, a traverse method was applied at the duct surface using hot wire anemometers to measure the discharged speed. The measurement points (9 points) were selected based on ASHRAE (American Society of Heating, Refrigerating and Air-Conditioning Engineers) standards [59]. While the wind speed of the exhaust fan for the case study was analysed by applying the ventilation flow rate provided by the manufacturer to the CFD model, in the case of the Korean standard pig house used for the case study, ventilation rates of the pig house were calculated from For the computation of the CFD model, a pressure-based solver, which is widely applicable in flow analysis, was performed and the gravitational acceleration was set to −9.81 m s −2 . A standard operating pressure of 101,325 Pa was applied. The solution was assumed to have converged when the sum of residuals for all of the cells in the computational domain was less than 1 × 10 −6 for continuity and energy, and 1 × 10 −3 for other variables such as X, Y, Z-velocity, k, and ε.   Turbulence models selected from the validation test (Standard k-ε, RNG k-ε, Realisable k-ε, LES) * field measured data; ** computed data.
As mentioned, odour dispersion in the atmosphere may vary depending on weather and structural variables. Therefore, wind speed (0.5, 1.5, and 2.5 m s −1 ), wind direction (north, east, south, and west), and atmospheric stability (stable, neutral, and unstable) were considered as weather factors, while the roof-chimney exhaust fan, sidewall exhaust fan, and conditions for using all the exhaust fans were considered as structural factors. The case study using 108 simulation cases for the standard Korean pig house was conducted using the determined grid size and turbulence model to evaluate the effect of various factors on odour dispersion distance. Data were calculated and analysed at the level of 1 OU, the minimum level of odour that humans could recognize at 1.5 m from the ground level.

Statistical Analysis for Identifying Key Factors Influencing Odour Dispersion
For the case study results, statistical analyses were conducted using the Statistical Package for the Social Sciences (SPSS) program (ver. 23.0, IBM Inc., Armonk, NY, USA). From the CFD-computed odour dispersion distances, a statistical analysis was performed to determine the key factors influencing odour dispersion distance. Before the analysis, a multicollinearity test was first carried out to determine which among the selected factors would be included in the multiple regression analysis. Multicollinearity between variables can cause a problem because it can increase the variance of the coefficient estimates and make the estimates very sensitive to minor changes in the model. A common criterion for multicollinearity problems includes calculation of the tolerance and variance inflation factor (VIF), as shown in Equations (12) and (13). Tolerance is associated with each independent variable and ranges from 0 to 1, while VIF with a value exceeding 10 is often regarded as indicating multicollinearity. Accordingly, when high multicollinearity was found, a highly correlated independent variable from the model was removed.
where R 2 i is the coefficient of determination between i-th independent variable and other independent variables, and VIF is the reciprocal of tolerance and quantifies the severity of multicollinearity in an ordinary least squares regression analysis.
Multiple regression analysis is a statistical approach used to predict the value of the dependent variable based on the value of two or more independent variables. This tool shows how the variation of each factor influences the total value of the predicted result. Using this analysis, the relative influences of the considered independent variable can be analysed. To supplement the result of the regression analysis, the statistical method called the relative weight analysis (RWA) method was implemented to identify which among the independent variables had the greatest contribution to odour dispersion. The RWA method is a tool that addresses the multicollinearity problem and helped to calculate the relative importance of an independent variable. In this method, the variables transformation approach was carried out to create a new set of predictors that are orthogonal to one another. The general procedure determining the RWA was to create and describe the relationship between the original variables and the orthogonal variables through the regression equation. Relative weight can be calculated using the following equations (Equations (14) and (15)) based on the procedure established by Tonidandel et al. [61]: Relative Weight = Raw relative weight where R 2 is the coefficient of determination; y i is the y value for observation i;ŷ i is the predicted value of y for observation i; y is the mean of y value.

Field Experiment Results (Yeo et al., 2019)
Since the dispersion of odour in the surrounding area was greatly influenced by environmental factors, field measurement of various meteorological conditions such as wind speed and the temperature was first conducted. The prevailing wind was northerly winds based on analysis from a weather station near the farm (analysis period: 1997.1.1-2017.12.31). This result was similar to a portable weather station installed in the experimental pig house. The result of the prevailing wind direction was greatly influenced by the topographical location of the pig houses as the selected farm was surrounded by hills along the eastern and western parts. This type of topographical orientation reduced the impact of wind originating from the east and west part of the computational domain. The atmospheric stability during the conduct of the field experiment was generally stable with a temperature ranging from 15.7 to 33.5 • C and wind was 0.5 m s −1 , as shown in Table 2. In terms of odour concentration measured, a detailed summary can be found in the study of Yeo et al. [18]. The amount of odour concentration inside the pig house had higher concentration during September, with a value ranging from 30.0 to 448.1 OU m −3 compared to August, with values ranging only from 5.5 to 44.8 OU m −3 . The variation in measured odour concentration inside the pig house was due to the large volume of manure present in the slurry pit during September. During this month, the farm owner no longer disposes of the manure, since the rearing period for that season was almost due. The amount of odour concentration measured outside the pig house was significantly lower compared to inside the pig house. This result accounted for the dilution of concentrated odour with clean air outside the pig house when the inside odour comes out from the pig house facilities. Considering that the amount of odour concentrations measured inside the pig houses was highest at the end of the rearing period, and a distinct odour concentration trend downwind of the pig houses was found, the field-measured odour concentration on September 6 was therefore used as a reference value to validate the CFD model by comparing the CFD-computed results and the field-measured data. Moreover, the field experiment result showed that the measured odour concentration at the boundary of the pig farm falls below the minimum level of allowable OU m −3 value set by the Korean government, which was 15 OU m −3 . This means that the odour emission complied with the maximum OU m −3 limit for Korean livestock houses.

Validation of the CFD Simulation Model
The validation of the designed CFD model included the grid independence test, time step test, and turbulence test. The size of the grids must be well identified, since it contributed to the accuracy of the CFD-computed results. Since the size of the computational domain has a large diameter of 5.0 km, the grid independence test of the CFD model was carried out using four grid sizes (4.0, 10.0, 20.0, and 40.0 m). According to the designed grid, the number of grids in the computational domain ranged from 3.76 million to 8.73 million combined hexahedral and tetrahedral elements. The number of grid elements decreased when a larger grid size was used, whereas a greater number of grid elements were created when smaller grids sizes were used. After convergence was achieved, the CFD-computed results showed considerable differences in the predicted odour between the four grid sizes as summarized in Table 3. When the RMSE was analysed, the model which designed the grid with a size of 4.0 m had the lowest error, while the error was significant when the grid was larger than 10.0 m. The CFD model, which had a relatively large grid size, such as 20.0 and 40.0 m, showed imprecise odour concentration results when compared with the field-measured data, while both the 4.0 and 10.0 m grids calculated the odour concentration relatively close to the field-measured data. These grid sizes obtained the highest R 2 value of 0.999 and 0.998 and lowest RSME value of 1.06 and 1.09, respectively. However, for the grid sizes of 4.0 and 10.0 m, there was no significant difference between the CFD-computed results and the field-measured data. Therefore, considering efficient computation time, the grid size for the CFD model was determined to be 10.0 m. The simulation results according to the four time steps are shown in Table 4. The use of 10and 20-s time intervals confirmed that the RSME (5.027 and 9.711) were relatively large compared to smaller time steps. In contrast, almost identical results were obtained between 1 and 5 s, with an RMSE level of 1.083 and 1.093, respectively. Based on this, it can be concluded that the accuracy of the simulation results using the simulation model is inversely proportional to the size of the time step. Considering the excessive calculation time when using the 1 s time step, the 5 s time step was selected as the most economic time step value to be used in the simulation. The results related to the turbulence model are summarized in Table 5. Considering the different turbulence models commonly used in odour dispersion studies, in this research, a total of four turbulence models were tested including the standard k-ε model, the RNG k-ε model, the standard k-ω model, and the LES model. The comparison of four turbulence models showed lower RMSE values for the standard k-ε and the RNG k-ε models. Statistical analysis with lower RSME values meant that the CFD-computed results had some discrepancy with the field-measured data. The results showed that the CFD-computed results using the standard k-ε model have the lowest RSME value of 1.093 compared with the other three models, with RSME values of 1.385, 1.696, and 2.165 for the RNG k-ε, standard k-ω, and LES models, respectively. Considering the statistical analysis results, the standard k-ε model was found to be the most appropriate turbulence model to be used.

Effect of Ventilation Types on Odour Dispersion
The odour dispersion distance when the sidewall exhaust fan or roof-chimney exhaust fan was used in the pig house was calculated using CFD. The CFD-computed results for the odour dispersion distance are summarized in Table 6. The CFD-computed results according to the cases ranged from a minimum of 78.9 m to a maximum of 1488.1 m. When using a sidewall exhaust fan, the odour dispersion ranged from 147.8 to 1488.1 m, depending on the case being simulated. On the other hand, the roof-chimney exhaust fan showed a minimum of 78.9 m to a maximum dispersion distance of 900.6 m. On average, the odour dispersion distance was 439.4 m when the exhaust fan on the sidewall was used, while the odour dispersion distance was 250.7 m when the roof-chimney exhaust fan was used. Specifically, the highest dispersion distance was computed when the atmospheric stability was neutral, the wind speed was low, the southerly wind was used, and the sidewall exhaust fan was operated. The lowest odour dispersion distance was computed when the atmospheric stability was unstable, the wind speed was high, with a westerly wind, and the roof-chimney exhaust fan was used. Under the same environmental conditions, there was a difference of 41.3% (maximum: 61.7%, minimum: 16.8%) dispersion distance, depending on the physical environment (ventilation method; sidewall exhaust fan; roof-chimney exhaust fan). The trend of odour dispersing outside the pig house according to the ventilation method under the condition of the north wind direction, wind speed of 1.5 m s −1 , and stable condition in the atmospheric stability is shown in Figure 9. Given that the initial amount of odour emission from exhaust fans is almost identical between the ventilation systems (577.5 OU m −3 ), it can be seen that the odour concentration emitted from the pig house using the sidewall exhaust fans was higher near the surrounding area of the pig houses. On the other hand, when the roof-chimney exhaust fan was used, the odour discharged through the exhaust fan has a lower odour concentration. Moreover, when the sidewall exhaust fan and the chimney exhaust fan were used simultaneously, the spreading distance of odour outside the pig house was 92.0 to 928.2 m depending on wind speed, wind direction, and atmospheric stability conditions. The odour dispersion distance was estimated to be 5.4% higher than that of the exhaust using only the roof-chimney exhaust fan. The emission of odour at the sidewall exhaust fan contributed to the increase in the distribution of the odour along downwind.

Effect of Atmospheric Stability in Odour Dispersion
The condition of the atmosphere is an important factor in determining the dispersion of odour in the atmosphere. The dispersion distance of odour emitted from pig houses when the condition of the atmosphere was stable was calculated from 112.3 to 1427.4 m. The average dispersion distance was 389.3 m. On the other hand, when the atmospheric stability was unstable, the odour emitted from the pig houses dispersed at a maximum distance of 485.3 m, a minimum of 78.9 m, and an average of 190.0 m. It decreased by 51.3% on average due to the change in atmospheric stability (stable to unstable). This was because, when the atmosphere was stable, the influence of the vertical momentum of air was relatively small resulting in higher odour dispersion distance. However, even if the atmospheric stability was stable, the dispersion distance of the odour was not extended. It was confirmed that the dispersion distance was greatly affected by the presence or absence of the obstacles toward the direction of the wind and wind speed. In contrast, when the atmospheric stability was unstable, vertical mixing and dispersion were very active, so the odour concentration was low near the ground, and the dispersion distance was short as well ( Figure 10). In addition, there was no significant difference between stable and neutral conditions in the stability of the atmosphere. The same trend result was found in the odour dispersion modelling result [20,45,60]. Representative figures showing odour dispersion according to the atmospheric stability are shown in Figure 10, where the sidewall exhaust fan was used as a ventilation system, wind speed of 0.5 m s −1 , and wind direction originating from the south part of the computational domain. As can be seen, the transport of odours was prominent and continuous under stable conditions, while the unstable condition was short. Furthermore, a maximum dispersion width of 155.24 m was calculated when the atmospheric stability was unstable. This was 17.42% higher than the maximum width when the atmospheric stability was stable (132.8 m).

Effect of Atmospheric Stability in Odour Dispersion
The condition of the atmosphere is an important factor in determining the dispersion of odour in the atmosphere. The dispersion distance of odour emitted from pig houses when the condition of the atmosphere was stable was calculated from 112.3 to 1427.4 m. The average dispersion distance was 389.3 m. On the other hand, when the atmospheric stability was unstable, the odour emitted from the pig houses dispersed at a maximum distance of 485.3 m, a minimum of 78.9 m, and an average of 190.0 m. It decreased by 51.3% on average due to the change in atmospheric stability (stable to unstable). This was because, when the atmosphere was stable, the influence of the vertical momentum of air was relatively small resulting in higher odour dispersion distance. However, even if the atmospheric stability was stable, the dispersion distance of the odour was not extended. It was confirmed that the dispersion distance was greatly affected by the presence or absence of the obstacles toward the direction of the wind and wind speed. In contrast, when the atmospheric stability was unstable, vertical mixing and dispersion were very active, so the odour concentration was low near the ground, and the dispersion distance was short as well ( Figure 10). In addition, there was no significant difference between stable and neutral conditions in the stability of the atmosphere. The same trend result was found in the odour dispersion modelling result [20,45,60]. Representative figures showing odour dispersion according to the atmospheric stability are shown in Figure 10, where the sidewall exhaust fan was used as a ventilation system, wind speed of 0.5 m s −1 , and wind direction originating from the south part of the computational domain. As can be seen, the transport of odours was prominent and continuous under stable conditions, while the unstable condition was short. Furthermore, a maximum dispersion width of 155.24 m was calculated when the atmospheric stability was unstable. This was 17.42% higher than the maximum width when the atmospheric stability was stable (132.8 m). wind speed, the dispersion distance was reduced by an average of 61.7% compared with 0.5 m s −1 . As the wind speed increased, the dispersion distance of odour emitted from pig houses decreased. This trend was mentioned in the studies of Chastain and Wolak [22], Guo et al. [62], Xing [63], and Heinemann and Wahanik [64]. Figure 11 represented CFD-computed results of odour dispersion for various wind speeds under the conditions that pig houses used the sidewall exhaust fan, the wind At 2.5 m s −1 wind speed, the dispersion distance was reduced by an average of 61.7% compared with 0.5 m s −1 . As the wind speed increased, the dispersion distance of odour emitted from pig houses decreased. This trend was mentioned in the studies of Chastain and Wolak [22], Guo et al. [62], Xing [63], and Heinemann and Wahanik [64]. Figure 11 represented CFD-computed results of odour dispersion for various wind speeds under the conditions that pig houses used the sidewall exhaust fan, the wind direction was southerly wind, and the atmospheric stability was stable. Indeed, the odour dispersed farther at high wind speeds, but the dispersed odour was more diluted with external air, so the odour dispersion distance seemed to decrease as the odour concentration dropped below the human perception limit. In addition, when the odour was continuously emitted from the pig houses and the wind was constantly blown at a specific wind speed and in a specific wind direction, the smaller the wind speed, the larger the horizontal dispersion width (left and right) from the outlet, as shown in Figure 11. On the contrary, when the wind speed around the pig houses increased, the odour dispersion to the horizontal direction did not occur much. The dispersion width at the lowest wind speed was 132.71 m; this is higher by 8.1% and 179.5% compared with the dispersion width when wind speed was 1.5 m s −1 (122.8 m) and 2.5 m s −1 (47.5 m), respectively.

Effect of Wind Direction in Odour Dispersion
Wind direction is also an important weather factor that determines the direction of odour dispersion. This wind direction is greatly influenced by the surrounding topography. Therefore, when considering the wind direction, it must be explained along with the complex terrain. The average dispersion distances for wind directions (north, east, south, and west) were calculated to be 369.0, 197.8, 470.5, and 235.1 m, respectively. The CFD-computed results showed that the odour dispersion distances were longer with a southerly wind and with the northerly wind at any atmospheric stability, ventilation type, and wind speed because there are no obstacles with high altitude, such as mountains and hills on the windward and leeward. When the northerly wind was generated, the average odour distance was 21.6% lower than that of the southerly wind. The reason for this was that there were no obstacles to the north from the pig house, but the mountain area was located at about 400 m away from the pig house in the south, which may affect the dispersion of odour. In the same way, the wind speed near the ground of southerly wind was lower than the speed of northerly wind because the wind did not recover from the height near the farm due to the mountain located at the southern part of the domain. Additionally, when the westerly wind was formed, the odour dispersion distance was increased by 18.8% compared with the easterly wind. This result was influenced by the steep slope of the mountain along the western direction (approximately 49.2%) compared with the slope of the mountain along the eastern direction (38.5%), which limits the transports of odours farther from the pig farm. Figure 12 showed the vector field around experimental pig houses under easterly wind conditions. As shown in the figure, when the wind originated from the eastern part of the computational domain, the wind passed through a mountain first before reaching the pig houses.
Agronomy 2020, 10, x 19 of 27 direction was southerly wind, and the atmospheric stability was stable. Indeed, the odour dispersed farther at high wind speeds, but the dispersed odour was more diluted with external air, so the odour dispersion distance seemed to decrease as the odour concentration dropped below the human perception limit. In addition, when the odour was continuously emitted from the pig houses and the wind was constantly blown at a specific wind speed and in a specific wind direction, the smaller the wind speed, the larger the horizontal dispersion width (left and right) from the outlet, as shown in Figure 11. On the contrary, when the wind speed around the pig houses increased, the odour dispersion to the horizontal direction did not occur much. The dispersion width at the lowest wind speed was 132.71 m; this is higher by 8.1% and 179.5% compared with the dispersion width when wind speed was 1.5 m s −1 (122.8 m) and 2.5 m s −1 (47.5 m), respectively.

Effect of Wind Direction in Odour Dispersion
Wind direction is also an important weather factor that determines the direction of odour dispersion. This wind direction is greatly influenced by the surrounding topography. Therefore, when considering the wind direction, it must be explained along with the complex terrain. The average dispersion distances for wind directions (north, east, south, and west) were calculated to be 369.0, 197.8, 470.5, and 235.1 m, respectively. The CFD-computed results showed that the odour dispersion distances were longer with a southerly wind and with the northerly wind at any atmospheric stability, ventilation type, and wind speed because there are no obstacles with high altitude, such as mountains and hills on the windward and leeward. When the northerly wind was generated, the average odour distance was 21.6% lower than that of the southerly wind. The reason for this was that there were no obstacles to the north from the pig house, but the mountain area was located at about 400 m away from the pig house in the south, which may affect the dispersion of odour. In the same way, the wind speed near the ground of southerly wind was lower than the speed of northerly wind because the wind did not recover from the height near the farm due to the mountain located at the southern part of the domain. Additionally, when the westerly wind was formed, the odour dispersion distance was increased by 18.8% compared with the easterly wind. This result was influenced by the steep slope of the mountain along the western direction (approximately 49.2%) compared with the slope of the mountain along the eastern direction (38.5%), which limits the transports of odours farther from the pig farm. Figure 12 showed the vector field around experimental pig houses under easterly wind conditions. As shown in the figure, when the wind originated from the eastern part of the computational domain, the wind passed through a mountain first before reaching the pig houses.

Statistical Analysis for Odour Dispersion
Multiple regression analyses were performed to estimate the odour dispersion distance depending on ventilation type, atmospheric stability, wind speed, and wind direction. In this study, the atmospheric stability, ventilation type, and wind direction were nominal variables, while both wind speed and the odour dispersion distance were scale variables. Therefore, a total of seven dummy variables were adopted in the statistical analysis of the CFD-computed results. In the analysis, dummy variable (a) was used to represent the statistical difference between the effect of sidewall exhaust fan and roof-chimney exhaust fan; dummy variable (b) was used to represent the statistical difference between the effect of sidewall exhaust fan and the combined sidewall-roofchimney exhaust fan. The dummy variable (c) was used to evaluate the statistical difference in the effect between stable atmospheric stability and neutral atmospheric stability; dummy variable (d) was used to evaluate the statistical difference in the effect between stable atmospheric stability and unstable atmospheric stability. To represent the statistical difference to north wind direction with other wind directions, dummy variables (e), (f), and (g) representing east, south, and west directions were used. The independent variable used includes speed, dummy variable (a), (b), (c), (d), (e), (f), and (g), while the dependent variable was the odour dispersion distance.
In the multicollinearity test, which was first conducted prior to the statistical analysis, the calculated VIF values were below 10 for each independent variable, indicating no significant relationship between variables. To determine the appropriate model to predict odour dispersion distance, a backward elimination process was implemented. As shown in Table 7, the R 2 and the adjusted R 2 value for the model with a lesser number of independent variables were lower compared with the models considering a greater number of independent variables. The linear regression equation derived is shown in Table 7 with an R 2 value and adjusted R 2 value of 0.685 and 0.663, respectively, for odour dispersion distance. In the regression analysis, the change between stable and neutral atmospheric stability (dummy variable c) showed a very small t-value of −0.409. In addition, the calculated p-value for the dummy variable (c) was found to be insignificant at p < 0.05, with a value of 0.683. This means that the change in dummy variable (c) is not associated with the changes in the dependent variable because no significant difference was found between the cases having an atmospheric stability of stable and neutral. Thus, the dummy variable representing the effect of stable and neutral conditions (dummy variable c) was excluded from the independent variables and regression analysis was recalculated.

Statistical Analysis for Odour Dispersion
Multiple regression analyses were performed to estimate the odour dispersion distance depending on ventilation type, atmospheric stability, wind speed, and wind direction. In this study, the atmospheric stability, ventilation type, and wind direction were nominal variables, while both wind speed and the odour dispersion distance were scale variables. Therefore, a total of seven dummy variables were adopted in the statistical analysis of the CFD-computed results. In the analysis, dummy variable (a) was used to represent the statistical difference between the effect of sidewall exhaust fan and roof-chimney exhaust fan; dummy variable (b) was used to represent the statistical difference between the effect of sidewall exhaust fan and the combined sidewall-roof-chimney exhaust fan. The dummy variable (c) was used to evaluate the statistical difference in the effect between stable atmospheric stability and neutral atmospheric stability; dummy variable (d) was used to evaluate the statistical difference in the effect between stable atmospheric stability and unstable atmospheric stability. To represent the statistical difference to north wind direction with other wind directions, dummy variables (e), (f), and (g) representing east, south, and west directions were used. The independent variable used includes speed, dummy variable (a), (b), (c), (d), (e), (f), and (g), while the dependent variable was the odour dispersion distance.
In the multicollinearity test, which was first conducted prior to the statistical analysis, the calculated VIF values were below 10 for each independent variable, indicating no significant relationship between variables. To determine the appropriate model to predict odour dispersion distance, a backward elimination process was implemented. As shown in Table 7, the R 2 and the adjusted R 2 value for the model with a lesser number of independent variables were lower compared with the models considering a greater number of independent variables. The linear regression equation derived is shown in Table 7 with an R 2 value and adjusted R 2 value of 0.685 and 0.663, respectively, for odour dispersion distance. In the regression analysis, the change between stable and neutral atmospheric stability (dummy variable c) showed a very small t-value of −0.409. In addition, the calculated p-value for the dummy variable (c) was found to be insignificant at p < 0.05, with a value of 0.683. This means that the change in dummy variable (c) is not associated with the changes in the dependent variable because no significant difference was found between the cases having an atmospheric stability of stable and neutral. Thus, the dummy variable representing the effect of stable and neutral conditions (dummy variable c) was excluded from the independent variables and regression analysis was recalculated. The results of the regression analysis in Table 8 showed that all the independent variables are significant at a confidence level of 95%. Furthermore, the wind speed has the highest t-value and lowest p-value of −9.057 and 1.15 × 10 −14 , respectively. This indicated that an abrupt change in wind speed can substantially affect the dispersion of odour along downwind. The influence of atmospheric stability was also found to have t and p-values of −6.367 and 5.90 × 10 −9 . The vertical movement of odorous parcels caused by the stability and instability of the atmospheric stability has a large impact on the dispersion of odour. Meanwhile, dummy variables (a) and (b) representing the statistical difference between the effect of using a sidewall exhaust fan and other ventilation systems showed almost identical t-values of −5.415 (dummy variable a) and −5.026 (dummy variable c). The odour dispersion distance can be predicted using the model shown below (Equation (16)  To supplement the regression analysis in identifying the key factors affecting the dispersion of odour outside the pig farm, a relative weight analysis on the independent variable was used. The calculated R 2 was 0.44 and the confidence level was set at 95%. The result of the analysis summarized in Table 9 showed a similar trend in terms of independent variables' ranking with regard to contribution to odour prediction. In the RWA method, the wind speed was the most influential factor that contributes to odour dispersion, followed by the atmospheric stability, ventilation type, and wind direction, with rescaled relative weight values of 58.2%, 23.2%, 17.9%, and 0.7%, respectively.

Discussion
Unlike typical results showing higher accuracy in a simulation using the LES turbulence model, the LES model showed a higher error compared with other turbulence models. This may be due to the y+ value related to mesh resolution. In the case of the LES turbulence model, the y+ value should be designed to be 1 or less. However, it was difficult to satisfy this criterion when considering the size of the target computational domain in this study. Therefore, unlike the RANS models (standard k-ε model, RNG k-ε model, and standard k-ω model), the LES model did not fully meet the criteria for the y+ value. Although the y+ design criterion was not met, the CFD-computed result using LES models was expected to have an accurate result. Thus, in this study, the LES model was used as a turbulence model for validating the three-dimensional model. However, contrary to the expected results, the LES model showed a rather low accuracy. If the y+ criterion is not met, it is thought that it can show lower accuracy than the RANS model in terms of accuracy. Among the RANS model, the standard k-ε turbulence model has been widely used in various CFD studies due to its relative simplicity, robustness, ease of usage, lesser computational requirements, and its applicability to various flows. Therefore, the standard k-ε turbulence model was applied in this study.
The difference in the odour dispersion distance according to the ways (ventilation method; sidewall exhaust fan; roof-chimney exhaust fan) of applying the exhaust fan of the pig house varied from 16.8% to 61.7%. The reason was that there was a difference in the odour concentration at the height of the receptor according to the effect of the shield installed on the roof-chimney exhaust fan and the height difference of the outlet. when the sidewall exhaust fan was used, the odorous air emitted directly travelled away from that source at a constant rate with no vortex formation (Figure 13a). In the case of a roof-chimney exhaust fan, a vortex was formed around the shield of the exhaust fan reducing the concentration of the odour emitted ( Figure 13b). The vortex created initiated the mixing of the air with the odorous air emitted from the pig house causing lower dispersion distance of odour.
There was no significant difference between stable and neutral conditions in the stability of the atmosphere. The same trend result was found in the odour dispersion modelling result [20,45,60]. The not significant difference between the result of neutral and stable atmospheric stability is due to equal or lower temperature of odorous parcel compared with the surrounding temperature allowing the pollutants to remain on the ground. For the dispersion width, the difference can be explained by the turbulent mixing that happened at unstable atmospheric stability commencing the movement of odorous compounds to move in all directions. Thus, the odorous air was easily diluted with air and became less detectable in the leeward direction. In addition, the odour dispersed farther at high wind speeds, but the dispersed odour was more diluted with external air, so the odour dispersion distance seemed to decrease as the odour concentration dropped below the human perception limit.
Surrounding complex terrain and mountain affected the dispersion distance; since the pig house was less than 100 m away from the foot of the hill along with the eastern domain, the wind was not strong enough to lift and transport the odorous substances away from the source. Specifically, when the wind from this direction flowed through the mountains, the wind accelerated uphill and decelerated near the downhill slopes. After the wind passed through the mountains, its momentum was reduced, resulting in lower wind speed. This causes lower dispersion at easterly and westerly wind directions.

Conclusions
The rapid intensification in pig production has led to an increase in livestock pollutant emissions. For this reason, attempts are being made to study the influence of various environmental conditions on odour dispersion. With the advancement of technology and computer software, dispersion models were one of the most common methods to visualize and estimate odour occurrences. Therefore, a three-dimensional model containing a complex terrain was developed in this study to realise odour dispersion originating from pig houses.
The CFD-computed results showed that odour dispersion distance at the level of 1 OU m −3 ranged from 78.9 to 1488.1 m. When using an exhaust fan on the sidewall, the odour dispersion ranged from 147.8 to 1488.1 m, depending on the simulated environmental conditions. On the other hand, the roof-chimney exhaust fan showed a minimum of 78.1 to 900.6 m. There was an average difference of 41.3% (maximum: 61.7%, minimum: 16.8%) dispersion distance, depending on the physical environment (ventilation method; sidewall exhaust fan; roof-chimney exhaust fan). Moreover, when the sidewall exhaust fan and the chimney exhaust fan were used simultaneously, the odour dispersion distance was estimated to be 5.4% larger than that of the exhaust using only the roof-chimney exhaust fan. It was shown that the application of a sidewall exhaust fan can

Conclusions
The rapid intensification in pig production has led to an increase in livestock pollutant emissions. For this reason, attempts are being made to study the influence of various environmental conditions on odour dispersion. With the advancement of technology and computer software, dispersion models were one of the most common methods to visualize and estimate odour occurrences. Therefore, a three-dimensional model containing a complex terrain was developed in this study to realise odour dispersion originating from pig houses.
The CFD-computed results showed that odour dispersion distance at the level of 1 OU m −3 ranged from 78.9 to 1488.1 m. When using an exhaust fan on the sidewall, the odour dispersion ranged from 147.8 to 1488.1 m, depending on the simulated environmental conditions. On the other hand, the roof-chimney exhaust fan showed a minimum of 78.1 to 900.6 m. There was an average difference of 41.3% (maximum: 61.7%, minimum: 16.8%) dispersion distance, depending on the physical environment (ventilation method; sidewall exhaust fan; roof-chimney exhaust fan). Moreover, when the sidewall exhaust fan and the chimney exhaust fan were used simultaneously, the odour