CN110232222A - Deposited tube flow field analysis method and system - Google Patents

Deposited tube flow field analysis method and system Download PDF

Info

Publication number
CN110232222A
CN110232222A CN201910437839.4A CN201910437839A CN110232222A CN 110232222 A CN110232222 A CN 110232222A CN 201910437839 A CN201910437839 A CN 201910437839A CN 110232222 A CN110232222 A CN 110232222A
Authority
CN
China
Prior art keywords
fluid
flow
flow field
geometrical model
field analysis
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
Application number
CN201910437839.4A
Other languages
Chinese (zh)
Other versions
CN110232222B (en
Inventor
赵彦琳
姚军
刘敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201910437839.4A priority Critical patent/CN110232222B/en
Publication of CN110232222A publication Critical patent/CN110232222A/en
Application granted granted Critical
Publication of CN110232222B publication Critical patent/CN110232222B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

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

Abstract

The present invention provides a kind of deposited tube flow field analysis method and system, which comprises obtains the geometrical model of fluid in the deposited tube formed there are the round tube of sedimentary;Define and determine the Reynolds number of the fluid;Numerical value calculating and flow field analysis are carried out to the geometrical model according to the Reynolds number, the present invention provides a kind of geometrical models of deposited tube with sedimentary, and calculate the fluid flow state analyzed in deposited tube with the method for flow field analysis by numerical value.

Description

Deposited tube flow field analysis method and system
Technical field
The present invention relates to deposited tube fluid analysis technical field more particularly to a kind of deposited tube flow field analysis method and it is System.
Background technique
Oil, gas and water conveying pipe are common in daily life and industrial production, their flow section is mostly circle.It is more Mutually flowing is important component therein, is answered extensively in fields such as environmental protection, pharmacy, metallurgy, building materials, the energy, chemical industry, space flight With.If there are solid phase particles, such as gas-solid or liquid fixed double phase flow in pipeline, solid phase particles may occur under the effect of gravity Deposition.
During pipeline transportation, the solid phase particles in fluid media (medium) can be moved downward and be deposited under gravity On pipeline wall surface, one layer of sedimentary is formed over time, so that area of section is become smaller, and generate between sedimentary and curved wall Angle.Currently, there is no the analytical technology for fluid flow state in deposited tube in the prior art.
Summary of the invention
It is an object of the present invention to provide a kind of deposited tube flow field analysis method, provide a kind of with sedimentary The geometrical model of deposited tube, and the fluid flow state analyzed in deposited tube with the method for flow field analysis is calculated by numerical value.This Another of invention is designed to provide a kind of deposited tube flow field analysis system.It is yet a further object of the present invention to provide one kind Computer equipment.Another purpose of the invention is that providing a kind of readable medium.
In order to reach the goals above, one aspect of the present invention discloses a kind of deposited tube flow field analysis method, comprising:
Obtain the geometrical model of fluid in the deposited tube formed there are the round tube of sedimentary;
Define and determine the Reynolds number of the fluid;
Numerical value calculating and flow field analysis are carried out to the geometrical model according to the Reynolds number.
Preferably, the geometrical model of the fluid is horizontally disposed column, and the cross section profile of the column includes The circular arc of bottom end opening and the horizontally disposed straight line being separately connected with the both ends of the opening, the straight line and the circular arc Distance is determined according to the thickness of the sedimentary.
Preferably, numerical value calculating is carried out to the geometrical model according to the Reynolds number to specifically include:
Grid dividing is carried out to the geometrical model of the fluid;
Geometrical model after grid dividing is arranged according to Reynolds number obtains analysis model;
Numerical value calculating is carried out to the analysis model.
Preferably, grid dividing is carried out to the geometrical model of the fluid to specifically include:
Selected at least one trellis-type;
Geometrical model is divided into multiple regions or grid stroke is carried out to geometrical model by solving elliptic differential equation method Point;
The grid of boundary layer region is encrypted.
Preferably, the geometrical model being arranged after grid dividing according to Reynolds number obtains analysis model and specifically includes:
The fluid flow state ginseng of the fluid is determined according to the sectional dimension of the geometrical model of the fluid and Reynolds number Number;
The fluid flow state parameter is set on the geometrical model and constraint condition obtains the analysis model.
Preferably, constraint condition is arranged on the geometrical model to specifically include:
Boundary condition type is set;
Fluid inlet, outlet and the flow direction of fluid are set;
It is arranged without slip boundary condition.
Preferably, the flow field analysis includes average flow field analysis and transient flow field analysis.
The invention also discloses a kind of deposited tube flow field analysis systems, comprising:
Geometrical model unit, for obtaining, there are the geometrical models of fluid in the deposited tube of the round tube of sedimentary formation;
Parameter determination unit, for defining and determining the Reynolds number of the fluid;
Analytical unit is solved, for carrying out numerical value calculating and flow field analysis to the geometrical model according to the Reynolds number.
The invention also discloses a kind of computer equipment, including memory, processor and storage are on a memory and can The computer program run on a processor,
The processor realizes method as described above when executing described program.
The invention also discloses a kind of computer-readable mediums, are stored thereon with computer program,
The program realizes method as described above when being executed by processor.
The present invention constructs the geometrical model for analyzing fluid flow fields according to the design feature of deposited tube, and passes through numerical value The flow regime numerical value of the fluid in deposited tube is calculated, further the flow regime numeric distribution of deposition pipe fluid is carried out Flow field analysis is, it can be achieved that theoretically analyze the flow feature of fluid in the deposited tube with different deposit thickness, such as flow The influence of Secondary Flow and Secondary Flow stream field in, and then in research engineering the problems such as the particle deposition of pipeline transportation.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.
Fig. 1 shows a kind of one of flow chart of one specific embodiment of deposited tube flow field analysis method of the present invention;
Fig. 2 shows the two of a kind of flow chart of one specific embodiment of deposited tube flow field analysis method of the invention;
Fig. 3 shows a kind of schematic diagram of one specific embodiment geometrical model of deposited tube flow field analysis method of the invention;
Fig. 4 shows a kind of design principle of one specific embodiment geometrical model of deposited tube flow field analysis method of the invention Figure;
Fig. 5 shows the three of the flow chart of one specific embodiment of a kind of deposited tube flow field analysis method of the invention;
Fig. 6 shows the four of the flow chart of one specific embodiment of a kind of deposited tube flow field analysis method of the invention;
Fig. 7 shows the five of the flow chart of one specific embodiment of a kind of deposited tube flow field analysis method of the invention;
Fig. 8 shows the six of the flow chart of one specific embodiment of a kind of deposited tube flow field analysis method of the invention;
Fig. 9 shows the seven of the flow chart of one specific embodiment of a kind of deposited tube flow field analysis method of the invention;
Figure 10 shows geometrical model grid dividing in a kind of one specific embodiment of deposited tube flow field analysis method of the present invention Schematic diagram;
Figure 11 a shows a kind of numerical result middle section of one specific example of deposited tube flow field analysis method of the present invention One of streamline schematic diagram;
Figure 11 b shows a kind of numerical result middle section of one specific example of deposited tube flow field analysis method of the present invention The two of streamline schematic diagram;
Figure 11 c shows a kind of numerical result middle section of one specific example of deposited tube flow field analysis method of the present invention The three of streamline schematic diagram;
Figure 12 a shows in the numerical result of one specific example of deposited tube flow field analysis method of the invention a kind of and flows to One of speed cloud atlas schematic diagram;
Figure 12 b shows in the numerical result of one specific example of deposited tube flow field analysis method of the invention a kind of and flows to The two of speed cloud atlas schematic diagram;
Figure 12 c shows in the numerical result of one specific example of deposited tube flow field analysis method of the invention a kind of and flows to The three of speed cloud atlas schematic diagram;
Figure 13 a shows secondary in a kind of numerical result of one specific example of deposited tube flow field analysis method of the invention Flow one of cloud atlas schematic diagram;
Figure 13 b shows secondary in a kind of numerical result of one specific example of deposited tube flow field analysis method of the invention Flow the two of cloud atlas schematic diagram;
Figure 13 c shows secondary in a kind of numerical result of one specific example of deposited tube flow field analysis method of the invention Flow the three of cloud atlas schematic diagram;
Figure 14 shows a kind of structure chart of one specific embodiment of deposited tube flow field analysis system of the invention;
Figure 15 shows the structural schematic diagram for being suitable for the computer equipment for being used to realize the embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
According to an aspect of the present invention, present embodiment discloses a kind of deposited tube flow field analysis methods.As shown in Figure 1, In the present embodiment, the flow field analysis method includes:
S100: the geometrical model of fluid in the deposited tube formed there are the round tube of sedimentary is obtained.
S200: define and determine the Reynolds number of the fluid.
S300: numerical value calculating and flow field analysis are carried out to the geometrical model according to the Reynolds number.
The present invention constructs the geometrical model for analyzing fluid flow fields according to the design feature of deposited tube, and passes through numerical value The flow regime numerical value of the fluid in deposited tube is calculated, further the flow regime numeric distribution of deposition pipe fluid is carried out Flow field analysis is, it can be achieved that theoretically analyze the flow feature of fluid in the deposited tube with different deposit thickness, such as flow The influence of Secondary Flow and Secondary Flow stream field in, and then in research engineering the problems such as the particle deposition of pipeline transportation.
In a preferred embodiment, as shown in Fig. 2, the method also includes:
S000: the geometrical model of fluid in the building deposited tube formed there are the round tube of sedimentary.
It, can be under conditions of not appreciably affecting flow field analysis result, according to one in the geometrical model construction process of fluid Set pattern then is simplified to obtain the geometrical model of fluid in deposited tube to the fluid flowing in the deposited tube with sedimentary, specifically Rule of simplification may include several following:
(1) ignore the particle complex behavior of precipitate in the duct in the attachment of deposited tube curved wall and medium.
(2) assume that sedimentary forms straight line interface under the washing away of fluid.
(3) assume that sedimentary does not flow with fluid and thickness, shape and the change of position occurs.
As a result, when constructing geometrical model, it can first draw the circle for indicating deposited tube and indicate sedimentary and fluid circle The straight line in face, straight line intersect with circle, and the region for clipping the expression sedimentary part that straight line and circle are crossed to form obtains geometrical model Section is left part and is axially stretched along deposited tube, i.e., stretches the geometry mould to form the fluid of needs along the direction vertical with section Type.
Specifically, as shown in Figure 3 and Figure 4, the geometrical model of the fluid is horizontally disposed column, the column Cross section profile include the circular arc of bottom end opening and the horizontally disposed straight line that is separately connected with the both ends of the opening, it is described straight Line determines at a distance from the circular arc according to the thickness of the sedimentary.The top in section is arc boundary, and downside is straight line Boundary, the two form angle, and with moving up and down for downside boundary, length, straight line and the arc boundary of upper and lower lateral boundaries are formed Angle, area of section, perimeter of section, section hydraulic diameter etc. corresponding change occurs.It is different according to the thickness of sedimentary, it can The geometrical model of different cross section is formed, the geometrical model in the present embodiment has fully considered sedimentary fluid flow in deposited tube Influence, to ensure that the accuracy of flow field analysis result.
In a preferred embodiment, as shown in figure 5, the S300 can include:
S310: grid dividing is carried out to the geometrical model of the fluid.
S320: the geometrical model after grid dividing is arranged according to Reynolds number obtains analysis model.
S330: numerical value calculating is carried out to the analysis model.
In a preferred embodiment, as shown in fig. 6, the S310 is specific can include:
S311: selected at least one trellis-type.For example, due in deposited tube geometrical model include arc boundary, can adopt Grid dividing is carried out to geometrical model with O-shaped grid, to adapt to arc boundary, obtains the grid of high quality.Due to below circular arc For straight border and straight border position is not fixed, and grid carries out different disposal below the center of circle, for example, can be used hexahedron and Tetrahedral grid improves mesh quality, improves digital simulation result precision.
S312: geometrical model is divided into multiple regions or net is carried out to geometrical model by solving elliptic differential equation method Lattice divide.The geometrical model of deposited tube includes arc boundary and straight border, geometrical model can be divided into multiple regions, according to The geometry feature of different zones selects different grids to be divided.
S313: the grid of boundary layer region is encrypted.It is understood that in hydrodynamics, boundary layer The range in boundary layer is selected near wall in calculating process by the region of near wall flowing for artificially divide, fluid, Mesh refinement processing to boundary layer, it is preferred that boundary layer can be used exponential form and be encrypted.
In a preferred embodiment, as shown in fig. 7, the S320 is specific can include:
S321: determine that the fluid of the fluid flows shape according to the sectional dimension of the geometrical model of the fluid and Reynolds number State parameter.
S322: being arranged the fluid flow state parameter on the geometrical model and constraint condition obtains the analysis mould Type.
Wherein, fluid flow state parameter may include that fluid flows to average speed.Average speed and several is flowed to according to fluid The section of what model defines Reynolds number, can be used to show flow regime.When fluid is flowed along pipeline, it is flat that Reynolds number is defined as flow direction The ratio between the product and fluid kinematic viscosity of equal speed and pipeline hydraulic diameter.Reynolds number is less than critical Reynolds number, and flow field is laminar flow; Reynolds number is greater than critical Reynolds number, and flow field is turbulent flow.When Reynolds number is in excessive section, flow field will also be in critical state.When To in deposited tube fluid carry out Turbulent Flow Field analysis when, by determine fluid formed turbulent flow when Reynolds number, according to Reynolds number with The relationship that fluid flows to average speed can obtain the value for flowing to average speed, for being used as fluid flow state parameter.
In a preferred embodiment, as shown in figure 8, setting constraint condition can on the geometrical model in the S322 Include:
S3221: setting boundary condition type.Periodic boundary condition or cyclic asynchronous communication can be used.In order to anti- The underdeveloped influence to calculated result of fluid stopping body preferably uses periodic boundary condition.
S3222: fluid inlet, outlet and the flow direction of fluid are set.Specifically, fluid inlet is arranged on geometrical model With flow, speed and the pressure of outlet etc., flows to direction and be set as duct orientation.
S3223: it is arranged without slip boundary condition.I.e. geometrical model is taken along the boundary of deposited tube extending direction without sliding side Boundary's condition.
In a preferred embodiment, as shown in figure 9, the S330 is specific can include:
S331: pressure correction mould is obtained using numerical simulator according to the VELOCITY DISTRIBUTION of the fluid and pressure distribution Type.In order to improve accuracy, numerical simulator must embody the anisotropy of eddy stress.Navier- can preferably be solved Stokes (NS) equation can be used turbulence model and solve average flow field when analyzing Turbulent Flow Field, can also be using sub- grid Stress model solves large scale eddy, can be with the direct numerical solution of full size.When solution Navier-Stokes equation When, Reynolds stress transformation model may be selected and dissipative shock wave model is closed.Wherein, can presuppose initial VELOCITY DISTRIBUTION and Pressure distribution.
S332: it is modified according to pressure and speed of the pressure correction model to the fluid.According to the first of hypothesis Beginning VELOCITY DISTRIBUTION and pressure distribution discrete equation constant term and further acquire pressure correction side in evaluation simulation model Journey further solves number according to revised VELOCITY DISTRIBUTION and pressure distribution to be modified to VELOCITY DISTRIBUTION and pressure distribution According to discretization equations other in simulation model.
S333: repeating above step until meeting preset iterated conditional obtains the numerical result of fluid flowing.It is logical It crosses and repeats the continuous VELOCITY DISTRIBUTION repaired in positive fluid flow of S410 and S420 step realization iterative calculation and pressure distribution until full The preset iterated conditional of foot obtains the numerical result of fluid flowing, include in the numerical result in geometrical model not The flow regimes parameter such as flowing velocity with position fluid, according to fluid flow in the flow regime parameter of each point can analyze to obtain Fluid state in deposited tube.
The present invention is further described below by a specific example.In a specific example, pass through this The flow field analysis method of embodiment carries out numerical simulation and theory analysis to the liquid-solid two-phase flow field in round tube.Fluid in pipeline For liquid fixed double phase flow, the solid precipitation in pipeline in fluid makes round tube become deposited tube, and the sedimentary in deposited tube forms secondary The occurrence condition of stream.Secondary Flow be flowing of the fluid on flow section, with flow to whirlpool there are related.Secondary Flow is divided into two Class, the first kind are pressure-actuated Secondary Flow, and the second class is the Secondary Flow of turbulent flow driving.First kind Secondary Flow such as bend pipe is secondary Stream, it may be possible to the transverse direction whirlpool deflection as caused by crooked pipeline and generate;Second class Secondary Flow such as square tube Secondary Flow, its generation It is related with the eddy stress gradient that complex section generates.
The generation of Secondary Flow is influenced by arc, angle and chord length, complexity with higher, is the prediction of Secondary Flow Bring larger difficulty.In addition, turbulent secondary flow meeting stream field has an impact, the solid, gas-solid multiphase flow for liquid, thus it is possible to vary particle Cross-sectional distribution.Therefore, the research of Secondary Flow is not only rich in challenge, or pipeline particle deposition problems, which provide, to be solved to think Road, therefore there is theory significance and practical value.
Firstly, the geometrical model of the turbulent fluid of building deposited tube.Geometrical model is horizontally disposed column, the column The cross section profile of shape body includes the circular arc of bottom end opening and the horizontally disposed straight line that is separately connected with the both ends of the opening.Table The radius for showing the circular arc of deposited tube is r, diameter D, indicates geometrical model the length L=20D, deposit thickness H of pipe rangeb, H can be enabledb=0.25r, 0.5r or r is to form three kinds of geometrical models with different deposit thickness, according to sedimentary HbThickness It can determine the straight border and circular arc center of circle distance d, i.e. d=r-H of sedimentaryb
Three-dimensional system of coordinate is set, x, y, z respectively correspond the radial level of deposited tube to, it is radial vertically to and axial, coordinate Origin is the circular arc center of circle, and using right-handed coordinate system, fluid is along z-axis forward flow.
Grid dividing is carried out using geometrical model of the O-shaped grid to the turbulent fluid, to adapt to arc boundary, and opposite side The grid in interlayer region is encrypted, and obtains high quality grid to guarantee the accuracy of calculated result, as shown in Figure 10.
Conditions setting, it is true according to the sectional dimension of the geometrical model of the turbulent fluid and the corresponding Reynolds number of turbulent flow The fluid flow state parameter of the fixed turbulent fluid.Wherein, reynolds number Re is defined as:
Wherein, UbTo flow to average speed, DHFor section hydraulic diameter, A is area of section, and l is perimeter of section, and ν is fluid Kinematic viscosity.
According to above formula, for different deposition thickness, selected reynolds number Re=80000 carry out can be calculated turbulent flow Flow to average speed.
Setting boundary condition type is periodic boundary condition, and fluid inlet, outlet and the flow direction of turbulent fluid is arranged And it is arranged and obtains analysis model without slip boundary condition.
Equal NS equation when being based on to the calculating of analysis model in this example, with model be reynolds stress model (RSM) Turbulence model carries out sunykatuib analysis to analysis model and respectively obtains numerical result.
Specifically, this example is studied in three-dimensional level.It is incompressible in deposited tube, equal Navier- when being based on Stokes equation, the quality and momentum conservation equation of fully developed turbulent flow are writeable are as follows:
Wherein, i=1,2 or 3, j=1,2 or 3, wherein 1,2 and 3 respectively corresponds x, y, the direction z, uiFor the corresponding direction i Fluid velocity, ujFor the fluid velocity in the corresponding direction j, xiFor the position in the corresponding direction i, xjFor the position in the corresponding direction i, τijFor Eddy stress, ρ are fluid density, and t is time, giFor the mass force in the corresponding direction i, p is pressure.
Using reynolds stress model (RSM) Lai Jinhang numerical simulation, isotropism eddy viscosity is abandoned it is assumed that by solving thunder The transmission equation of promise stress and the equation of dissipative shock wave close the Navier-Stokes equation of Reynolds average.
When fluid is flowed along pipeline, by adjusting the thickness and Reynolds number of sedimentary, available difference flow field.Sedimentary When with a thickness of 0, pipeline is round tube, without secondary flow;When Reynolds number is smaller, flow field is laminar flow, also without secondary flow. When sedimentary reaches certain thickness, Reynolds number increases, and flow field reaches critical state, is finally changed into turbulent flow, at this moment can generate two Secondary flowing.
The flow field analysis includes average flow field analysis and flow field analysis.The flow field analysis knot simulated by logarithm Fruit can analyze Secondary Flow velocity vector in the distribution in section, such as size, and direction or position equal distribution information can also analyze rapids Kinetic energy, eddy stress, wall shear stress analyze their relationships between Secondary Flow, the production to Secondary Flow in the distribution in section Raw and development mechanism is analyzed and predicted, and illustrates the turbulence characteristics in section.
Work as Hb=0.5r when reynolds number Re=80000, obtains secondary streamline as shown in fig. 11a according to flow field analysis result, DH is section hydraulic diameter in figure.The Secondary Flow that will be close to center is referred to as inside Secondary Flow (such as 1 in Figure 11 a- Figure 11 c With 2), the Secondary Flow on both sides is referred to as outside Secondary Flow (such as 3 and 4 in Figure 11 a- Figure 11 c).Straight line near wall generates two Secondary whirlpool.Curved wall nearby without there is sizable vortex, angular bisector only near corner as interior outside The tangent line of Secondary Flow.Outside is only producing fairly small Secondary Flow near corner.Inside Secondary Flow is from center Yan Jiao Corner is flowed on the outside of bisector, then returns to center from straight line wall surface bisector.Deposit thickness HbIt flows in the section of=0.75r and r Line chart is respectively such as Figure 11 b and Figure 11 c.Outside Secondary Flow increases with the increase of deposit thickness.
The fluid flowing simulated by logarithm can also analyze turbulent flow and flow to speed in the distribution in section, observe secondary Flow the influence to speed is flowed to.Obtaining Secondary Flow according to flow field analysis result makes to flow to speed constant velocity line such as Figure 12 a- Figure 12 c institute Show.In deposited tube near wall, constant velocity line is parallel with wall surface.Since outside does not generate biggish vortex, Secondary Flow is in arc Near wall mainstream speed along wall surface tangential direction, deposited tube wall surface does not occur obviously to deform.On the contrary, the Secondary Flow of inside will The high-velocity fluid at center is transmitted to the corner that circular arc and straight line are formed, and by the low-speed flow of near wall at straight line wall surface bisector Body is passed back to center, and it is downwardly convex in close corner to make mainstream speed line, and raises upward at straight line wall surface bisector.When Deposit thickness Hb is gradually increased, and is flowed to bending degree of the contoured velocity in section and is gradually decreased.Flow to contoured velocity Bending degree is related with Secondary Flow, this shows that Secondary Flow reduces with the increase of deposit thickness at this.
Figure 13 a- Figure 13 c is Secondary Flow size cloud atlas, using flowing to speed nondimensionalization.It can be found that for different In the case of deposit thickness, there is maximum at sedimentary, angular bisector and camber line in Secondary Flow.But as Hb=0.5r and When 0.75r, Secondary Flow reaches maximum value at bisector, this embodies superposition of the inside whirlpool at center;And in Hb=r, There is maximum value at straight line wall surface in Secondary Flow, shows that the superposition of Secondary Flow at center weakens.Secondary Flow under different cross section Maximum value 2% or so, it is identical as square tube.For different cross section, it is secondary that the size of outside Secondary Flow is always less than inside Stream, further demonstrates that arc to the inhibiting effect of Secondary Flow.Secondary Flow reaches maximum in straight line wall surface, it is therefore desirable to avoid wall The use of surface function.Since outside Secondary Flow is suppressed development, far from corner, inside and outside sidespin whirlpool is no longer mutual with diagonal line Symmetrically, inside Secondary Flow obtains bigger development space, promotes the development of inside Secondary Flow.As deposit thickness increases Add, the Secondary Flow on diagonal line is gradually increased.This is because outside Secondary Flow is gradually increased, corner phase adjacent vortex size connects gradually Closely, the tangent line of phase adjacent vortex moves closer to angular bisector, finally forms the biggish Secondary Flow from center to corner in diagonal line.It is right The minimum of line two sides is claimed to represent vortex core.The comparison of section top and the bottom Secondary Flow it can be found that Secondary Flow by corner Place generates, and lower half portion is mainly distributed on, in section-top very little.
In short, the turbulent flow in the circle straight tube with sedimentary produces the second class Secondary Flow in section.Secondary Flow exists Corner generates and develops, and causes the variation for flowing to speed etc. on section.Arc has inhibiting effect to Secondary Flow, reduces Secondary Flow speed near curved wall;Smaller angle can promote turbulent flow reciprocation, be conducive to the generation of Secondary Flow;It is interior There are restricting effect between the Secondary Flow of outside, forms them and change and shifting effect with deposit thickness.As a result table Bright, there are the consistency of mechanism of production and geometry between the deposited tube Secondary Flow with sedimentary and the Secondary Flow in other flow fields Factor otherness caused by Secondary Flow.The present invention proposes the geometrical model of deposited tube earliest, considers that different deposition thickness flow down The different characteristic of field, deeply dissect geometrical factor influences caused by Secondary Flow.It, can by the analysis to Effects of Secondary Flow factor To realize the control to Secondary Flow, solution is provided to solve particle deposition relevant issues in pipeline.The Secondary Flow meeting of generation Stream field generates other influences, changes the original distribution of physical quantity in flow field.
Since there are this secondary flows, there is higher application value in engineer application.When demand secondary flow When, it can be generated using this circular pipe shape.By adjusting deposit thickness and flow velocity, can be obtained at center Certain speed vertically upward, and the acquisition of this speed is not necessarily to pipe bending or other conditions.This is the general pipe of solution Road deposition, blockage problem and other pipeline transportation problems and provide brand-new resolving ideas.
Based on same principle, the present embodiment also discloses a kind of deposited tube flow field analysis system.As shown in figure 14, deposited tube Flow field analysis system includes geometrical model unit 11, parameter determination unit 12 and solution analytical unit 13.
The geometry mould of fluid in the deposited tube that the geometrical model unit 11 is used to obtain being formed there are the round tube of sedimentary Type.
The parameter determination unit 12 is used to define and determine the Reynolds number of the fluid.
The solution analytical unit 13 is used to carry out numerical value calculating and flow field to the geometrical model according to the Reynolds number Analysis.Wherein, the flow field analysis includes average flow field analysis and flow field analysis.
The present invention constructs the geometrical model for analyzing fluid flow fields according to the design feature of deposited tube, and passes through numerical value The flow regime numerical value of the fluid in deposited tube is calculated, further the flow regime numeric distribution of deposition pipe fluid is carried out Flow field analysis is, it can be achieved that theoretically analyze the flow feature of fluid in the deposited tube with different deposit thickness, such as flow The influence of Secondary Flow and Secondary Flow stream field in, and then in research engineering the problems such as the particle deposition of pipeline transportation.
In a preferred embodiment, the geometrical model unit 11 is further used for constructing that described there are the circles of sedimentary The geometrical model of fluid in the deposited tube that pipe is formed.
The geometrical model unit 11 can not appreciably affect flow field analysis knot in the geometrical model construction process of fluid Under conditions of fruit, the fluid flowing in the deposited tube with sedimentary is simplified to obtain in deposited tube according to certain rule and is flowed The geometrical model of body, specific rule of simplification may include several following:
(1) ignore the particle complex behavior of precipitate in the duct in the attachment of deposited tube curved wall and medium.
(2) assume that sedimentary forms straight line interface under the washing away of fluid.
(3) assume that sedimentary does not flow with fluid and thickness, shape and the change of position occurs.
As a result, when constructing geometrical model, it can first draw the circle for indicating deposited tube and indicate sedimentary and fluid circle The straight line in face, straight line intersect with circle, and the region for clipping the expression sedimentary part that straight line and circle are crossed to form obtains geometrical model Section is left part and is axially stretched along deposited tube, i.e., stretches the geometry mould to form the fluid of needs along the direction vertical with section Type.
Specifically, as shown in Figure 3 and Figure 4, the geometrical model of the fluid is horizontally disposed column, the column Cross section profile include the circular arc of bottom end opening and the horizontally disposed straight line that is separately connected with the both ends of the opening, it is described straight Line determines at a distance from the circular arc according to the thickness of the sedimentary.The top in section is arc boundary, and downside is straight line Boundary, the two form angle, and with moving up and down for downside boundary, length, straight line and the arc boundary of upper and lower lateral boundaries are formed Angle, area of section, perimeter of section, section hydraulic diameter etc. corresponding change occurs.It is different according to the thickness of sedimentary, it can The geometrical model of different cross section is formed, the geometrical model in the present embodiment has fully considered sedimentary fluid flow in deposited tube Influence, to ensure that the accuracy of flow field analysis result.
In a preferred embodiment, the analytical unit 13 that solves is specifically used for carrying out the geometrical model of the fluid Grid dividing, the geometrical model after grid dividing is arranged according to Reynolds number obtain analysis model, count to the analysis model Value calculates.
In a preferred embodiment, the solution analytical unit 13 is further used for selected at least one trellis-type, Geometrical model is divided into multiple regions or grid dividing is carried out to geometrical model by solving elliptic differential equation method, to boundary The grid of layer region is encrypted.
Wherein, at least one trellis-type is selected.For example, due in deposited tube geometrical model include arc boundary, can adopt Grid dividing is carried out to geometrical model with O-shaped grid, to adapt to arc boundary, obtains the grid of high quality.Due to below circular arc For straight border and straight border position is not fixed, and grid carries out different disposal below the center of circle, for example, can be used hexahedron and Tetrahedral grid improves mesh quality, improves digital simulation result precision.
Geometrical model is divided into multiple regions or grid stroke is carried out to geometrical model by solving elliptic differential equation method Point.The geometrical model of deposited tube includes arc boundary and straight border, geometrical model can be divided into multiple regions, according to difference The geometry feature in region selects different grids to be divided.
The grid of boundary layer region is encrypted.It is understood that boundary layer is artificial in hydrodynamics Divide, fluid selectes the range in boundary layer near wall, to boundary in calculating process by the region of near wall flowing The mesh refinement processing of layer, it is preferred that boundary layer can be used exponential form and be encrypted.
In a preferred embodiment, the analytical unit 13 that solves is specifically used for according to the geometrical model of the fluid Sectional dimension and Reynolds number determine the fluid flow state parameter of the fluid, and the fluid stream is arranged on the geometrical model Dynamic state parameter and constraint condition obtain the analysis model.
Wherein, fluid flow state parameter may include that fluid flows to average speed.Average speed and several is flowed to according to fluid The section of what model defines Reynolds number, can be used to show flow regime.When fluid is flowed along pipeline, it is flat that Reynolds number is defined as flow direction The ratio between the product and fluid kinematic viscosity of equal speed and pipeline hydraulic diameter.Reynolds number is less than critical Reynolds number, and flow field is laminar flow; Reynolds number is greater than critical Reynolds number, and flow field is turbulent flow.When Reynolds number is in excessive section, flow field will also be in critical state.When To in deposited tube fluid carry out Turbulent Flow Field analysis when, by determine fluid formed turbulent flow when Reynolds number, according to Reynolds number with The relationship that fluid flows to average speed can obtain the value for flowing to average speed, for being used as fluid flow state parameter.
In a preferred embodiment, constraint condition is arranged in the solution analytical unit 13 on the geometrical model, tool Body, it can be used for being arranged boundary condition type, fluid inlet, outlet and the flow direction of fluid be set, be arranged without slip boundary item Part.
Specifically, periodic boundary condition or cyclic asynchronous communication can be used in boundary condition type.It flows in order to prevent The underdeveloped influence to calculated result of body preferably uses periodic boundary condition.
Fluid inlet, outlet and the flow direction of fluid are set, fluid inlet and outlet can be set on geometrical model Flow, speed and pressure etc. flow to direction and are set as duct orientation.
It is arranged without slip boundary condition, i.e. geometrical model takes no slip boundary item along the boundary of deposited tube extending direction Part.
In a preferred embodiment, the solution analytical unit 13 is particularly used in the VELOCITY DISTRIBUTION according to the fluid Pressure correction model is obtained using numerical simulator with pressure distribution, according to the pressure correction model to the pressure of the fluid Power and speed are modified, and are computed repeatedly speed and are announced with pressure distribution until meeting preset iterated conditional obtains fluid flowing Numerical result.
Wherein, pressure correction mould is obtained using numerical simulator according to the VELOCITY DISTRIBUTION of the fluid and pressure distribution Type.In order to improve accuracy, numerical simulator must embody the anisotropy of eddy stress.Navier- can preferably be solved Stokes (NS) equation can be used turbulence model and solve average flow field when analyzing Turbulent Flow Field, can also be using sub- grid Stress model solves large scale eddy, can be with the direct numerical solution of full size.When solution Navier-Stokes equation When, Reynolds stress transformation model may be selected and dissipative shock wave model is closed.Wherein, can presuppose initial VELOCITY DISTRIBUTION and Pressure distribution.
It is modified according to pressure and speed of the pressure correction model to the fluid.According to the initial velocity of hypothesis Distribution and pressure distribution discrete equation constant term and further acquire pressure correction equation in evaluation simulation model, with right VELOCITY DISTRIBUTION and pressure distribution are modified, and further solve digital simulation mould according to revised VELOCITY DISTRIBUTION and pressure distribution Other discretization equations in type.
It computes repeatedly speed and announces and be distributed with pressure until meeting preset iterated conditional obtains the numerical value meter of fluid flowing Calculate result.Realize that the VELOCITY DISTRIBUTION in positive fluid flow is constantly repaired in iterative calculation and pressure distribution is straight by repeating to correct and calculate The numerical result of fluid flowing is obtained to preset iterated conditional is met, includes geometrical model in the numerical result The flow regimes parameters such as the flowing velocity of middle different location fluid, according to fluid flow in the flow regime parameter of each point can analyze Obtain the fluid state in deposited tube.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity, Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer equipment, specifically, computer is set It is standby for example can for personal computer, laptop computer, cellular phone, camera phone, smart phone, personal digital assistant, Media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or these equipment In any equipment combination.
Computer equipment specifically includes memory, processor and storage on a memory simultaneously in a typical example The computer program that can be run on a processor is realized when the processor executes described program and is held as described above by client Capable method, alternatively, the processor realizes the method executed as described above by server when executing described program.
Below with reference to Figure 15, it illustrates the structures for the computer equipment 600 for being suitable for being used to realize the embodiment of the present application to show It is intended to.
As shown in figure 15, computer equipment 600 includes central processing unit (CPU) 601, can be read-only according to being stored in Program in memory (ROM) 602 is loaded into random access storage device (RAM) from storage section 608) program in 603 And execute various work appropriate and processing.In RAM603, also it is stored with system 600 and operates required various program sum numbers According to.CPU601, ROM602 and RAM603 are connected with each other by bus 604.Input/output (I/O) interface 605 is also connected to Bus 604.
I/O interface 605 is connected to lower component: the importation 606 including keyboard, mouse etc.;It is penetrated including such as cathode The output par, c 607 of spool (CRT), liquid crystal ultramagnifier (LCD) etc. and loudspeaker etc.;Storage section 608 including hard disk etc.; And including such as LAN card, the communications portion 609 of the grid interface card of modem etc..Communications portion 609 via such as because The grid of spy's net executes communication process.Driver 610 is also connected to I/O interface 606 as needed.Detachable media 611, such as Disk, CD, magneto-optic disk, semiconductor memory etc. are mounted on as needed on driver 610, in order to read from thereon Computer program be mounted as needed such as storage section 608.
Particularly, according to an embodiment of the invention, may be implemented as computer above with reference to the process of flow chart description Software program.For example, the embodiment of the present invention includes a kind of computer program product comprising be tangibly embodied in machine readable Computer program on medium, the computer program include the program code for method shown in execution flow chart.At this In the embodiment of sample, which can be downloaded and installed from grid by communications portion 609, and/or from removable Medium 611 is unloaded to be mounted.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data. The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM), Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this The function of each unit can be realized in the same or multiple software and or hardware when application.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want There is also other identical elements in the process, method of element, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product. Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by Task is executed by the connected remote processing devices of grid communications.In a distributed computing environment, program module can be with In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal Replacement, improvement etc., should be included within the scope of the claims of this application.

