Realistic Animation of Fluids

Size: px
Start display at page:

Download "Realistic Animation of Fluids"

Transcription

1 Realistic Animation of Fluids p. 1/2 Realistic Animation of Fluids Nick Foster and Dimitri Metaxas

2 Realistic Animation of Fluids p. 2/2 Overview Problem Statement Previous Work Navier-Stokes Equations Explicit solving Results

3 Realistic Animation of Fluids p. 3/2 Problem Statement Simulate fluids for graphics applications Scalability Usability Easy to setup Controllable

4 Realistic Animation of Fluids p. 4/2 Related Work Non-physics based Paramtric functions Sinusoidal phase functions Physics based CFD Too expensive Not controllable Kass and Miller Shallow water Chen and Lobo Navier-Stokes in 2-D Fluid zero depth Use instantaneous pressure

5 Realistic Animation of Fluids p. 5/2 Navier-Stokes Equation Conservation of Momentum ( u 2 t + u u x + v u y + w u z = p x + g u x + µ x u y u z 2 ( 2 v v t + u v x + v v y + w v z = p y + g y + µ w t + u w x + v w y + w w z = p z + g z + µ Conservation of Mass u x + v y + w z = 0 ) ) x v y v z 2 ( 2 ) w x w y w z 2

6 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form ρ DV Dt = p + ρg + µ 2 V

7 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ ρ DV Dt = p + ρg + µ 2 V

8 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ DV Dt = p + g + µ 2 V

9 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ LHS Intertial force DV Dt = p + g + µ 2 V

10 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ LHS Intertial force DV Dt = p + g + µ 2 V RHS Surface forces due to pressure differences

11 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ LHS Intertial force DV Dt = p + g + µ 2 V RHS Surface forces due to pressure differences Gravity

12 Realistic Animation of Fluids p. 6/2 Conservation of Momentum Vector form Normalize for ρ LHS Intertial force DV Dt = p + g + µ 2 V RHS Surface forces due to pressure differences Gravity Sheer force due to viscosity

13 Realistic Animation of Fluids p. 7/2 LHS Begin with Newton s Second Law F = m dv dt

14 Realistic Animation of Fluids p. 7/2 LHS Begin with Newton s Second Law Using the Chain Rule F = m dv dt dv dt = V x dx dt + V y dy dt + V z dz dt + V t

15 Realistic Animation of Fluids p. 7/2 LHS Begin with Newton s Second Law Using the Chain Rule F = m dv dt dv dt = V x dx dt + V y dy dt + V z dz dt + V t = u V x + v V y + w V z + V t

16 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 dy F xx = µ 2 u x 2 dx

17 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

18 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

19 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

20 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

21 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

22 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2

23 Realistic Animation of Fluids p. 8/2 RHS Pressure Gravity Sheer forces F yx = µ 2 u y 2 F xx = µ 2 u x 2 In 3-D: ( 2 ) u F x = µ x u y u z 2

24 Realistic Animation of Fluids p. 9/2 Conservation of Mass Incompressible fluid Divergence should be zero

25 Realistic Animation of Fluids p. 9/2 Conservation of Mass Incompressible fluid Divergence be zero should ρ t + ρu x + ρv y + ρw z = 0

26 Realistic Animation of Fluids p. 9/2 Conservation of Mass Incompressible fluid Divergence be zero should ρ t + ρu x + ρv y + ρw z = 0 u x + v y + w z = 0

27 Realistic Animation of Fluids p. 10/2 Solving Navier-Stokes Specify obstacle/boundaries Discretize on a rough grid Solve equations Interpolate surface of fluid

28 Realistic Animation of Fluids p. 11/2 Discretization Rectangular Cartesian grid Finite differences Velocity defined on cell faces Pressure defined in cell Stability 1 > max [ u δt ] δt δt,v,w δx δy δz

29 Realistic Animation of Fluids p. 11/2 Discretization Rectangular Cartesian grid Finite differences Velocity defined on cell faces Pressure defined in cell Stability 1 > max [ u δt ] δt δt,v,w δx δy δz

30 Realistic Animation of Fluids p. 11/2 Discretization Rectangular Cartesian grid Finite differences Velocity defined on cell faces Pressure defined in cell Stability 1 > max [ u δt ] δt δt,v,w δx δy δz dy dx

