A Standardized Treatment Model for Head Loss of Farmland Filters Based on Interaction Factors

: A head loss model for pressureless mesh filters used in farmland irrigation was developed by integrating the four basic test factors: irrigation flow, filter cartridge speed, self-cleaning flow, and initial sand content. The model’s coefficient of determination was found to be 98.61%. Among the basic factors, the total irrigation flow accounted for only 17.20% of the relatively small self-cleaning flow. The contribution of initial sand content was found to be the smallest, with a coefficient of only 0.0166. Furthermore, the contribution rate of the flow term was significantly higher than that of the initial sand content, with a value of 159.73%. In terms of quadratic interaction, the difference between the interaction term of flushing flow and filter cartridge speed, and the interaction term of filter cartridge speed and self-cleaning flow was 38.42%. On the other hand, the difference within this level for the interaction term between initial sand content and filter cartridge speed, as well as the interaction term between irrigation flow and self-cleaning flow, was 2.82%. Finally, through joint optimization of the response surface and model, the optimal values for the irrigation flow rate, filter cartridge speed, self-cleaning flow rate, and initial sand content were determined to be 121.687 m 3 · h − 1 , 1.331 r · min − 1 , 19.980 m 3 · h − 1 , and 0.261 g · L − 1 ; the measured minimum head loss was found to be 21.671 kPa. These research findings can serve as a reference for enhancing the design of farmland filters and optimizing irrigation systems.


Introduction
Water scarcity poses a significant challenge to agricultural productivity in arid regions, necessitating the development of efficient irrigation systems.Among the vital components of these systems are irrigation filters-specifically mesh filters or screen filters-which facilitate the removal of debris and sediments from the water supply.The hydraulic performance of these filters directly impacts their ability to maintain a stable water flow and prevent clogging issues, ensuring the effective preservation of irrigation water by removing impurities [1,2].In Xinjiang, where irrigation water-source impurity levels are high, filters are essential to prevent serious internal blockages caused by impurities entering the pump [3].
In recent years, there has been a growing body of research focused on studying the hydraulic performance of filters in the context of arid agricultural irrigation.These studies have explored various aspects related to the filtration process, including flow characteristics, wear patterns, head loss, and pollutant-removal capabilities.By investigating these factors, researchers aim to enhance their efficiency, further enriching the knowledge surrounding the preservation of irrigation water and the role of filtration in arid regions.Some studies conducted [4,5] utilized numerical simulations to investigate the flow and wear characteristics of a Y-type mesh filter.The findings revealed a direct linkage between inlet velocity, vorticity, and the hydraulic performance of filters used in agricultural irrigation systems in arid regions.They provide insights into the relationship between flow characteristics and filter performance, contributing to the optimization of filter design and operation.The head loss of well screens was quantitatively studied [6,7] through the coupling of Darcy and turbulent flow numerical models.Their research yielded an empirical formula for quantifying head loss, contributing to our understanding of filter performance in arid agricultural irrigation.The study showed valuable insights into the flow characteristics at sub-millimeter scales and serves as a foundation for optimizing filter design and operation.Physical model tests on micro-pressure filters were conducted [8,9], and they identified the flow rate as the primary factor affecting head loss, followed by the water separator type.Their research provides crucial insights into the impact of these factors on the hydraulic performance of filters used in arid agricultural irrigation.By highlighting the importance of the flow rate and water-separator type, this study aids in optimizing filter design and operation.Some researchers studied [10,11] the hydraulic performance and flow behavior of an automatic flushing filter with the aim of reducing filtration pressure drop and expanding the operating flow range.Their research contributes to an improved understanding of filter shell design optimization and its relationship to water filtration in agricultural irrigation systems.This study directly contributes to the preservation of irrigation water quality and the prevention of clogging issues, both of which are crucial considerations in arid farming regions.Optimal discharge times for self-cleaning filter screens in micro-irrigation systems were explored [12].Through their experiments, they analyzed local head loss variation laws under clean and muddy water conditions, enabling a deeper understanding of the filtration process and its performance in arid agricultural irrigation.This research provides valuable insights into the factors influencing the hydraulic performance of irrigation filters and contributes to the optimization of filter operation in water-scarce environments.Some researchers studied [13,14] the hydraulic efficiency and pollutant-removal capabilities of a horizontal flow multimedia rainwater filter.Through column and batch studies, they developed corresponding isotherms for coconut fiber adsorption behavior, contributing to the understanding of the filtration process in arid agricultural irrigation systems.This research facilitates the adoption of effective mechanisms for pollutant removal and underscores the importance of hydraulic efficiency in preserving irrigation water quality.The hydraulic performance of a bottom flow screen was examined [15,16], evaluating the effect of geometric variables on capture flow and determining collection capacity and sediment removal efficiency.Their physical and numerical experiments shed light on the impact of filter geometry on hydraulic performance, specifically in capturing sediment and maintaining efficient water flow.This research serves as a valuable resource for optimizing filter design and operation in arid agricultural irrigation systems.
Although previous research has made strides in improving the hydraulic performance of filters through various optimizations, such as flow characteristics and wear patterns, there is a lack of studies that investigate the interaction between these influencing factors during filter operation, as well as the degree to which these factors impact the filter.Hence, this article focuses on a pressureless mesh filter as the subject of study and examines the effects of varying interactive factors on its hydraulic performance by considering factors such as filter flow rate, rotational speed of the cleaning cartridge, cleaning flow rate of the self-cleaning device, and initial sand addition rate.Through a systematic analysis, this study aims to provide a theoretical foundation and technical guidance for minimizing energy loss in pressureless mesh filters during operation, thereby promoting sustainable water use and supporting crop productivity in agricultural irrigation systems [17][18][19].

