Particle Systems. University of Texas at Austin CS384G - Computer Graphics
|
|
- Muriel Todd
- 5 years ago
- Views:
Transcription
1 Particle Systems University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
2 Reading Required: Witkin, Particle System Dynamics, SIGGRAPH 97 course notes on Physically Based Modeling. Witkin and Baraff, Differential Equation Basics, SIGGRAPH 01 course notes on Physically Based Modeling. Optional Hocknew and Eastwood. Computer simulation using particles. Adam Hilger, New York, Gavin Miller. The motion dynamics of snakes and worms. Computer Graphics 22: , University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 2
3 What are particle systems? A particle system is a collection of point masses that obeys some physical laws (e.g, gravity, heat convection, spring behaviors, ). Particle systems can be used to simulate all sorts of physical phenomena: University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 3
4 Particle in a flow field We begin with a single particle with: Position, " x = x % $ # y& y Velocity, v = x = d x dt " dx dt% = $ # dy dt& g(x,t) Suppose the velocity is actually dictated by some driving function g: x = g( x,t) x University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 4
5 Vector fields At any moment in time, the function g defines a vector field over x: How does our particle move through the vector field? University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 5
6 Diff eqs and integral curves The equation x = g( x,t) is actually a first order differential equation. We can solve for x through time by starting at an initial point and stepping along the vector field: Start Here This is called an initial value problem and the solution is called an integral curve. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 6
7 Euler s method One simple approach is to choose a time step, Δt, and take linear steps along the flow: x (t + Δt) = Writing as a time iteration: x (t) + Δt x(t) = x (t) + Δt g( x,t) x i+1 = x i + Δt v i This approach is called Euler s method and looks like: Properties: Simplest numerical method Bigger steps, bigger errors. Error ~ O(Δt 2 ). Need to take pretty small steps, so not very efficient. Better (more complicated) methods exist, e.g., Runge-Kutta and implicit integration. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 7
8 Particle in a force field Now consider a particle in a force field f. In this case, the particle has: Mass, m Acceleration, a x = d v dt = d 2 x dt 2 The particle obeys Newton s law: f = m a = mx The force field f can in general depend on the position and velocity of the particle as well as time. Thus, with some rearrangement, we end up with: f ( x,x,t) x = m University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 8
9 Second order equations This equation: x = f ( x,x,t) m is a second order differential equation. Our solution method, though, worked on first order differential equations. We can rewrite this as: x = v " f ( x, v % $,t) v = # $ m & where we have added a new variable v to get a pair of coupled first order equations. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 9
10 Phase space " x % $ # v & " x% $ # v& " x% " v % $ = $ # v& # f m& Concatenate x and v to make a 6- vector: position in phase space. Taking the time derivative: another 6-vector. A vanilla 1 st -order differential equation. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 10
11 Differential equation solver Starting with: Applying Euler s method: And making substitutions: Writing this as an iteration, we have: Again, performs poorly for large Δt. " x% " v % $ = $ # v& # f m& x (t + Δt) = x (t) + Δt x(t) x(t + Δt) = x(t) + Δt x(t) x (t + Δt) = x (t) + Δt v (t) x(t + Δt) = x(t) + Δt f ( x,x,t) m x i+1 = x i + Δt v i v i+1 = v f i i + Δt m University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 11
12 Verlet Integration Also called Størmer s Method Invented by Delambre (1791), Størmer (1907), Cowell and Crommelin (1909), Verlet (1960) and probably others More stable than Euler s method (timereversible as well) University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 12
13 Forces Each particle can experience a force which sends it on its merry way. Where do these forces come from? Some examples: Constant (gravity) Position/time dependent (force fields) Velocity-dependent (drag) Combinations (Damped springs) How do we compute the net force on a particle? University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 13
14 Gravity and viscous drag The force due to gravity is simply: f grav = m G p->f += p->m * F->G Often, we want to slow things down with viscous drag: f drag = k v p->f -= F->k * p->v University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 14
15 Damped spring Recall the equation for the force due to a spring: We can augment this with damping: f = k spring (Δ x r) [ ] f = k spring (Δ x r) + k damp v The resulting force equations for a spring between two particles become: * f 1 = k spring Δ x $ Δ v ( r) Δ x - + k damp & % Δ x Δ x, )/ +, (./ Δ x f 2 = " x 1 % p f 1 = $ 1 # v 1 & Δ x = x r = rest length 1 " x 2 % p 2 = $ # & v 2 x 2 University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 15
16 deriveval Clear forces Loop over particles, zero force accumulators Calculate forces Sum all forces into accumulators Return derivatives Loop over particles, return v and f/m " v 1 %" v $ 2 % " v $ n $ # f 1 m 1 &# f 2 m 2 & # f n m n % & Return derivatives to solver 1 Apply forces to particles 3 Clear force accumulators University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 16 2 " $ $ $ f 1 $ # F 1 " $ $ $ f 1 $ # x 1 v 1 x 1 v 1 m 1 m 1 F 2 F 3 F nf % & % & " $ $ $ $ # " $ $ $ $ # x 2 v 2 f 2 x 2 v 2 f 2 m 2 m 2 % " x n % $ v $ n $ f n $ & # & m n % " x n % $ v n $ $ f n $ & # & m n
17 Bouncing off the walls P N Add-on for a particle simulator For now, just simple point-plane collisions A plane is fully specified by any point P on the plane and its normal N. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 17
18 Collision Detection How do you decide when you ve crossed a plane? v P N x University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 18
19 Normal and tangential velocity To compute the collision response, we need to consider the normal and tangential components of a particle s velocity. v N v P N v vt v N = ( N v ) N v T = v v N University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 19
20 Collision Response v N k restitution v N v v T v T " v before v " = v T k restitution v N after Without backtracking, the response may not be enough to bring a particle to the other side of a wall. In that case, detection should include a velocity check: University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell 20
Particle Systems. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
Paticle Systems Univesity of Texas at Austin CS384G - Compute Gaphics Fall 2010 Don Fussell Reading Requied: Witkin, Paticle System Dynamics, SIGGRAPH 97 couse notes on Physically Based Modeling. Witkin
More informationReading. Required: Witkin, Particle System Dynamics, SIGGRAPH 97 course notes on Physically Based Modeling.
Particle Systems Readig Required: Witki, Particle System Dyamics, SIGGRAPH 97 course otes o Physically Based Modelig. Optioal Witki ad Baraff, Differetial Equatio Basics, SIGGRAPH 97 course otes o Physically
More informationReading. Particle Systems. What are particle systems? Overview
Readig Particle Systems Required: Witki, Particle System Dyamics, SIGGRAPH 97 course otes o Physically Based Modelig. Optioal Witki ad Baraff, Differetial Equatio Basics, SIGGRAPH 97 course otes o Physically
More informationPredicting the future with Newton s Second Law
Predicting the future with Newton s Second Law To represent the motion of an object (ignoring rotations for now), we need three functions x(t), y(t), and z(t), which describe the spatial coordinates of
More informationOptimising the Computational Simulation of Mechanical Models: Designing a Real-Time Mass-Aggregate Physics Engine
Optimising the Computational Simulation of Mechanical Models: Designing a Real-Time Mass-Aggregate Physics Engine Jordan Spoooner July 16, 2015 What is a Physics Engine? Engine: Computer software that
More informationDynamics. Kinematics and Dynamics. Passive vs. Active Dynamics. Passive Dynamics. Passive Dynamics. Particle Systems
Kinematics and Dynamics Dynamics Thomas Funkhouser Princeton University C0S 426, Spring 2004 Kinematics Considers only motion Determined by positions, velocities, accelerations Dynamics Considers underlying
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationPassive Dynamics and Particle Systems COS 426
Passive Dynamics and Particle Systems COS 426 Computer Animation Animation Make objects change over time according to scripted actions Simulation / dynamics Predict how objects change over time according
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationthe gamedesigninitiative at cornell university Lecture 17 Physics and Motion
Lecture 17 The Pedagogical Problem Physics simulation is a very complex topic No way I can address this in a few lectures Could spend an entire course talking about it CS 5643: Physically Based Animation
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationPhysically Based Modeling: Principles and Practice Differential Equation Basics
Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationCS229 Programming Assignment 3: Dynamics
: Assignment Out: Wednesday, Oct. 12, 1999 Assignment Due: Fri Oct 22, 1999 (11:59 pm) 1 Introduction Physically based modeling has been an extremely active area of research in computer graphics for over
More informationTime scales. Notes. Mixed Implicit/Explicit. Newmark Methods [ ] [ ] "!t = O 1 %
Notes Time scales! For using Pixie (the renderer) make sure you type use pixie first! Assignment 1 questions?! [work out]! For position dependence, characteristic time interval is "!t = O 1 % $ ' # K &!
More informationCHAPTER 1 INTRODUCTION TO NUMERICAL METHOD
CHAPTER 1 INTRODUCTION TO NUMERICAL METHOD Presenter: Dr. Zalilah Sharer 2018 School of Chemical and Energy Engineering Universiti Teknologi Malaysia 16 September 2018 Chemical Engineering, Computer &
More informationModeling and Solving Constraints. Erin Catto Blizzard Entertainment
Modeling and Solving Constraints Erin Catto Blizzard Entertainment Basic Idea Constraints are used to simulate joints, contact, and collision. We need to solve the constraints to stack boxes and to keep
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More information1.1 GRAPHS AND LINEAR FUNCTIONS
MATHEMATICS EXTENSION 4 UNIT MATHEMATICS TOPIC 1: GRAPHS 1.1 GRAPHS AND LINEAR FUNCTIONS FUNCTIONS The concept of a function is already familiar to you. Since this concept is fundamental to mathematics,
More informationComputer Animation. Animation Methods Keyframing Interpolation Kinematics Inverse Kinematics. Slides courtesy of Leonard McMillan and Jovan Popovic
Computer Animation Animation Methods Keyframing Interpolation Kinematics Inverse Kinematics Slides courtesy of Leonard McMillan and Jovan Popovic Lecture 13 6.837 Fall 2002 Office hours Administrative
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationOrdinary Differential Equations
Ordinary Differential Equations In this lecture, we will look at different options for coding simple differential equations. Start by considering bicycle riding as an example. Why does a bicycle move forward?
More informationPHYS 211 Lecture 9 - Examples of 3D motion 9-1
PHYS 211 Lecture 9 - Examples of 3D motion 9-1 Lecture 9 - Examples of 3D motion Text: Fowles and Cassiday, Chap. 4 In one dimension, the equations of motion to be solved are functions of only one position
More informationRigid Body Dynamics and Beyond
Rigid Body Dynamics and Beyond 1 Rigid Bodies 3 A rigid body Collection of particles Distance between any two particles is always constant What types of motions preserve these constraints? Translation,
More informationPSE Game Physics. Session (6) Angular momentum, microcollisions, damping. Oliver Meister, Roland Wittmann
PSE Game Physics Session (6) Angular momentum, microcollisions, damping Oliver Meister, Roland Wittmann 23.05.2014 Session (6)Angular momentum, microcollisions, damping, 23.05.2014 1 Outline Angular momentum
More informationCS 378: Computer Game Technology
CS 378: Computer Game Technolog 3D Engines and Scene Graphs Spring 202 Universit of Teas at Austin CS 378 Game Technolog Don Fussell Representation! We can represent a point, p =,), in the plane! as a
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More information(Refer Slide Time 1:25)
Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 2 Lecture - 24 Transient Response of Pressure Transducers
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationEuler s method for solving a differential equation (approximately)
Euler s method for solving a differential equation (approximately) Department of Mathematics, UW - Madison March 4, 2013 A chemical reaction A B A B A A A chemical reactor contains two kinds of molecules,
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationδ Substituting into the differential equation gives: x i+1 x i δ f(t i,x i ) (3)
Solving Differential Equations Numerically Differential equations are ubiquitous in Physics since the laws of nature often take on a simple form when expressed in terms of infinitesimal changes of the
More informationBoiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics. Derivation of Rayleigh and Rayleigh-Plesset Equations:
Boiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics Derivation of Rayleigh and Rayleigh-Plesset Equations: Let us derive the Rayleigh-Plesset Equation as the Rayleigh equation can
More informationIntroduction to First Order Equations Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationRead textbook CHAPTER 1.4, Apps B&D
Lecture 2 Read textbook CHAPTER 1.4, Apps B&D Today: Derive EOMs & Linearization undamental equation of motion for mass-springdamper system (1DO). Linear and nonlinear system. Examples of derivation of
More informationMatthew S. Norton PHY 112 Lab Report CL-4 Damped Oscillator: Feynman-Newton Method for Solving First-Order Differential Equations March 2, 2006
Matthew S. Norton PHY 112 Lab Report CL-4 Damped Oscillator: Feynman-Newton Method for Solving First-Order Differential Equations March 2, 2006 Norton 0 Norton 1 Abstract In this lab, I used the Feynman-Newton
More informationMATH 1231 MATHEMATICS 1B Calculus Section 3A: - First order ODEs.
MATH 1231 MATHEMATICS 1B 2010. For use in Dr Chris Tisdell s lectures. Calculus Section 3A: - First order ODEs. Created and compiled by Chris Tisdell S1: What is an ODE? S2: Motivation S3: Types and orders
More informationParticle Systems. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017
Particle Systems CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017 Particle Systems Particle systems have been used extensively in computer animation and special effects since their
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationCS 480/680: GAME ENGINE PROGRAMMING
CS 480/680: GAME ENGINE PROGRAMMING PHYSICS 2/14/2013 Santiago Ontañón santi@cs.drexel.edu https://www.cs.drexel.edu/~santi/teaching/2013/cs480-680/intro.html Outline Student Presentations Basics on Rigid
More informationPhysics 202 Laboratory 3. Root-Finding 1. Laboratory 3. Physics 202 Laboratory
Physics 202 Laboratory 3 Root-Finding 1 Laboratory 3 Physics 202 Laboratory The fundamental question answered by this week s lab work will be: Given a function F (x), find some/all of the values {x i }
More informationSimple Harmonic Motion
Simple Harmonic Motion (FIZ 101E - Summer 2018) July 29, 2018 Contents 1 Introduction 2 2 The Spring-Mass System 2 3 The Energy in SHM 5 4 The Simple Pendulum 6 5 The Physical Pendulum 8 6 The Damped Oscillations
More informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationThursday Simulation & Unity
Rigid Bodies Simulation Homework Build a particle system based either on F=ma or procedural simulation Examples: Smoke, Fire, Water, Wind, Leaves, Cloth, Magnets, Flocks, Fish, Insects, Crowds, etc. Simulate
More informationCS888 Advanced Topics in Computer Graphics: Physics-Based Animation. Jan. 7, 2014
CS888 Advanced Topics in Computer Graphics: Physics-Based Animation Jan. 7, 2014 Physics-Based Animation The use of physical simulation to generate animations of: Rigid bodies: Perfectly stiff or rigid
More informationAn Introduction to Physically Based Modeling: An Introduction to Continuum Dynamics for Computer Graphics
An Introduction to Physically Based Modeling: An Introduction to Continuum Dynamics for Computer Graphics Michael Kass Pixar Please note: This document is 1997 by Michael Kass. This chapter may be freely
More informationFigure 1: Doing work on a block by pushing it across the floor.
Work Let s imagine I have a block which I m pushing across the floor, shown in Figure 1. If I m moving the block at constant velocity, then I know that I have to apply a force to compensate the effects
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationIntroduction to Process Control
Introduction to Process Control For more visit :- www.mpgirnari.in By: M. P. Girnari (SSEC, Bhavnagar) For more visit:- www.mpgirnari.in 1 Contents: Introduction Process control Dynamics Stability The
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationGaussians Distributions, Simple Harmonic Motion & Uncertainty Analysis Review. Lecture # 5 Physics 2BL Summer 2015
Gaussians Distributions, Simple Harmonic Motion & Uncertainty Analysis Review Lecture # 5 Physics 2BL Summer 2015 Outline Significant figures Gaussian distribution and probabilities Experiment 2 review
More information3. What type of force is the woman applying to cart in the illustration below?
Name: Forces and Motion STUDY GUIDE Directions: Answer the following questions. 1. What is a force? a. A type of energy b. The rate at which an object performs work c. A push or a pull d. An object that
More informationRevision History Date Version Description Author 03/06/ Initial Version Sergio Miguel Martin. Table of Contents
Icarus Physics Engine Manual 1 - Basic Engine Operation Version 1.0 Revision History Date Version Description Author 03/06/2010 1.0 Initial Version Sergio Miguel Martin Table of Contents 1. Introduction
More informationIntegration of Differential Equations
Integration of Differential Equations Gabriel S. Cabrera August 4, 018 Contents 1 Introduction 1 Theory.1 Linear 1 st Order ODEs................................. 3.1.1 Analytical Solution...............................
More informationVisual Interactive Simulation, TDBD24, Spring 2006
Visual Interactive Simulation, TDBD24, Spring 2006 Lecture 6 Outline Cross product as matrix multiplication The collision matrix The Newton Coulomb impulse for rigid bodies N C updates of translational
More informationPhysics 132 3/31/17. March 31, 2017 Physics 132 Prof. E. F. Redish Theme Music: Benny Goodman. Swing, Swing, Swing. Cartoon: Bill Watterson
March 31, 2017 Physics 132 Prof. E. F. Redish Theme Music: Benny Goodman Swing, Swing, Swing Cartoon: Bill Watterson Calvin & Hobbes 1 Outline The makeup exam Recap: the math of the harmonic oscillator
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationCS-184: Computer Graphics
CS-184: Computer Graphics Lecture #25: Rigid Body Simulations Tobias Pfaff 537 Soda (Visual Computing Lab) tpfaff@berkeley.edu Reminder Final project presentations next week! Game Physics Types of Materials
More informationLecture 17 - Gyroscopes
Lecture 17 - Gyroscopes A Puzzle... We have seen throughout class that the center of mass is a very powerful tool for evaluating systems. However, don t let yourself get carried away with how useful it
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationMATH 312 Section 4.3: Homogeneous Linear Equations with Constant Coefficients
MATH 312 Section 4.3: Homogeneous Linear Equations with Constant Coefficients Prof. Jonathan Duncan Walla Walla College Spring Quarter, 2007 Outline 1 Getting Started 2 Second Order Equations Two Real
More informationPhysics 111. Tuesday, October 05, Momentum
ics Tuesday, ober 05, 2004 Ch 6: Ch 9: Last Example Impulse Momentum Announcements Help this week: Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468 Announcements This week s lab will be
More informationUpdated 2013 (Mathematica Version) M1.1. Lab M1: The Simple Pendulum
Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are
More informationHow to work out really complicated motion. Iteration and Problem Solving Strategies. Let s go. Vertical spring-mass.
Iteration and Problem Solving Strategies How to solve anything! How to work out really complicated motion Break it up into little tiny steps. Use an approximate method for each step. Add them all up. Vertical
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationIntroduction to Controls
EE 474 Review Exam 1 Name Answer each of the questions. Show your work. Note were essay-type answers are requested. Answer with complete sentences. Incomplete sentences will count heavily against the grade.
More informationChapter 3 Acceleration
Chapter 3 Acceleration Slide 3-1 Chapter 3: Acceleration Chapter Goal: To extend the description of motion in one dimension to include changes in velocity. This type of motion is called acceleration. Slide
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationPhysics 6010, Fall Relevant Sections in Text: Introduction
Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationDan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999)
Dan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999) This is intended to be a detailed by fairly-low-math explanation of Stam s Stable Fluids, one of the key papers in a recent series of advances
More informationMotion along a straight line. Physics 11a. 4 Basic Quantities in Kinematics. Motion
Physics 11a Motion along a straight line Motion Position and Average velocity and average speed Instantaneous velocity and speed Acceleration Constant acceleration: A special case Free fall acceleration
More information2.2 Separable Equations
2.2 Separable Equations Definition A first-order differential equation that can be written in the form Is said to be separable. Note: the variables of a separable equation can be written as Examples Solve
More informationspring mass equilibrium position +v max
Lecture 20 Oscillations (Chapter 11) Review of Simple Harmonic Motion Parameters Graphical Representation of SHM Review of mass-spring pendulum periods Let s review Simple Harmonic Motion. Recall we used
More informationProcedural Animation. D.A. Forsyth
Procedural Animation D.A. Forsyth Big points Two important types of procedural animation slice and dice data, like texture synthesis build (approximate) physical simulation Extremely powerful issues how
More informationNumerical methods for physics simulations.
Numerical methods for physics simulations. Claude Lacoursière HPC2N/VRlab, Umeå Universiy, Umeå, Sweden, and CMLabs Simulations, Montréal, Canada May 8, 26 Typeset by FoilTEX May 5th 26 Outline Typeset
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationChapter 3 Acceleration
Chapter 3 Acceleration Slide 3-1 PackBack The first answer gives a good physical picture. The video was nice, and worth the second answer. https://www.youtube.com/w atch?v=m57cimnj7fc Slide 3-2 Slide 3-3
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationFundamentals Physics. Chapter 15 Oscillations
Fundamentals Physics Tenth Edition Halliday Chapter 15 Oscillations 15-1 Simple Harmonic Motion (1 of 20) Learning Objectives 15.01 Distinguish simple harmonic motion from other types of periodic motion.
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationFeedback Basics. David M. Auslander Mechanical Engineering University of California at Berkeley. copyright 1998, D.M. Auslander
Feedback Basics David M. Auslander Mechanical Engineering University of California at Berkeley copyright 1998, D.M. Auslander 1 I. Feedback Control Context 2 What is Feedback Control? Measure desired behavior
More informationLecture 4. Magic Numbers
Lecture 4 In last lecture we discussed how one can solve the Schroedinger equation to determine the allowed energies of a particle trapped in a spherically symmetric potential. We will apply these methods
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
More informationDerivation of Kinematic Equations. View this after Motion on an Incline Lab
Derivation of Kinematic Equations View this after Motion on an Incline Lab Constant velocity Average velocity equals the slope of a position vs time graph when an object travels at constant velocity. v
More informationChapter 3 Acceleration
Chapter 3 Acceleration Slide 3-1 Chapter 3: Acceleration Chapter Goal: To extend the description of motion in one dimension to include changes in velocity. This type of motion is called acceleration. Slide
More informationMon Jan Improved acceleration models: linear and quadratic drag forces. Announcements: Warm-up Exercise:
Math 2250-004 Week 4 notes We will not necessarily finish the material from a given day's notes on that day. We may also add or subtract some material as the week progresses, but these notes represent
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationOscillations Simple Harmonic Motion
Oscillations Simple Harmonic Motion Lana Sheridan De Anza College Dec 1, 2017 Overview oscillations simple harmonic motion (SHM) spring systems energy in SHM pendula damped oscillations Oscillations and
More information