31 Realistic Animation of Fluids p. 11/2 Discretization Rectangular Cartesian grid Finite differences Velocity defined on cell faces Pressure defined in cell Stability 1 > max [ u δt ] δt δt,v,w δx δy δz

32 Realistic Animation of Fluids p. 12/2 Types of Cells Empty A cell containing no particles Surface A cell containing at least 1 particle that is adajacent to an Empty cell Full A cell containing at least 1 particle and not a Surface cell

33 Realistic Animation of Fluids p. 13/2 Calculating Velocities Velocities from momentum equations u t = p t + g x u u ( 2 ) x v u y w u z + µ u x u y u z 2 Pressure from continuity equations δp is proportional to D D = u x + v y + w z

34 Realistic Animation of Fluids p. 14/2 Consistency Velocity and pressure are updated Iterate until D < ǫ δp = βd u i+1/2,j,k u i 1/2,j,k v i,j+1/2,k v i,j 1/2,k w i,j,k+1/2 w i,j,k 1/2 p i,j,k = u i+1/2,j,k + (δt/δx)δp = u i 1/2,j,k (δt/δx)δp = u i,j+1/2,k + (δt/δy)δp = u i,j 1/2,k (δt/δy)δp = u i,j,k+1/2 + (δt/δz)δp = u i,j,k 1/2 (δt/δz)δp = p i,j,k + δp

35 Realistic Animation of Fluids p. 15/2 Boundary Conditions Stationary Obstacles Non-slip Free-slip Infow and outflow Free Surfaces u 0 = 0 p 0 = p 1 w 0 = w 1 w 0 = w 1

36 Realistic Animation of Fluids p. 16/2 Free Surface Free faces set so divergence zero Surface pressure set to atmospheric pressure Cases 3-D: 64 configurations

37 Realistic Animation of Fluids p. 16/2 Free Surface Free faces set so divergence zero Surface pressure set to atmospheric pressure Cases 3-D: 64 configurations

38 Realistic Animation of Fluids p. 16/2 Free Surface Free faces set so divergence zero Surface pressure set to atmospheric pressure Cases 3-D: 64 configurations

39 Realistic Animation of Fluids p. 17/2 Recap 1. Set boundary conditions 2. Classify cells 3. Calculate velocities for Full cells 4. Calculate pressures for Full cells 5. Update velocities for Surface cells 6. Iterate

40 Realistic Animation of Fluids p. 18/2 Problem Stability condition Limited spatial resolution Computational complexity Main contribution Interpolating methods for displaying the surface 3 methods for tracking fluid positions

41 Realistic Animation of Fluids p. 19/2 Marker Particles Convect massless particles with local fluid velocity Use weighted average of 4 nearest velocities Weighting based on area Multiplied by timestep Labeling done based on particles 1 particle test

42 Realistic Animation of Fluids p. 20/2 Free Surface Particles Only keep particles at surface Rules for adding and removing particles If far apart: Insert a particle If close: remove and connect neighbors Region growing algorithm

43 Realistic Animation of Fluids p. 21/2 Height Field Scalar function h t = w u h x v h y Cell configuration trivial Need other methods Overturning wave Spray Foam

44 Realistic Animation of Fluids p. 21/2 Height Field Scalar function h t = w u h x v h y Cell configuration trivial Need other methods Overturning wave Spray Foam

45 Realistic Animation of Fluids p. 22/2 Buoyancy Rigid object Force formula f ni = p i dv i + m i g f fluid = i f ni Lagrange equations of motions M q + D q = f q + g q Do not affect water flow Collisions treated separately

46 Realistic Animation of Fluids p. 23/2 Control Inflow and Outflow Instead of atmospheric pressure, apply surface pressure history Constant Time-depedent Height of surface P P P P P P p applied (z) = A + B cos(cz wt) δt

47 Results Realistic Animation of Fluids p. 24/2

48 Results Realistic Animation of Fluids p. 25/2

49 Results Realistic Animation of Fluids p. 26/2

50 Results Realistic Animation of Fluids p. 27/2

