Immersed Boundary Method and Chimera Method applied to Fluid-

Size: px
Start display at page:

Download "Immersed Boundary Method and Chimera Method applied to Fluid-"

Transcription

1 The numericsacademy Fixed Colloquium IBM on Moving Immersed IBM Boundary Applications Methods : Conclusion Current Status and Future Research Directions June 2009, Academy Building, Amsterdam, the Netherlands Immersed Boundary Method and Chimera Method applied to Fluid- Structure Interactions Y. Hoarau, V. Novacek, T. Deloze, A. Makradi Institut de Mécanique des Fluides et des Solides, Université de Strasbourg / CNRS Groupe Instabilité, Turbulence, Diphasique 2 rue Boussingault Strasbourg France hoarau@imfs.u-strasbg.fr

2 Motivations «Complex» geometries, FSI and structured grids?

3 NSMB solver (Navier-Stokes MultiBloc) - Navier-Stokes compressible and incompressible solver - Finite volumes, structured grids, multi-blocks and parallel (MPI) - Turbulence modelling : algebraic models, Spalart-Allmaras, two equations models, EARSM, NLEVM, RSM, OES, DES, DDES, WMLES, LES - ALE : moving grid, deforming grid (remeshing) - Chimera method (overlapping grid) - Immersed Boundary Method Compressible Version - Spatial schemes : centered with artificial dissipation (Jameson, Martinelli, Matrix), upwind (Roe, AUSM, AUSMP, AUSMDV), Harten, Van Leer, Riemman - Temporal scheme : Runge-Kutta, LU-SGS semi implicit, dual time stepping - Multi-grid - Preconditioning (Weiss Smith, Choi/Merkle, artificial compressibility) - Chemistry : Equilibrium/Non-Equilibrium i E i Air, Argon and N2 Chemistry for hypersonic and plasma flows application, diffusion flame combustion model, coupling with Chemkin II - Patch and sliding grids, AMR 3

4 NSMB solver (Navier-Stokes MultiBloc) Incompressible version - Pressure-velocity coupling : SIMPLE, SIMPLEC, PISO, Braza - Rhie & Chow stabilisation - Spatial schemes : centered (2 nd et 4 th order), upwind (1 st, 2 nd and 3 rd order) - Temporal scheme : 1 st,2 nd and 3 rd order backward Euler, Cranck-Nicolson - (Multi-grid) - Thermal convection - SPARSKIT library to solve sparse matrix 4

5 NSMB solver (Navier-Stokes MultiBloc) The Immersed Boundary Method in NSMB - Discrete forcing : Ghost-cell method - 2 layers of ghost-cells are interpolated - Interpolation based on inverse distance - Fully parallelised - No correction at the IB for the incompressibility constraint (for the moment!) Fluid domain Solid domain 5

6 NSMB solver (Navier-Stokes MultiBloc) The Chimera Method in NSMB + = 6

7 NSMB solver (Navier-Stokes MultiBloc) Limit of the sphere mesh Buffer layer cells Buffer layer cells Interpolated cells Interpolated cells ignored cells ignored cells Limit of the wall dominant layer near the tube wall mesh near wall dominant layer Couche dominante proche paroi 7

8 Validation of fixed IBM Trans-sonic/Supersonic sonic/supersonic flows around a NACA0012 2D steady cylinder Sphere 2D unsteady cylinder NACA0012 at 20 of incidence, Re=800 8

9 Validation of fixed IBM NACA0012, M=0.85, α=2, Re=100 Cx Cx IBM 4 th order Central 0, , rd order Roe 0, , rd order AUSM+ 0, , rd order Van Leer 0, ,

10 Validation of fixed IBM NACA0012, M=2, α=5, Re=500 3 rd order Roe scheme IBM : Cx = Cz = Body fitted grid Cx = 0.33 Cz =

11 Validation of fixed IBM 2D steady cylinder SIMPLE, 4 th order central scheme Re=20 Re=30 11

12 Validation of fixed IBM 2D steady cylinder SIMPLE, 4 th order central scheme 12

13 Validation of fixed IBM Steady sphere SIMPLE, 4 th order central scheme 13

14 Validation of fixed IBM Steady sphere SIMPLE, 4 th order central scheme 14

15 Validation of fixed IBM Steady sphere with thermal convection SIMPLE, 4 th order central scheme 15

16 Validation of fixed IBM Unsteady 2D cylinder 4 th order Braza central scheme Re=300 16