Claims (10)

1. a kind of deposited tube flow field analysis method characterized by comprising
Obtain the geometrical model of fluid in the deposited tube formed there are the round tube of sedimentary;
Define and determine the Reynolds number of the fluid;
Numerical value calculating and flow field analysis are carried out to the geometrical model according to the Reynolds number.
2. flow field analysis method according to claim 1, which is characterized in that the geometrical model of the fluid is horizontally disposed Column, the cross section profile of the column includes the circular arc of bottom end opening and the water that is separately connected with the both ends of the opening The straight line of flat setting, the straight line determine at a distance from the circular arc according to the thickness of the sedimentary.
3. flow field analysis method according to claim 1, which is characterized in that according to the Reynolds number to the geometrical model Numerical value calculating is carried out to specifically include:
Grid dividing is carried out to the geometrical model of the fluid;
Geometrical model after grid dividing is arranged according to Reynolds number obtains analysis model;
Numerical value calculating is carried out to the analysis model.
4. flow field analysis method according to claim 3, which is characterized in that carry out grid to the geometrical model of the fluid Division specifically includes:
Selected at least one trellis-type;
Geometrical model is divided into multiple regions or grid dividing is carried out to geometrical model by solving elliptic differential equation method;
The grid of boundary layer region is encrypted.
5. flow field analysis method according to claim 1, which is characterized in that after the setting grid dividing according to Reynolds number Geometrical model obtain analysis model and specifically include:
The fluid flow state parameter of the fluid is determined according to the sectional dimension of the geometrical model of the fluid and Reynolds number;
The fluid flow state parameter is set on the geometrical model and constraint condition obtains the analysis model.
6. flow field analysis method according to claim 3, which is characterized in that constraint condition is arranged on the geometrical model It specifically includes:
Boundary condition type is set;
Fluid inlet, outlet and the flow direction of fluid are set;
It is arranged without slip boundary condition.
7. flow field analysis method according to claim 1, which is characterized in that the flow field analysis includes average flow field analysis It is analyzed with transient flow field.
8. a kind of deposited tube flow field analysis system characterized by comprising
Geometrical model unit, for obtaining, there are the geometrical models of fluid in the deposited tube of the round tube of sedimentary formation;
Parameter determination unit, for defining and determining the Reynolds number of the fluid;
Analytical unit is solved, for carrying out numerical value calculating and flow field analysis to the geometrical model according to the Reynolds number.
9. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage Computer program, which is characterized in that
The processor is realized when executing described program such as any one of claim 1-7 the method.
10. a kind of computer-readable medium, is stored thereon with computer program, which is characterized in that
It realizes when the program is executed by processor such as any one of claim 1-7 the method.
CN201910437839.4A 2019-05-24 2019-05-24 Sediment tube flow field analysis method and system Active CN110232222B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910437839.4A CN110232222B (en) 2019-05-24 2019-05-24 Sediment tube flow field analysis method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910437839.4A CN110232222B (en) 2019-05-24 2019-05-24 Sediment tube flow field analysis method and system