51 Results Realistic Animation of Fluids p. 28/2

Realistic Animation of Fluids

Realistic Animation of Fluids 1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely

More information

Navier-Stokes & Flow Simulation

Navier-Stokes & Flow Simulation Last Time? Navier-Stokes & Flow Simulation Implicit Surfaces Marching Cubes/Tetras Collision Detection & Response Conservative Bounding Regions backtracking fixing Today Flow Simulations in Graphics Flow

More information

Navier-Stokes & Flow Simulation

Navier-Stokes & Flow Simulation Last Time? Navier-Stokes & Flow Simulation Pop Worksheet! Teams of 2. Hand in to Jeramey after we discuss. Sketch the first few frames of a 2D explicit Euler mass-spring simulation for a 2x3 cloth network

More information

Navier-Stokes & Flow Simulation

Navier-Stokes & Flow Simulation Last Time? Navier-Stokes & Flow Simulation Optional Reading for Last Time: Spring-Mass Systems Numerical Integration (Euler, Midpoint, Runge-Kutta) Modeling string, hair, & cloth HW2: Cloth & Fluid Simulation

More information

2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123

2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123 2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs

More information

CGT 581 G Fluids. Overview. Some terms. Some terms

CGT 581 G Fluids. Overview. Some terms. Some terms CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches

More information

Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)

Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final

More information

A Framework for Real-Time Animation of Liquid-Rigid Body Interaction

A Framework for Real-Time Animation of Liquid-Rigid Body Interaction A Framework for Real-Time Animation of Liquid-Rigid Body Interaction Melih KANDEMİR Tolga ÇAPIN Bülent ÖZGÜÇ Computer Engineering Department Bilkent University Abstract Both high-end training simulators

More information

Fluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]

Fluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11] Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t

More information

CS-184: Computer Graphics Lecture #21: Fluid Simulation II

CS-184: Computer Graphics Lecture #21: Fluid Simulation II CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves

More information

Chapter 1 - Basic Equations

Chapter 1 - Basic Equations 2.20 Marine Hydrodynamics, Fall 2017 Lecture 2 Copyright c 2017 MIT - Department of Mechanical Engineering, All rights reserved. 2.20 Marine Hydrodynamics Lecture 2 Chapter 1 - Basic Equations 1.1 Description

More information

