Sustainable Development Optimization of a Plant Factory for Reducing Tip Burn Disease

: It is generally believed that stable airﬂow can effectively reduce tip burn, a common lettuce plant disease in closed plant factories that severely restricts the sustainable development of these factories. This study aims to achieve stable airﬂow in the cultivator by zoning the seedling and growth stage crops and installing differential fans, while ensuring comprehensive quality. In this study, a three-dimensional simulation plant factory model was created to simulate the airﬂow inside the cultivator, taking crop shading and heat dissipation from LED light sources into account. Experiments on photosynthetic physiology and airﬂow were used to determine environmental thresholds for crop growth, which were then used as CFD boundary conditions. After adopting the optimized cultivation model, the comprehensive quality of lettuce increased by 22.28% during the seedling stage, and the tip burn rate decreased to 26.9%; during the growth stage, the comprehensive quality increased by 25.72%, and the tip burn rate decreased to 23.2%. The zoning optimization cultivation method and differential fan arrangement used in this study to improve the airﬂow ﬁeld of plant factories provide new ideas and reliable theoretical support for plant factories to combat lettuce tip burn disease.


Introduction
With the outbreak and spread of COVID-19 in recent years, sustainable development of agriculture, especially plant factories, has attracted increasing attention [1][2][3][4]. Plant factories can improve agricultural production efficiency by regulating environmental factors, enabling agricultural production to break free of the shackles of natural conditions and the effects of the outside world [5,6]. Lettuce is the most common crop in plant factories because the cultivation period is short, with less energy being required to produce the final product. Lettuce tip burn, which is also the most prevalent disease in plant factories, will result in the yellowing of the upper tip edges and the curling of leaves, significantly impacting crop quality and yield [7]. Scholars have conducted a significant amount of research on the pathogenesis of lettuce tip burn in the growth stage. Studies have shown that the pathogenesis of lettuce tip burn in the growth stage is due to the inhibition of transpiration in the high humidity environment of plant factories, which leads to the inability of Ca 2+ to be transported smoothly from the culture solution to the leaves through the roots [8][9][10][11].
Air conditioning systems are used to ventilate plant factories, which are highly enclosed structures. However, the airflow inlet and outlet locations are fixed for established plant factories. Shading of the airflow by the cultivation frame and plant cultivation location, along with heat dissipation from the light source, can affect the airflow field distribution [12][13][14]. After long-term planting, uneven airflow can result in issues such as tip burn, uneven production, and increased energy consumption [15,16]. When lettuce varieties and cultivation methods were identified, Kitaya and Ryohei et al. confirmed that stable and uniform airflow is an effective tool for resolving the problem of lettuce tip burn [17,18]. Computational fluid dynamics (CFD) is an efficient tool for reasonably simulating complex physical phenomena and analyzing environmental uniformity in controlled environments [19]. Ying et al. designed a three-row outlet duct with fans for downward ventilation, and the CFD simulations indicated that the device could effectively improve the airflow distribution in the crop canopy [20]. Benyamin et al. developed a computational fluid dynamics model to simulate the fluid flow in an indoor vertical farming system using photosynthesis. The case with the inlet at the bottom of the short side and the exit at the top of the long side, operating at 0.3 kg/s, demonstrated the desired lower average temperature compared to other vent locations operating at 0.4 or 0.5 kg/s [21]. However, these designs must be completed prior to the construction of the plant factories. For existing plant factories, it is evident that the location of the vents and the lack of monitoring and verification of the actual airflow field cannot be changed [22,23]. Kozai et al. concluded that LED heat dissipation affects the uniformity of the airflow field, particularly in the case of different light source arrangements in the plant factories [24]. Hui F et al. considered crop resistance and LED heat dissipation, designed ventilation ducts with different apertures, and simulated plant canopy airflow. The simulation results revealed that ventilation ducts with 12 mm apertures could increase the plant airflow comfort zone to 70.3% [25]. For plant factories with fixed air inlet and outlet locations, appropriate cultivation methods based on the inherent airflow distribution are more likely to reduce production costs [26,27]. Baek and Ahmed et al. proposed a CFD-based internal differential fan control simulation to improve the cultivation environment inside a plant factory. The uniformity of the airflow field and the comprehensive quality of the crop were significantly enhanced based on simulation results [28,29]. Li Kun used interlayer cool airflow (ILCA) to introduce room air into the internal plant canopies, resulting in a 41.7% reduction in root growth, without negatively impacting lettuce crop yield [30]. However, there is a lack of research regarding optimizing the crop canopy airflow field from the perspective of cultivation location distribution; most studies on optimizing such methods are focused on crop cultivation methods [31].
This study aimed to reduce lettuce tip burn by controlling airflow while maintaining comprehensive plant quality. This study utilized the analytic hierarchy process (Ahp) to analyze the results of lettuce photosynthetic physiology experiments to determine plant light source arrangement. Airflow treatment tests were performed to determine the airflow speed threshold for reducing the incidence of tip burns, and the results were used as simulation boundary conditions. Based on the above findings, we proposed zoning cultivation into high-and low-speed areas, based on the airflow distribution of the plant factory, and implemented CFD simulation. Considering the crop's shading, we conducted aerodynamic experiments to determine the parameters of the lettuce's porous medium during the seedling and growing stages. This study arranged differential fans at the side ends of the cultivators and adjusted the differential speed according to the airflow field of each cultivator to reduce the uneven airflow field caused by LED heat dissipation. The study's results indicated that the zoning cultivation and the arrangement of differential fans effectively mitigated tip burn disease in the plant factory and ensured the lettuce's comprehensive quality. For the study of sustainable development of plant factories, our study offers a novel concept for low-cost management of tip burn disease in plant factories, avoiding the modification of the plant factories' inherent building patterns.

