An added mass partitioned algorithm for rigid bodies and incompressible flows
|
|
- Amice Ramsey
- 5 years ago
- Views:
Transcription
1 An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216
2 Collaborators Bill Henshaw, Don Schwendeman, Qi Tang Department of Mathematical Sciences Rensselaer Polytechnic Institute Support Department of Energy Office of Advanced Scientific Computing Research Applied Mathematics Program NNSA Rensselaer Polytechnic Institute
3 In recent work we have focused on large deformation and/or displacement FSI using partitioned solvers Component solvers remain independent can use existing solvers no need to solve (or precondition) a coupled implicit system Can naturally take advantage of disparate time scales e.g. mixing implicit and explicit integration High levels of algorithmic concurrency maps well to modern and emerging computers
4 In recent work we have focused on large deformation and/or displacement FSI using partitioned solvers fluid solver interface solid solver Component solvers remain independent can use existing solvers no need to solve (or precondition) a coupled implicit system Can naturally take advantage of disparate time scales e.g. mixing implicit and explicit integration High levels of algorithmic concurrency maps well to modern and emerging computers
5 Traditional partitioned schemes have suffered from added-mass instabilities for solids which are sufficiently light when compared to the fluid Traditional partitioned FSI algorithms (Cirak, et. al. 27, Bungartz and Schafer 26) 1. advance fluid (using interface velocity/displacement from the solid) 2. advance solid (apply fluid forces to the solid) 3. possibly iterate with under-relaxation to convergence t (n) fluid force displ./vel. structure t (n+1) diagram from Keyes et. al. 212 underrelaxation interface quasi Newton Some analysis of added-mass instabilities can be found in the literature, for example Causin, Grebeau, and Nobile, 25 (stability with relaxation) Gretarsson, Kwatra, and Fedkiw 211 (semi-monolithic formulations)
6 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum force Body simply moves according to Newton s laws of motion
7 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum force Body simply moves according to Newton s laws of motion
8 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum in a fluid force force Body simply moves according to Newton s laws of motion
9 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum in a fluid fluid must be displaced force force Body simply moves according to Newton s laws of motion
10 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum in a fluid force force fluid must be entrained Body simply moves according to Newton s laws of motion
11 The origin of added-mass instabilities is that the effect of displaced fluid is not appropriately accounted for in the numerical algorithms in a vacuum in a fluid fluid contributing to added mass force force Body simply moves according to Newton s laws of motion Body must displace and entrain fluid to move and therefore appears more massive than in vacuum... the so called added mass
12 As a concrete motivating example consider a rising rigid body in counterflow Incompressible Navier-Stokes Light rigid body flow
13 This case has both strong added-mass and added-damping effects Added Mass Added Damping M a n p t m b G I z F Added mass relates to increased apparent mass owing to fluid displacement i.e. the local geometry occupied by solid changes Added damping relates to increased apparent inertia owing to viscous fluid drag i.e. the local geometry remains fixed
14 To understand these effects in isolation we derive extremely simple models by localizing and linearizing the problem near the interface y-translations relate to added-mass effects x-translations relate to added-damping effects
15 The resulting model is used to motivate our new AMP algorithms and discuss the performance of traditional partitioned (TP) schemes y = H fluid: y = interface: b rigid body: b y = H x = x = L Fluid: Rigid body: m b a u = m b a v + rp = v, x 2, r v =, x 2, Z L Z (x,,t) dx + g u(t), p(x,,t) dx + g v (t), Interface: v(x,,t)=v b (t), x 2 [,L], 2
16 An added-mass model problem is derived by considering vertical motions y = H y = y = H x = fluid: interface: b rigid body: b x = @y =, dv b m b dt = >: v(,t)=v b, Z L pdx, p(h, t) =p H (t), y 2 (,H), y 2 (,H), The AMP scheme matches the vertical accelerations at y= = dv b dt = a And applies a generalized Robin condition to the fluid pressure equation = m b a = Z L pdx
17 The resulting AMP scheme is stable for any finite mass, while the traditional scheme suffers y = H y = y = H x = fluid: interface: b rigid body: b x = @y =, dv b m b dt = >: v(,t)=v b, Z L pdx, p(h, t) =p H (t), y 2 (,H), y 2 (,H), The added mass for this problem is easily identified as M a = LH Thm: The 2nd order accurate AMP scheme for the added-mass model problem is stable provided m b + M a is bounded away from zero. Thm: The 2nd order traditional partitioned scheme is stable if and only if m b >M a
18 An added-damping model problem is derived by considering horizontal motions y = H y = fluid: interface: b rigid body: b = u 2, Z L du b m b dt = µ >: u(,t)=u b (,t) u(h, t) =u H (t), y 2 (,H), y = H x = x = L The AMP scheme uses the exact solution to form the discrete approximation µ Z (,t) Du b D µ Z L 1 e y = y p t/2 here is a ratio of viscous length scales, and D is an added-damping coefficient (in 3D these become added-damping tensors)
19 The AMP scheme with extra velocity projection is stable even for massless bodies y=h fluid: interface: b y= AMP VC, A=4.5, β=.25 1 n.5 Re(ξ ) 1 Re(A ) (scaled) rigid body: b j y= H > > = µ 2, y 2 (, H), > > Z L > mb =µ (, t) dx, > > AMP VC, A= i, β=.31 AMP VC, A= i, β= 6. dt > 1 : u(, t) = ub (t),iv u(h, t) = uh (t), II.5 x =L.5 Figure 3: The geometry for the rectangular geometry FSI.5 model problems. I Re(An) (scaled) j, n Re(ξj) 5 j, n 1 Re(An) (scaled) 1 Re(ξj) ic equation for the fluid pressure can be derived by taking the divergence of the momentum 1) and using the continuity equation in (2) to give p=, xstable 2. region (9) ty-pressure form of the Stokes equations, the momentum equation for the velocity in (1) is used place equation for the pressure in (9) instead of the divergence equation in (2). II n consistency with the original velocity-divergence form of the equations, zero divergence is an additional boundary condition, i.e. r v = for x see [23] for additional details. In n to follow we will use a fractional-step method based on the velocity-pressure form of the fluid define partitioned FSI algorithms. In Section 3 a traditional partitioned scheme, referred to as IV added mass partitioned scheme algorithm, is described. Then in Section 4 we present our new e referred to as the AMP-RB algorithm. The important di erence in the two approaches lies in of the fractional-step method for the fluid with the equations governing the motion of the body. I Stable region onal Partitioned/Rigid-Body algorithm TP-RB n our discussion of partitioned algorithms for the model FSI problem by considering a traditional 5 j, n 1
20 The AMP-RB scheme is implemented in Overture and is found to be stable against both added-mass and added damping instabilities without iteration b =.1 p t = p t = p t = G
21 .157 a G b.66 ab G8 ab G16 ab G32 By comparison, the traditional partitioned scheme requires ~85 under-relaxed iterations to provide comparable results t.7 ab G4 a G8.72 ab G16.74 b.5 a G32 b.1 2 b = t ab G4 ab G8 ab G16 a G32 b t Figure 11: Light rising body. Left: vertical acceleration of the light rising body. Right: zoom near t = 4.75 the perturbations to the accelerations caused by the changes in overlapping grid interpolation points as moves. the surface integrals defining the added-damping tensors (57)-(6), using the discrete geometry defining the rounded rectangle. For grid G (4), these computed added-damping coefficients for the rounded rectangle were found to be µ µ µ vv vv!! D , D22.917, D22.379, n n n p t = 4. p t are t = 2. = 1. withp other entries in the tensors being approximately zero. These above entries nearly equal to the actual coefficients for a square-cornered solid rectangle (??) of vv D11 = 2w µ µ =2, n n.141 vv D22 = 2h.444 µ µ =, n n.26!! D22 = wh(w + h)/2.719 µ 3 µ µ = =.375. n 8 n n Risplay ingabody, b =.1 For this problem, added-damping e ects can critical ρrole, especially on coarser grids p when the time.7 step is larger. Recall that added-damping is proportional to µ t/ n, where n = t/2 and thus p yb (4) added-damping e ects increase for larger values of t. Figure 13 shows that the simulation on grid Gfb.6 vb e ects when the added-damping coefficient d is chosen too small ( d =.5 is unstable due to added-damping in this case). The instability appears primarily in the rotation of the body, the time history of the angular ab acceleration,! 3, shows.5 a large high-frequency oscillation. This instability seems to become saturated due to yb TP SI counter-acting pressure forces (i.e. added-mass e ects). vb TP SI Figure 12 compares results from the AMP-RB scheme versus the TP-SI on grid G (4). The TP-SI scheme.4 required on average 85 sub-iteration per step with a relaxation parameter of! =.25. The results are ab TP SI nearly indistinguishable, indicating that the AMP scheme can achieve essentially the same solution as the.3 without sub-iterations. traditional scheme but Figure 9: Rising body, streamlines and contours of the pressure for b =.1 at times t = 2, t = 4 and t = 1 computed using grid Grid G (16) Rising body, ρ b =.1 yb G4 yb G8 yb G16.5 Ris ing body, ρ b = Ris ing body, ρ b =.1.4 Figure 8: yb G4 Ris ing body, ρ b =.1 6 t.14 yb G Figure 12: A comparison of the AMP-RB scheme versus the TP-SI for a light body. Results for grid G (4).1t = Rising this body: difficult Coarse composite grid Gthe, attp-si times 2, t = 4 and t = 1. Aon gridaverage, is located on top problem, scheme required, 85thesub-iterations per step. The results ar yb G8 yb G16 yb G v G4
22 A more challenging case is that of a light cylinder rising in counterflow The AMP-RB scheme is again stable without any iteration Traditional partitioned scheme needs ~2 sub iterations to stabilize b =.1 y x
23 A more challenging case is that of a light cylinder rising in counterflow b =.1 flow
24 Extensions to multiple bodies presents no particular challenges, and the scheme is robust even in the presence of both light and heavy bodies
25 Extensions to multiple bodies presents no particular challenges, and the scheme is robust even in the presence of both light and heavy bodies
26 Summary Careful analysis of simplified model problems motivates our stable and 2nd order accurate AMP-RB scheme for incompressible flows Stability analysis in simple geometries shows excellent stability properties even in the uniterated form Implementation within Overture illustrates the utility of the approach for both light and heavy bodies Future Work implement in 3D investigate the alternate added-mass potential formulation New FSI regimes (incompressible/incompressible, compressible/beams)
Solving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More 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 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 informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More 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 informationCharacteristic Aspects of SPH Solutions
Characteristic Aspects of SPH Solutions for Free Surface Problems: Source and Possible Treatment of High Frequency Numerical Oscillations of Local Loads. A. Colagrossi*, D. Le Touzé & G.Colicchio* *INSEAN
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationInteraction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 2015) Interaction of Fluid Simulation Based on PhysX Physics Engine Huibai Wang, Jianfei Wan, Fengquan Zhang College
More 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 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 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 informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More 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 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 informationPartitioned strongly coupled Fluid-Structure Interaction
Partitioned strongly coupled Fluid-Structure Interaction 7 th OpenFOAM Workshop Darmstadt, Germany Manuel Kosel * 1 and Ulrich Heck 2 1 Center for Computational Engineering Science, RWTH Aachen University,
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 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 informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More 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 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 informationA Review of the Implicit Motion Solver Algorithm in OpenFOAM to Simulate a Heaving Buoy
A Review of the Implicit Motion Solver Algorithm in OpenFOAM to Simulate a Heaving Buoy Devolder, B., Schmitt, P., Rauwoens, P., Elsaesser, B., & Troch, P. (2015). A Review of the Implicit Motion Solver
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
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 informationA THREE-DIMENSIONAL ADAPTIVE WAVELET METHOD FOR FLUID STRUCTURE INTERACTION
A THREE-DIMENSIONAL ADAPTIVE WAVELET METHOD FOR FLUID STRUCTURE INTERACTION N.K.-R. Kevlahan 1, O.V. Vasilyev 2, D. Goldstein 2, and A. Jay 1,3 kevlahan@mcmaster.ca 1 Department of Mathematics & Statistics,
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 informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationIsogeometric Analysis of Fluid-Structure Interaction
Isogeometric Analysis of Fluid-Structure Interaction Y. Bazilevs, V.M. Calo, T.J.R. Hughes Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA e-mail: {bazily,victor,hughes}@ices.utexas.edu
More informationOpenFOAM and Third Party Structural Solver for Fluid Structure Interaction Simulations
OpenFOAM and Third Party Structural Solver for Fluid Structure Interaction Simulations Robert L. Campbell rlc138@arl.psu.edu Fluids and Structural Mechanics Office Applied Research Laboratory The Pennsylvania
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 informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
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 informationT6: Position-Based Simulation Methods in Computer Graphics. Jan Bender Miles Macklin Matthias Müller
T6: Position-Based Simulation Methods in Computer Graphics Jan Bender Miles Macklin Matthias Müller Jan Bender Organizer Professor at the Visual Computing Institute at Aachen University Research topics
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Optional Reading for Last Time: Spring-Mass Systems Numerical Integration (Euler, Midpoint, Runge-Kutta) Modeling string, hair, & cloth HW2: Cloth & Fluid Simulation
More informationA singular value decomposition based generalized finite difference method for fluid solid interaction problems
Fluid Structure Interaction V 25 A singular value decomposition based generalized finite difference method for fluid solid interaction problems P. Yu, K. S. Yeo, X. Y. Wang & S. J. Ang Department of Mechanical
More informationRigid Body Dynamics, Collision Response, & Deformation
Rigid Body Dynamics, Collision Response, & Deformation Pop Worksheet! Teams of 2. SOMEONE YOU HAVEN T ALREADY WORKED WITH What are the horizontal and face velocities after 1, 2, and many iterations of
More informationFluid-Structure Interaction Modeling of High-Aspect Ratio Nuclear Fuel Plates using COMSOL
Fluid-Structure Interaction Modeling of High-Aspect Ratio Nuclear Fuel Plates using COMSOL Franklin Curtis 1 Kivanc Ekici 1 James Freels 2 1 University of Tennessee, Knoxville, TN 2 Oak Ridge National
More informationOverview and Recent Developments of Dynamic Mesh Capabilities
Overview and Recent Developments of Dynamic Mesh Capabilities Henrik Rusche and Hrvoje Jasak h.rusche@wikki-gmbh.de and h.jasak@wikki.co.uk Wikki Gmbh, Germany Wikki Ltd, United Kingdom 6th OpenFOAM Workshop,
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 informationProblem description. The FCBI-C element is used in the fluid part of the model.
Problem description This tutorial illustrates the use of ADINA for analyzing the fluid-structure interaction (FSI) behavior of a flexible splitter behind a 2D cylinder and the surrounding fluid in a channel.
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
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 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 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 informationDNV GL s 16th Technology Week
OIL & GAS DNV GL s 16th Technology Week Advanced Simulation for Offshore Application: Application of CFD for Computing VIM of Floating Structures 1 SAFER, SMARTER, GREENER OUTLINE Introduction Elements
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
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 information3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic
ISCS 2013 Selected Papers Dam-break effect on a Solid Wall 1 3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic Suprijadi a,b, F. Faizal b, C.F. Naa a and A.Trisnawan
More informationComputational Fluid Dynamics (CFD) using Graphics Processing Units
Computational Fluid Dynamics (CFD) using Graphics Processing Units Aaron F. Shinn Mechanical Science and Engineering Dept., UIUC Accelerators for Science and Engineering Applications: GPUs and Multicores
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 informationThree-dimensional numerical simulations of flapping wings at low Reynolds numbers
Three-dimensional numerical simulations of flapping wings at low Reynolds numbers OpenFOAM Workshop, Zagreb Frank Bos, Bas van Oudheusden, Hester Bijl 7-9 June 2007 1/22 Delft University of Technology
More informationFLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016
FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions
More informationFluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching
Fluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching Paul Groenenboom ESI Group Delft, Netherlands Martin Siemann German Aerospace Center (DLR) Stuttgart, Germany
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationIntroduction to the immersed boundary method
Introduction to the immersed boundary method Motivation. Hydrodynamics and boundary conditions The incompressible Navier-Stokes equations, ( ) u ρ + (u )u = p + ρν 2 u + f, () t are partial differential
More informationPost-processing utilities in Elmer
Post-processing utilities in Elmer Peter Råback ElmerTeam CSC IT Center for Science PATC course on parallel workflows Stockholm, 4-6.12.2013 Alternative postprocessors for Elmer Open source ElmerPost Postprocessor
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
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 informationExplicit and Implicit Coupling Strategies for Overset Grids. Jörg Brunswig, Manuel Manzke, Thomas Rung
Explicit and Implicit Coupling Strategies for s Outline FreSCo+ Grid Coupling Interpolation Schemes Implementation Mass Conservation Examples Lid-driven Cavity Flow Cylinder in a Channel Oscillating Cylinder
More informationPhysically Based Simulation
CSCI 480 Computer Graphics Lecture 21 Physically Based Simulation April 11, 2011 Jernej Barbic University of Southern California http://www-bcf.usc.edu/~jbarbic/cs480-s11/ Examples Particle Systems Numerical
More informationCFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence
CFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence Kavya H.P, Banjara Kotresha 2, Kishan Naik 3 Dept. of Studies in Mechanical Engineering, University BDT College of Engineering,
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Pop Worksheet! Teams of 2. Hand in to Jeramey after we discuss. Sketch the first few frames of a 2D explicit Euler mass-spring simulation for a 2x3 cloth network
More 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 informationPhysically Based Simulation
CSCI 420 Computer Graphics Lecture 21 Physically Based Simulation Examples Particle Systems Numerical Integration Cloth Simulation [Angel Ch. 9] Jernej Barbic University of Southern California 1 Physics
More informationIMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY DEFINED WALLS FOR PARTICLE BASED FLUID SIMULATION
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK IMPROVED WALL BOUNDARY CONDITIONS WITH IMPLICITLY
More informationFluid Structure Interaction Evaluation of two coupling techniques
Technical report, IDE1135, June 2011 Fluid Structure Interaction Evaluation of two coupling techniques Master s Thesis in Computational Science and Engineering Christoffer Andersson, Daniel Ahl School
More informationNumerical Estimation and Validation of Shallow Draft Effect on Roll Damping
The 14 th International Ship Stability Workshop (ISSW), 29 th September- 1 st October 2014, Kuala Lumpur, Malaysia Numerical Estimation and Validation of Shallow Draft Effect on Roll Damping Toru Katayama
More informationChapter 3: Computer Animation Reminder: Descriptive animation. Procedural animation : Examples. Towards methods that generate motion?
Chapter 3 : Computer Animation (continued) Chapter 3: Computer Animation Reminder: Descriptive animation Describes a single motion, with manual control Ex: direct kinematics with key-frames, inverse kinematics
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 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 informationFluidlrigid body interaction in complex industrial flows
Fluidlrigid body interaction in complex industrial flows D. ~bouri', A. parry1 & A. ~arndouni~ 1 Schlumberger - Riboud Product Center, Clamart, France 2 University of La Rochelle, LEPTAB, La Rochelle,
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationMarine Hydrodynamics Solver in OpenFOAM
Marine Hydrodynamics Solver in OpenFOAM p. 1/14 Marine Hydrodynamics Solver in OpenFOAM Hrvoje Jasak and Henrik Rusche h.jasak@wikki.co.uk, h.rusche@wikki.co.uk Wikki, United Kingdom and Germany 4 December
More informationAnalysis, extensions and applications of the Finite-Volume Particle Method (FVPM) PN-II-RU-TE Synthesis of the technical report -
Analysis, extensions and applications of the Finite-Volume Particle Method (FVPM) PN-II-RU-TE-2011-3-0256 - Synthesis of the technical report - Phase 1: Preparation phase Authors: Delia Teleaga, Eliza
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationImmersed Boundary Method and Chimera Method applied to Fluid-
The numericsacademy Fixed Colloquium IBM on Moving Immersed IBM Boundary Applications Methods : Conclusion Current Status and Future Research Directions 15-17 June 2009, Academy Building, Amsterdam, the
More informationCFD-1. Introduction: What is CFD? T. J. Craft. Msc CFD-1. CFD: Computational Fluid Dynamics
School of Mechanical Aerospace and Civil Engineering CFD-1 T. J. Craft George Begg Building, C41 Msc CFD-1 Reading: J. Ferziger, M. Peric, Computational Methods for Fluid Dynamics H.K. Versteeg, W. Malalasekara,
More informationEffect of initial turbulence intensity and velocity profile on liquid jets for IFE beamline protection
Effect of initial turbulence intensity and velocity profile on liquid jets for IFE beamline protection A. Konkachbaev, N.B. Morley and M. A. Abdou Mechanical and Aerospace Engineering Department, UCLA
More informationProgram: Advanced Certificate Program
Program: Advanced Certificate Program Course: CFD-Vehicle Aerodynamics Directorate of Training and Lifelong Learning #470-P, Peenya Industrial Area, 4th Phase Peenya, Bengaluru 560 058 www.msruas.ac.in
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 informationANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step
ANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Governing Equation Start-Up Geometry
More informationComputing Nearly Singular Solutions Using Pseudo-Spectral Methods
Computing Nearly Singular Solutions Using Pseudo-Spectral Methods Thomas Y. Hou Ruo Li January 9, 2007 Abstract In this paper, we investigate the performance of pseudo-spectral methods in computing nearly
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
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 informationNumerical Simulations of Fluid-Structure Interaction Problems using MpCCI
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI François Thirifay and Philippe Geuzaine CENAERO, Avenue Jean Mermoz 30, B-6041 Gosselies, Belgium Abstract. This paper reports
More informationNUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE
NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE Jungseok Ho 1, Hong Koo Yeo 2, Julie Coonrod 3, and Won-Sik Ahn 4 1 Research Assistant Professor, Dept. of Civil Engineering,
More 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 informationMiddle East Technical University Mechanical Engineering Department ME 485 CFD with Finite Volume Method Fall 2017 (Dr. Sert)
Middle East Technical University Mechanical Engineering Department ME 485 CFD with Finite Volume Method Fall 2017 (Dr. Sert) ANSYS Fluent Tutorial Developing Laminar Flow in a 2D Channel 1 How to use This
More informationEvaluation of hydrodynamic coefficients on riser floaters using CFD
Evaluation of hydrodynamic coefficients on riser floaters using CFD Erico Santos, Pedro Mendes, Bruno Luna PETROBRAS (CENPES/PDEP/TDUT) Ricardo Damian ESSS AGENDA MOTIVATION PROCEADURE OVERVIEW RISER GLOBAL
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationSystem Level Cooling, Fatigue, and Durability. Co-Simulation. Stuart A. Walker, Ph.D.
System Level Cooling, Fatigue, and Durability Analysis via Multiphysics Co-Simulation Stuart A. Walker, Ph.D. swalker@altair.com Outline Motivation Presentation of process Presentation of tools Presentation
More informationBest Practices Workshop: Overset Meshing
Best Practices Workshop: Overset Meshing Overview Introduction to Overset Meshes Range of Application Workflow Demonstrations and Best Practices What are Overset Meshes? Overset meshes are also known as
More informationSimulation of Overhead Crane Wire Ropes Utilizing LS-DYNA
Simulation of Overhead Crane Wire Ropes Utilizing LS-DYNA Andrew Smyth, P.E. LPI, Inc., New York, NY, USA Abstract Overhead crane wire ropes utilized within manufacturing plants are subject to extensive
More information