(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more

(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon

More information

Mass-Spring Systems. Last Time?

Mass-Spring Systems. Last Time? Mass-Spring Systems Last Time? Implicit Surfaces & Marching Cubes/Tetras Collision Detection & Conservative Bounding Regions Spatial Acceleration Data Structures Octree, k-d tree, BSF tree 1 Today Particle

More information

PHYSICALLY BASED ANIMATION

PHYSICALLY BASED ANIMATION PHYSICALLY BASED ANIMATION CS148 Introduction to Computer Graphics and Imaging David Hyde August 2 nd, 2016 WHAT IS PHYSICS? the study of everything? WHAT IS COMPUTATION? the study of everything? OUTLINE

More information

Chapter 5: Introduction to Differential Analysis of Fluid Motion

Chapter 5: Introduction to Differential Analysis of Fluid Motion Chapter 5: Introduction to Differential 5-1 Conservation of Mass 5-2 Stream Function for Two-Dimensional 5-3 Incompressible Flow 5-4 Motion of a Fluid Particle (Kinematics) 5-5 Momentum Equation 5-6 Computational

More information

AN ADAPTIVE SAMPLING APPROACH TO INCOMPRESSIBLE PARTICLE-BASED FLUID. A Dissertation WOO-SUCK HONG

AN ADAPTIVE SAMPLING APPROACH TO INCOMPRESSIBLE PARTICLE-BASED FLUID. A Dissertation WOO-SUCK HONG AN ADAPTIVE SAMPLING APPROACH TO INCOMPRESSIBLE PARTICLE-BASED FLUID A Dissertation by WOO-SUCK HONG Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements

More information

CS 231. Fluid simulation

CS 231. Fluid simulation CS 231 Fluid simulation Why Simulate Fluids? Feature film special effects Computer games Medicine (e.g. blood flow in heart) Because it s fun Fluid Simulation Called Computational Fluid Dynamics (CFD)

More information

Lecture 1.1 Introduction to Fluid Dynamics

Lecture 1.1 Introduction to Fluid Dynamics Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory

More information

Interaction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang

Interaction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang 4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 2015) Interaction of Fluid Simulation Based on PhysX Physics Engine Huibai Wang, Jianfei Wan, Fengquan Zhang College

More information

Realtime Water Simulation on GPU. Nuttapong Chentanez NVIDIA Research

Realtime Water Simulation on GPU. Nuttapong Chentanez NVIDIA Research 1 Realtime Water Simulation on GPU Nuttapong Chentanez NVIDIA Research 2 3 Overview Approaches to realtime water simulation Hybrid shallow water solver + particles Hybrid 3D tall cell water solver + particles

More information

Animation of Fluids. Animating Fluid is Hard

Animation of Fluids. Animating Fluid is Hard Animation of Fluids Animating Fluid is Hard Too complex to animate by hand Surface is changing very quickly Lots of small details In short, a nightmare! Need automatic simulations AdHoc Methods Some simple

More information

FLUID SIMULATION. Kristofer Schlachter

FLUID SIMULATION. Kristofer Schlachter FLUID SIMULATION Kristofer Schlachter The Equations Incompressible Navier-Stokes: @u @t = (r u)u 1 rp + vr2 u + F Incompressibility condition r u =0 Breakdown @u @t The derivative of velocity with respect

More information

Divergence-Free Smoothed Particle Hydrodynamics

Divergence-Free Smoothed Particle Hydrodynamics Copyright of figures and other materials in the paper belongs to original authors. Divergence-Free Smoothed Particle Hydrodynamics Bender et al. SCA 2015 Presented by MyungJin Choi 2016-11-26 1. Introduction

More information

An added mass partitioned algorithm for rigid bodies and incompressible flows

An added mass partitioned algorithm for rigid bodies and incompressible flows An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216 Collaborators Bill Henshaw,

More information

The 3D DSC in Fluid Simulation

The 3D DSC in Fluid Simulation The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations

More information

CHAPTER 3. Elementary Fluid Dynamics

CHAPTER 3. Elementary Fluid Dynamics CHAPTER 3. Elementary Fluid Dynamics - Understanding the physics of fluid in motion - Derivation of the Bernoulli equation from Newton s second law Basic Assumptions of fluid stream, unless a specific

More information

Water. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things:

Water. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things: Notes Added a 2D cross-section viewer for assignment 6 Not great, but an alternative if the full 3d viewer isn t working for you Warning about the formulas in Fedkiw, Stam, and Jensen - maybe not right

More information

MOTION CONTROL FOR FLUID ANIMATION: FLOW ALONG A CONTROL PATH

MOTION CONTROL FOR FLUID ANIMATION: FLOW ALONG A CONTROL PATH MOTION CONTROL FOR FLUID ANIMATION: FLOW ALONG A CONTROL PATH A Thesis Presented in Partial Fulfillment of the Requirements for Graduation with Distinction in the College of Engineering of The Ohio State

More information

Numerical Modeling Study for Fish Screen at River Intake Channel ; PH (505) ; FAX (505) ;

Numerical Modeling Study for Fish Screen at River Intake Channel ; PH (505) ; FAX (505) ; Numerical Modeling Study for Fish Screen at River Intake Channel Jungseok Ho 1, Leslie Hanna 2, Brent Mefford 3, and Julie Coonrod 4 1 Department of Civil Engineering, University of New Mexico, Albuquerque,

More information

A Novel Approach to High Speed Collision

A Novel Approach to High Speed Collision A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.

More information

Fluids in Games. Jim Van Verth Insomniac Games

Fluids in Games. Jim Van Verth Insomniac Games Fluids in Games Jim Van Verth Insomniac Games www.insomniacgames.com jim@essentialmath.com Introductory Bits General summary with some details Not a fluids expert Theory and examples What is a Fluid? Deformable

More information

Overview of Traditional Surface Tracking Methods

Overview of Traditional Surface Tracking Methods Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx

More information

ENERGY-224 Reservoir Simulation Project Report. Ala Alzayer

ENERGY-224 Reservoir Simulation Project Report. Ala Alzayer ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................

More information

Finite Volume Discretization on Irregular Voronoi Grids

Finite Volume Discretization on Irregular Voronoi Grids Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling

More information

3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic

3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic ISCS 2013 Selected Papers Dam-break effect on a Solid Wall 1 3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic Suprijadi a,b, F. Faizal b, C.F. Naa a and A.Trisnawan

More information

Rigid Body Dynamics, Collision Response, & Deformation

Rigid Body Dynamics, Collision Response, & Deformation Rigid Body Dynamics, Collision Response, & Deformation Pop Worksheet! Teams of 2. SOMEONE YOU HAVEN T ALREADY WORKED WITH What are the horizontal and face velocities after 1, 2, and many iterations of

More information

A Particle Cellular Automata Model for Fluid Simulations

A Particle Cellular Automata Model for Fluid Simulations Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 36(2), 2009, Pages 35 41 ISSN: 1223-6934 A Particle Cellular Automata Model for Fluid Simulations Costin-Radu Boldea Abstract. A new cellular-automaton

More information

Numerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1

Numerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1 Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns

More information

Math 690N - Final Report

Math 690N - Final Report Math 690N - Final Report Yuanhong Li May 05, 008 Accurate tracking of a discontinuous, thin and evolving turbulent flame front has been a challenging subject in modelling a premixed turbulent combustion.

More information

FEMLAB Exercise 1 for ChE366

FEMLAB Exercise 1 for ChE366 FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let

More information

Development of the Compliant Mooring Line Model for FLOW-3D

Development of the Compliant Mooring Line Model for FLOW-3D Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship

More information

Particle-Based Fluid Simulation. CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016

Particle-Based Fluid Simulation. CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016 Particle-Based Fluid Simulation CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016 Del Operations Del: = x Gradient: s = s x y s y z s z Divergence: v = v x + v y + v z x y z Curl: v = v z v

More information

Simulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method

Simulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method Simulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method June 21, 2011 Introduction Free Surface LBM Liquid-Gas-Solid Flows Parallel Computing Examples and More References Fig. Simulation

More information

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances

More information

Particle-based Fluid Simulation

Particle-based Fluid Simulation Simulation in Computer Graphics Particle-based Fluid Simulation Matthias Teschner Computer Science Department University of Freiburg Application (with Pixar) 10 million fluid + 4 million rigid particles,

More information

Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods

Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,

More information

Permeable and Absorbent Materials in Fluid Simulations

Permeable and Absorbent Materials in Fluid Simulations Permeable and Absorbent Materials in Fluid Simulations Nate Andrysco Bedrich Benes Department of Computer Science Slide 1 Motivation Fluid simulations are geared toward impermeable materials What if you

More information

Simulating Fluid-Solid Interaction

Simulating Fluid-Solid Interaction Simulating Fluid-Solid Interaction Olivier Génevaux Arash Habibi Jean-Michel Dischler LSIIT UMR CNRS-ULP 7005 Boulevard Sébastien Brant, F67400 Illkirch, France {genevaux, habibi, dischler}@dpt-info.u-strasbg.fr

More information

Solidification using Smoothed Particle Hydrodynamics

Solidification using Smoothed Particle Hydrodynamics Solidification using Smoothed Particle Hydrodynamics Master Thesis CA-3817512 Game and Media Technology Utrecht University, The Netherlands Supervisors: dr. ir. J.Egges dr. N.G. Pronost July, 2014 - 2

More information

Fluid Simulation. Dhruv Kore, Giancarlo Gonzalez, and Jenny Sum CS 488: Introduction to Computer Graphics Professor Angus Forbes

Fluid Simulation. Dhruv Kore, Giancarlo Gonzalez, and Jenny Sum CS 488: Introduction to Computer Graphics Professor Angus Forbes Fluid Simulation Dhruv Kore, Giancarlo Gonzalez, and Jenny Sum CS 488: Introduction to Computer Graphics Professor Angus Forbes 1) Introduction Figure 1.1: Realistic water and ice simulation using Blender,