Structure of Pressure-Free Net Filter
The pressure-free filter is composed of a pontoon, and the pressure-free filter is composed of a float, a rotating motor, a dirt baffle, a filter cylinder, a main bracket, and a self-cleaning device.The overall structure and parameters (unit: mm) are shown in Figure 1.Silt-containing water for irrigation is deposited in the sedimentation tank and enters the filter from outside to inside due to suction from the water pump.The water is filtered by the filter screen, with impurities larger than the pore diameter forming accumulation.Meanwhile, water filtered by the Y-shaped filter on the component enters the filter through the spray pipe, and silt or organic impurities on the filter element are carried out by high-pressure spray holes.

Structure of Pressure-Free Net Filter
The pressure-free filter is composed of a pontoon, and the pressure-free filter is composed of a float, a rotating motor, a dirt baffle, a filter cylinder, a main bracket, and a selfcleaning device.The overall structure and parameters (unit: mm) are shown in Figure 1.Silt-containing water for irrigation is deposited in the sedimentation tank and enters the filter from outside to inside due to suction from the water pump.The water is filtered by the filter screen, with impurities larger than the pore diameter forming accumulation.Meanwhile, water filtered by the Y-shaped filter on the component enters the filter through the spray pipe, and silt or organic impurities on the filter element are carried out by high-pressure spray holes.

Test Device Structure
The test machine's entire structure primarily consists of six parts: water return device, electromagnetic flowmeter, electromagnetic pump, stirring device, filtration system, and high-precision negative pressure gauge.The equipment utilized in this test forms a comprehensive filtration circulation system.The diagram of the test system device can be found in Figure 2.

Test Device Structure
The test machine's entire structure primarily consists of six parts: water return device, electromagnetic flowmeter, electromagnetic pump, stirring device, filtration system, and high-precision negative pressure gauge.The equipment utilized in this test forms a comprehensive filtration circulation system.The diagram of the test system device can be found in Figure 2.

Test Flow and Materials
In accordance with the actual situation, a 48-mesh filter screen was utilized for the research.The prototype muddy water test was employed, with the filter flow rate, rotation speed of the self-cleaning net cylinder, self-cleaning device flow rate, and initial sand amount as the factors considered.The use of a specific filter screen allows for a consistent factor in the test setup.The prototype muddy water test was conducted to simulate real-world conditions and evaluate the performance of the pressure-free mesh filter.A fourfactor, three-level division of test groups was designed based on the pressure-free filter's operating conditions.This division allows for a comprehensive evaluation of the filter's performance under different operating conditions.The test flow was selected according to the previous tests, indicating that the motor-driven filter cartridge speed of 4 r/min yielded the best pressure-free filter performance.Consequently, the self-cleaning filter cartridge speed was set at (1, 2.5, and 4 r•min −1 ) to examine the impact of different rotation speeds on the filter's efficiency.

Test Flow and Materials
In accordance with the actual situation, a 48-mesh filter screen was utilized for the research.The prototype muddy water test was employed, with the filter flow rate, rotation speed of the self-cleaning net cylinder, self-cleaning device flow rate, and initial sand amount as the factors considered.The use of a specific filter screen allows for a consistent factor in the test setup.The prototype muddy water test was conducted to simulate realworld conditions and evaluate the performance of the pressure-free mesh filter.A fourfactor, three-level division of test groups was designed based on the pressure-free filter's operating conditions.This division allows for a comprehensive evaluation of the filter's performance under different operating conditions.The test flow was selected according to the previous tests, indicating that the motor-driven filter cartridge speed of 4 r/min yielded the best pressure-free filter performance.Consequently, the self-cleaning filter cartridge speed was set at (1, 2.5, and 4 r•min −1 ) to examine the impact of different rotation speeds on the filter's efficiency.
The pipeline flow calculation was based on the micro-irrigation engineering standard: the pipeline operation scenario under the most unfavorable hydraulic conditions was adopted, and the calculation accuracy of all flow variations in emi ers was not to exceed 20%.Therefore, according to the actual situation, the flow factor was determined to set the test flow, resulting in three intermediate values (120, 140, and 160 m 3 •h −1 ).Considering the principle that sediment concentration should not be too high to affect the test sample and not too low to prolong the experiment, the initial sediment concentration factor was set as (0.2, 0.5, and 0.8 kg•m −3 ).This range allows for a variation in sediment concentration without compromising the representativeness of the experimental samples.Additionally, the total flow extracted by the self-cleaning device was determined by combining the test parameter combinations with the calculation of the diffusion section of the self-cleaning nozzle.Three flows (1, 5, and 20 m 3 •h −1 ) were selected to evaluate the efficiency of the selfcleaning mechanism under different flow rates.
During the experiments, the filter system was equipped with a pressure gauge and a flow meter to monitor the performance of the pressure-free mesh filter.To ensure accurate readings, the pressure gauge and flow meter readings were taken three times after stabilizing.
In this study, natural bed sand from Site No. 146 of the Xinjiang Corps was utilized as the sediment material.The sediment particle grading curve of the sand is presented in The pipeline flow calculation was based on the micro-irrigation engineering standard: the pipeline operation scenario under the most unfavorable hydraulic conditions was adopted, and the calculation accuracy of all flow variations in emitters was not to exceed 20%.Therefore, according to the actual situation, the flow factor was determined to set the test flow, resulting in three intermediate values (120, 140, and 160 m 3 •h −1 ).Considering the principle that sediment concentration should not be too high to affect the test sample and not too low to prolong the experiment, the initial sediment concentration factor was set as (0.2, 0.5, and 0.8 kg•m −3 ).This range allows for a variation in sediment concentration without compromising the representativeness of the experimental samples.Additionally, the total flow extracted by the self-cleaning device was determined by combining the test parameter combinations with the calculation of the diffusion section of the self-cleaning nozzle.Three flows (1, 5, and 20 m 3 •h −1 ) were selected to evaluate the efficiency of the self-cleaning mechanism under different flow rates.
During the experiments, the filter system was equipped with a pressure gauge and a flow meter to monitor the performance of the pressure-free mesh filter.To ensure accurate readings, the pressure gauge and flow meter readings were taken three times after stabilizing.
In this study, natural bed sand from Site No. 146 of the Xinjiang Corps was utilized as the sediment material.The sediment particle grading curve of the sand is presented in Figure 3.This information provides insights into the characteristics of the sediment used in the experimental tests.