17 Validation of fixed IBM Unsteady 2D cylinder 4 th order Braza central scheme 17

18 Validation of fixed IBM NACA0012, α =20, Re=800 4 th order Braza central scheme Nombre de Strouhal NSMB IBM Icare 0.55 Ventikos 0.52 Pulliam 0.5 Bouhadji

19 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 Forced oscillation motion : z=z z 0 +Asin( A.sin(ω t) F = f osc / f 0 Lock-in zone by Koopman, JFM 1967 Comparison with Placzek & al, Comp. & Flds

20 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 (A,F)=(0.25,1) (A,F)=(0.25,0.9) 20

21 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 (A,F)=(0.25,1.5) 21

22 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 22

23 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 (A,F)=(1,0.9) Mode 2P 23

24 Validation of moving IBM Unsteady oscillating 2D cylinder 4 th order Braza central scheme, Re=100 (A,F)=(1.25,1.5) Mode P+S 24

25 Validation of moving IBM Z = z d Unsteady free oscillating 2D cylinder 4 th order Braza central scheme, Re=100 2 d Z dz * * * * m. + b. + k. Z = C ( t ) 2 Z dt dt U m b * * * *, t = t, m =, b =, k = 1 1 d 2 ρ U d 2 ρ U d 1 2 k ρ U 2 k* m* A f* Cx Shiels & al NSMB Shiels & al NSMB Shiels & al NSMB

26 Laminar tube array 26

27 Icicle melting. Q dt = λ ds dn dm = dt. Q L 27

28 Turek FSI test-case 28

29 Conclusion - Chimera and Immersed boundary methods successfully implemented in NSMB - Very powerful tools to handle complex geometries and FSI - IBM : need to check convergence order and incompressibility correction - Physical applications : add free movement of the sphere, add a solid solver Thank You for your attention ti 29

Development 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 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 information

Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)

Investigation 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 information

Implementation of a new discrete Immersed Boundary Method in OpenFOAM

Implementation of a new discrete Immersed Boundary Method in OpenFOAM Implementation of a new discrete Immersed Boundary Method in OpenFOAM CONSTANT Eddy*, COLIN-BELLOT Clothilde, FAVIER Julien, MELIGA Philippe, SERRE Éric Aix-Marseille Université, CNRS, Ecole Centrale Marseille,

More information

The Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids

The Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning

More information

Direct Numerical Simulation and Turbulence Modeling for Fluid- Structure Interaction in Aerodynamics

Direct Numerical Simulation and Turbulence Modeling for Fluid- Structure Interaction in Aerodynamics Available online at www.prace-ri.eu Partnership for Advanced Computing in Europe Direct Numerical Simulation and Turbulence Modeling for Fluid- Structure Interaction in Aerodynamics Thibaut Deloze a, Yannick

More information

Auto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial

Auto 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 information

Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić

Coupling 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 information

Debojyoti 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 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 information

Studies 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 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 information

A 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 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 information

Recent applications of overset mesh technology in SC/Tetra

Recent applications of overset mesh technology in SC/Tetra Recent applications of overset mesh technology in SC/Tetra NIA CFD Seminar October 6, 2014 Tomohiro Irie Software Cradle Co., Ltd. 1 Contents Introduction Software Cradle SC/Tetra Background of Demands

More information

EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES

EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary

More information

Overview and Recent Developments of Dynamic Mesh Capabilities

Overview 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 information

Computational Fluid Dynamics for Engineers

Computational Fluid Dynamics for Engineers Tuncer Cebeci Jian P. Shao Fassi Kafyeke Eric Laurendeau Computational Fluid Dynamics for Engineers From Panel to Navier-Stokes Methods with Computer Programs With 152 Figures, 19 Tables, 84 Problems and

More information

cuibm A GPU Accelerated Immersed Boundary Method

cuibm 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 information

Solving Partial Differential Equations on Overlapping Grids

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 information

Turbulence Modeling. Gilles Eggenspieler, Ph.D. Senior Product Manager

Turbulence Modeling. Gilles Eggenspieler, Ph.D. Senior Product Manager Turbulence Modeling Gilles Eggenspieler, Ph.D. Senior Product Manager 1 Overview The Role of Steady State (RANS) Turbulence Modeling Overview of Reynolds-Averaged Navier Stokes (RANS) Modeling Capabilities

More information

SPH: Why and what for?