More information

Smoke Simulation using Smoothed Particle Hydrodynamics (SPH) Shruti Jain MSc Computer Animation and Visual Eects Bournemouth University

Smoke Simulation using Smoothed Particle Hydrodynamics (SPH) Shruti Jain MSc Computer Animation and Visual Eects Bournemouth University Smoke Simulation using Smoothed Particle Hydrodynamics (SPH) Shruti Jain MSc Computer Animation and Visual Eects Bournemouth University 21st November 2014 1 Abstract This report is based on the implementation

More information

MODELLING-Representing space and time in a numerical model

MODELLING-Representing space and time in a numerical model MODELLING-Representing space and time in a numerical model Simulation parameters When setting up any simulation, whether computer-based or to be computed manually, it is important to identify the parameters

More information

Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray

Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray 39th Dayton-Cincinnati Aerospace Sciences Symposium Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray Murat Dinc Prof. Donald D. Gray (advisor), Prof. John M. Kuhlman, Nicholas L. Hillen,

More information

1.2 Numerical Solutions of Flow Problems

1.2 Numerical Solutions of Flow Problems 1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian

More information

IMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY DEFINED WALLS FOR PARTICLE BASED FLUID SIMULATION

IMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY DEFINED WALLS FOR PARTICLE BASED FLUID SIMULATION 6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK IMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY

More information

Simple and Fast Fluids

Simple and Fast Fluids Simple and Fast Fluids Martin Guay, Fabrice Colin, Richard Egli To cite this version: Martin Guay, Fabrice Colin, Richard Egli. Simple and Fast Fluids. GPU Pro, A.K. Peters, Ltd., 2011, GPU Pro, pp.433-444.

More information

The Jello Cube Assignment 1, CSCI 520. Jernej Barbic, USC

The Jello Cube Assignment 1, CSCI 520. Jernej Barbic, USC The Jello Cube Assignment 1, CSCI 520 Jernej Barbic, USC 1 The jello cube Undeformed cube Deformed cube The jello cube is elastic, Can be bent, stretched, squeezed,, Without external forces, it eventually

More information

Moving Interface Problems: Methods & Applications Tutorial Lecture II

Moving Interface Problems: Methods & Applications Tutorial Lecture II Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,

More information

Solving Partial Differential Equations on Overlapping Grids

Solving Partial Differential Equations on Overlapping Grids **FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific

More information

NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE

NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE Jungseok Ho 1, Hong Koo Yeo 2, Julie Coonrod 3, and Won-Sik Ahn 4 1 Research Assistant Professor, Dept. of Civil Engineering,

More information

Vector Visualisation 1. global view

Vector Visualisation 1. global view Vector Field Visualisation : global view Visualisation Lecture 12 Institute for Perception, Action & Behaviour School of Informatics Vector Visualisation 1 Vector Field Visualisation : local & global Vector

More information

Two-Phase flows on massively parallel multi-gpu clusters

Two-Phase flows on massively parallel multi-gpu clusters Two-Phase flows on massively parallel multi-gpu clusters Peter Zaspel Michael Griebel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität Bonn Workshop Programming of Heterogeneous

More information

The Simulation of Liquids

The Simulation of Liquids The Simulation of Liquids Mario Ausseloos mario.ausseloos@hotmail.com In partial fulfillment of the requirements for the degree MSc Computer Animation National Centre for Computer Animation Bournemouth

More information

Hardware Accelerated Real-Time Fluid Surface Simulation

Hardware Accelerated Real-Time Fluid Surface Simulation Hardware Accelerated Real-Time Fluid Surface Simulation Björn Hellstrandh bjornh@cs.umu.se Jesper Byström c99jbm@cs.umu.se July 1, 2008 Master s Thesis in Computing Science, 2*30 ECTS credits Supervisor

More information

CUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata

CUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata CUDA Fluid simulation Lattice Boltzmann Models Cellular Automata Please excuse my layout of slides for the remaining part of the talk! Fluid Simulation Navier Stokes equations for incompressible fluids

More information

Driven Cavity Example

Driven Cavity Example BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square

More information

6 Fluid. Chapter 6. Fluids. Department of Computer Science and Engineering 6-1

6 Fluid. Chapter 6. Fluids. Department of Computer Science and Engineering 6-1 Fluids 6-1 Among the most difficult graphical objects to model and animate are those that are not defined by a static, rigid, topological simple structure. Many of these complex forms are found in nature.

