WO2017147728A1 - 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品 - Google Patents

滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品 Download PDF

Info

Publication number
WO2017147728A1
WO2017147728A1 PCT/CN2016/000339 CN2016000339W WO2017147728A1 WO 2017147728 A1 WO2017147728 A1 WO 2017147728A1 CN 2016000339 W CN2016000339 W CN 2016000339W WO 2017147728 A1 WO2017147728 A1 WO 2017147728A1
Authority
WO
WIPO (PCT)
Prior art keywords
flow
flow channel
drip irrigation
flow path
emitter
Prior art date
Application number
PCT/CN2016/000339
Other languages
English (en)
French (fr)
Inventor
李云开
周博
冯吉
杨培岭
Original Assignee
中国农业大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 中国农业大学 filed Critical 中国农业大学
Priority to US15/773,948 priority Critical patent/US10750683B2/en
Priority to PCT/CN2016/000339 priority patent/WO2017147728A1/zh
Publication of WO2017147728A1 publication Critical patent/WO2017147728A1/zh

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
    • A01G25/00—Watering gardens, fields, sports grounds or the like
    • A01G25/02—Watering arrangements located above the soil which make use of perforated pipe-lines or pipe-lines with dispensing fittings, e.g. for drip irrigation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
    • Y02A40/22—Improving land use; Improving water use or availability; Controlling erosion

