CN114184239A - Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove - Google Patents

Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove Download PDF

Info

Publication number
CN114184239A
CN114184239A CN202111485577.2A CN202111485577A CN114184239A CN 114184239 A CN114184239 A CN 114184239A CN 202111485577 A CN202111485577 A CN 202111485577A CN 114184239 A CN114184239 A CN 114184239A
Authority
CN
China
Prior art keywords
rectangular groove
shrinkage
flow
optimal
central cylindrical
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
CN202111485577.2A
Other languages
Chinese (zh)
Inventor
朱晓明
姚莉
李金波
顾信钦
陈浩
刘志峰
周意恒
张思金
石先罗
任长江
斯静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Institute of Technology
Original Assignee
Nanchang Institute of Technology
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 Nanchang Institute of Technology filed Critical Nanchang Institute of Technology
Priority to CN202111485577.2A priority Critical patent/CN114184239A/en
Publication of CN114184239A publication Critical patent/CN114184239A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/002—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow wherein the flow is in an open channel
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • 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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于最优收缩度中心圆柱矩形槽的渠道测流法，主要包括通过CFD方法对给定中心圆柱矩形槽在不同收缩度条件下进行模拟，通过水跃的形态来判断柱前是否存在有效水深，当某一最小收缩度在任意流量条件下均存在有效水深时即得到最优收缩度；借助SPSS和Origin对不同收缩度及流量条件下的水深和收缩过流宽度进行数据处理，得到基于不同收缩度的中心圆柱矩形槽的流量预测公式，将最优收缩度带入即可得到基于最优收缩度中心圆柱矩形槽的流量预测公式，在渠道测流时依据渠道尺寸将最优收缩度的中心圆柱矩形槽安置渠道中，待水流稳定后测得水深带入流量公式即可得到实时的渠道流量。本发明能够准确地测试渠道流量，操作便捷，制作廉价，适应性强。

Figure 202111485577

The invention discloses a channel flow measurement method based on a central cylindrical rectangular groove with an optimal shrinkage degree. Whether there is an effective water depth before, when a certain minimum shrinkage has an effective water depth under any flow conditions, the optimal shrinkage is obtained; with the help of SPSS and Origin, the water depth and shrinkage overflow width under different shrinkage and flow conditions are analyzed. Processing, the flow prediction formula of the central cylindrical rectangular groove based on different shrinkage degrees is obtained, and the optimal shrinkage degree is brought into the flow prediction formula of the central cylindrical rectangular groove based on the optimal shrinkage degree. The central cylindrical rectangular groove with the optimal shrinkage is placed in the channel. After the water flow is stable, the measured water depth is brought into the flow formula to obtain the real-time channel flow. The invention can accurately test the channel flow, is convenient to operate, inexpensive to manufacture and has strong adaptability.

Figure 202111485577

