当前位置:首页 >> >> Simulation, Modelling and Rendering of Incompressible Fluids in Real Time

Simulation, Modelling and Rendering of Incompressible Fluids in Real Time


Simulation, Modelling and Rendering of Incompressible Fluids in Real Time
Thomas Klein, Mike Eissele, Daniel Weiskopf, Thomas Ertl Universit¨ t Stuttgart a Institut f¨ r Visualisierung und Interaktive Systeme (VIS) u Universit¨ tsstra?e 38, 70569 Stuttgart, Germany a Email: {klein,eissele,weiskopf,ertl}@informatik.uni-stuttgart.de

Abstract
In this paper we present techniques for the real-time simulation and rendering of liquids. Appropriate approximations to a full 3D simulation are applied to reduce the numerical complexity. Fluid ?ow is described by the 2D Navier-Stokes equations; wave effects are represented according to the wave equation for shallow water waves on a height-?eld. A novel mechanism for coupling ?ow and wave behaviour is introduced to ef?ciently handle transport along a ?uid ?ow and accurate wave propagation on the ?uid surface. An additional noise-based animation of detailed ?uid structures further improves the realistic appearance. Rendering makes use of fragment operations on consumer-level GPUs to allow for a physically guided representation of re?ection, refraction and Fresnel blending of light.

1 Introduction
Water is present in our everyday lives. Therefore the realism of any virtual environment is greatly enhanced by including water phenomena, and the animation and rendering of liquids is an important and challenging topic in computer graphics. In recent years much progress has been made with the physically based simulation of water effects. Most techniques are based on direct simulation methods to solve the full 3D Navier-Stokes equations for incompressible ?uids, e.g. [6]. Recent publications concerning ?uids mostly make use of a semi-Lagrangian treatment introduced to computer graphics by Stam [18]. This method avoids numerical oscillations in convection-dominated ?ows and allows large time steps. State-of-the-art animation and rendering of water using 3D Navier-Stokes is VMV 2003

presented by Enright et al. [4]. The major disadvantage of these accurate simulation methods are the required computation and rendering times. In this work we do not focus on a physically correct simulation but rather on real-time animation and rendering of realistic looking water effects. There is a large, commercially interesting market for such applications, reaching from the game industry to visual simulation applications and interactive controlling and modelling of liquids in animation tools. Typical performance should be some 5– 60 frames per second depending on the type of application, i.e., between two to three orders of magnitude faster than a full-bodied simulation. Even with the increasing performance of PC hardware in the years to come, a real-time application of a complete 3D liquid simulation seems to be impossible on PCs in the foreseeable future. Therefore we propose an approach that allows for a restricted, yet important class of water phenomena. The main idea is to reduce a full 3D Navier-Stokes simulation to an appropriate 2D approximation. A twodimensional computation is less expensive than a three-dimensional and, in particular, scales much better with increasing numerical grid resolution. In this context one class of previous and related animation approaches focuses on a reduced model representation of ?uid surfaces by means of parametric modelling [7, 14, 20]. Another approach [13], guided by oceanographic observations, is based upon the ?ltering of a noise ?eld by a wave spectrum using the fast Fourier transform. In recent work [10] an adaptive scheme for the real-time animation and display of ocean waves is proposed. Common to all above methods is that they contain no concept of the ?uid ?ow or mass transport, but only a pure description of the ?uids surface behaviour. Munich, Germany, November 19–21, 2003

Another class of techniques is based on 2D approximations of the full 3D Navier-Stokes equations. The approach presented here follows this line of reasoning and is closely related to the following two papers. Kass and Miller [12] obtain a height?eld representation of the ?uid surface by means of a linearised form of the 2D shallow water wave equations. Chen and Lobo [2] use the 2D incompressible Navier-Stokes equations to compute a 2D ?ow and directly drive a height-?eld by the pressure scalar ?eld obtained from the ?ow simulation. An important contribution of this paper is to introduce a coupling between ?ow and wave behaviour. In our system, ?ow properties can induce wave effects in order to consider, e.g., an in?ow into a lake that produces both streaming and wave effects. Using a height-?eld restricts the possible classes of phenomena that can be achieved, e.g., splashing or pouring cannot be reproduced. However, visually important effects of wave propagation are taken into account and moreover, transport along the ?ow of the liquid is considered, by means of the ?ow simulation based upon the 2D NavierStokes equations. A typical application is the simulation of the ?ow and waves for a river or a lake. Additional noise-based animation techniques are used to further improve the appearance of the ?uid surface. Noise-based animation has been found useful for approximating a wide range of natural phenomena. Stam and Fiume [19] use an FFT method to add a small scale component to turbulent wind ?elds. Also, Schneider and Westermann [16] generate a water surface by completely restricting themselves to noise in the form of fractional Brownian motion. We use this approach to add some additional detail structure to our simulated ?uid surface. Besides the animation of ?uids the visual representation of their surfaces is extremely important for the realism of a graphical application. There are many rendering algorithms that can properly simulate complex light transport, but due to their computational complexity most of these algorithms are not suitable for real-time applications. We propose a hardware-accelerated approach to allow for physics-based re?ection and refraction effects, which is faster than a related approach [11]. The computations are performed with per-pixel accuracy and exploiting the capabilities of low-cost graphics hardware. 666

