(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
|
|
- Rosalyn Bruce
- 5 years ago
- Views:
Transcription
1 Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more Lehrstuhl für Informatik 10 (Systemsimulation Systemsimulation) Universität Erlangen-Nürnberg www10.informatik.uni-erlangen.de February 24,
2 Outline Free-Surface LBM Boundary condition Gas volume pressure Surface tension Curvature computation for surface tension Free-Surface in motion images Further extensions Coupling with moving objects Disjoining pressure for foam simulation Capillary force and contact angles 2
3 Free-Surface LBM Cell types and gas phase model Free surface is two phase, however: Compute only liquid phase Model gas phase by ideal gas equation (no flow but only pressure is modeled) T. Pohl, High Performance Simulation of Free Surface Flows Using the Lattice Boltzmann Method, PhD Thesis, 2008, Pohl pdf N. Thürey, Physically based Animation of Free Surface Flows with the Lattice Boltzmann Method, PhD Thesis, 2007, 3
4 Free-Surface LBM Fraction of fluid volume in Interface cells Interface cells store fluid fraction (similar to VoF) Fluid advection: Mass transfer between interface and liquid id cells Mass is conserved For two interface cells a weighting factor depends on the fill levels 4
5 Outline Free-Surface LBM Boundary condition Gas volume pressure Surface tension Curvature computation for surface tension Free-Surface in motion images Further extensions Coupling with moving objects Disjoining pressure for foam simulation Capillary force and contact angles 5
6 Free-Surface LBM Interface boundary condition Simplification of stress tensor at interface Navier-Stokes for incompressible fluid with stress tensor Treating gas as inert phase with vanishing density reduces stress to sole influence of pressure 6
7 Free-Surface LBM Interface boundary condition: gas pressure Gas volume pressure Initial state with initial gas volume Tracking gas volume changes in all interface cells leads to current volume Pressure is ratio of current volume to initial volume Remarks: Tracking gas volume changes throughout a parallel simulation involves means of all-to-all communication Gas pressure depends on accurate mass tracking in interface cells Supports continuous inflation of bubbles 7
8 Free-Surface LBM Interface boundary condition: surface tension Surface tension Energy balance relates pressure to surface tension: Young-Laplace equation (for 3 dimensions) Takes only mean curvature into account, which is the mean value of the two principal curvatures 8
9 Free-Surface LBM Interface boundary condition in LBM Reconstruction of distribution functions at interface Simplified stress tensor means in LBM: Velocity of fluid and gas phase have to be equal at interface Force of the gas has to be balanced with force by fluid Reconstruct distribution functions pointing in opposite direction of surface normal Density of gas phase is computed from gas pressure 9
10 Outline Free-Surface LBM Boundary condition Gas volume pressure Surface tension Curvature computation for surface tension Free-Surface in motion images Further extensions Coupling with moving objects Disjoining pressure for foam simulation Capillary force and contact angles 10
11 Free-Surface LBM Surface curvature computation Calculating the surface curvature Needs information of a 5x5x5 neighborhood For parallelization: Split in 3 parts: Surface normal Reconstruction of surface point Surface curvature Only one layer has to be communicated for parallelization 11
12 Free-Surface LBM Surface curvature computation Calculating surface normal Treat fluid fraction of cells as scalar field Compute the gradient and its normal by Parker-Youngs approximation Similar to Center-of-Mass method as used in VoF, but with coefficients different from 1 12
13 Free-Surface LBM Surface curvature computation Calculating the surface point For simplification: Surface fraction in a cell is assumed to be planar, fulfilling two requirements: Plane s normal is equal to surface normal Plane bisects cell into two volumes corresponding to fluid fraction of cell 13
14 Free-Surface LBM Surface curvature computation Calculating the surface point (cont.) Direct calculation of offset a is difficult Solve the inverse problem by bisection algorithm Calculate the fluid volume with a Correct a and proceed Calculation of fluid volume difficult in 3D: 14
15 Free-Surface LBM Surface curvature computation Calculating the surface curvature (of 3x3x3 neighborhood s interface cells) Select appropriate subset of neighboring surface points (to avoid degenerated triangles) Local triangulation of chosen surface points and calculation of surface curvature (similar algorithm as described in: G. Taubin. Estimating the Tensor of Curvature of a Surface from a Polyhedral Approximation. In Proceedings of the Fifth International Conference on Computer Vision, pages , 1995) 15
16 Free-Surface LBM Bubble coalescence Bubble coalescence (Bubble merging) If two bubbles touch each other: Add initial and current volume from bubble with higher index to bubble with lower index Delete bubble with higher index Connect all cells belonging to old bubble to new bubble If target bubble merges with other bubble, do process iterativelyti 16
17 Outline Free-Surface LBM Boundary condition Gas volume pressure Surface tension Curvature computation for surface tension Free-Surface in motion images Further extensions Coupling with moving objects Disjoining pressure for foam simulation Capillary force and contact angles 17
18 Free-Surface LBM in motion images Drop falling on free surface 18
19 Free-Surface LBM in motion images Rising and coalescing bubbles 19
20 Free-Surface LBM in motion images Many rising and coalescing bubbles 20
21 Outline Free-Surface LBM Boundary condition Gas volume pressure Surface tension Curvature computation for surface tension Free-Surface in motion images Further extensions Coupling with moving objects Disjoining pressure for foam simulation Capillary force and contact angles 21
22 Free-Surface LBM Extension: Coupling with moving objects Coupling with FSI-LBM New cell conversion rules Adapted force calculation for interface cells Adapted interface boundary condition for particle cells S. Bogner, Simulation of Floating Objects in Free-Surface Flows, Diploma Thesis, 2009, 22
23 Free-Surface LBM Extension: Disjoining pressure Disjoining i i pressure for foam modeling Add pressure component depending on distance to next bubble C. Körner, M. Thies, T. Hofmann, N. Thürey, and U. Rüde, Lattice Boltzmann Model for Free Surface Flow for Modeling Foaming, Journal of Statistical Physics, Vol. 121, Nos. 1/2, 2005, 23
24 Free-Surface LBM Extension: Capillary force and contact angles Different approaches: Define additional pressure depending on contact angle θ Contact Angle Enforce a curvature at triple line by manipulating curvature computation Introduce additional force based on Young s equation S. Donath, to be published 2010, 24
25 Free-Surface LBM Extension: Capillary force and contact angles Drop on Plate Defined volume, prescribed contact angle, measure height h = V ( ) 1 cosϑ 3 π D Θ=10,, h=3.28 Θ=90,, h=15.12 Θ=150,, h=
26 The End Thank you for your attention! 26
27 Lattice Boltzmann Method Boltzmann equation Evolved from Lattice-Gas-Cellular-Automata Lattice, Rules, Local Mesoscopic - works with distributions of particles Boltzmann equation f f + ξ f = + f F r r ( x, p, t) f = ξ F...external forces t p Fulfills the Maxwell distribution Ω Solves the Navier-Stokes equations ξ Ω... particle velocity... collision operator 27
28 Lattice Boltzmann Method BGK collision i model Collision i Operator 1 )... particle velocity (0 (0) f + f = [ f f ] f... equilibrium distribution function t ξ τ τ... relaxation time ξ f ( 0 ) α Equilibrium distribution function = wα ρ[ 1 + e u ( e u ) 2 α + 4 α 4 c 2c 2c Macroscopic quantities ( u u )] ρ = fα ρu = α α e α w a... weighting factors dx c... lattice speed dt u... macroscopic velocity f α 28
29 Lattice Boltzmann Method Discretization of Velocity Space Discretization of particle velocity space t (finite set of discrete velocities) f + ξ f f α α α = 1 (eq f f ) α α [ τ Different discretization schemes Numerical accuracy and stability Computational speed and simplicity ] f α r r ( x, t) = f ( x, ξα, t) r (0) r ( x, t) = f ( x,, t) ( ) f eq α ξ α D3Q15 D3Q19 D3Q27 29
30 Lattice Boltzmann Method Spatial and Temporal Discretization Discretization in space x and dtime t: collision step: streaming step: f α ~ ( eq f ) α fα ( xi, t) = fα ( xi, t) ω [ fα ( xi, t) ~ ( x + e r δ t, t + δt) = f ( x, t) i α α i ( x i, t)] 30
31 Lattice Boltzmann Method Boundary Conditions Boundary conditions at walls No slip / Free slip / Part slip boundary conditions May be enough for qualitative examination (without contact angle modeling) 31
32 Lattice Boltzmann Method Parameterization ti All physical quantities have to be scaled such that the following conditions hold in lattice space: δt * * * * 1 1 = 1; δx = 1; ρ = 1; c s = ; 0.5 < < 1.95; u 3 τ Physical quantities are scaled and dimensionlessi * < 0.1 Kinematic viscosity m υ s 2 : v * = v δt x 6 υ τ = 2 δx 2 * + 1 Velocity m u : s u * = δt u δx Surface Tension kg : 2 s δt ρ δx 2 * σ σ = σ 3 32
33 Parallelization Why parallelize? li Because of the complicated geometry 3.0 x 0.15x 0.16mm³ resolution dx = 0.1µm Domain of 7.2 * cells 425 Bytes per cell Approx. 28 TB memory Subset of 0.5 x 0.1 x 0.1mm³ results in 2 TB of memory Code is parallel, yes. But not suitable for large-scale parallel simulations! Portation to walberla framework. 33
34 walberla Widelyy applicable lattice Boltzmann from Erlangen pp g CFD project based on lattice Boltzmann method Modular software concept Various applications under development: Blood flow Moving particles and agglomerates Charged colloids Free surfaces Species transport Multi-component multi-phase Brownian motion Integration in efficient massive-parallel environment Project involving currently 5 scientists and many diploma theses 34
35 walberla Patch Concept Supports Fast Integration ti of FreeSurfaces Optimized i for Large-Scale Parallelization li and Porous Media 35
36 Outlook: Contact Angles In our model: Fill levels in cells influence normal Normals influence curvature and hence Fill levels influence force exerted by surface tension Idea: Control contact angle at walls by fill level in wall cells Körner et al. (at FAU Material Sciences) successfully showed relation between physical contact angle and fill levels for 2D 36
37 Outlook: Contact Angles Problems and Open Questions Interface modeling Calculations in 3D Which fill levels taken into account? Which cells for curvature computation? 37
38 Outlook: Contact Angles Open Questions Will contact angle induce capillary effects? Capillaries? Drops? Will capillary effects be necessary? Will this contact angle model work with the thin fibres in the substrate? 38
A Contact Angle Model for the Parallel Free Surface Lattice Boltzmann Method in walberla Stefan Donath (stefan.donath@informatik.uni-erlangen.de) Computer Science 10 (System Simulation) University of Erlangen-Nuremberg
More informationSimulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method
Simulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method June 21, 2011 Introduction Free Surface LBM Liquid-Gas-Solid Flows Parallel Computing Examples and More References Fig. Simulation
More informationFree Surface Lattice-Boltzmann fluid simulations. with and without level sets.
ree Surface Lattice-Boltzmann fluid simulations with and without level sets Nils Thürey, Ulrich Rüde University of Erlangen-Nuremberg System Simulation roup Cauerstr. 6, 91054 Erlangen, ermany Email: Nils.Thuerey@cs.fau.de
More informationComputational Fluid Dynamics with the Lattice Boltzmann Method KTH SCI, Stockholm
Computational Fluid Dynamics with the Lattice Boltzmann Method KTH SCI, Stockholm March 17 March 21, 2014 Florian Schornbaum, Martin Bauer, Simon Bogner Chair for System Simulation Friedrich-Alexander-Universität
More informationRealistic Animation of Fluids
Realistic Animation of Fluids p. 1/2 Realistic Animation of Fluids Nick Foster and Dimitri Metaxas Realistic Animation of Fluids p. 2/2 Overview Problem Statement Previous Work Navier-Stokes Equations
More informationThe walberla Framework: Multi-physics Simulations on Heterogeneous Parallel Platforms
The walberla Framework: Multi-physics Simulations on Heterogeneous Parallel Platforms Harald Köstler, Uli Rüde (LSS Erlangen, ruede@cs.fau.de) Lehrstuhl für Simulation Universität Erlangen-Nürnberg www10.informatik.uni-erlangen.de
More informationPerformance Analysis of the Lattice Boltzmann Method on x86-64 Architectures
Performance Analysis of the Lattice Boltzmann Method on x86-64 Architectures Jan Treibig, Simon Hausmann, Ulrich Ruede Zusammenfassung The Lattice Boltzmann method (LBM) is a well established algorithm
More informationSoftware and Performance Engineering for numerical codes on GPU clusters
Software and Performance Engineering for numerical codes on GPU clusters H. Köstler International Workshop of GPU Solutions to Multiscale Problems in Science and Engineering Harbin, China 28.7.2010 2 3
More informationLehrstuhl für Informatik 10 (Systemsimulation)
FRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG INSTITUT FÜR INFORMATIK (MATHEMATISCHE MASCHINEN UND DATENVERARBEITUNG) Lehrstuhl für Informatik 10 (Systemsimulation) On the Resource Requirements of
More informationPerformance and Software-Engineering Considerations for Massively Parallel Simulations
Performance and Software-Engineering Considerations for Massively Parallel Simulations Ulrich Rüde (ruede@cs.fau.de) Ben Bergen, Frank Hülsemann, Christoph Freundl Universität Erlangen-Nürnberg www10.informatik.uni-erlangen.de
More informationHigh Performance Simulation of Free Surface Flows Using the Lattice Boltzmann Method
High Performance Simulation of Free Surface Flows Using the Lattice Boltzmann Method Hochleistungssimulation von Strömungen mit freien Oberflächen basierend auf der Lattice-Boltzmann-Methode Der Technischen
More informationAdarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)
OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final
More informationPrinciples and Goals
Simulations of foam rheology Simon Cox Principles and Goals Use knowledge of foam structure to predict response when experiments are not easy to isolate causes of certain effects to save time in designing
More informationLattice Boltzmann with CUDA
Lattice Boltzmann with CUDA Lan Shi, Li Yi & Liyuan Zhang Hauptseminar: Multicore Architectures and Programming Page 1 Outline Overview of LBM An usage of LBM Algorithm Implementation in CUDA and Optimization
More informationIntroducing a Cache-Oblivious Blocking Approach for the Lattice Boltzmann Method
Introducing a Cache-Oblivious Blocking Approach for the Lattice Boltzmann Method G. Wellein, T. Zeiser, G. Hager HPC Services Regional Computing Center A. Nitsure, K. Iglberger, U. Rüde Chair for System
More informationFree Surface Flows with Moving and Deforming Objects for LBM
Free Surface Flows with Moving and Deforming Objects for LBM Nils Thürey, Klaus Iglberger, Ulrich Rüde Institute for System Simulation (LSS), University of Erlangen-Nuremberg, Cauerstr. 6, 91058 Erlangen,
More informationFluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]
Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationLattice Boltzmann Method: A New Alternative for Loss Estimation Juan C. Díaz Cruz* a, Sergio E. Rodríguez a, Sergio A. Ordóñez b, Felipe Muñoz a a
799 VOL. 31, 2013 CHEMICAL ENGINEERING TRANSACTIONS A publication of The Italian Association of Chemical Engineering Online at: www.aidic.it/cet Guest Editors: Eddy De Rademaeker, Bruno Fabiano, Simberto
More informationA 3D VOF model in cylindrical coordinates
A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved
More informationAn accurate numerical method for computing surface tension forces in CFD codes
1 An accurate numerical method for computing surface tension forces in CFD codes Numerical experiments with surface tension M. Coquerelle 1 S. Glockner 1 1 Numerical Fluid Mechanics (MFN), département
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 informationRealistic Animation of Fluids
1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely
More informationMassively Parallel Phase Field Simulations using HPC Framework walberla
Massively Parallel Phase Field Simulations using HPC Framework walberla SIAM CSE 2015, March 15 th 2015 Martin Bauer, Florian Schornbaum, Christian Godenschwager, Johannes Hötzer, Harald Köstler and Ulrich
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 informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationHigh Scalability of Lattice Boltzmann Simulations with Turbulence Models using Heterogeneous Clusters
SIAM PP 2014 High Scalability of Lattice Boltzmann Simulations with Turbulence Models using Heterogeneous Clusters C. Riesinger, A. Bakhtiari, M. Schreiber Technische Universität München February 20, 2014
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More informationsimulation framework for piecewise regular grids
WALBERLA, an ultra-scalable multiphysics simulation framework for piecewise regular grids ParCo 2015, Edinburgh September 3rd, 2015 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler
More informationCUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata
CUDA Fluid simulation Lattice Boltzmann Models Cellular Automata Please excuse my layout of slides for the remaining part of the talk! Fluid Simulation Navier Stokes equations for incompressible fluids
More informationENERGY-224 Reservoir Simulation Project Report. Ala Alzayer
ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................
More informationPerformance Optimization of a Massively Parallel Phase-Field Method Using the HPC Framework walberla
Performance Optimization of a Massively Parallel Phase-Field Method Using the HPC Framework walberla SIAM PP 2016, April 13 th 2016 Martin Bauer, Florian Schornbaum, Christian Godenschwager, Johannes Hötzer,
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Christian Wollblad Copyright 2017 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of
More informationTechnical Report TR
Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,
More informationSimulating Sinkage & Trim for Planing Boat Hulls. A Fluent Dynamic Mesh 6DOF Tutorial
Simulating Sinkage & Trim for Planing Boat Hulls A Fluent Dynamic Mesh 6DOF Tutorial 1 Introduction Workshop Description This workshop describes how to perform a transient 2DOF simulation of a planing
More informationThe Lattice Boltzmann Method used for fluid flow modeling in hydraulic components
The 15th Scandinavian International Conference on Fluid Power, SICFP 17, June 7-9, 2017, Linköping, Sweden The Lattice Boltzmann Method used for fluid flow modeling in hydraulic components Bernhard Manhartsgruber
More informationWebinar #3 Lattice Boltzmann method for CompBioMed (incl. Palabos)
Webinar series A Centre of Excellence in Computational Biomedicine Webinar #3 Lattice Boltzmann method for CompBioMed (incl. Palabos) 19 March 2018 The webinar will start at 12pm CET / 11am GMT Dr Jonas
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Implicit Surfaces Marching Cubes/Tetras Collision Detection & Response Conservative Bounding Regions backtracking fixing Today Flow Simulations in Graphics Flow
More informationFree Surface Flow Simulations
Free Surface Flow Simulations Hrvoje Jasak h.jasak@wikki.co.uk Wikki Ltd. United Kingdom 11/Jan/2005 Free Surface Flow Simulations p.1/26 Outline Objective Present two numerical modelling approaches for
More informationLattice Boltzmann Methods on the way to exascale
Lattice Boltzmann Methods on the way to exascale Ulrich Rüde (LSS Erlangen, ulrich.ruede@fau.de) Lehrstuhl für Simulation Universität Erlangen-Nürnberg www10.informatik.uni-erlangen.de HIGH PERFORMANCE
More informationFRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG. Lehrstuhl für Informatik 10 (Systemsimulation)
FRIEDRICH-ALEXANDER-UNIVERSITÄT ERLANGEN-NÜRNBERG INSTITUT FÜR INFORMATIK (MATHEMATISCHE MASCHINEN UND DATENVERARBEITUNG) Lehrstuhl für Informatik 10 (Systemsimulation) walberla: Visualization of Fluid
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationNumerical Algorithms on Multi-GPU Architectures
Numerical Algorithms on Multi-GPU Architectures Dr.-Ing. Harald Köstler 2 nd International Workshops on Advances in Computational Mechanics Yokohama, Japan 30.3.2010 2 3 Contents Motivation: Applications
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationInternational Supercomputing Conference 2009
International Supercomputing Conference 2009 Implementation of a Lattice-Boltzmann-Method for Numerical Fluid Mechanics Using the nvidia CUDA Technology E. Riegel, T. Indinger, N.A. Adams Technische Universität
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 information13. Learning Ballistic Movementsof a Robot Arm 212
13. Learning Ballistic Movementsof a Robot Arm 212 13. LEARNING BALLISTIC MOVEMENTS OF A ROBOT ARM 13.1 Problem and Model Approach After a sufficiently long training phase, the network described in the
More informationNUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)
University of West Bohemia» Department of Power System Engineering NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV) Publication was supported by project: Budování excelentního
More informationAuto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial
Auto Injector Syringe A Fluent Dynamic Mesh 1DOF Tutorial 1 2015 ANSYS, Inc. June 26, 2015 Prerequisites This tutorial is written with the assumption that You have attended the Introduction to ANSYS Fluent
More informationLARGE-SCALE FREE-SURFACE FLOW SIMULATION USING LATTICE BOLTZMANN METHOD ON MULTI-GPU CLUSTERS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationTwo-Phase flows on massively parallel multi-gpu clusters
Two-Phase flows on massively parallel multi-gpu clusters Peter Zaspel Michael Griebel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität Bonn Workshop Programming of Heterogeneous
More informationDrop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationHIGH PERFORMANCE COMPUTATION (HPC) FOR THE
HIGH PERFORMANCE COMPUTATION (HPC) FOR THE DEVELOPMENT OF FLUIDIZED BED TECHNOLOGIES FOR BIOMASS GASIFICATION AND CO2 CAPTURE P. Fede, H. Neau, O. Simonin Université de Toulouse; INPT, UPS ; IMFT ; 31400
More informationA Particle Cellular Automata Model for Fluid Simulations
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 36(2), 2009, Pages 35 41 ISSN: 1223-6934 A Particle Cellular Automata Model for Fluid Simulations Costin-Radu Boldea Abstract. A new cellular-automaton
More informationReal-time Thermal Flow Predictions for Data Centers
Real-time Thermal Flow Predictions for Data Centers Using the Lattice Boltzmann Method on Graphics Processing Units for Predicting Thermal Flow in Data Centers Johannes Sjölund Computer Science and Engineering,
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 informationLattice Boltzmann Method for Simulating Turbulent Flows
Lattice Boltzmann Method for Simulating Turbulent Flows by Yusuke Koda A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Applied Science
More informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
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 informationChapter 2 Water and Bubbles
Chapter 2 Water and Bubbles 2.1 Animation of Bubbles in Liquid Abstract This section introduces a new fluid animation technique in which liquid and gas interact with each other, using the example of bubbles
More informationPhysically based Animation of Free Surface Flows with the Lattice Boltzmann Method
Physically based Animation of Free Surface Flows with the Lattice Boltzmann Method Physikalische Animation von Strömungen mit freien Oberflächen mit der Lattice-Boltzmann-Methode Der Technischen Fakultät
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationDynamic Mode Decomposition analysis of flow fields from CFD Simulations
Dynamic Mode Decomposition analysis of flow fields from CFD Simulations Technische Universität München Thomas Indinger Lukas Haag, Daiki Matsumoto, Christoph Niedermeier in collaboration with Agenda Motivation
More informationPreliminary Spray Cooling Simulations Using a Full-Cone Water Spray
39th Dayton-Cincinnati Aerospace Sciences Symposium Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray Murat Dinc Prof. Donald D. Gray (advisor), Prof. John M. Kuhlman, Nicholas L. Hillen,
More informationThe gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last
Parallel Implementation of Gas-Kinetic BGK Scheme on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer Middle East Technical
More informationSailfish: Lattice Boltzmann Fluid Simulations with GPUs and Python
Sailfish: Lattice Boltzmann Fluid Simulations with GPUs and Python Micha l Januszewski Institute of Physics University of Silesia in Katowice, Poland Google GTC 2012 M. Januszewski (IoP, US) Sailfish:
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationMarker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 1, pp. 1-41 Commun. Comput. Phys. January 2009 REVIEW ARTICLE Marker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
More informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationInvestigation of the critical submergence at pump intakes based on multiphase CFD calculations
Advances in Fluid Mechanics X 143 Investigation of the critical submergence at pump intakes based on multiphase CFD calculations P. Pandazis & F. Blömeling TÜV NORD SysTec GmbH and Co. KG, Germany Abstract
More informationDevelopment of Hybrid Fluid Jet / Float Polishing Process
COMSOL Conference - Tokyo 2013 Development of Hybrid Fluid Jet / Float Polishing Process A. Beaucamp, Y. Namba Dept. of Mechanical Engineering, Chubu University, Japan Zeeko LTD, United Kingdom Research
More information1. Carlos A. Felippa, Introduction to Finite Element Methods,
Chapter Finite Element Methods In this chapter we will consider how one can model the deformation of solid objects under the influence of external (and possibly internal) forces. As we shall see, the coupled
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationSimulation of moving Particles in 3D with the Lattice Boltzmann Method
Simulation of moving Particles in 3D with the Lattice Boltzmann Method, Nils Thürey, Christian Feichtinger, Hans-Joachim Schmid Chair for System Simulation University Erlangen/Nuremberg Chair for Particle
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationMath 690N - Final Report
Math 690N - Final Report Yuanhong Li May 05, 008 Accurate tracking of a discontinuous, thin and evolving turbulent flame front has been a challenging subject in modelling a premixed turbulent combustion.
More informationLATTICE-BOLTZMANN AND COMPUTATIONAL FLUID DYNAMICS
LATTICE-BOLTZMANN AND COMPUTATIONAL FLUID DYNAMICS NAVIER-STOKES EQUATIONS u t + u u + 1 ρ p = Ԧg + ν u u=0 WHAT IS COMPUTATIONAL FLUID DYNAMICS? Branch of Fluid Dynamics which uses computer power to approximate
More informationCFD Simulation of a dry Scroll Vacuum Pump including Leakage Flows
CFD Simulation of a dry Scroll Vacuum Pump including Leakage Flows Jan Hesse, Rainer Andres CFX Berlin Software GmbH, Berlin, Germany 1 Introduction Numerical simulation results of a dry scroll vacuum
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More informationMeasurement Techniques. Digital Particle Image Velocimetry
Measurement Techniques Digital Particle Image Velocimetry Heat and Mass Transfer Laboratory (LTCM) Sepideh Khodaparast Marco Milan Navid Borhani 1 Content m Introduction m Particle Image Velocimetry features
More informationCoupled Analysis of FSI
Coupled Analysis of FSI Qin Yin Fan Oct. 11, 2008 Important Key Words Fluid Structure Interface = FSI Computational Fluid Dynamics = CFD Pressure Displacement Analysis = PDA Thermal Stress Analysis = TSA
More informationNumerical Modeling Study for Fish Screen at River Intake Channel ; PH (505) ; FAX (505) ;
Numerical Modeling Study for Fish Screen at River Intake Channel Jungseok Ho 1, Leslie Hanna 2, Brent Mefford 3, and Julie Coonrod 4 1 Department of Civil Engineering, University of New Mexico, Albuquerque,
More informationMONTEREY, CALIFORNIA DISSERTATION LATTICE BOLTZMANN METHODS FOR FLUID STRUCTURE INTERACTION. Stuart R. Blair. September Dissertation Supervisor:
MONTEREY, CALIFORNIA DISSERTATION LATTICE BOLTZMANN METHODS FOR FLUID STRUCTURE INTERACTION by Stuart R. Blair September 2012 Dissertation Supervisor: Young Kwon Approved for public release; distribution
More informationLeast Square Curvature Calculation Method for VOF Schemes
ILASS Americas 28th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2016 Least Square Curvature Calculation Method for VOF Schemes E. Cauble, M. Owkes Department of Mechanical
More informationCFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe
CFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe Adib Zulhilmi Mohd Alias, a, Jaswar Koto, a,b,* and Yasser Mohamed Ahmed, a a) Department of Aeronautics, Automotive and Ocean
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 informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationNumerical Simulation of Fuel Filling with Volume of Fluid
Numerical Simulation of Fuel Filling with Volume of Fluid Master of Science Thesis [Innovative and Sustainable Chemical Engineering] Kristoffer Johansson Department of Chemistry and Bioscience Division
More informationREGULARITIES OF INK-WATER BALANCE STABILITY IN OFFSET PRINTING
REGULARITIES OF INK-WATER BALANCE STABILITY IN OFFSET PRINTING Olha Blagodir o.blagodir@kpi.ua Каteryna Zolotukhina k.zolotukhina@kpi.ua Bogdan Kushlyk b.kushlyk@kpi.ua Olena Velychko o.velychko@kpi.ua
More information