CN110232222A - Deposited tube flow field analysis method and system - Google Patents
Deposited tube flow field analysis method and system Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power 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
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.
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)
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)
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 |
-
2019
- 2019-05-24 CN CN201910437839.4A patent/CN110232222B/en active Active
Patent Citations (3)
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)
Title |
---|
YANZHI WANG, ET AL.: "Large eddy simulation of particle deposition and resuspension in turbulent duct flows", 《ADVANCED POWDER TECHNOLOGY》 * |
Cited By (2)
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 |