Analysis Methods
To comprehensively investigate the interactive effects among the parameters of the pressure-free net filter device, a four-factor, three-level response surface methodology (RSM) was employed.The four critical factors under consideration were filter flow rate, self-cleaning net drum rotation speed, self-cleaning device flow rate, and initial sand addition rate.

Analysis Methods
To comprehensively investigate the interactive effects among the parameters of the pressure-free net filter device, a four-factor, three-level response surface methodology (RSM) was employed.The four critical factors under consideration were filter flow rate, self-cleaning net drum rotation speed, self-cleaning device flow rate, and initial sand addition rate.
(1) Establishment of the response surface function: The response surface function, denoted as (Iij), was formulated to capture the relationships between the factors and the response variable.This function encompassed constant, linear, and interaction terms.
where Iij is the construct response surface function, a0 is the constant term coefficient, aj is the first term coefficient, and aij is the second term coefficient.The author of the response equation is Nikolai Vasilievich Ermolaev.The response surface methodology was introduced by George E. P. Box and William G. Hunter in 1957.In the response surface methodology, the response surface function (Iij) is used to capture the relationship between factors and response variables, including constant terms, linear terms, and interaction terms.
(2) Construction matrix A design matrix was constructed based on the chosen levels of the factors.This matrix facilitated the organization and analysis of the experimental data.(1) Establishment of the response surface function: The response surface function, denoted as (I ij ), was formulated to capture the relationships between the factors and the response variable.This function encompassed constant, linear, and interaction terms.
where I ij is the construct response surface function, a 0 is the constant term coefficient, a j is the first term coefficient, and a ij is the second term coefficient.The author of the response equation is Nikolai Vasilievich Ermolaev.The response surface methodology was introduced by George E. P. Box and William G. Hunter in 1957.In the response surface methodology, the response surface function (I ij ) is used to capture the relationship between factors and response variables, including constant terms, linear terms, and interaction terms.
(2) Construction matrix A design matrix was constructed based on the chosen levels of the factors.This matrix facilitated the organization and analysis of the experimental data.

3) Determine structural coefficients
To determine the construction coefficient, β j , m groups of independent tests are conducted in this paper, with each test requiring different parameters for each variable.Using this approach, we can obtain m response values corresponding to the points.
Substitute the above m sample points, x(j) (j = 0, 1, . .., m − 1), into Equation ( 4) to obtain the response surface function value as follows: The response surface function (I ij ) served as an approximation of the actual performance function, introducing an error term.
To minimize and obtain the response surface closest to the test data point, the least squares method is used.
The necessary condition for obtaining the minimum value of Equation ( 8) is as follows: (4) Combining Equations ( 1)-( 9) yields the following: Next, it is converted to matrix form as The final response surface equation was obtained by combining the equations derived from the previous steps.This equation provided a mathematical representation of the relationships between the factors and the response variable, facilitating a further analysis and optimization, wherein we have the following:

Standardized Processing
For the purpose of standardization processing, z-score standardization is employed as the initial step.This technique ensures uniformity in the scale of the data and mitigates biases and misunderstandings resulting from different dimensions.Additionally, z-score normalization centers the data on zero, removing any inherent bias and enhancing data comprehensibility and comparability.Simultaneously, this standardization method preserves the distribution shape and relative positioning of the original data, thereby retaining the information contained within.Furthermore, the utilization of z-score standardized data facilitates the process of data analysis and modeling, as the consistent data scale and eliminated biases enhance the accuracy and reliability of the model.Consequently, z-score normalization is employed to address the issue of data comparability in our study.
where µ is the mean (average) of the raw data, σ is the standard deviation of the raw data, and z is the standardized value.

Results and Discussion
According to the test parameter combination and the Box-Behnken [20][21][22] model center principle, a four-factor, three-stage response surface test was conducted using the RSM theoretical design, as illustrated in the response surface basic factor table (Table 1).