SPH: 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 information

1.2 Numerical Solutions of Flow Problems

1.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 information

High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering

High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson

More information

Numerical Methods in Aerodynamics. Fluid Structure Interaction. Lecture 4: Fluid Structure Interaction

Numerical Methods in Aerodynamics. Fluid Structure Interaction. Lecture 4: Fluid Structure Interaction Fluid Structure Interaction Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark

More information

ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects

ALE 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 information

A higher-order finite volume method with collocated grid arrangement for incompressible flows

A higher-order finite volume method with collocated grid arrangement for incompressible flows Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.

More information

Best Practices Workshop: Overset Meshing

Best 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 information

EVALUATE SHOCK CAPTURING CAPABILITY WITH THE NUMERICAL METHODS IN OpenFOAM

EVALUATE SHOCK CAPTURING CAPABILITY WITH THE NUMERICAL METHODS IN OpenFOAM THERMAL SCIENCE: Year 2013, Vol. 17, No. 4, pp. 1255-1260 1255 Open forum EVALUATE SHOCK CAPTURING CAPABILITY WITH THE NUMERICAL METHODS IN OpenFOAM by Reza KHODADADI AZADBONI a*, Mohammad Rahim MALEKBALA

More information

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Reports of Research Institute for Applied Mechanics, Kyushu University No.150 (71 83) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 3: For the Case

More information

High-Fidelity Simulation of Unsteady Flow Problems using a 3rd Order Hybrid MUSCL/CD scheme. A. West & D. Caraeni

High-Fidelity Simulation of Unsteady Flow Problems using a 3rd Order Hybrid MUSCL/CD scheme. A. West & D. Caraeni High-Fidelity Simulation of Unsteady Flow Problems using a 3rd Order Hybrid MUSCL/CD scheme ECCOMAS, June 6 th -11 th 2016, Crete Island, Greece A. West & D. Caraeni Outline Industrial Motivation Numerical

More information

Explicit and Implicit Coupling Strategies for Overset Grids. Jörg Brunswig, Manuel Manzke, Thomas Rung

Explicit 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 information

Introduction To FLUENT. David H. Porter Minnesota Supercomputer Institute University of Minnesota

Introduction To FLUENT. David H. Porter Minnesota Supercomputer Institute University of Minnesota Introduction To FLUENT David H. Porter Minnesota Supercomputer Institute University of Minnesota Topics Covered in this Tutorial What you can do with FLUENT FLUENT is feature rich Summary of features and

More information

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Reports of Research Institute for Applied Mechanics, Kyushu University, No.150 (60-70) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 2: For the Case

More information

Partitioned strongly coupled Fluid-Structure Interaction

Partitioned 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 information

Math 690N - Final Report

Math 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 information

HPC Usage for Aerodynamic Flow Computation with Different Levels of Detail

HPC Usage for Aerodynamic Flow Computation with Different Levels of Detail DLR.de Folie 1 HPCN-Workshop 14./15. Mai 2018 HPC Usage for Aerodynamic Flow Computation with Different Levels of Detail Cornelia Grabe, Marco Burnazzi, Axel Probst, Silvia Probst DLR, Institute of Aerodynamics

More information

CFD in COMSOL Multiphysics

CFD 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 information

CFD VALIDATION FOR SURFACE COMBATANT 5415 STRAIGHT AHEAD AND STATIC DRIFT 20 DEGREE CONDITIONS USING STAR CCM+

CFD VALIDATION FOR SURFACE COMBATANT 5415 STRAIGHT AHEAD AND STATIC DRIFT 20 DEGREE CONDITIONS USING STAR CCM+ CFD VALIDATION FOR SURFACE COMBATANT 5415 STRAIGHT AHEAD AND STATIC DRIFT 20 DEGREE CONDITIONS USING STAR CCM+ by G. J. Grigoropoulos and I..S. Kefallinou 1. Introduction and setup 1. 1 Introduction The

More information

System Level Cooling, Fatigue, and Durability. Co-Simulation. Stuart A. Walker, Ph.D.

System 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 information

Overview of Traditional Surface Tracking Methods

Overview 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 information

LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.

LS-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 information

A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS

A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.

More information

DNV GL s 16th Technology Week

DNV 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 information

A fully implicit Navier-Stokes algorithm for unstructured grids incorporating a two-equation turbulence model

