An Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE)
|
|
- Byron Fox
- 5 years ago
- Views:
Transcription
1 An Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE) Shingyu Leung Department of Mathematics, Hong Kong University of Science and Technology May, Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
2 Outline Introduction to FTLE ECS vs. LCS LCS through FTLE Lagrangian Formulation to FTLE 3 Part : Eulerian Formulation to FTLE 4 Part : Propagating the FTLE for many time-levels 5 Advantages and Limitations 6 Part 3: FTLE on a Codimension One Manifold 7 Examples 8 Conclusion Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
3 Introduction to FTLE Eulerian Coherent Structure (ECS) To segment the domain into different regions with similar behavior according to an Eulerian quantity such as the strain, kinetic energy, or the vorticity. Lagrangian Coherent Structure (LCS) To partition the space-time domain into different regions according to a Lagrangian quantity advected along with passive tracers. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
4 One such Lagrangian quantity is the FTLE, rate of separation between adjacent particles over a finite time interval with an infinitesimal perturbation in the initial location. x (t;x,t ) = u(x(t;x,t ),t) with the initial condition x(t ) = x and u : R d R R d. Flow map Φ : Ω R R R d, Φ(x;t,T) = x(t;x,t ) The leading order in the perturbation δx(t) = δx(),[dφ(x;t,t)] DΦ(x;t,T)δx() Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 4 / 33
5 Denoting (x;t,t) the Cauchy-Green deformation tensor (x;t,t) = [DΦ(x;t,T)] DΦ(x;t,T), the largest strength deformation max δx(t) = λ max [ (x;t,t)] e() = exp[σ T (x,t ) T ] e() δx() FTLE σ T (x,t ) = T ln λ max [ (x;t,t)] forward FTLE if T > and the backward FTLE if T <. LCS: (roughly speaking) Ridge of the FTLE. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 5 / 33
6 Lagrangian formulation to FTLE Lagrangian approach to FTLE Compute the flow map Φ(x;t,T) by solving the ODE system x (t;x,t ) = u(x(t;x,t ),t) on an initial Cartesian mesh at time t = t. Construct the deformation tensor (x;t,t) by finite differencing the flow map Φ(x i ;t,t) and then find it s largest eigenvalue. Remarks The velocity field might be defined discretely on a mesh. The numerical accuracy when using a high order adaptive time integrator such as ode45 in MATLAB. 3 For computing the FTLE at the next time step t = t = t + t, all rays obtained on the previous time step t = t are all discarded. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 6 / 33
7 Part : Eulerian formulation to FTLE for one time level Goal: Lagrangian computations Eulerian formulation Define Ψ = (Ψ,Ψ,,Ψ d ) : Ω R R d with the initial condition Ψ(x,t ) = x = (x,x,,x d ). Following the particle trajectory with x = x, any particle identity should be preserved in the Lagrangian framework, i.e. DΨ(x, t) Dt = Ψ(x,t) t +(u )Ψ(x,t) = Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 7 / 33
8 These level set functions defined on a uniform Cartesian mesh give the backward flow map from t = t +T to t = t, i.e. Φ(y;t +T, T) = Ψ(y,t +T) Figure: Lagrangian and Eulerian interpretations of the function Ψ. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 8 / 33
9 Computing the backward FTLE σ T (x,t +T) Discretize the computational domain Initialize Ψ (x i,y j,t ) = x i and Ψ (x i,y j,t ) = y j 3 Solve the Liouville equations Ψ m t +(u )Ψ m =, m =, for up to t = t K = t +T. 4 Compute the Cauchy-Green deformation tensor with Ψ = (Ψ,Ψ ) T, (x i,y j ;t K, T) = [DΨ(x i,y j,t K )] DΨ(x i,y j,t K ) 5 Determine the backward FTLE at (x i,y j,t +T) σ T (x i,y j,t K ) = T ln λ max [ (x i,y j ;t K, T)] Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 9 / 33
10 Computing the forward FTLE σ T (x,t ) Initialize the level set functions at t = t +T by Ψ(x,t +T) = x. Solve the corresponding level set equations backward in time. 3 Determining the Jacobian of the resulting flow map and then computing the largest eigenvalue of the deformation tensor (x;t,t). 4 The forward FTLE is formed by σ T (x,t ) = T ln λ max [ (x;t,t)]. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
11 Part : Propagating the FTLE for many time-levels Goal: Propagate the backward FTLE σ T (x,t +T) forward in time from t +T to T f in order to approximate σ T (x,t) or σ (t t ) (x,t) for t > t +T. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
12 Theorem (Theorem 3. in Shadden, Lekien and Marsden 5) The traditional (forward) Lyapunov exponent is constant along trajectories. Also, the (forward) finite-time Lyapunov exponent becomes constant along trajectories for large integration times T. The difference between σ T (x,t ) and σ T (Φ s (x;t ),t +s) for some arbitrary but fixed s, σ T (x,t ) σ T (Φ s (x;t ),t +s) s T max t σ t t (x,t ) = O( T ) The material derivative of the (forward) FTLE is bounded by Dσ T (x,t) Dt = lim σ T (Φ s (x;t),t +s) σ T (x,t) = O( T ) s s If T is large enough, we obtain the approximation Dσ T (x,t) Dt = Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
13 If T is large enough, we can approximate the solution σ T (x,t) for t > t +T by simply solving the following Liouville equation σ T (x,t) t +u σ T (x,t) = with the initial condition at t = t +T. Let S be the time difference from t = t +T. Since σ T (Φ S (x;t +T),t +T +S) σ T (x,t +T) S ( ) S T max σ t (x,t) = O, t T the error in the approximation is linear in S. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
14 Advantages and limitations of the Eulerian formulation Timestep constraint Lagrangian: t, stiffness of the velocity field Eulerian: t = O( x), CFL condition Computational Complexity d: the dimension of the domain N = O( x ): number of mesh points in each spacial dimension M = T/ t = O( t ): number of time steps for t = T K = (T f t )/ t = O( t ): number of time steps until t = T f (number of time levels of FTLE to be computed) Lagrangian: O(N d MK) = O( x d t t ) Eulerian: O(N d [(d )M +K +]) = O( x d t ) Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 4 / 33
15 Velocity Field CFD solver Real life measurements Storage for the velocity field Lagrangian approach: O(dN d+ ), N d+ mesh points in the space-time domain and a d-vector per node. Eulerian approach: O(dN d ), N d mesh points in the space domain and a d-vector per node. Storage for the flow map and the FTLE Lagrangian approach: O(d) Eulerian approach: O(dN d ) Numerical dissipation in the Eulerian formulation An order of O(S/T) error is introduced in the approximation Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 5 / 33
16 Part 3: FTLE on a Codimension One Manifold Goal: A simple approach to compute the FTLE on a codimension one manifold. A codimension one manifold M(t) in R d evolving in time under the velocity field ũ : M(t) t R d. A possible approach: triangulate the surface, evolve all vertices according to the flow, compute the Jacobian matrix based on a local coordinate system on the surface,... Our proposed approach: the evolving manifold represented implicitly using a level set function. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 6 / 33
17 Computing the backward FTLE σ T (x,t +T) for x M(t +T) With φ(x i,y j,t ) at t = t so that M = {φ = }, solve φ t +(u )φ = for up to t = t +T With Ψ (x i,y j,t ) = x i, Ψ (x i,y j,t ) = y j, solve the Liouville equations Ψ m t +(u )Ψ m =, m =, from t = t to t = t +T 3 Extend the map off the manifold by Ψ m τ +sgn(φ)(n )Ψ m =, m =, 4 Compute the deformation tensor = (DΨ) DΨ 5 Determine the backward FTLE σ T = lnλ max ( ) Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 7 / 33
18 Example: Double-Gyre flow The flow is modeled by the following stream-function ψ(x,y,t) = Asin[πf(x,t)]sin(πy), where f(x,t) = a(t)x +b(t)x, a(t) = ǫsin(ωt), b(t) = ǫsin(ωt) with A =., ω = π/ and ǫ =. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 8 / 33
19 (a) (b) (c) (d) (a) The backward FTLE using x = y = /8 at t = 5, T = 5 and S =. (b) The backward FTLE using x = y = /5 at t = 5, T = 5 and S =. (c) The forward FTLE using x = y = /8 at t =, T = 5 and S =. (d) The forward FTLE using x = y = /5 at t =, T = 5 and S =. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 9 / 33
20 (a) (b) (c) (d) (a) The backward FTLE with x = y = /5 at t = 5. using the Lagrangian approach (MATLAB function ode45) with T = 5. (b) The backward FTLE with x = y = /5 at t = 5. using the Lagrangian approach (RK4 with a fixed timestep) with T = 5. (c) The backward FTLE at t = 5 using the proposed Eulerian approach with x = y = /5 and T = 5 and S =. (d) The backward FTLE at t = 5 using the proposed Eulerian approach with x = y = /5 and T = 5. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
21 (a) (b) (c) (d) (a) the forward FTLE with x = y = /5 at t =. using the Lagrangian approach (MATLAB function ode45) with T = 5. (b) The forward FTLE with x = y = /5 at t =. using the Lagrangian approach (RK4 with a fixed timestep) with T = 5. (d) The forward FTLE at t = using the proposed Eulerian approach with x = y = /5, T = 5 and S =. (c) The forward FTLE at t = using the proposed Eulerian approach with x = y = /5 and T = 5. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
22 x Eul. Rate Lag. ( t = ) Rate Lag. ( t = x) Rate / / / / CPU times (in seconds) for the Eulerian approach with T = 5 and S =, and the Lagrangian approach with t = and t = x. Eulerian approach O( x 3 ) Lagrangian approach O( x 3 t ) Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, / 33
23 Example: Rayleigh-Bénard convection The velocity field is given by u = A ξsinkr cosz, k v = A y sinkr cosz, k w = Asinz (r coskr + k ) sinkr, where A =.4, k =., r = ξ +y, ξ = x g(t), and g(t) =.cos(πt). Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
24 (a) (b) (a) The backward FTLE using x = y = /3 at t = 5, T = and S = 38. (b) The forward FTLE using x = y = /3 at t =, T = and S = 38. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 4 / 33
25 Example: FTLE on a Sphere (a) (b) The forward FTLE on a sphere under the rigid body rotation. The exact solution is zero everywhere on the sphere. The underlying mesh size (a) x = 3/3 and (b) x = 3/8. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 5 / 33
26 (a) (b) The error in the forward FTLE on a sphere under the motion in the outward normal direction. The underlying mesh size (a) x = 3/3 and (b) x = 3/8. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 6 / 33
27 (a) (b) The error in the forward FTLE on a sphere under the motion in the inward normal direction. The underlying mesh size (a) x = 3/3 and (b) x = 3/8. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 7 / 33
28 Example: Four point vortices on a sphere The velocity of any marker particle satisfies the motion x = π 4 i= x i x ( x x i ), with the point vortices centered at (/ 3, / 3,/ 3), (/ 3, / 3, / 3), ( /3,/ 3,) and ( / 3, / 3,/ 3). Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 8 / 33
29 (a) (b) (a) The backward FTLE at t = with T =. (b) The forward FTLE at t =. with T =. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 9 / 33
30 Example: Sphere under a vortex flow The sphere evolves under the vortex motion governed by u = sin (πx)sin(πy)sin(πz) v = sin (πy)sin(πx)sin(πz) w = sin (πz)sin(πy)sin(πx). (a) (b) (a) The backward FTLE at t =.5 with T =.5. (b) The forward FTLE at t =. with T =.5. Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
31 Conclusion Eulerian framework for approximating the FTLE based on the level set formulation Compute the related flow map by solving Liouville equations Determine the FTLE on an evolving manifold implicitly represented by a level set function Future work Extract the Lagrangian coherent structure from the FTLE Combine the current approaches with a CFD solver Incorporate our algorithms with flow fields from real life measurements Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
32 References Shingyu Leung () An Eulerian Approach for Computing the Finite Time Lyapunov Exponent. Journal of Computational Physics, Volume 3, Issue 9, May, Pages Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 3 / 33
33 Thank you Shingyu Leung (HKUST) An Eulerian Approach to the FTLE May, 33 / 33
Hot Topics in Visualization
Hot Topic 1: Illustrative visualization 12 Illustrative visualization: computer supported interactive and expressive visualizations through abstractions as in traditional illustrations. Hot Topics in Visualization
More informationThree-Dimensional Oceanic Flows from Eulerian Velocity Data
Second-year Ph.D. student, Applied Math and Scientific Computing Project Advisor: Kayo Ide Department of Atmospheric and Oceanic Science Center for Scientific Computation and Mathematical Modeling Earth
More informationPACS numbers: Valid PACS appear here Keywords: Suggested keywords
Integrated computation of finite-time Lyapunov exponent fields during direct numerical simulation of unsteady flows Justin Finn a) and Sourabh V. Apte b) School of Mechanical, Industrial and Manufacturing
More informationHot Topics in Visualization. Ronald Peikert SciVis Hot Topics 12-1
Hot Topics in Visualization Ronald Peikert SciVis 2007 - Hot Topics 12-1 Hot Topic 1: Illustrative visualization Illustrative visualization: computer supported interactive and expressive visualizations
More informationA Grid Based Particle Method for Evolution of Open Curves and Surfaces
A Grid Based Particle Method for Evolution of Open Curves and Surfaces Shingyu Leung a,, Hongkai Zhao b a Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Hong
More informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationCooperative Control of a Team of Unmanned Vehicles Using Smoothed Particle Hydrodynamics
Cooperative Control of a Team of Unmanned Vehicles Using Smoothed Particle Hydrodynamics Doug Lipinski and Kamran Mohseni A procedure is outlined for trajectory planning based on Lagrangian coherent structures
More informationInvestigation of Lagrangian Coherent Structures. - To Understand and Identify Turbulence JOHAN JAKOBSSON
Investigation of Lagrangian Coherent Structures - To Understand and Identify Turbulence JOHAN JAKOBSSON Department of Chemical and Biological Engineering Division of Chemical Reaction Engineering CHALMERS
More informationVector Visualization. CSC 7443: Scientific Information Visualization
Vector Visualization Vector data A vector is an object with direction and length v = (v x,v y,v z ) A vector field is a field which associates a vector with each point in space The vector data is 3D representation
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationChapter 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 informationlibcfd2lcs Manual A work in progress... Justin R. Finn June 17, 2016
libcfd2lcs Manual A work in progress... Justin R. Finn June 17, 2016 1 Overview libcfd2lcs is a general purpose numerical library that performs the calculations needed to extract Lagrangian Coherent Structures
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationNumerical 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 informationLevel 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 informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationA gradient-augmented level set method with an optimally local, coherent advection scheme
A gradient-augmented level set method with an optimally local, coherent advection scheme The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationSurfaces and Integral Curves
MODULE 1: MATHEMATICAL PRELIMINARIES 16 Lecture 3 Surfaces and Integral Curves In Lecture 3, we recall some geometrical concepts that are essential for understanding the nature of solutions of partial
More informationChapter 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 informationApplications. Human and animal motion Robotics control Hair Plants Molecular motion
Multibody dynamics Applications Human and animal motion Robotics control Hair Plants Molecular motion Generalized coordinates Virtual work and generalized forces Lagrangian dynamics for mass points
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More informationCS-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 informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationStructure and Synthesis of Robot Motion
Structure and Synthesis of Robot Motion Dynamics: Constraints, Continua, etc. Subramanian Ramamoorthy School of Informatics 5 February, 2009 Recap Last time, we discussed two major approaches to describing
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationReview Initial Value Problems Euler s Method Summary
THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More information10. Cartesian Trajectory Planning for Robot Manipulators
V. Kumar 0. Cartesian rajectory Planning for obot Manipulators 0.. Introduction Given a starting end effector position and orientation and a goal position and orientation we want to generate a smooth trajectory
More informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationPHYSICALLY 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 informationA splitting method for numerical simulation of free surface flows of incompressible fluids with surface tension
A splitting method for numerical simulation of free surface flows of incompressible fluids with surface tension Kirill D. Nikitin Maxim A. Olshanskii Kirill M. Terekhov Yuri V. Vassilevski Abstract The
More informationA brief introduction to fluidstructure. O. Souček
A brief introduction to fluidstructure interactions O. Souček Fluid-structure interactions Important class computational models Civil engineering Biomechanics Industry Geophysics From http://www.ihs.uni-stuttgart.de
More informationVorticity dynamics & the small scales of turbulence Particle Tracking Velocimetry 3DPTV & DNS
Vorticity dynamics & the small scales of turbulence Particle Tracking Velocimetry 3DPTV & DNS Luthi et al. JFM (528) pp 87, 2005 Guala et al. JFM (533) pp 339, 2005 Hoyer et al., Exp. in fluids (39) pp
More informationLecture 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 informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationData Visualization. Fall 2017
Data Visualization Fall 2017 Vector Fields Vector field v: D R n D is typically 2D planar surface or 2D surface embedded in 3D n = 2 fields tangent to 2D surface n = 3 volumetric fields When visualizing
More informationLevel Set Method in a Finite Element Setting
Level Set Method in a Finite Element Setting John Shopple University of California, San Diego November 6, 2007 Outline 1 Level Set Method 2 Solute-Solvent Model 3 Reinitialization 4 Conclusion Types of
More informationVerification of Moving Mesh Discretizations
Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions
More informationLecture overview. Visualisatie BMT. Vector algorithms. Vector algorithms. Time animation. Time animation
Visualisatie BMT Lecture overview Vector algorithms Tensor algorithms Modeling algorithms Algorithms - 2 Arjan Kok a.j.f.kok@tue.nl 1 2 Vector algorithms Vector 2 or 3 dimensional representation of direction
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationMass-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 informationMesh-free Semi-Lagrangian Methods for Transport on a Sphere Using Radial Basis Functions
Mesh-free Semi-Lagrangian Methods for Transport on a Sphere Using Radial Basis Functions Varun Shankar a,, Grady B. Wright b a Department of Mathematics, University of Utah, UT, USA b Department of Mathematics,
More informationM.Sc. in Computational Science. Fundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationContextual Geometric Structures
Contextual Geometric Structures Modeling the fundamental components of cultural behavior Bradly Alicea http://www.msu.edu/~aliceabr/ What is the essence of culture? Heredity (beliefs that propagate)? Plasticity?
More informationFLOWING FLUIDS AND PRESSURE VARIATION
Chapter 4 Pressure differences are (often) the forces that move fluids FLOWING FLUIDS AND PRESSURE VARIATION Fluid Mechanics, Spring Term 2011 e.g., pressure is low at the center of a hurricane. For your
More informationRegistration: Rigid vs. Deformable
Lecture 20 Deformable / Non-Rigid Registration ch. 11 of Insight into Images edited by Terry Yoo, et al. Spring 2017 16-725 (CMU RI) : BioE 2630 (Pitt) Dr. John Galeotti The content of these slides by
More informationSemi-Lagrangian Advection. The Basic Idea. x 1.
3.2.6. Semi-Lagrangian Advection We have studied the Eulerian leapfrog scheme and found it to be conditionally stable. The criterion for stability was the CFL condition µ c t x 1. For high spatial resolution
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationSimulation in Computer Graphics. Introduction. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Introduction Matthias Teschner Computer Science Department University of Freiburg Contact Matthias Teschner Computer Graphics University of Freiburg Georges-Koehler-Allee
More informationLagrangian Data-Driven Reduced Order Modeling of Finite Time Lyapunov Exponents
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 Lagrangian Data-Driven Reduced Order Modeling of Finite Time Lyapunov Exponents Xuping Xie, Peter J. Nolan, Shane D. Ross, and Traian Iliescu Abstract. This
More informationChapter 1 - Basic Equations
2.20 - Marine Hydrodynamics, Sring 2005 Lecture 2 2.20 Marine Hydrodynamics Lecture 2 Chater 1 - Basic Equations 1.1 Descrition of a Flow To define a flow we use either the Lagrangian descrition or the
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationChapter 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 informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationAn adaptive discretization of compressible flow using a multitude of moving Cartesian grids
An adaptive discretization of compressible flow using a multitude of moving Cartesian grids Linhai Qiu, Wenlong Lu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationAnalyzing fluid flows via the
Analyzing fluid flows via the ergodicity defect Sherry E. Scott FFT 2013, Norbert Wiener Center February 22, 2013 Support: ONR-MURI Ocean 3D+1(N00014-11-1-0087) Outline Background & Motivation tion General
More informationConservative Semi-Lagrangian solvers on mapped meshes
Conservative Semi-Lagrangian solvers on mapped meshes M. Mehrenberger, M. Bergot, V. Grandgirard, G. Latu, H. Sellama, E. Sonnendrücker Université de Strasbourg, INRIA Grand Est, Projet CALVI, ANR GYPSI
More informationOverview 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 informationOptimal Segmentation and Understanding of Motion Capture Data
Optimal Segmentation and Understanding of Motion Capture Data Xiang Huang, M.A.Sc Candidate Department of Electrical and Computer Engineering McMaster University Supervisor: Dr. Xiaolin Wu 7 Apr, 2005
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More information1. Introduction 1 2. Mathematical Representation of Robots
1. Introduction 1 1.1 Introduction 1 1.2 Brief History 1 1.3 Types of Robots 7 1.4 Technology of Robots 9 1.5 Basic Principles in Robotics 12 1.6 Notation 15 1.7 Symbolic Computation and Numerical Analysis
More informationLECTURE 8: SMOOTH SUBMANIFOLDS
LECTURE 8: SMOOTH SUBMANIFOLDS 1. Smooth submanifolds Let M be a smooth manifold of dimension n. What object can be called a smooth submanifold of M? (Recall: what is a vector subspace W of a vector space
More informationLagrangian Coherent Structures near a. subtropical jet stream
Lagrangian Coherent Structures near a subtropical jet stream Wenbo Tang School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85283 Manikandan Mathur, George Haller Department
More informationarxiv: v1 [math.na] 21 Feb 2019
A Fully Lagrangian Meshfree Framework for PDEs on Evolving Surfaces Pratik Suchde a,, Jörg Kuhnert a a Fraunhofer ITWM, 67663 Kaiserslautern, Germany arxiv:1902.08107v1 [math.na] 21 Feb 2019 Abstract We
More informationThe 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 information2011 IEEE Visualization Contest Winner: Visualizing Unsteady Vortical Behavior of a Centrifugal Pump
2011 IEEE Visualization Contest Winner: Visualizing Unsteady Vortical Behavior of a Centrifugal Pump Mathias Otto 1, Alexander Kuhn 1, Wito Engelke 1, 2, and Holger Theisel 1 1 University of Magdeburg,
More informationInverse and Implicit functions
CHAPTER 3 Inverse and Implicit functions. Inverse Functions and Coordinate Changes Let U R d be a domain. Theorem. (Inverse function theorem). If ϕ : U R d is differentiable at a and Dϕ a is invertible,
More informationEulerian and Lagrangian Pictures of Mixing
Eulerian and Lagrangian Pictures of Mixing Jean-Luc Thiffeault Department of Mathematics Imperial College London with Steve Childress Courant Institute of Mathematical Sciences New York University http://www.ma.imperial.ac.uk/
More informationComputer Science Senior Thesis
Computer Science Senior Thesis Matt Culbreth University of Colorado May 4, 2005 0.1 Abstract Two methods for the identification of coherent structures in fluid flows are studied for possible combination
More informationA Grid Based Particle Method for Evolution of Open Curves and Surfaces
A Grid Based Particle Method for Evolution of Open Curves and Surfaces Shingyu Leung a,, Hongkai Zhao a a Department of Mathematics, University of California at Irvine, Irvine, CA 92697-3875. Abstract
More informationLONG-TIME SIMULATIONS OF VORTEX SHEET DYNAMICS
LONG-TIME SIMULATIONS OF VORTEX SHEET DYNAMICS BY DANIEL SOLANO 1, LING XU 2, AND SILAS ALBEN 2 1 Rutgers University, New Brunswick, NJ 891, USA. 2 University of Michigan, Ann Arbor, MI 4819, USA. Abstract.
More informationParticle-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 informationDetection of Lagrangian Coherent Structures in 3D Turbulence
Under consideration for publication in J. Fluid Mech. 1 Detection of Lagrangian Coherent Structures in 3D Turbulence By M. A. G R E E N 1, C. W. R O W L E Y 1 A N D G. H A L L E R 2 1 Department of Mechanical
More informationAn adaptive discretization of incompressible flow using a multitude of moving Cartesian grids
An adaptive discretization of incompressible flow using a multitude of moving Cartesian grids R. Elliot English, Linhai Qiu,YueYu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA
More information3D vector fields. Contents. Introduction 3D vector field topology Representation of particle lines. 3D LIC Combining different techniques
3D vector fields Scientific Visualization (Part 9) PD Dr.-Ing. Peter Hastreiter Contents Introduction 3D vector field topology Representation of particle lines Path lines Ribbons Balls Tubes Stream tetrahedra
More informationFast Multipole Method on the GPU
Fast Multipole Method on the GPU with application to the Adaptive Vortex Method University of Bristol, Bristol, United Kingdom. 1 Introduction Particle methods Highly parallel Computational intensive Numerical
More informationCloth Hair. and. soft bodies
Cloth Hair Lesson 11 and soft bodies Lesson 08 Outline Problem definition and motivations Modeling deformable solids with mass-spring model Position based dynamics Modeling cloths with mass-spring model
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More informationSUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS
SUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS M. Hofer and H. Pottmann Institute of Geometry Vienna University of Technology, Vienna, Austria hofer@geometrie.tuwien.ac.at, pottmann@geometrie.tuwien.ac.at
More informationThis version of flowvc is capable of reading in 2D or 3D velocity data and
1 Capabilities This version of flowvc is capable of reading in 2D or 3D velocity data and 1. Computing FTLE fields, 2. Computing tracer/particle trajectories, 3. Interpolating the velocity data onto another
More informationA Lagrangian analysis of a developing and non-developing disturbance observed during the PREDICT experiment
Manuscript prepared for Atmos. Chem. Phys. with version 3.2 of the L A TEX class copernicus.cls. Date: 5 September 11 A Lagrangian analysis of a developing and non-developing disturbance observed during
More informationImplementation of PML in the Depth-oriented Extended Forward Modeling
Implementation of PML in the Depth-oriented Extended Forward Modeling Lei Fu, William W. Symes The Rice Inversion Project (TRIP) April 19, 2013 Lei Fu, William W. Symes (TRIP) PML in Extended modeling
More informationEssay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions
Table of Contents Essay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions and Scaling Analysis 1 About dimensional analysis 4 1.1 Thegoalandtheplan... 4 1.2 Aboutthisessay... 5 2 Models
More informationAn Embedded Boundary Method with Adaptive Mesh Refinements
An Embedded Boundary Method with Adaptive Mesh Refinements Marcos Vanella and Elias Balaras 8 th World Congress on Computational Mechanics, WCCM8 5 th European Congress on Computational Methods in Applied
More informationMOTION OF A LINE SEGMENT WHOSE ENDPOINT PATHS HAVE EQUAL ARC LENGTH. Anton GFRERRER 1 1 University of Technology, Graz, Austria
MOTION OF A LINE SEGMENT WHOSE ENDPOINT PATHS HAVE EQUAL ARC LENGTH Anton GFRERRER 1 1 University of Technology, Graz, Austria Abstract. The following geometric problem originating from an engineering
More informationAn adaptive discretization of incompressible flow using a multitude of moving Cartesian grids
An adaptive discretization of incompressible flow using a multitude of moving Cartesian grids R. Elliot English, Linhai Qiu, Yue Yu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford,
More informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationVolumes of Solids of Revolution Lecture #6 a
Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationCGT 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 informationarxiv:math/ v1 [math.ap] 18 Jul 2006
arxiv:math/6747v [math.ap] 8 Jul 6 A linear dispersive mechanism for numerical error growth: spurious caustics Claire David, Pierre Sagaut, Tapan Sengupta Université Pierre et Marie Curie-Paris 6 Laboratoire
More informationDD2429 Computational Photography :00-19:00
. Examination: DD2429 Computational Photography 202-0-8 4:00-9:00 Each problem gives max 5 points. In order to pass you need about 0-5 points. You are allowed to use the lecture notes and standard list
More informationParticle Tracing Module
Particle Tracing Module Particle Tracing Module Released with version 4.2a in October 2011 Add-on to COMSOL Multiphysics Combines with any COMSOL Multiphysics Module Particle Tracing Particle tracing
More information