2 Numerical Description of Fluid Flow
The dynamic behaviour of a viscous ?uid, like water, is completely described by the so-called NavierStokes equations. For incompressible Newtonian ?uids they state that, in every point, the following conditions must be ful?lled, ?v 1 + (v · ?) v = Re ?2 v ? ?p + f ?t ?·v =0 , , (1a) (1b)

where v is the ?uid velocity, p the pressure, f the external body forces and Re the so-called Reynolds number, which accounts for the overall behaviour of the ?ow. For a detailed derivation and description of the Navier-Stokes equations we refer to textbooks on ?uid dynamics, such as [1, 8]. This system of partial differential equations has to be solved for the velocity v and pressure p in every point of the computational domain. In general there is no analytical solution, so numerical methods have to be employed. In this work a ?nitedifferences approach based on the marker-and-cell (MAC) method by Harlow and Welch [9] is followed. Spatial derivatives are discretised using central differences on a staggered grid. But employing central differences for discretising all spatial derivatives leads to problems, like numerical oscillations for convection-dominated ?ows [1, 8]. Therefore we propose a combination of central differences and a donor-cell scheme to discretise the convective term (v · ?) v. The donorcell scheme is well-know in the CFD community, cf. [8], and has been used, in a slightly different form, for example in [21]. The qualitative behaviour of the donor-cell scheme can be compared to the semi-Lagrangian methodology [18]. Combining the donor-cell scheme and central differences has the nice property, that it provides a user-adjustable factor that can be used to ?nd a compromise between the numerical damping due to ?rst order interpolation and the second order accuracy of central differences. The development in time due to Eqns. (1) is treated by a standard projection method [3, 8]. According to that, the numerical algorithm for solving Eqns. (1) can be split into three major steps. First, a l+1 temporary velocity ?eld vtmp is computed by applying an explicit Euler method, while neglecting the

coupling to the pressure. Second, a Poisson equation for the pressure ?eld is derived and solved by treating the pressure derivatives implicitly in time and enforcing conservation of mass by means of Eqn. (1b). Third, the new velocity ?eld v l+1 is coml+1 puted by updating vtmp using the new pressure values obtained in the previous step. For a visual real-time simulation, a constant computation time per frame should be reached to facilitate a constant application ?ow. The computation time for steps one and three is of O(n), where n equals the number of grid cells. However, computation time for the second step is rather problematic. A numerical approach to the Poisson equation leads to a sparse system of linear equations which could be solved by a number of different methods. With respect to performance, iteration methods are better suited than direct methods for such sparse systems. A popular method is SOR iteration (successive over-relaxation) [22], for which the number of iteration cycles depends on the error level and is not constant. In our approach, on the other hand, only a ?xed number of steps of the iterative algorithm is computed regardless of the achieved convergence. Although this does not lead to an accurate solution, because the continuity condition is not completely satis?ed, the achieved accuracy is still suf?cient for a qualitatively correct ?uid behaviour in a visual simulation. Since the SOR iteration has a preferred direction of iterating through the solution vector, artifacts can occur if the iteration is terminated before convergence is reached. These artifacts can be reduced by alternating the direction of updates in every iteration step, which leads to the symmetric SOR (SSOR) method [22]. From experience, six to eight iteration steps are suf?cient to achieve results with suf?cient accuracy for typical visual simulations of ?uids. The ?xed number of iteration steps guarantees a constant execution time. Another issue is numerical stability. Because we use a semi-explicit time integration scheme, our method exhibits the inherent instability associated with explicit methods. The maximum time step is restricted by various factors, including the Reynolds number, the cell size, the donor-cell damping factor and last but not least the applied boundary conditions. We do not present a thorough stability analysis of the solution scheme; the interested reader may refer to the CFD literature, for example [9, 1, 8] for 666