A fully implicit Navier-Stokes algorithm for unstructured grids incorporating a two-equation turbulence model Copyright 1996, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1996 A9618376, AIAA Paper 96-0414 A fully implicit Navier-Stokes algorithm for unstructured

More information

Solver Settings. Introductory FLUENT Training ANSYS, Inc. All rights reserved. ANSYS, Inc. Proprietary

Solver Settings. Introductory FLUENT Training ANSYS, Inc. All rights reserved. ANSYS, Inc. Proprietary Solver Settings Introductory FLUENT Training 2006 ANSYS, Inc. All rights reserved. 2006 ANSYS, Inc. All rights reserved. 5-2 Outline Using the Solver Setting Solver Parameters Convergence Definition Monitoring

More information

On the high order FV schemes for compressible flows

On the high order FV schemes for compressible flows Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech

More information

Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent

Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent MEGR 7090-003, Computational Fluid Dynamics :1 7 Spring 2015 Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent Rahul R Upadhyay Master of Science, Dept of Mechanical

More information

Turbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics

Turbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics Turbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics Rob J.M Bastiaans* Eindhoven University of Technology *Corresponding author: PO box 512, 5600 MB, Eindhoven, r.j.m.bastiaans@tue.nl

More information

Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak

Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society

More information

Numerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1

Numerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1 Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns

More information

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software

Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Reports of Research Institute for Applied Mechanics, Kyushu University No.150 (47 59) March 2016 Reproducibility of Complex Turbulent Using Commercially-Available CFD Software Report 1: For the Case of

More information

Numerical Analysis of Shock Tube Problem by using TVD and ACM Schemes

Numerical 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 information

The Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations

The 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 information

An 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 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 information

Two-dimensional laminar shock wave / boundary layer interaction

Two-dimensional laminar shock wave / boundary layer interaction Two-dimensional laminar shock wave / boundary layer interaction J.-Ch. Robinet (), V. Daru (,) and Ch. Tenaud () () SINUMEF Laboratory, ENSAM-PARIS () LIMSI-CNRS 5, Bd. de l Hôpital, PARIS 753, France

More information

Use of CFD in Design and Development of R404A Reciprocating Compressor

Use of CFD in Design and Development of R404A Reciprocating Compressor Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2006 Use of CFD in Design and Development of R404A Reciprocating Compressor Yogesh V. Birari

More information

Evaluation of Flow Solver Accuracy using Five Simple Unsteady Validation Cases

Evaluation of Flow Solver Accuracy using Five Simple Unsteady Validation Cases Evaluation of Flow Solver Accuracy using Five Simple Unsteady Validation Cases Bradford E. Green, Ryan Czerwiec Naval Air Systems Command, Patuxent River, MD, 20670 Chris Cureton, Chad Lillian, Sergey

More information

Flow in an Intake Manifold

Flow in an Intake Manifold Tutorial 2. Flow in an Intake Manifold Introduction The purpose of this tutorial is to model turbulent flow in a simple intake manifold geometry. An intake manifold is a system of passages which carry

More information

Multigrid Algorithms for Three-Dimensional RANS Calculations - The SUmb Solver

Multigrid Algorithms for Three-Dimensional RANS Calculations - The SUmb Solver Multigrid Algorithms for Three-Dimensional RANS Calculations - The SUmb Solver Juan J. Alonso Department of Aeronautics & Astronautics Stanford University CME342 Lecture 14 May 26, 2014 Outline Non-linear

More information

A THREE-DIMENSIONAL ADAPTIVE WAVELET METHOD FOR FLUID STRUCTURE INTERACTION

A 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 information

Comparison Between Different Immersed Boundary Conditions for Simulation of Complex Fluid Flows

Comparison Between Different Immersed Boundary Conditions for Simulation of Complex Fluid Flows Copyright 2011 Tech Science Press FDMP, vol.7, no.3, pp.241-258, 2011 Comparison Between Different Immersed Boundary Conditions for Simulation of Complex Fluid Flows A. Mark 1 2, R. Rundqvist 1 and F.

More information

FLOWVISION CFD FREQUENTLY ASKED QUESTIONS

FLOWVISION CFD FREQUENTLY ASKED QUESTIONS FLOWVISION CFD FREQUENTLY ASKED QUESTIONS 1. Installation and Licensing 1.1. Does FlowVision have floating licenses? 1.1.1. Actually all FlowVision licenses have floating capability and no extra fees are

More information

Introduction to C omputational F luid Dynamics. D. Murrin