Materials and Conditions
The photosynthetic physiology experiment was conducted in an artificial climate incubator (MRC-550C). The climate incubator was a three-layer spacer-type construction, with adjustable spacing between the spacers. The environmental conditions were controlled, including airflow, temperature, humidity, and light cycle. The light source radiation form was vertical irradiation, the ventilation mode was horizontal or vertical, the temperature range was 0-50 • C, the maximum light intensity was 500 µmol/m 2 s, and the humidity adjustment range was 30-95% [32,33]. The LED light used a Guixiang full-spectrum white LED light source and a ZPDT812 adjustable red and blue LED light source, with a light period of 16 h/d. The lettuce cultivar chosen was Italian bolting lettuce. After soaking for 7 h, the seeds were wrapped in wet gauze for 3-4 days. The seedlings were cultured in water with a seedling sponge after budding. The seedlings with the same growth rate were selected and transplanted to each treatment group (n = 24) once they exhibited three true leaves. The nutrient solution follows the Japanese garden test nutrient solution (pH 6.0-6.5, EC 1.0 ms cm −1 ), the ambient temperature was maintained at 21 • C, and the relative humidity was set at 80%.

Test Methods
The photosynthetic physiology experiment was conducted on seedling and growth stage lettuce with the test schedule illustrated in Figure 1. There were six treatment groups for the seedling and growth stage lettuce photosynthetic physiology experiments, respectively. The number of samples analyzed per treatment group and time point was n = 24. The photosynthetic physiology experiments were treated with each light mode, as shown in Table 1, where WL is the full-spectrum white LED light source, and RB is the red and blue light source, and the different coefficients of the RB prefix represent the ratio of red and blue light sources [34]. The lettuce's growth and physiological indexes in the seedling stage were measured 22 days after sowing. After the seedling stage, lettuce with consistent growth without tip burn disease was selected for the growing stage treatment. Meanwhile, lettuce's growth and physiological indexes in the growth stage were measured 36 days after sowing.

Plant Growth and Physiological Index Measurements
In this section, eight physiological growth indexes (leaf area, tip burn rate, height, fresh weight, Spad, soluble sugar, ascorbic acid (AsA), and nitrate) were sele to constitute the evaluation index system of lettuce. Plant fresh weight for quantit analysis was measured by a ten-thousandth electronic analytical balance. Plant heigh measured by an electronic vernier caliper. Lettuce leaves with more than three br spots are considered to have tip burn disease, and the tip burn rate is the percentage o burn plants in the treatment group. The leaf area of the plants was measured by a Ya leaf area instrument, the AsA content was determined by molybdenum blue colorim the soluble sugars content was measured using the anthrone method, and the nitrate tent was determined by salicylic sulfuric acid colorimetry [35].