Definitions

  • the present invention proposes a flow channel design method for drip irrigation emitters based on the disadvantages of existing product development theory and process.
  • the method is designed in the design process.
  • the present invention develops a fractal flow drip emitter product with both hydraulic and anti-clogging properties, which is excellent in both sheet and cylindrical shape.
  • the anti-clogging performance evaluation parameter P is used to evaluate the anti-clogging performance of the drip emitter, and the higher the value of P, the better the anti-blocking performance;
  • the flow path configuration with the largest turbulence intensity is selected as the optimal flow path configuration.
  • step 6 the prediction model of the drip irrigation water flow coefficient k is formula (3): Among them, a, b, and c are undetermined coefficients, and L, W, and D are respectively length, width, and depth in the structural parameters of the flow channel.
  • the d in the calculation formula (4) of the anti-blocking evaluation index P is -1.345, and f is 3.489, that is, the formula (7):
  • One port of the fractal flow channel is in communication with the water inlet of the drip emitter, and the other port is connected to the water outlet of the drip emitter
  • the length of the fractal flow channel is 214.4 mm, the width is 1.27 mm, and the depth is 0.745 mm.
  • the tooth tip water-facing area and the heel-facing water-storing area of the inner wall of the flow channel It is an arc with a radius of 0.423mm.
  • the literature for the disclosure of the modified fractal-M flow path is: Li Yunkai, 2005, Experimental study and numerical simulation of the design of the flow pattern of the dripper and its flow characteristics (dissertation), Chapter 4, page 38), the amendment The fractal-M flow path is based on the fractal flow path designed on the Minkowski curve, that is, the fractal-M flow path is prototyped (as shown in (a) of Fig.
  • a drip irrigation emitter prototype primary circulation design method (ie step 5)-7)) is proposed, with the help of the optimal numerical simulation model (ie under real roughness conditions)
  • the simulation model determines the optimal flow path configuration and structural parameters (including length, width and depth) of the drip irrigation emitter, and clarifies the design method of the primary prototype structure of the drip irrigation emitter.
  • Figure 1 is a flow chart of the design method of the water channel structure of the emitter for flow design requirements; wherein, the multi-scenario simulation is a solid-liquid-gas three-phase flow simulation;
  • Figure 2 is a fractal-M flow path prototype and a modified fractal-M flow path
  • Figure 4 is a distribution diagram of the flow shear force near the inner wall of the flow channel of the fractal sheet type FE38# drip irrigation emitter;
  • Figure 5 is a distribution diagram of the flow shear force near the inner wall of the flow channel of the fractal cylindrical CE91# drip irrigation emitter;
  • Fig. 6 is a distribution diagram of the flow shear force near the inner wall of the flow passage after the secondary fine structure optimization of the drip irrigation emitter; wherein, the left figure is the optimized fractal sheet type FE38#, and the right picture shows the optimized fractal cylinder type CE91#. The figure shows the magnitude of the flow shear force on the inner wall of the inner wall of the flow channel with different colors;
  • Figure 8 is a three-dimensional stereotype of the optimized fractal cylindrical CE91# product
  • W, L and D in the formulas of the following examples represent the flow path width, length and depth of the drip emitter, respectively (as shown in Figure 9), in mm;
  • Q rating represents the nominal design flow of the drip emitter (ie flow design requirements)
  • the unit is L/h.
  • the water emitters of the present invention are drip irrigation emitters.
  • the simulation model under the initial roughness condition is to use both the RNG (renormalization group) k- ⁇ model and the VOF (fluid volume function) model to characterize the water flow and the particulate matter solid-liquid-gas inside the emitter flow channel configuration.
  • RNG normalization group
  • VOF fluid volume function
  • the inlet is the pressure inlet (pressure is 0.1 MPa)
  • the solid phase volume fraction is 0.9%
  • the gas phase volume fraction is 1.1%
  • the liquid phase volume fraction is 98%
  • the outlet is the pressure outlet (pressure It is 0 MPa).
  • all other fluid and solid contact surfaces are set to have no slip boundary; solved by standard wall function; the mean wall roughness is 869 nm; the numerical calculation uses the finite volume method to discretize the governing equation.
  • the pressure term adopts the second-order upwind style; the coupling of pressure and velocity is solved by SIMPLE algorithm; the residual value is used as the basis for convergence. When the exit flow is basically stable and the residual value is lower than 104, the iterative calculation is considered to be convergent.
  • the flow path configuration proposed by the plurality of users is an existing fractal flow path, a toothed flow path, a triangular flow path, a rectangular flow path, and a trapezoidal flow path;
  • the physical prototypes ie, the simulated objects have the following dimensions:
  • Fractal flow channel width 1mm, length 35mm, depth 0.73mm;
  • Toothed flow path width 1mm, length 35mm, depth 0.73mm, tooth height 0.84mm, The tooth angle is 54° and the tooth pitch is 1.42 mm;
  • Trapezoidal flow channel width 1mm, length 35mm, depth 0.73mm, tooth height 0.84mm, tooth angle 54°, tooth spacing 1.8mm;
  • Triangular flow path width 1mm, length 35mm, depth 0.73mm, tooth height 1.2mm, tooth angle 54°, tooth spacing 1.42mm;
  • Rectangular flow path width 1mm, length 35mm, depth 0.73mm.
  • the fractal flow channel is a modified fractal-M flow channel
  • the literature for disclosing the modified fractal-M flow channel is: Li Yunkai, 2005, experimental research and numerical simulation of the dripper fractal flow channel design and its flow characteristics ( The thesis), Chapter 4, page 38)
  • the fractal flow path ie the fractal-M flow path
  • the design of the flow channel requires that the fractal flow path be simplified according to the principle that the number and sequence of the energy dissipating units are substantially consistent with the prototype, that is, the modified fractal, while maintaining the width and depth of the flow path of the flow dissipating unit.
  • -M flow path (as shown in (b) of Figure 2).
  • the mean value of the flow shear force near the inner wall of the inner wall is obtained.
  • the simulation result is 0.5Pa, which is substituted into formula (1) to obtain the mean value of the roughness of the plugged material.
  • the simulation model under the initial roughness condition is substituted to obtain the simulation model under the true roughness condition of each flow channel configuration; the subsequent simulation is performed.
  • Figure 3 is a color map of the turbulence intensity (ie, anti-blocking performance) of different drip irrigation emitter flow paths directly output using FLUENT software.
  • the color bars of different colors represent different values of turbulence intensity.
  • Figure 3 shows that the turbulence intensity of the drip irrigation water channel configuration is the largest when the fractal flow channel is used (the turbulence intensity is an indicator of the anti-blocking performance, and the turbulence intensity represents The size of the anti-clogging performance, the greater the turbulence intensity, the higher the anti-clogging performance), so it is determined that the optimal flow path configuration of the drip irrigation emitter is a fractal flow channel (ie, the configuration of the emitter flow channel in FIG. 1 is completed).
  • step 5 In the range of values obtained in step 5), take 20 representative values in the range of the width W, the range of the length L, and the range of the depth D, and the representative values are equalized. Pitch selection;
  • the flow coefficient k of the drip irrigation emitter represented by each combination is obtained, and the flow channel structural parameters and flow coefficient k represented by each combination are fitted by excel, and the prediction models of the drip irrigation emitter flow coefficient k are obtained as follows:
  • step 6 According to the method of step 6), set 100 representative values, and combine the representative values from small to large in order, and arrange the combinations according to the flow path structure parameters from small to large, according to formula (5) and ( 6) Calculate the anti-blocking evaluation parameter P of each combination (the results are shown in Table 2 and Table 3), and then calculate the deviation between the P values of the adjacent two combinations, and the deviation between the P values of the adjacent two combinations is smaller than 1% (ie s is 1) of several combinations of smaller flow path structure parameters (ie In the ith combination of the time, i represents a different combination), take a certain number of combinations, and substitute the values of the structural parameters into the corresponding formula (8) or formula (9) in step 6) to calculate the flow coefficient k.
  • the combined flow path structure parameters with the smallest k value are determined as the optimal flow path structure parameters of the drip emitter;
  • the fractal flow channel Determined as the optimal flow channel structure parameter, the fractal flow channel has a length of 214.4 mm, a width of 1.27 mm, and a depth of 0.745 mm;
  • the FE38# and CE91# were determined as the original shape of the optimal runner structure of the emitter, and the second-level fine structure design was carried out.
  • Fig. 4 and Fig. 5 are the results of color maps of the flow shear force distribution on the inner wall near the wall of different drip irrigation emitters using FLUENT software.
  • the different colors of the color column represent the values of the flow shear force of the inner wall near the wall. .
  • the flow shear force of the inner wall near the wall is in the water-carrying area and the heel-watering area at the inner wall of the flow channel, and the shearing force is between 0.2 and 0.4 Pa (marked by an arrow, The inventors of the present application have found that the flow shear force of the inner wall near the wall is in the interval of 0.2-0.4 Pa, and the blockage material grows fastest, so the occurrence of this interval is avoided. Therefore, the inner wall is optimized by the vortex wall optimization design method. The location is optimized as follows:
  • the water-pointing area of the tooth tip of the flow channel and the water-contact area of the tooth heel are designed as circular arcs with the same or the same shape and size as the outer edge of the vortex.
  • the tip of the flow channel is optimized for the water-facing area and the tooth-to-water area.
  • the optimized drip irrigation emitter is subjected to solid-liquid-gas three-phase flow simulation, and the flow shear force of the inner wall of the drip irrigation emitter is near wall surface.
  • the optimized flow path structure is determined as the optimal flow path structure of the drip irrigation emitter.
  • the FE38# drip emitter and the CE91# drip emitter are respectively optimized by arcs with circular arc radius of 1, 2, and 1/3 of the flow path width, respectively, and the fractal flow obtained by step 3) is used.
  • the simulation model under the true roughness condition of the road simulates the solid-liquid-gas three-phase flow of the internal flow and the movement of the particulate matter in the optimized FE38# and CE91# drip irrigation emitters respectively;
  • the shear force of the inner wall of the inner wall of the CE91# drip irrigation emitter is not 0.2- Within the 0.4 Pa interval (the right image in Figure 6), the optimized flow path structure of the CE91# drip emitter at this time is determined as the optimal flow path structure of the drip emitter.
  • the optimal flow path structure of the two drip irrigation emitters obtained in step 8) is the physical prototype, that is, the simulation object, and the simulation model under the real roughness condition of the fractal flow channel obtained in the step 3) is simulated, and the flow is directly outputted in the FLUENT.
  • the optimized flow value after simulation is 1.58L/h, which is 1.25% deviation from the rated design flow of drip irrigation emitters of 1.6L/h, indicating that the boundary optimization (ie, the optimal design of the vortex wall-washing channel structure) has no flow rate. Have an impact.
  • the optimal flow path structure of the two drip irrigation emitters obtained in step 8) is UG NX series software for 3D shaping, high-precision mold development, and finally high-precision mold is developed according to user application requirements (accuracy is not less than ⁇ 5 ⁇ m) ), select other materials such as drip irrigation pipe material, wall thickness, etc. to realize the industrialization of new products of drip irrigation emitters, and obtain two kinds of drip irrigation emitter products:
  • Drip irrigation emitter product A the type is sheet type, the length of the fractal flow channel is 39.567mm, the width is 0.824mm, the depth is 0.7555mm, the depth of the inner wall of the flow channel is the tip of the water and the heel of the water is 0.412mm.
  • Arc as shown in Figure 7;
  • Drip irrigation emitter product B type is cylindrical, fractal flow channel each flow channel structural unit has a length of 214.4mm, a width of 1.27mm and a depth of 0.745mm; the tip of the flow channel has a tip-to-water area and teeth
  • the arc is a circular arc with a radius of 0.423 mm; as shown in Figure 8.
  • the drip emitter products A and B shown in Fig. 7 and Fig. 8 were applied in actual field, and it was found through testing that the two drip emitters verified that the drip emitters developed by the method provided by the invention have extremely high Hydraulic performance (product flow index between 0.50-0.52) and anti-clogging performance (system continuous running time of 680-840h).

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental Sciences (AREA)
  • Architecture (AREA)
  • Soil Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Civil Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Software Systems (AREA)
  • Nozzles (AREA)

Abstract

一种滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品，该方法在设计过程中综合考虑了滴灌灌水器的水力性能和抗堵塞性能，主要包括模拟方法建立、流道构型选择、结构参数确定、流道边界优化等阶段；借助最优数值模拟模型确定了滴灌灌水器流道的最优流道构型及结构参数取值范围，明确了滴灌灌水器初级雏形结构设计方法；提出了一种滴灌灌水器流道边界的漩涡洗壁优化设计方法，确定了边界优化设计的控制阈值范围，确定了二级精细结构设计方法；应用该方法结合分形几何理论，设计了片式和圆柱式两种分形流道灌水器产品，具有极高的水力性能和抗堵塞性能。