Publications (2)

Publication Number Publication Date
CN110232222A true CN110232222A (en) 2019-09-13
CN110232222B CN110232222B (en) 2020-12-08

Family

ID=67861553

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910437839.4A Active CN110232222B (en) 2019-05-24 2019-05-24 Sediment tube flow field analysis method and system

Country Status (1)

Country Link
CN (1) CN110232222B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113111610A (en) * 2021-05-10 2021-07-13 中国空气动力研究与发展中心计算空气动力研究所 Novel sub-lattice scale model establishing method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102257407A (en) * 2008-12-18 2011-11-23 埃克森美孚上游研究公司 Sediment transport by fully developed turbulent flows
CN107261714A (en) * 2017-07-31 2017-10-20 盛志方 A kind of fine particle thermophoretic deposition dust removal method and its system
CN109726495A (en) * 2019-01-07 2019-05-07 暨南大学 A method of simulation and prediction Kandelia candelforest deposit total phosphorus concentration variation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102257407A (en) * 2008-12-18 2011-11-23 埃克森美孚上游研究公司 Sediment transport by fully developed turbulent flows
CN107261714A (en) * 2017-07-31 2017-10-20 盛志方 A kind of fine particle thermophoretic deposition dust removal method and its system
CN109726495A (en) * 2019-01-07 2019-05-07 暨南大学 A method of simulation and prediction Kandelia candelforest deposit total phosphorus concentration variation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YANZHI WANG, ET AL.: "Large eddy simulation of particle deposition and resuspension in turbulent duct flows", 《ADVANCED POWDER TECHNOLOGY》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113111610A (en) * 2021-05-10 2021-07-13 中国空气动力研究与发展中心计算空气动力研究所 Novel sub-lattice scale model establishing method
CN113111610B (en) * 2021-05-10 2022-10-14 中国空气动力研究与发展中心计算空气动力研究所 Sub-lattice scale model establishing method

Also Published As

Publication number Publication date
CN110232222B (en) 2020-12-08

Similar Documents

Publication Publication Date Title
Wan et al. Direct numerical simulation of particulate flow via multigrid FEM techniques and the fictitious boundary method
Miller et al. Globular dynamics: A connected particle system for animating viscous fluids
Monaghan Smoothed particle hydrodynamic simulations of shear flow
CN105069184B (en) A kind of stirred autoclave analogy method based on immersed Boundary Method
Zhang et al. An optimized Eulerian–Lagrangian method for two-phase flow with coarse particles: Implementation in open-source field operation and manipulation, verification, and validation
Wen et al. Lattice-type-dependent momentum-exchange method for moving boundaries
CN102436550B (en) Self-adaptive simulative method of dam break flood on complex border and actual landform
Dietzel et al. Numerical calculation of flow resistance for agglomerates with different morphology by the Lattice–Boltzmann Method
US20240160817A1 (en) Computer system for simulating physical processes using surface algorithm
Fang et al. Calculations of nonsubmerged groin flow in a shallow open channel by large-eddy simulation
CN109783935B (en) Implementation method for improving splash fluid stability based on ISPH
Ramos et al. Flat vs. curved rigid-lid LES computations of an open-channel confluence
Apsley et al. Bed-load sediment transport on large slopes: model formulation and implementation within a RANS solver
Abergo et al. Aerodynamic shape optimization based on discrete adjoint and RBF
Song et al. Skin friction reduction characteristics of nonsmooth surfaces inspired by the shapes of barchan dunes
CN109635405A (en) A kind of Multiphase Flow transient state calculation method and system based on space conservation
Winchenbach et al. Semi-analytic boundary handling below particle resolution for smoothed particle hydrodynamics
Chiang et al. Finite volume analysis of spiral motion in a rectangular lid‐driven cavity
CN110232222A (en) Deposited tube flow field analysis method and system
Naghian et al. Numerical simulation of turbulent flows using a least squares based meshless method
Zhang et al. Coarse particle-laden flows and energy dissipation in inclined hydraulic conveying pipes
Chai et al. Improving the accuracy of fictitious domain method using indicator function from volume intersection
Romanus et al. An immersed boundary-lattice Boltzmann framework for fully resolved simulations of non-spherical particle settling in unbounded domain
CN111931306A (en) Method and system for regulating and controlling interaction process of main vortex pair based on influence of auxiliary vortex pair
Yokoi Numerical method for interaction among multi-particle, fluid and arbitrary shape structure

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