Data Processing and Analysis
Data management and analysis were performed using SPSS 26.0. Significant d ence analysis was performed using Duncan's multiple range test, and differences considered statistically significant at the 5% probability level.   CK1 is the seedling control group; CK2 is the growth stage control group; WL is the white full-spectrum LED lamp; RB is the red and blue light source, and the different coefficients of RB prefix represent the ratio of red and blue light sources.

Plant Growth and Physiological Index Measurements
In this section, eight physiological growth indexes (leaf area, tip burn rate, plant height, fresh weight, Spad, soluble sugar, ascorbic acid (AsA), and nitrate) were selected to constitute the evaluation index system of lettuce. Plant fresh weight for quantitative analysis was measured by a ten-thousandth electronic analytical balance. Plant height was Sustainability 2023, 15, 5607 4 of 16 measured by an electronic vernier caliper. Lettuce leaves with more than three brown spots are considered to have tip burn disease, and the tip burn rate is the percentage of tip burn plants in the treatment group. The leaf area of the plants was measured by a Ya-Xin leaf area instrument, the AsA content was determined by molybdenum blue colorimetry, the soluble sugars content was measured using the anthrone method, and the nitrate content was determined by salicylic sulfuric acid colorimetry [35].

Data Processing and Analysis
Data management and analysis were performed using SPSS 26.0. Significant difference analysis was performed using Duncan's multiple range test, and differences were considered statistically significant at the 5% probability level.

Analytic Hierarchy Process (Ahp)
The comprehensive evaluation system of the comprehensive quality of the lettuce is composed of the experimental indexes in Figure 2 (leaf area, tip burn rate, plant height, fresh weight, Spad, soluble sugar, ascorbic acid (AsA) and nitrate) [36][37][38]. The comprehensive quality of lettuce is evaluated by Ahp fuzzy comprehensive evaluation, and the Wi of each index is shown in to constitute the evaluation index system of lettuce. Plant fresh weight for quantitat analysis was measured by a ten-thousandth electronic analytical balance. Plant height w measured by an electronic vernier caliper. Lettuce leaves with more than three brow spots are considered to have tip burn disease, and the tip burn rate is the percentage of burn plants in the treatment group. The leaf area of the plants was measured by a Ya-X leaf area instrument, the AsA content was determined by molybdenum blue colorimet the soluble sugars content was measured using the anthrone method, and the nitrate co tent was determined by salicylic sulfuric acid colorimetry [35].

Data Processing and Analysis
Data management and analysis were performed using SPSS 26.0. Significant diff ence analysis was performed using Duncan's multiple range test, and differences w considered statistically significant at the 5% probability level.

Analytic Hierarchy Process (Ahp)
The comprehensive evaluation system of the comprehensive quality of the lettuce composed of the experimental indexes in Figure 2 (leaf area, tip burn rate, plant heig fresh weight, Spad, soluble sugar, ascorbic acid (AsA) and nitrate) [36][37][38]. The comp hensive quality of lettuce is evaluated by Ahp fuzzy comprehensive evaluation, and t Wi of each index is shown in Table 2

Airflow Experiment
The seedling lettuce and growth stage lettuce were transferred to an artificial climate incubator (MRC-550C) for a different airflow. The environmental factors were set in the same way as in the photosynthetic physiology experiment in Section 2.1. The lighting pattern setting for the airflow test needed to be based on the results of the plant light and physiology experiments in Section 2.1. The airflow test time period was the same as the photoperiod. The airflow direction in the climate incubator is horizontal or vertical, from top to bottom, respectively. Our previous research experience showed airflow speeds ranging from 0.1 to 0.9 m/s in the seedling stage and from 0.5 to 1.4 m/s in the growth stage [39]. The tip burn rate of seedling stage lettuce was measured at 22 days after sowing; the tip rate of growth stage lettuce was measured at 36 days after sowing.