Description

Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove
Technical Field
The invention relates to a channel flow measurement method based on an optimal shrinkage center cylindrical rectangular groove, and belongs to the technical field of hydraulic engineering.
Background
The water management of the irrigation area is crucial to the construction of the irrigation area, and mainly refers to the reasonable adjustment and distribution of irrigation water quantity, flow and irrigation time, so that the benefits of irrigation and drainage engineering are brought into play, irrigation water sources are effectively utilized, the irrigation quality and efficiency are improved, and the stable and high yield of agriculture is promoted. The current commonly used method for water measurement in irrigation areas comprises the following steps: using water in canal system building; water is measured by using a special water measuring facility; water is measured with a current meter, etc. The method has high measurement requirement, the building cannot be damaged, deformed, leaked water, silted, blocked and the like, and internal test equipment is difficult to calibrate and maintain, and the test method is complicated and is not beneficial to application. In addition, the first task of using canal system buildings to measure water in engineering is to construct corresponding buildings, which requires a lot of manpower, material resources and financial resources. (II) the water is supplied by utilizing a special water supply setting, and the method mainly comprises four modes: firstly, a U-shaped channel flat-bottom parabola throat-section-free rectangular groove is used for measuring water quantity, the method has certain defects, the accuracy of flow prediction is poor, the scale effect is not considered during coefficient check, the flow coefficient checked on the model scale is applied to the prototype scale and is necessarily inaccurate, the speed coefficient has great experience, the application range is small, and the method can only be used for a flat-bottom channel and cannot be moved once being built; secondly, the Parshall rectangular groove is also a water measuring facility formed by contraction of open channels, and has harsh flow measuring conditions, high design and installation requirements and large measuring error; thirdly, the triangular water measuring weir has the advantages of simple structure, low manufacturing cost, convenient test and relatively high precision, but has small water passing capacity, complex calculation and poor universality; the trapezoidal measuring weir has the advantages that the structure is simple, the manufacturing cost is low, the test is convenient, when the water level in front of the weir is high, silting and congestion can be caused in a channel with poor water quality, the channel can only be used for flow measurement in clear water, the technical requirements on users during installation and use are high, and the popularization is poor; the current measuring method by using the current meter has the disadvantages of complicated current measuring process, long time consumption, poor measuring precision and difficult popularization and application in irrigation areas. With the acceleration of the modernization process, ultrasonic flow meters, electromagnetic flow meters and the like can also be applied to irrigation areas, but the ultrasonic flow meters, the electromagnetic flow meters and the like are expensive in manufacturing cost and high in maintenance cost, and the wide application of the ultrasonic flow meters, the electromagnetic flow meters and the like in the irrigation areas is limited. It becomes a difficult problem how to simply and accurately test the channel flow of the irrigation area.
Scholars at home and abroad set obstacles in the rectangular grooves to form critical flows in the rectangular grooves, wherein the critical flows comprise trihedrons, cylinders, semicylinders and the like, a flow measurement theoretical formula is provided through dimension analysis or column writing energy equations, and the flow measurement can be achieved by measuring corresponding parameters. The method for measuring the flow through the rectangular groove of the central cylinder is a simple and convenient flow measuring method with high precision, has strong applicability, and is more suitable for measuring the flow of all levels of channels in irrigation areasMore importantly, compared with obstacles such as three-surface-shaped obstacles, semi-cylinders on two sides and the like, the central cylindrical rectangular groove is simple in structure and simple in theoretical formula, but the relation coefficient of the flow and the water depth in front of the column needs to be determined in advance, and the relation between the flow and the water depth in front of the column can be obtained through dimensional analysis: q ═ a (B)c 2.5)(g0.5)(h/Bc)bAt present, students usually determine coefficients a and b through a physical model test, but the workload of the model test is large, the limiting conditions are many, errors hardly meet the precision requirement, and more importantly, when the coefficients are calibrated, the researchers usually think that the coefficients a and b obtained through 2-3 shrinkage degrees can be suitable for predicting the flow rate of a rectangular groove with any shrinkage degree, which is obviously unreasonable, and the popularization and the application of the flow measuring method of the central cylindrical rectangular groove are seriously restricted. In addition, the relation between the flow and the water depth in front of the column is usually determined by adopting a model test method at present, the calibration method has large workload, a plurality of limiting conditions and errors, the precision requirement is difficult to meet, and the influence of the shrinkage is difficult to consider in the derivation of a theoretical formula, so that the rectangular groove flow measurement method of the central cylinder is difficult to popularize and apply. In addition, the shrinkage of the rectangular groove of the central cylinder also influences the diameter of the central cylinder, which determines the amount of consumable materials of the central cylinder, and if a minimum and universal shrinkage (optimal shrinkage) is determined, the energy conservation and environmental protection are facilitated. With the rapid development of fluid dynamics (CFD) and cloud computing technologies, numerical simulation becomes another effective research means due to the advantages of convenient modeling, direct simulation of a prototype, good test repeatability and the like. However, no method for predicting the flow of the central cylindrical rectangular groove by using a numerical test means is available in the prior art. Chinese patent CN201910015788.6 discloses a method for predicting critical submergence depth of a wet room type pump station by CFD, which uses CFD to predict critical submergence depth of a pump station, but the method cannot be used to predict flow of a central cylindrical rectangular groove, and cn201811325431.x discloses a method for predicting flow of a large low-lift water pump, which uses CFD in combination with a differential pressure flow measurement method to predict flow of a large low-lift water pump, and also the method cannot be used for channel flow test of the present application. Thus, how to utilize the cloudThe flow measurement precision of the central cylindrical rectangular groove is improved by a calculation and CFD method, the influence of the shrinkage on the flow measurement precision is fully considered, and a determination method of the optimal shrinkage is provided on the basis, so that the irrigation district channel flow can be accurately and economically tested, and the method is a very worthy of research, and the construction process of realizing intelligent irrigation in the irrigation district is promoted.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects in the prior art, the invention provides a channel flow measuring method based on an optimal shrinkage degree central cylindrical rectangular groove, which can predict the flow of the central cylindrical rectangular groove by using a numerical test means, can accurately and economically test the channel flow of an irrigation area under the condition of fully considering the shrinkage degree and the influence of the optimal shrinkage degree on the flow measuring precision, effectively avoids the waste of irrigation water sources, and improves the irrigation quality and efficiency.
The technical scheme is as follows: in order to solve the technical problem, the channel flow measurement method based on the optimal shrinkage degree central cylindrical rectangular groove comprises the following steps of:
step 1: establishing a physical test model of the central cylindrical rectangular groove and carrying out a physical model test, collecting the height of the submerged water depth at the front end of the cylinder after the water flow is stable, and determining the optimal shrinkage of the central cylindrical rectangular groove by a CFD (computational fluid dynamics) method;
in step 1, the determination of the optimal shrinkage r comprises the following steps:
step 1-1: carrying out a physical model test on the central cylindrical rectangular groove, manufacturing a central cylindrical rectangular groove model for fixing a certain shrinkage degree, wherein the material is transparent organic glass, the bottom of the rectangular groove is horizontally placed, the overflowing section of the rectangular groove is rectangular, the bottom width is B, the height is M, the diameter of the central cylinder is D, the ratio of the diameter of the central cylinder to the width of the rectangular groove is defined as the shrinkage degree r, and the test is carried out under different flow conditions to test the water depth in front of the central cylindrical groove;
step 1-2: determining an optimal numerical simulation method, modeling and performing numerical simulation calculation on a central cylindrical rectangular groove under a given shrinkage condition in a physical model test, meshing the established model by adopting the sizes of meshes with different scales, importing an output calculation file into Fluent, performing numerical simulation by adopting different turbulence models, dispersion methods, numerical algorithms and time step lengths, respectively deriving calculation results, performing Post-processing through CFD-Post and collecting data, finding out a calculation condition close to the test water depth of the central cylindrical rectangular groove with the given shrinkage in the step 1 from a plurality of numerical calculation data, and outputting the corresponding mesh size, turbulence models, dispersion methods, numerical algorithms and time step lengths;
step 1-3: determining the optimal shrinkage of the central cylindrical rectangular groove, keeping the width of the rectangular groove unchanged, determining the shrinkage of the central cylindrical rectangular groove by changing the diameter of the central cylinder, respectively performing geometric modeling on the central cylindrical rectangular groove under different shrinkage r conditions, and performing geometric modeling on the central cylindrical rectangular groove under different dimensionless numbers F0＝QB-1M-1.5g-0.5Under the condition, the optimal numerical simulation method obtained in the step 1-2 is adopted to carry out CFD calculation, and the dimensionless number F0The value interval is 0.01, the bottom width of the rectangular groove is kept the same during modeling, different shrinkage degrees are determined by changing the diameter of the central cylinder, the value range of the shrinkage degree r is more than 0.0 and less than 1.0, and the value interval of r is 0.01, namely 0.01, 0.02, 0.03, … …, 0.97, 0.98 and 0.99 when the geometric model is established; the order of calculation of the degree of shrinkage r is chosen as: when the calculation sequence of the shrinkage r is determined by adopting a center mean value method, namely, the flow field of the central cylindrical rectangular groove when r is 0.99 is calculated, and all dimensionless numbers F when r is 0.99 are calculated0The water jump shape of (1) is that the rectangular groove has no optimal shrinkage degree when the water jump on the column is carried out, and when r is all dimensionless numbers F of 0.990When the hydraulic jump is normal (the hydraulic jump does not occur on the column), the CFD calculation is performed by re-determining the shrinkage degree downward at a value of r ═ 0.01+0.99)/2 ═ 0.50; when r is 0.50, there is a dimensionless number F0When the hydraulic jump is on-column hydraulic jump, the upward redetermined shrinkage degree is r ═ 0.50+0.99)/2 ═ 0.745, rounded to 0.75, and when r ═ 0.50, all dimensionless numbers F0When the hydraulic jump is normal, the shrinkage is re-determined downward at a value of r ═ 0.50+0.01)/2 ═ 0.255, rounded to 0.26, and the water jump is thus characterizedCalculating the flow field of the central cylindrical rectangular groove with different shrinkage degrees until the water jump of the central cylindrical rectangular groove with a certain shrinkage degree r is a normal water jump and the water jump of the central cylindrical rectangular groove with the shrinkage degree smaller than the normal shrinkage degree by 0.01 is an on-column water jump, stopping calculating and outputting the shrinkage degree, wherein the shrinkage degree can be determined as the optimal shrinkage degree of the target rectangular groove and is recorded as r0Observing the form of the front hydraulic jump of the central cylinder in the water tank under the conditions of any shrinkage and any flow, wherein the judgment standard of the hydraulic jump form is as follows: when the hydraulic jump is on-column hydraulic jump, the shrinkage of the central cylindrical rectangular groove is ineffective shrinkage, and when the hydraulic jump is not on-column hydraulic jump (normal hydraulic jump), the shrinkage is effective shrinkage of the central cylindrical rectangular groove;
step 2: for the optimal shrinkage r which is more than or equal to the output of the step 1-30(r≥r0) 6 shrinkage degrees, namely r1, r2, r3, r4, r5 and r6 are selected, and r 1-r 6 are uniformly distributed on r0R is less than or equal to 1, the central cylindrical rectangular groove with the selected shrinkage degree is subjected to numerical calculation under different flow conditions, corresponding calculated flow Q, water depth h at the front end of the cylinder and flow width Bc after the rectangular groove is shrunk are output, and SPSS software is adopted to carry out numerical calculation on submergence water depth h at the front end of the cylinder and flow width B after the flow passage is shrunk under the conditionscAnd the corresponding calculated flow Q is combined with a theoretical formula Q ═ a (B)c 2.5)(g0.5)(h/Bc)bCarrying out nonlinear processing, respectively outputting coefficients a and b under different shrinkage r conditions, then adopting Origin to carry out fitting analysis on the coefficients a, b and r and outputting a functional relation of the coefficients, wherein g is gravity acceleration;
and step 3: substituting the fitting relational expression of the coefficients a, b and r into a theoretical formula to obtain an expression of real-time flow in the rectangular groove considering the influence of shrinkage;
and 4, step 4: will optimize the degree of shrinkage r0And corresponding contracted rectangular groove overflowing width Bc0Substituting the expression of the real-time flow in the rectangular groove considering the influence of the shrinkage in the step 3 to obtain the expression of the real-time flow in the central cylindrical rectangular groove based on the optimal shrinkage;
and 5: surveying the channel needing flow measurement in field, and measuring the channel bottom width WAmplifying the central cylindrical rectangular groove with the optimal shrinkage in the step 1 to obtain a prototype central cylindrical rectangular groove with the optimal shrinkage, wherein the amplification factor is lambda0Arranging a central cylinder rectangular groove in the channel, and measuring the depth h of the submerged water at the front end of the central cylinder after the water flow is stable0And a shrinkage width W-Wr0And (4) substituting the expression based on the real-time flow in the central cylindrical rectangular groove with the optimal shrinkage obtained in the step (4) to obtain the real-time channel flow.
Further, in the step 1-1, the length L of the rectangular groove is not less than 10D, and the length L of the central cylinder from the overflow section to the inlet section is the diameter perpendicular to the positive direction of the incoming flowinNot less than 7D, length L from outlet sectionoutNot less than 3D, wherein D is central cylinder diameter, and central cylinder overlaps with rectangular channel central symmetry face along the cross section that rivers positive direction diameter place.
Further, in the step 1-1, model tests are performed on the central cylindrical rectangular groove under different flow conditions, the flow value is selected from small to large, the water depth in the rectangular groove is gradually deepened along with the increase of the flow, the test is stopped when the water depth at the inlet of the rectangular groove is lower than the height of the rectangular groove by 0.2D, D is the diameter of the cylinder, the flow Qt in the rectangular groove is tested through an electromagnetic flowmeter in the test process, the submergence depth ht at the front end of the corresponding central cylinder is read through a water ruler arranged at the front end of the cylinder, subscript t of Qt and ht represents the test conditions, the unit of Q is L/s, and the unit of h is cm.
Further, hexahedral structured grids of different scales are adopted in grid division in the step 1-2, an external expansion structure is adopted in central cylindrical grid division, so that a plurality of different grid schemes are obtained, a VOF method is adopted in Fluent calculation for each grid scheme, different turbulence models and different numerical algorithms are selected for the VOF method, flow adopted in numerical calculation corresponds to output flow tested in the step 1 one to one, boundary conditions adopted in calculation are that an inlet is a speed inlet, an outlet is a pressure outlet, a surface is a pressure inlet, a solid boundary is set to wall, an initial condition is that air proportion is 100%, water depth hs obtained through numerical simulation is water depth hs under the condition that water and air in Post-processing CFD-Post respectively occupy 50% of free liquid level, the unit is cm, and subscript s is a digital-analog condition.
Further, error analysis is carried out on the water depth ht under the test condition obtained in the step 1-2 and the water depth hs under the digital-analog condition, the calculation formula is that | ht-hs | and/ht are less than or equal to k%, wherein | | | represents an absolute value, and when all the k values under a certain numerical calculation scheme are not more than 7, a grid division scheme and a Fluent setting scheme are output.
Further, in step 3, the relationship between the coefficient a and r is expressed by a quadratic term a ═ a1r2+B1r+C1The relation between the coefficient b and r is expressed in the form of a quadratic term b ═ a2r2+B2r+C2Wherein A is1，A2，B1，B2，C1，C2Are column term coefficients.
Further, in the step 4, based on the expression of real-time flow in the central cylindrical rectangular slot with the optimal degree of contraction
Figure BDA0003397407160000051
r0And Bc0The flow width of the rectangular groove corresponding to the optimal shrinkage degree and the optimal shrinkage degree is obtained.
Further, in the step 5, the real-time channel flow obtained based on the expression of the real-time flow in the central cylindrical rectangular groove with the optimal contraction degree is
Figure BDA0003397407160000052
Furthermore, cloud computing technology is introduced into both the model and prototype CFD computing in the invention, so that the numerical computing precision and the computing efficiency can be effectively improved.
The invention tests the real-time flow in irrigation district channels by a channel flow measuring method based on a central cylindrical rectangular groove with optimal shrinkage, which comprises the steps of verifying a numerical simulation method, carrying out numerical calculation and physical tests on a rectangular groove model, and analyzing and comparing the submergence depth at the front end of the cylinder to obtain a reasonable CFD numerical scheme comprising a mesh subdivision mode, a turbulence model and a numerical algorithm; then simulating a plurality of shrinkage degrees (fixing the width of the water tank and changing the diameter of the central cylinder) by obtaining a CFD numerical scheme, judging the pre-column water depth when the water jump on the column occurs in the central cylindrical rectangular groove as an invalid water depth, judging the pre-column water depth when the water jump on the column does not occur in the central cylindrical rectangular groove as an effective water depth, and outputting the minimum shrinkage degree when the effective water depth exists under all flow working conditions as the optimal shrinkage degree; then, the SPSS carries out nonlinear analysis on the water depth at the front end of the cylinder, the width of the corresponding rectangular groove after contraction and the calculated flow under the condition of different degrees of contraction which are output and are more than or equal to the optimal degree of contraction, and is combined with a theoretical formula to obtain a relational expression of a variable coefficient and the degree of contraction in the theoretical formula, and then the variable coefficient expressed by the degree of contraction is brought back to the theoretical formula to obtain a flow prediction formula of the rectangular groove with the center cylinder under the condition of any degree of contraction which is more than or equal to the optimal degree of contraction; and finally, substituting the obtained optimal shrinkage and the shrinkage width corresponding to the optimal shrinkage into a flow prediction formula under any shrinkage condition with the optimal shrinkage being more than or equal to the optimal shrinkage to obtain an expression based on the flow of the central cylindrical rectangular groove with the optimal shrinkage. According to the invention, the accuracy of CFD numerical simulation is improved by comparing a physical model test with a corresponding model CFD numerical calculation and using the model test to constrain a numerical calculation method, prototype data are predicted by calculating a prototype rectangular groove by an obtained high-precision numerical method, the influence of the shrinkage on the precision of a prediction result is considered, and how to reduce the manufacturing cost of a central cylinder to the minimum is further considered, so that the optimal shrinkage of any central cylinder rectangular groove is provided, and finally the real-time flow of a channel under the condition of the optimal shrinkage in an irrigation area can be accurately predicted.
Has the advantages that: the channel flow measurement method based on the central cylindrical rectangular groove with the optimal contraction degree can simply, accurately and economically predict the real-time flow of the channel in the irrigation area, so that the utilization efficiency of irrigation water is improved, the phenomenon that return water is wasted due to overlarge flow or irrigation is insufficient due to undersize flow is avoided, the optimal contraction degree of the central cylindrical rectangular groove is provided, the manufacturing consumables of the central cylinder are effectively reduced, and the manufacturing cost is saved.
Drawings
FIG. 1 is a flow chart of a channel flow measurement method based on an optimal shrinkage degree central cylindrical rectangular groove;
FIG. 2 is a schematic diagram of a central cylindrical rectangular slot model;
FIG. 3 is a split view of a central cylindrical rectangular slot grid;
FIG. 4 is a diagram of verification of numerical simulation accuracy of a central cylindrical rectangular groove;
FIG. 5 is a diagram of a central cylindrical rectangular channel numerically simulated free liquid level;
FIG. 6 is a schematic view of a normal hydraulic jump of a central cylindrical rectangular tank;
FIG. 7 is a schematic view of a hydraulic jump on a central cylindrical rectangular groove column;
FIG. 8 is a diagram illustrating the calculation of the optimal shrinkage;
FIG. 9 is a drawing showing data acquisition of the depth of water submerged at the front end of a cylinder in a central cylindrical rectangular groove with various degrees of contraction;
FIG. 10 is a graph of a quadratic polynomial fit of r to coefficient a;
FIG. 11 is a graph of a quadratic polynomial fit of r to coefficient b.
Detailed Description
The present invention will be further described with reference to the accompanying drawings.
A flow chart of the channel flow measurement method based on the optimal shrinkage degree central cylindrical rectangular groove in this embodiment is shown in fig. 1, and the specific steps are as follows:
step 1: establishing a physical test model of the central cylindrical rectangular groove and carrying out a physical model test, collecting the height of the submerged water depth at the front end of the cylinder after the water flow is stable, and determining the optimal shrinkage of the central cylindrical rectangular groove by a CFD (computational fluid dynamics) method;
in the physical model in the step 1-1, the central cylindrical rectangular groove model is made of transparent organic glass, the bottom of the rectangular groove is horizontally placed, the flow cross section of the rectangular groove is rectangular, the bottom width is 17cm, the height is 30cm, the diameter of the central cylinder is 11.4cm, and the corresponding shrinkage ratio is 0.670.
The length of the rectangular groove is 300cm, the length of the overflow section where the diameter of the central cylinder perpendicular to the positive direction of the incoming flow is located from the inlet section is 200cm, the length of the overflow section from the outlet section is 100cm, and the section where the diameter of the central cylinder along the positive direction of the incoming flow is located is overlapped with the central symmetry plane of the rectangular groove.
In the test process, the flow rates are 7.026L/s, 8.080L/s, 9.016L/s, 10.070L/s, 11.241L/s, 11.475L/s, 12.178L/s, 14.520L/s, 20.258L/s and 22.482L/s from small to large respectively. When the shrinkage degree is measured by test, the water depth is 15.785cm, 16.861cm, 17.523cm, 18.441cm, 18.578cm, 19.254cm, 19.844cm, 20.920cm, 25.020cm, 25.439cm and 26.961 cm.
Step 1-2: carrying out numerical simulation calculation on a central cylindrical rectangular groove of a physical model test, establishing a geometric model of the central cylindrical rectangular groove as shown in figure 2, dividing grids as shown in figure 3, outputting a calculation file, importing Fluent to carry out numerical calculation, exporting a calculation result, carrying out Post-processing through CFD-Post, and collecting data;
the method adopts a VOF method in Fluent calculation, adopts boundary conditions that an inlet is a speed inlet, an outlet is a pressure outlet, a surface is a pressure inlet, a solid boundary is set to be wall, and the initial condition is that the air proportion is 100 percent;
the numerical simulation flow rate and the model test flow rate are same and respectively 7.026L/s, 8.080L/s, 9.016L/s, 10.070L/s, 11.241L/s, 11.475L/s, 12.178L/s, 14.520L/s, 20.258L/s and 22.482L/s;
finding out a calculation working condition close to the water depth tested in the step 1-1 from a plurality of numerical calculation results, and outputting a corresponding grid size, a turbulence model and a numerical algorithm;
when the water depth and height obtained by numerical simulation are the water depth under the condition that water and air respectively account for 50% of free liquid level in the Post-treatment CFD-Post, and the water depth obtained by simulation is 14.873cm, 15.893cm, 17.145cm, 18.417cm, 19.244cm, 19.702cm, 20.143cm, 21.275cm, 25.607cm, 26.098cm and 27.678cm, the relative error between the water depth obtained by numerical simulation and the water depth obtained by model test is not more than 7% (as shown in figure 4), the grid subdivision at the moment is output, the grid is 1cm multiplied by 1cm, the turbulence model is Realicablk-epsilon, and the discrete mode of the control equation is finite volume method FVM; the diffusion term adopts a second-order central difference format, the convection term adopts a Quick format, the coupling of pressure and speed adopts a SIMPLEC algorithm, and the calculation mode adopts parallel calculation.
And a specific numerical calculation step:
1. modeling a central cylindrical rectangular groove of the model by using three-dimensional modeling software SolidWorks (as shown in FIG. 2), and then outputting a suffix IGS file;
2. importing the suffix name IGS file into mesh subdivision software Ansys ICEM, carrying out mesh subdivision on a central cylindrical rectangular groove shown in figure 2 (shown in figure 3), wherein the mesh size is 1cm multiplied by 1cm, the mesh type is hexahedron structured mesh, and then generating an Ansys Fluent numerical value to calculate the suffix name of the file;
3. importing the calculation file into an Ansys Fluent to calculate, wherein the specific setting mode is as follows:
reducing x, y and z of Scale in General by 0.001 times in a same ratio, and enabling the modeling and the numerical calculation to be consistent in geometric dimension;
opening multiphase flow Volume of Fluid setting in Models, selecting Implicit calculation Implicit, and taking 1 × 10 for Curian number-8Considering the volume Force Implicit Body Force, selecting Realizblek-epsilon by a turbulence model;
adding a medium water into Matricals, wherein the main term in the multiphase flow is the water, and the secondary term in the multiphase flow is air;
selecting an inlet as a speed inlet from Boundary Conditions, adopting a side wall condition for wall, adopting an outlet as a pressure outlet, and adopting a pressure inlet on the upper surface;
coupling of pressure and speed in Solution Methods adopts a Coupled algorithm, and the discrete mode of a control equation is FVM; the diffusion item adopts a second-order central difference format, and the convection item adopts a Quick format;
the Moinors sets the convergence condition of the calculation, and in the calculation, in order to improve the calculation precision, the consistency, x, y and z are all set to reach 1 multiplied by 10-6Stopping the calculation;
setting an Initialization method as Standard Initialization in Solution Initialization, and setting 100% of an initial item as air;
setting Time Step Size to be 0.01s, number Time Steps to be 200000, Max operations to be 20, automatically saving Step Size to be 100 and 1s and saving a Calculation result in Run Calculation;
storing the Fluent calculation result after the calculation is stopped, and importing CFD-Post processing, wherein the method comprises the following specific steps:
the Isosurface is set, the water volume Fraction is set to be 50%, the free liquid level (as shown in figure 5) can be obtained, and the water depth at the front end of the cylinder can be measured.
Step 1-3: determining the optimal shrinkage of the central cylindrical rectangular groove, keeping the width of the rectangular groove unchanged, determining the shrinkage of the central cylindrical rectangular groove by changing the diameter of the central cylinder, respectively performing geometric modeling on the central cylindrical rectangular groove under different shrinkage r conditions, and performing geometric modeling on the central cylindrical rectangular groove under different dimensionless numbers F0＝QB-1M-1.5g-0.5Carrying out CFD calculation by adopting the optimal numerical simulation method obtained in the step 1-2 under the condition, keeping the bottom width of the rectangular groove the same during modeling, determining different shrinkage degrees by changing the diameter of the central cylinder, wherein the value range of the shrinkage degree r is more than 0.0 and less than 1.0, and the value interval of r is 0.01, namely 0.01, 0.02, 0.03, … …, 0.97, 0.98 and 0.99 during establishing the geometric model; the order of calculation of the degree of shrinkage r is chosen as: when the calculation sequence of the shrinkage r is determined by adopting a center mean value method, namely, the flow field of the central cylindrical rectangular groove when r is 0.99 is calculated, and all dimensionless numbers F of r is 0.99 are calculated0The hydraulic jump (c) of (a) is a normal hydraulic jump (as shown in fig. 6), and CFD calculation is performed with the degree of contraction re-determined downward as r ═ 0.01+0.99)/2 ═ 0.50, and r ═ 0.50 for all dimensionless numbers F0The water jump form (2) is normal water jump, the shrinkage is re-determined downwards as r ═ 0.50+0.01)/2 ═ 0.255, the value of r ═ 0.26 is rounded off, and CFD calculation is carried out, when the dimensionless number F is0When the water jump is less than 0.173, the water jump is on-column water jump (as shown in fig. 7), CFD calculation is carried out when the upward redetermined shrinkage degree value is r ═ 0.26+0.50)/2 ═ 0.38, the water jump is normal water jump, CFD calculation is carried out when the downward redetermined shrinkage degree value is r ═ 0.38+0.26)/2 ═ 0.32, and when the dimensionless number F is smaller than0The water jump state is water jump on column when < 0.122, the upward redetermined shrinkage degree is r ═ 0.32+0.38)/2 ═ 0.35, and all dimensionless numbers F thereof are0Water jumpFor normal hydraulic jump, the downward redetermined shrinkage has a value of r ═ 0.35+0.32)/2 ═ 0.34, when the dimensionless number F is0If the water jump is on-column water jump < 0.106 and the calculation results of all the above shrinkage degrees are shown in FIG. 8, the optimal shrinkage degree of the central cylindrical rectangular groove is r0＝0.35。
Step 2: respectively carrying out geometric modeling on the rectangular groove with the center cylinder in the step 1-1 under the condition of different shrinkage degrees, carrying out CFD calculation on the split grid after the modeling is finished, carrying out Post-processing on the calculation result by adopting CFD-Post (the result is shown in figure 9), wherein the grid splitting mode, the Fluent setting method and the CFD-Post-processing method are the same as the numerical method output in the step 1-2, SPSS software is adopted to carry out submerged water depth at the front end of the output cylinder, the overflow width after the runner is shrunk and the corresponding calculated flow, and the result theoretical formula Q is a (B)c 2.5)(g0.5)(h/Bc)bAnd performing nonlinear processing, respectively outputting coefficients a and b under different gradient conditions, performing fitting analysis on the coefficients a, b and r by adopting Origin, and outputting a functional relationship of the coefficients, wherein the expression form of the relationship between the coefficients a and r is shown in figure 10 by quadratic fitting, and the expression form of the relationship between the coefficients b and r is shown in figure 11 by quadratic fitting.
And step 3: substituting the fitting relation of the coefficients a, B and r into a theoretical formula Q ═ a (B)c 2.5)(g0.5)(h/Bc)bThe expression of the real-time flow of the rectangular groove can be obtained as
Figure BDA0003397407160000091
And 4, step 4: substituting the optimal contraction degree of 0.35 into the step 3 to obtain an expression Q which is 0.606 (B) based on the real-time flow in the central cylindrical rectangular groove of the optimal contraction degreec 2.5)(g0.5)(h/Bc)1.571
And 5: measuring the bottom width of a channel in an irrigation area to be 1.7m, amplifying the bottom width of the channel by 10 times compared with the bottom width of the rectangular groove of the test model in the step 1, wherein the shrinkage r is the optimal shrinkage of 0.35, adaptively arranging the rectangular groove of the prototype center cylinder with the amplification of 10 times and the shrinkage of 0.35 in the channel, and measuring the right center cylinder after the water flow is stableThe front end submerged depth is 0.5m, the corresponding contracted rectangular groove overflowing width is 1.7 x (1-0.35) ═ 1.105m, and the real-time channel flow is 0.701m by substituting the parameters into a real-time flow expression3/s。
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.

