Description of AdjointShapeOptimizationFoam and how to implement new cost functions
|
|
- Simon McCarthy
- 5 years ago
- Views:
Transcription
1 Description of AdjointShapeOptimizationFoam and how to implement new cost functions A project work in the course CFD with OpenSource software, taught by Håkan Nilsson. Chalmers University of Technology, Gothenburg, Sweden December 10, 2013 December 10, /
2 Overview Introduction Theory Implementation Modifications Results December 10, /
3 Introduction Optimization in CFD Increased importance Optimization algorithm - Simplex - Evolution Strategies - Gradient based -... Demanding - Number of variables - Constraints December 10, /
4 Introduction Adjoint approach Gradient based method Tool to compute sensitivity w.r.t. design variables Well behaved objective function Large number of design variables Constraints, number of functions December 10, /
5 Optimization problem Objective function minimize J = J(v, p, α) (1) such that (v )v + p (2νD(v)) + αv = 0, (2) v = 0. (3) Gerneral cost function Incompressible, steady state RANS equation Porosity α Topology optimization December 10, /
6 Optimization problem Lagrangian relaxation minimize L := J + Lagrange multipliers, (u, q) A penalty to violate the constraints Adjoint velocity, u, and pressure, q No physical meaning Ω (u, q)rdω. (4) December 10, /
7 Optimization problem Variations w.r.t. flow and design variables δl = δ α L + δ v L + δ p L, (5) δ v L + δ p L = 0. (6) Total variation. Choose adjoint variables so that the variation with respect to flow variables vanishes. Equation of the adjoint variables and an expression for the sensitivity of L with respect to the porosity of each cell, derivations out of the scope of the project. Approximations needed to arrive at the adjoint equations. - Forzen turbulence. - Ducted flow specialization. Integration by parts to arrive at volume contribution and boundary contribution. December 10, /
8 Resulting equations Sensitivity Expression for cell i V i the cell volume. L α i = u i v i V i, (7) Equation of the adjoint variables and an expression for the sensitivity of L with respect to the porosity of each cell December 10, /
9 Resulting equations Adjoint equations 2D(u)v = q + (2νD(u)) αu J Ω v, (8) u = J Ω p. (9) Contribution from volumetric part of the cost function, J Ω. December 10, /
10 Resulting equations Adjoint boundary conditions Wall and inlet boundary conditions Outlet boundary conditions u t = 0, (10) u n = J Ω p, (11) n q = 0. (12) q = u v + u n v n + ν(n )u n + J Γ v n, (13) 0 = v n u t + ν(n )u t + J Γ v t. (14) December 10, /
11 Resulting equations Adjoint boundary conditions Wall Inlet Outlet v No-slip Prescribed value Zero gradient p Zero gradient Zero gradient p = 0 Table: Boundary conditions for the primal quantities v and p. Only valid for specific primal boundary condition. Contribution from both the volume and boundary part of the cost function, J Ω and J Γ respectively. Eq. 13 and 14 used to calculate the adjoint pressure and tangential part of the adjoint velocity at the outlet. December 10, /
12 Resulting equations Summary of the theory Introducing a Lagrange relaxation, with multipliers u and q. Choice of adjoint variables gives adjoint equations and an expression for the sensititivity. Lengthy derivations. Approximations and requirements to ensure a valid method. If the cost function satisfies J ω = 0 the dependence of J enters the solver only through the boundary conditions. The optimization problem reduces to solving two flow equations, the primal and the adjoint. When the sensitivity is calculated an algorithm is needed to update the porosity. December 10, /
13 Steepest descent algorithm Steepest descent p k = f(x k ), (15) α n+1 = α n u i v i V i δ. (16) General search direction in the opposite direction of the gradient of the objective function. Considering a linear problem a minimum will be found for small enough steps in the search direction. Using the sensitivity, L/ α to update. The step length δ. An optimization problem of its own. Under relaxation and limits in the final implementation. December 10, /
14 Solver Solution procedure Main-loop Update the porosity Steepest descent Solve the primal system (v, p) Solve the adjoint system (u, q) Frozen tubulence Convergence assessment End time or convergence condition Converged End Figure 1: Block scheme of the solution procedure used in. December 10, /
15 Solver $FOAM SOLVERS/incompressible/ /.C Let us study the implementation together, trying to identify the underlying theory we have covered so far. In the order it is implemented in the code. Porosity update - Implementation of the steepest descent using under relaxation. - Correct? The primal pressure-velocity SIMPLE corrector - How is the sink term implemented? The adjoint pressure-velocity SIMPLE corrector - Is the current implementation done for objective functions with volumteric contribution (J Ω 0)? December 10, /
16 Solver Porosity update 94 a l p h a += 95 mesh. f i e l d R e l a x a t i o n F a c t o r ( "alpha" ) 96 ( min ( max( alpha + lambda (Ua & U), zeroalpha ), alphamax ) alpha ) ; Listing 1: file.c α n+1 = α n (1 γ) + γmin (max ((α n + u i v i V i δ).0), α max ), (17) December 10, /
17 Solver Porosity update However the there seems to be an incorrect sign in the steepest descent implementation. The correct eq. and implementation should be α n+1 = α n (1 γ) + γmin (max ((α n u i v i V i δ), 0), α max ). (18) 94 a l p h a += 95 mesh. f i e l d R e l a x a t i o n F a c t o r ( "alpha" ) 96 ( min ( max( alpha lambda (Ua & U), zeroalpha ), alphamax ) alpha ) ; Listing 2: file.c December 10, /
18 Solver Primal SIMPLE loop 103 // Momentum p r e d i c t o r tmp<f v V e c t o r M a t r i x > UEqn 106 ( 107 fvm : : d i v ( phi, U) t u r b u l e n c e >d i v D e v R e f f (U) fvm : : Sp ( alpha, U) 110 ) ; UEqn ( ). r e l a x ( ) ; s o l v e (UEqn ( ) == f v c : : grad ( p ) ) ; Listing 3: file.c December 10, /
19 Solver Primal SIMPLE loop fvm::sp(alpha, U) Implicit implementation of the extra source term. December 10, /
20 Solver Adjoint SIMPLE loop 156 // A d j o i n t Momentum p r e d i c t o r volvectorfield adjointtransposeconvection ( ( fvc : : grad (Ua) & U ) ) ; 159 // v o l V e c t o r F i e l d a d j o i n t T r a n s p o s e C o n v e c t i o n 160 // ( 161 // fvc : : reconstruct 162 // ( 163 // mesh. magsf ( ) ( f v c : : sngrad (Ua) & f v c : : i n t e r p o l a t e (U) ) 164 // ) 165 // ) ; z e r o C e l l s ( a d j o i n t T r a n s p o s e C o n v e c t i o n, i n l e t C e l l s ) ; tmp<f v V e c t o r M a t r i x > UaEqn 170 ( 171 fvm : : d i v ( phi, Ua) 172 a d j o i n t T r a n s p o s e C o n v e c t i o n t u r b u l e n c e >d i v D e v R e f f (Ua) fvm : : Sp ( alpha, Ua) 175 ) ; UaEqn ( ). r e l a x ( ) ; s o l v e ( UaEqn ( ) == f v c : : grad ( pa ) ) ; Listing 4: file.c December 10, /
21 Solver Adjoint SIMPLE loop Similar to the primal system. No additional terms containing information about the cost function, i.e. current implementation treats cost functions of the type J Γ = 0 December 10, /
22 Boundary conditions of the adjoint variables Adjoint boundary conditions The bad results when the sign in the steepest descent algorithm is changed indicates additional problems with the current solver. Since the actual solver seems to agree with the theory, what remains to be examined is the boundary conditions. First the boundary condition for a specific cost function need to be calculated. December 10, /
23 Boundary conditions of the adjoint variables Power dissipation J := Only boundary contribution. Γ (p v2 )v n dγ. (19) Energy loss, due to pressure drop and kinetic energy. December 10, /
24 Boundary conditions of the adjoint variables Inlet boundary conditions u t = 0, (20) u n = J Ω p, (21) n q = 0. (22) u t = { 0, u n = n q = 0. 0 at wall, v n at inlet, (23) December 10, /
25 Boundary conditions of the adjoint variables Inlet boundary conditions Use existing conditions - Prescribed value - No slip - Zero-gradient Could be useful to implement a inlet boundary condition for the adjoint velocity, if the velocity profile of the primal velocity at the inlet is not easily described. December 10, /
26 Boundary conditions of the adjoint variables Outlet boundary conditions q = u v + u n v n + ν(n )u n + J Γ v n, (24) 0 = v n u t + ν(n )u t + J Γ v t. (25) { q = u v + u n v n + ν(n )u n 1 2 v2 v 2 n, 0 = v n (u t v t ) + ν(n )u t. (26) December 10, /
27 Boundary conditions of the adjoint variables Outlet boundary conditions Use eq. 26 to prescribe a value to the adjoint pressure Use eq. 26 to prescribe a value to the tangential component of the adjoint velocity u t = ν u neigh,t v p,n v p,t v p,n + ν, (27) December 10, /
28 Boundary conditions of the adjoint variables Implementation of pressure condition 86 // Member F u n c t i o n s // void Foam : : adjointoutletpressurefvpatchscalarfield : : updatecoeffs ( ) 89 { 90 i f ( updated ( ) ) 91 { 92 return ; 93 } const fvspatchfield <scalar>& phip = 96 patch ( ). lookuppatchfield <s u r f a c e S c a l a r F i e l d, s c a l a r >("phi" ) ; const fvspatchfield <scalar>& phiap = 99 patch ( ). lookuppatchfield <s u r f a c e S c a l a r F i e l d, s c a l a r >("phia" ) ; const f v P a t c h F i e l d <v e c t o r >& Up = 102 patch ( ). lookuppatchfield <v o l V e c t o r F i e l d, vector >("U" ) ; const f v P a t c h F i e l d <v e c t o r >& Uap = 105 patch ( ). lookuppatchfield <v o l V e c t o r F i e l d, vector >("Ua" ) ; operator==((phiap /patch ( ). magsf ( ) 1.0) phip /patch ( ). magsf ( ) + (Up & Uap ) ) ; fixedvaluefvpatchscalarfield : : updatecoeffs ( ) ; 110 } Listing 5: file adjointoutletpressurefvpatchscalarfield.c December 10, /
29 Boundary conditions of the adjoint variables Implementation of velocity condition 83 // Member F u n c t i o n s // // Update the coefficients associated with the patch f i e l d 86 void Foam : : adjointoutletvelocityfvpatchvectorfield : : updatecoeffs ( ) 87 { 88 i f ( updated ( ) ) 89 { 90 return ; 91 } const fvspatchfield <scalar>& phiap = 94 patch ( ). lookuppatchfield <s u r f a c e S c a l a r F i e l d, s c a l a r >("phia" ) ; const f v P a t c h F i e l d <v e c t o r >& Up = 97 patch ( ). lookuppatchfield <v o l V e c t o r F i e l d, vector >("U" ) ; s c a l a r F i e l d Un(mag( patch ( ). n f ( ) & Up ) ) ; 100 vectorfield UtHat ( (Up patch ( ). nf ( ) Un ) / (Un + SMALL ) ) ; v e c t o r F i e l d Uan ( patch ( ). nf ( ) ( patch ( ). n f ( ) & p a t c h I n t e r n a l F i e l d ( ) ) ) ; vectorfield : : operator=(phiap patch ( ). Sf ( ) / sqr ( patch ( ). magsf ( ) ) + UtHat ) ; 105 // vectorfield : : operator=(uan + UtHat ) ; fixedvaluefvpatchvectorfield : : updatecoeffs ( ) ; 108 } Listing 6: file adjointoutletvelocityfvpatchvectorfield.c December 10, /
30 Boundary conditions of the adjoint variables Analyzing the implementation Listing 5 gives an implementation of the pressure condition according to compare to q = (u n 1)v n + u v. (28) q = u v + u n v n + ν(n )u n 1 2 v2 v 2 n, (29) Listing 6 gives an implementation of the velocity condition according to compare to u p,t = v p v p,n u p,n + SMALL. (30) u t = ν u neigh,t v p,n v p,t v p,n + ν, (31) December 10, /
31 Boundary conditions of the adjoint variables Analyzing the theory Not equal to the theory. Another cost function, total pressure loss. Approximations or derivations done makes the cost function hard to identify. December 10, /
32 Analyzing the theory Changing the sign of the steepest descent implementation. Implementing the cost functions according to derived equations. December 10, /
33 Adjoint pressure boundary condition Similar to the existing implementation A few additional terms - Velocity in the neighbouring node - Effective viscosity const scalarfield& deltainv = patch ( ). deltacoeffs ( ) ; // distance ˆ( 1) between patch and neighbouring node. // create the object needed to get the v i s c o u s i t y : const incompressible : : RASModel& rasmodel = db ( ). lookupobject<incompressible : : RASModel>("RASProperties" ) ; // n u e f f from the t u r u l e n c e model, rasmodel : s c a l a r F i e l d n u e f f = rasmodel. n u E f f ( ) ( ). b o u n d a r y F i e l d ( ) [ patch ( ). i n d e x ( ) ] ; // Neighbouring node s v e l o c i t y ( normal component ) : s c a l a r F i e l d Uaneigh n = (Uap. p a t c h I n t e r n a l F i e l d ( ) & patch ( ). nf ( ) ) ; Listing 7: file myadjointoutletvelocityfvpatchvectorfield.c December 10, /
34 Adjoint velocity boundary condition Vector fields instead of scalar fields // patchinternalfield to get the adjoint v e l o c i t y of neighbouring node. vectorfield Uaneigh = Uap. patchinternalfield ( ) ; // I t s tangential and normal components vectorfield Uaneigh n = ( Uaneigh & normal ) normal ; v e c t o r F i e l d Uaneigh t = Uaneigh Uaneigh n ; Listing 8: file myadjointoutletvelocityfvpatchvectorfield.c December 10, /
35 Case settings The boundary conditions of the primal and adjoint variables are set according to the table below. More detailed information and m4 script of the box example case can be found in the provided files. Wall Inlet Outlet u Prescribed value, 0 Prescribed value, u n = v n adjointoutletvelocitypower q Zero gradient Zero gradient adjointoutletpressurepower December 10, /
36 Power dissipation pitzdaily porosity field December 10, /
37 Power dissipation pitzdaily value of the objective function valueobjfunc.xy J Iteration step December 10, /
38 Power dissipation pitzdaily value of the objective function, without porosity update valueobjfunc.xy J Iteration step December 10, /
39 Power dissipation Pipe bend example December 10, /
40 Power dissipation Thank you for listening! Questions December 10, 2013 /
Rhie-Chow interpolation in OpenFOAM 1
Rhie-Chow interpolation in OpenFOAM 1 Fabian Peng Kärrholm Department of Applied Mechanics Chalmers Univesity of Technology Göteborg, Sweden, 2006 1 Appendix from Numerical Modelling of Diesel Spray Injection
More informationINDUSTRIAL APPLICATION OF CONTINUOUS ADJOINT FLOW SOLVERS FOR THE OPTIMIZATION OF AUTOMOTIVE EXHAUST SYSTEMS. C. Hinterberger, M.
CFD & OPTIMIZATION 2011-069 An ECCOMAS Thematic Conference 23-25 May 2011, Antalya TURKEY INDUSTRIAL APPLICATION OF CONTINUOUS ADJOINT FLOW SOLVERS FOR THE OPTIMIZATION OF AUTOMOTIVE EXHAUST SYSTEMS C.
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationRANS turbulence treatment for continuous adjoint optimization
Turbulence, Heat and Mass Transfer 8 c 2015 Begell House, Inc. RANS turbulence treatment for continuous adjoint optimization M. Popovac 1, C. Reichl 1, H. Jasa 2,3 and H. Rusche 3 1 Austrian Institute
More informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationStratified scavenging in two-stroke engines using OpenFOAM
Stratified scavenging in two-stroke engines using OpenFOAM Håkan Nilsson, Chalmers / Applied Mechanics / Fluid Dynamics 1 Acknowledgements I would like to thank Associate Professor Håkan Nilsson at the
More informationLarge Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes
Large Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes Thomas Gallinger * Alexander Kupzok Roland Wüchner Kai-Uwe Bletzinger Lehrstuhl für Statik Technische Universität München
More informationLecture 2: Fundamentals 2/3: Navier-Stokes Equation + Basics of Discretization
Lecture 2: Fundamentals 2/3: Navier-Stokes Equation Basics of Discretization V.Vuorinen Aalto University School of Engineering CFD Course, Spring 2018 January 15th 2018, Otaniemi ville.vuorinen@aalto.fi
More informationBoundary Conditions - Inlet
Boundary Conditions - Inlet Boundary Conditions Need? What is physical and mathematical significance of a boundary condition? The boundary conditions of any problem is used to define the upper and lower
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationHow to implement your own turbulence model(1/3)
How to implement your own turbulence model(1/3) The implementations of the turbulence models are located in $FOAM_SRC/turbulenceModels Copythesourceoftheturbulencemodelthatismostsimilartowhatyouwanttodo.Inthis
More informationMicro-Scale CFD Modeling of Packed-Beds
Micro-Scale CFD Modeling of Packed-Beds Daniel P. Combest and Dr. P.A. Ramachandran and Dr. M.P. Dudukovic Chemical Reaction Engineering Laboratory (CREL) Department of Energy, Environmental, and Chemical
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationHow to implement your own turbulence model (1/3)
How to implement your own turbulence model 1/3) The implementations of the turbulence models are located in $FOAM_SRC/turbulenceModels Copy the source of the turbulence model that is most similar to what
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationstudies Maryse Page Hydro-Québec, Research Institute Håkan Nilsson Chalmers University of Technology Omar Bounous Chalmers University of Technology
OpenFOAM Turbomachinery Working Group: ERCOFTAC conical diffuser case-studies studies Maryse Page Hydro-Québec, Research Institute Olivier Petit Chalmers University of Technology Håkan Nilsson Chalmers
More informationCoupled calculations in OpenFOAM -
Coupled calculations in OpenFOAM - Multiphysics handling, structures and solvers, Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg Chalmers University of Technology November 14, 2012 Outline
More informationCFDOFAIRFLOWINHYDROPOWERGENERATORS FOR CONVECTIVE COOLING, USING OPENFOAM
CFDOFAIRFLOWINHYDROPOWERGENERATORS FOR CONVECTIVE COOLING, USING OPENFOAM ECCOMAS CFD 21 Pirooz Moradnia, Håkan Nilsson Lisbon-Portugal 21-6-16 Pirooz Moradnia, Chalmers/ Applied Mechanics/ Fluid Dynamics
More informationviscousheatingsolver Martin Becker
viscousheatingsolver Adding temperature transport and viscous heating to simplefoam Solver and Test Case Martin Becker martin_becker@dhcae-tools.de 11/09/2012 Abstract The modifications of the simplefoam
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationCHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM
CHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM Suheni and Syamsuri Department of Mechanical Engineering, Adhi Tama Institute of Technology Surabaya, Indonesia E-Mail:
More informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
More informationUnsteady Rotor-Stator Simulation of the U9 Kaplan Turbine Model
Unsteady Rotor-Stator Simulation of the U9 Kaplan Turbine Model Olivier Petit, Håkan Nilsson Chalmers University 6 th OpenFOAM workshop Penn State University, USA 13-16 June 2011 Develop OpenFOAM as the
More informationKeywords - Gas Turbine, Exhaust Diffuser, Annular Diffuser, CFD, Numerical Simulations.
Numerical Investigations of PGT10 Gas Turbine Exhaust Diffuser Using Hexahedral Dominant Grid Vaddin Chetan, D V Satish, Dr. Prakash S Kulkarni Department of Mechanical Engineering, VVCE, Mysore, Department
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationA Thorough Description Of How Wall Functions Are Implemented In OpenFOAM
Cite as: Fangqing Liu.: A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM. In Proceedings of CFD with OpenSource Software, 2016, Edited by Nilsson. H., http://www.tfd.chalmers.se/~hani/kurser/os_cfd_2016
More informationPredictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical and CFD Method
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869, Volume-3, Issue-5, May 2015 Predictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationSimplified flow around a propeller for the course CFD with OpenSource Software
Simplified flow around a propeller for the course CFD with OpenSource Software Gonzalo Montero Applied Mechanics/Fluid Dynamics, Chalmers University of Technology, Gothenburg, Sweden 2015-12-08 Gonzalo
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationMAE 598 Project #1 Jeremiah Dwight
MAE 598 Project #1 Jeremiah Dwight OVERVIEW A simple hot water tank, illustrated in Figures 1 through 3 below, consists of a main cylindrical tank and two small side pipes for the inlet and outlet. All
More informationConservation of Angular Momentum
10 March 2017 Conservation of ngular Momentum Lecture 23 In the last class, we discussed about the conservation of angular momentum principle. Using RTT, the angular momentum principle was given as DHo
More informationDraft Tube calculations using OpenFOAM-1.5dev and validation with FLINDT data
6th OpenFOAM Workshop, June 13-16 2011, PennState University, USA Draft Tube calculations using OpenFOAM-1.5dev and validation with FLINDT data C. Devals, Y. Zhang and F.Guibault École Polytechnique de
More informationDevelopment of adjoint-based optimization methods for ducted flows in vehicles
THESIS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY IN THERMO AND FLUID DYNAMICS Development of adjoint-based optimization methods for ducted flows in vehicles EYSTEINN HELGASON Department of Applied Mechanics
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationTOPOLOGY OPTIMIZATION OF NAVIER-STOKES FLOW IN MICROFLUIDICS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationLevel Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationGENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION
25 th National Symposium on Wind Engineering, Tokyo, Japan, 3-5 December 2018 第 25 回風工学シンポジウム (2018) GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION Takafumi
More informationINTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA)
INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) Title 3 Column (full page) 2 Column What is Finite Element Analysis? 1 Column Half page The Finite Element Method The Finite Element Method (FEM)
More informationFluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers
Fluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers R. Sanchez 1, D. Thomas 2, R. Palacios 1, V. Terrapon 2 1 Department of Aeronautics, Imperial College London 2 Department
More informationDEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.
DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M. Sahu 1, Kishanjit Kumar Khatua and Kanhu Charan Patra 3, T. Naik 4 1, &3 Department of Civil Engineering, National Institute of technology,
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationCFD with OpenSource software A course at Chalmers University of Technology Taught by HÅKAN NILSSON. Project work: intersettlingfoam
Taught by HÅKAN NILSSON Project work: intersettlingfoa Developed for OpenFOAM-2.2.0 Author: Pedra Rain Deceber, 2013 2/21/2014 1 Presentation Outline: The objective of this presentation Introduction to
More informationDevelopment of a consistent and conservative Eulerian - Eulerian algorithm for multiphase flows
Development of a consistent and conservative Eulerian - Eulerian algorithm for multiphase flows Ana Cubero Alberto Sanchez-Insa Norberto Fueyo Numerical Fluid Dynamics Group University of Zaragoza Spain
More informationA numerical investigation of tip clearance flow in Kaplan water turbines
Published in the proceedings of HYDROPOWER INTO THE NEXT CENTURY - III, 1999. ISBN 9522642 9 A numerical investigation of tip clearance flow in Kaplan water turbines M.Sc. H. Nilsson Chalmers University
More informationOpenFOAM SIMULATION OF THE FLOW IN THE HÖLLEFORSEN DRAFT TUBE MODEL
Turbine-99 III Proceedings of the third IAHR/ERCOFTAC workshop on draft tube flow 8-9 December 2005, Porjus, Sweden Paper No. 8 OpenFOAM SIMULATION OF THE FLOW IN THE HÖLLEFORSEN DRAFT TUBE MODEL Nilsson
More informationTransient Thermal Flow and Thermal Stress Analysis Coupled NASTRAN and SC/Tetra
Transient Thermal Flow and Thermal Stress Analysis Coupled NASTRAN and SC/Tetra Qin Yin Fan Software CRADLE Co., Ltd. ABSTRACT In SAE paper 2004-01-1345, author focused on how to use a steady state temperature
More informationEvaluation of OpenFOAM for CFD of turbulent flow in
!" #$&% ' (*)+,, -. /12 Evaluation of OpenFOAM for CFD of turbulent flow in water turbines Håkan NILSSON Chalmers Univ. of Technology, Göteborg, Sweden hani@chalmers.se Key words: CFD, Validation, Kaplan,
More informationComputational Fluid Dynamics
Computational Fluid Dynamics A Practical Approach Jiyuan Tu RMIT University, Australia Guan Heng Yeoh Australian Nuclear Science and Technology Organisation Chaoqun Liu University of Texas, Arlington ^fl
More informationCFD with OpenSource software
CFD with OpenSource software The Harmonic Balance Method in foam-extend Gregor Cvijetić Supervisor: Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture University of Zagreb gregor.cvijetic@fsb.hr
More informationProcess Chemistry Toolbox - Mixing
Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws
More informationvector H. If O is the point about which moments are desired, the angular moment about O is given:
The angular momentum A control volume analysis can be applied to the angular momentum, by letting B equal to angularmomentum vector H. If O is the point about which moments are desired, the angular moment
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 3B: Conservation of Mass C3B: Conservation of Mass 1 3.2 Governing Equations There are two basic types of governing equations that we will encounter in this course Differential
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationTurbulence Modeling of Air Flow in the Heat Accumulator Layer
PIERS ONLINE, VOL. 2, NO. 6, 2006 662 Turbulence Modeling of Air Flow in the Heat Accumulator Layer I. Behunek and P. Fiala Department of Theoretical and Experimental Electrical Engineering Faculty of
More informationreport: Computational Fluid Dynamics Modelling of the Vortex Ventilator MK4 Rev 2 Ventrite International
report: Computational Fluid Dynamics Modelling of the Vortex Ventilator MK4 Rev 2 PrePAreD for: Ventrite International D O C U men t C O N trol Any questions regarding this document should be directed
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationPrediction of Performance Characteristics of Orifice Plate Assembly for Non-Standard Conditions Using CFD
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869, Volume-3, Issue-5, May 2015 Prediction of Performance Characteristics of Orifice Plate Assembly for Non-Standard Conditions
More informationMagnetic induction and electric potential solvers for incompressible MHD flows
Magnetic induction and electric potential solvers for incompressible MHD flows Alessandro Tassone DIAEE - Nuclear Section, Sapienza Università di Roma, Rome, Italy 05/12/2016 Alessandro Tassone CMHD Solvers
More informationCOMPUTATIONAL FLUID DYNAMIC ANALYSIS ON THE EFFECT OF PARTICLES DENSITY AND BODY DIAMETER IN A TANGENTIAL INLET CYCLONE HEAT EXCHANGER
THERMAL SCIENCE: Year 2017, Vol. 21, No. 6B pp. 2883-2895 2883 COMPUTATIONAL FLUID DYNAMIC ANALYSIS ON THE EFFECT OF PARTICLES DENSITY AND BODY DIAMETER IN A TANGENTIAL INLET CYCLONE HEAT EXCHANGER by
More informationHybrid RANS/LES simulations of a cavitating flow in Venturi
Hybrid RANS/LES simulations of a cavitating flow in Venturi TSFP-8, 28-31 August 2013 - Poitiers - France Jean Decaix jean.decaix@hevs.ch University of Applied Sciences, Western Switzerland Eric Goncalves
More informationDescription and validation of the rotordisksource class for propeller performance estimation
Description and validation of the rotordisksource class for propeller performance estimation Alexandre Capitao Patrao Department of Mechanics and Maritime Sciences Division of Fluid Dynamics Chalmers University
More informationEnhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins
Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins Jubin Jose 1, Reji Mathew 2 1Student, Dept. of Mechanical Engineering, M A
More informationInitial and Boundary Conditions
Initial and Boundary Conditions Initial- and boundary conditions are needed For a steady problems correct initial conditions is important to reduce computational time and reach convergence Boundary conditions
More informationA Simplified Numerical Analysis for the Performance Evaluation of Intercooler
, pp.34-38 http://dx.doi.org/10.14257/astl.2013.41.09 A Simplified Numerical Analysis for the Performance Evaluation of Intercooler Vashahi Foad, Myungjae Lee, Jintaek Kim, Byung Joon Baek* School of Mechanical
More informationCFD analysis of the transient flow in a low-oil concentration hydrocyclone
CFD analysis of the transient flow in a low-oil concentration hydrocyclone Paladino, E. E. (1), Nunes, G. C. () and Schwenk, L. (1) (1) ESSS Engineering Simulation and Scientific Software CELTA - Rod SC-41,
More informationViscoelastic Fluid Simulation with OpenFOAM
Viscoelastic Fluid Simulation with OpenFOAM Joan Marco Rimmek and Adrià Barja Romero Universitat Politècnica de Catalunya. BarcelonaTECH. (Dated: May 31, 2017) Viscoelastic fluids do not behave like newtonian
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationEXAMPLE 115 NON-ISOTHERMAL FLOW IN A SINGLE SCREW EXTRUDER DESCRIPTION KEYWORDS FILENAMES. Example 115
EXAMPLE 115 NON-ISOTHERMAL FLOW IN A SINGLE SCREW EXTRUDER DESCRIPTION In this example, we evaluate the flow in the last part of a single screw extruder, including the last few conveying screw elements,
More informationPrepared by: Simanto. Date: Sunday, August 17, 2014
C Date: Sunday, August 17, 2014 Report Contents Design 1 2 Scenario 1... 2 Materials... 2 boundary conditions... 3 Initial Conditions... 4 mesh... 4 Physics... 5 Solver Settings... 5 Convergence... 5 Results...
More informationECE580 Fall 2015 Solution to Midterm Exam 1 October 23, Please leave fractions as fractions, but simplify them, etc.
ECE580 Fall 2015 Solution to Midterm Exam 1 October 23, 2015 1 Name: Solution Score: /100 This exam is closed-book. You must show ALL of your work for full credit. Please read the questions carefully.
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationDrag Coefficient of Tall Building by CFD Method using ANSYS
Drag Coefficient of Tall Building by CFD Method using ANSYS Sharma P K 1, Dr. Parekar S R 2 1ME Student, Department of Civil Engineering, AISSM S, Pune, Maharashtra, India 2HOD, Department of Civil Engineering,
More information2Dimensional CFD Simulation of Gas Fired Furnace during Heat Treatment Process
2Dimensional CFD Simulation of Gas Fired Furnace during Heat Treatment Process Vivekanand D. Shikalgar 1, Umesh C.Rajmane 2, Sanjay S. Kumbhar 3, Rupesh R.bidkar 4, Sagar M. Pattar 5, Ashwini P. Patil
More informationSolving PDEs with OpenFOAM
Solving PDEs with OpenFOAM ThePDEswewishtosolveinvolvederivativesoftensorfieldswith respect to time and space The PDEs must be discretized in time and space before we solve them We will start by having
More informationDescription of poroussimplefoam and adding the Brinkmann model to the porous models Developed for OpenFOAM-2.2.x
CFD with OpenSource Software A course at Chalmers University of Technology Taught by Hakan Nilsson Project work: Description of poroussimplefoam and adding the Brinkmann model to the porous models Developed
More informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationJMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman
JMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman The exercises will be carried out on PC s in the practicum rooms. Several (Matlab and Fortran) files are required. How these can be obtained
More informationA Two-Fluid/DQMOM Methodology For Condensation In Bubbly Flow
A Two-Fluid/DQMOM Methodology For Condensation In Bubbly Flow Gothenburg Region OpenFOAM User Group Meeting, November 11, 2015 Klas Jareteg Chalmers University of Technology Division of Nuclear Engineering,
More informationABSTRACT I. INTRODUCTION
2017 IJSRSET Volume 3 Issue 3 Print ISSN: 2395-1990 Online ISSN : 2394-4099 Themed Section: Engineering and Technology CFD Analysis of Flow through Single and Multi Stage Eccentric Orifice Plate Assemblies
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationAngular momentum equation
Angular momentum equation For angular momentum equation, B =H O the angular momentum vector about point O which moments are desired. Where β is The Reynolds transport equation can be written as follows:
More informationCFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University
CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University email: bengt.sunden@energy.lth.se CFD? CFD = Computational Fluid Dynamics;
More informationABSTRACT I. INTRODUCTION
2016 IJSRSET Volume 2 Issue 4 Print ISSN : 2395-1990 Online ISSN : 2394-4099 Themed Section: Engineering and Technology Analysis of Compressible Effect in the Flow Metering By Orifice Plate Using Prasanna
More informationConstrained optimization
Constrained optimization In general, the formulation of constrained optimization is as follows minj(w), subject to H i (w) = 0, i = 1,..., k. where J is the cost function and H i are the constraints. Lagrange
More information