Porous Medium Model of Lettuce
The simulation of plant factory fluid must account for plant occlusion of airflow. Therefore, aerodynamic tests on the crop are needed to obtain the parameters of the crop porous medium. After establishing plants as porous media, the resistance coefficients in all directions can be artificially defined to replace the resistance of solids to fluids in porous media, i.e., a momentum sink related to velocity is added, and its expression equation is: where v is the velocity value, D is the viscous drag coefficient, which in turn can be rewritten as 1 k , k is the permeability, C 2 is the inertial drag coefficient, momentum sink S i acts on the fluid to produce a pressure gradient ∇P, ∇P = −S i * l, and l is the thickness of the porous media domain [40].
From the above relationship between the pressure drop ∇P and the flow velocity v, we obtain a quadratic equation Y = Ax 2 + Bx, and Equation (1) can be rewritten as: The aerodynamic test was conducted in Jilin University's Bionic Key Laboratory's low-speed airflow tunnel. The seedling lettuce porous medium consisted of 3 × 4 crops; the growth period lettuce porous medium consisted of 2 × 3 crops. The monitoring points of pressure drop and the direction of tunnel airflow are shown in Figure 3 [41]. The dimension of the xyz vector direction of lettuce is 360 mm × 250 mm × 200 mm, and the dimension of the xyz vector direction of lettuce in seedling period is 360 mm × 250 mm × 100 mm. l is the length of the lettuce porous medium along the airflow direction in the airflow test, i.e., the xyz direction vector value. The pressure drop of the lettuce canopy was measured in three xyz directions at different airflow speeds. For aerodynamic tests on porous media, velocity ranges of 0.1-1.3 m/s were set for the seedling stage and 0.5-1.5 m/s for the growth stage. According to Equation (2), the polynomial coefficients D and C 2 could be solved. The leaf area density of the seedling stage lettuce was 8.25 m 2 /m 3 , and the leaf area density of growth stage lettuce was 32.3 m 2 /m 3 . The air density ρ is 1.205 kg/m 3 , and the air viscosity µ is 17.9 × 10 −6 Pa·s. The three-dimensional model was imported into the star ccm+ module, and the spatial discretization of the flow field in the greenhouse was carried out using a cut-body grid (trim). At a spatial grid volume of 1850 w, the volume change quality of the grid cells was 0.938, thus meeting the solution requirements. on porous media, velocity ranges of 0.1-1.3 m/s were set for the seedling stage and 0.5 1.5 m/s for the growth stage. According to Equation (2), the polynomial coefficients D an C2 could be solved. The leaf area density of the seedling stage lettuce was 8.25 m 2 /m 3 , an the leaf area density of growth stage lettuce was 32.3 m 2 /m 3 . The air density is 1.20 kg/m 3 , and the air viscosity is 17.9 × 10 −6 Pa·s.

Geometric Model of the Plant Factory
This experiment was carried out at Jilin University's artificial light plant factory. Th plant factory's specification is 5575 mm × 3950 mm × 3175 mm, as shown in the Figure 4 The three-dimensional model was imported into the star ccm+ module, and the spatia discretization of the flow field in the greenhouse was carried out using a cut-body gri (trim). At a spatial grid volume of 1850 w, the volume change quality of the grid cells wa 0.938, thus meeting the solution requirements.

Geometric Model of the Plant Factory
This experiment was carried out at Jilin University's artificial lig plant factory's specification is 5575 mm × 3950 mm × 3175 mm, as sh The three-dimensional model was imported into the star ccm+ mod discretization of the flow field in the greenhouse was carried out u (trim). At a spatial grid volume of 1850 w, the volume change quality 0.938, thus meeting the solution requirements.

Control Equation
The factory LED light source is the only heat source, and the k-ε turbulence model was used for momentum equation. The governing equations include the continuity equation, the momentum equation, the energy equation, the K turbulent kinetic energy equation, and the ε dissipation rate equation, which can be expressed by general Equation (3) ρ is the density, kg/m 3 ,