Introduction 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 information

Possibility 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 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 information

ITU/FAA Faculty of Aeronautics and Astronautics

ITU/FAA Faculty of Aeronautics and Astronautics S. Banu YILMAZ, Mehmet SAHIN, M. Fevzi UNAL, Istanbul Technical University, 34469, Maslak/Istanbul, TURKEY 65th Annual Meeting of the APS Division of Fluid Dynamics November 18-20, 2012, San Diego, CA

More information

A Hybrid Cartesian Grid and Gridless Method for Compressible Flows

A Hybrid Cartesian Grid and Gridless Method for Compressible Flows rd AIAA Aerospace Sciences Meeting and Exhibit,, January 5, Reno, Nevada A Hybrid Cartesian Grid and Gridless Method for Compressible Flows Hong Luo and Joseph D. Baum Science Applications International

More information

CS205b/CME306. Lecture 9

CS205b/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 information

EFFECT OF ARTIFICIAL DIFFUSION SCHEMES ON MULTIGRID CONVERGENCE

EFFECT OF ARTIFICIAL DIFFUSION SCHEMES ON MULTIGRID CONVERGENCE AIAA Paper 95-1670 EFFECT OF ARTIFICIAL DIFFUSION SCHEMES ON MULTIGRID CONVERGENCE Seokkwan Yoon *, Antony Jameson t, and Dochan Kwak $ NASA Ames Research Center Moffett Field, California 94035 Abstract

More information

Algorithmic Developments in TAU

Algorithmic Developments in TAU Algorithmic Developments in TAU Ralf Heinrich, Richard Dwight, Markus Widhalm, and Axel Raichle DLR Institute of Aerodynamics and Flow Technology, Lilienthalplatz 7, 38108, Germany ralf.heinrich@dlr.de,

More information

FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD

FAST 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 information

Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders

Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.

More information

Adaptive Mesh Astrophysical Fluid Simulations on GPU. San Jose 10/2/2009 Peng Wang, NVIDIA

Adaptive Mesh Astrophysical Fluid Simulations on GPU. San Jose 10/2/2009 Peng Wang, NVIDIA Adaptive Mesh Astrophysical Fluid Simulations on GPU San Jose 10/2/2009 Peng Wang, NVIDIA Overview Astrophysical motivation & the Enzo code Finite volume method and adaptive mesh refinement (AMR) CUDA

More information

Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach

Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for

More information

A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws

A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)

More information

Usage of CFX for Aeronautical Simulations

Usage of CFX for Aeronautical Simulations Usage of CFX for Aeronautical Simulations Florian Menter Development Manager Scientific Coordination ANSYS Germany GmbH Overview Elements of CFD Technology for aeronautical simulations: Grid generation

More information

An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid

An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid J. V. Lassaline Ryerson University 35 Victoria St, Toronto, ON, M5B 2K3, Canada D. W. Zingg University

More information

An added mass partitioned algorithm for rigid bodies and incompressible flows

An added mass partitioned algorithm for rigid bodies and incompressible flows An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216 Collaborators Bill Henshaw,

More information

Metafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction

Metafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for

More information

2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123

2.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 information

Introducing Overdecomposition to Existing Applications: PlasComCM and AMPI

Introducing Overdecomposition to Existing Applications: PlasComCM and AMPI Introducing Overdecomposition to Existing Applications: PlasComCM and AMPI Sam White Parallel Programming Lab UIUC 1 Introduction How to enable Overdecomposition, Asynchrony, and Migratability in existing

More information

OverBlown : A Fluid Flow Solver For Overlapping Grids, User Guide, Version 1.0

OverBlown : A Fluid Flow Solver For Overlapping Grids, User Guide, Version 1.0 1 OverBlown : A Fluid Flow Solver For Overlapping Grids, User Guide, Version 1.0 William D. Henshaw Centre for Applied Scientific Computing Lawrence Livermore National Laboratory Livermore, CA, 94551.

More information

DES Turbulence Modeling for ICE Flow Simulation in OpenFOAM

DES Turbulence Modeling for ICE Flow Simulation in OpenFOAM 2 nd Two-day Meeting on ICE Simulations Using OpenFOAM DES Turbulence Modeling for ICE Flow Simulation in OpenFOAM V. K. Krastev 1, G. Bella 2 and G. Campitelli 1 University of Tuscia, DEIM School of Engineering

More information

The 3D DSC in Fluid Simulation

