CN105760655B - A method of obtaining the fully developed flow boundary condition of channel entrance - Google Patents
A method of obtaining the fully developed flow boundary condition of channel entrance Download PDFInfo
- Publication number
- CN105760655B CN105760655B CN201610070282.1A CN201610070282A CN105760655B CN 105760655 B CN105760655 B CN 105760655B CN 201610070282 A CN201610070282 A CN 201610070282A CN 105760655 B CN105760655 B CN 105760655B
- Authority
- CN
- China
- Prior art keywords
- fully developed
- entrance
- type groove
- groove road
- channel entrance
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
The invention belongs to Computational Mechanics technical fields, are related to a kind of method for obtaining the fully developed flow boundary condition of channel entrance.This method comprises: firstly, before former channel entrance add one section of plane conduit, obtain modified conduit;Secondly, being based on set monitoring point, judge whether the length of additional plane conduit meets former channel entrance boundary and obtain the requirement of fully developed flowing velocity distribution, if being unsatisfactory for being adjusted;Judge whether modified channel entrance boundary amount of turbulence meets the requirement that former channel entrance boundary obtains fully developed flowing tubulence energy distribution, is adjusted if being unsatisfactory for;If satisfied, then completing a kind of method for obtaining the fully developed flow boundary condition of channel entrance.The mentioned method of the present invention can effectively reduce acquisition experimental data point and do curve matching bring error, simplify the operating process for obtaining fully developed distribution, fast and effeciently can obtain fully developed flowing entrance boundary condition in channel entrance.
Description
Technical field
The invention belongs to Computational Mechanics technical field, it is related to a kind of fully developed flow boundary condition of channel entrance that obtains
Method.
Background technique
U-shaped open channel flow (U-bend flow) is a kind of flowing with high current line curvature characteristic, and can be sent out in inner wall of the pipe
Raw a large amount of boundary layer separations and again phenomenon.These flow characteristics, which calculate numerical value, has very big challenge, therefore U-type groove road
Stream becomes the classical example of verifying numerical method.
Monson and Seegmiller in 1989 has carried out experimental study to U-type groove road stream and numerical value calculates.Numerical value calculates
In, the VELOCITY DISTRIBUTION and turbulence in the bend upstream somewhere that inlet flowing all assumes that as fully developed flowing, and experiment is measured
It can be distributed and be used as entrance boundary condition.Subsequent Shur and Strelets et al. propose to obtain by calculating plane conduit stream
Then these are distributed entering as U-type groove road stream calculation domain again by the VELOCITY DISTRIBUTION and tubulence energy distribution for obtaining fully developed flowing
Mouth boundary condition.However the acquisition computational domain entrance boundary of Monson and Seegmiller or Shur and Strelets proposition
There is some problems for VELOCITY DISTRIBUTION and the method for tubulence energy distribution:
(1) method that Monson and Seegmiller is proposed needs to carry out experimental data minimum quadratic fit, but this side
Method needs to obtain multiple data points, realizes that process is complex, and the curve and true velocity distribution and tubulence energy that fitting obtains
There are errors for distribution;
(2) method that Shur and Strelets are proposed needs Calculation Plane open channel flow example, then extraction rate distribution and
Tubulence energy distribution, and these distributions are assigned to U-type groove road inflow entrance boundary, it may be assumed that it needs to extract data and enters to U-type groove road
Mouth boundary assignment, realizes that process complexity is higher, causes Calculation Plane open channel flow and follow-up data to handle and brings additional work amount.
Thinking of the present invention is similar with (2) method, by the side for adding one section of plane conduit directly before U-type groove road inflow entrance
Method is distributed to obtain fully developed flowing velocity distribution required for entrance section and tubulence energy.Compared with (2) method, simplify
Obtain the operating process of fully developed flow distribution;Compared with (1) method, solves unnecessary mistake caused by curve matching
Difference.
Summary of the invention
The purpose of the invention is to reduce the operation complexity and reduction error that obtain fully developed flow distribution, mention
A kind of method for obtaining the fully developed flow boundary condition of channel entrance is gone out.
A method of the fully developed flow boundary condition of channel entrance being obtained, the specific implementation steps are as follows:
Step 1 adds one section of plane conduit before former channel entrance, obtains modified conduit;
Wherein, preferred former conduit is U-type groove road;
Preferred additional plane channel lengths range is 10 times to 15 times of former channel width;
Step 2 gives the entrance boundaries conditions such as conduit average speed, amount of turbulence in modified channel entrance, and in former slot
Monitoring point is set on road entrance section;
Step 3 judges whether the VELOCITY DISTRIBUTION on former channel entrance section meets fully developed flowing velocity distribution, and root
The case where according to whether satisfaction, carries out following corresponding operating:
3.1 if not satisfied, then former channel entrance section does not obtain required fully developed flowing velocity distribution, skips to step
1, continue the length for adjusting plane conduit;
3.2 if satisfied, fully developed flowing velocity distribution needed for then former channel entrance section obtains, skips to step 4;
Step 4, after step 3, i.e., the VELOCITY DISTRIBUTION on former channel entrance section meets fully developed flowing velocity distribution,
The amount of turbulence such as tubulence energy, the tubulence energy dissipative shock wave on modified channel entrance boundary are adjusted again;
Step 5 judges whether the tubulence energy distribution on former channel entrance section meets fully developed flowing velocity distribution, and
The case where according to whether satisfaction, carries out following corresponding operating:
5.1 if not satisfied, then skip to step 4, continues to adjust the tubulence energy on modified channel entrance boundary, tubulence energy dissipates
The amount of turbulence such as rate;
5.2 if satisfied, the fully developed flowing tubulence energy needed for then former channel entrance section obtains is distributed;
After step 5, i.e., the tubulence energy distribution on former channel entrance section has met fully developed flowing velocity distribution, from
This, a kind of method for obtaining the fully developed flow boundary condition of channel entrance is realized from step 1 to step 5.2.
Beneficial effect
A kind of method for obtaining the fully developed flow boundary condition of channel entrance of the present invention, comparison prior art have as follows
Advantage:
(1) acquisition experimental data point can be effectively reduced and do curve matching bring error;
(2) operating process for obtaining fully developed distribution is simplified, is extracted after Calculation Plane open channel flow can be effectively prevented from
VELOCITY DISTRIBUTION and tubulence energy distribution, and these distributions are assigned to the additional complex works such as U-type groove road inflow entrance boundary.
Detailed description of the invention
Fig. 1 is a kind of method flow diagram for obtaining the fully developed flow boundary condition of channel entrance of the present invention;
Fig. 2 is the tradition original U-type groove road geometric shape ABCD that compares with the mentioned method of the present invention;
Fig. 3 is modified U-type groove road geometric shape ABCDEF in the present invention and embodiment;
Fig. 4 is the former U-type groove road entrance section for obtain after numerical value calculating using modified U-type groove road geometric shape
The comparison of fully developed flowing velocity distribution and tubulence energy distribution and experimental result;
Wherein, " experiment " carries out the VELOCITY DISTRIBUTION and tubulence energy at the section CD obtained after numerical value calculating for former U-type groove road
Distribution;" present invention " is the section that modified U-type groove road geometric shape obtain after numerical value calculating in the present invention and embodiment
Fully developed flowing velocity distribution and tubulence energy distribution at CD;
Fig. 5 is the coefficient of friction and pressure that conduit inner and outer wall is calculated using modified U-type groove road geometric shape
The distribution curve of coefficient, " Shur and Strelets and Monson and the mentioned method of Seegmiller " and to former U-type groove road into
The Comparative result of row experiment three;
Wherein, " experiment " is the coefficient of friction on former U-type groove road ABCD inner and outer wall and pressure coefficient distribution;" this hair
It is bright " it is the method for obtaining the fully developed flow boundary condition of channel entrance a kind of in the present invention and embodiment using modified U-shaped
Conduit geometric shape carries out the coefficient of friction obtained after numerical value calculating and pressure coefficient distribution;" Shur and Strelets " is Shur
The result being calculated with Strelets and Monson and the Seegmiller method proposed.
Specific embodiment
It elaborates with reference to the accompanying drawings and examples to the present invention.
Embodiment 1
As Fig. 1 gives a kind of method flow diagram for obtaining the fully developed flow boundary condition of channel entrance of the present invention;Fig. 2
Give the present embodiment Central Plains U-type groove road geometric shape ABCD;Fig. 3 is modified U-type groove road geometric shape in the present embodiment
ABCDEF。
A kind of method for obtaining the fully developed flow boundary condition of channel entrance of the present embodiment, specifically includes the following steps:
Step 1 adds the plane conduit that a segment length is 15 times of channel widths before former U-type groove road ABCD entrance CD
CDEF obtains modified U-type groove road ABFE, as shown in Figure 3;
Step 2 gives the entrance boundaries conditions such as conduit average speed, amount of turbulence in modified U-type groove road ABFE entrance EF,
And 150 monitoring points are set on former U-type groove road ABCD entrance section CD;
Step 3 judges whether the VELOCITY DISTRIBUTION on former U-type groove road ABCD entrance section CD meets fully developed flowing velocity
Distribution, and according to meet whether the case where carry out following corresponding operating:
3.1 if not satisfied, then original U-type groove road ABCD entrance section CD does not obtain required fully developed flowing velocity point
Cloth skips to step 1, continues the length for adjusting plane conduit CDEF;
3.2 if satisfied, then original U-type groove road ABCD entrance section CD obtains required fully developed flowing velocity distribution, such as
Shown in Fig. 4 (a), step 4 is skipped to;
The amount of turbulence such as step 4, tubulence energy, the tubulence energy dissipative shock wave of the modified U-type groove road entrance boundary the ABFE EF of adjustment;
Step 5 judges whether the tubulence energy distribution on former U-type groove road ABCD entrance section CD meets fully developed flowing speed
Degree distribution, and according to meet whether the case where carry out following corresponding operating:
5.1 if not satisfied, then skip to step 4, continue to adjust the tubulence energy of the modified U-type groove road entrance boundary ABFE EF,
The amount of turbulence such as tubulence energy dissipative shock wave;
5.2 if satisfied, then original U-type groove road ABCD entrance section CD obtain needed for fully developed flowing tubulence energy distribution,
As shown in Fig. 4 (b);
Since then, from step 1 to step 5.2, final determine obtains fully developed flowing velocity and tubulence energy distribution in section CD
The length and the required given amount of turbulence on the EF of entrance boundary of required additional plane conduit CDEF, realizes a kind of acquisition
The method of the fully developed flow boundary condition of channel entrance.
Fig. 4 is the former U-type groove road entrance section for obtain after numerical value calculating using modified U-type groove road geometric shape
The comparison of fully developed flowing velocity distribution and tubulence energy distribution and experimental result;It is calculated by numerical value, respectively obtains former U-type groove
Coefficient of friction on road ABCD inner and outer wall and pressure coefficient distribution and Shur and Strelets and Monson and
The result that the method that Seegmiller is proposed is calculated.
Fig. 5 is the coefficient of friction and pressure that conduit inner and outer wall is calculated using modified U-type groove road geometric shape
The distribution curve of coefficient, " Shur and Strelets and Monson and the mentioned method of Seegmiller " and to former U-type groove road into
The Comparative result of row experiment three.
As can be seen from Figure 5 method of the invention has effectively obtained sending out completely required for U-type groove road entrance boundary
Open up flowing velocity distribution and tubulence energy distribution, and then the coefficient of friction for the U-type groove road inner wall and outer wall being calculated and pressure system
The result that number distribution is obtained with method used by Shur and Strelets et al. meets very good, and meets with experimental result
It obtains preferably, and the external wall friction coefficient for the conduit endpiece that method used by Monson and Seegmiller et al. calculates relatively is tested
As a result higher 19.5%.
Embodiment 2
Conduit class in the present embodiment method for obtaining the fully developed flow boundary condition of channel entrance a kind of to the present invention
Type is extended, and other than U-type groove road, is also possible to V-groove road, C-shaped channel road and M type groove road.
Embodiment described above is a kind of method for obtaining the fully developed flow boundary condition of channel entrance of the present invention
Preferred embodiment, the present invention should not be limited to the embodiment and attached drawing disclosure of that.It is all not depart from the present invention
The equivalent or modification completed under disclosed spirit, both falls within the scope of protection of the invention.
Claims (2)
1. a kind of method for obtaining the fully developed flow boundary condition of channel entrance, it is characterised in that: comprising steps are as follows:
Step 1 adds the plane conduit that a segment length is 10 to 15 times of channel widths before former U-type groove road ABCD entrance CD
CDEF obtains modified U-type groove road ABFE;
Step 2 gives the entrance boundary item based on conduit average speed and amount of turbulence in modified U-type groove road ABFE entrance EF
Part, and 150 monitoring points are set on former U-type groove road ABCD entrance section CD;
Step 3 judges whether the VELOCITY DISTRIBUTION on former U-type groove road ABCD entrance section CD meets fully developed flowing velocity point
Cloth, and according to meet whether the case where carry out following corresponding operation:
3.1 if not satisfied, then U-shaped original conduit ABCD entrance section CD does not obtain required fully developed flowing velocity distribution, jumps
To step 1, continue the length for adjusting plane conduit CDEF;
3.2 if satisfied, then original U-type groove road ABCD entrance section CD obtain needed for fully developed flowing velocity distribution, skip to step
Rapid 4;
Step 4, the tubulence energy of the modified U-type groove road entrance boundary the ABFE EF of adjustment, the amount of turbulence based on tubulence energy dissipative shock wave;
Step 5 judges whether the tubulence energy distribution on former U-type groove road ABCD entrance section CD meets fully developed flowing velocity point
Cloth, and according to meet whether the case where carry out corresponding operating, specifically:
5.1 if not satisfied, then skip to step 4, continues the tubulence energy and rapids that adjust the modified U-type groove road entrance boundary ABFE EF
Amount of turbulence based on kinetic energy dissipation rate;
5.2 if satisfied, the fully developed flowing tubulence energy needed for then original U-type groove road ABCD entrance section CD is obtained is distributed;
Since then, from step 1 to step 5.2, final determine obtains needed for fully developed flowing velocity and tubulence energy distribution in section CD
The length and the required given amount of turbulence on the EF of entrance boundary for wanting additional plane conduit CDEF, realize a kind of acquisition conduit
The method of the fully developed flow boundary condition of entrance.
2. a kind of method for obtaining the fully developed flow boundary condition of channel entrance as described in claim 1, feature also exist
In: the U-type groove road in step 1 to step 5.2 is also possible to V-groove road, C-shaped channel road and M type groove road.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610070282.1A CN105760655B (en) | 2016-02-02 | 2016-02-02 | A method of obtaining the fully developed flow boundary condition of channel entrance |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610070282.1A CN105760655B (en) | 2016-02-02 | 2016-02-02 | A method of obtaining the fully developed flow boundary condition of channel entrance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105760655A CN105760655A (en) | 2016-07-13 |
CN105760655B true CN105760655B (en) | 2019-10-11 |
Family
ID=56343011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610070282.1A Expired - Fee Related CN105760655B (en) | 2016-02-02 | 2016-02-02 | A method of obtaining the fully developed flow boundary condition of channel entrance |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105760655B (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104634810A (en) * | 2013-11-08 | 2015-05-20 | 国核华清(北京)核电技术研发中心有限公司 | High Rayleigh (Ra) number coupling heat-transfer characteristic measuring and evaluating device |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7954544B2 (en) * | 2007-11-28 | 2011-06-07 | Uop Llc | Heat transfer unit for high reynolds number flow |
-
2016
- 2016-02-02 CN CN201610070282.1A patent/CN105760655B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104634810A (en) * | 2013-11-08 | 2015-05-20 | 国核华清(北京)核电技术研发中心有限公司 | High Rayleigh (Ra) number coupling heat-transfer characteristic measuring and evaluating device |
Non-Patent Citations (2)
Title |
---|
《Turbulence measurements and mass transfer in fully developed flow in a triangular duct with a narrow apex angle》;Hiromoto Usui 等;《Int.J.Heat Mass Transfer》;19820531;第25卷(第5期);全文 * |
《反映强流动曲率效应的非线性湍流模型》;徐晶磊 等;《应用数学和力学》;20080115;第29卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105760655A (en) | 2016-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106202935B (en) | A kind of bearing calibration of Watershed Runoff forecast and its system | |
CN106547972B (en) | Fluid line method of topological optimization design based on parametrization Level Set Method | |
Li et al. | Existence of solutions to nonlinear Hammerstein integral equations and applications | |
Kalisch et al. | Instability of solitary waves for a nonlinearly dispersive equation | |
JP2017004103A (en) | Turbulence simulation method | |
CN105760655B (en) | A method of obtaining the fully developed flow boundary condition of channel entrance | |
CN104810791B (en) | A kind of differential protecting method based on FPGA | |
CN106897789B (en) | Wind power active power climbing event detection method and device | |
CN105117609A (en) | Dynamic weighing method based on generalized K-Means classification decision | |
CN101599871A (en) | A kind of SFARIMA network flow prediction method | |
CN103106332B (en) | A kind of analytical approach of uncertainty of measurement | |
CN109740963A (en) | A method of identification urban road network's key intersection | |
Kang et al. | Design of microirrigation laterals on nonuniform slopes | |
CN106294913B (en) | The method for improving components CALCULATION OF THERMAL result reliability | |
CN106991436A (en) | Noise spot detection method and device | |
CN105719300A (en) | Riverway main stream line detection method based on SNE manifold learning | |
KR100739083B1 (en) | Performance analysis method of mixed-flow impellers covering the low flow rate characteristics | |
CN107135126B (en) | Flow online identification method based on sub-flow fractal index | |
CN112131698B (en) | Method and device for obtaining optimized value of pipe-strip radiator | |
CN103926021A (en) | Method for determining electron temperature of plasma | |
CN107063370A (en) | A kind of time response compensation method of heat distribution type gas mass flow meter | |
CN105740628B (en) | A kind of method that boundary-layer marginal position is determined according to Flow Field Numerical Calculation result | |
CN117119465B (en) | Broadband signal parallel spectrum sensing method based on Adaboost | |
CN105846956B (en) | Block code code weight analysis recognition methods based on sliding window filtering | |
CN111951555B (en) | Traffic flow speed critical phase change interval estimation method based on RFID data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191011 Termination date: 20210202 |