CN106530411A - GPGPU-rendering-based simulation method and simulation system for 3D sea water - Google Patents
GPGPU-rendering-based simulation method and simulation system for 3D sea water Download PDFInfo
- Publication number
- CN106530411A CN106530411A CN201611013190.6A CN201611013190A CN106530411A CN 106530411 A CN106530411 A CN 106530411A CN 201611013190 A CN201611013190 A CN 201611013190A CN 106530411 A CN106530411 A CN 106530411A
- Authority
- CN
- China
- Prior art keywords
- gpgpu
- wave
- seawater
- sea
- computing
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a general-purpose-graphics-processing-unit (GPGPU)-rendering-based simulation method and simulation system for 3D sea water. The simulation method comprises: a CPU is used for initialization and an ocean-surface grid is generated dynamically; parallel computing is carried out on the ocean-surface grid by using the GPGPU; and the GPGPU executes a 3D Kochin algorithm to obtain a marine navigation track of a ship and calculates a mixing effect of sea water and a coastline. According to the method provided by the application, all computing and rendering are arranged in the GPGPU and sea water module computing and rendering form an integrated flow in the GPGPU, thereby reducing communication between an internal memory and a video memory and thus realizing the GPGPU performances fully and improving the simulation efficiency. On the basis of the method, rendering, reflecting, and refraction of sea water and waves as well as good coordination of the visual effect reality sense of the ship track and simulation real-time performance are realized.
Description
Technical field
The application is related to data processing field, more particularly to a kind of emulation mode of the 3D seawater rendered based on GPGPU and
Analogue system.
Background technology
A difficult problem showing of present video card when seawater is emulated, due to its reflection, refraction, metamorphosis etc. so that it is extra large
The emulation of water is extremely difficult.Prior art its fidelity when seawater is emulated is not high enough, it is impossible to realize the real-time simulation in big marine site,
Than relatively low, operational efficiency is low etc. for the heaving of the sea fidelity of wave.As simulation algorithm is not ideal enough, ship rides the sea
Without obvious wave effect, or wave effect is realized by special efficacy, and emulation sense is weaker.It is not carried out similar helicopter
The ripple effect that sea occurs when sea is hanging.And cannot be cross-platform, can only run on single platform.
The content of the invention
In view of this, this application provides the emulation mode and analogue system of a kind of 3D seawater rendered based on GPGPU, with
Solve the problems, such as that prior art fidelity is inadequate.
Specifically, the application is achieved by the following technical solution：
The invention provides a kind of emulation mode of the 3D seawater rendered based on GPGPU, the emulation mode includes：
Initialized by central processor CPU, dynamic generates sea grid；
The parallelization computing of the sea grid is carried out by general-purpose computations graphic process unit GPGPU；
3D Kochin algorithms are performed by the GPGPU and obtains the track that ship rides the sea, and seawater is calculated with sea
The mixed effect of water front.
The present invention is also provided based on a kind of analogue system of the 3D seawater rendered based on GPGPU, and the analogue system includes：
Sea mess generation unit, for being initialized by central processor CPU, dynamic generates sea grid；
Parallelization arithmetic element, for carrying out the parallel of the sea grid by general-purpose computations graphic process unit GPGPU
Change computing；
Track and mixed effect computing unit, perform 3D Kochin algorithms acquisition ship by the GPGPU and at sea navigate
Capable track, and calculate the mixed effect of seawater and coastline.
The embodiment of the present invention, by all of calculating and renders to be arranged in GPGPU and carries out, and seawater model is calculated and rendered
An overall flow is constituted in GPGPU, the communication between internal memory and video memory is reduced, can so be given full play to GPGPU performances,
Simulation efficiency is substantially improved.Method by more than, reaches seawater, wave and renders, reflect, reflect, hull track it is visual
Change the good coordination of the effect sense of reality and emulation real-time.
Description of the drawings
Fig. 1 is a kind of emulation mode of the 3D seawater rendered based on GPGPU shown in one exemplary embodiment of the application
Flow chart；
Fig. 2 is a kind of analogue system of the 3D seawater rendered based on GPGPU shown in one exemplary embodiment of the application
Structure chart.
Specific embodiment
Here in detail exemplary embodiment will be illustrated, its example is illustrated in the accompanying drawings.Explained below is related to
During accompanying drawing, unless otherwise indicated, the same numbers in different accompanying drawings represent same or analogous key element.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistent with the application.Conversely, they be only with as appended by
The example of consistent apparatus and method in terms of some described in detail in claims, the application.
It is the purpose only merely for description specific embodiment in term used in this application, and is not intended to be limiting the application.
" one kind ", " described " and " being somebody's turn to do " of singulative used in the application and appended claims is also intended to include majority
Form, unless context clearly shows that other implications.It is also understood that term "and/or" used herein is referred to and is wrapped
Containing one or more associated any or all possible combinations for listing project.
It will be appreciated that though various information, but this may be described using term first, second, third, etc. in the application
A little information should not necessarily be limited by these terms.These terms are only for same type of information is distinguished from each other out.For example, without departing from
In the case of the application scope, the first information can also be referred to as the second information, and similarly, the second information can also be referred to as
One information.Depending on linguistic context, word as used in this " if " can be construed to " ... when " or " when ...
When " or " in response to determining ".
A kind of emulation of the 3D seawater rendered based on GPGPU being illustrated in figure 1 shown in one exemplary embodiment of the application
The flow chart of method, the emulation mode include：
Step S101, is initialized by central processor CPU, and dynamic generates sea grid.
In embodiments of the present invention, analogue system carries out initial work by central processor CPU first, applies being used for
Emulation internal memory, and dynamic generate sea grid, wherein, initialization is included but is not limited to：Initialization unified calculation equipment frame
Structure CUDA (Compute Unified Device Architecture) equipment, application video memory, establishment grid sponge, establishment
VBO, LOD mess generation.
Step S102, carries out the parallelization computing of the sea grid by general-purpose computations graphic process unit GPGPU.
In embodiments of the present invention, analogue system dynamic generate sea grid after, in GPGPU (General
Purpose Graphics Processing Unit, general-purpose computations graphic process unit) in carry out the parallelization of the sea grid
Computing, wherein, parallelization computing is included but is not limited to：Calculated by the innovatory algorithm Gerstner waves equations based on FFT
The wave disturbance field of feature time point, the height field by Fourier transform and statistics wave pattern calculating wave generation, calculating are more
The generation of any wave is realized in the height field modeling of individual FFT magnanimity side slope vector wave, calculating Fourier amplitudes wave field, calculating is refined can
The simulation spray, foam and splashing are realized than Matrix Formula.The result of above-mentioned parallelization computing will be exported in video memory.
Wherein, each concurrent operation is specifically included：
(1), the wave of the innovatory algorithm Gerstner waves equation calculation features time points by based on FFT is disturbed
Dynamic field, including：
Fluctuation superposition effect of the hypothesis time t in any point (x, y) of seawater surface, the following institute of computing formula
Show：
X=x0-(k/κ)Asin(k·x0-wt)
Y=Acos (k x0-wt)
The horizontal vector that what wherein parameter k was represented be wave vector, be a ripple.What κ was represented is magnitude, its computing formula
It is as follows：
κ=2 π/λ
What wherein λ was represented is the length of ripple, and what w was represented is the frequency relative to wave vector, wave be it is yo-yo, and
Gerstner ripples are a single sine wave, it is impossible to which satisfaction realizes the demand of seawater, is folded below by iterative formula
Add up formula as follows：
kiRepresent be wave vector set, AiWhat is represented is amplitude, wiWhat is represented is frequency, φiWhat is represented is the rank of ripple
Section.The fluctuating behavior of its medium wave is by frequency wiDetermine, wherein w0Formula it is as follows：
What wherein T was represented be in seawater sometime, wiFrequency formula it is as follows：
WhereinBe quantify seawater dispersive surface [[]] represents is to round part, what w (k) was represented is the color of seawater
Scattered relation.Its calculation expression is as follows：
w2(k)=gk (1+k2L2)
What wherein g was represented is gravitational constant, the unit length of the ripple that L is represented.By the calculating wave disturbance of above formula
It is formed.
(2) it is, described that the height field that wave is produced is calculated by Fourier transform and statistics wave pattern, including：
The computing formula of the height field is as follows：
What wherein t was represented is the time, and what x=(x, z) was represented is the horizontal level of seawater grid, and k is a bivector,
K=(κx,κz), wherein κxAnd κzFormula be expressed as follows respectively shown in：
κx=2 π n/Lx
κz=2 π m/Lz
Parameter m, n are integers, and its scope is respectively：
- N/2≤n≤N/2,
- M/2≤m≤M/2,
The discrete point of height field is generated using fft algorithm process：
X=(nLx/N,mLz/M)。
Above formula completes the calculating of statistics wave pattern height field, next solves the Karen Phillips POP side of modified version
Journey solves height field：
Wherein, L=V2/ g, what parameter V was represented are the speed of wind, and g is gravitational constant,It is the direction of wind, what A was represented is
Digital constantThe movement for being perpendicular to wind direction for representing, wherein parameter l≤L.
(3), the height field modeling for calculating multiple FFT magnanimity side slope vector waves, including：
The accounting equation of the height field modeling is as follows：
What parameter, Δ x was represented be horizontal separation 2d it is vectorial.
(4) it is, described to calculate the generation that Fourier amplitudes wave field realizes any wave, including：
Wave in reality is all that random milli is irregular, is also required to do which any ripple life in ocean waves simulation algorithm
Into, next proceed to solve Fourier transform equation：
Wherein parameter ξrAnd ξiWhat is represented is random numbers of Gaussian distribution, and random numbers of Gaussian distribution tends to follow test
Accurate sea wave data, certain other can use random-number distribution.Solving equation is next proceeded to, in Fu of time t
Leaf amplitude wave field realizes that equation is as follows：
WhereinThe rolling of the ripple of expression calculates still to the right, finally the height of wave to the left.
The calculating of formula is just completed by above equation, the calculating for being applied to wave can just produce random wave height true to nature
Degree.
(5), the calculating Jacobian matrix formula realizes the simulation spray, foam and splashing, including：
The overlapping region of its wave is judged by the ripple of Jacobian matrix signal overlap：
J (x)=JxxJyy-JxyJyx
Above-mentioned equation decomposition is obtained into below equation：
Wherein D=(Dx,Dy), continue solving equation as follows：
Solve for convenience, the characteristic value and characteristic vector of above-mentioned equation matrix are expressed as into formula：
Wherein J-And J+It is two eigenvalue matrix,WithCorresponding respectively is normal orthogonal characteristic vector.Wherein
J=J-J+, it is most fast to calculate characteristic value and eigenvectors matrix using equation below, and formula is as follows：
Wherein shown in being calculated as follows of characteristic value：
Wherein：
Step S103, performs 3D Kochin algorithms by the GPGPU and obtains the track that ship rides the sea, and count
Calculate the mixed effect of seawater and coastline.
In embodiments of the present invention, after concurrent operation is carried out to sea grid, you can to the afloat navigation of ship
Track is calculated, and its circular is as follows：Assume that ship is travelled in the x-direction with speed v, the free surface wave edge of generation
The direction at an angle with x-axis to propagate, surface wave height is expressed as：
Wherein k0=g/v2, g is acceleration of gravity, and R is curve curvature radius, and A (α) represents the freedom of ship parameter attribute
Surface wave spectrum, this can be tried to achieve by Kochin equations.Current around ship are produced by the source that the central plane of ship is distributed, source
Intensity be directly proportional to the local gradient of hull.The expression formula of source strength is：
Wherein f (x, z) is the equation for representing hull parameters.For the tail in far field, in the equation above A (α) expression
Formula is as follows：
In above formula, H (k0, α) and for Kochin equations：
The model of wherein hull be function f (x, z) can with the calculating of simplified model, then by above formula recursion go out ζ (x,
Y), so as to solving 3D Kelvin equations.
It is pointed out that hull is in the marine spray and ship trajectory for navigating by water and having splashing, and hull can be with sea
Wave fluctuates, and its code is realized as follows：
WakeGeneratorParameters parameters；// produce waves splash about structure
Parameters.sprayEffects=true；// whether produce
Parameters.bowSprayOffset=100.0；The cheap value that // track produces
Parameters.beamWidth=20.0；// track width
Parameters.length=100.0；The length of // track
// produce hull track
WakeGenerator*ship=new WakeGenerator (ocean, parameters)；
Need to update per frame in a program, renewal function is as follows：
ship->Update(Vector3(shipX,shipY,shipZ),shipDirection,shipVelocity,
now)；
Parameter shipX, shipY, shipZ represent the position of hull, and shipDirection represents hull navigation direction,
Whether shipVelocity hull headways, now render.
Meanwhile, the particle trajectory that hull is navigated by water in the seawater is rendered by GPGPU dynamics, renders code as follows：
Vec4 color=texColor*mix (vec4 (fogColor, 1.0), lightColor, fogFactor) *
decay*transparency；
Wherein texColor is by function vec4 texColor=texture (particleTexture, tex)；
Arrive, wherein particleTexture is particle grain, and tex is gl_PointCoord point coordinates, and fogColor is mist effect face
Color, lightColor are light colors, and fogFactor is mist effect coefficient, and decay is applique coefficient, and it is by code float
Decay=clamp (exp (- 2.0*elapsed) * 2.0*sin (3.0*elapsed), 0.0,1.0)；Get, elased
Expression time, what transparency was represented are coefficients of transparency.It is the built-in mixed function of language that wherein mix is GLSL.clamp
It is GLSL interpolating functions, texture is texture shan.
The embodiment of the present invention, by all of calculating and renders to be arranged in GPGPU and carries out, and seawater model is calculated and rendered
An overall flow is constituted in GPGPU, the communication between internal memory and video memory is reduced, can so be given full play to GPGPU performances,
Simulation efficiency is substantially improved.Method by more than, reaches seawater, wave and renders, reflect, reflect, hull track it is visual
Change the good coordination of the effect sense of reality and emulation real-time.
Used as an alternative embodiment of the present invention, the emulation mode is further comprising the steps of：
The reflecting effect and refraction effect of seawater are rendered by the GPGPU.
In embodiments of the present invention, reflecting effect and refraction effect, the reflecting effect and folding are also included in real seawater
Penetrate effect and the realization of Fresnel algorithm is run by GPGPU generally.
Used as an alternative embodiment of the present invention, the emulation mode is further comprising the steps of：
The floating thing of seawater surface is rendered by the GPGPU using 3D applique principles.
In embodiments of the present invention, various floating things are included toward contact in seawater, such as：Greasy dirt etc., can be made by GPGPU
The floating thing of seawater surface is rendered with 3D applique principles.
A kind of emulation of the 3D seawater rendered based on GPGPU being illustrated in figure 2 shown in one exemplary embodiment of the application
The structure chart of system, the analogue system, including：
Sea mess generation unit 201, for being initialized by central processor CPU, dynamic generates sea grid.
In embodiments of the present invention, sea mess generation unit 201 carries out initial chemical industry by central processor CPU first
Make, apply for the internal memory of emulation being used for, and dynamic generates sea grid, wherein, initialization is included but is not limited to：Initialization unification
Computing device framework CUDA (Compute Unified Device Architecture) equipment, application video memory, establishment grid sea
Continuous, establishment VBO, LOD mess generation.
Parallelization arithmetic element 202, for carrying out the sea grid simultaneously by general-purpose computations graphic process unit GPGPU
Rowization computing.
In embodiments of the present invention, after dynamic generates sea grid, parallelization arithmetic element 202 is controlled analogue system
GPGPU (General Purpose Graphics Processing Unit, general-purpose computations graphic process unit) carries out the sea
The parallelization computing of grid, wherein, parallelization computing is included but is not limited to：By the innovatory algorithm Gerstner based on FFT
The wave disturbance field of waves equation calculation features time points, by Fourier transform and statistics wave pattern calculate wave produce
Height field, the height field modeling for calculating multiple FFT magnanimity side slope vector waves, calculating Fourier amplitudes wave field realize any wave
Generate, calculate the Jacobian matrix formula realization simulation spray, foam and splashing.The result of above-mentioned parallelization computing will be exported
In video memory.
Wherein, each concurrent operation includes：
(1), the wave of the innovatory algorithm Gerstner waves equation calculation features time points by based on FFT is disturbed
Dynamic field, including：
Fluctuation superposition effect of the hypothesis time t in any point (x, y) of seawater surface, the following institute of computing formula
Show：
X=x0-(k/κ)Asin(k·x0-wt)
Y=Acos (k x0-wt)
The horizontal vector that what wherein parameter k was represented be wave vector, be a ripple.What κ was represented is magnitude, its computing formula
It is as follows：
κ=2 π/λ
What wherein λ was represented is the length of ripple, and what w was represented is the frequency relative to wave vector, wave be it is yo-yo, and
Gerstner ripples are a single sine wave, it is impossible to which satisfaction realizes the demand of seawater, is folded below by iterative formula
Add up formula as follows：
kiRepresent be wave vector set, AiWhat is represented is amplitude, wiWhat is represented is frequency, φiWhat is represented is the rank of ripple
Section.The fluctuating behavior of its medium wave is by frequency wiDetermine, wherein w0Formula it is as follows：
What wherein T was represented be in seawater sometime, wiFrequency formula it is as follows：
WhereinBe quantify seawater dispersive surface [[]] represents is to round part, what w (k) was represented is the color of seawater
Scattered relation.Its calculation expression is as follows：
w2(k)=gk (1+k2L2)
What wherein g was represented is gravitational constant, the unit length of the ripple that L is represented.By the calculating wave disturbance of above formula
It is formed.
(2) it is, described that the height field that wave is produced is calculated by Fourier transform and statistics wave pattern, including：
The computing formula of the height field is as follows：
What wherein t was represented is the time, and what x=(x, z) was represented is the horizontal level of seawater grid, and k is a bivector,
K=(κx,κz), wherein κxAnd κzFormula be expressed as follows respectively shown in：
κx=2 π n/Lx
κz=2 π m/Lz
Parameter m, n are integers, and its scope is respectively：
- N/2≤n≤N/2,
- M/2≤m≤M/2,
The discrete point of height field is generated using fft algorithm process：
X=(nLx/N,mLz/M)。
Above formula completes the calculating of statistics wave pattern height field, next solves the Karen Phillips POP side of modified version
Journey solves height field：
Wherein, L=V2/ g, what parameter V was represented are the speed of wind, and g is gravitational constant,It is the direction of wind, what A was represented is
Digital constantThe movement for being perpendicular to wind direction for representing, wherein parameter l≤L.
(3), the height field modeling for calculating multiple FFT magnanimity side slope vector waves, including：
The accounting equation of the height field modeling is as follows：
What parameter, Δ x was represented be horizontal separation 2d it is vectorial.
(4) it is, described to calculate the generation that Fourier amplitudes wave field realizes any wave, including：
Wave in reality is all that random milli is irregular, is also required to do which any ripple life in ocean waves simulation algorithm
Into, next proceed to solve Fourier transform equation：
Wherein parameter ξrAnd ξiWhat is represented is random numbers of Gaussian distribution, and random numbers of Gaussian distribution tends to follow test
Accurate sea wave data, certain other can use random-number distribution.Solving equation is next proceeded to, in Fu of time t
Leaf amplitude wave field realizes that equation is as follows：
WhereinThe rolling of the ripple of expression calculates still to the right, finally the height of wave to the left.
The calculating of formula is just completed by above equation, the calculating for being applied to wave can just produce random wave height true to nature
Degree.
(5), the calculating Jacobian matrix formula realizes the simulation spray, foam and splashing, including：
The overlapping region of its wave is judged by the ripple of Jacobian matrix signal overlap：
J (x)=JxxJyy-JxyJyx
Above-mentioned equation decomposition is obtained into below equation：
Wherein D=(Dx,Dy), continue solving equation as follows：
Solve for convenience, the characteristic value and characteristic vector of above-mentioned equation matrix are expressed as into formula：
Wherein J-And J+It is two eigenvalue matrix,WithCorresponding respectively is normal orthogonal characteristic vector.Wherein
J=J-J+, it is most fast to calculate characteristic value and eigenvectors matrix using equation below, and formula is as follows：
Wherein shown in being calculated as follows of characteristic value：
Wherein：
Track and mixed effect computing unit 203, perform 3D Kochin algorithms by the GPGPU and obtain ship in sea
The track of upper navigation, and calculate the mixed effect of seawater and coastline.
In embodiments of the present invention, after concurrent operation is carried out to sea grid, track and mixed effect computing unit
203rd, the afloat ship trajectory of ship is calculated, its circular is as follows：Assume ship with speed v in the x-direction
Traveling, the free surface wave of generation are propagated along the direction at an angle with x-axis, and surface wave height is expressed as：
Wherein k0=g/v2, g is acceleration of gravity, and R is curve curvature radius, and A (α) represents the freedom of ship parameter attribute
Surface wave spectrum, this can be tried to achieve by Kochin equations.Current around ship are produced by the source that the central plane of ship is distributed, source
Intensity be directly proportional to the local gradient of hull.The expression formula of source strength is：
Wherein f (x, z) is the equation for representing hull parameters.For the tail in far field, in the equation above A (α) expression
Formula is as follows：
In above formula, H (k0, α) and for Kochin equations：
The model of wherein hull be function f (x, z) can with the calculating of simplified model, then by above formula recursion go out ζ (x,
Y), so as to solving 3D Kelvin equations.
It is pointed out that hull is in the marine spray and ship trajectory for navigating by water and having splashing, and hull can be with sea
Wave fluctuates, and its code is realized as follows：
WakeGeneratorParameters parameters；// produce waves splash about structure
Parameters.sprayEffects=true；// whether produce
Parameters.bowSprayOffset=100.0；The cheap value that // track produces
Parameters.beamWidth=20.0；// track width
Parameters.length=100.0；The length of // track
// produce hull track
WakeGenerator*ship=new WakeGenerator (ocean, parameters)；
Need to update per frame in a program, renewal function is as follows：
ship->Update(Vector3(shipX,shipY,shipZ),shipDirection,shipVelocity,
now)；
Parameter shipX, shipY, shipZ represent the position of hull, and shipDirection represents hull navigation direction,
Whether shipVelocity hull headways, now render.
Meanwhile, the particle trajectory that hull is navigated by water in the seawater is rendered by GPGPU dynamics, renders code as follows：
Vec4 color=texColor*mix (vec4 (fogColor, 1.0), lightColor, fogFactor) *
decay*transparency；
Wherein texColor is by function vec4 texColor=texture (particleTexture, tex)；
Arrive, wherein particleTexture is particle grain, and tex is gl_PointCoord point coordinates, and fogColor is mist effect face
Color, lightColor are light colors, and fogFactor is mist effect coefficient, and decay is applique coefficient, and it is by code float
Decay=clamp (exp (- 2.0*elapsed) * 2.0*sin (3.0*elapsed), 0.0,1.0)；Get, elased
Expression time, what transparency was represented are coefficients of transparency.It is the built-in mixed function of language that wherein mix is GLSL.clamp
It is GLSL interpolating functions, texture is texture shan.
The embodiment of the present invention, by all of calculating and renders to be arranged in GPGPU and carries out, and seawater model is calculated and rendered
An overall flow is constituted in GPGPU, the communication between internal memory and video memory is reduced, can so be given full play to GPGPU performances,
Simulation efficiency is substantially improved.Method by more than, reaches seawater, wave and renders, reflect, reflect, hull track it is visual
Change the good coordination of the effect sense of reality and emulation real-time.
Used as an alternative embodiment of the present invention, the analogue system also includes：
Reflection and refraction effect rendering unit, for the reflecting effect and refraction effect of seawater are rendered by the GPGPU.
In embodiments of the present invention, reflecting effect and refraction effect, the reflecting effect and folding are also included in real seawater
Penetrate effect and the realization of Fresnel algorithm is run by GPGPU generally.
Used as another alternative embodiment of the invention, the analogue system also includes：
Floating thing rendering unit, for rendering the floating thing of seawater surface by the GPGPU using 3D applique principles.
In embodiments of the present invention, various floating things are included toward contact in seawater, such as：Greasy dirt etc., can be made by GPGPU
The floating thing of seawater surface is rendered with 3D applique principles.
Device embodiment described above is only schematic, wherein the unit as separating component explanation can
To be or may not be physically separate, as the part that unit shows can be or may not be physics list
Unit, you can local to be located at one, or can also be distributed on multiple NEs.Which is selected according to the actual needs can
In some or all of module realizing the purpose of application scheme.Those of ordinary skill in the art are not paying creative labor
In the case of dynamic, you can to understand and implement.
The preferred embodiment of the application is the foregoing is only, not to limit the application, all essences in the application
Within god and principle, any modification, equivalent substitution and improvements done etc. are should be included within the scope of the application protection.
Claims (10)
1. a kind of emulation mode of the 3D seawater rendered based on GPGPU, it is characterised in that the emulation mode includes：
Initialized by central processor CPU, dynamic generates sea grid；
The parallelization computing of the sea grid is carried out by general-purpose computations graphic process unit GPGPU；
3D Kochin algorithms are performed by the GPGPU and obtains the track that ship rides the sea, and calculate seawater and coastline
Mixed effect.
2. emulation mode as claimed in claim 1, it is characterised in that the emulation mode also includes：
The reflecting effect and refraction effect of seawater are rendered by the GPGPU.
3. emulation mode as claimed in claim 1, it is characterised in that the emulation mode also includes：
The floating thing of seawater surface is rendered by the GPGPU using 3D applique principles.
4. the emulation mode as described in any one of claims 1 to 3, it is characterised in that initialization includes：Initialization unified calculation
Equipment framework CUDA equipment, application video memory, establishment grid sponge, establishment VBO, LOD mess generation.
5. the emulation mode as described in any one of claims 1 to 3, it is characterised in that the parallelization computing, including：Pass through
Based on the wave disturbance field of the innovatory algorithm Gerstner waves equation calculation features time points of FFT, by Fourier transform
The height field of wave generation is calculated, the height field modeling of multiple FFT magnanimity side slope vector waves is calculated, is calculated Fu with statistics wave pattern
Vertical leaf amplitude wave field is realized the generation of any wave, calculates the Jacobian matrix formula realization simulation spray, foam and splashing.
6. a kind of analogue system of the 3D seawater rendered based on GPGPU, it is characterised in that the management system includes：
Sea mess generation unit, for being initialized by central processor CPU, dynamic generates sea grid；
Parallelization arithmetic element, the parallelization for the sea grid is carried out by general-purpose computations graphic process unit GPGPU are transported
Calculate；
Track and mixed effect computing unit, perform 3D Kochin algorithms by the GPGPU and obtain what ship rode the sea
Track, and calculate the mixed effect of seawater and coastline.
7. analogue system as claimed in claim 6, it is characterised in that the analogue system also includes：
Reflection and refraction effect rendering unit, for the reflecting effect and refraction effect of seawater are rendered by the GPGPU.
8. analogue system as claimed in claim 6, it is characterised in that the analogue system also includes：
Floating thing rendering unit, for rendering the floating thing of seawater surface by the GPGPU using 3D applique principles.
9. the analogue system as described in any one of claim 6～8, it is characterised in that the initialization includes：Initialization unification
Computing device framework CUDA equipment, application video memory, establishment grid sponge, establishment VBO, LOD mess generation.
10. the analogue system as described in any one of claim 6～8, it is characterised in that the parallelization computing includes：Pass through
Based on the wave disturbance field of the innovatory algorithm Gerstner waves equation calculation features time points of FFT, by Fourier transform
The height field of wave generation is calculated, the height field modeling of multiple FFT magnanimity side slope vector waves is calculated, is calculated Fu with statistics wave pattern
Vertical leaf amplitude wave field is realized the generation of any wave, calculates the Jacobian matrix formula realization simulation spray, foam and splashing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611013190.6A CN106530411A (en) | 2016-11-17 | 2016-11-17 | GPGPU-rendering-based simulation method and simulation system for 3D sea water |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611013190.6A CN106530411A (en) | 2016-11-17 | 2016-11-17 | GPGPU-rendering-based simulation method and simulation system for 3D sea water |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106530411A true CN106530411A (en) | 2017-03-22 |
Family
ID=58352444
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611013190.6A Pending CN106530411A (en) | 2016-11-17 | 2016-11-17 | GPGPU-rendering-based simulation method and simulation system for 3D sea water |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106530411A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416831A (en) * | 2018-03-14 | 2018-08-17 | 苏州蜗牛数字科技股份有限公司 | A kind of method of real-time simulation spray and foam |
| CN111950157A (en) * | 2020-08-17 | 2020-11-17 | 北京世冠金洋科技发展有限公司 | Simulation output method and device and electronic equipment |
| CN113345066A (en) * | 2021-05-26 | 2021-09-03 | 网易（杭州）网络有限公司 | Method, device, equipment and computer-readable storage medium for rendering sea waves |
| CN114792360A (en) * | 2022-05-07 | 2022-07-26 | 北京领为军融科技有限公司 | Dynamic ocean simulation method based on fast Fourier transform |
| WO2023093231A1 (en) * | 2021-11-26 | 2023-06-01 | 腾讯科技（深圳）有限公司 | Virtual scene-based rendering method and apparatus, electronic device, computer-readable storage medium and computer program product |
| CN117292038A (en) * | 2023-10-23 | 2023-12-26 | 上海臻图信息技术有限公司 | Rendering method, system, equipment and storage medium for sea surface model |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102789650A (en) * | 2012-07-19 | 2012-11-21 | 中国科学院软件研究所 | Sea surface track parallel simulation method based on particle system |
| CN103577656A (en) * | 2013-11-25 | 2014-02-12 | 哈尔滨工业大学 | Three-dimensional dynamic simulation method for water outlet process of submarine-launched missiles |
-
2016
- 2016-11-17 CN CN201611013190.6A patent/CN106530411A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102789650A (en) * | 2012-07-19 | 2012-11-21 | 中国科学院软件研究所 | Sea surface track parallel simulation method based on particle system |
| CN103577656A (en) * | 2013-11-25 | 2014-02-12 | 哈尔滨工业大学 | Three-dimensional dynamic simulation method for water outlet process of submarine-launched missiles |
Non-Patent Citations (3)
| Title |
|---|
| 於文: "基于OSG的航海模拟器视景系统关键技术研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 * |
| 李波: "复杂环境下的海面实时建模与仿真研究", 《中国博士学位论文全文数据库 (信息科技辑)》 * |
| 陈丽宁 等: "表面波法海浪绘制综述", 《计算机工程与设计》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416831A (en) * | 2018-03-14 | 2018-08-17 | 苏州蜗牛数字科技股份有限公司 | A kind of method of real-time simulation spray and foam |
| CN111950157A (en) * | 2020-08-17 | 2020-11-17 | 北京世冠金洋科技发展有限公司 | Simulation output method and device and electronic equipment |
| CN111950157B (en) * | 2020-08-17 | 2023-08-08 | 北京世冠金洋科技发展有限公司 | Simulation output method and device and electronic equipment |
| CN113345066A (en) * | 2021-05-26 | 2021-09-03 | 网易（杭州）网络有限公司 | Method, device, equipment and computer-readable storage medium for rendering sea waves |
| CN113345066B (en) * | 2021-05-26 | 2022-07-22 | 网易（杭州）网络有限公司 | Method, device, equipment and computer-readable storage medium for rendering sea waves |
| WO2023093231A1 (en) * | 2021-11-26 | 2023-06-01 | 腾讯科技（深圳）有限公司 | Virtual scene-based rendering method and apparatus, electronic device, computer-readable storage medium and computer program product |
| CN114792360A (en) * | 2022-05-07 | 2022-07-26 | 北京领为军融科技有限公司 | Dynamic ocean simulation method based on fast Fourier transform |
| CN117292038A (en) * | 2023-10-23 | 2023-12-26 | 上海臻图信息技术有限公司 | Rendering method, system, equipment and storage medium for sea surface model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106530411A (en) | GPGPU-rendering-based simulation method and simulation system for 3D sea water | |
| Zhang et al. | A study on prediction of ship maneuvering in regular waves | |
| US20210232733A1 (en) | Systems and methods for computer simulation of detailed waves for large-scale water simulation | |
| Kim et al. | Estimation of added resistance and ship speed loss in a seaway | |
| Liu et al. | Prediction of added resistance of ships in waves | |
| Park et al. | Experimental and numerical studies on added resistance of ship in oblique sea conditions | |
| Chen et al. | Numerical ship navigation based on weather and ocean simulation | |
| Zhao et al. | Nonlinear sloshing in rectangular tanks under forced excitation | |
| Ueng et al. | A ship motion simulation system | |
| Kim et al. | Development of a numerical simulation tool for efficient and robust prediction of ship resistance | |
| CN101706972A (en) | Three-dimensional visualized algorithm for maritime oil spill | |
| CN103984793A (en) | Method for predicting movement of FLNG (floating liquefied natural gas) by considering rocking influence of liquid cabin | |
| Wu et al. | Nonlinear wave forces on large-scale submerged tunnel element | |
| Pinkster et al. | A real-time simulation technique for ship-ship and ship-port interactions | |
| Wang et al. | A Gardner evolution equation for topographic Rossby waves and its mechanical analysis | |
| Deusen et al. | The elements of nature: interactive and realistic techniques | |
| Liu et al. | Time domain simulation of nonlinear ship motions using an impulse response function method | |
| Xu et al. | Hydrodynamic analysis of a tethered underwater robot with control equipment subjected to cnoidal waves | |
| Xu et al. | Experimental study on inertial hydrodynamic behaviors of a complex remotely operated vehicle | |
| Hassan et al. | Numerical solution of the rotating shallow water flows with topography using the fractional steps method | |
| Lo | Numerical simulation of hydrodynamic interaction produced during the overtaking and the head-on encounter process of two ships | |
| US20140081612A1 (en) | Adaptive Physics Engine for Rendering Rigid Body and or Soft Body Physics for Virtual Objects in Contact with Voxelized Fluid | |
| Isaacson | Ice mass motions near an offshore structure | |
| Li | Application of 3D virtual ocean image simulation in ship vision | |
| Liu et al. | On nonlinear simulation methods and tools for evaluating the performance of ships and offshore structures in waves |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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: 20170322 |