Claims (8)

1.一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于，包括以下步骤：1. a channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove, is characterized in that, comprises the following steps: 步骤1：建立中心圆柱矩形槽的物理试验模型并进行物理模型试验，当水流流动稳定后采集圆柱正前端淹没水深的高度，并针对所述中心圆柱矩形槽通过CFD方法确定其最优收缩度；Step 1: Establish a physical test model of the central cylindrical rectangular groove and carry out the physical model test. When the water flow is stable, collect the height of the submerged water depth at the front end of the cylindrical groove, and determine its optimal shrinkage degree by CFD method for the central cylindrical rectangular groove; 所述步骤1中，最优收缩度r的确定包含以下步骤：In the step 1, the determination of the optimal shrinkage r includes the following steps: 步骤1-1：对中心圆柱矩形槽进行物理模型试验，制作某一固定收缩度的中心圆柱矩形槽模型，其材质采用透明有机玻璃，矩形槽的底部为水平放置，矩形槽过流断面为矩形，底宽为B，高度为M，中心圆柱直径为D，定义中心圆柱直径为D与矩形槽宽度B之比为收缩度r，在不同流量Q条件下进行试验，测试中心圆柱的柱前水深；Step 1-1: Carry out a physical model test on the central cylindrical rectangular groove, and make a central cylindrical rectangular groove model with a fixed shrinkage degree. The material is transparent plexiglass, the bottom of the rectangular groove is placed horizontally, and the flow section of the rectangular groove is rectangular , the bottom width is B, the height is M, the diameter of the central cylinder is D, and the ratio between the diameter of the central cylinder D and the width B of the rectangular groove is defined as the shrinkage r, and the test is carried out under different flow conditions Q, and the water depth in front of the column of the central cylinder is tested. ; 步骤1-2：确定最优的数值模拟方法，对物理模型试验中给定收缩度条件下的中心圆柱矩形槽进行建模和数值模拟计算，采用不同尺度的网格大小对建立的模型进行网格剖分，输出计算文件导入Fluent中，采用不同的湍流模型、离散方法、数值算法、时间步长进行数值模拟并分别导出计算结果，通过CFD-Post进行后处理并采集数据，从若干数值计算数据中找出与步骤1中给定收缩度中心圆柱矩形槽测试水深相近的计算工况，输出对应的网格尺寸、湍流模型、离散方法、数值算法、时间步长；Step 1-2: Determine the optimal numerical simulation method, carry out modeling and numerical simulation calculation of the central cylindrical rectangular groove under the condition of given shrinkage in the physical model test, and use the mesh size of different scales to mesh the established model. Lattice division, the output calculation file is imported into Fluent, different turbulence models, discrete methods, numerical algorithms, and time steps are used for numerical simulation and the calculation results are exported separately, and post-processing and data collection are performed through CFD-Post. Find out the calculation conditions that are similar to the test water depth of the given shrinkage center cylindrical rectangular groove in step 1 from the data, and output the corresponding grid size, turbulence model, discrete method, numerical algorithm, and time step; 步骤1-3：确定中心圆柱矩形槽的最优收缩度，保持矩形槽宽度不变，通过改变中心圆柱直径来确定中心圆柱矩形槽的收缩度，对不同收缩度r条件下的中心圆柱矩形槽分别进行几何建模，在不同无量纲数F0＝QB-1M-1.5g-0.5，条件下采用步骤1-2得到的最优数值模拟方法进行CFD计算，其中底宽为B，高度为M，g为重力加速度，Q为中心圆柱矩形槽模拟时给定的流量；无量纲数F0取值间隔为0.01，建模时保持矩形槽底宽相同，通过改变中心圆柱的直径来确定不同收缩度，收缩度r的取值范围为0.0＜r＜1.0，建立几何模型时r取值间隔为0.01，即0.01、0.02、0.03、……、0.97、0.98、0.99；收缩度r的计算顺序的选择为：当收缩度r的计算顺序采用中心均值法进行确定，即先计算r＝0.99时的中心圆柱矩形槽的流场，当r＝0.99的所有无量纲数F0的水跃形态为柱上水跃时该矩形槽无最优收缩度，当r＝0.99的所有无量纲数F0水跃形态为正常水跃时，向下重新确定收缩度的值为r＝(0.01+0.99)/2＝0.50进行CFD计算；当r＝0.50的存在无量纲数F0水跃形态为柱上水跃时，向上重新确定收缩度的值为r＝(0.50+0.99)/2＝0.745，四舍五入为0.75，当r＝0.50的所有无量纲数F0水跃形态为正常水跃时，向下重新确定收缩度的值为r＝(0.50+0.01)/2＝0.255，四舍五入为0.26，以此类推对不同收缩度的中心圆柱矩形槽流场进行计算，直至某收缩度r的中心圆柱矩形槽的水跃为正常水跃，而比该收缩度小0.01的中心圆柱矩形槽的水跃为柱上水跃时停止计算并输出该收缩度，该收缩度可以确定为目标矩形槽的最优收缩度记为r0，在任意收缩度和任意流量条件下观察水槽中中心圆柱前水跃的形态，水跃形态的判断标准为：当水跃形态为柱上水跃时的中心圆柱矩形槽收缩度为无效收缩度，当水跃形态为不为柱上水跃时的收缩度为中心圆柱矩形槽有效收缩度；Step 1-3: Determine the optimal shrinkage degree of the central cylindrical rectangular groove, keep the width of the rectangular groove unchanged, and determine the shrinkage degree of the central cylindrical rectangular groove by changing the diameter of the central cylinder. Geometric modeling is carried out respectively, and CFD calculation is carried out using the optimal numerical simulation method obtained in steps 1-2 under the condition of different dimensionless numbers F 0 =QB -1 M -1.5 g -0.5 , where the bottom width is B and the height is M, g are the gravitational acceleration, Q is the flow rate given in the simulation of the central cylindrical rectangular groove; the dimensionless number F 0 is set at an interval of 0.01, the bottom width of the rectangular groove is kept the same during modeling, and the difference is determined by changing the diameter of the central cylindrical Shrinkage, the value range of shrinkage r is 0.0 < r < 1.0, and the interval between r values when building a geometric model is 0.01, that is, 0.01, 0.02, 0.03, ..., 0.97, 0.98, 0.99; the calculation sequence of shrinkage r The choice is: when the calculation order of the shrinkage r is determined by the center mean method, that is, the flow field of the central cylindrical rectangular groove when r=0.99 is calculated first, and the hydraulic jump form of all dimensionless numbers F 0 when r=0.99 is When the hydraulic jump on the column has no optimal shrinkage degree, when all the dimensionless F 0 hydraulic jumps with r=0.99 are normal hydraulic jumps, the value of the shrinkage degree is re-determined downward as r=(0.01+0.99) /2=0.50 for CFD calculation; when r=0.50 exists dimensionless number F 0 hydraulic jump form is on-column hydraulic jump, the value of the upward re-determination of the contraction degree is r=(0.50+0.99)/2=0.745, rounded up is 0.75, when all the dimensionless numbers F 0 hydraulic jumps with r=0.50 are normal hydraulic jumps, the downward re-determination value of the contraction degree is r=(0.50+0.01)/2=0.255, which is rounded up to 0.26, so that By analogy, the flow field of the central cylindrical rectangular groove with different shrinkage degrees is calculated until the hydraulic jump of the central cylindrical rectangular groove with a certain shrinkage degree r is a normal hydraulic jump, and the hydraulic jump of the central cylindrical rectangular groove with a shrinkage degree less than 0.01 is a column Stop the calculation and output the shrinkage degree when the hydraulic jump is started. The shrinkage degree can be determined as the optimal shrinkage degree of the target rectangular tank and is recorded as r 0 . Observe the shape of the hydraulic jump in front of the central cylinder in the tank under the conditions of any shrinkage degree and arbitrary flow rate. , the judgment standard of hydraulic jump form is: when the hydraulic jump form is an on-column hydraulic jump, the shrinkage degree of the central cylindrical rectangular groove is the invalid shrinkage degree, and when the hydraulic jump form is not an on-column hydraulic jump, the shrinkage degree is the central cylindrical rectangle The effective shrinkage of the groove; 步骤2：对大于等于步骤1-3输出的最优收缩度r0，即r≥r0，中选取6个收缩度即r1、r2、r3、r4、r5和r6，r1～r6均匀分布于r0≤r＜1之间，对所选收缩度的中心圆柱矩形槽在不同流量条件下进行数值计算，输出对应的计算流量Q、圆柱正前端水深h及矩形槽收缩后过流宽度Bc，再采用SPSS软件对条件下圆柱正前端淹没水深h、流道收缩后的过流宽度Bc及对应的计算流量Q，结合理论公式Q＝a(Bc 2.5)(g0.5)(h/Bc)b进行非线性处理，分别输出不同收缩度r条件下的系数a和b，再采用Origin对系数a、b与r进行拟合分析并输出其函数关系，g为重力加速度；Step 2: For the optimal shrinkage degree r 0 greater than or equal to the output of steps 1-3, that is, r≥r 0 , select 6 shrinkage degrees, namely r1, r2, r3, r4, r5 and r6, and r1 to r6 are evenly distributed in the Between r 0 ≤ r < 1, carry out numerical calculation for the central cylindrical rectangular groove of the selected shrinkage degree under different flow conditions, and output the corresponding calculated flow Q, the water depth h at the front of the cylinder, and the overflow width Bc after the rectangular groove shrinks, Then use SPSS software to calculate the submerged water depth h at the front end of the cylinder, the overflow width B c after the channel shrinks and the corresponding calculated flow Q under the conditions, combined with the theoretical formula Q=a(B c 2.5 )(g 0.5 )(h/B c ) b is subjected to nonlinear processing, and the coefficients a and b under different shrinkage degrees r are output respectively, and then Origin is used to fit and analyze the coefficients a, b and r and output their functional relationship, and g is the acceleration of gravity; 步骤3：将系数a、b与r的拟合关系式带入理论公式，即可得到考虑收缩度影响的矩形槽中实时流量的表达式；Step 3: Bring the fitting relationship of the coefficients a, b and r into the theoretical formula, and then the expression of the real-time flow in the rectangular slot considering the effect of shrinkage can be obtained; 步骤4：将最优收缩度r0和对应的收缩后矩形槽过流宽度Bc0带入步骤3中考虑收缩度影响的矩形槽中实时流量的表达式，即可得到基于最优收缩度中心圆柱矩形槽中实时流量的表达式；Step 4: Bring the optimal shrinkage degree r 0 and the corresponding rectangular groove overflow width B c0 after shrinkage into the expression of the real-time flow rate in the rectangular groove considering the effect of the shrinkage degree in step 3, and then the center based on the optimal shrinkage degree can be obtained. Expression of real-time flow in cylindrical rectangular groove; 步骤5：对需要测流的渠道进行实地勘测，量测渠道底宽W，对步骤1中最优收缩度中心圆柱矩形槽进行放大得到原型最优收缩度中心圆柱矩形槽，放大倍数为λ0＝W/B，将中心圆柱矩形槽安置在渠道中，当水流稳定后，测得中心圆柱正前端淹没水深h0和收缩宽度W-W*r0带入步骤4中得到的基于最优收缩度中心圆柱矩形槽中实时流量的表达式即可得到实时的渠道流量。Step 5: Conduct on-the-spot survey of the channel that needs to measure the flow, measure the bottom width W of the channel, and enlarge the central cylindrical rectangular groove with the optimal shrinkage degree in step 1 to obtain the prototype optimal shrinkage degree central cylindrical rectangular groove, and the magnification is λ 0 =W/B, place the central cylindrical rectangular groove in the channel, when the water flow is stable, the submerged water depth h 0 and shrinkage width WW*r 0 at the front end of the central cylinder are measured and brought into the center based on the optimal shrinkage degree obtained in step 4. The real-time channel flow can be obtained by the expression of the real-time flow in the cylindrical rectangular groove. 2.根据权利要求1所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于：所述步骤1-1中，矩形槽长度L不小于10D，中心圆柱垂直于来流正方向的直径所在过流断面距进口断面的长度Lin不小于7D、距出口断面的长度Lout不小于3D，其中D为中心圆柱直径，中心圆柱顺水流正方向直径所在断面与矩形槽中心对称面重叠。2. a kind of channel flow measurement method based on the optimal shrinkage degree central cylindrical rectangular groove according to claim 1, is characterized in that: in described step 1-1, the rectangular groove length L is not less than 10D, the central cylinder is perpendicular to The length L in of the flow section where the diameter in the positive direction of the incoming flow is located from the inlet section is not less than 7D, and the length L out from the outlet section is not less than 3D, where D is the diameter of the central cylinder, and the section where the diameter of the central cylinder is located in the positive direction of the water flow is the same as the rectangle. Slot center symmetry planes overlap. 3.根据权利要求1所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于：所述步骤1-1中，在不同流量条件下对中心圆柱矩形槽进行模型试验，流量的取值由小至大选取，随着流量的增大，矩形槽中水深逐渐加深，当矩形槽进口水深低于矩形槽高度0.2D时停止试验，D为圆柱直径，试验过程中通过电磁流量计测试矩形槽中的流量Qt，相应的中心圆柱正前端的淹没水深ht通过在圆柱正前端安装水尺进行读数，Qt与ht的下标t表示试验条件，Qt的单位为L/s，ht的单位为cm。3. a kind of channel flow measurement method based on the optimal shrinkage central cylindrical rectangular groove according to claim 1, is characterized in that: in described step 1-1, under different flow conditions, the central cylindrical rectangular groove is modeled In the test, the value of the flow rate is selected from small to large. As the flow rate increases, the water depth in the rectangular tank gradually deepens. When the inlet water depth of the rectangular tank is lower than the height of the rectangular tank by 0.2D, the test is stopped. D is the diameter of the cylinder. The flow rate Qt in the rectangular tank is tested by an electromagnetic flowmeter, and the submerged water depth ht at the front end of the corresponding central cylinder is read by installing a water gauge at the front end of the cylinder. The subscript t of Qt and ht represents the test condition, and the unit of Qt is L/ The unit of s and ht is cm. 4.根据权利要求1所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于：所述步骤1-2中在网格划分中采用不同尺度的六面体结构化网格，中心圆柱网格剖分时采用外拓扑结构，从而得到若干不同的网格方案，针对各个网格网格方案在Fluent计算时采用VOF方法，针对VOF方法选取不同的湍流模型及不同的数值算法，数值计算中采用的流量与步骤1中试验的输出流量一一对应，计算中采用的边界条件进口为速度进口、出口为压力出口、表面为压力进口、固体边界设置为wall，计算初始条件为空气占比为100％，在数值模拟求得的水深高度为后处理CFD-Post中水和空气各占50％自由液面条件下的水深hs，单位为cm，下标s为数模条件。4. a kind of channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove according to claim 1, is characterized in that: in described step 1-2, adopt the hexahedral structured net of different scales in mesh division The outer topology structure is used when the central cylindrical grid is divided, so that several different grid schemes are obtained. For each grid scheme, the VOF method is used in the Fluent calculation, and different turbulence models and different numerical values are selected for the VOF method. Algorithm, the flow rate used in the numerical calculation corresponds to the output flow rate of the test in step 1, the boundary conditions used in the calculation are the velocity inlet, the outlet is the pressure outlet, the surface is the pressure inlet, and the solid boundary is set to wall, and the initial conditions are calculated. The water depth obtained by numerical simulation is the water depth hs under the condition that water and air each account for 50% of the free liquid surface in the post-processing CFD-Post, the unit is cm, and the subscript s is the digital-analog condition. . 5.根据权利要求1所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于：将步骤1-2中得到试验条件下的水深ht与数模条件下的水深hs进行误差分析，其计算式为|ht-hs|/ht≤k％，式中||表示绝对值，当某个数值计算方案下所有k值均不大于7时输出网格划分方案及Fluent设置方案。5. a kind of channel flow measurement method based on the optimal shrinkage center cylindrical rectangular groove according to claim 1, is characterized in that: the water depth ht under the test condition and the water depth under the digital-analog condition are obtained in step 1-2 hs is used for error analysis. The calculation formula is |ht-hs|/ht≤k%, where || represents the absolute value. When all k values under a certain numerical calculation scheme are not greater than 7, the grid division scheme and Fluent are output. Set up scenarios. 6.根据权利要求1所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于，所述步骤3中，系数a与r的关系表现形式为二次项a＝A1r2+B1r+C1，系数b与r的关系表现形式为二次项b＝A2r2+B2r+C2，其中A1，A2，B1，B2，C1，C2为列项系数。6. a kind of channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove according to claim 1, is characterized in that, in described step 3, the relation expression form of coefficient a and r is quadratic term a= A 1 r2+B 1 r+C 1 , the relationship between coefficient b and r is expressed as quadratic term b=A 2 r2+B 2 r+C 2 , where A 1 , A 2 , B 1 , B 2 , C 1 , C2 is the column term coefficient. 7.根据权利要求6所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于，所述步骤4中，基于最优收缩度中心圆柱矩形槽中实时流量的表达式
Figure FDA0003397407150000031
r0和Bc0为最优收缩度和最优收缩度对应的矩形槽过流宽度。7. a kind of channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove according to claim 6, is characterized in that, in described step 4, based on the expression of real-time flow in the optimal shrinkage center cylindrical rectangular groove Mode
Figure FDA0003397407150000031
r 0 and B c0 are the optimal shrinkage and the corresponding rectangular slot overcurrent width. 8.根据权利要求6所述的一种基于最优收缩度中心圆柱矩形槽的渠道测流法，其特征在于，所述步骤5中，基于最优收缩度中心圆柱矩形槽中实时流量的表达式得到的实时渠道流量为：
Figure FDA0003397407150000032
8. a kind of channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove according to claim 6, is characterized in that, in described step 5, based on the expression of real-time flow in the optimal shrinkage center cylindrical rectangular groove The real-time channel traffic obtained by the formula is:
Figure FDA0003397407150000032
CN202111485577.2A 2021-12-07 2021-12-07 Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove Pending CN114184239A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111485577.2A CN114184239A (en) 2021-12-07 2021-12-07 Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111485577.2A CN114184239A (en) 2021-12-07 2021-12-07 Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove

Publications (1)

Publication Number Publication Date
CN114184239A true CN114184239A (en) 2022-03-15

Family

ID=80603659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111485577.2A Pending CN114184239A (en) 2021-12-07 2021-12-07 Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove

Country Status (1)

Country Link
CN (1) CN114184239A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116384289A (en) * 2023-06-05 2023-07-04 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) Method for predicting pier block type fishway flow through computational fluid dynamics
CN117829036A (en) * 2024-03-05 2024-04-05 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) Method for determining flow of downstream branch channels of T-shaped branch channels through CFD

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116384289A (en) * 2023-06-05 2023-07-04 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) Method for predicting pier block type fishway flow through computational fluid dynamics
CN116384289B (en) * 2023-06-05 2023-08-08 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) Method for predicting pier block type fishway flow through computational fluid dynamics
CN117829036A (en) * 2024-03-05 2024-04-05 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) Method for determining flow of downstream branch channels of T-shaped branch channels through CFD
CN117829036B (en) * 2024-03-05 2024-05-14 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) A method for determining the flow rate of downstream branch channels of T-shaped branched open channels using CFD

Similar Documents

Publication Publication Date Title
CN111369059A (en) Urban waterlogging prediction method and system based on waterlogging fast simulation coupling model
CN109460605B (en) A method for predicting the flow rate of a large-scale low-lift pump
CN110135069B (en) Method and device for acquiring silt characteristics during water delivery of water delivery tunnel and computer equipment
CN114184239A (en) Channel flow measurement method based on optimal shrinkage center cylindrical rectangular groove
CN111241757A (en) Dam break three-dimensional numerical simulation method for uranium tailing pond based on computational fluid mechanics
CN107145671A (en) A kind of numerical reservoir simulation method and system
Ran et al. Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D
CN114117965A (en) A Method for Measuring Flow in Irrigation Areas Based on the Gradient Effect of Half-cylindrical Rectangular Slots on Both Sides
CN114330156A (en) A channel flow measurement method based on the scale effect of the central cylindrical rectangular slot
Zeng et al. Applications of computational fluid dynamics to flow ratings at prototype spillways and weirs. I: Data generation and validation
CN114357913A (en) Open channel junction hydrodynamic simulation method based on momentum conservation and Riemann solution
Ben Said et al. Performance of rectangular labyrinth weir–an experimental and numerical study
CN104778356B (en) A kind of method for numerical simulation of Convention diffusion mass transport process
Sun et al. Study on the applicability of four flumes in small rectangular channels
Wang et al. Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area
CN109614757B (en) A method of moist chamber type pumping plant critical submergence depth is predicted by CFD
CN114925624B (en) A numerical simulation method for three-dimensional water flow in natural rivers
Li et al. Wind tunnel study on influences of morphological parameters on drag coefficient of horizontal non-uniform buildings
CN118171603B (en) A method for determining flow through rectangular vegetation weirs using CFD
CN113971381A (en) Method for predicting local flow field influence range of sudden-expansion pipeline through CFD
CN117574716A (en) A flow coefficient correction method for low-temperature vortex flowmeter based on numerical simulation
CN117350194A (en) A flow velocity distribution simulation model for double-layer rigid vegetation open channel water flow
CN110728064A (en) Modeling method of hydrothermal coupling in riparian undercurrent zone based on COMSOL math module
CN118261073A (en) Distributed two-dimensional coupling simulation method and system suitable for rain and flood disasters
CN107808021A (en) The drag computation method of fluid device based on CFD

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20220315