More information

Microwell Mixing with Surface Tension

Microwell Mixing with Surface Tension Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical

More information

Coastal impact of a tsunami Review of numerical models

Coastal impact of a tsunami Review of numerical models Coastal impact of a tsunami Review of numerical models Richard Marcer 2 Content Physics to simulate Different approaches of modelling 2D depth average Full 3D Navier-Stokes 3D model Key point : free surface

More information

Technical Report TR

Technical Report TR Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,

More information

Chapter 6. Semi-Lagrangian Methods

Chapter 6. Semi-Lagrangian Methods Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called

More information

Literature Report. Daniël Pols. 23 May 2018

Literature Report. Daniël Pols. 23 May 2018 Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +

More information

2.29 Numerical Marine Hydrodynamics Spring 2007

2.29 Numerical Marine Hydrodynamics Spring 2007 Numerical Marine Hydrodynamics Spring 2007 Course Staff: Instructor: Prof. Henrik Schmidt OCW Web Site: http://ocw.mit.edu/ocwweb/mechanical- Engineering/2-29Spring-2003/CourseHome/index.htm Units: (3-0-9)

More information

Chapter 2 Water and Bubbles

Chapter 2 Water and Bubbles Chapter 2 Water and Bubbles 2.1 Animation of Bubbles in Liquid Abstract This section introduces a new fluid animation technique in which liquid and gas interact with each other, using the example of bubbles

More information

A higher-order finite volume method with collocated grid arrangement for incompressible flows

A higher-order finite volume method with collocated grid arrangement for incompressible flows Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.

More information

A High Quality, Eulerian 3D Fluid Solver in C++ A Senior Project. presented to. the Faculty of the Computer Science Department of

A High Quality, Eulerian 3D Fluid Solver in C++ A Senior Project. presented to. the Faculty of the Computer Science Department of A High Quality, Eulerian 3D Fluid Solver in C++ A Senior Project presented to the Faculty of the Computer Science Department of California Polytechnic State University, San Luis Obispo In Partial Fulfillment

More information

Simulation of Offshore Wave Impacts with a Volume of Fluid Method. Tim Bunnik Tim Bunnik MARIN

Simulation of Offshore Wave Impacts with a Volume of Fluid Method. Tim Bunnik Tim Bunnik MARIN Simulation of Offshore Wave Impacts with a Volume of Fluid Method Tim Bunnik Tim Bunnik MARIN Outline Part I: Numerical method -Overview Part II: Applications - Dambreak - Wave run-up - Sloshing loads

More information

More Animation Techniques

More Animation Techniques CS 231 More Animation Techniques So much more Animation Procedural animation Particle systems Free-form deformation Natural Phenomena 1 Procedural Animation Rule based animation that changes/evolves over

More information

NUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.

NUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document

More information

SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH

SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH Fifth International Conference on CFD in the Process Industries CSIRO, Melbourne, Australia 3-5 December 006 SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH Manasa Ranjan Behera and Murali.

More information

Modeling and simulation the incompressible flow through pipelines 3D solution for the Navier-Stokes equations

Modeling and simulation the incompressible flow through pipelines 3D solution for the Navier-Stokes equations Modeling and simulation the incompressible flow through pipelines 3D solution for the Navier-Stokes equations Daniela Tudorica 1 (1) Petroleum Gas University of Ploiesti, Department of Information Technology,

More information

Solution for Euler Equations Lagrangian and Eulerian Descriptions

Solution for Euler Equations Lagrangian and Eulerian Descriptions Solution for Euler Equations Lagrangian and Eulerian Descriptions Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com Abstract We find an exact solution for the system of Euler equations, following

More information

Auto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial

Auto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial Auto Injector Syringe A Fluent Dynamic Mesh 1DOF Tutorial 1 2015 ANSYS, Inc. June 26, 2015 Prerequisites This tutorial is written with the assumption that You have attended the Introduction to ANSYS Fluent

More information

Physically Based Simulation and Animation of Gaseous Phenomena in a Periodic Domain