Description

滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品 技术领域
本发明涉及高效节水灌溉技术领域，特别是一种面向流量设计需求的滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品。
背景技术
滴灌系统的滴灌灌水器(即滴头)堵塞问题解决的好坏程度直接影响系统灌水均匀度、使用寿命和运行效益。尤其是目前全球水资源紧缺和水污染严重并存，再生水、地表微污染水、高含沙水、微咸水等非常规水源也常常回用于农田灌溉，其复杂的水质情况极大的增加了滴灌灌水器堵塞的风险，也使得堵塞机理更为复杂。
解决滴灌灌水器堵塞问题，关键就是提升自身抗堵塞能力。对滴灌灌水器结构进行优化设计，提升滴灌灌水器对颗粒物的输移能力及壁面的自清洗能力是提升滴灌灌水器自身抗堵塞能力的有效途径。滴灌灌水器结构设计关键在于消能流道设计。迷宫流道是消能流道设计中最为常见的形式，但由于迷宫流道的复杂性，消能机理与设计理论研究的缺失，导致传统的滴灌灌水器设计开发过程缺乏设计知识的积累，主要依托注塑模具设计与制造工艺，仿造国外先进的滴头结构，一个产品的定型需要若干次反复修改，尤其在模具加工阶段，并且若零件试验性能不能满足产品要求，必须从模具修改开始，甚至模具还有报废的风险。这样的设计开发周期一般为4-5个月，成本一般在5万元以上，不仅费时费力，造价极高，最为重要的是产品精度得不到保证。
也有部分专家学者提出了滴灌灌水器的优化设计方法。如：中国农业大学杨培岭等提出了一种抗堵塞滴灌灌水器设计方法，通过减小流道特征部位结构尺寸对流道结构进行了优化(中国专利，申请号为CN200710063794.6)；中国水利水电科学研究院王建东等提出一种低 压滴灌锯齿形灌水器流道结构优化设计方法，该方法通过试验与模拟相结合的方式对灌水器结构进行优化设计(中国专利，申请号为CN201220669437.0)。上述方法均只是着重于对灌水器流道边界进行优化设计，但是并未确定流道边界设计的控制阈值范围，更未提出面向流量设计需求的一整套滴灌灌水器流道结构设计方法。
发明内容
针对现有技术中存在的缺陷，本发明在已有产品开发理论与过程中存在弊端的基础上，提出了一种面向流量设计需求的滴灌灌水器流道结构设计方法，该方法在设计过程中综合考虑了滴灌灌水器的水力性能和抗堵塞性能，本发明应用该方法开发了片式和圆柱式两种水力性能与抗堵塞性能俱优的分形流道滴灌灌水器产品，本发明使得在国际范围内滴灌灌水器设计开发产生质的提升。
为达到以上目的，本发明采取的技术方案是：
一种面向流量设计需求的滴灌灌水器流道结构设计方法，包括如下步骤：
步骤1，根据用户所提出的流量设计需求，在若干用户所提出的流道构型中，确定最优流道构型；
步骤2，通过试算抗堵塞性能评估参数P的值和流量系数k的值，来确定最优流道构型的最优流道结构参数；
所述流道结构参数包括；宽度W、长度L和深度D；
所述抗堵塞性能评估参数P用于评估滴灌灌水器的抗堵塞性能的优劣，P的值越高，抗堵塞性能越好；
所述流量系数k用于评估滴灌灌水器的出流流量随进口压力波动的敏感程度，k的值越小，滴灌灌水器的出流流量随进口压力波动的敏感程度越小，水力性能越好；
步骤3，通过最优流道结构参数，得到最优流道结构原形。
所述结构参数的单位为mm。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，试算抗堵塞性能评估参数P的值和流量系数k的值时，选择若干组流道结构参数参与试算。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，所述步骤1包括如下步骤：
1)使用初始粗糙度条件下的模拟模型，分别以若干用户所提出的流道构型(如分形流道、齿形流道、三角形流道、矩形流道和梯形流道)为模拟对象，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟，得到各不同的流道构型的内壁近壁面流动剪切力均值
Figure PCTCN2016000339-appb-000001
(单位为Pa)模拟结果；
所述初始粗糙度条件下的模拟模型中的边壁粗糙度均值为默认值；
2)将步骤1)得到的
Figure PCTCN2016000339-appb-000002
的模拟结果分别代入公式(1)：
Figure PCTCN2016000339-appb-000003
Figure PCTCN2016000339-appb-000004
中，得到各不同的流道构型的堵塞物质粗糙度均值
Figure PCTCN2016000339-appb-000005
(单位为μm)；
所述公式(1)为本发明发明人通过对再生水、高含沙水、地表微咸水、高含沙水和地表微咸水混配水等多水源滴灌条件以及常用的片式、圆柱式、单翼迷宫式、贴条式等多类型灌水器抗堵塞试验，得到的在滴灌灌水器堵塞程度为5％时滴灌灌水器流道的内壁近壁面流动剪切力均值
Figure PCTCN2016000339-appb-000006
和堵塞物质粗糙度均值
Figure PCTCN2016000339-appb-000007
的相关关系；
所述公式(1)适用于不同类型(如片式、圆柱式、单翼迷宫式)的滴灌灌水器及滴灌灌水器不同的流道构型；
3)将步骤2)得到的各
Figure PCTCN2016000339-appb-000008
作为对应的流道构型的边壁粗糙度均值，代入初始粗糙度条件下的模拟模型，得到各流道构型的真实粗糙度条件下的模拟模型；
4)使用步骤3)得到的各流道构型的真实粗糙度条件下的模拟模型，分别以若干用户所提出的流道构型为模拟对象，对流道内部的 水流及颗粒物运动进行固-液-气三相流动模拟，得到各流道构型内部的湍流强度；
选择湍流强度最大的流道构型确定为最优流道构型。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，所述步骤2包括如下步骤：
5)使用最优流道构型对应的真实粗糙度条件下的模拟模型，在满足流道消能需求(即不产生射流)和流态指数x为0.50-0.55，且流道消能需求和流态指数x随所述最优流道构型的流道结构参数变化不敏感的前提下，确定所述最优流道构型的流道结构参数的取值范围(即控制阈值)，具体包括：宽度W的取值范围、长度L的取值范围和深度D的取值范围；
6)在步骤5)得到的取值范围内，分别在宽度W的取值范围、长度L的取值范围和深度D的取值范围中，取相同个数的代表值(代表值的个数越多，最后得到的最优流道结构参数越准确)，且各代表值按均等间距选取；
使用最优流道构型对应的真实粗糙度条件下的模拟模型，在不同进口压力H条件下，分别模拟各代表值从小到大依次组合所代表的滴灌灌水器出流流量Q，得到公式(2)：Q＝kHx，公式(2)中的Q为出流流量，k为流量系数，H为进口压力，x为流态指数；
得到各组合所代表的滴灌灌水器的流量系数k，将各组合所代表的流道结构参数与流量系数k进行拟合，得到滴灌灌水器流量系数k的预报模型；
7)将步骤6)各组合或按照步骤6)的方法设置的各组合按照流道结构参数从小到大的顺序排列，计算各组合抗堵塞评估指数P值， 再计算相邻两个组合的P值间的偏差(即
Figure PCTCN2016000339-appb-000009
)，从相邻两个组合的P值间的偏差小于s％(s为小于100的正数，可取s为1)的流道结构参数较小的若干组合(即
Figure PCTCN2016000339-appb-000010
时的第i个组合，i代表不同组合)中，选取一定数量的组合，将其结构参数的值分别代入步骤6)滴灌灌水器流量系数k的预报模型中，计算流量系数k，将k值最小的组合中的流道结构参数确定为滴灌灌水器最优流道结构参数。
在最优流道构型的流道结构参数的取值范围(即控制阈值)内，水力性能参数即流态指数x仅受流道构型的影响，对于流道结构参数的改变较为不敏感，但对滴灌灌水器抗堵塞性能产生显著影响，因此本发明主要通过试算抗堵塞性能评估参数P值来确定最优的滴灌灌水器流道结构参数，进而得到滴灌灌水器流道最优结构原形。具体为选择多组长度、宽度、深度组合进行抗堵塞评估参数P的试算。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，步骤6)滴灌灌水器流量系数k的预报模型为公式(3)：
Figure PCTCN2016000339-appb-000011
其中，a、b、c为待定系数，L、W、D分别为流道结构参数中的长度、宽度和深度。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，步骤7)中抗堵塞评估指数P值的计算方法如下：
Figure PCTCN2016000339-appb-000012
公式(4)中的d和f为待定系数；L、W、D分别为流道的长度、 宽度和深度，单位为mm；Q额定代表滴灌灌水器额定设计流量即用户所提出的流量设计需求，单位为L/h。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，步骤7)中，当滴灌灌水器的类型为片式时，所述抗堵塞评估指数P的计算公式(4)中的d为-15.723，f为21.093，即为式(5)：
Figure PCTCN2016000339-appb-000013
Figure PCTCN2016000339-appb-000014
当滴灌灌水器的类型为圆柱式时，所述抗堵塞评估指数P的计算公式(4)中的d为-0.167，f为1.243，即为式(6)：
Figure PCTCN2016000339-appb-000015
Figure PCTCN2016000339-appb-000016
当滴灌灌水器的类型为单翼迷宫式时，所述抗堵塞评估指数P的计算公式(4)中的d为-1.345，f为3.489，即为式(7)：
Figure PCTCN2016000339-appb-000017
Figure PCTCN2016000339-appb-000018
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，还包括步骤4：
对步骤3所述最优流道结构原形进行二级精细结构设计，将所述最优流道结构原形的流道中内壁近壁面流动剪切力处在0.2-0.4Pa间的内壁位置优化为内壁近壁面流动剪切力不在0.2-0.4Pa间。
所述内壁位置分为齿尖迎水区、齿跟迎水区、齿尖背水区和齿跟背水区。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，步骤4所述优化使用漩涡洗壁优化设计方法：
分析步骤3的所述最优流道结构原形流道内漩涡分布，根据漩涡 外边缘的形状和大小，将流道中内壁近壁面流动剪切力处在0.2-0.4Pa间的内壁位置设计为与漩涡外边缘形状和大小相近或相同的圆弧。
该漩涡洗壁优化设计方法可以使漩涡充分发展，提升水流对流道壁面的自清洗能力，进而提升滴灌灌水器的抗堵塞能力。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，所述流道中内壁近壁面流动剪切力，是使用最优流道结构原形对应的真实粗糙度条件下的模拟模型或最优流道构型的真实粗糙度条件下的模拟模型，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟得到的。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，所述模拟使用FLUENT软件进行。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，所述初始粗糙度条件下的模拟模型为同时使用RNG(重整化群)k-ε模型和VOF(流体体积函数)模型表征灌水器流道构型内部的水流及颗粒物运动固-液-气三相流动湍流模拟模型。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，初始粗糙度条件下的模拟模型中的边壁粗糙度均值以经验值作为默认值；所述经验值具体可为869nm。
在上述面向流量设计需求的滴灌灌水器流道结构设计方法中，使用所述初始粗糙度条件下的模拟模型进行模拟时，固相体积分数、气相体积分数和液相体积分数根据用户实际使用的水源水质的不同而不同，三者之和为1；
在实施例1中，固相体积分数、气相体积分数和液相体积分数分别为0.9％，1.1％，98％。
本发明保护上述任一所述面向流量设计需求的滴灌灌水器流道结构设计方法在生产滴灌灌水器产品中的应用，即将上述任一所述面向流量设计需求的灌水器流道结构设计方法得到的最优流道构型的最优流道结构，采用软件如UG NX系列软件进行3维定型，最终依据用户应用需求，开发高精度模具(精度不低于±5μm)，选定滴灌管材料、壁厚等其它参数，实现灌水器新产品的产业化，获得滴灌灌水器产品。
本发明保护所述应用中的滴灌灌水器产品。
所述滴灌灌水器的流道构型为分形流道，
所述分形流道为修正的分形-M流道，
所述分形流道的一端口与所述滴灌灌水器的进水口相通，另一端口与所述滴灌灌水器的出水口相通，
当所述滴灌灌水器的类型为片式时，所述分形流道的长度为39.567mm，宽度为0.824mm，深度为0.7555mm，流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.412mm的圆弧；
当所述滴灌灌水器的类型为圆柱式时，所述分形流道的长度为214.4mm，宽度为1.27mm，深度为0.745mm，流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.423mm的圆弧。
公开所述修正的分形-M流道的文献为：李云开，2005，滴头分形流道设计及其流动特性的试验研究与数值模拟(学位论文)，第四章第38页)，所述修正的分形-M流道是以Minkowski曲线为基础设计的分形流道即分形-M流道为原型(如图2中的(a)图所示)，考虑到灌水器流道的设计要求，在保持流道消能单元流道的宽、深不变的条件下，根据消能单元数目与顺序排列基本符合原型的原则对分形流道进行了简化，得到所述修正的分形-M流道(如图2中的(b)图 所示)。
本发明所述的流道长度为流道中心线长度，即指流道边壁垂直距离中点的连线。
本发明的有益效果如下：
本发明可以有效解决以下几方面的问题：
(1)提出了一种面向流量设计需求的滴灌灌水器流道结构循环逐级优化设计方法。该方法在设计过程中综合考虑了滴灌灌水器的水力性能(流态指数x)和抗堵塞性能(湍流强度、抗堵塞评估指数P)，主要包括模拟方法建立(即步骤1)-3))、流道构型选择(即步骤4))、结构参数确定(即步骤5)-7))、流道边界优化(即步骤4)等阶段。
(2)依据数值模拟和大量的试验测试统计分析结果，提出了一种滴灌灌水器雏形初级循环设计方法(即步骤5)-7))，借助最优数值模拟模型(即真实粗糙度条件下的模拟模型)确定了滴灌灌水器最优流道构型及结构参数(包括长度、宽度和深度)取值范围，明确了滴灌灌水器初级雏形结构设计方法。
(3)提出了一种滴灌灌水器流道边界的漩涡洗壁优化设计方法，确定了滴灌灌水器边界优化设计的控制阈值范围(即内壁近壁面流动剪切力不在0.2-0.4Pa)，确定了滴灌灌水器二级精细结构设计方法。
(4)应用本发明所提出的设计方法结合分形几何理论，设计了片式和圆柱式两种分形流道灌水器产品(即产品A和B)，具有极高的水力性能(产品流态指数介于0.50-0.52)和抗堵塞性能(系统持续运行时间可达680-840h)。
附图说明
本发明有如下附图：
图1为面向流量设计需求的灌水器流道结构设计方法流程图；其中，多情境模拟即为固-液-气三相流动模拟；
图2为分形-M流道原型及修正的分形-M流道；
图3为滴灌灌水器不同流道构型的湍流强度(即抗堵塞性能)结果；
图4为分形片式FE38#滴灌灌水器流道内壁近壁面流动剪切力分布图；
图5为分形圆柱式CE91#滴灌灌水器流道内壁近壁面流动剪切力分布图；
图6为滴灌灌水器二级精细结构优化后流道内壁近壁面流动剪切力分布图；其中，左图为优化后的分形片式FE38#，右图为优化后的分形圆柱式CE91#，图中用不同颜色代表流道内壁近壁面流动剪切力的大小；
图7为优化后的分形片式FE38#产品3维定型图；
图8为优化后的分形圆柱式CE91#产品3维定型图；
图9为本发明滴灌灌水器流道结构参数示意图，其中，W、L和D分别代表滴灌灌水器流道宽度、长度和深度。
具体实施方式
以下结合附图对本发明作进一步详细说明。
以下实施例公式中的W、L和D分别代表滴灌灌水器流道宽度、长度和深度(如图9所示)，单位为mm；Q额定代表滴灌灌水器额定设计流量(即流量设计需求)，单位为L/h。
本发明所述灌水器均为滴灌灌水器。
实施例1、应用面向流量设计需求的滴灌灌水器流道结构设计方法进行滴灌灌水器设计
用户所提出的Q额定为1.6L/h，滴灌灌水器类型为片式和圆柱式。设计过程如图1所示，具体如下：
一、考虑边壁粗糙度条件下灌水器内部固-液-气三相流动模拟模型建立
1)使用初始粗糙度条件下的模拟模型，分别以若干用户所提出 的流道构型为模拟对象，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟，得到各不同的流道构型的内壁近壁面流动剪切力均值
Figure PCTCN2016000339-appb-000019
(单位为Pa)模拟结果；
2)将步骤1)得到的将
Figure PCTCN2016000339-appb-000020
的模拟结果分别代入公式(1)：
Figure PCTCN2016000339-appb-000021
Figure PCTCN2016000339-appb-000022
中，得到各不同的流道构型的堵塞物质粗糙度均值
Figure PCTCN2016000339-appb-000023
(单位为μm)；
3)将步骤2)得到的各
Figure PCTCN2016000339-appb-000024
(单位为μm)作为对应的流道构型的边壁粗糙度均值，代入初始粗糙度条件下的模拟模型，得到各流道构型的真实粗糙度条件下的模拟模型；
所述初始粗糙度条件下的模拟模型为同时使用RNG(重整化群)k-ε模型和VOF(流体体积函数)模型表征灌水器流道构型内部的水流及颗粒物运动固-液-气三相流动湍流模拟模型；
使用所述初始粗糙度条件下的模拟模型进行模拟时，具体设置如下：
在流场计算中，初始条件设置：进口为压力进口(压力为0.1MPa)，固相体积分数为0.9％，气相体积分数为1.1％，液相体积分数为98％，出口为压力出口(压力为0MPa)。除了计算域的进水口与出水口，其它所有流体和固体接触的面均设置为无滑移边界；通过标准壁面函数来求解；边壁粗糙度均值为869nm；数值计算采用有限体积法离散控制方程；压力项采用二阶迎风格式；压力与速度的耦合采用SIMPLE算法求解；以残差值作为是否收敛的依据，当出口流量基本稳定且残差值低于104时，认为迭代计算达到收敛。
以上模型模拟使用FLUENT软件(软件版本号为6.3，软件开发商为ansys公司)进行。
所述若干用户所提出的流道构型为现有的分形流道、齿形流道、三角形流道、矩形流道和梯形流道；其物理原型(即模拟对象)的尺寸分别如下：
分形流道：宽度为1mm，长度为35mm，深度为0.73mm；
齿形流道：宽度为1mm，长度为35mm，深度0.73mm，齿高为0.84mm， 齿角度为54°，齿间距为1.42mm；
梯形流道：宽度为1mm，长度为35mm，深度0.73mm，齿高为0.84mm，齿角度为54°，齿间距为1.8mm；
三角形流道：宽度为1mm，长度为35mm，深度0.73mm，齿高为1.2mm，齿角度为54°，齿间距为1.42mm；
矩形流道：宽度为1mm，长度为35mm，深度为0.73mm。
其中，所述分形流道为修正的分形-M流道(公开该修正的分形-M流道的文献为：李云开，2005，滴头分形流道设计及其流动特性的试验研究与数值模拟(学位论文)，第四章第38页)，是以Minkowski曲线为基础设计的分形流道即分形-M流道为原型(如图2中的(a)图所示)，考虑到滴灌灌水器流道的设计要求，在保持流道消能单元流道的宽度、深度不变的条件下，根据消能单元数目与顺序排列基本符合原型的原则对分形流道进行了简化，即修正的分形-M流道(如图2中的(b)图所示)。
通过分析分形流道内流动特征，得到内壁近壁面流动剪切力均值
Figure PCTCN2016000339-appb-000025
的模拟结果为0.5Pa，代入公式(1)，得到堵塞物质粗糙度均值
Figure PCTCN2016000339-appb-000026
为1.4μm，代入初始粗糙度条件下的模拟模型，得到各流道构型的真实粗糙度条件下的模拟模型；进行后续模拟。
二、滴灌灌水器初级雏形结构设计
4)使用步骤3)得到的各流道构型的真实粗糙度条件下的模拟模型，分别以若干用户所提出的流道构型为模拟对象，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟，得到各流道构型内部的湍流强度；结果如图3所示；
图3为使用FLUENT软件直接输出的不同滴灌灌水器流道构型的湍流强度(即抗堵塞性能)彩图结果，不同颜色的色柱代表不同的湍流强度的数值。
图3显示，滴灌灌水器流道构型为所述分形流道时的湍流强度最大(湍流强度是一个体现抗堵塞性能的指标，湍流强度的大小就代表 抗堵塞性能的大小，湍流强度越大，抗堵塞性能越高)，因此确定滴灌灌水器最优流道构型为分形流道(即完成图1中的灌水器流道构型选择)。
5)使用步骤3)得到的分形流道的真实粗糙度条件下的模拟模型，在满足流道消能需求(即不产生射流)和流态指数x为0.50-0.55时，且流道消能需求和流态指数x随结构参数(流道的长度L、宽度W、深度D)变化不敏感的前提下，确定分形片式和分形圆柱式流道构型的流道结构参数控制阈值(即流道宽度W的取值范围、长度L的取值范围和深度D的取值范围)，结果如表1所示。
6)在步骤5)得到的取值范围内，分别在宽度W的取值范围、长度L的取值范围和深度D的取值范围中，取20个的代表值，且各代表值按均等间距选取；
使用分形流道对应的真实粗糙度条件下的模拟模型，在不同进口压力H(具体为0.01、0.03、0.05、0.07、0.09、0.1、0.11、0.13、0.15Mpa)条件下，分别模拟各代表值从小到大依次组合所代表的滴灌灌水器出流流量Q，得到公式(2)：Q＝kHx，公式(2)中的Q为出流流量，k为流量系数，H为进口压力，x为流态指数；
得到各组合所代表的滴灌灌水器的流量系数k，将各组合所代表的流道结构参数与流量系数k采用excel进行拟合，得到滴灌灌水器流量系数k的预报模型分别为：
分形片式滴灌灌水器：
Figure PCTCN2016000339-appb-000027
分形圆柱式滴灌灌水器：
Figure PCTCN2016000339-appb-000028
7)按照步骤6)的方法设置设100个代表值，将各代表值从小到大依次进行参数间组合，将各组合按照流道结构参数从小到大的顺序排列，根据公式(5)和(6)计算各组合的抗堵塞评估参数P(结果如表2和表3所示)，再计算相邻两个组合的P值间的偏差，从相邻 两个组合的P值间的偏差小于1％(即s为1)的流道结构参数较小的若干组合(即
Figure PCTCN2016000339-appb-000029
时的第i个组合，i代表不同组合)中，取一定数量的组合，将其结构参数的值分别代入步骤6)相应的公式(8)或公式(9)中，计算流量系数k，将k值最小的组合的流道结构参数确定为滴灌灌水器最优流道结构参数；
分形片式滴灌灌水器：
Figure PCTCN2016000339-appb-000030
分形圆柱式滴灌灌水器：
Figure PCTCN2016000339-appb-000031
分形片式滴灌灌水器当试算至第38个组合(即i＝38)至第99个组合(即i＝99)时，
Figure PCTCN2016000339-appb-000032
均小于1％，将第38个组合至第42个组合结构参数的值分别代入公式(8)中计算kFE即k值，将k值最小的第38个组合即FE38#的流道结构参数确定为最优流道结构参数，其分形流道的长度为39.567mm，宽度为0.824mm，深度为0.7555mm；
分形圆柱式滴灌灌水器当试算至第91个组合(即i＝91)至第99个组合(即i＝99)时，
Figure PCTCN2016000339-appb-000033
均小于1％，将第91个组合至第96个组合结构参数的值分别代入公式(9)中计算kCE即k值，将k值最小的第91个组合即CE91#的流道结构参数确定为最优流道结构参数，其分形流道的长度为214.4mm，宽度为1.27mm，深度为0.745mm；
将FE38#和CE91#确定为灌水器最优流道结构原形，进行下一步的二级精细结构设计。
表1.滴灌灌水器流道结构参数阈值
类型 长度(mm) 宽度(mm) 深度(mm)
分形片式 32.5-51.6 0.75-0.95 0.70-0.85
分形圆柱式 128.0-224.0 1.00-1.30 0.70-0.75
表2.分形片式滴灌灌水器抗堵塞性能特征值结果
Figure PCTCN2016000339-appb-000034
表3.分形圆柱式滴灌灌水器抗堵塞性能特征值结果
Figure PCTCN2016000339-appb-000035
Figure PCTCN2016000339-appb-000036
三、灌水器二级精细结构设计
8)漩涡洗壁流道结构优化设计
使用步骤3)得到的分形流道的真实粗糙度条件下的模拟模型，分别以FE38#和CE91#为模拟对象，对滴灌灌水器流道内部水流及颗粒物运动进行固-液-气三相流动模拟，得出滴灌灌水器内壁近壁面流动剪切力分布情况，结果如图4和图5所示。图4和图5为使用FLUENT软件直接输出的不同滴灌灌水器流道构型的内壁近壁面流动剪切力分布的彩图结果，色柱的不同颜色代表不同内壁近壁面流动剪切力的数值。
从图4和图5中可以看出，在流道内壁位置的齿尖迎水区和齿跟迎水区，内壁近壁面流动剪切力处在0.2-0.4Pa间(已用箭头标出，本申请的发明人经研究发现，内壁近壁面流动剪切力在0.2-0.4Pa区间内，堵塞物质生长最快，所以要避免这个区间的出现)，故采用漩涡洗壁优化设计方法对该内壁位置进行优化，具体方法如下：
分析流道内漩涡分布，根据漩涡外边缘的形状和大小，将流道的齿尖迎水区和齿跟迎水区设计为与漩涡外边缘形状和大小相近或相同的圆弧，对滴灌灌水器流道的齿尖迎水区和齿跟迎水区进行优化，
然后使用步骤3)得到的分形流道的真实粗糙度条件下的模拟模型对该优化后的滴灌灌水器进行固-液-气三相流动模拟，当滴灌灌水器的内壁近壁面流动剪切力均不处于0.2-0.4Pa区间内时，将该优化后的流道结构确定为滴灌灌水器最优流道结构，
当滴灌灌水器的内壁近壁面流动剪切力还存在0.2-0.4Pa区间时，按照上述方法进行再次优化和模拟，直到滴灌灌水器的内壁近壁面流 动剪切力均不处于0.2-0.4Pa区间内。
本实施例首先分别采用圆弧半径为流道宽度的1、1/2、1/3的圆弧对FE38#滴灌灌水器和CE91#滴灌灌水器分别进行优化，使用步骤3)得到的分形流道的真实粗糙度条件下的模拟模型对本次优化的FE38#和CE91#滴灌灌水器流道内部水流及颗粒物运动分别进行固-液-气三相流动模拟；
结果发现：当采用圆弧半径为流道宽度的1/2的圆弧对流道的齿尖迎水区和齿跟迎水区进行优化时，FE38#滴灌灌水器内壁近壁面流动剪切力均不处于0.2-0.4Pa区间内(图6中的左图)，将此时的FE38#滴灌灌水器经过优化后的流道结构确定为滴灌灌水器最优流道结构；
当采用圆弧半径为流道宽度的1/3的圆弧对流道的齿尖迎水区和齿跟迎水区进行优化时，CE91#滴灌灌水器内壁近壁面剪切力均不处于0.2-0.4Pa区间内(图6中的右图)，将此时的CE91#滴灌灌水器经过优化的流道结构确定为滴灌灌水器最优流道结构。
9)滴灌灌水器设计流量校核
以步骤8)得到的两种滴灌灌水器最优流道结构为物理原型即模拟对象，采用步骤3)得到的分形流道的真实粗糙度条件下的模拟模型进行模拟，直接在FLUENT中输出流量，得到优化后模拟的流量值，为1.58L/h，与滴灌灌水器额定设计流量1.6L/h的偏差为1.25％，说明边界优化(即漩涡洗壁流道结构优化设计)对流量基本没有产生影响。
10)滴灌灌水器产品定型
分别以步骤8)获得的两种滴灌灌水器最优流道结构，采用UG NX系列软件进行3维定型，开发高精度模具，最终依据用户应用需求，开发高精度模具(精度不低于±5μm)，选定滴灌管材料、壁厚等其它参数，实现滴灌灌水器新产品的产业化，获得两种滴灌灌水器产品：
滴灌灌水器产品A：类型为片式，分形流道的长度为39.567mm，宽度为0.824mm，深度为0.7555mm流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.412mm的圆弧；如图7所示；
滴灌灌水器产品B：类型为圆柱式，分形流道每个流道结构单元的长度为长度为214.4mm，宽度为1.27mm，深度为0.745mm；流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.423mm的圆弧；如图8所示。
将图7和图8所示的滴灌灌水器产品A和B在实际大田中进行应用，通过测试发现这两个滴灌灌水器验证了应用本发明所提供的方法开发的滴灌灌水器具有极高的水力性能(产品流态指数介于0.50-0.52)和抗堵塞性能(系统持续运行时间达680-840h)。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (10)

  1. 一种面向流量设计需求的滴灌灌水器流道结构设计方法，包括如下步骤：
    步骤1，根据用户所提出的流量设计需求，在若干用户所提出的流道构型中，确定最优流道构型；
    步骤2，通过试算抗堵塞性能评估参数P的值和流量系数k的值，来确定最优流道构型的最优流道结构参数；
    所述流道结构参数包括：宽度W、长度L和深度D；
    所述抗堵塞性能评估参数P用于评估滴灌灌水器的抗堵塞性能的优劣，P的值越高，抗堵塞性能越好；
    所述流量系数k用于评估滴灌灌水器的出流流量随进口压力波动的敏感程度，k的值越小，滴灌灌水器的出流流量随进口压力波动的敏感程度越小，水力性能越好；
    步骤3，通过最优流道结构参数，得到最优流道结构原形。
  2. 如权利要求1所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：所述步骤1包括如下步骤：
    1)使用初始粗糙度条件下的模拟模型，分别以若干用户所提出的流道构型为模拟对象，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟，得到各不同的流道构型的内壁近壁面流动剪切力均值
    Figure PCTCN2016000339-appb-100001
    拟结果；
    所述初始粗糙度条件下的模拟模型中的边壁粗糙度均值为默认值；
    2)将步骤1)得到的
    Figure PCTCN2016000339-appb-100002
    的模拟结果分别代入公式(1)：
    Figure PCTCN2016000339-appb-100003
    Figure PCTCN2016000339-appb-100004
    中，得到各流道构型的堵塞物质粗糙度均值
    Figure PCTCN2016000339-appb-100005
    3)将步骤2)得到的
    Figure PCTCN2016000339-appb-100006
    为对应的流道构型的边壁粗糙度均值，代入初始粗糙度条件下的模拟模型，得到各流道构型的真实粗糙度条件下的模拟模型；
    4)使用步骤3)得到的各流道构型的真实粗糙度条件下的模拟模型，分别以若干用户所提出的流道构型为模拟对象，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟，得到各流道构型内部的湍流强度；
    选择湍流强度最大的流道构型确定为最优流道构型。
  3. 如权利要求1所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：所述步骤2包括如下步骤：
    5)使用最优流道构型对应的真实粗糙度条件下的模拟模型，在满足流道消能需求和流态指数x为0.50-0.55，且流道消能需求和流态指数x随所述最优流道构型的流道结构参数变化不敏感的前提下，确定所述最优流道构型的流道结构参数的取值范围，具体包括：宽度W的取值范围、长度L的取值范围和深度D的取值范围；
    6)在步骤5)得到的取值范围内，分别在宽度W的取值范围、长度L的取值范围和深度D的取值范围中，取相同个数的代表值，且各代表值按均等间距选取；
    使用最优流道构型对应的真实粗糙度条件下的模拟模型，在不同进口压力H条件下，分别模拟各代表值从小到大依次组合所代表的滴灌灌水器出流流量Q，得到公式(2)：Q＝kHx，公式(2)中的Q为出流流量，k为流量系数，H为进口压力，x为流态指数；
    得到各组合所代表的滴灌灌水器的流量系数k，将各组合所代表的流道结构参数与流量系数k进行拟合，得到滴灌灌水器流量系数k的预报模型；
    7)将步骤6)各组合或按照步骤6)的方法设置的各组合按照流道结构参数从小到大的顺序排列，计算各组合抗堵塞评估指数P值，再计算相邻两个组合的P值间的偏差，从相邻两个组合的P值间的偏 差小于s％的流道结构参数较小的若干组合中，选取一定数量的组合，将其结构参数的值分别代入步骤6)滴灌灌水器流量系数k的预报模型中，计算流量系数k，将k值最小的组合的流道结构参数确定为滴灌灌水器最优流道结构参数。
  4. 如权利要求1所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：还包括步骤4：
    对步骤3所述最优流道结构原形进行二级精细结构设计，将所述最优流道结构原形的流道中内壁近壁面流动剪切力处在0.2-0.4Pa间的内壁位置优化为内壁近壁面流动剪切力不在0.2-0.4Pa间。
  5. 如权利要求4所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：步骤4所述优化使用漩涡洗壁优化设计方法：
    分析步骤3的所述最优流道结构原形流道内漩涡分布，根据漩涡外边缘的形状和大小，将流道中内壁近壁面流动剪切力处在0.2-0.4Pa间的内壁位置设计为与漩涡外边缘形状和大小相近或相同的圆弧。
  6. 如权利要求4或5所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：所述流道中内壁近壁面流动剪切力，是使用最优流道结构原形对应的真实粗糙度条件下的模拟模型或最优流道构型的真实粗糙度条件下的模拟模型，对流道内部的水流及颗粒物运动进行固-液-气三相流动模拟得到的。
  7. 如权利要求2-5中任一所述的面向流量设计需求的滴灌灌水器流道结构设计方法，其特征在于：所述初始粗糙度条件下的模拟模型为同时使用RNG k-ε模型和VOF模型表征滴灌灌水器流道构型内部的水流及颗粒物运动固-液-气三相流动湍流模拟模型；
    和/或，初始粗糙度条件下的模拟模型中的边壁粗糙度均值以经验值作为默认值；
    和/或，所述内壁位置分为齿尖迎水区、齿跟迎水区、齿尖背水区和齿跟背水区。
  8. 权利要求1-7中任一所述的面向流量设计需求的滴灌灌水器流道结构设计方法在生产滴灌灌水器产品中的应用。
  9. 权利要求8所述应用中生产的滴灌灌水器产品。
  10. 如权利要求9所述的滴灌灌水器产品，其特征在于；所述滴灌灌水器的流道构型为分形流道，
    所述分形流道为修正的分形-M流道，
    所述分形流道的一端口与所述滴灌灌水器的进水口相通，另一端口与所述滴灌灌水器的出水口相通，
    当所述滴灌灌水器的类型为片式时，所述分形流道的长度为39.567mm，宽度为0.824mm，深度为0.7555mm，流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.412mm的圆弧；
    当所述滴灌灌水器的类型为圆柱式时，所述分形流道的长度为214.4mm，宽度为1.27mm，深度为0.745mm，流道内壁位置的齿尖迎水区和齿跟迎水区为半径为0.423mm的圆弧。
