Smoothed Particle Hydrodynamics (SPH) 4. May 2012
|
|
- Lynne Ball
- 6 years ago
- Views:
Transcription
1 Smoothed Particle Hydrodynamics (SPH) 4. May 2012
2 Calculating density
3 SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension of inverse volume. h scale parameter determining fall-off rate of W. Conservation of mass implies W normalisation: V W (r r b, h)dv = 1
4 The smoothing kernel (Weight function) Positive weighting, monotonically decreasing with relative distance, smooth derivatives. Symmetry with respect to (r r ) i.e., W (r r, h) W ( r r, h) Flat central portion (No strong influence by small change in near neighbour position change) A natural choice that satisfies all of the properties is the Gaussian: d number of dimensions. σ normalisation factor. ] W (r r, h) = σ exp [ (r r ) 2 h d h 2
5 The smoothing kernel (Weight function) Problem with Gaussian: Interaction with all particles in the domain necessary, although the contribution of particles quickly becomes negligible with increasing distance. Better choice: Gaussian shaped kernel, but truncated at some finite distance (a few times the scale length h). A popular choice is the Schoenberg (1946) B-Spline: w(q) = σ 1 4 (2 q)3 (1 q) 3, 0 q < 1; 1 4 (2 q)3, 1 q < 2; 0. q 2, There is a large number of alternative Kernels.
6 The smoothing length (scale length) h Specifies the fall-off of the kernel weighting with respect the to particle separation. The number of neighbours (inside the smoothing kernel) can change without changing h (Higher order B-splines) To resolve clustered and sparse regions evenly a natural choice is to relate h to the local number density of particles: h(r) n(r) 1/d n(r) = b W [r r b, h(r)]. For equal mass particles this is equivalent to making h proportional to density itself. Since in turn density itself is a function of h, this leads to the idea of iterative summation to simultaneous obtain density and h. This is done for the two equations: ρ(r a ) = b m bw (r a r b, h a ) h(r a ) = η ( ma ρ a ) 1/d
7 From density to equations of motion How to calculate arbitrary quantities A(r) (which could be velocity, energy, etc.) using the density estimate? First use the identity A(r) = A(r )δ(r r )dr and approximate using the smoothing kernel instead of the Dirac delta function. A(r) = A(r )W (r r, h)dr + O(h 2 ) Now discretise onto a finite set of interpolation points. A(r) = A(r ) ρ(r ) W (r r, h)ρ(r )dr N neigh b=1 m b A b ρ b W (r r b, h)
8 From density to equations of motion The summation interpolant forms the base of all SPH formalisms. Gradient terms can now be straightforwardly calculated. A(r) b m b A b ρ b W (r r b, h) A(r) b m b A b ρ b W (r r b, h) A(r) b m b A b ρ b W (r r b, h) j A i (r) b m b Ai b ρ b j W (r r b, h) In case of a non-constant smoothing length also the derivative of the smoothing length has to be taken into account.
9 From density to equations of motion (Example: continuity equation) Start with the base expression for ρ: Time derivative: ρ a = b m bw (r a r b, h a ) dρ a dt = b m b (v a v b ) a W ab (h). Rewrite: dρ a dt = v a b m b ρ b ρ b a W ab m b b ρ b (ρ b v b ) a W ab Which is an approximate expression for (derivatives from the last slide): v a ρ (ρv) ρ a ( v) a
10 Non constant smoothing length In the case of a non-constant smoothing length, also the derivative of the smoothing length needs to be taken into account. Using the time derivative of the density again: dρ a dt = b m b (v a v b ) a W ab (h). Changes to: With the correction term Ω a : dρ a dt = 1 Ω a b m b (v a v b ) a W ab (h a ) Ω a [ ] 1 ha ρ a b m b W ab(h a) h a
11 Generalised first derivate operators The gradient operators can be generalised using the expressions A = 1 φ [ (φa) A φ] b m b ρ b φ b φ a (A b A a ) a W ab and [ ( )] A = φ A φ + A φ 2 φ ( ) m b φb b ρ b φ a A a + φa φ b A b a W ab where φ is any arbitrary, differentiable scalar quantity defined on the particles. Various alternative SPH formulations are based on different choices for φ. Two examples using the expression for acceleration: (Hernquist, Katz, 1989) dv dt = b m b (Morris, 1996) dv dt = b m b ( 2 PaP b ρ aρ b ) a W ab ( Pa P b ρ aρ b ) a W ab
12 Choice of the derivative The paradox of SPH is that choosing a less accurate derivate estimator can actually be a good decision. Using the accurate derivate estimator of Morris for a constant density box of particles (randomly distributed, which is a bad choice) the particle-particle force is: m a m b [ Pa P b ρ aρ b ] F ab (ˆr a ˆr b ). Here, the force vanishes when the pressure is constant regardless of the particle distribution. Using a less accurate estimate: [ ] m a m P b F ρ 2 ab (ˆr a ˆr b ). a + P ρ 2 b there will be a net (repulsive) force between the particles until a certain arrangement of particles is reached (locally isotropic and regular).
13 Alternative derivation of the SPH formalism Price derives the equations of motion for SPH from first principles using only ( the) Lagrangian L = T V inserted into the Euler-Lagrange equations d L dt v a L r a = 0 From this the whole SPH equations of motion can also be derived. The full equations of motion are (including the treatment of magnetic fields): dv i a dt = b m b [ S ij a Ω aρ 2 a ( ) with S ij P + 1 2µ 0 B 2 δ ij + 1 µ 0 B i B j ] j aw ab (h a ) + Sij b j Ω b ρ aw 2 ab (h b ) b Induction equation for the evolution of the magnetic field: [ db i dt = 1 N [ Ωρ i j=1 m Bi j ( v ij i W i ) v ij ( B i i W i )]]
14 Higher order derivatives With some work one can introduce higher order derivates to calculate quantities such as 2 A or ( A) Those higher order derivatives can be used to discretise additional physics (depending on higher order derivates) like viscosity or thermal conductivity. Taking viscosity as an example, the standard formulation of artificial viscosity is given by: ( dv ) dt diss = b m b α c s,abµ ab +βµ 2 ab ρ ab ˆr ab F ab µ ab = hv ab r ab r ab +ɛ h ab 2 Barred quantities correspond to an average c s is the sound speed α and β are dimensionless parameters (α = 1, β = 2 normally) ɛ 0.01 is a small parameter to prevent divergences.
15 Example for particle rearrangement and (artificial) viscosity
16 Treatment of discontinuities using artificial viscosity and artificial conductivity du dt diss = mb b ρ ab P h 1 2 αvsig (vab i u (u u ) F r ab )2 + αu vsig a b ab
17 SPH Advantage No geometry constraints Auto resolution Perfect advection Energy conservation Easy NBody coupling Disadvantage Shocks broadened Mixing Conduction by hand Derivative estimator Visualisation
18 SPH Codes Gadget Phantom NDSPMHD SPHNG Seren VINE GCSPH SPARTACUS
Fluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationarxiv: v3 [astro-ph] 2 Jun 2008
Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH Daniel J. Price arxiv:0709.2772v3 [astro-ph] 2 Jun 2008 Abstract School of Physics, University of Exeter, Exeter EX4 4QL, UK In this
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationCSCI1950V Project 4 : Smoothed Particle Hydrodynamics
CSCI1950V Project 4 : Smoothed Particle Hydrodynamics Due Date : Midnight, Friday March 23 1 Background For this project you will implement a uid simulation using Smoothed Particle Hydrodynamics (SPH).
More informationSmooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520
Smooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520 OUTLINE Ø Introduction and Historical Perspective: Ø General Principles: Ø Governing
More informationNumerical Simulations. Duncan Christie
Numerical Simulations Duncan Christie Motivation There isn t enough time to derive the necessary methods to do numerical simulations, but there is enough time to survey what methods and codes are available
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationDynamics of Galaxies: Basics. Frontiers in Numerical Gravitational Astrophysics July 3, 2008
Dynamics of Galaxies: Basics Frontiers in Numerical Gravitational Astrophysics July 3, 2008 c = G = 1 h = 0 Outline (a) Galaxies as physical systems. (b) The Collisionless Boltzmann Equation. (c) N-body
More informationThe Evolution of SPH. J. J. Monaghan Monash University Australia
The Evolution of SPH J. J. Monaghan Monash University Australia planetary disks magnetic fields cosmology radiation star formation lava phase change dam break multiphase waves bodies in water fracture
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationImprovement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method *
Materials Transactions, Vol. 58, No. 3 (2017) pp. 479 to 484 2017 Japan Foundry Engineering Society Improvement of Calculation Stability for Slow Fluid Flow Analysis Using Particle Method * Naoya Hirata
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationSmoothed Particle Hydrodynamics (SPH) Huamin Wang
Smoothed Particle Hydrodynamics (SPH) Huamin Wang Fluid Representation Fluid is represented using a set of particles. Particle i has position x i, velocity v i, and mass m i. Animation Example Using 10
More informationSmoothed Particle Hydrodynamics
I went on to test the program in every way I could devise. I strained it to expose its weaknesses. I ran it for high-mass stars and low-mass stars, for stars orn exceedingly hot and those orn relatively
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
More informationToy stars in one dimension
Mon. Not. R. Astron. Soc. 350, 1449 1456 (004) doi:10.1111/j.1365-966.004.07748.x Toy stars in one dimension J. J. Monaghan 1 and D. J. Price 1 School of Mathematical Sciences, Monash University, Clayton
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationProblem Set #3: 2.11, 2.15, 2.21, 2.26, 2.40, 2.42, 2.43, 2.46 (Due Thursday Feb. 27th)
Chapter Electrostatics Problem Set #3:.,.5,.,.6,.40,.4,.43,.46 (Due Thursday Feb. 7th). Coulomb s Law Coulomb showed experimentally that for two point charges the force is - proportional to each of the
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationarxiv:astro-ph/ v2 22 Nov 1996
Force Error Optimization in Smooth Particle Hydrodynamics 1 R. Capuzzo Dolcetta, 1,2 R. Di Lisio arxiv:astro-ph/9611177v2 22 Nov 1996 1 Istituto Astronomico Università di Roma La Sapienza Via G.M. Lancisi
More informationAbel Inversion using the Maximum Entropy Method
Abel Inversion using the Maximum Entropy Method Danilo R. Neuber Wolfgang von der Linden 3rd October 2003 Inst. für Theoretische Physik, Tel.: +43/3 16/8 73-81 88; neuber@itp.tu-graz.ac.at Inst. für Theoretische
More informationImproving Boundary Conditions in Time- Domain Self-Force Calculations
Improving Boundary Conditions in Time- Domain Self-Force Calculations Carlos F. Sopuerta Institute of Space Sciences National Spanish Research Council Work in Collaboration with Anil Zenginoğlu (Caltech)
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationNumber-Flux Vector and Stress-Energy Tensor
Massachusetts Institute of Technology Department of Physics Physics 8.962 Spring 2002 Number-Flux Vector and Stress-Energy Tensor c 2000, 2002 Edmund Bertschinger. All rights reserved. 1 Introduction These
More informationCosmic Structure Formation on Supercomputers (and laptops)
Cosmic Structure Formation on Supercomputers (and laptops) Lecture 4: Smoothed particle hydrodynamics and baryonic sub-grid models Benjamin Moster! Ewald Puchwein 1 Outline of the lecture course Lecture
More informationarxiv: v1 [astro-ph.co] 9 Dec 2009
Astronomy & Astrophysics manuscript no. EvoL I arxiv c ESO 29 December 9, 29 EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics Emiliano Merlin,
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
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 informationDaniel Price Monash University Melbourne, Australia. Guillaume Laibe University of St Andrews, Scotland
Dust+gas in SPH en.wikipedia.org Daniel Price Monash University Melourne, Australia Guillaume Laie University of St Andrews, Scotland Image: Gemini Oservatory/ AURA Artwork y Lynette Cook HST/ 10th SPHERIC
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationLecture IX: Field equations, cosmological constant, and tides
Lecture IX: Field equations, cosmological constant, and tides Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: October 28, 2011) I. OVERVIEW We are now ready to construct Einstein
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationIntroduction to Electrostatics
Chapter 1 Introduction to Electrostatics Problem Set #1: 1.5, 1.7, 1.12 (Due Monday Feb. 11th) 1.1 Electric field Coulomb showed experimentally that for two point charges the force is -proportionaltoeachofthecharges,
More informationSink particle accretion test
Sink particle accretion test David A. Hubber & Stefanie Walch 1 Objectives Simulate spherically-symmetric Bondi accretion onto a sink particle for an isothermal gas. Calculate the accretion rate onto a
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationChapter 3. Forces, Momentum & Stress. 3.1 Newtonian mechanics: a very brief résumé
Chapter 3 Forces, Momentum & Stress 3.1 Newtonian mechanics: a very brief résumé In classical Newtonian particle mechanics, particles (lumps of matter) only experience acceleration when acted on by external
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationIntroduction and Review Lecture 1
Introduction and Review Lecture 1 1 Fields 1.1 Introduction This class deals with classical electrodynamics. Classical electrodynamics is the exposition of electromagnetic interactions between the develoment
More informationClassical and Quantum Mechanics of a Charged Particle Moving in Electric and Magnetic Fields
Classical Mechanics Classical and Quantum Mechanics of a Charged Particle Moving in Electric and Magnetic Fields In this section I describe the Lagrangian and the Hamiltonian formulations of classical
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationarxiv: v1 [astro-ph.im] 23 Feb 2011
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 25 February 2011 (MN LATEX style file v2.2) Divergence-free Interpolation of Vector Fields From Point Values - Exact B = 0 in Numerical Simulations
More informationTutorial on the smooth particle-mesh Ewald algorithm
Tutorial on the smooth particle-mesh Ewald algorithm These notes explain the smooth particle-mesh Ewald algorithm (SPME) in detail. The accompanying library libpme6 is a fully-featured implementation of
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
More informationGeneral Relativity and Cosmology Mock exam
Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers
More informationSPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology
SPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology G. Bilotta 1,2 1 Dipartimento di Matematica e Informatica, Università di Catania, Italy 2 Istituto Nazionale Geofisica e Vulcanologia,
More informationTruncation Error Estimates for Mesh Refinement in Lagrangian Hydrocodes Sarah Lianne Cole
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Truncation Error Estimates for Mesh Refinement in Lagrangian Hydrocodes Sarah Lianne Cole Thesis submitted for the degree of Doctor of Philosophy
More informationTopics in Relativistic Astrophysics
Topics in Relativistic Astrophysics John Friedman ICTP/SAIFR Advanced School in General Relativity Parker Center for Gravitation, Cosmology, and Astrophysics Part I: General relativistic perfect fluids
More informationModeling of Water Flows around a Circular Cylinder with the SPH Method
Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 1 2, pp. 39 60 DOI: 10.1515/heem-2015-0003 IBW PAN, ISSN 1231 3726 Modeling of Water Flows around a Circular Cylinder with
More informationApplications and Analysis
SPH and α-sph: Applications and Analysis John Anthony Mansour B.Sc.(Hons) Thesis sumitted for the degree of Doctor of Philosophy School of Mathematical Sciences Monash University Octoer 2007 for my parents
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationOn Fluid Maxwell Equations
On Fluid Maxwell Equations Tsutomu Kambe, (Former Professor, Physics), University of Tokyo, Abstract Fluid mechanics is a field theory of Newtonian mechanics of Galilean symmetry, concerned with fluid
More informationSPH development at Cranfield University. Prof. Rade Vignjevic. IV SPHERIC Nantes, May 2009
SPH development at Cranfield University Prof. Rade Vignevic IV SPHERIC Nantes, May 2009 SPH development at Cranfield University Presentation outline 1. Introduction (CU, Motivation) 2. Normalised Corrected
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More informationAstronomy, Astrophysics, and Cosmology
Astronomy, Astrophysics, and Cosmology Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson VI March 15, 2016 arxiv:0706.1988 L. A. Anchordoqui (CUNY)
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationDivergence-free interpolation of vector fields from point values exact B = 0 in numerical simulations
Mon. Not. R. Astron. Soc. 413, L76 L80 (2011) doi:10.1111/j.1745-3933.2011.01037.x Divergence-free interpolation of vector fields from point values exact B = 0 in numerical simulations Colin P. McNally
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationPhysics 505 Fall Homework Assignment #9 Solutions
Physics 55 Fall 25 Textbook problems: Ch. 5: 5.2, 5.22, 5.26 Ch. 6: 6.1 Homework Assignment #9 olutions 5.2 a) tarting from the force equation (5.12) and the fact that a magnetization M inside a volume
More information5.1 Fluid momentum equation Hydrostatics Archimedes theorem The vorticity equation... 42
Chapter 5 Euler s equation Contents 5.1 Fluid momentum equation........................ 39 5. Hydrostatics................................ 40 5.3 Archimedes theorem........................... 41 5.4 The
More informationWe start with some important background material in classical and quantum mechanics.
Chapter Basics We start with some important background material in classical and quantum mechanics.. Classical mechanics Lagrangian mechanics Compared to Newtonian mechanics, Lagrangian mechanics has the
More information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationNotes. Multi-Dimensional Plasticity. Yielding. Multi-Dimensional Yield criteria. (so rest state includes plastic strain): #=#(!
Notes Multi-Dimensional Plasticity! I am back, but still catching up! Assignment is due today (or next time I!m in the dept following today)! Final project proposals: I haven!t sorted through my email,
More informationMaxwell s equations for electrostatics
Maxwell s equations for electrostatics October 6, 5 The differential form of Gauss s law Starting from the integral form of Gauss s law, we treat the charge as a continuous distribution, ρ x. Then, letting
More informationIntroduction to Fluid Dynamics
Introduction to Fluid Dynamics Roger K. Smith Skript - auf englisch! Umsonst im Internet http://www.meteo.physik.uni-muenchen.de Wählen: Lehre Manuskripte Download User Name: meteo Password: download Aim
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationComputational challenges in Numerical Weather Prediction
Computational challenges in Numerical Weather Prediction Mike Cullen Oxford 15 September 2008 Contents This presentation covers the following areas Historical background Current challenges Why does it
More informationConstrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds Terrence S. Tricco a,b,, Daniel J. Price c, Matthew R. Bate b a Canadian Institute for
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationElastic Wave Theory. LeRoy Dorman Room 436 Ritter Hall Tel: Based on notes by C. R. Bentley. Version 1.
Elastic Wave Theory LeRoy Dorman Room 436 Ritter Hall Tel: 4-2406 email: ldorman@ucsd.edu Based on notes by C. R. Bentley. Version 1.1, 20090323 1 Chapter 1 Tensors Elastic wave theory is concerned with
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 (Many figures from C. M. Bishop, "Pattern Recognition and ") 1of 254 Part V
More informationSlide05 Haykin Chapter 5: Radial-Basis Function Networks
Slide5 Haykin Chapter 5: Radial-Basis Function Networks CPSC 636-6 Instructor: Yoonsuck Choe Spring Learning in MLP Supervised learning in multilayer perceptrons: Recursive technique of stochastic approximation,
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationMeasuring the effects of artificial viscosity in SPH simulations of rotating fluid flows
Mon. Not. R. Astron. Soc. 426, 1687 1700 (2012) doi:10.1111/j.1365-2966.2012.21747.x Measuring the effects of artificial viscosity in SPH simulations of rotating fluid flows P. A. Taylor 1,2,3 and J. C.
More informationConvergence of AMR and SPH simulations I. Hydrodynamical resolution and convergence tests
MNRAS 432, 711 727 (2013) Advance Access publication 2013 April 22 doi:10.1093/mnras/stt509 Convergence of AMR and SPH simulations I. Hydrodynamical resolution and convergence tests D. A. Hubber, 1,2,3
More informationConstrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds Terrence S. Tricco a,b,, Daniel J. Price c, Matthew R. Bate b a Canadian Institute for
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationWAVE-STRUCTURE INTERACTION USING SMOOTHED PARTICLE HYDRODYNAMICS
Fifth International Conference on CFD in the Process Industries CSIRO, Melbourne, Australia 13-1 December 26 WAVE-STRUCTURE INTERACTION USING SMOOTHED PARTICLE HYDRODYNAMICS Todd B. SILVESTER 1 and Paul
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationPhysics-Based Animation
CSCI 5980/8980: Special Topics in Computer Science Physics-Based Animation 13 Fluid simulation with grids October 20, 2015 Today Presentation schedule Fluid simulation with grids Course feedback survey
More informationCreated by T. Madas VECTOR OPERATORS. Created by T. Madas
VECTOR OPERATORS GRADIENT gradϕ ϕ Question 1 A surface S is given by the Cartesian equation x 2 2 + y = 25. a) Draw a sketch of S, and describe it geometrically. b) Determine an equation of the tangent
More informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More information