Simulation Boundary Conditions of the Plant Factory
According to the photosynthetic physiology experiment, T S and T G are the LED light temperatures under the optimal light treatment mode during the seedling and growth stages, respectively. The gradient term used the least square method; the pressure term used the standard algorithm; and the momentum, energy, and viscosity terms all used the first-order up airflow scheme to achieve faster convergence. Boundary conditions of CFD simulations are shown in Table 3, the relaxation factor for the energy term was set to 10 −6 , and the rest of the terms are set to 10 −3 . We set the inlet as the velocity inlet, with 8 m/s as the airflow speed. Meanwhile, we set the outlet as the pressure outlet, with 0 Pa as the pressure. Wall Insulation T S is lamp temperature for optimal light treatment at the seedling stage; T G is lamp temperature for optimal light treatment at growth stage.

Plant Physiological Indexes
The analytic hierarchy process (Ahp) method was used to obtain the comprehensive evaluation set of the photosynthetic physiology experiment. The original indices for seedling and growth stages in Table 4 were normalized and multiplied by the Wi weight coefficient. Because the tip burn rate and nitrate levels were both negative indicators of lettuce quality, the reciprocal calculation was used after normalization. The comprehensive evaluation subset of each treatment was obtained based on the growth and physiological indexes in Table 4 and the weight coefficients of each index. In the seedling stage, the comprehensive evaluation subset B1 = (0.5847, 0.917, 0.929, 0.621, 0.983, 0.979), while in the growth stage, the comprehensive evaluation subset B2 = (0.618, 0.7976, 0.8875, 0.6648, 0.9254, 1.077). According to the B1 and B2 evaluation sets, the light pattern of the seedling stage was ranked T4 > T5 > T2 > T1 > T3 > CK1, while the light pattern of the growing stage was ranked W5 > W4 > W2 > W1 > W3 > CK2. The best light mode, T S , T4 (2R:1B, 100 µmol/m 2 ·s), was selected during the seedling period; the best light mode, T G , W5 (4R:1B, 200 µmol/m 2 ·s), was selected during the growth stage. Consistent with the findings of Haaak et al., low light intensity, according to the experimental results, can significantly reduce the incidence of tip burn disease while inhibiting growth. Red and blue LEDs, with a ratio between 5 and 20% for the blue light, are optimal to improve lettuce plant growth and light use efficiency. However, red and blue light sources do not significantly alleviate tip burn disease when compared to full-spectrum LED lamps [42,43].

Airflow Experiment Results and Analysis
As shown in Figure 5, for the seedling stage lettuce, when the airflow speed is lower than 0.5 m/s, the tip burn rate decreases to 31 However, the airflow speed should be, at most, the critical value of 1.4 m/s; otherwise, the lettuce leaves will be broken, reducing plant yield and quality. In determining airflow thresholds unfavorable to the development of tip burn in lettuce, based on airflow treatment experiments, the seedling stage airflow thresholds were 0.1-0.5 m/s, and the growth stage airflow threshold was 1.1-1.4 m/s. This airflow experiment showed that appropriate airflow speed could reduce lettuce tip burn disease.

Porous Medium Parameters
The average pressure drop values at six monitoring points were obtained based on aerodynamic tests conducted in a low-speed airflow tunnel for seedlings and growing lettuce. The pressure drop curves (∇P−v) were plotted, as shown in Figure 6. Based on th pressure drop curves, we deduced the inertial and viscous drag values in three vecto directions (x, y, and z) for the lettuce seedling and growth stage, as shown in Table 5.

Porous Medium Parameters
The average pressure drop values at six monitoring points were obtained based on aerodynamic tests conducted in a low-speed airflow tunnel for seedlings and growing lettuce. The pressure drop curves (∇P−v) were plotted, as shown in Figure 6. Based on the pressure drop curves, we deduced the inertial and viscous drag values in three vector directions (x, y, and z) for the lettuce seedling and growth stage, as shown in Table 5.
The average pressure drop values at six monitoring points were obtained based on aerodynamic tests conducted in a low-speed airflow tunnel for seedlings and growing lettuce. The pressure drop curves (∇P−v) were plotted, as shown in Figure 6. Based on the pressure drop curves, we deduced the inertial and viscous drag values in three vector directions (x, y, and z) for the lettuce seedling and growth stage, as shown in Table 5.