PCT/CN2016/000339 2016-06-24 2016-06-24 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品 WO2017147728A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US15/773,948 US10750683B2 (en) 2016-06-24 2016-06-24 Drip irrigation emitter flow channel structural design method and fractal flow channel drip irrigation emitter product therefor
PCT/CN2016/000339 WO2017147728A1 (zh) 2016-06-24 2016-06-24 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2016/000339 WO2017147728A1 (zh) 2016-06-24 2016-06-24 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品

Publications (1)

Publication Number Publication Date
WO2017147728A1 true WO2017147728A1 (zh) 2017-09-08

Family

ID=59742352

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/000339 WO2017147728A1 (zh) 2016-06-24 2016-06-24 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品

Country Status (2)

Country Link
US (1) US10750683B2 (zh)
WO (1) WO2017147728A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109542426A (zh) * 2018-11-30 2019-03-29 广州市网拓信息技术有限公司 一种图像传感行为仿真模型用生成软件及方法
CN110583433A (zh) * 2019-09-27 2019-12-20 夏志强 一种用于农业的防堵型节水灌溉装置
CN112818447A (zh) * 2021-01-29 2021-05-18 中国地质大学（武汉） 一种城市空间结构分析方法
CN113455344A (zh) * 2021-07-27 2021-10-01 山东大学 一种星形迷宫流道滴灌灌水器

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3485215B1 (en) 2016-07-12 2023-06-07 Alexander Poltorak System and method for maintaining efficiency of a heat sink

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454459A (zh) * 2002-04-29 2003-11-12 中国农业大学 一种滴灌灌水器
CN101021881A (zh) * 2007-02-09 2007-08-22 中国农业大学 一种抗堵塞滴灌灌水器设计方法
CN101551828A (zh) * 2009-03-25 2009-10-07 西北农林科技大学 一种迷宫流道结构灌水器的抗堵塞流道结构设计方法
CN101667218A (zh) * 2009-10-22 2010-03-10 北京市水利科学研究所 削除流道内悬浮颗粒集聚位置的灌水器抗堵塞设计方法
CN101694679A (zh) * 2009-10-21 2010-04-14 西北农林科技大学 一种抗堵塞滴灌灌水器流道结构设计方法
CN102609569A (zh) * 2012-01-16 2012-07-25 中国水利水电科学研究院 一种低压滴灌锯齿型灌水器流道结构优化设计方法
US20150106237A1 (en) * 2013-10-11 2015-04-16 Steele Lorenz System for design of drip irrigation system
CN106096179A (zh) * 2016-06-24 2016-11-09 中国农业大学 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202958354U (zh) 2012-12-06 2013-06-05 中国水利水电科学研究院 一种低压滴灌专用锯齿型灌水器

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454459A (zh) * 2002-04-29 2003-11-12 中国农业大学 一种滴灌灌水器
CN101021881A (zh) * 2007-02-09 2007-08-22 中国农业大学 一种抗堵塞滴灌灌水器设计方法
CN101551828A (zh) * 2009-03-25 2009-10-07 西北农林科技大学 一种迷宫流道结构灌水器的抗堵塞流道结构设计方法
CN101694679A (zh) * 2009-10-21 2010-04-14 西北农林科技大学 一种抗堵塞滴灌灌水器流道结构设计方法
CN101667218A (zh) * 2009-10-22 2010-03-10 北京市水利科学研究所 削除流道内悬浮颗粒集聚位置的灌水器抗堵塞设计方法
CN102609569A (zh) * 2012-01-16 2012-07-25 中国水利水电科学研究院 一种低压滴灌锯齿型灌水器流道结构优化设计方法
US20150106237A1 (en) * 2013-10-11 2015-04-16 Steele Lorenz System for design of drip irrigation system
CN106096179A (zh) * 2016-06-24 2016-11-09 中国农业大学 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109542426A (zh) * 2018-11-30 2019-03-29 广州市网拓信息技术有限公司 一种图像传感行为仿真模型用生成软件及方法
CN109542426B (zh) * 2018-11-30 2022-06-21 广州市网拓信息技术有限公司 一种图像传感行为仿真模型生成系统及方法
CN110583433A (zh) * 2019-09-27 2019-12-20 夏志强 一种用于农业的防堵型节水灌溉装置
CN112818447A (zh) * 2021-01-29 2021-05-18 中国地质大学（武汉） 一种城市空间结构分析方法
CN113455344A (zh) * 2021-07-27 2021-10-01 山东大学 一种星形迷宫流道滴灌灌水器
CN113455344B (zh) * 2021-07-27 2022-07-01 山东大学 一种星形迷宫流道滴灌灌水器