stability conditions that have to be satis?ed. But it should be noted that due to our speci?c simpli?cations in the pressure correction step these conditions may not be very reliable. Generally it can be said that small time steps, lower Reynolds numbers, small in?ow velocities, high donor-cell damping factors and good-natured boundary conditions lead to stable solutions.

3

Shallow Water Waves

For simulating the ?uid surface we use a method proposed by Kass and Miller [12] to compute the wave propagation on a height-?eld. This method is based on the so-called shallow water equations, a simpli?ed form of the previously described NavierStokes equations. Introducing several further assumptions these equations can be reduced to a single wave equation: ? ? 2 ?2h ?2h ? h , (2) = gd + ?t2 ?x2 ?y 2 where h is the deviation of the surface from the ?uid at rest, d the water depth and g the gravitational constant. This wave equation can be ef?ciently solved using the ADI (alternating-direction implicit) method [12]. The advantage of the ADI method over a simple explicit method is its unconditional stability which results from an implicit discretisation. Moreover, computational costs are linear to the number of grid cells. Therefore a constant simulation time can be achieved for real-time applications. We introduce additional damping via a Lax-Wendroff dissipative interface as described in [1] to account for the viscosity of the ?uid.

4

Coupling Flow and Wave Behaviour

So far two distinct methods that describe some partial properties of the ?uid behaviour have been presented: Wave propagation on a height-?eld and 2D ?ow. The fundamental problem one encounters at this point is how to couple the pure 2D ?ow with the height-?eld in 3D. Considering real ?uid motion one observes that ?uid ?ow and surface waves are two, although not independent, yet loosely coupled ?uid properties. Waves have virtually no effect on the ?ow itself because they exhibit no mass transportation. The underlying ?uid ?ow instead

is—apart from other sources like wind or bodies falling into or moving through the ?uid—the main reason for the existence of surface waves. So the surface waves can be considered as a minor effect compared to the ?uid ?ow and therefore one can focus on a mapping from 2D ?ow to wave effects on the height-?eld. In [2] this problem is solved by linearly mapping pressure to height values. This approach has the drawback of a direct ?ow-to-wave relationship that does not allow for effects like superposed waves originating from other locations of the ?ow or outside sources (objects falling into the liquid, waves re?ected at boundaries, etc.). Instead we propose an indirect coupling between the 2D ?ow and the height-?eld. The solution of the ?ow simulation is utilised to generate initial disturbances in the height-?eld that are propagated by means of the wave equation. We use the local divergence of the temporary l+1 velocity ?eld vtmp generated by the Navier-Stokes solver to disturb the height-?eld. This follows a very simple heuristic: divergence indicates a change in mass, i.e., mass has been moved to or from neighbouring cells. This mass divergence must result in a change of the surface height. If there is a mass surplus (negative divergence) in the cell, the ?uid surface has to rise; vice versa if the divergence is positive, the surface height must fall. The disturbed height-?eld before the wave propagation step is then given by
l+1 ? hl+1 = hl + fh (?vtmp ) ,

takes into account that moving further away from the equilibrium state is considerably harder than moving towards the ?uid height at rest. The combination of a 2D Navier-Stokes based ?ow simulation and a wave equation driven height?eld allows us to bene?t from the advantages of both methods. Using the wave equation enables one to take advantage of a wide range of wave phenomena. This includes, e.g., wave refraction due to ground unevennesses, re?ection of waves at obstacles or superposition of waves originating from wind force or ?uid-user interactions. In contrast, the ?ow simulation can be employed for material transport, e.g., by tracing massless particles—which can be used for simulating leafs ?owing on water (cf. Fig. 4). Another application is the interaction between the ?uid and rigid bodies ?owing within it.