Optimization of Three-Dimensional Cultivation
In the airflow field of the plant factory with no load, as shown in Figure 7, the high airflow speed area is concentrated in the H (high airflow speed) area of the cultivation rack, while the low speed area is concentrated in the L (low airflow speed) area. According to the simulation results of the airflow field in the no-load plant factory, lettuce was

Optimization of Three-Dimensional Cultivation
In the airflow field of the plant factory with no load, as shown in Figure 7, the high airflow speed area is concentrated in the H (high airflow speed) area of the cultivation rack, while the low speed area is concentrated in the L (low airflow speed) area. According to the simulation results of the airflow field in the no-load plant factory, lettuce was planted in the L (low airflow speed) area during the seedling stage and the H (high airflow speed) area during the growth stage. Figure 8 depicts a planting density distribution that ensured that the number of seedlings and growth stage plants correspond to 1:1, thus facilitating the rapid cycling cultivation of lettuce under controlled environments.     . Cultivation pattern diagram. The seedling zone (low airflow speed zone) planting density is twice as high as the growing period (high airflow speed zone) planting density, and the seedling zone (low airflow speed zone) planting area is 0.5 times the growing period (high airflow speed zone) planting area.

Optimization of Airflow Field
In the H (high airflow speed) area, a porous media model for the growth stage was input and the LED lights were arranged in W5 mode (4R:1B, 200 μmol/m 2 .s), with a heat source temperature of 307.15 K. In the L (low airflow speed) area, a porous media model for the lettuce seedling stage was input, the LEDs were arranged in the T4 mode (2R:1B, 100μmol/m 2 .s), and the heat source input temperature was 304.15K. The rest of the boundary conditions were the same as for the no-load CFD simulation. Therefore, creating vortex areas in the back area where the air supply was opposite the heat source was simple. The dominant indoor airflow after being blocked and attenuated by the supply heat source is the updraft caused by the heat source. The axial speed of the z direction changes very little in the vertical direction of the upper part of the heat source. The updraft changed direction and was discharged from the return air outlet only after it reached the top and was subjected to the combined action of the ceiling and the return air outlet. To address the problem of the uneven distribution of airflow field zoning in the unloaded plant, a fan was installed on the left side of each cultivation frame, with the fan position Figure 8. Cultivation pattern diagram. The seedling zone (low airflow speed zone) planting density is twice as high as the growing period (high airflow speed zone) planting density, and the seedling zone (low airflow speed zone) planting area is 0.5 times the growing period (high airflow speed zone) planting area.

Optimization of Airflow Field
In the H (high airflow speed) area, a porous media model for the growth stage was input and the LED lights were arranged in W5 mode (4R:1B, 200 µmol/m 2 s), with a heat source temperature of 307.15 K. In the L (low airflow speed) area, a porous media model for the lettuce seedling stage was input, the LEDs were arranged in the T4 mode (2R:1B, 100 µmol/m 2 s), and the heat source input temperature was 304.15 K. The rest of the boundary conditions were the same as for the no-load CFD simulation. Therefore, creating vortex areas in the back area where the air supply was opposite the heat source was simple. The dominant indoor airflow after being blocked and attenuated by the supply heat source is the updraft caused by the heat source. The axial speed of the z direction changes very little in the vertical direction of the upper part of the heat source. The updraft changed direction and was discharged from the return air outlet only after it reached the top and was subjected to the combined action of the ceiling and the return air outlet. To address the problem of the uneven distribution of airflow field zoning in the unloaded plant, a fan was installed on the left side of each cultivation frame, with the fan position shown in Figure 9 and the cultivation frames numbered as shown in Figure 10. The differential speed of the fans was adjusted according to the flow field simulation, with the specific fan speed settings shown in Table 6. shown in Figure 9 and the cultivation frames numbered as shown in Figure 10. The differential speed of the fans was adjusted according to the flow field simulation, with the specific fan speed settings shown in Table 6.     Following the differential speed adjustment of the fan at the side end of each cultivation frame, the flow field of each cultivation frame plant canopy is shown in Figure 11. The airflow comfort zone values for each cultivation frame are shown in Table 7. The average growth comfort zone as a percentage of the total low airflow speed cultivation frame at the seedling stage (low airflow speed zone) was 84.5%. At the growth stage (high airflow speed zone), the average growth comfort zone as a percentage of the total high airflow speed cultivation frame was 94.1%.  Following the differential speed adjustment of the fan at the side end of each cultivation frame, the flow field of each cultivation frame plant canopy is shown in Figure 11. The airflow comfort zone values for each cultivation frame are shown in Table 7. The average growth comfort zone as a percentage of the total low airflow speed cultivation frame at the seedling stage (low airflow speed zone) was 84.5%. At the growth stage (high airflow speed zone), the average growth comfort zone as a percentage of the total high airflow speed cultivation frame was 94.1%.  S is seedling stage lettuce; G is growth stage lettuce. Figure 11. Vector diagram of cultivation frame canopy airflow field. The dotted line represents the low airflow speed zone; the rest of the area is the high airflow speed zone.

