CN106874539A - For the method and device of the filter disc flow passage structure design of laminated filter - Google Patents
For the method and device of the filter disc flow passage structure design of laminated filter Download PDFInfo
- Publication number
- CN106874539A CN106874539A CN201611267307.3A CN201611267307A CN106874539A CN 106874539 A CN106874539 A CN 106874539A CN 201611267307 A CN201611267307 A CN 201611267307A CN 106874539 A CN106874539 A CN 106874539A
- Authority
- CN
- China
- Prior art keywords
- runner
- model
- parameter
- flow
- structure model
- 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.)
- Granted
Links
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3323—Design verification, e.g. functional simulation or model checking using formal methods, e.g. equivalence checking or property checking
Abstract
The present invention relates to agricultural water-saving irrigation equipment technical field, more particularly to the method and device that a kind of filter disc flow passage structure for laminated filter is designed.Wherein, method is mainly included the following steps that：The first structure model of runner between adjacent filter disc is set up, numerical simulation calculation is carried out to first structure model, to determine the first flow parameter of runner；According to first flow parameter, runner correction term is preset；According to runner correction term, amendment is iterated to the computational fields border of first structure model by fractal algorithm, to obtain the second structural model；Reliability demonstration is carried out to the second structural model.Device mainly includes modeling unit, amending unit and authentication unit.The method and device can shorten and manufacture and design the cycle, reduce design cost, and the filter disc flow passage structure that its optimization design goes out has good filter effect, effectively reduce the loss of flood peak that filtering is produced；The method and device can also reliably, low consumption verify the validity of structure design simultaneously.
Description
Technical field
The present invention relates to agricultural water-saving irrigation equipment technical field, more particularly to a kind of filter disc stream for laminated filter
The method and device of road structure design.
Background technology
Slight irrigation is the Water Saving Irrigation Mode more popularized at present, and it significantly reduces Water Consumption in Agriculture, but this kind of
The irrigation technique of type is more strict for the requirement of water quality, if all kinds of impurity in irrigation water can not be removed effectively,
The quality of water can be had a strong impact on, the blocking of equipment at different levels is caused, serious meeting cause whole drip irrigation or irrigation system paralysis and
Scrap.
For the medium and small agricultural of China is irrigated, especially the Hetao Irrigation District of Huanghe water, the master in its diversion are being quoted
It is a large amount of fine sediments to want impurity.Therefore, in the irrigation system of slight irrigation, the filter efficiency of silt just becomes irrigation system is
No successful key.
At present, used as preferred filter in slight irrigation irrigation system, its core component is lamination, lamination to laminated filter
The plastic laminate upper and lower surface of formula filter is carved with a large amount of longitudinal direction microchannels and is laminated in inner support, by spring and fluid
Pressure compaction, the groove between lamination intersects, and forms a series of in-depth filtration unit of unique filtration channels.Laminated filter
During work, current are to enter lamination filter disc from filter core outer shroud to be flowed out by inner ring, and lamination is compressed by spring and Fluid pressure, pressure difference
Bigger, lamination pressure must be tighter, forms lamination self-locking performance high efficiency filter.
Currently for the research of laminated filter filter efficiency, the analysis to its automatic backwash is focused primarily upon, its
Middle field experiment and platform test have obtained more wide because of the guarantee of its accuracy in the optimization design of laminated filter
General application.But, the design processing for laminated filter needs to make corresponding mould, and the cycle that it manufactures and designs is very
It is long, and using experimental test to the inadequate economical rationality of method that the flow passage structure of filter disc in filter is optimized.
The content of the invention
(1) technical problem to be solved
The technical problem to be solved in the present invention there is provided a kind of filter disc flow passage structure design for laminated filter
Method and device, can shorten and manufacture and design the cycle, and the filter disc flow passage structure that its optimization design goes out has good filter effect,
Effectively reduce the loss of flood peak that filtering is produced.
(2) technical scheme
In order to solve the above-mentioned technical problem, the invention provides a kind of filter disc flow passage structure design for laminated filter
Method, comprise the following steps：
The first structure model of runner between adjacent filter disc is set up, numerical simulation calculation is carried out to the first structure model,
To determine the first flow parameter of the runner；
According to the first flow parameter, runner correction term is preset；
According to the runner correction term, the computational fields border of the first structure model is iterated by fractal algorithm
Amendment, to obtain the second structural model；
Reliability demonstration is carried out to second structural model.
Further, the described first structure model for setting up runner between filter disc, line number is entered to the first structure model
Value simulation is calculated, and to determine the first flow parameter of the runner, is further included：
The first structure model is built by Geometric Modeling；
Mesh generation is carried out to the first structure model；
To the computational domain setting boundary condition of the first structure model；
According to the boundary condition, the first structure model after mesh generation is carried out by computational fluid dynamics algorithm
Numerical simulation calculation, to determine the first flow parameter of the runner.
Further, it is described according to boundary condition, after computational fluid dynamics algorithm is to mesh generation first
Structural model carries out numerical simulation calculation, to determine the first flow parameter of the runner, further includes：
Determining the first flow parameter of the runner includes the turbulence model of fluid, pressure term and pressure x velocity coupling in runner
Close item；
The turbulence model is calculated by half permanent implicit Decoupled algorithm；
The pressure term is calculated using Pressure solution device, by Second-order Up-wind form；
The pressure x velocity coupling terms are calculated by SIMPLE algorithms.
Further, it is described according to first flow parameter, runner correction term is preset, further include：
According to first flow parameter, the filtering accuracy and suspen ded sediment load use parameter in the runner are preset.
Further, it is described according to runner correction term, by fractal algorithm to the computational fields of the first structure model
Border is iterated amendment, to obtain the second structural model, further includes：
According to the runner correction term, the computational fields border of the first structure model is iterated using KOCH curves
Amendment.
Further, it is described according to runner correction term, using KOCH curves to the computational fields of the first structure model
Border is iterated amendment, further includes：
Determine that the filtering accuracy in the runner includes cross section of fluid channel Internal tangential circular areas；
Determine whole corners on the computational fields border of the first structure model；
Keep the cross section of fluid channel Internal tangential circular areas constant, be by any one described corner iterated revision using KOCH curves
Three corners；
Chamfering is carried out to the revised whole corner, to obtain second structural model.
Further, described holding cross section of fluid channel Internal tangential circular areas are constant, using KOCH curves by any one side
Angle iterated revision is three corners, is further included：
Determine the triangular-section of the corner；
According to the filtering accuracy in the runner, respectively the predeterminated position iteration increase by two on two waists of the corner
Triangle wedge angle, so that the corner iterated revision is three corners；
According to the predeterminated position, the overflow section of the checking runner is repeatedly calculated, to determine that the runner excessively stream is cut
The optimum position of wedge angle when face is maximum.
Further, it is described that reliability demonstration is carried out to the second structural model, further include：
Numerical simulation calculation is carried out to second structural model, is joined with the second flow channel for determining the revised runner
Number；
Contrast the first flow parameter and second flow channel parameter；
Wherein, the second flow channel parameter includes the turbulence model of fluid, pressure term and pressure speed in revised runner
Degree coupling terms.
Further, it is described that reliability demonstration is carried out to the second structural model, further include：
Fluid dynamic simulation calculating is carried out to second structural model, is joined with the filtering for determining the revised runner
Number；
Contrast the filtration parameter and runner correction term；
Wherein, the filtration parameter includes the filtering accuracy and suspen ded sediment load use parameter in the revised runner.
Present invention also offers a kind of device of the filter disc flow passage structure design for laminated filter, including：
Modeling unit, for setting up filter disc between runner first structure model, numerical value is carried out to the first structure model
Simulation is calculated, to determine the first flow parameter of the runner；
Amending unit, for according to the first flow parameter, presetting runner correction term, according to the runner correction term,
Amendment is iterated to the computational fields border of the first structure model by fractal algorithm, to obtain the second structural model；
Authentication unit, for carrying out reliability demonstration to second structural model.
(3) beneficial effect
Above-mentioned technical proposal of the invention has the advantages that：Filter disc runner for laminated filter of the invention
In the method and device of structure design, method is mainly included the following steps that：The first structure model of runner between adjacent filter disc is set up,
Numerical simulation calculation is carried out to first structure model, to determine the first flow parameter of runner；According to first flow parameter, preset
Runner correction term；According to runner correction term, amendment is iterated to the computational fields border of first structure model by fractal algorithm,
To obtain the second structural model；Reliability demonstration is carried out to the second structural model.Device mainly includes modeling unit, amending unit
And authentication unit.The method and device can shorten and manufacture and design the cycle, reduce design cost, the filter disc stream that its optimization design goes out
Road structure has good filter effect, effectively reduces the loss of flood peak that filtering is produced；The method and device can be with can simultaneously
By ground, verify to low consumption the validity of structure design.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the method for the embodiment of the present invention；
Fig. 2 is the module connection diagram of the device of the embodiment of the present invention；
Fig. 3 is the structural section figure of the first structure model of experimental example of the present invention；
Fig. 4 is the structural section figure of the second structural model of the embodiment of the present invention；
Fig. 5 is the section speed vector figure of the first structure model of the embodiment of the present invention；
Fig. 6 is the section speed vector figure of the second structural model of the embodiment of the present invention；
Fig. 7 is the pressure-plotting of the first structure model of the embodiment of the present invention；
Fig. 8 is the pressure-plotting of the second structural model of the embodiment of the present invention；
Fig. 9 a~Fig. 9 j are the sediment charge distribution map of the second structural model of the embodiment of the present invention；
Wherein, Fig. 9 a are the runner interception schematic diagram of the second structural model of the embodiment of the present invention；
Fig. 9 b are the section sediment charge distribution map of target runner in Fig. 9 a；
Fig. 9 c are the enlarged diagram of a-quadrant in Fig. 9 b；
Fig. 9 d~Fig. 9 j are the section sediment charge distribution map of interception runner in Fig. 9 a.
Wherein, 100, modeling unit；200th, amending unit；300th, authentication unit；110th, corner region；111st, wedge angle；120、
Inscribed circle；130th, stringcourse；140th, target runner；150th, runner is intercepted.
Specific embodiment
Embodiments of the present invention are described in further detail with reference to the accompanying drawings and examples.Following examples are used for
The present invention is illustrated, but can not be used for limiting the scope of the present invention.
Embodiment one
The method and device of the filter disc flow passage structure design for laminated filter that the present embodiment one is provided, can shorten
Manufacture and design the cycle, reduce design cost, the filter disc flow passage structure that its optimization design goes out has good filter effect, effectively drop
The loss of flood peak that low filtering is produced；The method and device reliably, low consumption can also verify the effective of structure design simultaneously
Property.
Specifically, as shown in figure 1, the method includes at least including step S1, step S2, tetra- steps of step S3 and step S4
Suddenly.
S1, the first structure model for setting up runner between adjacent filter disc, numerical simulation calculation is carried out to first structure model, with
Determine the first flow parameter of runner, by the calculating of step S1, pressure field, velocity field equal distribution parameter can be obtained, it is determined that
Between laminated filter lamination in runner the flow feature of current, VELOCITY DISTRIBUTION, pressure distribution etc. influence sediment settlement factor with
And influence of the flow passage structure to filter efficiency.
Wherein, step S1 further comprises following steps：
S110, by Geometric Modeling build first structure model, the several of runner are preferably realized using PRO/E softwares during modeling
What is modeled, and channel size is measured using electron microscope and slide measure.
S120, mesh generation is carried out to first structure model, it is preferred to use GAMBIT carries out mesh generation to whole model,
GAMBIT be to aid in analyst and designer set up and gridding Fluid Mechanics Computation (CFD) model and design one is soft
Part bag, grid sets certain mesh refinement preferably using the tetrahedral grid of 0.04mm near wall region boundary.
S130, the computational domain setting boundary condition to first structure model, preferably by the flow inlet of runner and outlet
Boundary condition is respectively set to flow velocity entrance and free discharge.
S140, according to boundary condition, the first structure model after mesh generation is entered by computational fluid dynamics algorithm
Row numerical simulation calculation, to determine the first flow parameter of runner；Preferably, using FLUENT softwares by calculating fluid dynamic
Learning algorithm carries out numerical simulation calculation to the first structure model after mesh generation, and FLUENT softwares are a kind of CFD software bags, with
The relevant industrial circle such as fluid, heat transfer and chemical reaction can be used.
In step S140, for the ease of calculating, following steps have been still further comprised：
S141, determine the turbulence model of fluid in the first flow parameter including but not limited to runner of runner, pressure term and
Pressure x velocity coupling terms.
S142, turbulence model is calculated by permanent half implicit Decoupled algorithm, preferably turbulence model is Realizable
K- ε models, the model is one of default turbulence model in FLUENT softwares；In the half implicit Decoupled algorithm for using, fluid
Parameters be constant constant.
S143, using Pressure solution device, by Second-order Up-wind form calculate pressure term, Pressure solution device refers to based on pressure
Solver.
S144, by SIMPLE algorithms calculate pressure x velocity coupling terms, wherein, according to pressure term, by measurement stream track forces
The preset flow rate of the pressure term is calculated after very little；SIMPLE algorithms are that a kind of numerical value of flowing that can calculate any flow velocity is calculated
Method, the algorithm is applied to compressible flow field computation and incompressible fluid flows are calculated.
In above-mentioned step S140, calculated based on CFD technical modellings, velocity flow profile cloud atlas, pressure distribution cloud can be obtained
The analog results such as figure, Turbulent Kinetic intensity distribution cloud atlas, as shown in figure 3, can be with the cross section of fluid channel figure of VELOCITY DISTRIBUTION cloud atlas
Find out, the VELOCITY DISTRIBUTION of whole cross section of fluid channel is：The kernel of section speed of runner is maximum, more smaller to four circular velocities, especially exists
Each corner region 110, can form larger low regime；The pressure distribution of the overall runner of observation, it can be seen that by water inlet to
Delivery port is in the pressure value of uniform downward trend, and the cross sectional flow rate near exit is obviously more than near porch, and its turbulence
Intensity is stronger.By the analysis to above-mentioned first flow parameter, it can be seen that the preferable still phase of the grains of sand effect of settling of low regime
It is smaller to area.
S2, according to first flow parameter, preset runner correction term, wherein, runner correction term is included but is not limited in runner
Filtering accuracy and suspen ded sediment load use parameter, in the present embodiment, with the area table of each corner inscribed circle 120 of cross section of fluid channel
Show filtering accuracy.
So in step S2, further include：
S210, according to first flow parameter, preset the filtering accuracy and suspen ded sediment load use parameter in runner.
S3, according to runner correction term, amendment is iterated to the computational fields border of first structure model by fractal algorithm,
To obtain the second structural model, by changing the boundary shape of runner in first structure model, runner is increased with fractal algorithm
Overflow section area, to strengthen the efficiency that current pass through, while reduce the loss of flood peak in flow field, at the same time increase stream
The low regime in road section, increased the possibility of sediment settlement, effectively strengthen the efficiency of sediment settlement.
Specifically, by taking triangle runner common at present as an example, according to runner correction term, using KOCH curves to the first knot
The computational fields border of structure model is iterated amendment.
Wherein, step S3 is further included：
S310, the filtering accuracy determined in runner include the area of cross section of fluid channel inscribed circle 120；
S320, determine first structure model computational fields border whole corners；
S330, the constant area of holding cross section of fluid channel inscribed circle 120, be by any one corner iterated revision using KOCH curves
Three corners.
In other words, any one corner of runner is considered as a triangle, it is modified, repairing timing introduces a point shape
Algorithm, is iterated using KOCH curves to the border of computational fields, and alternative manner is preferably：KOCH curves are directly used in three
On two sides of angular corner, it is three that the quantity in the region of complete wedge angle 111 on corner is increased by one, when being improved
Keep the radius of inscribed circle 120 of cross section of fluid channel constant to ensure original filtering accuracy, due to the increase of wedge angle 111 so that section
Area increased, so as to increased low speed grains of sand decanting zone, and increased cross section of fluid channel product, reduce the loss of flood peak.
Specifically, in step S330, it is further comprising the steps：
S331, the triangular-section for determining corner；
S332, according to the filtering accuracy in runner, the predeterminated position iteration on two waists of corner increases by two three respectively
Angular wedge angle 111, so that corner iterated revision is three corners；Wherein, the predeterminated position on two waists of corner is respectively corner
Stringcourse 130 on four fourth class partial nodes in any one node, i.e. the predeterminated position of wedge angle 111 is：Two stringcourses of corner
On 130 1/3 at, at 1/2 and at 2/3 in one or two combinations.
The overflow section of runner is verified in the selected predeterminated position of wedge angle 111 of S333, basis, repeated multiple times calculating, to determine
The optimum position of wedge angle 111 when runner overflow section is maximum.
S340, chamfering is carried out to revised whole corners, to obtain the second structural model, preferred pair each corner rounding
Angle, rounding treatment is carried out by by the wedge angle 111 in original section, it is to avoid because easy accumulation is a large amount of at the wedge angle 111 of runner
Silt, it is difficult to clean up and by long-term accumulation, and causes the problem for blocking or even scrapping of filter, effectively improves
The efficiency of backwash, while can also be brought convenience for manufacture, reduces manufacturing cost.
S4, reliability demonstration is carried out to the second structural model, the reliability demonstration described in the embodiment of the present invention is included but not
It is limited to the checking of runner parameter comparison and the checking of runner filter effect.
Wherein, runner parameter comparison checking is the parameter by contrasting the runner flow field before and after optimization, is used to determine
The raising of the strainability of filter, determines the pressure loss of filter, it was demonstrated that the fact that it reduces the loss of flood peak.
The checking of runner parameter comparison is mainly included the following steps that：
S411, numerical simulation calculation is carried out to second structural model, to determine the second of the revised runner
Runner parameter, wherein, second flow channel parameter includes the turbulence model of fluid, pressure term and pressure x velocity coupling in revised runner
Close item.
Specifically, by changing the relevant parameter in the second structural model, such as predeterminated position of wedge angle 111 repeatedly, it is repeatedly right
Second structural model carries out numerical simulation calculation, so that it is determined that most preferred second structural model, according to most preferred second knot
Structure model determines the second flow channel parameter of revised runner.
Wherein, the filter disc of the laminated filter after the determination of most preferred second structural model is according to optimization is in difference
The measure of working condition, obtains corresponding data and determines under flow condition, specially：
The preset boundary conditions when determining, with flow velocity entrance as boundary condition, the outlet of runner is with freedom for the entrance of runner
It is boundary condition to go out stream, when the boundary condition of flow channel entry point is set into flow velocity entrance, is accumulated according to flow and cross section of fluid channel
Size is calculated, and the boundary condition of runner exit positioned at filter interior due to that cannot measure any numerical value, therefore is set to
Free discharge
S412, contrast first flow parameter and second flow channel parameter, such as contrast optimization before and after runner in turbulence model,
Pressure term and pressure x velocity coupling terms, its result of calculation are relative value, by the ginseng for contrasting the runner flow field before and after optimization
Number, to determine the raising of the strainability of filter, by the pressure loss for determining filter, it was demonstrated that it reduces the loss of flood peak
It is true.
The checking of runner filter effect is mainly included the following steps that：
S421, fluid dynamic simulation calculating is carried out to the second structural model, to determine the filtration parameter of revised runner；
Wherein, filtration parameter includes the filtering accuracy and suspen ded sediment load use parameter in revised runner.
S422, contrast filtration parameter and runner correction term.
Preferably, flow of fluid is calculated in the runner using two-phase simulated flow algorithm to the second structural model, for
Water sand mixed flow is calculated using Euler-Euler model, so that silt distribution in designed product is verified, to confirm sediment settlement
The raising of level.
Specifically, Eulerian equation of the filter disc according to the laminated filter after optimization under silt two-phase simulated flow algorithm,
Simulation calculates distribution situation of the silt in runner.Wherein, the same Realizable k- ε models from turbulence model, calculate
Method uses half implicit Decoupled algorithm of parameter constant, and pressure term uses Second-order Up-wind form, and used is based on pressure
Pressure solution device, pressure x velocity coupling terms are calculated using SIMPLE algorithms.In the setting of Euler's model, solid phase particles are set
It is the grains of sand (sand) that the numerical value such as corresponding density, particle phase volume fraction are set.Obtain result of calculation, by result of calculation with it is default
Runner correction term is contrasted, and the distribution according to sediment concentration judges whether the second model structure produces a desired effect.
The optimization design of the method filter disc structure suitable for various types of laminated filters described in the embodiment of the present invention,
Filter disc its strainability after optimization has raising largely.Due to the change of flow channel shape, the filter disc stream after optimization
Road structure has increased its area of section, can reduce the loss of flood peak in filter process, meanwhile, increased wedge angle 111 is
Filtering provides more sediment settlement regions, enhances filter effect, by the wedge angle 111 that rounding is designed, beneficial to recoil
Wash and be more easily manufactured, largely reduce the hidden danger of runner blocking, make it to reduce backwash time with long-play
Number and efficient before can keeping after backwashed.
From the foregoing, the method makes full use of the Visual calculation advantage of CFD, design cost greatly reduces, contract
The short construction cycle and improve ability to develop independently；The method is for improving the Optimization Design of laminated filter filter disc
Theory, making up the deficiency in laminated filter design theory has significant meaning and wide prospect.
Embodiment two
Based on the method described in embodiment one, the present embodiment two provides a kind of filter disc runner knot for laminated filter
The device of structure design, can be very good to implement the method as described in embodiment one, so as to realize to laminated filtering using the device
The flow passage structure that the filter disc of device is formed carries out reliable optimization, can shorten and manufacture and design the cycle, reduces design cost, and its optimization sets
The filter disc flow passage structure counted out has good filter effect, effectively reduces the loss of flood peak that filtering is produced；Simultaneously the method and
Device can also reliably, low consumption verify the validity of structure design.
As shown in Fig. 2 the device mainly includes modeling unit 100, amending unit 200 and authentication unit 3003, wherein, build
Form unit 100 is used to set up the first structure model of runner between filter disc, numerical simulation calculation is carried out to first structure model, with true
The first flow parameter in constant current road；Amending unit 200, for according to first flow parameter, runner correction term being preset, according to runner
Correction term, amendment is iterated by fractal algorithm to the computational fields border of first structure model, to obtain the second structural model；
Authentication unit 3003, for carrying out reliability demonstration to the second structural model.
The specific of the method for embodiment one and the device of embodiment two is further described below by way of experimental example
Structure and principle.
Experimental example
As described in embodiment one, the experimental provision of use is as described in embodiment two for the method being related in this experimental example.Below
By taking the structure of runner between the two neighboring filter disc of agricultural laminated filter as an example, the method for illustrating the optimization of its flow passage structure
And device.
1st, the first structure model of runner between adjacent filter disc is set up by modeling unit 100, first structure model is carried out
Numerical simulation calculation, to determine the first flow parameter of the runner.
As shown in figure 3, realizing the Geometric Modeling of runner using PRO/E softwares, the cross sectional shape of runner is triangle, altogether
Three corners, channel size includes：A diameter of 217.32mm of cross section of fluid channel inscribed circle 120, overflow section area is
0.0614mm2；Mesh generation is carried out to first structure model using GAMBIT softwares, it is the four sides of 0.04mm that grid uses yardstick
Volume mesh, and the near wall region boundary of the entrance, outlet and computational fields to runner sets mesh refinement.The boundary condition of entrance
Flow velocity entrance is set to, the boundary condition of outlet is set to free discharge, and ingress flow rate is with by the flow rate calculation of filter
Draw.
2nd, according to boundary condition, using modeling unit 100 after computational fluid dynamics algorithm is to mesh generation
One structural model carries out numerical simulation calculation, to determine the first flow parameter of runner.
In experiment, the flowing of fluid in the runner between two neighboring filter disc can be considered as typical VISCOUS FLOW, and can not
Compression, meets the continuity equation and momentum conservation equation of fluid motion；Using FLUENT simulation softwards to first structure model
Water movement is simulated, and selected turbulence model is by the Realizable k- ε models of low reynolds number amendment, computational methods
Using half implicit Decoupled algorithm of constant parameter, pressure term passes through Second-order Up-wind format algorithms using the solver based on pressure
Calculate and obtain, pressure x velocity coupling terms are calculated using SIMPLE algorithms.
3rd, the first flow parameter of runner is determined by modeling unit 100, specially：Calculate the flow field of first structure model
The parameter such as speed, pressure distribution, determine sediment settlement area, the determinant of sediment quantities deposited on the bed and flow path boundary structure for
The influence of filter efficiency and blockage problem.
As shown in Figure 5, Figure 7, the simulation result of calculation based on CFD software, has obtained velocity flow profile cloud atlas and pressure distribution
The parameter schematic diagram such as cloud atlas, the gradual change numerical value bar of lower section illustrates flow speed value, the pressure value representated by various color depths in figure
Isoparametric scope.
As seen from Figure 5, velocity flow profile has very strong regularity, the larger area of speed (figure middle section center in flow field
Locate to be light area, be expressed as region I) and speed compared with cell, (speed is shown as dark area at nearly wall, is expressed as region II, region
III and region IV).From center to external diffusion, flow velocity is gradually reduced, and middle main flow area flow velocity is higher, such as region I, mainly by
In water is crossed, the blending of current is taken into account；The flow velocity of near wall region is minimum, such as region II, region III and region IV, beneficial to the heavy of the grains of sand
Drop, forms low flow velocity area, and scouring effect is poor, the grains of sand is largely settled, and is the main utilization region of filtering, should give attention
And increase its region area as far as possible, especially as the region IV in Fig. 5, the region flow velocity are small, area is relatively large, strengthens silt
The effect of sedimentation, should give maximum reinforcing in optimization design.
4th, according to first flow parameter, default runner correction term is input into amending unit 200, according to runner correction term,
Amendment is iterated to the computational fields border of first structure model by fractal algorithm using amending unit 200, to obtain second
Structural model, the runner to the second structural model is modified analysis.
It can be seen from above-mentioned analysis, to improve the strainability of laminated filter, should try one's best the low flow velocity increased at nearly wall
Area, while also to ensure the filtering accuracy of filter：Should give reservation and strengthened in the corner of low flow velocity so that sediment settlement
Region further increase, and ensure filtering accuracy on the premise of (i.e. ensure overflow section the diameter of inscribed circle 120 not
Become), the area of overflow section is improved as far as possible, increase the inflow-rate of water turbine under identical flow velocity, the head reduced in filter process is damaged
Lose.
For the desilting advantage in the low flow velocity area of corner, flow passage structure is optimized：
As shown in figure 4, the triangle first to the runner of first structure model is modified, fractal algorithm is introduced, utilized
Koch curves are iterated amendment to cross section of fluid channel, and Koch curves are directly used on two stringcourses 130 of triangle, will be complete
It is three that the quantity of full wedge angle 111 is increased by one, keeps the radius of inscribed circle 120 of cross section of fluid channel constant to ensure original during amendment
There is filtering accuracy；Due to the increase of wedge angle 111 so that the area in section has increased, so as to increased low speed grains of sand decanting zone,
And the overflow section area of runner is increased, reduce the loss of flood peak.
Wherein, the predeterminated position of wedge angle 111 can be arranged on the midpoint of the triangle stringcourse 130 of cross section of fluid channel, three points
One of at place or 2/3rds etc., it is contemplated that effectively expand cross section of fluid channel and accumulate and increase low regime area, according under different pieces of information
Result, preferably by wedge angle 111 be arranged on the midpoint of stringcourse 130 be wedge angle 111 optimum position.
Then, original corner and newly-increased wedge angle 111 are carried out into rounding treatment respectively, preferably using size is
The circular arc of 0.1mm carries out rounding, carries out after rounding treatment, can be prevented effectively from wedge angle 111 and stockpile silt not easy cleaning
The drawbacks of, therefore the efficiency of backwash can be improved, while can also be brought convenience for manufacture.
As shown in figure 4, the cross section of fluid channel diameter of the second structural model that Fig. 4 is represented keeps 217.32mm constant, point shape is calculated
In method, the predeterminated position of wedge angle 111 is arranged at two 1/2nd of stringcourse 130 of optimum position, i.e. triangular-section,
As seen from Figure 4, now runner has three wedge angles 111, and the area of overflow section is 0.0656mm.
5th, carry out reliability demonstration using 300 pairs of the second structural models of authentication unit, described reliability demonstration include but
The checking of runner parameter comparison and the checking of runner filter effect are not limited to, and performance to runner after optimization is analyzed.
Wherein, the analysis result of runner parameter comparison checking is as follows：
Fig. 3 and Fig. 4 is contrasted respectively, Fig. 5 and Fig. 6 contrasts be can be seen that into the first structure model before optimization and optimization
The second structural model afterwards is compared, the area increase of the runner overflow section after optimization, and ensure that filtering accuracy is constant simultaneously,
At the same time, the area of section of the corner region 110 of the second structural model, the corner region than the first structure model before optimization
110 area of section increased more than twice, it is ensured that the expansion in effective sediment settlement area.
By Fig. 7 and Fig. 8 contrasts as can be seen that by after amendment, the pressure drop of the whole flow passage area of the second structural model has
It is obvious to reduce, do not changed significantly in the uniformity for dropping under stress, but the pressure differential of start-stop significantly decreases,
This represents that the loss of flood peak in whole flow process is obviously reduced.
The analysis result of runner filter effect checking is as follows：
Flow of fluid is calculated in runner by using two-phase simulated flow algorithm to the second structural model, husky for water
Mixed flow is calculated using Euler-Euler model, so that silt distribution in designed product is verified, to confirm sediment settlement level
Raising.
Simulation based on suspen ded sediment load use, using governing equation group of the SIMPLE Algorithm for Solving after discrete, each solution becomes
The calculation of amount uses Second-order Up-wind format algorithms, and setting solid phase particles are the grains of sand, and density is set to 2650kg/m3；It is determined that
Solid phase particles size uniform and it is shaped as spherical, (solid phase particles diameter can be according to required applicable situation for a diameter of 0.05mm
It is adjusted change)；Solid phase particles concentration distribution is configured again, it is 0.2% to set concentration, is simulated and calculates.
Fig. 9 a~Fig. 9 j are the sediment charge distribution map of the second structural model of the embodiment of the present invention, are flowed according to after optimization
From the point of view of the silt distribution situation in road section, due to reducing from flow inlet to the overflow section of outlet product is linear, cause unit plane
Flow velocity in product increases；And runner center belongs to main flow area, velocity ratio near wall region is big.
The entrance sediment concentration for setting runner is 0.2%, in the result of calculation of two-phase simulated flow algorithm, the body of silt phase
Product concentration distribution scope can be seen that the volume ratio of silt phase according to Fig. 9 a~Fig. 9 j middle sections silt distribution 0~0.26%
It is distributed and uneven in whole cross section, main flow area sediment concentration is smaller compared with corner region 110；According to Fig. 9 a~Fig. 9 j it can be pointed out that
Suspen ded sediment load use in single flow channel is stablized relatively, when lower flow channel coincidence process is run into motion process, shared by superposition section
Ratio is bigger, and the concentration of sand grain phase is higher in runner, wherein, superposition section refers to that neighbouring two runner bases mutually overlap
The part of conjunction；Similarly with the reduction of superposition section ratio, the concentration of sand grain phase also accordingly reduces；In its sand grain phase concentration
When higher, good condition has been created in the corner low regime of near wall region for the sedimentation of sand grain phase.
Specifically, the center of target runner is have chosen as shown in Fig. 9 a~Fig. 9 j, any interception with adjacent two filter disc
Staggeredly 7 of superposition intercept the sections of runner 150 between runner, by the volumetric concentration ratio of sand grain phase on 7 sections, point
Not Fen Xi silt staggeredly before, interleaved processes neutralize staggeredly after silt phase concentration the regularity of distribution.
When the porch sediment concentration be given in Fig. 9 a~Fig. 9 j is 0.2%, the silt in 7 sections being truncated to respectively
Phase volume concentration profile, interception has 7 schematic cross-sections of the interception runner 150 of section a~g respectively in Fig. 9 a, 7 in figure
The interception section a~g for intercepting runner 150 sequentially corresponds to shown in Fig. 9 d~Fig. 9 j respectively；In Fig. 9 a below interception runner 150
Runner be defined as target runner 140, as shown in Fig. 9 b and Fig. 9 c；Be given in Fig. 9 a color gradient contrast legend be Fig. 9 b~
The corresponding pictorial representation of each region labeling in Fig. 9 j, includes that 11 color scales are so that contrast is used altogether respectively in the legend,
Respectively yardstick 1~11.
Observation understands, in the section suspen ded sediment load use figure of the interception runner 150 shown in Fig. 9 d~Fig. 9 j, Fig. 9 d and Fig. 9 e
The runner of middle top is same runner, referred to as first runner；Fig. 9 f~Fig. 9 i are same runner, referred to as second runner；
Fig. 9 j are the 3rd runner.Wherein, because data shown in second runner are more comprehensive, so by taking second runner as an example.Work as mesh
It is higher away from sediment concentration at second wedge angle of runner 111 before mark runner 140 overlaps with second runner, and this feature
Runner rear end is always maintained at, is illustrated during being overlapped with another lamination groove runner, water in external flow passage
Stream can cause sand grain to move in the opposite direction in the turbulent flow that itself runner is caused, and be offset to point shape wedge angle in trapezoidal waist
At 111, the filter capacity of laminated filter can be effectively increased.
With the motion of current, the part that target runner 140 overlaps with second flow channel is more and more, according to the figure in Fig. 9 a
Example contrast understands, very unstable in the silt phase concentration distribution of overlapping position, increases yardstick 9 from yardstick 3, i.e., from 0.18%
Increase 0.32%；As intersection is more and more, silt phase concentration is also increasing.Illustrate sand grain in single runner
The volumetric concentration of phase is more stable, but when current run into the runner of other coincidences in runner, the concentration of sand grain phase will
Change, when overlapping more, concentration is higher；When current continue to move, superposition section again can fewer and feweri, silt phase concentration
Can reduce accordingly.
Can be seen that the wedge angle 111 that silt createed after optimization by above-mentioned reliability demonstration has obvious really
Aggregation situation, its silt aggregation extent is obviously improved than the runner before optimization, has reached the performance requirement of enhancing sediment settlement.
To sum up, in the method and device for the filter disc flow passage structure design of laminated filter of the present embodiment, method master
Comprise the following steps：The first structure model of runner between adjacent filter disc is set up, numerical simulation meter is carried out to first structure model
Calculate, to determine the first flow parameter of runner；According to first flow parameter, runner correction term is preset；According to runner correction term, lead to
Cross fractal algorithm and amendment is iterated to the computational fields border of first structure model, to obtain the second structural model；To the second knot
Structure model carries out reliability demonstration.Device mainly includes modeling unit 100, amending unit 200 and authentication unit 3003.The method
And device can shorten and manufacture and design the cycle, design cost is reduced, the filter disc flow passage structure that its optimization design goes out has good
Filter effect, effectively reduces the loss of flood peak that filtering is produced；The method and device can also reliably, low consumption be verified simultaneously
The validity of structure design.
One of ordinary skill in the art will appreciate that：Realize the complete of the method that the related each method embodiments of Fig. 1 are provided
Portion or part steps can be completed by the related hardware of programmed instruction, and foregoing program can be stored in a computer-readable
Take in storage medium, the program upon execution, performs the step of including above method embodiment；And foregoing storage medium bag
Include：ROM, RAM, magnetic disc or CD etc. are various can be with the medium of store program codes.
Through the above description of the embodiments, those skilled in the art can be understood that each implementation method can
Realized by the mode of software plus required general hardware platform, naturally it is also possible to by hardware.Based on such understanding, on
Stating the part that technical scheme substantially contributes to prior art in other words can be embodied in the form of software product, should
Computer software product can be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disc, CD, including some fingers
Order is used to so that a computer equipment (can be personal computer, server, or network equipment etc.) performs each implementation
Method described in some parts of example or embodiment.
Finally it should be noted that：Various embodiments above is only used to illustrate the technical scheme of embodiments of the invention, rather than right
Its limitation；Although being described in detail to embodiments of the invention with reference to foregoing embodiments, the ordinary skill of this area
Personnel should be understood：It can still modify to the technical scheme described in foregoing embodiments, or to which part
Or all technical characteristic carries out equivalent；And these modifications or replacement, do not make the essence disengaging of appropriate technical solution
The scope of each embodiment technical scheme of embodiments of the invention.
Claims (10)
1. a kind of method that filter disc flow passage structure for laminated filter is designed, it is characterised in that comprise the following steps：
The first structure model of runner between adjacent filter disc is set up, numerical simulation calculation is carried out to the first structure model, with true
The first flow parameter of the fixed runner；
According to the first flow parameter, runner correction term is preset；
According to the runner correction term, the computational fields border of the first structure model is iterated by fractal algorithm is repaiied
Just, obtaining the second structural model；
Reliability demonstration is carried out to second structural model.
2. method according to claim 1, it is characterised in that the described first structure model for setting up runner between filter disc,
Numerical simulation calculation is carried out to the first structure model, to determine the first flow parameter of the runner, is further included：
The first structure model is built by Geometric Modeling；
Mesh generation is carried out to the first structure model；
To the computational domain setting boundary condition of the first structure model；
According to the boundary condition, numerical value is carried out to the first structure model after mesh generation by computational fluid dynamics algorithm
Simulation is calculated, to determine the first flow parameter of the runner.
3. method according to claim 2, it is characterised in that described according to boundary condition, by calculating fluid dynamic
Learn algorithm carries out numerical simulation calculation to the first structure model after mesh generation, is joined with the first flow for determining the runner
Number, further includes：
Determining the first flow parameter of the runner includes the turbulence model of fluid, pressure term and pressure x velocity coupling in runner
；
The turbulence model is calculated by half permanent implicit Decoupled algorithm；
The pressure term is calculated using Pressure solution device, by Second-order Up-wind form；
The pressure x velocity coupling terms are calculated by SIMPLE algorithms.
4. method according to claim 1, it is characterised in that described according to first flow parameter, presets runner amendment
, further include：
According to first flow parameter, the filtering accuracy and suspen ded sediment load use parameter in the runner are preset.
5. method according to claim 4, it is characterised in that described according to runner correction term, by fractal algorithm pair
The computational fields border of the first structure model is iterated amendment, to obtain the second structural model, further includes：
According to the runner correction term, the computational fields border of the first structure model is iterated using KOCH curves is repaiied
Just.
6. method according to claim 4, it is characterised in that described according to runner correction term, using KOCH curves pair
The computational fields border of the first structure model is iterated amendment, further includes：
Determine that the filtering accuracy in the runner includes cross section of fluid channel Internal tangential circular areas；
Determine whole corners on the computational fields border of the first structure model；
Keep the cross section of fluid channel Internal tangential circular areas constant, using KOCH curves by any one described corner iterated revision be three
Corner；
Chamfering is carried out to the revised whole corner, to obtain second structural model.
7. method according to claim 6, it is characterised in that described holding cross section of fluid channel Internal tangential circular areas are constant, profit
With KOCH curves by any one described corner iterated revision be three corners, further include：
Determine the triangular-section of the corner；
According to the filtering accuracy in the runner, the predeterminated position iteration on two waists of the corner increases by two triangles respectively
Shape wedge angle, so that the corner iterated revision is three corners；
According to the predeterminated position, the overflow section of the checking runner is repeatedly calculated, to determine the runner overflow section most
The optimum position of wedge angle when big.
8. method according to claim 1, it is characterised in that described to carry out reliability demonstration to the second structural model, enters
One step includes：
Numerical simulation calculation is carried out to second structural model, to determine the second flow channel parameter of the revised runner；
Contrast the first flow parameter and second flow channel parameter；
Wherein, the second flow channel parameter includes the turbulence model of fluid, pressure term and pressure x velocity coupling in revised runner
Close item.
9. method according to claim 1, it is characterised in that described to carry out reliability demonstration to the second structural model, enters
One step includes：
Fluid dynamic simulation calculating is carried out to second structural model, to determine the filtration parameter of the revised runner；
Contrast the filtration parameter and runner correction term；
Wherein, the filtration parameter includes the filtering accuracy and suspen ded sediment load use parameter in the revised runner.
10. the device that a kind of filter disc flow passage structure for laminated filter is designed, it is characterised in that including：
Modeling unit, for setting up filter disc between runner first structure model, numerical simulation is carried out to the first structure model
Calculate, to determine the first flow parameter of the runner；
Amending unit, for according to the first flow parameter, presetting runner correction term, according to the runner correction term, passes through
Fractal algorithm is iterated amendment to the computational fields border of the first structure model, to obtain the second structural model；
Authentication unit, for carrying out reliability demonstration to second structural model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611267307.3A CN106874539B (en) | 2016-12-31 | 2016-12-31 | Method and device for designing filter sheet flow passage structure of laminated filter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611267307.3A CN106874539B (en) | 2016-12-31 | 2016-12-31 | Method and device for designing filter sheet flow passage structure of laminated filter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106874539A true CN106874539A (en) | 2017-06-20 |
| CN106874539B CN106874539B (en) | 2020-06-19 |
Family
ID=59164718
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611267307.3A Active CN106874539B (en) | 2016-12-31 | 2016-12-31 | Method and device for designing filter sheet flow passage structure of laminated filter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106874539B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763638A (en) * | 2018-04-20 | 2018-11-06 | 中国农业大学 | A kind of centrifugal type filter structural optimization method towards regulation and control internal flow |
| JP7422649B2 (en) | 2020-12-03 | 2024-01-26 | 大成建設株式会社 | Optimal design equipment and method for purification walls and water-stop walls |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5888399A (en) * | 1997-07-21 | 1999-03-30 | Rutledge; Dwight Dean | Water-removing funnel insert and method of application thereof |
| JP2002001027A (en) * | 2000-06-22 | 2002-01-08 | Tennex Corp | Multilayered dust collecting filter |
| CN1656311A (en) * | 2002-03-22 | 2005-08-17 | 克里萨里斯技术公司 | Fuel system for an internal combustion engine |
| CN104978456A (en) * | 2015-06-25 | 2015-10-14 | 北京信息科技大学 | General valve design optimization method |
| CN106082469A (en) * | 2016-06-14 | 2016-11-09 | 中国农业大学 | The silt of a kind of drip irrigation system regulates and controls method step by step |
-
2016
- 2016-12-31 CN CN201611267307.3A patent/CN106874539B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5888399A (en) * | 1997-07-21 | 1999-03-30 | Rutledge; Dwight Dean | Water-removing funnel insert and method of application thereof |
| JP2002001027A (en) * | 2000-06-22 | 2002-01-08 | Tennex Corp | Multilayered dust collecting filter |
| CN1656311A (en) * | 2002-03-22 | 2005-08-17 | 克里萨里斯技术公司 | Fuel system for an internal combustion engine |
| CN104978456A (en) * | 2015-06-25 | 2015-10-14 | 北京信息科技大学 | General valve design optimization method |
| CN106082469A (en) * | 2016-06-14 | 2016-11-09 | 中国农业大学 | The silt of a kind of drip irrigation system regulates and controls method step by step |
Non-Patent Citations (2)
| Title |
|---|
| 李云开 等: "《分形流道设计及几何参数对滴头水力性能的影响》", 《机械工程学报》 * |
| 李浩 等: "《微灌用叠片过滤器全流场数值模拟与优化分析》", 《灌溉排水学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108763638A (en) * | 2018-04-20 | 2018-11-06 | 中国农业大学 | A kind of centrifugal type filter structural optimization method towards regulation and control internal flow |
| CN108763638B (en) * | 2018-04-20 | 2022-08-19 | 中国农业大学 | Centrifugal filter structure optimization method for regulating and controlling internal flow |
| JP7422649B2 (en) | 2020-12-03 | 2024-01-26 | 大成建設株式会社 | Optimal design equipment and method for purification walls and water-stop walls |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106874539B (en) | 2020-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Johnson et al. | Normalized smoothing functions for SPH impact computations | |
| CN103886187B (en) | A kind of urban river water sand real-time predicting method based on data assimilation | |
| CN103761583B (en) | A kind of reservoir sand discharge generating Multiobjective Optimal Operation method | |
| CN109255185A (en) | An a kind of two-dimentional hydrodynamics coupling analytical method based on urban surface underground pipe network | |
| CN112257352A (en) | Coupling method and system of one-dimensional hydrodynamic model and two-dimensional hydrodynamic model | |
| CN106874539A (en) | For the method and device of the filter disc flow passage structure design of laminated filter | |
| CN107145678A (en) | A kind of rating method of Two Dimensional Plane Flow in Rivers model roughness | |
| CN104265652A (en) | Collaborative optimization design method for large-flow low-lift pump station | |
| CN101930483A (en) | Method for simplifying numerical map settlement place polygon by utilizing parametric design model | |
| CN103020361B (en) | A kind of no-checkerboard topological diagram from compliant mechanism extracting method | |
| CN103294898A (en) | Method for calculating single rod power of overall reactor core | |
| CN105512417A (en) | Particle tracking based three-dimensional migration simulation method for pore underground water pollutants | |
| CN101892836A (en) | Method for preparing large-scale porous percolating medium | |
| CN102129517A (en) | High-precision two-phase fluid interface capturing method | |
| CN103810326A (en) | Fishpond flow state digital analysis and construction method | |
| US20190110413A1 (en) | Drip Irrigation Emitter Flow Channel Structural Design Method and Fractal Flow Channel Drip Irrigation Emitter Product Therefor | |
| CN106844963A (en) | Excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process | |
| CN111177945B (en) | Optimization design method of centrifugal filter | |
| Xiao et al. | 2D numerical modelling of meandering channel formation | |
| Li et al. | Multi-objective optimal allocation of sediment resources under multiple uncertainties | |
| Asadsangabi et al. | Two phase flow modeling in shaft-spillways using volume of fluid (VOF) method | |
| CN110414075B (en) | CFD-based inversion calculation method for optimizing site selection of urban contradiction area | |
| CN111553010A (en) | BIM-based intelligent design method for gravity dam | |
| Chang et al. | Lattice Boltzmann simulation of fluid flows with fractal geometry: An unknown-index algorithm | |
| CN115659859A (en) | Flood simulation calculation method for river network structure containing flood area |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |