STCE. An (more) effective Discrete Adjoint Model for OpenFOAM
|
|
- Harry Norman
- 6 years ago
- Views:
Transcription
1 An (more) effective Discrete Adjoint Model for OpenFOAM Markus Towara, Uwe Naumann Software and Tools for Computational Engineering Science RWTH Aachen University EuroAD 2013, Oxford, 10. December 2013
2 A. Sen (aboutflow Project, discrete OpenFOAM) A. Daztouri (aboutflow Project, discrete ACE+) M. Foerster (fortissifoam Project, OpenFOAM in the cloud) me (discrete OpenFOAM) RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
3 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
4 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
5 OpenFOAM Open Field Operation and Manipulation Open-Source (GPLv3) CFD solver developed by OpenCFD Ltd., currently at version 2.2.x > 1 Mio. Lines of Code includes tools for meshing, pre-, post-processing rising adoption in industry and academia due to lack of licence costs well suited for parallel architectures RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
6 CFD Basics Navier Stokes equations for incompressible steady flow: v v = ν 2 v 1 ρ p v = 0 momentum conservation mass conservation Or in three Dimensions: vx vx vx vx + vy + vz = p x y z x vy vy vy vx + vy + vz = p x y z y vz vz vz vx + vy + vz = p x y z z u x + v + w = 0 y z + ν 2 vx x 2 + ν 2 vy x 2 + ν 2 vz x vx y vy y vz y vx z vy z vz z 2 RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
7 CFD Basics (Decoupled) Solution of the partial differential equations (SIMPLE-Algorithm): discretize / linearize momentum conservation equations solve momentum equations for velocity v, assume pressure p as known obtained velocity field fulfills momentum equation but not mass conservation equation as the pressure field was guessed and not correct discretize mass conservation equation use mass conservation equation to correct pressure field use new pressure field to correct velocity field Loop... RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
8 Topology Optimization no optimization added material reconstruction Add penalty term α to Navier-Stokes equation 1 : (v ) v = ν 2 v p αv 1 C. Othmer: A continuous adjoint formulation for the computation of topological and surface sensitivities of ducted flows. Intern. J. f. Num. Meth. in Fluids. p , RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
9 How to find appropriate α Define Cost Function J, e.g. total pressure loss between inlet and outlet: J = p v2 n dγ Γ Calculate sensitivity of the Cost function w.r.t. parameters α i J α i =??? Calculate updated porosity field α n+1, e.g.: α n+1 i = αi n λ J n, while insuring 0 < α i < α max Loop until α converged... α n i RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
10 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
11 Motivation for AD in Context of Shape Optimization AD promises: greater flexibility w.r.t new objectives easy adaption to new solver types / generations calculated derivatives are exact w.r.t the used discretization higher order derivatives available Problems: memory requirements how to retain parallelism? RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
12 How to implement n m, so Adjoint Mode is method of choice source code transformation operator overloading due to complexity of OpenFOAM code (templates, massively object oriented) the only suitable option operator overloading tool dco/c++ every floating point value in OpenFOAM is replaced with custom adjoint datatype dco::a1s::type needed intermediate values are stored in a structure called tape RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
13 How to implement put dco into src/openfoam include dco.hpp replace doubles with active datatype from dco OpenFOAM has own typedef for scalar floating point values just one substitution in theory we now just need to recompile OpenFOAM and are ready to go RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
14 OpenFOAM a1s mode in src/openfoam/primitives/scalar/doublescalar/doublescalar.h: replace: namespace Foam { } with: typedef double doublescalar;... #include dco.hpp namespace Foam { } typedef dco::a1s::type doublescalar;... RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
15 OpenFOAM t1s mode in src/openfoam/primitives/scalar/doublescalar/doublescalar.h: replace: namespace Foam { } with: typedef double doublescalar;... #include dco.hpp namespace Foam { } typedef dco::t1s::type doublescalar;... RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
16 in Reality some minor changes have to be made in the OpenFOAM code: unions don t support active datatypes some casts are missing and have to be done by hand some functions (pow,max,min) have to be adapted RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
17 usage of discrete OpenFOAM - t1s Black-Box tangent-linear Version of simplefoam, calculates J α i : double sens = 0; dco::t1s::set(alpha[i],1,1); for (runtime++;!runtime.end(); runtime++) {... // solve for U,p } // Sum pressure over inlet faces scaled with face area doublescalar J = gsum( p.boundaryfield() patch.magsf() ); dco::t1s::get(j,sens,1); Need to do this N-times to get full sensitivity field! RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
18 usage of discrete OpenFOAM - a1s Black-Box adjoint Version of simplefoam, calculates gradient of J: dco::a1s::static tape tape(tapesize); double sens = new double[alpha.size()]; for(int i=0; i<alpha.size(); i++) tape.register variable(alpha[i]); for (runtime++;!runtime.end(); runtime++) {... // solve for U,p } // Sum pressure over inlet faces scaled with face area doublescalar J = gsum( p.boundaryfield() patch.magsf() ); dco::a1s::set(j,1, 1); tape.interpret adjoint(); for(int i = 0; i<alpha.size(); i++) dco::a1s::get(alpha[i],sens[i], 1); RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
19 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
20 First Results Black-Box approach severely bound by memory bandwith! (Need to allocate around 20 GB) By doing Checkpointing we can actually get faster... RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
21 How to improve? Knowledge and Profiling reveals that most of calculation time and tape memory is spent in (iterative) linear solvers Analytical insight allows us to treat the adjoints of linear solvers analytically Lemma For a Linear Equation Sytem Ax = b we can calculate the adjoints for A and b by: 2 b = A T x A T b = x Ā = b x T This gives us an additional Linear Equation System which we have to solve during the gathering of the adjoints 2 M. Giles, Collected Matrix Derivative Results for Forward and Reverse Mode Algorithmic Differentiation, Advances in Automatic Differentiation 2008 RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
22 Treatment of Linear Solvers Black-Box Continuous Linear Solver when we encounter an linear solver during the augmented forward run we can stop taping when we encounter the gap in the tape during the interpretation we have to fill in the gap by hand have to remember A and x! RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
23 Application to OpenFOAM Abstract Base Class for all Linear Solvers exists can tackle all implemented solvers at once Discretization matrices are stored in a coordinate sparse matrix format transposing matrix straightforward take care not to mess up symmetry of matrices RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
24 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
25 Results we are seeing a nice improvement in both runtime and memory usage RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
26 Verification We calculate a reference solution with solvers GAMG for U and smoothsolver for p, both with tolerance 1e 10 After that we calculate solutions with different tolerances for the forward and adjoint linear solvers and compare to the reference solution 11k Cells in mesh Velocity Sensitivity RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
27 Verification Runtime Tape-Memory RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
28 Error RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
29 Runtime RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
30 Outline OpenFOAM Introduction Discrete OpenFOAM Improvements Verification and Results Work in Progress - PISO Foam RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
31 Work in Progress - Unsteady Adjoints Unsteady PISO solver with unsteady cost function (potentially depends on every timestep) RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
32 Work in Progress - Unsteady Adjoints RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
33 Work in Progress - Unsteady Adjoints RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
34 Summary and Outlook Fully discrete mode uses loads of memory Treatment of linear solvers greatly reduces this Its also faster Application to unsteady flow, heat transfer, turbulence in progress RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
35 Thank you! Questions? RWTH-Aachen, Discrete OpenFOAM, EuroAD / 35
Combining source transformation and operator overloading approaches for adjoint CFD solvers
Combining source transformation and operator overloading approaches for adjoint CFD solvers Azar Dastouri Sinan Gezgin Prof. Uwe Naumann Software and Tools for Computational Engineering, RWTH Aachen University,
More informationEfficiency of second-order differentiation schemes by algorithmic differentiation: Case Study
Efficiency of second-order differentiation schemes by algorithmic differentiation: Case Study Thorsten Lajewski Supervisor: Johannes Lotz STCE, RWTH Aachen May 6, 2013 1 Introduction In numerical calculations
More informationTopology optimization of heat conduction problems
Topology optimization of heat conduction problems Workshop on industrial design optimization for fluid flow 23 September 2010 Misha Marie Gregersen Anton Evgrafov Mads Peter Sørensen Technical University
More informationFlow and Heat Transfer in a Mixing Elbow
Flow and Heat Transfer in a Mixing Elbow Objectives The main objectives of the project are to learn (i) how to set up and perform flow simulations with heat transfer and mixing, (ii) post-processing and
More informationMulti-objective adjoint optimization of flow in duct and pipe networks
Multi-objective adjoint optimization of flow in duct and pipe networks Eugene de Villiers Thomas Schumacher 6th OPENFOAM Workshop PennState University, USA 13-16 June, 2011 info@engys.eu Tel: +44 (0)20
More informationAdjoint Solver Advances, Tailored to Automotive Applications
Adjoint Solver Advances, Tailored to Automotive Applications Stamatina Petropoulou s.petropoulou@iconcfd.com 1 Contents 1. Icon s Principal Work in FlowHead 2. Demonstration Cases 3. Icon s Further Development
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 informationAdjoint Solver Workshop
Adjoint Solver Workshop Why is an Adjoint Solver useful? Design and manufacture for better performance: e.g. airfoil, combustor, rotor blade, ducts, body shape, etc. by optimising a certain characteristic
More informationI. Introduction. Optimization Algorithm Components. Abstract for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany
An Aerodynamic Optimization Framework for the Automotive Industry, based on Continuous Adjoint and OpenFOAM E. Papoutsis-Kiachagias 1, V. Asouti 1, K. Giannakoglou 1, K. Gkagkas 2 1) National Technical
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 informationIntroduction to fluid mechanics simulation using the OpenFOAM technology
Introduction to fluid mechanics simulation using the OpenFOAM technology «Simulation in porous media from pore to large scale» Part II: Mesh complex geometries, application to the evaluation of permeability,
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 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 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 informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 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 informationShape optimisation using breakthrough technologies
Shape optimisation using breakthrough technologies Compiled by Mike Slack Ansys Technical Services 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Introduction Shape optimisation technologies
More informationDevelopment of a Consistent Discrete Adjoint Solver for the SU 2 Framework
Development of a Consistent Discrete Adjoint Solver for the SU 2 Framework Tim Albring, Max Sagebaum, Nicolas Gauger Chair for Scientific Computing TU Kaiserslautern 16th Euro-AD Workshop, Jena December
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 informationAerodynamic Study of a Realistic Car W. TOUGERON
Aerodynamic Study of a Realistic Car W. TOUGERON Tougeron CFD Engineer 2016 Abstract This document presents an aerodynamic CFD study of a realistic car geometry. The aim is to demonstrate the efficiency
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 informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
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 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 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 informationStudies of the ERCOFTAC Centrifugal Pump with OpenFOAM
Title 1/20 Studies of the ERCOFTAC Centrifugal Pump with OpenFOAM Olivier Petit Håkan Nilsson Outline Outline Geometry Boundary conditions Method and cases Comparison of numerical results with the available
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 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 informationAutomatic Hex-Dominant Mesh Generation for CFD Analysis of Formula One Car with cfmeshpro
Automatic Hex-Dominant Mesh Generation for CFD Analysis of Formula One Car with cfmeshpro Alen Cukrov and Franjo Juretić Creative Fields Ltd, X Vrbik 4, 10000 Zagreb, Croatia 1 Introduction This report
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 informationPerformance of Implicit Solver Strategies on GPUs
9. LS-DYNA Forum, Bamberg 2010 IT / Performance Performance of Implicit Solver Strategies on GPUs Prof. Dr. Uli Göhner DYNAmore GmbH Stuttgart, Germany Abstract: The increasing power of GPUs can be used
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 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 informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationThe Continuous Adjoint Method for the Computation of First- and Higher-Order Sensitivities
NATIONAL TECHNICAL UNIVERSITY OF ATHENS Parallel CFD & Optimization Unit Laboratory of Thermal Turbomachines The Continuous Adjoint Method for the Computation of First- and Higher-Order Sensitivities (Activities,
More informationSemi-automatic transition from simulation to one-shot optimization with equality constraints
Semi-automatic transition from simulation to one-shot optimization with equality constraints Lisa Kusch, Tim Albring, Andrea Walther, Nicolas Gauger Chair for Scientific Computing, TU Kaiserslautern, www.scicomp.uni-kl.de
More informationChart 1 Application of AD in Turbomachinery Design 19 th European Workshop on Automatic Differentiation Jan Backhaus DLR Cologne
www.dlr.de Chart 1 Application of AD in Turbomachinery Design 19 th European Workshop on Automatic Differentiation Jan Backhaus DLR Cologne www.dlr.de Chart 2 CFD based Optimization CRISP 1 rig Gradient-free
More informationA 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 informationTeam 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud
Team 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud CFD Support s OpenFOAM and UberCloud Containers enable efficient, effective, and easy access and use of MEET THE TEAM End-User/CFD
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 informationUse 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 informationComputational Fluid Dynamics using OpenCL a Practical Introduction
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 2011 http://mssanz.org.au/modsim2011 Computational Fluid Dynamics using OpenCL a Practical Introduction T Bednarz
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 informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationMultiphysics simulations of nuclear reactors and more
Multiphysics simulations of nuclear reactors and more Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg klasjareteg@chalmersse Division of Nuclear Engineering Department of Applied Physics Chalmers
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationPressure Correction Scheme for Incompressible Fluid Flow
AALTO UNIVERSITY School of Chemical Technology CHEM-E7160 Fluid Flow in Process Units Pressure Correction Scheme for Incompressible Fluid Flow Ong Chin Kai 620503 Lee De Ming Benedict 620448 Page 1 Abstract
More informationPublic Release: Native Overset Mesh in FOAM-Extend
Public Release: Native Overset Mesh in FOAM-Extend Vuko Vukčević and Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture, Uni Zagreb, Croatia Wikki Ltd. United Kingdom 5th OpenFOAM UK
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationFinite Volume Discretization on Irregular Voronoi Grids
Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling
More informationStreamlining Aircraft Icing Simulations. D. Snyder, M. Elmore
Streamlining Aircraft Icing Simulations D. Snyder, M. Elmore Industry Analysis Needs / Trends Fidelity Aircraft Ice Protection Systems-Level Modeling Optimization Background Ice accretion can critically
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 informationOpenFOAM: Open Platform for Complex Physics Simulations
OpenFOAM: Open Platform for Complex Physics Simulations Hrvoje Jasak h.jasak@wikki.co.uk, hrvoje.jasak@fsb.hr FSB, University of Zagreb, Croatia Wikki Ltd, United Kingdom 18th October 2007 OpenFOAM: Open
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 informationModeling and Simulation of Single Phase Fluid Flow and Heat Transfer in Packed Beds
Modeling and Simulation of Single Phase Fluid Flow and Heat Transfer in Packed Beds by:- Balaaji Mahadevan Shaurya Sachdev Subhanshu Pareek Amol Deshpande Birla Institute of Technology and Science, Pilani
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
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 informationOpenFOAM Programming the basic classes
OpenFOAM Programming the basic classes Prof Gavin Tabor Friday 25th May 2018 Prof Gavin Tabor OpenFOAM Programming the basic classes Friday 25th May 2018 1 / 30 OpenFOAM : Overview Overview : programming
More informationLab 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 informationOptimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston Optimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells James A. Peitzmeier *1, Steven Kapturowski 2 and Xia Wang
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
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 informationInverse Design of an Indoor Environment by Using a CFD-based Adjoint Method with Adaptive Step Size for Adjusting Design Parameters
Zhao, X., Liu, W., Liu, S., Zou, Y., and Chen, Q. 2017. Inverse design of an indoor environment by using a CFD-based adjoint method with adaptive step size for adjusting design parameter Numerical Heat
More informationComputational 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 informationRBF Morph An Add-on Module for Mesh Morphing in ANSYS Fluent
RBF Morph An Add-on Module for Mesh Morphing in ANSYS Fluent Gilles Eggenspieler Senior Product Manager 1 Morphing & Smoothing A mesh morpher is a tool capable of performing mesh modifications in order
More informationThe Adjoint Method Hits the Road: Applications in Car Aerodynamics
Symposium in Honor of Antony Jameson s 80th Birthday, Stanford University 2014 The Adjoint Method Hits the Road: Applications in Car Aerodynamics Dr. Carsten Othmer, Volkswagen AG, Corporate Research,
More informationStudy Of Overloading Effects, In A Refrigerated Display Case
Study Of Overloading Effects, In A Refrigerated Display Case Sandeep Palaksha Senior CAE Engineer HUSSMANN Basavanagudi Bangalore 04, India p.sandeep@hussmann.com Narasimhamurthy CAE Engineer HUSSMANN
More informationUsing Multiple Rotating Reference Frames
Tutorial 10. Using Multiple Rotating Reference Frames Introduction Many engineering problems involve rotating flow domains. One example is the centrifugal blower unit that is typically used in automotive
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 informationCo-Simulation von Flownex und ANSYS CFX am Beispiel einer Verdrängermaschine
Co-Simulation von Flownex und ANSYS CFX am Beispiel einer Verdrängermaschine Benoit Bosc-Bierne, Dr. Andreas Spille-Kohoff, Farai Hetze CFX Berlin Software GmbH, Berlin Contents Positive displacement compressors
More informationMesh Morphing and the Adjoint Solver in ANSYS R14.0. Simon Pereira Laz Foley
Mesh Morphing and the Adjoint Solver in ANSYS R14.0 Simon Pereira Laz Foley 1 Agenda Fluent Morphing-Optimization Feature RBF Morph with ANSYS DesignXplorer Adjoint Solver What does an adjoint solver do,
More informationDirect Numerical Simulation of a Low Pressure Turbine Cascade. Christoph Müller
Low Pressure NOFUN 2015, Braunschweig, Overview PostProcessing Experimental test facility Grid generation Inflow turbulence Conclusion and slide 2 / 16 Project Scale resolving Simulations give insight
More informationChapter 1 Introduction
Chapter 1 Introduction Ibrahim Sezai Department of Mechanical Engineering Eastern Mediterranean University Fall 2009-2010 What is CFD? CFD is the simulation of fluids engineering systems using modeling
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 informationUsing Multiple Rotating Reference Frames
Tutorial 9. Using Multiple Rotating Reference Frames Introduction Many engineering problems involve rotating flow domains. One example is the centrifugal blower unit that is typically used in automotive
More informationReal-Time Shape Editing using Radial Basis Functions
Real-Time Shape Editing using Radial Basis Functions, Leif Kobbelt RWTH Aachen Boundary Constraint Modeling Prescribe irregular constraints Vertex positions Constrained energy minimization Optimal fairness
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 informationAerodynamic Topology Optimization: Some Observations on Hysteresis in Separated Flows
Aerodynamic Topology Optimization: Some Observations on Hysteresis in Separated Flows Carsten Othmer Electronics Research Laboratory, Volkswagen Group of America Inc., Belmont, CA 94002, USA David E. Manosalvas-Kjono
More informationHigh Performance Computing
High Performance Computing ADVANCED SCIENTIFIC COMPUTING Dr. Ing. Morris Riedel Adjunct Associated Professor School of Engineering and Natural Sciences, University of Iceland Research Group Leader, Juelich
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 informationFlow 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 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 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 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 informationOptimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers
Optimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers Jeffrey A. Fike and Juan J. Alonso Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305,
More informationIsotropic Porous Media Tutorial
STAR-CCM+ User Guide 3927 Isotropic Porous Media Tutorial This tutorial models flow through the catalyst geometry described in the introductory section. In the porous region, the theoretical pressure drop
More informationComputational Fluid Dynamics Simulation of a Rim Driven Thruster
Computational Fluid Dynamics Simulation of a Rim Driven Thruster Aleksander J Dubas, N. W. Bressloff, H. Fangohr, S. M. Sharkh (University of Southampton) Abstract An electric rim driven thruster is a
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 informationSolving systems of nonlinear equations with arbitrary precision
Solving systems of nonlinear equations with arbitrary precision CES Seminar January 2017 Presented by Supervisor Gregor Heiming Dr. Ralf Hannemann-Tamás Systemverfahrenstechnik AVT, RWTH Aachen University
More informationLES of Vortex Formation in a Spherical Dimple and Genetic Optimization of Heat Transfer Surfaces
LES of Vortex Formation in a Spherical Dimple and Genetic Optimization of Heat Transfer Surfaces J.Turnow, N. Kornev, E. Hassel Department of Technical Thermodynamics University of Rostock Germany 5 th
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 informationUse of numerical flow simulations (CFD) for optimising heat exchangers
www.guentner.eu Technical article from 03.04.2017 Author Dr. Andreas Zürner Research Güntner GmbH & Co. KG Use of numerical flow simulations (CFD) for optimising heat exchangers Numerical flow simulations
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 informationmidas NFX 2017R1 Release Note
Total Solution for True Analysis-driven Design midas NFX 2017R1 Release Note 1 midas NFX R E L E A S E N O T E 2 0 1 7 R 1 Major Improvements Midas NFX is an integrated finite element analysis program
More informationVerification of Laminar and Validation of Turbulent Pipe Flows
1 Verification of Laminar and Validation of Turbulent Pipe Flows 1. Purpose ME:5160 Intermediate Mechanics of Fluids CFD LAB 1 (ANSYS 18.1; Last Updated: Aug. 1, 2017) By Timur Dogan, Michael Conger, Dong-Hwan
More informationAM119: Yet another OpenFoam tutorial
AM119: Yet another OpenFoam tutorial Prof. Trask April 11, 2016 1 Todays project Today we re going to implement a projection method for the Navier-Stokes, learn how to build a mesh, and explore the difference
More informationPermeable and Absorbent Materials in Fluid Simulations
Permeable and Absorbent Materials in Fluid Simulations Nate Andrysco Bedrich Benes Department of Computer Science Slide 1 Motivation Fluid simulations are geared toward impermeable materials What if you
More informationFree Convection Cookbook for StarCCM+
ME 448/548 February 28, 2012 Free Convection Cookbook for StarCCM+ Gerald Recktenwald gerry@me.pdx.edu 1 Overview Figure 1 depicts a two-dimensional fluid domain bounded by a cylinder of diameter D. Inside
More informationParallel Computing and AD. Tools and Issues
Parallel Computing and AD Tools and Issues Graz, 27.06.2013 Matthias Meuter Parallel Computing 2 Parallel Computing Multicore Desktop PCs Large Clusters 3 Parallel Computing Multicore Desktop PCs Large
More information