Verification of Simulation Results
According to Figure 12, 33 monitoring points were evenly spaced at 20 cm above the cultivation frame and measured using an RA620 anemometer. The average absolute error was 6.342%, and the simulation model was accurate and reliable.

Verification of Simulation Results
According to Figure 12, 33 monitoring points were evenly spaced at 20 cm above the cultivation frame and measured using an RA620 anemometer. The average absolute error was 6.342%, and the simulation model was accurate and reliable. Figure 11. Vector diagram of cultivation frame canopy airflow field. The dotted line represents the low airflow speed zone; the rest of the area is the high airflow speed zone.

Verification of Simulation Results
According to Figure 12, 33 monitoring points were evenly spaced at 20 cm above the cultivation frame and measured using an RA620 anemometer. The average absolute erro was 6.342%, and the simulation model was accurate and reliable.

Comparison of Tip Burn Rate of Plants before and after Optimization
Lettuce was cultivated according to the planting pattern after the HL (high airflow area and low airflow area) division. As shown in Table 8, the tip burn rate in the seedling period decreased to 26.9%, and that in the growth period, it decreased to 23.3%. Compared with the tip burn rate before optimizing the airflow field, the tip burn rate decreased by 64.3% in the seedling stage and by 69.1% in the growth stage, and the growth quality was significantly improved.

Comparison of Tip Burn Rate of Plants before and after Optimization
Lettuce was cultivated according to the planting pattern after the HL (high airflow area and low airflow area) division. As shown in Table 8, the tip burn rate in the seedling period decreased to 26.9%, and that in the growth period, it decreased to 23.3%. Compared with the tip burn rate before optimizing the airflow field, the tip burn rate decreased by 64.3% in the seedling stage and by 69.1% in the growth stage, and the growth quality was significantly improved. The analytic hierarchy process (Ahp) weight coefficient was used to evaluate the comprehensive quality of lettuce after optimization, and the comprehensive evaluation set of seedling stage lettuce b1 was (0.983, 1.202), and the comprehensive evaluation set of growth stage lettuce b2 was (1.077, 1.354). According to the comprehensive evaluation value, the quality of lettuce in the seedling stage and growth stage improved after optimization. Comparing the lettuce growth quality and tip burn before and after zoning, the overall growth quality at the seedling stage increased by 22.28%, and the integrated quality at the growth stage increased by 25.72%. Therefore, the optimization and improvement described in this paper could improve the problem of tip burn in plant factories.