Physically Based Simulation and Animation of Gaseous Phenomena in a Periodic Domain Physically Based Simulation and Animation of Gaseous Phenomena in a Periodic Domain Andrew Nealen Department of Computer Science, University of British Columbia CPSC 533b: Algorithmic Animation, Project

More information

NUMERICAL SIMULATION OF SHALLOW WATERS EFFECTS ON SAILING SHIP "MIRCEA" HULL

NUMERICAL SIMULATION OF SHALLOW WATERS EFFECTS ON SAILING SHIP MIRCEA HULL NUMERICAL SIMULATION OF SHALLOW WATERS EFFECTS ON SAILING SHIP "MIRCEA" HULL Petru Sergiu ȘERBAN 1 1 PhD Student, Department of Navigation and Naval Transport, Mircea cel Batran Naval Academy, Constanța,

More information

Volcanic Smoke Animation using CML

Volcanic Smoke Animation using CML Volcanic Smoke Animation using CML Ryoichi Mizuno Yoshinori Dobashi Tomoyuki Nishita The University of Tokyo Tokyo, Japan Hokkaido University Sapporo, Hokkaido, Japan {mizuno,nis}@nis-lab.is.s.u-tokyo.ac.jp

More information

Computer animation for fluid simulation of a high viscous fluid melting

Computer animation for fluid simulation of a high viscous fluid melting 2 nd WIETE Annual Conference on Engineering and Technology Education 2011 WIETE Pattaya, Thailand, 25-28 January 2011 Computer animation for fluid simulation of a high viscous fluid melting S. Lertkulvanich

More information

A Contact Angle Model for the Parallel Free Surface Lattice Boltzmann Method in walberla Stefan Donath (stefan.donath@informatik.uni-erlangen.de) Computer Science 10 (System Simulation) University of Erlangen-Nuremberg

More information

A 3D VOF model in cylindrical coordinates

A 3D VOF model in cylindrical coordinates A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved

More information

Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić

Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume

More information

Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow

Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training

More information

Animation. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 4/23/07 1

Animation. Computer Graphics COMP 770 (236) Spring Instructor: Brandon Lloyd 4/23/07 1 Animation Computer Graphics COMP 770 (236) Spring 2007 Instructor: Brandon Lloyd 4/23/07 1 Today s Topics Interpolation Forward and inverse kinematics Rigid body simulation Fluids Particle systems Behavioral

More information

WAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD

WAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD Proceedings of the ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering OMAE2011 June 19-24, 2011, Rotterdam, The Netherlands OMAE2011-49593 WAVE PATTERNS, WAVE INDUCED FORCES

More information

Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi

Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow

More information

Animation of Bubbles in Liquid

Animation of Bubbles in Liquid EUROGRAPHICS 2003 / P. Brunet and D. Fellner (Guest Editors) Volume 22 (2003), Number 3 Animation of Bubbles in Liquid Jeong-Mo Hong and Chang-Hun Kim Department of Computer Science, Korea University Abstract

More information

Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)

Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,

More information

1 Mathematical Concepts

1 Mathematical Concepts 1 Mathematical Concepts Mathematics is the language of geophysical fluid dynamics. Thus, in order to interpret and communicate the motions of the atmosphere and oceans. While a thorough discussion of the

More information

Georgia Institute of Technology. Atlanta, GA [obrienj j dynamic models of the uid and the impacting. object.

Georgia Institute of Technology. Atlanta, GA [obrienj j dynamic models of the uid and the impacting. object. Dynamic Simulation of Splashing Fluids James F. O'Brien and Jessica K. Hodgins College of Computing Georgia Institute of Technology Atlanta, GA 30332-0280 [obrienj j jkh]@cc.gatech.edu To appear in Proceedings

More information

A MULTIPHASE CFD METHOD FOR PREDICTION OF FLOODWATER DYNAMICS

A MULTIPHASE CFD METHOD FOR PREDICTION OF FLOODWATER DYNAMICS 10 th International Conference 307 A MULTIPHASE CFD METHOD FOR PREDICTION OF FLOODWATER DYNAMICS Zhiliang Gao 1, Dracos Vassalos 1, Qiuxin Gao 1 1 Department of Naval Architecture and Marine Engineering,

More information

Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field

Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi

More information