Head Loss Model Calculation and Variance Analysis
A multi-factor analysis of variance (ANOVA) was conducted to analyze the standardized test results.The ANOVA results are presented in Table 3.The number '2' in the "Divisor" column represents the quadratic term.
The standardized model is presented in Table 4.The number '2' in the "Divisor" column represents the quadratic term.
The index function is influenced by numerous factors, and their combined impact on the target is not independent.Instead, there is a common influence on the dependent variable.First, employing multi-factor variance analysis [23][24][25], the joint effect of multiple factors on dependent variables is analyzed.Quadratic interaction terms with various parameters and factors are subject to variance analysis, resulting in a table with 14 degrees of freedom.
The coefficient of determination (R 2 ) for the model is 0.9861, indicating that the model explains a significant portion of the variance in the response variable.The p-values obtained from the ANOVA analysis were used to test the significance of each factor and their interactions.A p-value less than the significance level (α = 0.05) indicates that the factor has a statistically significant effect on the response variable.
In this study, the head loss model was analyzed in terms of the initial influencing factors: flow rate, rotation speed, self-cleaning flow rate, initial sediment concentration, their respective quadratic terms, and the interaction of the four factors.The results showed that the p-values of factors such as flow rate, rotation speed, and self-cleaning flow rate were significant, indicating their crucial role in determining head loss.However, the p-value of the initial sediment concentration factor was found to be non-significant, suggesting a limited direct influence on head loss.
The model's goodness-of-fit was further evaluated by comparing the R 2 value with the adjusted R 2 and predicted R 2 values.A close agreement between these values indicates the model's reliability and predictive capability.Additionally, the lack of fit test was performed to assess whether the model adequately represents the data.A non-significant lack of fit (p-value > 0.05) indicates that the model is suitable for predicting the response variable.
In conclusion, this study successfully applied the Box-Behnken design and response surface methodology to analyze the head loss model in terms of various influencing factors.The results provide valuable insights into the factors that significantly impact head loss and can be used to optimize system performance.Future research could explore additional factors or refine the existing model to further enhance its predictive accuracy.

Factor Contribution to Head Loss
Table 5 presents the calculation of the contribution of each factor to the model.The number '2' in the "Divisor" column represents the quadratic term.
The RSM factor's contribution to head loss reflects the ratio of each factor's useful results to resource consumption and occupation, demonstrating the extent of its impact on index head loss.
We analyzed the factors influencing head loss, including flow rate, rotation speed, selfcleaning flow rate, initial sediment concentration, their quadratic terms, and the pairwise joint calculation contribution degree of these four influencing factors.The reference index is derived from the absolute value estimated by the system [26].
Table 5 reveals the contribution rate order for head-loss influencing factors: flow rate term, self-cleaning flow rate term, self-cleaning flow 2 term, rotation speed term, flow 2 term, speed × self-cleaning flow term, speed × initial sediment concentration term, flow × self-cleaning flow term, initial sediment concentration 2 term, self-cleaning flow × initial sediment concentration term, initial sediment concentration term, and flow × initial sediment concentration term.The estimation suggests that a better-fitting model corresponds to a larger value and improved explanation of the independent variable to the dependent variable.
Figure 4 presents the contribution parameter diagram of filter factors.In this diagram, Figure 4a depicts the estimated value of the model system and the VIF (Variance Inflation Factor) value of the factor, while Figure 4b illustrate the lower and upper limits of the confidence interval associated with the factor.ment concentration term.The estimation suggests that a be er-fi ing model corresponds to a larger value and improved explanation of the independent variable to the dependent variable.
Figure 4 presents the contribution parameter diagram of filter factors.In this diagram, Figure 4a depicts the estimated value of the model system and the VIF (Variance Inflation Factor) value of the factor, while Figure 4b illustrate the lower and upper limits of the confidence interval associated with the factor.Analyzing the data presented in Figure 4, it can be observed that, in the basic factor group, flow and self-cleaning flow exert dominant contributions, with absolute contribution values of 1.02 and 0.8445, respectively.However, the contribution of total irrigation flow to head loss is significantly smaller compared to flow.The contribution of self-cleaning flow is 17.20% higher, whereas the contribution of initial sand content is the smallest, merely 0.0166.Moreover, the contribution of self-cleaning flow surpasses the contribution of rotational speed by 46.18%, and the contribution of rotational speed exceeds the contribution of initial sand content by 96.34%.Notably, the largest disparity between groups lies in the contribution of the flow term, which is 159.73% higher than the contribution of the initial sand content.
In the quadratic interaction term factor group, the interaction term of flow rate and rotational speed prevails with an absolute value of 0.2080.Conversely, the contribution of the interaction term between flow rate and initial sand content is the smallest, merely 0.036.This quadratic interaction term factor group can be divided into three levels based on their values.The first level consists of the interaction term between flow rate and rotational speed, as well as the interaction term between rotational speed and self-cleaning flow, exhibiting a 38.42% difference within the level.The second level includes the interaction term between initial sand content and rotational speed, as well as the interaction term between initial sand content and self-cleaning flow, with a 2.82% difference within the level.Lastly, the third level comprises the interaction term between flow rate and initial sand content, as well as the interaction term between self-cleaning flow and initial Analyzing the data presented in Figure 4, it can be observed that, in the basic factor group, flow and self-cleaning flow exert dominant contributions, with absolute contribution values of 1.02 and 0.8445, respectively.However, the contribution of total irrigation flow to head loss is significantly smaller compared to flow.The contribution of self-cleaning flow is 17.20% higher, whereas the contribution of initial sand content is the smallest, merely 0.0166.Moreover, the contribution of self-cleaning flow surpasses the contribution of rotational speed by 46.18%, and the contribution of rotational speed exceeds the contribution of initial sand content by 96.34%.Notably, the largest disparity between groups lies in the contribution of the flow term, which is 159.73% higher than the contribution of the initial sand content.
In the quadratic interaction term factor group, the interaction term of flow rate and rotational speed prevails with an absolute value of 0.2080.Conversely, the contribution of the interaction term between flow rate and initial sand content is the smallest, merely 0.036.This quadratic interaction term factor group can be divided into three levels based on their values.The first level consists of the interaction term between flow rate and rotational speed, as well as the interaction term between rotational speed and self-cleaning flow, exhibiting a 38.42% difference within the level.The second level includes the interaction term between initial sand content and rotational speed, as well as the interaction term between initial sand content and self-cleaning flow, with a 2.82% difference within the level.Lastly, the third level comprises the interaction term between flow rate and initial sand content, as well as the interaction term between self-cleaning flow and initial sand content, with an 87.72% difference within the level.Furthermore, the largest disparity between groups is observed in the contribution of the interaction term between flow rate and rotational speed, which is 218.25% higher than the contribution of the interaction term between flow rate and initial sand content.
In the quadratic self-interaction factor group, the contribution of the self-cleaning flow interaction dominates, with an absolute contribution value of 0.6896.Conversely, the contribution of the self-interaction term of the rotational speed is the smallest, only 0.0289.Notably, the largest disparity between the groups lies in the contribution of the self-cleaning flow self-interaction, which is 180.24% higher than the contribution of the self-interaction term of the speed term.