The 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 information

Solution of 2D Euler Equations and Application to Airfoil Design

Solution of 2D Euler Equations and Application to Airfoil Design WDS'6 Proceedings of Contributed Papers, Part I, 47 52, 26. ISBN 8-86732-84-3 MATFYZPRESS Solution of 2D Euler Equations and Application to Airfoil Design J. Šimák Charles University, Faculty of Mathematics

More information

Numerical 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 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 information

SIMULATION OF FLOW AROUND KCS-HULL

SIMULATION OF FLOW AROUND KCS-HULL SIMULATION OF FLOW AROUND KCS-HULL Sven Enger (CD-adapco, Germany) Milovan Perić (CD-adapco, Germany) Robinson Perić (University of Erlangen-Nürnberg, Germany) 1.SUMMARY The paper describes results of

More information

Modelling of Levitation Melting using a Fixed Mesh Method

Modelling of Levitation Melting using a Fixed Mesh Method International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, October 27-29, 2008 Modelling of Levitation Melting using a Fixed Mesh Method D. Hectors, E. Toorman, K. Van Reusel

More information

02_CRADLE 03_CSSRC 04_HSVA

02_CRADLE 03_CSSRC 04_HSVA A Computational Domain Domain Topology A1 1 steady domain 1 rotating domain Multiple domains Grid coupling technique A2 Sliding Overset Multiple ref. Frames B Propeller Representation Resolution B1 Geometrically

More information

CFD Simulation of a dry Scroll Vacuum Pump including Leakage Flows

CFD 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 information

RANS COMPUTATION OF RIBBED DUCT FLOW USING FLUENT AND COMPARING TO LES

RANS COMPUTATION OF RIBBED DUCT FLOW USING FLUENT AND COMPARING TO LES RANS COMPUTATION OF RIBBED DUCT FLOW USING FLUENT AND COMPARING TO LES Máté M., Lohász +*& / Ákos Csécs + + Department of Fluid Mechanics, Budapest University of Technology and Economics, Budapest * Von

More information

Problem description. The FCBI-C element is used in the fluid part of the model.

Problem 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 information

This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:

This 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 information

Numerical and theoretical analysis of shock waves interaction and reflection

Numerical and theoretical analysis of shock waves interaction and reflection Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for

More information

NUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.

NUMERICAL 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 information

Multigrid 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 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 information

Simulating 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 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 information

Investigation of the Influence of the Turbulent Transition on the Transonic Periodic Flow

Investigation of the Influence of the Turbulent Transition on the Transonic Periodic Flow Investigation of the Influence of the Turbulent Transition on the Transonic Periodic Flow C. Tulita *, E. Turkbeyler and S. Raghunathan Department of Aeronautical Engineering, The Quenn s University of

More information

Profile Catalogue for Airfoil Sections Based on 3D Computations

Profile Catalogue for Airfoil Sections Based on 3D Computations Risø-R-58(EN) Profile Catalogue for Airfoil Sections Based on 3D Computations Franck Bertagnolio, Niels N. Sørensen and Jeppe Johansen Risø National Laboratory Roskilde Denmark December 26 Author: Franck

More information

An 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 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 information

Second Symposium on Hybrid RANS-LES Methods, 17/18 June 2007

Second Symposium on Hybrid RANS-LES Methods, 17/18 June 2007 1 Zonal-Detached Eddy Simulation of Transonic Buffet on a Civil Aircraft Type Configuration V.BRUNET and S.DECK Applied Aerodynamics Department The Buffet Phenomenon Aircraft in transonic conditions Self-sustained

More information

TAU mesh deformation. Thomas Gerhold

TAU mesh deformation. Thomas Gerhold TAU mesh deformation Thomas Gerhold The parallel mesh deformation of the DLR TAU-Code Introduction Mesh deformation method & Parallelization Results & Applications Conclusion & Outlook Introduction CFD

More information

Keywords: flows past a cylinder; detached-eddy-simulations; Spalart-Allmaras model; flow visualizations

Keywords: flows past a cylinder; detached-eddy-simulations; Spalart-Allmaras model; flow visualizations A TURBOLENT FLOW PAST A CYLINDER *Vít HONZEJK, **Karel FRAŇA *Technical University of Liberec Studentská 2, 461 17, Liberec, Czech Republic Phone:+ 420 485 353434 Email: vit.honzejk@seznam.cz **Technical

More information

An 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 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 information