Discussion
The photosynthesis experiments showed that the best comprehensive quality of seedling lettuce was achieved under T4 lighting treatment (2R:1B, 100 µmol/m 2 s), and the best comprehensive quality of growing lettuce was achieved under W5 (4R:1B, 200 µmol/m 2 s) lighting treatment. Although photosynthesis is integral to tip burn disease, this issue has yet to be well addressed. Meanwhile, low light intensity reduces tip burn disease and significantly impacts plant growth quality; this result is consistent with the findings of Carotti et al. [44]. Furthermore, the increase in lamp temperature, as the heat source of LED, will cause the canopy airflow to move upward. However, we expect the airflow to move into the plant canopy [45]. Therefore, appropriate LED lighting can reduce heat dissipation disturbances to the plant canopy's airflow field, save energy, and improve plant quality and yield.
The airflow experiment showed that the lettuce tip burn rate decreased when the seedling stage's airflow speed was 0.1-0.5 m/s and the airflow speed in the growth stage was 1.1-1.4 m/s. In contrast to other studies, such as that of Luuk et al., which merely increased the airflow speed to reduce the resistance of the plant boundary layer, we found that the plant canopy at the seedling stage required a lower airflow speed to maintain the water potential of the leaves [46]. In addition, experiments showed that, in the growth stage, vertical airflow from top to bottom was more effective than horizontal blowing in treating tip burn disease. This could be because vertical airflow could increase airflow in the plant's canopy and relieve tip burn disease in the first leaf position. In contrast to the demand trend of airflow speed in the seedling stage and growth stage, this may be because too high an airflow speed will accelerate the leaf water potential decrease in the seedling stage, resulting in tip burn; increasing airflow speed in the growth stage can promote high humidity airflow around plants, promote plant transpiration power, and subsequently, promote plant Ca 2+ transport.
The aerodynamic tests in the three xyz vector directions for lettuce in the seedling and growth stage, respectively, as well as the porous media model of lettuce in different periods, were conducted. This study showed that the viscous drag coefficient D and inertial drag coefficient C 2 of porous media in lettuce at the seedling and growth stage were highly variable; thus, it was necessary to test them separately. The significant differences in porous media parameters compared to those of tomatoes and other vegetables referenced in previous studies occur because lettuce has different physical characteristics from other plants and because of the high density of lettuce cultivation in plant factories [45].
This study cultivated lettuce in the seedling stage and the growth stage in different regions. We optimized the three-dimensional cultivation structure of the plant factory based on the flow field distribution of an empty plant factory; one side's airflow speed was high, and the other side's airflow speed was low. In the upper high and low airflow speed areas, we set the LED illumination to 2R:1B (100 µmol/m 2 s) and 4R:1B (200 µmol/m 2 s), respectively. Compared with the traditional plant factory cultivation mode, our light source arrangement pattern reduces light source heat dissipation and plant factories operating costs. To further optimize the uniformity of the airflow between the various cultivation racks, the side fans of the cultivation frame are adjusted for differential speed, according to the simulation results. The airflow field can be adjusted according to the required environment of the cultivated crop, due to the versatility of the fan unit.

Conclusions
Intending to solve the problem of tip burn caused by non-uniform airflow distribution in existing plant factories, this study proposed a method for zoning plants based on photosynthetic physiology experiments and airflow field distribution in unladen plant factories, and further optimizing the airflow field by arranging differential fans in cultivation frames. The best comprehensive quality of seedling lettuce was achieved when the airflow speed was 0.1-0.5 m/s in T4 light mode (2R:1B, 100 µmol/m 2 s), and the best comprehensive quality of growing lettuce was achieved when the airflow speed was 1.1-1.4 m/s in W5 light mode (4R:1B, 200 µmol/m 2 s). The zoning optimization cultivation method and differential fan arrangement were used in this study to improve the airflow field of plant factories. After setting the differential fan to regulate the airflow field, the percentage of the suitable growth area reached 84.5% in the seedling stage and 94.1% in the growth stage. Compared with the lettuce growth quality and tip burn rate before and after zoning, the integrated quality at the seedling stage increased by 22.28% and the tip burn decreased to 26.9%; the integrated quality at the growth stage increased by 25.72%, and the tip burn decreased to 23.2%. Based on the analytic hierarchy process (Ahp) comprehensive evaluation values, it can be seen that the comprehensive quality of both the seedling and growth stages were improved after optimized cultivation. In this study, the optimization of plant cultivation methods and airflow distribution can effectively reduce the tip burn rate in the plant factory, which has a good guiding significance for the sustainable development of plant factories. Future research can quantify the energy costs and crop benefit output after zone-optimized cultivation to provide more favorable support for plant factory sustainability studies.

Conflicts of Interest:
The authors declare no conflict of interest.