Response Surface Analysis
Except for the quadratic interaction term with a VIF value greater than 1, a response effect analysis was conducted on the index head loss parameter using other interaction terms.

Interaction Term of Filter Flow Rate and Filter Cartridge Speed
(1) The response surface for the interaction factor filter flow rate and filter cartridge speed under the influence of a significant item is established in Figure 5.The response surface analysis reveals that, with a constant filter cartridge speed, the head loss increases as flow rate increases.Similarly, with a constant flow rate, the negative pressure value and head loss relative to speed increase and decrease.

Filter Flow and Self-Cleaning Flow Interaction
(1) This interaction item represents the interaction of filter flow and self-cleaning flow.The response surface, as illustrated in Figure 6, demonstrates that flow rate has a direct impact on negative pressure values.When the self-cleaning device's flow rate remains constant, increasing flow rate leads to decreased negative pressure values and increased head loss.Conversely, when flow rate is constant, increasing self-cleaning flow rate results in increased negative pressure values and decreased head loss.
(2) As illustrated in Figure 6, when the flow rate is from 130 m 3 /h to 160 m 3 /h, the curvature changes more sharply.The upper response surface exhibits a wider range of high negative pressure values, with rapid response value changes and larger curvature.The optimal depth curvature follows, followed by a region of intense change behind the violent change area.The region in front is more inclined towards a mixed state, indicating that irrigation flow and self-cleaning flow are relatively smaller under low flow rates, but more prominent under high flow rates.Consequently, head-loss changes are more pronounced.
(3) Figure 6 reveals that the flow rate positively affects head loss, while self-cleaning flow rate has a suppressing effect.In comparison to the violent change area, both selfcleaning and flow rates are higher, leading to a more mixed area and intense head-loss changes.This should be avoided.The figure demonstrates that this significant term has a negative interaction effect on head loss and plays a significant role in irrigation equipment from an energy perspective.The self-cleaning flow term is the primary influence factor, with the steepest curvature corresponding to a large effect.(2) Up to a flow rate of 150 m 3 /h, the response surface's upper part exhibits a large high-pressure region with a small curvature and broad area.Optimal curvature results in a significant decrease in the negative pressure value and a sharp increase in head loss.This indicates that an increased irrigation flow rate leads to a more turbulent flow pattern within the filter, thus increasing head loss.However, the contact surface between the filter cartridge and water flow is large, leaving little air above to create compaction effects on the water flow layer.At low rotation speeds, the filter screen clearance exerts a small force on water flow, allowing for easy passage through the filter screen device in the impurity layer.This reduces the collision times between sediment particles and the filter cartridge, resulting in a small head loss.
(3) From the response surface analysis to negative pressure in Figure 5, it is observed that this interaction item exhibits a positive interaction effect on head loss.The energy term has a minor influence on irrigation devices, with the flow term being the primary effect factor for this interaction item.The steepest curvature is small, suggesting that this interaction should be slightly reduced under irrigation-condition satisfaction.

Filter Flow and Self-Cleaning Flow Interaction
(1) This interaction item represents the interaction of filter flow and self-cleaning flow.The response surface, as illustrated in Figure 6, demonstrates that flow rate has a direct impact on negative pressure values.When the self-cleaning device's flow rate remains constant, increasing flow rate leads to decreased negative pressure values and increased head loss.Conversely, when flow rate is constant, increasing self-cleaning flow rate results in increased negative pressure values and decreased head loss.
(2) As illustrated in Figure 6, when the flow rate is from 130 m 3 /h to 160 m 3 /h, the curvature changes more sharply.The upper response surface exhibits a wider range of high negative pressure values, with rapid response value changes and larger curvature.The optimal depth curvature follows, followed by a region of intense change behind the violent change area.The region in front is more inclined towards a mixed state, indicating that irrigation flow and self-cleaning flow are relatively smaller under low flow rates, but more prominent under high flow rates.Consequently, head-loss changes are more pronounced.

Filter Flow and Initial Sediment Concentration Interaction
(1) This influential term represents the interaction of filter flow and the initial sediment concentration term.The response surface, as depicted in Figure 7, demonstrates that the flow rate has a direct impact on the negative pressure value when the initial sediment concentration remains constant.An increased flow rate leads to decreased negative pressure values and an increased head loss.In contrast, when the flow rate is constant, changes in the head loss relative to initial sediment concentration are minimal.
(2) Figure 7 reveals a small interaction effect of this interaction item on head loss.Compared to the intense change area, the influence of both high initial sediment concentration and discharge is particularly pronounced.This area tends to be more marginal, and head-loss changes are more significant.However, the interaction effect of this interaction item on head loss is not particularly evident.

