CN111767684B - Modeling method of optimized friction source item implicit format two-dimensional shallow water equation - Google Patents
Modeling method of optimized friction source item implicit format two-dimensional shallow water equation Download PDFInfo
- Publication number
- CN111767684B CN111767684B CN202010612654.5A CN202010612654A CN111767684B CN 111767684 B CN111767684 B CN 111767684B CN 202010612654 A CN202010612654 A CN 202010612654A CN 111767684 B CN111767684 B CN 111767684B
- Authority
- CN
- China
- Prior art keywords
- unit
- flow
- time
- friction
- source item
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 230000004907 flux Effects 0.000 claims abstract description 23
- 230000003746 surface roughness Effects 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 6
- 108010001267 Protein Subunits Proteins 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 3
- 238000001764 infiltration Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 210000003205 muscle Anatomy 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 7
- 239000002352 surface water Substances 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an optimized friction source item implicit format two-dimensional shallow water equation modeling method, which comprises the steps of firstly dispersing terrain data of a river basin by adopting rectangular grids in a structural grid, and assigning a value to each parameter according to boundary conditions; adopting a central-format space discrete two-dimensional shallow water equation by a finite volume method; dividing a solving area into a plurality of small misaligned subareas based on a limited volume method of Godunov format, adopting an approximation solving Riemann problem method proposed by Harten, lax and Leer, namely an HLLC Riemann operator, and controlling flux calculation on a unit interface; and then calculating a source item and a stability condition, and finally updating the single-width flow, so that the modeling method of the optimized friction source item implicit format two-dimensional shallow water equation is completed. The invention improves the simulation efficiency and the calculation precision of the surface water dynamic process and provides effective technical support for the prediction of the associated process based on the surface water.
Description
Technical Field
The invention belongs to the technical field of surface water numerical simulation methods, and particularly relates to an optimized friction source term implicit format two-dimensional shallow water equation modeling method.
Background
The two-dimensional shallow water equation numerical model based on dynamic waves can simulate surface runoff and associated processes thereof, such as urban rainfall flood process, sub-flood process, pollutant diffusion process and sediment transport process. The method can provide effective technical support for prediction and early warning of urban flood disasters, river flood control, non-point source pollution and water and soil loss.
The high-efficiency accurate simulation calculation result is a basis for guaranteeing the accuracy and timeliness of prediction early warning, the high-resolution terrain data can truly represent the terrain, and the model simulation precision is improved. Meanwhile, the high-resolution topographic data restricts the operation efficiency of the model, and one of the restriction factors is the calculation of friction source items in the model operation. Proper treatment of friction source terms is one of the important conditions to ensure the computational stability of the model at shallow waters. Currently, the active and effective friction source treatment is based on an implicit method. However, the iteration of implicit calculations greatly increases the computational burden of the model, especially for high resolution terrain over large scale areas. Therefore, development of a new scheme for treating friction source item with positive and efficient is needed. Provides technical support for disaster prevention, disaster reduction, water pollution prevention and ecological civilization construction.
Disclosure of Invention
The invention aims to provide an optimized friction source item implicit format two-dimensional shallow water equation modeling method, which improves the simulation efficiency and calculation accuracy of a surface water power process and provides effective technical support for prediction of a surface water-based associated process.
The technical scheme adopted by the invention is that an optimized friction source item implicit format two-dimensional shallow water equation modeling method is implemented according to the following steps:
step 1, dispersing terrain data of a river basin by adopting rectangular grids in a structural grid, and assigning a value to each parameter according to boundary conditions;
step 2, adopting a central-format space discrete two-dimensional shallow water equation by a finite volume method;
Step 3, dividing the solving area into a plurality of small non-coincident subareas based on a limited volume method of Godunov format, adopting a method for solving the Riemann problem by approximation proposed by Harten, lax and Leer, namely an HLLC Riemann operator, and controlling flux calculation on a unit interface;
Step 4, source item calculation: the source item part comprises a bottom slope item and a friction item, the bottom slope source item processing adopts integration of decomposing integration on a cell into sub-units, and the bottom slope source item of the grid is converted into flux on the surface of the grid by assuming linear change of the height and the water level of the bed bottom in the grid; the friction source item treatment adopts a friction source item calculation method;
step 5, calculating stability conditions: the stability condition is satisfied by the time step and the space step;
And 6, updating the single wide flow, and thus, completing the modeling method of the two-dimensional shallow water equation with the implicit format of the optimized friction source item.
The present invention is also characterized in that,
The two-dimensional shallow water equation in step 2 is expressed as follows:
Wherein: t represents a time step, in s; x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q represents a conservation variable, q x represents a single-width flow in the x direction, m 2/s,qx=uh;qy represents a single-width flow in the y direction, and m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; f represents the flux vector in the x-direction; g represents the flux vector in the y-direction; s represents a source item vector; s b represents a bottom slope source term; s f represents a friction source item; z b represents the height of the river bed, and the unit is m; c f represents the bed surface roughness coefficient, C f=gn2/h1 /3; n represents a Manning coefficient, the unit is s/m -1/3, eta represents the water level, and the unit is m; η=h+z b; i represents static rainfall in mm/h, i=i r-ii;ir represents rainfall intensity in mm/h; i i represents the infiltration rate in mm/h.
The step 3 is specifically as follows:
Step 3.1, firstly, performing spatial second-order precision reconstruction on the interface variables to form a Riemann problem, performing numerical reconstruction on the variables at the left side and the right side of the unit interface by adopting a MUSCL state interpolation method, and simultaneously limiting gradients by combining a Min-Mod limiter to effectively inhibit numerical oscillation;
and 3.2, introducing a non-negative water reconstruction method to ensure harmony and conservation of mass of numerical calculation at a dry-wet interface.
The friction resistance source term calculation formula in the step4 is specifically as follows:
Wherein: x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3.
The step 5 is specifically as follows:
under the space step length deltax, the calculation formula of the value of the time step length deltat is as follows:
Wherein: Δx min represents the minimum distance from the center of the grid cell to the corresponding interface, in m; u is the flow rate in m/s; g is the gravity acceleration, and the value is 9.8m/s 2; h is the depth of water in m.
The step6 is specifically as follows:
Step 6.1, the single-wide traffic is updated at a new time step:
qn+1=qn+Δq+ΔtSf (8)
Wherein: q n+1 represents the single wide flow of n+1, the unit is m 2/s,qn represents the single wide flow at the time of n, the unit is m 2/s;Δq=F(qn,hn)+ΔtSb, and the flow term is calculated through HLLC; f represents interface flux; s b represents the bottom slope source term flux; s f represents friction source term flux; Δt represents the time step in s;
And 6.2, substituting the calculation formula of the friction resistance source term in the step 4 into a formula (8) to obtain the following two-dimensional form:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)And/>The single wide flow in the x direction at the time of n+1 and n is calculated by HLLC in the y direction, wherein the unit is m 2/s;Δqy; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s.
X represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3.
Step 6.3, solving the steps (9) and (10), and calculating to obtain the single wide flow under the new time step as follows:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)Is the single wide flow rate at the time of n+1 in the y direction, and the unit is m 2/s,/>The unit is m 2/s;Δqy, which is the single wide flow in the y direction at the time n, and is the flow term calculated by HLLC in the y direction; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s; θ represents a parameter introduced to avoid denominator of 0, in m 2/s,θ=1.0e-12.
The invention has the beneficial effects that the modeling method of the optimized friction source item implicit format two-dimensional shallow water equation aims at solving the problem of unavoidable redundant iteration of the implicit scheme adopted by the existing numerical model in the aspect of processing friction source item calculation, and provides a new processing scheme for avoiding iteration, which simulates the surface runoff process based on dynamic waves under a structure limited volume numerical discrete frame, can accurately ensure the conservation of water quantity in the whole calculation process, and efficiently and stably simulate the motion process of the surface runoff under the conditions of large gradient and small water depth of a river basin.
Drawings
FIG. 1 is a schematic view of a "V" shaped basin;
FIG. 2 is a comparison of simulated and analyzed values of the "V" watershed slope flow process;
fig. 3 is a comparison of the simulated value and the analytical solution value of the "V" basin exit flow process.
Detailed Description
The invention will be described in detail below with reference to the drawings and the detailed description.
The invention discloses an optimized friction source item implicit format two-dimensional shallow water equation modeling method, which is implemented according to the following steps:
step 1, dispersing terrain data of a river basin by adopting rectangular grids in a structural grid, and assigning a value to each parameter according to boundary conditions;
step 2, adopting a central-format space discrete two-dimensional shallow water equation by a finite volume method;
the two-dimensional shallow water equation in step 2 is expressed as follows:
Wherein: t represents a time step, in s; x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q represents a conservation variable, q x represents a single-width flow in the x direction, m 2/s,qx=uh;qy represents a single-width flow in the y direction, and m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; f represents the flux vector in the x-direction; g represents the flux vector in the y-direction; s represents a source item vector; s b represents a bottom slope source term; s f represents a friction source item; z b represents the height of the river bed, and the unit is m; c f represents the bed surface roughness coefficient, C f=gn2/h1 /3; n represents a Manning coefficient, the unit is s/m -1/3, eta represents the water level, and the unit is m; η=h+z b; i represents static rainfall in mm/h, i=i r-ii;ir represents rainfall intensity in mm/h; i i represents the infiltration rate in mm/h.
Step 3, dividing the solving area into a plurality of small non-coincident subareas based on a limited volume method of Godunov format, adopting a method for solving the Riemann problem by approximation proposed by Harten, lax and Leer, namely an HLLC Riemann operator, and controlling flux calculation on a unit interface;
the step 3 is specifically as follows:
step 3.1, firstly, performing spatial second-order precision reconstruction on the interface variables to form a Riemann problem, performing numerical reconstruction on the variables at the left side and the right side of the unit interface by adopting a MUSCL (Monotonic Upstream Scheme for Conservation Laws) state interpolation method, and simultaneously limiting gradients by combining a Min-Mod limiter to effectively inhibit numerical oscillation;
And 3.2, introducing a non-negative water reconstruction method in the process of calculating the boundary flux by using the HLLC Riemann operator so as to ensure harmony and conservation of mass of numerical calculation at a dry-wet interface.
Step 4, source item calculation: the source item part comprises a bottom slope item and a friction item, the bottom slope source item processing adopts integration of decomposing integration on a cell into sub-units, and the bottom slope source item of the grid is converted into flux on the surface of the grid by assuming linear change of the height and the water level of the bed bottom in the grid; the friction source item treatment adopts a friction source item calculation method;
the friction resistance source term calculation formula in the step4 is specifically as follows:
Wherein: x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3.
Step 5, calculating stability conditions: the stability condition is satisfied by the time step and the space step;
the step 5 is specifically as follows:
under the space step length deltax, the calculation formula of the value of the time step length deltat is as follows:
Wherein: Δx min represents the minimum distance from the center of the grid cell to the corresponding interface, in m; u is the flow rate in m/s; g is the gravity acceleration, and the value is 9.8m/s 2; h is the depth of water in m.
And 6, updating single wide flow, namely finishing the modeling method of the two-dimensional shallow water equation with the implicit format of the optimized friction resistance source item, wherein the method comprises the following steps of:
Step 6.1, the single-wide traffic is updated at a new time step:
qn+1=qn+Δq+ΔtSf (8)
Wherein: q n+1 represents the single wide flow of n+1, the unit is m 2/s,qn represents the single wide flow at the time of n, the unit is m 2/s;Δq=F(qn,hn)+ΔtSb, and the flow term is calculated through HLLC; f represents interface flux; s b represents the bottom slope source term flux; s f represents friction source term flux; Δt represents the time step in s;
And 6.2, substituting the calculation formula of the friction resistance source term in the step 4 into a formula (8) to obtain the following two-dimensional form:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)And/>The single wide flow in the x direction at the time of n+1 and n is calculated by HLLC in the y direction, wherein the unit is m 2/s;Δqy; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s.
X represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3.
Step 6.3, solving the steps (9) and (10), and calculating to obtain the single wide flow under the new time step as follows:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)Is the single wide flow rate at the time of n+1 in the y direction, and the unit is m 2/s,/>The unit is m 2/s;Δqy, which is the single wide flow in the y direction at the time n, and is the flow term calculated by HLLC in the y direction; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s; θ represents a parameter introduced to avoid denominator of 0, in m 2/s,θ=1.0e-12.
Examples
Discrete and conventional parameter assignment of the drainage basin terrain grid: the watershed topography was discretized into a structured uniform grid with 5m cell side length for a total of 64,800 cells. The length and width of the slope are 1000m and 800m respectively, and the gradient is 0.05. The channel width was 20m, the riverbed slope was 0.02, and the sidewall design of the trench was high enough to prevent the effects of backwater, as shown in FIG. 1. Then, a rainfall event with constant rainfall intensity of 10.8mm/h and rainfall duration of 1.5h was simulated. The Manning coefficients of the slope and the channel were taken to be 0.015s/m 1/3 and 0.15s/m 1/3, respectively.
The numerical calculation effect is solved for comparing the optimized friction source item implicit format, in addition, the friction source item solving method of the explicit type and split point implicit format is respectively modeled in a numerical mode and is respectively simulated. The rainfall-runoff process of the V-shaped basin is simulated based on the model for processing the friction source item by the 3 different algorithms, wherein FIG. 2 is a comparison of a simulation value and an analytic solution value of the slope flow process of the V-shaped basin, and FIG. 3 is a comparison of a simulation value and an analytic solution value of the outlet flow process of the V-shaped basin. As is evident from a comparison of the flow process lines in fig. 2 and 3, the model for processing the friction source item based on the implicit and explicit algorithms of the split point cannot truly reflect the flow process of the slope and the drainage basin outlet. The algorithm for processing the friction source term by applying the new format of the implicit algorithm in the model can better simulate the flow process of the slope and the channel outlet, and the result is very close to the value of the analytic solution.
Claims (3)
1. An optimized friction source term implicit format two-dimensional shallow water equation modeling method is characterized by comprising the following steps:
step 1, dispersing terrain data of a river basin by adopting rectangular grids in a structural grid, and assigning a value to each parameter according to boundary conditions;
step 2, adopting a central-format space discrete two-dimensional shallow water equation by a finite volume method;
Step 3, dividing the solving area into a plurality of small non-coincident subareas based on a limited volume method of Godunov format, adopting a method for solving the Riemann problem by approximation proposed by Harten, lax and Leer, namely an HLLC Riemann operator, and controlling flux calculation on a unit interface;
Step 4, source item calculation: the source item part comprises a bottom slope item and a friction item, the bottom slope source item processing adopts integration of decomposing integration on a cell into sub-units, and the bottom slope source item of the grid is converted into flux on the surface of the grid by assuming linear change of the height and the water level of the bed bottom in the grid; the friction source item treatment adopts a friction source item calculation method;
the friction resistance source term calculation formula in the step 4 is specifically as follows:
Wherein: x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3;
step 5, calculating stability conditions: the stability condition is satisfied by the time step and the space step;
The step 5 specifically comprises the following steps:
under the space step length deltax, the calculation formula of the value of the time step length deltat is as follows:
Wherein: Δx min represents the minimum distance from the center of the grid cell to the corresponding interface, in m; u is the flow rate in m/s; g is the gravity acceleration, and the value is 9.8m/s 2; h is the depth of water, and the unit is m;
step 6, updating single wide flow, so far, completing the modeling method of the two-dimensional shallow water equation with the implicit format of the optimized friction resistance source item,
The step 6 specifically comprises the following steps:
Step 6.1, the single-wide traffic is updated at a new time step:
qn+1=qn+Δq+ΔtSf (8)
Wherein: q n+1 represents the single wide flow of n+1, the unit is m 2/s,qn represents the single wide flow at the time of n, the unit is m 2/s;Δq=F(qn,hn)+ΔtSb, and the flow term is calculated through HLLC; f represents interface flux; s b represents the bottom slope source term flux; s f represents friction source term flux; Δt represents the time step in s;
And 6.2, substituting the calculation formula of the friction resistance source term in the step 4 into a formula (8) to obtain the following two-dimensional form:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)And/>The single wide flow in the x direction at the time of n+1 and n is calculated by HLLC in the y direction, wherein the unit is m 2/s;Δqy; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s;
x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q x denotes the single-width flow in x direction, the unit is m 2/s,qx=uh;qy denotes the single-width flow in y direction, the unit is m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; s fx represents a friction source item in the x direction; s fy represents a friction source term in the y direction; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, and the unit is s/m -1/3;
Step 6.3, solving the steps (9) and (10), and calculating to obtain the single wide flow under the new time step as follows:
Wherein: For single wide flow in x direction at time n+1, the unit is m 2/s,/> For single wide flow in the x direction at time n, the unit is m 2/s;Δqx which is a flow term calculated by HLLC in the x direction; /(I)Is the single wide flow rate at the time of n+1 in the y direction, and the unit is m 2/s,/>The unit is m 2/s;Δqy, which is the single wide flow in the y direction at the time n, and is the flow term calculated by HLLC in the y direction; c f represents the roughness coefficient of the bed surface; h n represents the water depth at time n, and the unit is m; Δt represents the time step in s; θ represents a parameter introduced to avoid denominator of 0, in m 2/s,θ=1.0e-12.
2. The modeling method of an optimized friction source term implicit format two-dimensional shallow water equation according to claim 1, wherein the two-dimensional shallow water equation in step 2 is represented as follows:
Wherein: t represents a time step, in s; x represents a horizontal transverse coordinate, and y represents a longitudinal coordinate; q represents a conservation variable, q x represents a single-width flow in the x direction, m 2/s,qx=uh;qy represents a single-width flow in the y direction, and m 2/s,qy =vh; h represents the water depth, and the unit is m; u represents the average flow velocity in the x direction in m/s; v represents the average flow velocity in the y direction in m/s; f represents the flux vector in the x-direction; g represents the flux vector in the y-direction; s represents a source item vector; s b represents a bottom slope source term; s f represents a friction source item; z b represents the height of the river bed, and the unit is m; c f represents the bed surface roughness coefficient, C f=gn2/h1/3; n represents a Manning coefficient, the unit is s/m -1/3, eta represents the water level, and the unit is m; η=h+z b; i represents static rainfall in mm/h, i=i r-ii;ir represents rainfall intensity in mm/h; i i represents the infiltration rate in mm/h.
3. The modeling method of the two-dimensional shallow water equation with the implicit format of the friction drag source term according to claim 2, wherein the step 3 is specifically as follows:
Step 3.1, firstly, performing spatial second-order precision reconstruction on the interface variables to form a Riemann problem, performing numerical reconstruction on the variables at the left side and the right side of the unit interface by adopting a MUSCL state interpolation method, and simultaneously limiting gradients by combining a Min-Mod limiter to effectively inhibit numerical oscillation;
and 3.2, introducing a non-negative water reconstruction method to ensure harmony and conservation of mass of numerical calculation at a dry-wet interface.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010612654.5A CN111767684B (en) | 2020-06-30 | 2020-06-30 | Modeling method of optimized friction source item implicit format two-dimensional shallow water equation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010612654.5A CN111767684B (en) | 2020-06-30 | 2020-06-30 | Modeling method of optimized friction source item implicit format two-dimensional shallow water equation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111767684A CN111767684A (en) | 2020-10-13 |
CN111767684B true CN111767684B (en) | 2024-04-19 |
Family
ID=72723069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010612654.5A Active CN111767684B (en) | 2020-06-30 | 2020-06-30 | Modeling method of optimized friction source item implicit format two-dimensional shallow water equation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111767684B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111767684B (en) * | 2020-06-30 | 2024-04-19 | 西安理工大学 | Modeling method of optimized friction source item implicit format two-dimensional shallow water equation |
CN112257313B (en) * | 2020-10-21 | 2024-05-14 | 西安理工大学 | GPU acceleration-based high-resolution numerical simulation method for pollutant transportation |
CN117010232B (en) * | 2023-06-27 | 2024-06-14 | 西安理工大学 | Urban non-point source pollution whole process high-resolution simulation method based on GPU (graphic processing Unit) acceleration technology |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108561A (en) * | 2017-12-22 | 2018-06-01 | 广州地理研究所 | Mud-rock flow integrated disaster reduction method based on dynamic process and energy spectrum analysis |
CN108446502A (en) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | A method of obtaining basin unit line using full two-dimensional Shallow Water Equations |
CN109271672A (en) * | 2018-08-21 | 2019-01-25 | 长江水资源保护科学研究所 | A kind of river-lake-pumping plant influences each other the river water surface curve calculation method under acting on |
CN110147646A (en) * | 2019-06-13 | 2019-08-20 | 中国水利水电科学研究院 | A kind of method for stream processing excessively of the non-structural flood numerical simulation frame lower linear liquid-retaining structure of two dimension |
CN110414088A (en) * | 2019-07-10 | 2019-11-05 | 武汉大学 | In conjunction with the wader habitat suitability ambiguity of space angle evaluation method of hydrodynamic model |
CN110929443A (en) * | 2019-12-02 | 2020-03-27 | 中国水利水电科学研究院 | Two-dimensional flood simulation method based on high-precision terrain generalization |
CN110990926A (en) * | 2019-12-02 | 2020-04-10 | 中国水利水电科学研究院 | Urban surface building hydrodynamic simulation method based on area correction rate |
CN111768502A (en) * | 2020-07-08 | 2020-10-13 | 西安理工大学 | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology |
CN111767684A (en) * | 2020-06-30 | 2020-10-13 | 西安理工大学 | Optimized friction resistance source term implicit format two-dimensional shallow water equation modeling method |
CN113656926A (en) * | 2021-08-26 | 2021-11-16 | 河海大学 | Pipeline transient flow simulation method based on Schohl convolution approximation |
CN116451611A (en) * | 2023-04-20 | 2023-07-18 | 武汉大学 | Method and system for constructing and simulating bright full-flow simulation model |
CN117010232A (en) * | 2023-06-27 | 2023-11-07 | 西安理工大学 | Urban non-point source pollution whole process high-resolution simulation method based on GPU (graphic processing Unit) acceleration technology |
-
2020
- 2020-06-30 CN CN202010612654.5A patent/CN111767684B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108108561A (en) * | 2017-12-22 | 2018-06-01 | 广州地理研究所 | Mud-rock flow integrated disaster reduction method based on dynamic process and energy spectrum analysis |
CN108446502A (en) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | A method of obtaining basin unit line using full two-dimensional Shallow Water Equations |
CN109271672A (en) * | 2018-08-21 | 2019-01-25 | 长江水资源保护科学研究所 | A kind of river-lake-pumping plant influences each other the river water surface curve calculation method under acting on |
CN110147646A (en) * | 2019-06-13 | 2019-08-20 | 中国水利水电科学研究院 | A kind of method for stream processing excessively of the non-structural flood numerical simulation frame lower linear liquid-retaining structure of two dimension |
CN110414088A (en) * | 2019-07-10 | 2019-11-05 | 武汉大学 | In conjunction with the wader habitat suitability ambiguity of space angle evaluation method of hydrodynamic model |
CN110929443A (en) * | 2019-12-02 | 2020-03-27 | 中国水利水电科学研究院 | Two-dimensional flood simulation method based on high-precision terrain generalization |
CN110990926A (en) * | 2019-12-02 | 2020-04-10 | 中国水利水电科学研究院 | Urban surface building hydrodynamic simulation method based on area correction rate |
CN111767684A (en) * | 2020-06-30 | 2020-10-13 | 西安理工大学 | Optimized friction resistance source term implicit format two-dimensional shallow water equation modeling method |
CN111768502A (en) * | 2020-07-08 | 2020-10-13 | 西安理工大学 | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology |
CN113656926A (en) * | 2021-08-26 | 2021-11-16 | 河海大学 | Pipeline transient flow simulation method based on Schohl convolution approximation |
CN116451611A (en) * | 2023-04-20 | 2023-07-18 | 武汉大学 | Method and system for constructing and simulating bright full-flow simulation model |
CN117010232A (en) * | 2023-06-27 | 2023-11-07 | 西安理工大学 | Urban non-point source pollution whole process high-resolution simulation method based on GPU (graphic processing Unit) acceleration technology |
Non-Patent Citations (5)
Title |
---|
A deep learning technique-based automatic monitoring method for experimental urban road inundation;Hao Han, Jingming Hou, Ganggang Bai, Bingyao Li, Tian Wang, Xuan Li, Xujun Gao, Feng Su, Zhaofeng Wang, Qiuhua Liang and Jiahui Gong;《Journal of Hydroinformatics》;20210403;全文 * |
二维水动力及污染物输移过程耦合高分辨率数值模型研究;石宝山;《中国优秀硕士学位论文全文数据库》;20220115;全文 * |
基于动力波的高性能分布式水文模型及应用;郭凯华;《中国优秀硕士学位论文全文数据库》;20190815;全文 * |
基于自适应网格的水文水力耦合 模型研究;于汪洋;《中国优秀硕士学位论文全文数据库》;20220815;全文 * |
考虑矿山堆渣影响的山洪形成及洪水演进淹没过程数值模拟;李桂伊;《中国优秀硕士学位论文全文数据库》;20181115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111767684A (en) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111767684B (en) | Modeling method of optimized friction source item implicit format two-dimensional shallow water equation | |
Afshar et al. | Particle swarm optimization for automatic calibration of large scale water quality model (CE-QUAL-W2): Application to Karkheh Reservoir, Iran | |
Afshar et al. | Application of cellular automata to sewer network optimization problems | |
CN114117848B (en) | Basin water environment simulation prediction method and device based on multi-model coupling | |
Murillo et al. | A conservative 2D model of inundation flow with solute transport over dry bed | |
CN110765676B (en) | Watershed water quality simulation method based on stable flow field | |
Salmasi et al. | Discharge coefficients for rectangular broad-crested gabion weirs: experimental study | |
CN102521886A (en) | Three-dimensional simulation method for chemical vapor deposition process | |
CN117010232B (en) | Urban non-point source pollution whole process high-resolution simulation method based on GPU (graphic processing Unit) acceleration technology | |
Ding et al. | Water quality safety prediction model for drinking water source areas in Three Gorges Reservoir and its application | |
van Loon | Numerical methods in smog prediction | |
CN112464584A (en) | Method for estimating water level and flow of free surface flow | |
CN110990926B (en) | Urban surface building hydrodynamic simulation method based on area correction rate | |
Yang et al. | Exploration of the dynamic water resource carrying capacity of the Keriya River Basin on the southern margin of the Taklimakan Desert, China | |
CN103605859A (en) | Optimum design method for secondary sedimentation tank of sewage treatment work | |
Sridharan et al. | Local-inertial shallow water model on unstructured triangular grids | |
Van den Bout et al. | A breakthrough in fast flood simulation | |
Wu et al. | Simulation of lake-groundwater interaction under unsteady-state flow using the sloping lakebed method | |
CN110837684B (en) | Basin water quality real-time simulation method based on general calculation | |
Ding et al. | Spatial-temporal variation of permanganate index under accident conditions during different water periods in Zigui drinking water source area, Three Gorges reservoir area, China | |
Hou et al. | Cyclic feedback updating approach and uncertainty analysis for the source identification of DNAPL-contaminated aquifers | |
CN103870699A (en) | Hydrodynamics flood routing analogy method based on double-deck asynchronous iteration strategy | |
Khrapov et al. | Lagrange-Eulerian method for numerical integration of the gas dynamics equations: parallel implementation on GPUs | |
CN114003977B (en) | Electrode configuration optimization method for constructed wetland-microbial fuel cell system | |
CN106650293A (en) | AM nest sampling algorithm based groundwater model evaluation method |
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 |