Also Published As

Publication number Publication date
US10750683B2 (en) 2020-08-25
US20190110413A1 (en) 2019-04-18

Similar Documents

Publication Publication Date Title
WO2017147728A1 (zh) 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品
CN106096179B (zh) 滴灌灌水器流道结构设计方法及其分形流道滴灌灌水器产品
CN105677964B (zh) 一种基于cfd仿真和网格自适应的阀门流量系数计算方法
CN101551828B (zh) 一种迷宫流道结构灌水器的抗堵塞流道结构设计方法
CN105781626B (zh) 一种大子午扩张涡轮的三维正交非定常设计方法
CN110534776B (zh) 一种基于fluent仿真的燃料电池多堆分配内腔结构优化的方法
CN204612895U (zh) 一种小型流场模拟水槽
CN105550481B (zh) 一种基于水重力热管的烟气换热器的优化设计方法
CN102630426B (zh) 文丘里施肥器的结构优化方法
CN112818495A (zh) 一种新型管道压降测算算法参数的动态修正方法
CN105045987B (zh) 孔板厚度对孔板能量损失系数的影响关系计算方法
CN110196987B (zh) 基于代理模型的风道结构尺寸优化方法
CN108256262A (zh) 一种轴对称射流稳压腔参数设计的数值模拟方法
CN106874539B (zh) 用于叠片过滤器的滤片流道结构设计的方法及装置
CN100567082C (zh) 一种用于构造进气道斜切进口的方法
CN108446452A (zh) 一种混流泵叶轮鲁棒优化设计方法
CN102609569A (zh) 一种低压滴灌锯齿型灌水器流道结构优化设计方法
Qingsong et al. Study of hydraulic performance of the eddy channel for drip emitters
CN106777589A (zh) 一种基于流量需求的锯齿型灌水器逆向设计方法
CN111177945A (zh) 一种离心式过滤器的优化设计方法
CN104537170A (zh) 用于泵站三维流场计算的网格划分方法及装置
CN103077317B (zh) 一种用于稳定激波的泄漏槽的流量计算方法
CN109141542A (zh) 基于临界水深槽工作原理的矩形渠道无槽量水方法
Guo et al. Evaluation of numerical simulation accuracy for two-ways mixed flow drip irrigation emitter based on CFD
Koren Euler flow solutions for a transonic windtunnel section

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16891921

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16891921

Country of ref document: EP

Kind code of ref document: A1