Interaction Term of Cartridge Speed and Self-Cleaning Flow Rate
(1) For the interaction factors of filter cartridge speed and self-cleaning flow rate, its response surface is shown in Figure 8.The response surface analysis of this interaction item to negative pressure value in Figure 7 reveals that when the self-cleaning flow rate remains constant, the negative pressure value decreases with the increase in the cartridge (3) Figure 6 reveals that the flow rate positively affects head loss, while self-cleaning flow rate has a suppressing effect.In comparison to the violent change area, both selfcleaning and flow rates are higher, leading to a more mixed area and intense head-loss changes.This should be avoided.The figure demonstrates that this significant term has a negative interaction effect on head loss and plays a significant role in irrigation equipment from an energy perspective.The self-cleaning flow term is the primary influence factor, with the steepest curvature corresponding to a large effect.

Filter Flow and Initial Sediment Concentration Interaction
(1) This influential term represents the interaction of filter flow and the initial sediment concentration term.The response surface, as depicted in Figure 7, demonstrates that the flow rate has a direct impact on the negative pressure value when the initial sediment concentration remains constant.An increased flow rate leads to decreased negative pressure values and an increased head loss.In contrast, when the flow rate is constant, changes in the head loss relative to initial sediment concentration are minimal.

Filter Flow and Initial Sediment Concentration Interaction
(1) This influential term represents the interaction of filter flow and the initial sediment concentration term.The response surface, as depicted in Figure 7, demonstrates that the flow rate has a direct impact on the negative pressure value when the initial sediment concentration remains constant.An increased flow rate leads to decreased negative pressure values and an increased head loss.In contrast, when the flow rate is constant, changes in the head loss relative to initial sediment concentration are minimal.
(2) Figure 7 reveals a small interaction effect of this interaction item on head loss.Compared to the intense change area, the influence of both high initial sediment concentration and discharge is particularly pronounced.This area tends to be more marginal, and head-loss changes are more significant.However, the interaction effect of this interaction item on head loss is not particularly evident.(1) For the interaction factors of filter cartridge speed and self-cleaning flow rate, its response surface is shown in Figure 8.The response surface analysis of this interaction item to negative pressure value in Figure 7 reveals that when the self-cleaning flow rate remains constant, the negative pressure value decreases with the increase in the cartridge (2) Figure 7 reveals a small interaction effect of this interaction item on head loss.Compared to the intense change area, the influence of both high initial sediment concentration and discharge is particularly pronounced.This area tends to be more marginal, and head-loss changes are more significant.However, the interaction effect of this interaction item on head loss is not particularly evident.
3.3.4.Interaction Term of Cartridge Speed and Self-Cleaning Flow Rate (1) For the interaction factors of filter cartridge speed and self-cleaning flow rate, its response surface is shown in Figure 8.The response surface analysis of this interaction item to negative pressure value in Figure 7 reveals that when the self-cleaning flow rate remains constant, the negative pressure value decreases with the increase in the cartridge rotation speed term, and the head loss increases accordingly.Conversely, when the cartridge rotation speed term is constant, the negative pressure value increases with the increase in the self-cleaning flow rate, and the head loss decreases.
crease in the self-cleaning flow rate, and the head loss decreases.
(2) When the cartridge rotation speed ranges from 2 r/min and the self-cleaning flow rate is approximately 15 m 3 /h, 20 m 3 /h, and 28 m 3 /h, the curvature variation range indicates that the response value changes more sharply compared to Figures 5 and 7.However, the response value variation range is gentler compared to Figure 6.The curvature of the initial optimal path of the response surface graph changes rapidly, with a large area.As the optimal path curvature increases, it signifies a gradual decrease in negative pressure value.
This indicates that the amplitude range of head-loss changes more evenly, the mixing phenomenon is less prominent, and the annular state appears.For the initial sediment concentration and self-cleaning flow rate, the response value change is small, while for the initial sediment concentration and self-cleaning flow rate, the response value change is more intense.This suggests that the head loss decreases with the increase in self-cleaning device flow rate and increases with the increase in cartridge rotation speed.This is because the contact surface between the filter cartridge and water flow is not fully closed, and there is no complete obstruction to the water flow layer.The head loss depends on the collision times and entry probability of sediment particles with the filter cartridge.
(3) This significant term in Figure 7 has a negative interaction effect on head loss, with a large impact on irrigation equipment from an energy perspective.The self-cleaning flow term is the primary influence factor in this interaction item, with the steepest curvature.The interaction term has a negative effect on head loss.(1) A response surface was established to analyze the interactive effects of filter cartridge speed and initial sand content, as depicted in Figure 9. Additionally, another response surface was developed to examine the interactive effects of self-cleaning flow and initial sand content, as illustrated in Figure 10.According to the response surface analysis, it is observed that when the initial sediment concentration remains constant, the negative pressure value decreases with an increase in filter cartridge rotation speed, leading to an increase in head loss.Similarly, when the initial sediment is constant, the negative pressure value increases as the self-cleaning flow rate increases.However, when either the (2) When the cartridge rotation speed ranges from 2 r/min and the self-cleaning flow rate is approximately 15 m 3 /h, 20 m 3 /h, and 28 m 3 /h, the curvature variation range indicates that the response value changes more sharply compared to Figures 5  and 7.However, the response value variation range is gentler compared to Figure 6.The curvature of the initial optimal path of the response surface graph changes rapidly, with a large area.As the optimal path curvature increases, it signifies a gradual decrease in negative pressure value.
This indicates that the amplitude range of head-loss changes more evenly, the mixing phenomenon is less prominent, and the annular state appears.For the initial sediment concentration and self-cleaning flow rate, the response value change is small, while for the initial sediment concentration and self-cleaning flow rate, the response value change is more intense.This suggests that the head loss decreases with the increase in self-cleaning device flow rate and increases with the increase in cartridge rotation speed.This is because the contact surface between the filter cartridge and water flow is not fully closed, and there is no complete obstruction to the water flow layer.The head loss depends on the collision times and entry probability of sediment particles with the filter cartridge.
(3) This significant term in Figure 7 has a negative interaction effect on head loss, with a large impact on irrigation equipment from an energy perspective.The self-cleaning flow term is the primary influence factor in this interaction item, with the steepest curvature.The interaction term has a negative effect on head loss.

