The content of the invention
For Shortcomings in the prior art, the invention provides a kind of adaptive mesh dynamic weighting of two-phase flow pump mistake
Poor appraisal procedure, using the mark mode for being gradually reduced threshold value, i.e. encryption every time, the quantity of computing unit is with the reduction of threshold value
And reduce, the computations repeatedly just for some regions are so effectively prevent, thus reduce the number of iteration, Neng Gouzheng
Really and error larger unit completely is evaluated, so as to instruct being smoothed out for encrypted work.
The present invention is to realize above-mentioned technical purpose by following technological means.
A kind of adaptive mesh dynamic weighting error evaluation method of two-phase flow pump, comprises the following steps:
(1) by carrying out external characteristics experiment to two-phase flow pump, the lift and efficiency data of two-phase flow pump, generation model are obtained
Pump threedimensional model simultaneously carries out initial mesh division;
(2) CFD numerical computations are carried out;
(3) the velocity gradient value for the unit i that the CFD numerical results according to step (2) obtainWith
Pressure gradientThe Grad of the unit i isWherein λ1For the power of velocity gradient
Weight coefficient, λ2For the weight coefficient of barometric gradient, λ1+λ2=1;
(4) gradient difference value of the unit i isPlace is normalized in unit i gradient difference value
Reason, obtains normalized gradient differenceWhereinFor the maximum of all unit gradient difference values;
(5) to describedUnit more than given error threshold is demarcated and encrypted, the given mistake
Poor threshold value is between 0~1;
(6) repeat step (2)~(5), and the error threshold that step (5) gives every time is gradually reduced, until what is obtained raises
Untill journey and efficiency meet required precision compared with the lift described in step (1) and efficiency.
Further, step (3) λ1>λ2。
For two-phase flow pump, the error that flow field is characterized compared to barometric gradient velocity gradient is more effective (especially not
Steady flow region), therefore of the invention when establishing speed and barometric gradient error evaluation function, the weight coefficient of velocity gradient
More than barometric gradient.
In such scheme, the barometric gradient of the unit i is:
WhereinThe barometric gradients of respectively described unit i in the X, Y, Z direction.
In such scheme, the velocity gradient of the unit i is:
WhereinThe velocity gradients of respectively described unit i in the X, Y, Z direction.
Unit i velocity gradient and barometric gradient computational methods is in addition to the method provided in the present invention, in prior art
In also have other computational methods, do not list one by one herein.
Beneficial effects of the present invention:
(1) the drawbacks of using unitary variant assessment strategy for existing appraisal procedure, the present invention use speed and pressure
The mode that gradient is combined is assessed, therefore can preferably react the cell position being had a great influence to computational accuracy, so as to solve
Existing evaluation method of having determined can not consider the deficiency of speed and barometric gradient simultaneously.
(2) in order to solve different example situations can not accurate assigned error threshold value situation, the present invention using calculate it is each
Normalized mode is used during the gradient difference value of individual unit so that all unit gradient difference values are all between 0~1, so as to add
The fast ciphering process of adaptive mesh.
(3) unit is obtained more than set-point to gradient difference value by way of dynamic alignment error threshold value in addition to be demarcated simultaneously
Encryption, i.e., give a larger error threshold (close to maximum 1), the error threshold in subsequent iterative process in first time
Be gradually reduced until the adaptive polo placement of two-phase flow pump lift and efficiencies with experiment compared with meet given accuracy untill, enter
One step improves the efficiency and validity of the process of adaptive refinement, and efficiently avoid different examples can not accurate assigned error threshold
The problem of value.
(4) according to the present invention program, it is possible to achieve the carry out error evaluation to all kinds grid, and assess efficiency high.
Embodiment
Below in conjunction with the accompanying drawings and specific embodiment the present invention is further illustrated, but protection scope of the present invention is simultaneously
Not limited to this.
Velocity gradient and barometric gradient of the present invention be in unit center of gravity (unit can be that tetrahedron can also
It is hexahedron).
If Fig. 1 is adaptive mesh generating process flow chart of the present invention.
Step 1, builds two-phase flow pump external characteristics testing stand, and the lift of pump has inlet and outlet pressure table measurement to obtain;Using electricity
The power of survey method measuring pump, the optimal operating condition point flow that external characteristics result obtains are 34.48m3/ h lifts are 4.76, and efficiency is
57.3%, three-dimensional modeling is carried out using Pro/E according to the hydraulic model of two-phase flow pump, and using ICEM generation initial mesh.
Step 2, calculating use CFX, standard k-ε turbulence model, SIMPLE algorithms, and non-coupled implicit aspect is solved.
Inlet boundary condition:Using speed import, assume to determine import axial velocity by mass conservation law and irrotationality.
Export boundary condition:Free discharge, it is assumed that flowing fully develops at outlet border, and exit region is apart from recirculating zone
Farther out.
Wall condition:Solid wall surface uses side wall non-slip condition.
Lift:
Efficiency:
In formula:PoutRepresent the pressure of pump discharge, PinPump inlet pressure is represented, ρ represents the density of liquid in pump, and Q represents pump
Flow, H represent pump lift, g is acceleration of gravity;M ' is the power of front side of vane, the back side and front and rear cover plate inner and outer surfaces
Square sum;ω represents the angular speed of pump;η ' is to include pump calculation domain forecasting efficiency value after volumetric loss, disc friction losses;
The loss of bearing and sealing takes 3%;Therefore efficiency eta=η ' × (1-3%) of pump is predicted.
Step 3, the velocity gradient value of unit is obtained according to result of calculationAnd pressure gradientAnd according to formulaCalculate the value of unit, wherein λ1For the weight system of velocity gradient
Number, λ2For the weight coefficient of barometric gradient, λ1+λ2=1, i be unit label, λ1Recommendation is 0.6, λ2Recommendation is 0.4;
1st, unit i pressure gradientComputational methods are as follows:
(1) unit i barometric gradient in the X direction:
Pi+1For the pressure of label i+1 units, Pi-1For the pressure of label i-1 units, Δ xI+1, iFor unit i+1 and unit i
In the distance (distance between unit center of gravity in X-direction) of X-direction, Δ xI, i-1For unit i and unit i-1 X-direction distance
(distance between unit center of gravity in X-direction);
Barometric gradients of the unit i in Y, Z-direction and the barometric gradient computational methods of the unit i in the X direction
Unanimously;
(2) unit i pressure gradient is:
WhereinThe barometric gradients of respectively described unit i in the X, Y, Z direction.
2nd, unit i velocity gradient valueComputational methods are as follows:
(1) present invention defines X, and Y, the speed of Z-direction is respectively u, v, w, then unit i velocity gradient calculates in X-direction
Method is as follows:
Wherein ui+1For the speed of label i+1 unit X-directions, ui-1For the speed of label i-1 unit X-directions, Δ xi+1,iFor
Unit i+1 and unit i is in the distance (distance between unit center of gravity in X-direction) of X-direction, Δ xi,i-1For unit i and unit i-1
In the distance (distance between unit center of gravity in X-direction) of X-direction.
(2) the velocity gradient size of X-direction is:
It is consistent with X-direction computational methods for the computational methods in other directions (such as Y, Z);
(3) it is in unit i velocity gradient size:
WhereinThe velocity gradients of respectively described unit i in the X, Y, Z direction.
3rd, on border unit barometric gradient computational methods:
Unit (such as i=0) on left margin, unit i barometric gradient computational methods are as follows in the X direction:
It is consistent with X-direction computational methods for the computational methods in other directions, therefore, in unit 0
Barometric gradient size be:
Unit (such as i=i on right marginmax, imaxRepresent largest unit number), unit i in the X directionmaxPressure ladder
It is as follows to spend computational methods,Computational methods and X-direction computational methods one for other directions
Cause, therefore, in unit imaxPoint barometric gradient size be:
4th, on border unit velocity gradient computational methods:
Unit (such as i=0) on left margin, the velocity gradient computational methods of unit 0 are as follows in the X direction:
The velocity gradient size of X-direction is:
It is consistent with X-direction computational methods for the computational methods in other directions (such as Y, Z), therefore, in the speed ladder of unit 0
Spending size is:
Unit (such as i=i on right marginmax, imaxRepresent largest unit number), unit i in the X directionmaxSpeed ladder
It is as follows to spend computational methods:
The velocity gradient size of X-direction is:
It is consistent with X-direction computational methods for the computational methods in other directions (such as Y, Z), therefore, in unit imaxSpeed
Gradient magnitude is:
Step 4:According to formulaThe gradient difference value of unit is calculated, and to the gradient of unit
Difference is normalized, i.e.,WhereinFor the maximum of all unit gradient difference values;
Step 5:The unit more than the value is demarcated according to given error threshold, if giving one every time
Fixed value can increase the time of whole adaptive process, therefore be pushed away in first time adaptive process to a larger error threshold
Recommend using 0.95, that is, demarcate all unitsValue be more than 0.95, and in subsequent calculating iterative process not
It is disconnected to reduce the value;
Step 6:The unit of all demarcation is encrypted, i.e., the midpoint on each side of demarcation unit adds an encryption
Node, then continue to carry out CFD calculating on this basis;
Step 7:Repeat step two arrives step 6, and the error threshold given in step 5 be gradually reduced (as the second time to
Fixed 0.9, third time given 0.88, by that analogy), until the CFD obtained lifts calculated and efficiency are less than 2% with test error
Untill (recommend calculation error be 2%), adaptive iteration process can be terminated.
The present embodiment is 860,000 in initial mesh, and the grid number obtained after adaptive refinement is the adaptometer after 1,830,000
Result and the external characteristics Comparative result of experiment are calculated, lift error is 1.2%, and efficiency error is 1.6%, and it is 6 hours to calculate the time
26 minutes;And if initial mesh is 1,850,000, the result of calculating is that lift error is 4.2%, and efficiency error is 5.1%, and 4 is small
When 14 minutes;And grid number reaches 3,980,000, can just obtain with adaptive mesh identical computational accuracy, and calculate the time be 16
43 minutes hours, it is seen that the method that this patent proposes can preferably improve computational efficiency and computational accuracy.
The embodiment is preferred embodiment of the invention, but the present invention is not limited to above-mentioned embodiment, not
Away from the present invention substantive content in the case of, those skilled in the art can make it is any it is conspicuously improved, replace
Or modification belongs to protection scope of the present invention.