CN109960836A - A kind of method, apparatus and terminal device of turbulent wake emulation - Google Patents
A kind of method, apparatus and terminal device of turbulent wake emulation Download PDFInfo
- Publication number
- CN109960836A CN109960836A CN201711429916.9A CN201711429916A CN109960836A CN 109960836 A CN109960836 A CN 109960836A CN 201711429916 A CN201711429916 A CN 201711429916A CN 109960836 A CN109960836 A CN 109960836A
- Authority
- CN
- China
- Prior art keywords
- particle
- sph
- vortex filament
- sampling
- turbulent wake
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000002245 particle Substances 0.000 claims abstract description 104
- 238000005070 sampling Methods 0.000 claims abstract description 45
- 239000012530 fluid Substances 0.000 claims abstract description 41
- 230000004888 barrier function Effects 0.000 claims abstract description 23
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 23
- 230000000903 blocking effect Effects 0.000 claims abstract description 19
- 238000004590 computer program Methods 0.000 claims description 19
- 238000003860 storage Methods 0.000 claims description 12
- 238000003786 synthesis reaction Methods 0.000 claims description 9
- 238000009792 diffusion process Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000004907 flux Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 abstract description 5
- 238000000926 separation method Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 241000283973 Oryctolagus cuniculus Species 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000001869 rapid Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 235000020069 metaxa Nutrition 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention is suitable for computer graphics techniques field, provides the method, apparatus and terminal device of a kind of turbulent wake emulation, this method comprises: choosing the SPH particle of blocking surfaces as burble point；Vortex filament particle sampler is carried out at the burble point, forms vortex filament；The vortex filament of sampling is discharged into SPH molecular SPH fluids and forms turbulent wake；Processing is diffused to the turbulent wake of formation.The present invention, based on this separation point sampling vortex filament, can effectively solve the Sampling of the problem of being difficult treatment fluid and barrier interaction under the prior art and vortex filament by choosing burble point in blocking surfaces.
Description
Technical field
A kind of emulated the invention belongs to computer graphics techniques field more particularly to turbulent wake method, apparatus and
Terminal device.
Background technique
In computer graphics, fluid simulation has developed for many years, and is widely used.
In computer graphics, using the method progress fluid simulation of Eulerian mesh, they make by Foster and Metaxas
Liquid is tracked in simulation process with label particles, particle density is transferred on grid in rendering, to utilize volume drawing
Technology realizes fluid simulation.
But the method for Eulerian mesh is difficult the interaction between treatment fluid and barrier at present, especially for
Barrier with small―gap suture and wedge angle, needs to be continuously updated and adapts to mobile object with tessellated mesh.
Summary of the invention
In view of this, the embodiment of the invention provides the method, apparatus and terminal device of a kind of emulation of turbulent wake, with solution
The problem of being certainly difficult the interaction between treatment fluid and barrier in the prior art.
The first aspect of the embodiment of the present invention provides a kind of method of turbulent wake emulation, comprising:
The SPH particle of blocking surfaces is chosen as burble point；
Vortex filament particle sampler is carried out at the burble point, forms vortex filament；
The vortex filament of sampling is discharged into SPH molecular SPH fluids and forms turbulent wake；
Processing is diffused to the turbulent wake of formation.
The second aspect of the embodiment of the present invention provides a kind of device of turbulent wake emulation characterized by comprising
Burble point selection unit, for choosing the SPH particle of blocking surfaces as burble point；
Vortex filament sampling unit forms vortex filament for carrying out vortex filament particle sampler at the burble point；
Wake flow synthesis unit forms turbulent flow for the vortex filament of sampling to be discharged into SPH molecular SPH fluids
Wake flow；
Wake flow diffusion unit, for being diffused processing to the turbulent wake of formation.
The third aspect of the embodiment of the present invention provides a kind of terminal device of turbulent wake emulation, including memory, place
The computer program managing device and storage in the memory and can running on the processor, which is characterized in that described
The step of processor realizes method as described in relation to the first aspect when executing the computer program.
The fourth aspect of the embodiment of the present invention provides a kind of computer readable storage medium, the computer-readable storage
Media storage has computer program, and the step of method as described in relation to the first aspect is realized when the computer program is executed by processor
Suddenly.
Existing beneficial effect is the embodiment of the present invention compared with prior art:
The present invention carries out vortex filament sampling by choosing burble point in blocking surfaces, and the vortex filament of sampling is discharged into SPH
In fluid, turbulent wake is generated, then processing is diffused to the turbulent wake.Therefore, the present invention chooses in blocking surfaces
Burble point, based on this separation point sampling vortex filament, can effectively solve be difficult the problem for the treatment of fluid and barrier interact and
The Sampling of vortex filament.
Detailed description of the invention
It to describe the technical solutions in the embodiments of the present invention more clearly, below will be to embodiment or description of the prior art
Needed in attached drawing be briefly described, it should be apparent that, the accompanying drawings in the following description is only of the invention some
Embodiment for those of ordinary skill in the art without any creative labor, can also be according to these
Attached drawing obtains other attached drawings.
Fig. 1 is the implementation process schematic diagram of turbulent wake emulation mode provided in an embodiment of the present invention；
Fig. 2 is the flow diagram of the step S101 of turbulent wake emulation mode provided in an embodiment of the present invention；
Fig. 3 A, Fig. 3 B, Fig. 4 and Fig. 5 are that rabbit barrier provided in an embodiment of the present invention forms flame in SPH fluid
Vortex filament samples schematic diagram during stream；
Fig. 6 is the flow diagram of the step S104 of turbulent wake emulation mode provided in an embodiment of the present invention；
Fig. 7 is the schematic diagram of the device of turbulent wake emulation provided in an embodiment of the present invention；
Fig. 8 is the schematic diagram of the terminal device of turbulent wake emulation provided in an embodiment of the present invention.
Specific embodiment
In being described below, for illustration and not for limitation, the tool of such as particular system structure, technology etc is proposed
Body details, to understand thoroughly the embodiment of the present invention.However, it will be clear to one skilled in the art that there is no these specific
The present invention also may be implemented in the other embodiments of details.In other situations, it omits to well-known system, device, electricity
The detailed description of road and method, in case unnecessary details interferes description of the invention.
In order to illustrate technical solutions according to the invention, the following is a description of specific embodiments.
As shown in Figure 1, being the schematic flow diagram for a kind of method that the embodiment of the present invention provides turbulent wake emulation.The turbulent flow
The method of wake flow emulation suitable for fluid simulation to the situation of turbulent wake analog simulation, by device that turbulent wake emulates Lai
It executes.The device of turbulent wake emulation is configured at the terminal device of turbulent wake emulation, can be by hardware and/or software Lai real
It is existing.As shown in Figure 1, the method for turbulent wake emulation can include: step S101 to S104.
S101 chooses the SPH particle of blocking surfaces as burble point.
Smoothed Particle Hydrodynamics Method (Smoothed Particle Hydrodynamics, SPH) is widely used in
The fluid simulation of field of Computer Graphics.Basin is divided into one group of discrete particle by SPH, then applies a smooth kernel function
Carry out the dynamic characteristic of approximate particle, and with these mobile particles of its local velocity.Due to this smoothing technique, the space of particle is led
Number can be obtained from the interaction of particle adjacent thereto.
In embodiments of the present invention, using weak compressible SPH method as basic framework；Barrier uses the surface triangulation network
Lattice model indicates；And fluid particles, then it is indicated using SPH particle.
In the embodiment of the present invention, first by being fixed on the SPH particle on its surface, i.e. side in blocking surfaces uniform sampling
Boundary's particle, SPH fluid can then be transformed into fluid SPH particle with the reciprocation of barrier and interact work between boundary particle at this time
With.Therefore, SPH particle of the present invention by choosing blocking surfaces is solved as burble point, and based on the separation point sampling vortex filament
The Sampling of the problem of being difficult treatment fluid and barrier interaction under the prior art of having determined and vortex filament
Optionally, as shown in Fig. 2, step S101 includes S1011 and S1012.
S1011, calculates the Rayleigh coefficient of boundary particle, and the boundary particle is the SPH particle of blocking surfaces.
Boundary particle is the SPH particle of blocking surfaces, can be by being fixed on its surface in blocking surfaces uniform sampling
SPH particle obtain.
In embodiments of the present invention, the Rayleigh coefficient of boundary particle is calculated using the method under the prior art, and the present invention is not
It repeats again.
S1012, if the Rayleigh coefficient of the boundary particle is greater than preset threshold, using the boundary particle as burble point.
If the Rayleigh coefficient of the boundary particle is greater than preset threshold, using the boundary particle as burble point.Due to rapids
Stream wake flow results from barrier rear, and therefore, this method calculates the Rayleigh coefficient at boundary particle, and it is big then to choose Rayleigh coefficient
In preset threshold boundary particle as burble point, vortex filament is sampled at this burble point, synthesize for subsequent turbulent wake based on
It calculates.The embodiment of the present invention reduces sampling by selecting the boundary particle for meeting preset condition in boundary particle as burble point
While quantity improves simulation efficiency, required real simulation effect is also reached.It should be noted that preset threshold is herein
Empirical value, those skilled in the art, which can according to need, to be configured, and the present invention is not particularly limited this.
S102 carries out vortex filament particle sampler at the burble point, forms vortex filament.
As previously described, because turbulent wake results from barrier rear, therefore, this method calculates the Rayleigh system of boundary particle
Number then chooses Rayleigh coefficient and is greater than the boundary particle of preset threshold as burble point.According to barrier at this burble point
Configuration sampling one encloses vortex filament particle, forms vortex filament, is used for subsequent turbulent wake composite calulation, can obtain emulation more true to nature
Effect.
Optionally, step S102 includes: to make perpendicular by the burble point, and barrier synthesizes region in turbulent wake
Interior profile and the perpendicular form boundary line, in one circle vortex filament particle of boundary line up-sampling, form vortex filament.
Turbulent wake results from barrier rear, but not all areas at barrier rear can generate flame
Stream, for example, in the barrier and part that is not disposed in SPH fluid, can't generate turbulent wake.Therefore, turbulent flow will be generated
The region of wake flow, referred to as turbulent wake synthesize region.Under normal conditions, the β plane for being parallel to horizontal plane with one is partitioned into
Turbulent wake synthesizes region.Specifically, β plane is below partly belongs to turbulent wake synthesis region.
In an embodiment of the invention, using rabbit as barrier, in SPH fluid for progress vortex filament sampling.Such as
Shown in Fig. 3 A- Fig. 5, vortex filament samples schematic diagram during forming turbulent wake in SPH fluid for rabbit barrier.Wherein, scheme
3A is schematic side view, and Fig. 3 B is schematic top plan view.Fig. 4's and Fig. 5 is respectively sampling vortex filament particle and the amplified figure of vortex filament.
In Fig. 3 A and Fig. 3 B, pi, pjIt is two burble points of rabbit blocking surfaces, it should be noted that other separation
Point is not shown in Fig. 3 A and Fig. 3 B.αsIt is the perpendicular by burble point, in Fig. 3 A and Fig. 3 B, burble point pi, pjIt is proper
It is located at the same perpendicular α wellsIt is interior.β plane is partitioned into turbulent wake synthesis region, below to partly belong to turbulent wake conjunction
At region.
In conjunction with shown in Fig. 3 A and Fig. 4, rabbit barrier is formed in β plane profile below and the perpendicular to have a common boundary
Line, in one circle vortex filament particle of boundary line up-samplingAlso referred to as vortex filament vertex, nsIt is number of samples, this
One circle vortex filament particle obviously includes burble point piAnd pj.Referring back to shown in Fig. 5, this circle vortex filament particleAccording to
It is secondary to be connected, form vortex filament γs.Specifically, the sampling interval of vortex filament particle, the present invention is not specifically limited.In the same sampling
Between, vortex filament sampling is carried out in all burble points.For the vortex filament sampling of other burble points, and take identical vortex filament sampling side
Formula, details are not described herein again.
The vortex filament of sampling is discharged into SPH molecular SPH fluids and forms turbulent wake by S103.
After vortex filament sampling, the vortex filament of sampling is discharged respectively according to the preset time interval according to its precedence generated
To among SPH molecular SPH fluids, turbulent wake is formed.For example, in an embodiment of the invention, being arranged every 10
Time step carries out vortex filament sampling, and the vortex filament of sampling is discharged into SPH fluid, makes it along with one cocurrent flow of SPH fluid
It is dynamic.
SPH particle cuts down under your rule (Biot-SavartLaw) acts in Biot-Sa and moves, i.e., for each production
Raw vortex filament unit, all generates a vorticity to each SPH particle in entire fluid mass, and SPH particle completes all whirlpools
After the superposition of silk effect, toroidal vortex flow wake flow, calculation are generated are as follows:
Wherein, the v (x) indicates that position coordinates are the fluid velocity of the SPH particle at x, and m is the quantity of the vortex filament,
ΓiIt is the vortex flux of i-th of vortex filament, γi(t) be i-th of vortex filament t-th of vortex filament particle position coordinates.
S104 is diffused processing to the turbulent wake of formation.
Turbulent wake under true environment can be spread apart after its formation with the time to dribble, and therefore, the present invention is real
It applies in example, after synthesizing turbulent wake, processing is diffused to the turbulent wake of formation, to obtain emulation effect more true to nature
Fruit.
Optionally, as shown in fig. 6, in embodiments of the present invention, step S104 includes S1041 and S1042.
S1041 distinguishes surface SPH particle and inside SPH particle using the SPH particle of color field flag SPH fluid.
Wherein, using the SPH particle of color field flag SPH fluid:
Wherein, ciIndicate the color field scalar value of particle, mjIt is the quality of j-th of neighborhood particle of i-th of particle, ρjIt is
The density of j-th of neighborhood particle of i-th of particle, WijIt is the smoothing kernel function in SPH method.
By color field flag, SPH particle in SPH fluid is marked, distinguishes surface SPH particle and inside SPH
Son.Surface SPH particle has surface tension, and inside SPH particle is without surface tension.
S1042 introduces surface tension model to surface SPH particle, carries out at discretization to the turbulent wake of formation
Reason.
Wherein, for surface SPH particle, continuum surface force can be calculated by the following formula:
Wherein, σ is surface tension coefficient, niIt is surface particle normal vector.
Since surface tension is a kind of internal power in Free Surface, liquid surface area may make to minimize, smooth liquid
Gap between face can prevent the untrue phenomenon generation on surface during fluid motion, and therefore, the present invention passes through the company of introducing
Continued face tension model carries out sliding-model control using SPH method, realizes turbulent wake diffusion effect true to nature.
It should be understood that the size of the serial number of each step is not meant that the order of the execution order in above-described embodiment, each process
Execution sequence should be determined by its function and internal logic, the implementation process without coping with the embodiment of the present invention constitutes any limit
It is fixed.
It is a kind of schematic block diagram of the device of turbulent wake emulation provided in an embodiment of the present invention referring to Fig. 7, Fig. 7.The rapids
It is not described in detail place in the Installation practice of stream wake flow emulation, refers to method Fig. 1-6 implementation of turbulent wake emulation
Example.The device of the turbulent wake emulation of the present embodiment includes: burble point selection unit 71, vortex filament sampling unit 72, wake flow synthesis
Unit 73 and wake flow diffusion unit 74.
Wherein, burble point selection unit 71, for choosing the SPH particle of blocking surfaces as burble point；
Vortex filament sampling unit 72 forms vortex filament for carrying out vortex filament particle sampler at the burble point；
Wake flow synthesis unit 73 forms rapids for the vortex filament of sampling to be discharged into SPH molecular SPH fluids
Flow wake flow；
Wake flow diffusion unit 74, for being diffused processing to the turbulent wake of formation.
Optionally, burble point selection unit 71 includes coefficients calculation block 711 and burble point determining module 712.
Wherein, coefficients calculation block 711, for calculating the Rayleigh coefficient of boundary particle, the boundary particle is barrier
The SPH particle on surface.
Burble point determining module 712, if the Rayleigh coefficient for the boundary particle is greater than preset threshold, by the boundary
Particle is as burble point.
Optionally, wake flow diffusion unit 74 includes that color field mark module 741 and surface tension introduce module 742.
Wherein, color field mark module 741 distinguishes surface SPH for using the SPH particle of color field flag SPH fluid
Particle and inside SPH particle.
Surface tension introduces module 742, for introducing surface tension model to surface SPH particle, to the rapids of formation
Stream wake flow is diffused processing.
Fig. 8 is the schematic diagram of the terminal device for the turbulent wake emulation that one embodiment of the invention provides.As shown in figure 8, should
The terminal device 8 of embodiment includes: processor 80, memory 81 and is stored in the memory 81 and can be in the processing
The computer program 82 run on device 80, such as the program of turbulent wake emulation.The processor 80 executes the computer journey
The step in the embodiment of the method for above-mentioned each turbulent wake emulation is realized when sequence 82, such as step S101 shown in FIG. 1 is extremely
S104.Alternatively, the processor 80 realizes each unit/module in above-mentioned each Installation practice when executing the computer program 82
Function, such as the function of unit 71 to 74 shown in Fig. 7.
Illustratively, the computer program 82 can be divided into one or more units/modules, it is one or
Multiple units/modules are stored in the memory 81, and are executed by the processor 80, to complete the present invention.Described one
A or multiple units/modules can be the series of computation machine program instruction section that can complete specific function, which is used for
Implementation procedure of the computer program 82 in the terminal device 8 that the turbulent wake emulates is described.For example, the computer
It is (empty that program 82 can be divided into burble point selection unit, vortex filament sampling unit, wake flow synthesis unit and wake flow diffusion unit
Unit in quasi- device), each unit concrete function is as follows:
Burble point selection unit, for choosing the SPH particle of blocking surfaces as burble point；
Vortex filament sampling unit forms vortex filament for carrying out vortex filament particle sampler at the burble point；
Wake flow synthesis unit forms turbulent flow for the vortex filament of sampling to be discharged into SPH molecular SPH fluids
Wake flow；
Wake flow diffusion unit, for being diffused processing to the turbulent wake of formation.
The terminal device 8 of the turbulent wake emulation can be desktop PC, notebook, palm PC and cloud clothes
Business device etc. calculates equipment.The terminal device may include, but be not limited only to, processor 80, memory 81.Those skilled in the art
It is appreciated that Fig. 8 is only the example of terminal device 8, the restriction to terminal device 8 is not constituted, may include than illustrating more
More or less component perhaps combines certain components or different components, such as the terminal device can also include input
Output equipment, network access equipment, bus etc..
Alleged processor 80 can be central processing unit (Central Processing Unit, CPU), can also be
Other general processors, digital signal processor (Digital Signal Processor, DSP), specific integrated circuit
(Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor is also possible to any conventional processor
Deng.
The memory 81 can be the internal storage unit of the terminal device 8, such as the hard disk or interior of terminal device 8
It deposits.The memory 81 is also possible to the External memory equipment of the terminal device 8, such as be equipped on the terminal device 8
Plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital (Secure Digital, SD) card dodge
Deposit card (Flash Card) etc..Further, the memory 81 can also both include the storage inside list of the terminal device 8
Member also includes External memory equipment.The memory 81 is for storing needed for the computer program and the terminal device 8
Other programs and data.The memory 81 can be also used for temporarily storing the data that has exported or will export.
It is apparent to those skilled in the art that for convenience of description and succinctly, only with above-mentioned each function
Can unit, module division progress for example, in practical application, can according to need and by above-mentioned function distribution by different
Functional unit, module are completed, i.e., the internal structure of described device is divided into different functional unit or module, more than completing
The all or part of function of description.Each functional unit in embodiment, module can integrate in one processing unit, can also
To be that each unit physically exists alone, can also be integrated in one unit with two or more units, it is above-mentioned integrated
Unit both can take the form of hardware realization, can also realize in the form of software functional units.In addition, each function list
Member, the specific name of module are also only for convenience of distinguishing each other, the protection scope being not intended to limit this application.Above system
The specific work process of middle unit, module, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In the above-described embodiments, it all emphasizes particularly on different fields to the description of each embodiment, is not described in detail or remembers in some embodiment
The part of load may refer to the associated description of other embodiments.
Those of ordinary skill in the art may be aware that list described in conjunction with the examples disclosed in the embodiments of the present disclosure
Member and algorithm steps can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually
It is implemented in hardware or software, the specific application and design constraint depending on technical solution.Professional technician
Each specific application can be used different methods to achieve the described function, but this realization is it is not considered that exceed
The scope of the present invention.
In embodiment provided by the present invention, it should be understood that disclosed device/terminal device and method, it can be with
It realizes by another way.For example, device described above/terminal device embodiment is only schematical, for example, institute
The division of module or unit is stated, only a kind of logical function partition, there may be another division manner in actual implementation, such as
Multiple units or components can be combined or can be integrated into another system, or some features can be ignored or not executed.Separately
A bit, shown or discussed mutual coupling or direct-coupling or communication connection can be through some interfaces, device
Or the INDIRECT COUPLING or communication connection of unit, it can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated module/unit be realized in the form of SFU software functional unit and as independent product sale or
In use, can store in a computer readable storage medium.Based on this understanding, the present invention realizes above-mentioned implementation
All or part of the process in example method, can also instruct relevant hardware to complete, the meter by computer program
Calculation machine program can be stored in a computer readable storage medium, the computer program when being executed by processor, it can be achieved that on
The step of stating each embodiment of the method.Wherein, the computer program includes computer program code, the computer program generation
Code can be source code form, object identification code form, executable file or certain intermediate forms etc..The computer-readable medium
It may include: any entity or device, recording medium, USB flash disk, mobile hard disk, magnetic that can carry the computer program code
Dish, CD, computer storage, read-only memory (Read-Only Memory, ROM), random access memory (Random
Access Memory, RAM), electric carrier signal, telecommunication signal and software distribution medium etc..It should be noted that the meter
The content that calculation machine readable medium includes can carry out increase and decrease appropriate according to the requirement made laws in jurisdiction with patent practice,
It such as does not include electric carrier signal and telecommunications according to legislation and patent practice, computer-readable medium in certain jurisdictions
Signal.
Embodiment described above is merely illustrative of the technical solution of the present invention, rather than its limitations；Although referring to aforementioned reality
Applying example, invention is explained in detail, those skilled in the art should understand that: it still can be to aforementioned each
Technical solution documented by embodiment is modified or equivalent replacement of some of the technical features；And these are modified
Or replacement, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution should all
It is included within protection scope of the present invention.
Claims (11)
1. a kind of method of turbulent wake emulation characterized by comprising
The SPH particle of blocking surfaces is chosen as burble point；
Vortex filament particle sampler is carried out at the burble point, forms vortex filament；
The vortex filament of sampling is discharged into SPH molecular SPH fluids and forms turbulent wake；
Processing is diffused to the turbulent wake of formation.
2. the method as described in claim 1, which is characterized in that it is described that vortex filament sampling is carried out at the burble point, form whirlpool
Silk, comprising:
Vortex filament particle is enclosed according to the configuration sampling of barrier one at the burble point, forms vortex filament.
3. method according to claim 2, which is characterized in that it is described at the burble point according to the configuration sampling of barrier
One circle vortex filament particle, forms vortex filament, comprising:
Make perpendicular, profile and the perpendicular shape of the barrier in turbulent wake synthesis region by the burble point
Vortex filament is formed in one circle vortex filament particle of boundary line up-sampling at boundary line.
4. the method according to claim 1, which is characterized in that the SPH particle for choosing blocking surfaces is made
For burble point, comprising:
The Rayleigh coefficient of boundary particle is calculated, the boundary particle is the SPH particle of blocking surfaces；
If the Rayleigh coefficient of the boundary particle is greater than preset threshold, using the boundary particle as burble point.
5. the method according to claim 1, which is characterized in that described that the vortex filament of sampling is discharged into SPH
Turbulent wake is formed in the molecular SPH fluid of grain, comprising:
The vortex filament of sampling is discharged into SPH particle group according to its precedence generated respectively according to the preset time interval
At SPH fluid in form turbulent wake.
6. method as claimed in claim 5, which is characterized in that the calculation formula of the turbulent wake are as follows:
Wherein, the v (x) indicates that position coordinates are the fluid velocity of the SPH particle at x, and m is the quantity of the vortex filament, ΓiIt is
The vortex flux of i-th of vortex filament, γi(t) be i-th of vortex filament t-th of vortex filament particle position coordinates.
7. the method according to claim 1, which is characterized in that the turbulent wake of described pair of formation expands
Dissipate processing, comprising:
Using the SPH particle of color field flag SPH fluid, surface SPH particle and inside SPH particle are distinguished；
Surface tension model is introduced for surface SPH particle, sliding-model control is carried out to the turbulent wake of formation.
8. the method for claim 7, which is characterized in that the surface tension model, are as follows:
Wherein, σ is surface tension coefficient, niIt is the normal vector of the surface SPH particle of SPH fluid, ciIndicate surface SPH particle
Color field scalar value.
9. a kind of device of turbulent wake emulation characterized by comprising
Burble point selection unit, for choosing the SPH particle of blocking surfaces as burble point；
Vortex filament sampling unit forms vortex filament for carrying out vortex filament particle sampler at the burble point；
Wake flow synthesis unit forms flame for the vortex filament of sampling to be discharged into SPH molecular SPH fluids
Stream；
Wake flow diffusion unit, for being diffused processing to the turbulent wake of formation.
10. a kind of terminal device of turbulent wake emulation, including memory, processor and storage are in the memory and can
The computer program run on the processor, which is characterized in that the processor is realized when executing the computer program
Such as the step of any one of claim 1 to 8 the method.
11. a kind of computer readable storage medium, the computer-readable recording medium storage has computer program, and feature exists
In when the computer program is executed by processor the step of any one of such as claim 1 to 8 of realization the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711429916.9A CN109960836A (en) | 2017-12-26 | 2017-12-26 | A kind of method, apparatus and terminal device of turbulent wake emulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711429916.9A CN109960836A (en) | 2017-12-26 | 2017-12-26 | A kind of method, apparatus and terminal device of turbulent wake emulation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109960836A true CN109960836A (en) | 2019-07-02 |
Family
ID=67021885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711429916.9A Pending CN109960836A (en) | 2017-12-26 | 2017-12-26 | A kind of method, apparatus and terminal device of turbulent wake emulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109960836A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110617238A (en) * | 2019-09-26 | 2019-12-27 | 浙江理工大学 | Optimization design method of centrifugal pump impeller |
| CN110929450A (en) * | 2019-11-22 | 2020-03-27 | 北京航空航天大学 | Realistic turbulence simulation method based on SPH |
| CN111824193A (en) * | 2020-07-29 | 2020-10-27 | 西南交通大学 | High-speed train wake passive control method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050109879A1 (en) * | 2003-11-21 | 2005-05-26 | Patterson Robert A. | Method and apparatus for quantum vortex implosion propulsion and species |
| CN103810332A (en) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | Helicopter ship surface wake flow influence simulation method |
| WO2014088680A2 (en) * | 2012-09-20 | 2014-06-12 | Stephen Heppe | Systems and methods for long endurance airship operations |
| CN104537175A (en) * | 2014-12-30 | 2015-04-22 | 中国科学院深圳先进技术研究院 | SPH-algorithm-based fluid simulation method and device |
-
2017
- 2017-12-26 CN CN201711429916.9A patent/CN109960836A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050109879A1 (en) * | 2003-11-21 | 2005-05-26 | Patterson Robert A. | Method and apparatus for quantum vortex implosion propulsion and species |
| WO2014088680A2 (en) * | 2012-09-20 | 2014-06-12 | Stephen Heppe | Systems and methods for long endurance airship operations |
| CN103810332A (en) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | Helicopter ship surface wake flow influence simulation method |
| CN104537175A (en) * | 2014-12-30 | 2015-04-22 | 中国科学院深圳先进技术研究院 | SPH-algorithm-based fluid simulation method and device |
Non-Patent Citations (6)
| Title |
|---|
| XIANGYUN LIAO； WEIXIN SI；BIAO XU等: "Filament-based realistic turbulent wake synthesis", 《COMPUTER ANIMATION & VIRTUAL WORLDS》 * |
| 刘彦文 * |
| 游振江,林建忠: "气固两相流叶轮机械磨损问题研究综述", 《流体机械 》 * |
| 潘永琛; 姚建伟; 李昌烽: "高速列车近尾流区湍流涡旋的数值计算及动力学讨论", 《力学季刊》 * |
| 肖凯涛等: "粒子系统在复杂风场及湍流场中的实现", 《北京交通大学学报》 * |
| 范全林; 王希麟; 张会强等: "圆湍射流拟序结构研究进展", 《力学进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110617238A (en) * | 2019-09-26 | 2019-12-27 | 浙江理工大学 | Optimization design method of centrifugal pump impeller |
| CN110617238B (en) * | 2019-09-26 | 2021-01-05 | 浙江理工大学 | Optimization design method of centrifugal pump impeller |
| CN110929450A (en) * | 2019-11-22 | 2020-03-27 | 北京航空航天大学 | Realistic turbulence simulation method based on SPH |
| CN110929450B (en) * | 2019-11-22 | 2021-09-21 | 北京航空航天大学 | Realistic turbulence simulation method based on SPH |
| CN111824193A (en) * | 2020-07-29 | 2020-10-27 | 西南交通大学 | High-speed train wake passive control method |
| CN111824193B (en) * | 2020-07-29 | 2021-10-26 | 西南交通大学 | High-speed train wake passive control method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Golovashkin et al. | Solving diffractive optics problems using graphics processing units | |
| CN110392902A (en) | Use the operation of sparse volume data | |
| Spencer et al. | Evenly spaced streamlines for surfaces: An image‐based approach | |
| CN109960836A (en) | A kind of method, apparatus and terminal device of turbulent wake emulation | |
| CN107067045A (en) | Data clustering method, device, computer-readable medium and electronic equipment | |
| Edmunds et al. | Automatic stream surface seeding: A feature centered approach | |
| CN107122243A (en) | Heterogeneous Cluster Environment and CFD computational methods for CFD simulation calculations | |
| Mirzaee et al. | Smoothness-increasing accuracy-conserving filters for discontinuous Galerkin solutions over unstructured triangular meshes | |
| Murthy et al. | Buffer merging—a powerful technique for reducing memory requirements of synchronous dataflow specifications | |
| Wu et al. | A heterogeneous platform with GPU and FPGA for power efficient high performance computing | |
| CN106776023B (en) | Task load balancing method for self-adaptive GPU unified dyeing array | |
| Xu et al. | GPU-accelerated evaluation platform for high fidelity network modeling | |
| Thulasidasan et al. | Explicit spatial scattering for load balancing in conservatively synchronized parallel discrete event simulations | |
| Chen et al. | Scalability of Algebraic Multigrid in Computer Science | |
| Zhao et al. | Flow simulation with locally-refined LBM | |
| Wang et al. | GPU‐accelerated finite‐difference time‐domain method for dielectric media based on CUDA | |
| CN110120959A (en) | Big data method for pushing, device, system, equipment and readable storage medium storing program for executing | |
| CN109670158A (en) | A kind of method and apparatus for according to information data generation content of text | |
| CN106874925A (en) | object grouping method, model training method and device | |
| CN103440395A (en) | Design system and method for realizing BlinnPhong illumination model based on lookup table | |
| Szymczak et al. | Nearly recurrent components in 3D piecewise constant vector fields | |
| Mi et al. | A novel and adaptive mesh method for arbitrary discrete fracture networks DFN simulations | |
| CN109246606A (en) | Extended method, device, terminal device and the storage medium of robot localization network | |
| US20230394768A1 (en) | System and method of generating smooth spline surface model preserving feature of physical object | |
| CN109286798A (en) | A kind of bezel locations recognition methods, system and the terminal device of projection screen |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190702 |