Interaction Term between Filter Cartridge Speed and Initial Sediment Concentration, and Interaction Term between Self-Cleaning Flow Rate and Initial Sediment Concentration
(1) A response surface was established to analyze the interactive effects of filter cartridge speed and initial sand content, as depicted in Figure 9. Additionally, another response surface was developed to examine the interactive effects of self-cleaning flow and initial sand content, as illustrated in Figure 10.According to the response surface analysis, it is observed that when the initial sediment concentration remains constant, the negative pressure value decreases with an increase in filter cartridge rotation speed, leading to an increase in head loss.Similarly, when the initial sediment is constant, the negative pressure value increases as the self-cleaning flow rate increases.However, when either the filter cartridge rotation speed or self-cleaning flow rate is held constant, the change in head loss relative to the initial sediment concentration is not significant.Figure 9

Interaction Correction
To account for interaction effects, a response surface analysis was employed to modify the regression equation compared to previous head loss analyses.The results are presented in Table 6.

Interaction Correction
To account for interaction effects, a response surface analysis was employed to modify the regression equation compared to previous head loss analyses.The results are presented in Table 6.(2) Figures 9 and 10 demonstrate that the corresponding interaction term has a minimal impact on head loss.In comparison to the intense change area, the curvature of the response surface tends to be higher in the edge state, indicating a more substantial change in head loss.

Interaction Correction
To account for interaction effects, a response surface analysis was employed to modify the regression equation compared to previous head loss analyses.The results are presented in Table 6.By integrating the response surfaces in Figures 5-10 and employing the optimized model, a total of 100 solutions were obtained by utilizing the method of determining the optimal slope of the response surface.The optimization plan is presented in Figure 11 as a data diagram.From this diagram, the best solution was identified by removing the standardization.It was found that, at a flow rate of 121.687 m 3 •h −1 , a rotation speed of 1.331 r•min −1 , a self-cleaning flow rate of 19.980 m 3 •h −1 , and an initial sand content of 0.261 g•L −1 , the minimum head loss obtained was 21.671 kPa.This result signifies the optimal combination of parameters that leads to the lowest head loss.The optimization plan outlined in Figure 11 serves as a reference for future research in this area.Researchers and practitioners can utilize this plan to guide their decision-making process when designing and operating filtration systems, with the aim of minimizing head loss and maximizing system performance.It should be noted that the methodology employed in this study, which involved combining response surfaces and an optimized model, demonstrates the effectiveness of such an approach in understanding complex relationships between multiple factors.By visualizing the data in Figures 5-10 and applying the method of finding the optimal slope of the response surface, a comprehensive set of solutions was obtained, leading to the identification of the best possible combination of parameters.By integrating the response surfaces in Figures 5-10 and employing the optimized model, a total of 100 solutions were obtained by utilizing the method of determining the optimal slope of the response surface.The optimization plan is presented in Figure 11 as a data diagram.From this diagram, the best solution was identified by removing the standardization.It was found that, at a flow rate of 121.687 m 3 •h −1 , a rotation speed of 1.331 r•min −1 , a self-cleaning flow rate of 19.980 m 3 •h −1 , and an initial sand content of 0.261 g•L −1 , the minimum head loss obtained was 21.671 kPa.This result signifies the optimal combination of parameters that leads to the lowest head loss.The optimization plan outlined in Figure 11 serves as a reference for future research in this area.Researchers and practitioners can utilize this plan to guide their decision-making process when designing and operating filtration systems, with the aim of minimizing head loss and maximizing system performance.It should be noted that the methodology employed in this study, which involved combining response surfaces and an optimized model, demonstrates the effectiveness of such an approach in understanding complex relationships between multiple factors.By visualizing the data in Figures 5-10 and applying the method of finding the optimal slope of the response surface, a comprehensive set of solutions was obtained, leading to the identification of the best possible combination of parameters.
The different combinations of lines in Figure 11 represent different optimal solutions obtained by adjusting parameters through model algorithms when applying the response surface method.These lines demonstrate the relationship between parameters and system performance, and help researchers determine the possible optimal combination of parameters.

Discussion
Although numerous studies have explored the factors influencing filter performance, few have comprehensively considered the interactions among multiple factors [27][28][29][30][31].By introducing interaction terms, this study not only unveils the complex relationships The different combinations of lines in Figure 11 represent different optimal solutions obtained by adjusting parameters through model algorithms when applying the response surface method.These lines demonstrate the relationship between parameters and system performance, and help researchers determine the possible optimal combination of parameters.