5 Controlling the Flow
The behaviour of the ?uid ?ow is mainly controlled by setting boundary conditions. Three different types of boundary conditions are supported: Each cell on the computational grid can be set to be an in?ow, an out?ow or an obstacle cell. These boundary conditions are applied during each simulation cycle, i.e., values for velocity and pressure are set according to the respective conditions. Setting boundary conditions can be done in two different ways. First, the simulation system itself can set obstacle cells automatically depending on the geometry of the simulation environment. The data describing the geometrical properties of the simulation environment is speci?ed in the form of a 3D scene ?le (in the current implementation a simpli?ed variant of the AutoDesk 3DS ?le format). It must contain the geometry of the scene surrounding the ?uid and one or more planar meshes that de?ne the geometries of the ?uid surfaces. When the scene is loaded, each ?uid surface gives rise to an appropriate rectangular grid whose cells are classi?ed either to be ?uid or obstacle cells. Fig. 1 illustrates how the grid is generated from a triangulated ?uid surface. Here, the cells completely covered by the original mesh are regarded as ?uid cells. The other possibility is to set the cell type interactively before or even during the simulation run. For this purpose our system provides a cell-based editor that can be used to manipulate the types of the grid cells and to set their properties, e.g., time-dependent in?ow ve666

where the function fh de?nes how the current height-?eld hl is affected by additional disturbance. ` ?2 De?ning fh (δ) = sδ/(1 + hl ) has proved to be useful because this mapping takes into account that the greater the deviation between the actual ?uid surface and the ?uid surface at rest is, the smaller is the effect of any additional disturbance. The strength of the perturbation can be controlled via the scaling factor s. Other choices for the mapping function fh are possible as well. In particular, the direction of the deviation compared to the direction of the newly introduced disturbance can be taken into account in addition to the magnitude of deviation. For example, fh (δ) = ( sδ
sδ 2 1+(hl )

for sign(δ) = sign(hl ) for sign(δ) = sign(hl )

in surface-local space, we have to transform them to world space. The so calculated re?ection vector points into a cube texture map to provide the re?ection colour according to the environment mapping approach. Unfortunately, refraction cannot be imY

Z

N V
X

Figure 1: Automatic boundary condition setup. locity. Moreover the editor allows the user to add sources for massless particles or rigid bodies.
?N

?N T T

Figure 2: Refraction space.

6 Rendering
High-quality rendering of the ?uid surface is crucial for a realistic looking simulation. Important optical phenomena related to the interaction between light and a ?uid surface have to be taken into account by a physical model. We present a solution that allows for effects of re?ection, refraction and Fresnel blending. All computations are based on the limited per-fragment operations of today’s pixel shader 1.0 compliant consumer graphics hardware such as the GeForce3. The geometry of the ?uid surface is built from the height-?eld over the uniform grid that is calculated in the simulation part. The calculation of the aforementioned effects take place in several different coordinate systems. Object space and world space are known from the standard rendering pipeline. We need two additional spaces: surface-local space [5] and refraction space. The surface-local space, also known as tangent space or texture space, is de?ned for every point on the surface such that the point lies in the origin and the normal vector points along the positive z axis. The x and y axes lie in the tangent plane and are orthogonal to each other. The setup of the refraction space is illustrated in Fig. 2. The y axis is oriented in a way to make the incoming line of sight V point along the negative y axis. Fragment operations allow to calculate the re?ection vector for a given viewing vector and surface normal. Since we apply the normals per fragment 666 plemented in a similar, straightforward way. The aforementioned class of graphics hardware imposes rather sharp restrictions onto the number and types of possible fragment operations. Due to these limitations, the refraction vector cannot be computed and used for a lookup in an environment map. To overcome this problem we propose an approximation based on the refraction space as illustrated in Fig. 2. The negative of the normal vector, ?N , represents the transmission vector for a refraction index of in?nity; the vector V represents the transmission vector for a refraction index n = 1 (no refraction). The transmission vector for a refraction index n ∈ [1, ∞) has to lie between ?N and V . This observation leads to the following qualitatively correct approximation of the transmission vector: the negative normal ?N is scaled down in its x and z components to yield the approximated transmission vector. This non-uniform scale forces the negative normal to converge to the transmission vector. The amount of the scale corresponds to the refraction index. The approximated transmission vector is used for lookup in the cubic environment map to collect the refraction colour. The advantage of this approach is the missing computation of complicated numerical terms that are replaced by a linear transformation via a combination of several matrices M = (R)(T)(Rt )(S). The matrix S transforms the normals from surface-local space to world space. The

matrix R transforms vectors from refraction space to world space. A non-uniform scale of the normal in refraction space is performed through the matrix T. The ?nal matrix M transforms the surface-local normal vectors to world-space transmission vectors. A ?uid surface that is only represented through the height-?eld from a simulation often misses to display ?ne details on the surface. A typical resolution of the computational grid is some 1282 cells for a real-time simulation on a current PC. This resolution is roughly one to two orders of magnitude smaller than a typical screen resolution. However, ?ne surface details contribute a huge amount of realism to the appearance of the ?uid. A realistic and ef?cient possibility to add detail to the surface makes use of a noise-based approach. The perfragment normals of the ?uid surface are perturbed, in order to achieve a realistic appearance. We choose a turbulence function consisting of different frequencies of Perlin noise [15]. Blending different scales of the same original noise ?eld yields the different frequencies. This method is well suited for a real-time implementation because all necessary computational steps are supported by graphics hardware. So far only a static perturbation of the ?uid height-?eld is possible which would look totally unrealistic in an animation. To overcome this problem the surface detail is animated with a time-dependent spatial offset for the different frequency patterns. The combination of these patterns is computed for each frame with changing offset values. The result is an animated random surface detail where the viewer is no longer able to recognise a single random pattern. Moreover, a computation of the overall noise texture for each frame allows us to include timedependent local distortions that only show up in a spatially limited region. For example, a small gust of wind can have a local effect on the appearance of the ?uid surface because wave ripples are introduced by wind dragging. These local effects massively contribute to a realistic appearance of the whole ?uid. The local features very often concern high frequency noise that is modulated upon the random surface detail, as illustrated in Fig. 6. A static shape of the local distortions with sharp edges would lead to an unrealistic impression. To circumvent this behaviour the intensity of the effect is faded out at the border and the shape is animated 666

using a weighting factor represented through a 2D texture. The position, size, motion and shape of local effects can easily be managed by changing this texture. An ef?cient approach to animate the shape of the effect is to blend several static shapes with different blending factors. Therefore the surface appearance can be controlled by outside effects and react to conditions of the surrounding environment. Finally, a lighting term is evaluated per fragment to modulate the refraction colour. This allows us to simulate the lighting of the ground through the ?uid. Afterwards the blending of the re?ection and refraction terms is achieved via a approximated Fresnel term F = (1 ? cos θi )4 , since an evaluation of the correct Fresnel term using only fragment operations is not possible on the targeted graphics hardware. The value θi describes the angle between the surface normal and the viewing direction. In most cases, the viewer cannot distinguish the true solution of the Fresnel term from this rough approximation, whereas the calculation of this simple approximation uses only a few fragment operations. A single cubic environment map is used for the re?ection and refraction computations. This environment map has to be updated whenever the viewer’s position is changed. Otherwise spatial distortions may become apparent if the position of the viewer is not equal to the position from where the environment map was created. Note that this dynamic environment map is generated once per frame and can be re-used for other rendering effects of a 3D engine.

7 Implementation and Results
A key element of the implementation is the partitioning into simulation and rendering parts. The implementation of the simulation part is based on C++ and makes use of the CPU only. Rendering is based on a GeForce3 GPU using OpenGL with several extensions. This separation leads to a parallelisation where the CPU calculates the simulation while the GPU renders the previous step. The rendering part is divided into four steps. The ?rst step for every frame is to update the dynamic environment map. For fast updates of dynamic textures any data transfer between the GPU and the CPU should be avoided. Since the functionality to directly render into textures is only available for some operating systems, we use a PBuffer to ren-

der into an offscreen area. In the second step the surface detail texture is updated. Two stages of the static noise texture are blended with different scaling and offset values. With two additional stages a spatially limited change of the detail texture is introduced. The ?rst of these two texture stages is setup to blend three static shapes of a local feature which are encoded in the colour channels of a texture. The following stage uses the blended shape to mask the high frequency noise texture for the local effect. The surface detail now contains perturbations of the normals for each fragment of the ?uid surface. This normal texture is used in the following rendering passes for re?ections and refraction. The re?ection and refraction rendering passes use vertex programmes to ef?ciently calculate the normals of the height-?eld. OpenGLs vertex pointers provide a fast access to the scalar height-?eld in order to evaluate central differences at the vertices. The x, y and z components in the vertex pointers hold three, shifted versions of the height?eld, such that the y components of the ?rst and third pointer and the x and z components of the second represent the height values for the direct neighbours of the current vertex. The y component of the second pointer provides the height of the vertex itself. Thus, the vertex programmes have access to all data that is required. In both programmes the normal vector, in combination with the tangent and the binormal vector, is used to calculate the appropriate transformation matrix that is then applied to the surface-detail normal in the fragment operations. For re?ections four texture shader stages are setup to perform a transformation of each detail normal from surface-local space to world coordinates. The normalization of the transformed normal is achieved by a lookup in a normalisation cube map. The texture stages also compute the re?ection vector based on the transformed normal. This re?ection vector points into the environment map to obtain the re?ection colour. A register combiner setup evaluates the Fresnel approximation based on the normalised normal and stores the result in the α portion of the frame buffer. For refractions the normals of the height-?eld are calculated identically to the re?ection case, but the transformation matrix is set to the matrix M. The texture shaders transform the normal from surfacelocal space into the approximate transmission vec666

tor with respect to world coordinates. The transmission colour is obtained by a texture lookup in the environment map. Since all texture stages are already consumed, the surface normal vector cannot be transformed. Therefore the lighting calculation has to take place in surface-local space, utilising the register combiners. Fig. 3 shows the results of different rendering effects. Fig. 3 (a) displays re?ections only. Image (b) shows refractions of the ?uid surface. The fabric on the pool ground appears distorted through the waves on the ?uid surface. The Fresnel blending in (c) clearly exhibits that, in the area close to the viewer, one can only see the refractions while further away from the viewer the re?ections play a major role. Fine details of the ?uid surface, introduced by the per fragment surface detail, provides a natural appearance of the surface, as shown in (d). In order to visualise the ?ow of the ?uid, particles and other objects can be inserted in the ?uid. In Fig. 4 the stream is visualised through leafs ?oating in the water. Fig. 5 shows an example for the coupling between ?ow and wave excitation. Fluid is streaming into the simulation grid from the top-left border. Behind the rectangular pillar the ?ow exhibits the typical vortices forming the wake behind an obstacle in the ?ow. Further material in the form of videos is available for download1 . The achieved frame rates depend on the CPU for the simulation part and on the GPU for the rendering part. Table 1 shows results for an Athlon 1200 system with a GeForce3 graphics card. All test scenes were rendered in full quality. The bottom row of the table represents the mere rendering performance without any simulation. The case where the CPU tends to be the limiting factor can be seen in the ?rst row. Here the frame rates increase only a few frames by decreasing the percentage of the displayed ?uid. Grids about the size 642 provide a good trade-off between the load of the CPU and the load of the GPU.

8

Conclusion and Future Work

We have introduces a novel technique for real-time simulation and rendering of ?uids for computer graphics applications. Combining a Navier-Stokes based 2D ?ow computation with a wave propagation technique on a height-?eld, very fast simula1

http://wwwvis.informatik.uni-stuttgart.de/?weiskopf/?uid

simulation grid size 128x128 64x64 16x16 none

screen area of displayed ?uid 100% 75% 50% 25% 0% 15 16 16 17 18 40 45 52 61 62 50 56 70 99 130 50 56 70 100 131

nology, 1989. [2] J. X. Chen and N. D. V. Lobo. Toward interactive-rate simulation of ?uids with moving obstacles using Navier-Stokes equations. Graphical models and image processing: GMIP, 57(2):107–116, 1995. [3] A. J. Chorin. Numerical solution of the Navier-Stokes equations. Mathematics of Computation, 22(104):745–762, 1968. [4] D. P. Enright, S. R. Marschner, and R. P. Fedkiw. Animation and rendering of complex water surfaces. In SIGGRAPH 2002 Conference Proceedings, pages 736–744, 2002. [5] C. Everitt. Mathematics of per-pixel lighting. http://developer.nvidia.com, 2001. NVidia Corporation. [6] N. Foster and D. Metaxas. Realistic animation of liquids. In Graphics Interface, pages 204–212, 1996. [7] A. Fournier and W. T. Reeves. A simple model of ocean waves. In SIGGRAPH 1986 Conference Proceedings, pages 75–84, 1986. [8] M. Griebel, T. Dornseifer, and T. Neunhoeffer. Numerical Simulation in Fluid Dynamics – A Practical Introduction. SIAM, Philadelphia, 1994. [9] F. H. Harlow and J. E. Welch. Numerical calculation of timedependent viscous incompressible ?ow of ?uid with free surface. The Physics of Fluids, 8(12):2182–2189, 1965. [10] D. Hinsinger, F. Neyret, and M.-P. Cani. Interactive animation of ocean waves. In Proceedings of the ACM SIGGRAPH Symposium on Computer Animation, pages 161–166, 2002. [11] K. Iwasaki, Y. Dobashi, and T. Nishita. Ef?cient rendering of optical effects within water using graphics hardware. In Proceedings of Paci?c Graphics 2001, pages 374–383, 2001. [12] M. Kass and G. Miller. Rapid, stable ?uid dynamics for computer graphics. In SIGGRAPH 1990 Conference Proceedings, pages 49–57, 1990. [13] G. A. Mastin, P. A. Watterberg, and J. F. Mareda. Fourier synthesis of ocean scenes. IEEE Computer Graphics and Applications, 7(3):16–23, 1987. [14] D. R. Peachey. Modeling waves and surf. In SIGGRAPH 1986 Conference Proceedings, pages 65–74, 1986. [15] K. Perlin. An image synthesizer. In SIGGRAPH 1985 Conference Proceedings, pages 287–296, 1985. [16] J. Schneider and R. Westermann. Towards real-time visual simulation of water surfaces. In Proceedings of the Vision, Modeling, and Visualization Conference 2001 (VMV01), pages 211–218, 2001. [17] J. Stam. Random caustics: Natural textures and wave theory revisited. In SIGGRAPH 1996 Visual Proceedings, page 151, 1996. Technical Sketch SIGGRAPH’96. [18] J. Stam. Stable ?uids. In SIGGRAPH 1999 Conference Proceedings, pages 121–128, 1999. [19] J. Stam and E. Fiume. Turbulent wind ?elds for gaseous phenomena. In SIGGRAPH 1993 Conference Proceedings, pages 369–376, 1993. [20] P. Y. Ts’o and B. A. Barsky. Modeling and rendering waves: Wave-tracing using beta-splines and re?ective and refractive texture mapping. ACM Transactions on Graphics, 6(3):191– 214, 1987. [21] G. D. Yngve, J. F. O’Brien, and J. K. Hodgins. Animating explosions. In SIGGRAPH 2000 Conference Proceedings, pages 29–36, 2000. [22] D. M. Young. Iterative Solution of Large Linear Systems. Academic Press, New York, 1971.

Table 1: Performance (fps) for an Athlon 1200 with a GeForce3 at a resolution of 800 × 600 pixels. tions of realistic looking ?uid behaviour are possible. Our method is well suited for modelling gentle motions of water in rivers or lakes. The highquality rendering of the ?uid surface supports re?ections, refractions and Fresnel blending at fragment level. We have presented a method to exploit current consumer-level graphics hardware to evaluate these effects completely within the GPU. Since our hardware requirements are extremely low, we serve a wide range of systems and provide a very good compromise of speed versus quality and portability. To further enhance the visual apperance of the rendered ?uid surface, other properties of ?uid light interaction could be considered. This includes, e.g., adding caustics to the rendering using a simpli?ed model such as by Stam [17] which is also appropriate for consumer-level graphics hardware. Additionally several enhancements are possible for the simulation part. E.g., the coupling between ?ow and wave propagation can be extended by advecting the height values with the ?ow. This would produce a more realistic transport of waves. Similarly the motion of the surface detail texture should be affected by the velocity ?eld. Texture advection could be used to implement this feature on graphics hardware. Also an approach to automatically sample the depth of the ?uid from the surrounding scene geometry should be investigated. In the rendering part the obtained depth information will help to attenuate the light through the ?uid, leading to a more sophisticated lighting model. The simulation part could also take advantage of these automatically generated depth information. For example the refraction of waves due to the real shape of the underlying scene geometry could be modelled.

References
[1] M. Abbot and D. Basco. Computational Fluid Dynamics – An Introduction For Engineers. Longman Science & Tech-

666

(a)

(b)

(c)

(d)

Figure 3: Comparison of rendering effects. Image (a) shows re?ections only, (b) refractions only, (c) Fresnel blending of re?ections and refractions, (d) all rendering effects including a noise-based surface detail.

Figure 4: Visualisation of the stream through leafs in the ?uid.

Figure 5: The typical wake, forming behind an obstacle in the ?ow.

Figure 6: Modulating the surface detail with a local effects mask. 666


更多相关文档:

...Rendering of Incompressible Fluids in Real Time_....pdf

Simulation, Modelling and Rendering of Incompressible Fluids in Real Time Thomas Klein, Mike Eissele, Daniel Weiskopf, Thomas Ertl Universit¨ t Stuttgart ...

Navier-Stokes方程组驱动的虚拟人群.pdf

arealtimepedestriansimulationmodelbasedon fluid...whichdescribe2Dincompressiblefluid,avelocityfieldwhich...balancebetweencomputationaltimeandrenderingquality....

保证时间一致性的风格化液体动画_图文.pdf

can not simulate fluid motion as in real word....3 Liquid Simulation We use the incompressible ...Real-Time rendering of cartoon smoke and clouds....

Fluids1_图文.pdf

Incompressible velocity field Incompressibility ...//www.cs.ubc.ca/~rbridson/fluidsimulation/ ? ...“Stable Fluids” [Stam 1999] ? “Real-Time ...

Multiphase CFD Simulation of a Solid Bowl Centrifuge.pdf

Multiphase CFD Simulation of a Solid Bowl ... and water as an incompressible Newtonian fluid. ...between the simulated model and the real ...

Large-Eddy Simulation of Unsteady Flow in a Mixed-Flow Pump_....pdf

ow of an incompressible ?uid, as represented in...The time increment was set so that 8000 time ...Modelling & Simulation in Engineering Volume 2014 ...

CAPABILITIES AND CHALLENGES OF CFD MULTIPHASE SIMULATION OF ....pdf

In contrast, simulation with computational fluid ...In nature there are little real flows with only...These are modeled as spherical, incompressible ...

Thermal stress simulation and interface destabilisation in ....pdf

Thermal stress simulation and interface de...Since real experiments are very expensive and di?...gas is considered as a Newtonian, incompressible ...

2011-5681 Method for CFD Simulation of Propellant Slosh in a ....pdf

to inaccurately model viscosity in real fluid. ...simulation treats all the fluids as incompressible,...The time step size is chosen because it captures...

DEM simulation of particle attrition in dilute-phase ....pdf

DEM simulation of particle attrition in dilute-...(DEM) was developed for modelling and simulating ... particles are spherical and incompressible before ...

Fluent-FSI_14.5_lect-05_CoSimulation_Convergence_图文.pdf

This may be different than the real-world ...planned FSI simulation (same time step size, etc...The fluid is (or is almost) incompressible, ...

A Real Time Transient Hydraulic Model For Leak Detection.pdf

incompressible liquids but becomes questionable for ...In particular, Real Time Transient Modelling (RTT...simulation are: al Conservation of Mass: dlAr) ...

INCOMPRESSIBLE NAVIER-STOKES CALCULATIONS FOR THE DEVELOPMENT....pdf

INCOMPRESSIBLE NAVIER-STOKES CALCULATIONS FOR THE ...ow simulation capability developed for liquid rocket...In the steady-state formulation, the time ...

Numerical Simulation of Two-Dimensional Flow and.pdf

Numerical Simulation of Two-Dimensional Flow and_...ber is incompressible. Therefore, the bulk strain...It was also called the ‘‘real three-...

Gpu simulation and rendering of volumetric effects for ....pdf

Modelling with Flow Patterns Ideally, the simulation of incompressible non-...4. Real-Time Fluid Effects Figure 6: The user interface to de?ne ...

...Culling for Efficient Fluid Flow Simulation and Rendering,....pdf

Fluid Flow Simulation and Rendering,” ATI_专业...as “early-z culling” in real-time rendering....algorithm, which simulates incompressible fluid flow...

Clouds and Atmospheric Phenomena Simulation in Real-Time 3D ....pdf

incompressible flow [1], Lattice Boltzmann Model [2] or cellular automaton ...Harris, M. J.: Real-Time Cloud Simulation and Rendering. Technical Report...

SUMMARY_图文.pdf

cloudiness contour surface in a climate simulation....rendering, and give techniques for handling time-...For incompressible 2D ?ows, stretching in one ...

Simulation of Three-Dimensional Symmetric and Asymmetric ....pdf

For time marching, a time-splitting procedure was...Simulation," Journal of Fluid Mechanics Vol. 163... \The Stability of Incompressible Attachment-Line ...

whitepaperRF4_图文.pdf

incompressible liquid using a collection of ...bitmaps that can be used at the rendering stage...simulation and the exporting phase in REALFLOW4....

更多相关标签:
网站地图

文档资料共享网 nexoncn.com copyright ©right 2010-2020。
文档资料共享网内容来自网络,如有侵犯请联系客服。email:zhit325@126.com