Discussion
Although numerous studies have explored the factors influencing filter performance, few have comprehensively considered the interactions among multiple factors [27][28][29][30][31].By introducing interaction terms, this study not only unveils the complex relationships among these factors but also enhances the predictive accuracy of the model.This aligns with the recent emphasis in academic journals on model complexity and predictive capability [32].
Furthermore, the findings of this study corroborate some of the latest literature, further validating the model's reliability.Several studies have pointed out the significant impact of irrigation flow rate on filter performance, which concurs with the predominant role of irrigation flow rate observed in this study.Additionally [31][32][33], this study identifies the important influence of self-cleaning flow rate and filter cartridge speed on head loss, echoing discussions in the literature about filter self-cleaning functions and cartridge design optimization [32].
However, compared to the existing literature, this study reveals some novel findings [34,35].For example, the interaction between irrigation flow rate and filter cartridge speed has a particularly significant effect on head loss, highlighting the need to pay closer attention to the combination of these two factors during irrigation system design.Moreover, the strongest contribution from the self-cleaning flow rate underscores the importance of self-cleaning functionality in filter design.
Lastly, this study acknowledges certain limitations compared to the existing literature.Specifically, it focuses on four primary test factors, while acknowledging that actual irrigation systems may involve other potential influencing factors.Future research could expand the model's applicability by considering additional factors and their interactions on filter performance.Moreover, validating the model's universal applicability and reliability across various types of filters would be a worthwhile pursuit [36].

•
Integrating four key test factors, including irrigation flow, filter cartridge speed, selfcleaning flow, and initial sand content, we developed a head loss model for pressureless mesh filters used in farmland irrigation.After refining the response surface, the model showed a significant improvement with a coefficient of determination of 98.61%.• The total irrigation flow had a 17.20% higher contribution than the self-cleaning flow rate, while the initial sand content had the smallest contribution at 0.0166.Additionally, the self-cleaning flow rate had a 46.18% higher contribution than the filter cartridge speed, and the filter cartridge speed had a 96.34% higher contribution than the initial sand content.• The interaction term between irrigation flow and filter cartridge speed had a contribution 218.25% higher than the interaction term between irrigation flow and initial sand content.The self-cleaning flow rate's self-interaction had the highest contribution, while the filter cartridge speed's self-interaction had the smallest contribution.

•
By optimizing the response surface and the model, the optimal parameters were determined to be an irrigation flow rate of 121.687 m 3 •h −1 , a filter cartridge speed of 1.331 r•min −1 , a self-cleaning flow rate of 19.980 m 3 •h −1 , and an initial sand content of 0.261 g•L −1 , resulting in a minimum head loss of 21.671 kPa.

Figure 3 .
Figure 3.This information provides insights into the characteristics of the sediment used in the experimental tests.

Figure 4 .
Figure 4. Analysis of the contribution parameters and confidence intervals of filter factors.

Figure 4 .
Figure 4. Analysis of the contribution parameters and confidence intervals of filter factors.

Agriculture 2024 , 20 Figure 5 .
Figure 5. Normalization of filter flow and cartridge speed interaction item to negative pressure value.

Figure 5 .
Figure 5. Normalization of filter flow and cartridge speed interaction item to negative pressure value.

Agriculture 2024 , 20 Figure 6 .
Figure 6.Normalization of filter flow and self-cleaning flow interaction to negative pressure values.

Figure 7 .
Figure 7. Normalization of filter flow and initial sediment interaction to negative pressure values.

Figure 6 .
Figure 6.Normalization of filter flow and self-cleaning flow interaction to negative pressure values.

Agriculture 2024 , 20 Figure 6 .
Figure 6.Normalization of filter flow and self-cleaning flow interaction to negative pressure values.

Figure 7 .
Figure 7. Normalization of filter flow and initial sediment interaction to negative pressure values.

Figure 7 .
Figure 7. Normalization of filter flow and initial sediment interaction to negative pressure values.

Figure 8 .
Figure 8. Normalization of cartridge speed and self-cleaning flow interaction term to the negative pressure value.

Figure 8 .
Figure 8. Normalization of cartridge speed and self-cleaning flow interaction term to the negative pressure value.
exhibits a gentle change, while Figure 10 reveals a more inclined response value influenced by the self-cleaning flow rate.Agriculture 2024, 14, x FOR PEER REVIEW 16 of 20 filter cartridge rotation speed or self-cleaning flow rate is held constant, the change in head loss relative to the initial sediment concentration is not significant.Figure 9 exhibits a gentle change, while Figure 10 reveals a more inclined response value influenced by the self-cleaning flow rate.(2)Figures 9 and 10 demonstrate that the corresponding interaction term has a minimal impact on head loss.In comparison to the intense change area, the curvature of the response surface tends to be higher in the edge state, indicating a more substantial change in head loss.

Figure 9 .
Figure 9. Normalization of cartridge speed and initial sediment interaction to negative pressure value.

Figure 10 .
Figure 10.Normalization of self-cleaning flow and initial sediment interaction to negative pressure value.

Figure 9 .
Figure 9. Normalization of cartridge speed and initial sediment interaction to negative pressure value.

Figure 9 .
Figure 9. Normalization of cartridge speed and initial sediment interaction to negative pressure value.

Figure 10 .
Figure 10.Normalization of self-cleaning flow and initial sediment interaction to negative pressure value.

Figure 10 .
Figure 10.Normalization of self-cleaning flow and initial sediment interaction to negative pressure value.
. The results highlight a substantial improvement in the determination coefficient, indicating a stronger correlation between the variables.
. The results highlight a substantial improvement in the determination coefficient, indicating a stronger correlation between the variables.

Figure 11 .
Figure 11.Cartridge speed and initial sediment interaction to negative pressure value.

Figure 11 .
Figure 11.Cartridge speed and initial sediment interaction to negative pressure value.

Table 1 .
Response surface base factor.The quantities of standardized factors for the test and the statistics of test results are presented in Table2.

Table 2 .
Standardized response surface experimental design and result statistics table.

Table 4 .
The standardized model.

Table 5 .
Contribution response of factors to indicator head loss.

Table 6 .
Modify the model.

Table 6
presents the optimized model derived from the response surfaces depicted in Figures 5-10

Table 6 .
Modify the model.

Table 6
presents the optimized model derived from the response surfaces depicted in Figures 5-10