Flow Control of Curved Air Ducts using Topological Derivatives
|
|
- Clifford Townsend
- 6 years ago
- Views:
Transcription
1 8 th World Congress on Structural and Multidisciplinary Optimization June 1-5, 2009, Lisbon, Portugal Flow Control of Curved Air Ducts using Topological Derivatives Jonathan Chetboun 1,2, Grégoire Allaire 1 1 CMAP - École Polytechnique Palaiseau, France 2 Dassault-Aviation DGT/DTIAE/AERAV 78, quai Marcel Dassault, Saint Cloud, France Jonathan.Chetboun@dassault-aviation.com 1. Abstract A method for flow control of curved air ducts using mechanical vortex generators (VG) is described. This method includes the creation of vortex generators in a reference flow using the notion of topological derivative. A classical procedure of parametric optimization of the created VG by an adjoint state approach is also presented. Three-dimensional numerical simulations are performed illustrating the control of an s-shaped duct by VG. 2. Keywords: Topological derivative, Parametric optimization, Flow control, Vortex generators. 3. Introduction The control of aerodynamic flows is a crucial challenge in the aircraft and automotive industries. On military aircraft, the engine is supplied in air by a duct ; its length and shape are dictated by the global design of the airplane. The stealth constraints yield very curved ducts in order to hide the engine rotating blades, largely responsible for the radar signature of the aircraft. These levels of curvature cause flow separation harmful to the quality of the air feeding the engine, therefore imposing the use of control devices. Figure 1 shows the project of a European unmanned combat air vehicle (UCAV) called neuron and a generic u-shaped air duct. Figure 1: neuron and u-shaped air duct We are interested here in controling the flow inside a curved air duct by mechanical vortex generators. These are typically small aerodynamic surfaces attached to the duct wall and causing a swirl in their wake, thus delaying the appearance of separations. Rather than finely mesh the VG in numerical simulations, we prefer to use a model for their action on the fluid through the addition of a source term to the Navier- Stokes equations. Our goal is to determine whether it is relevant to introduce an additional VG into a reference flow in order to minimize a given cost function J. Note that one or several VG may already be present in this reference flow field. We use the idea of topological derivative, first introduced in the framework of structural optimum design [1], [2], [3], [4], [5], and subsequently applied, among other areas, to the incompressible Navier-Stokes equations [6]. In section 4, emphasis is made on the computation of the topological derivative for the creation of a source term in the case of a simple model problem. As a basis for our numerical simulations, we use Aether, a finite element industrial code developed at Dassault-Aviation, which solves the turbulent compressible Navier-Stokes equations [7]. Section 5 1
2 briefly describes the numerical method used in the code and its adjoint version obtained by automatic differentiation. We also deal with the modeling of vortex generators by source term. Section 6 is devoted to numerical simulations of the flow control of an s-shaped air duct using topological derivatives. In particular, we show the reduction of a given cost function consecutive to the placement of the VG according to the topological derivative map. Finally, in section 7, we apply a classical parametric optimization procedure, using the same adjoint state as in the topological derivative, in order to achieve a complete optimization of the created vortex generator. 4. Topological derivative of a source term As a vortex generator is modeled by a source term in numerical simulations, we need to adapt the concept of topological derivative developed in shape optimization to the case of the creation of a source term. The aim of this section is thus to compute the topological derivative of a source term in the simple case of the Laplace equation. The (formal) computation of this topological derivative in the more interesting case of the compressible Navier-Stokes system will be found in [8]. Let then Ω be an open subset of R N (N 2) and let us consider the following problem : { u = χf in Ω (1) u = 0 on Ω, with χ the characteristic function of the source term f. We define the cost function J : J(χ) = j(u)dx. (2) Ω Our goal is to evaluate the sensibility of the cost function J when a small perturbation χ is added to the characteristic function χ of the source term f. We thus define the following perturbed problem : { u = (χ + χ)f in Ω u (3) = 0 on Ω. The inclusion χ is defined from an elementary inclusion B (for example B = { x R N : x 1 } ) as follows : ( ) x x0 χ(x) = χ B. (4) Moreover, we assume that the supports of the original source term and of the inclusion do not overlap : χ(x)χ(x) = 0. (5) In the sequel, we need to evaluate the difference between the initial solution u and the perturbed one u. Thus, we define : v = u u, (6) which is solution of : { v = χf in Ω v = 0 on Ω. (7) We introduce the microscopic space variable y = x x0 in order to rescale the previous problem. This rescaling, centered on the inclusion, transforms the problem posed on Ω into a problem posed on R N, in the limit as tends to zero. Consequently, let us define w the solution of the canonical problem : { w = χb f(x 0 ) on R N (8) w 0 at, and the rescaled solution : ( ) x w (x) = 2 x0 w. (9) Then we have : { v = χ (f(x) f(x 0 )) w in Ω v = 0 on Ω. (10) 2
3 Using the decay properties at infinity of w, we can show that : v 2 L 2 (Ω) O ( N+1). (11) This estimate is not optimal, but sufficient for our purpose. One can refer to [8] for the detailed calculus. Introducing the adjoint state p solution of : { p = j (u) in Ω (12) p = 0 on Ω, and assuming that j (u) is bounded (for example by taking j quadratic), we can prove that : J(χ + χ) = J(χ) + N B f(x 0 )p(x 0 ) + O ( N+1). (13) In conclusion, the topological derivative G consecutive to the addition of a source term f at a point x 0 simply reads : G(x 0 ) = f(x 0 )p(x 0 ). (14) Consequently, in order to decrease the cost function J, and for a small, it is relevant to create a new source term f at the points where G is the more negative. The same kind of result can be obtained in the case of a source term f depending on the state u and when the perturbation lies at a distance of the boundary (which is the case for the source term modeling a vortex generator attached to the wall). We can also consider the cases of the Stokes problem and the incompressible Navier-Stokes equations. In the case of the compressible Navier-Stokes equations, only a formal calculus has been made. 5. Implementation of the topological derivative in a CFD industrial software 5.1. The Navier-Stokes solver Our numerical simulations are made using the software Aether, a compressible Navier-Stokes finite element solver developed at Dassault-Aviation. Let us first recall the expression of the compressible Navier-Stokes equations. Let ρ, u and e be respectively the density, the velocity and the total energy per unit mass of fluid. The Navier-Stokes equations read : ρ t + (ρu) = 0 ρu + (ρu u) = t σ + ρf (15) ρe + (ρeu) = (σu q) + ρf u, t where σ is the Cauchy stress tensor, q the heat-flux vector and f the source term. The numerical method used in the Aether solver is a Galerkin / least-squares formulation of the system written in terms of entropy variables and includes a dicontinuity capturing operator. The code also includes turbulent flow modeling through the Reynolds Averaged Navier-Stokes equations (RANS). The closure of the Reynolds stress tensor and heat flux is obtained using a classical Boussinesq hypothesis and the concept of eddyviscosity. The steady numerical simulations presented in the sequel are performed using a two-layer k ε turbulence model Modeling vortex generators In order to control the fluid inside the duct, we use passive devices called mechanical vortex generators. These are small blades attached to the duct wall and causing a swirl in their wake. This swirl prevents flow separation by giving energy to the boundary layer. As said in the introduction, vortex generators are not meshed in numerical simulations. Indeed, meshing an actual VG would require a too high level of mesh refinement. Moreover, one wants to test several configurations of vortex generators without the obligation of performing a different mesh. That is why we rather use a source term model in the code. The source term added to the Navier-Stokes equations Eq.(15) is defined as follows [9] : f = c S V G ρ (n u)(b u), (16) where c is a modeling constant, S V G is the surface of the vortex generator and vectors b and n are deduced from the orientation t of the vortex generator according to Figure 2. Therefore, the force f is 3
4 b x V G n t Figure 2: Geometrical description of a vortex generator perpendicular to the flow and tends to deviate it into a direction parallel to the vortex generator. In the code, the force is multiplied by an exponential decreasing with the distance from the vortex generator, in such a way that its action is distributed over a small volume surrounding the VG. Note that, as u, f is equal to zero at the wall Computation of the topological derivative We formally rewrite the compressible Navier-Stokes equations as : E(U, f(u, l)) = 0 (17) where U is the vector of entropy variables and l is the vector of the vortex generator parameters (its size and orientation). Then, we can write the adjoint system as follows [10], [11] : ( ) T ( ) T E J P = (18) U U This system is obtained by automatic differentiation of the Aether solver using the software Tapenade developed at INRIA [12]. It is clearly a discrete approach to the computation of the adjoint state rather than a discretization of the continuous equations. Note that only the Navier-Stokes part and not the turbulence part of the code has been differentiated. Numerical simulations show that this is sufficient to compute accurate sensitivities. The adjoint version of the code is validated by comparison with finite differences. As the source term modeling a vortex generator lies on a small volume around the VG, we rather use in our code an averaged version of the expression Eq. (14) : G(x 0 ) = f P dx. (19) V G The volume V G is taken as the support of the small source term we want to add. 6. Numerical simulations 6.1. The test case Our numerical simulations are performed on the generic s-shaped air duct with circular section presented on Figure 3. The duct is 1020 mms long. In the sequel, we make comparisons between two meshes : the coarse one includes approximately nodes and the fine one approximately nodes. First, we compute a steady reference flow using the Aether solver and a k ε turbulence model on half a mesh with symmetry conditions. The cost function considered here is the swirl at the duct output. It is the ratio between the velocity on the output plane and the velocity perpendicular to this plane Creation of a first vortex generator Figure 4 and Figure 5 show the topological derivative map on the two meshes for the creation of a 5-mm-long and 2-mm-high vortex generator. The created vortex generators are oriented at 45 degrees apart from flow direction. Circles represent results of numerical simulations with actual VG and show the decrease of cost function. Each circle comes with its corresponding value of cost function (Swirl) and topological derivative (TD). We can notice that, in the case of the coarse mesh, the localization VG2 presents a swirl smaller than the localization VG1 corresponding to the minimum of topological derivative. Nevertheless, this default is no longer visible in the case of the fine mesh. Besides, the swirl of localizations VG2 and VG3 escapes the classification predicted by the topological derivative. However, the topological derivative gives us a good idea of where to create a vortex generator (areas in blue and green on the two images, the cost function 4
5 Figure 3: Geometry of the test case decreases) and where not to create a vortex generator (areas in red, the cost function increases). Nevertheless, in order to keep decreasing the cost function, it is useful to implement a method of parametric optimization to adjust the size and orientation of the created vortex generator Addition of a second vortex generator Figure 6 and Figure 7 show the topological derivative map on the two meshes for the addtition of a second vortex generator. The first vortex generator has been put on the localization VG1 of Figure 4 and Figure 5. In this case, the topological derivative map is different between the two meshes. Unfortunately, the swirl resulting from numerical simulations with actual VG does not match the value of the topological derivative. The cost function increases wherever we put the second VG. However, in the case of the fine mesh, we can notice that this increase is more significant at the points where the topological derivative is positive, whereas it is very small where the topological derivative is negative. In this test case, where the flow separation is not massive at all, we can assume that a single vortex generator is sufficient to control the flow. It seems that the addition of a second VG is not suitable. 7. Parametric optimization In the previous section, we saw that the topological derivative allows us to know where it is relevant to create a vortex generator. Rather than a precise localization of its minimum, the topological derivative gives us an area where the cost function can be decreased. In order to find this minimum, it is useful to perform a classical method of parametric optimization, involving the same adjoint state that appeared in the computation of the topological derivative. Therefore, we optimize parameters such as the size or the orientation of the VG, which were fixed in the computation of the topological derivative. Let us first recall that, formally, the Aether software solves the following system of equations : and that the adjoint state equations can be written as : E(U, f(u, l)) = 0, (20) ( ) T ( ) T E J P =. (21) U U Then, it is well-known that the variation of cost function J resulting from the variation l of the parameters is given by : J = P T E l. (22) l The derivative in Eq. (22) is also computed using automatic differentiation of the code. We performed a parametric optimization of the vortex generator named VG1 given by the topological derivative on the coarse mesh using an algorithm of sequential quadratic programming. The five parameters we want 5
6 Figure 4: Topological derivative for the creation of a first VG - Coarse mesh to optimize are the two variables of the VG position on the surface of the duct, its angle (orientation compared with flow direction) and its size (height and length). Now, the cost function considered is J = 500 Swirl 2. Figure 8 shows the cost function and parameters evolution during the optimization process. In this picture, we present each step of the optimization process, which includes gradient computation and line search. That is why we can notice a brief increase in the cost function as the algorithm searches for the optimal step. In term of actual swirl, the creation of the first VG using topological derivative allows to reduce the cost function of 7 % whereas the parametric optimization produces a decrease of 3 % more. At the end of the optimization process, the position is not very different from the one given by the topological derivative. The orientation changes from 45 to approximately 30 degrees. The size of the VG is reduced from 5 x 2 mms to approximately 4 x 1.8 mms. 8. Conclusion An efficient method for the flow control inside a curved air duct using mechanical vortex generators has been presented. The idea of topological derivative has been used in order to decide where it is relevant to create a new VG in a reference flow. A classical method of parametric optimization by adjoint state has been applied successfully to the created vortex generator. Further developments are under way concerning the use of this method in the case of unsteady flows. Indeed, in order to achieve more complex simulations with shorter or more curved air ducts, it will be essential to take into account the unsteady nature of the flow with methods such as Detached Eddy Simulation. In particular, the computation of the adjoint state in this framework will be a crucial challenge, both from a mathematical and a numerical point of view. 9. References [1] H.A. Eschenauer, V.V. Kobelev and A. Schumacher, Bubble method for topology and shape optimization of structures, Structural Optimization, 8, 42-51, [2] J. Sokolowski and A. Zochowski, On the topological derivative in shape optimization, SIAM J. Control 6
7 Figure 5: Topological derivative for the creation of a first VG - Fine mesh Optim., 37, , [3] J. Céa, S. Garreau, P. Guillaume and M. Masmoudi, The shape and topological optimizations connection, Comput. Methods Appl. Mech. Engrg., 188, , [4] S. Garreau, P. Guillaume and M. Masmoudi, The topological asymptotic for PDE systems : the elasticity case, SIAM J. Control Optim., 39, , [5] P. Guillaume and K. Sid Idris, The topological asymptotic expansion for the Dirichlet problem, SIAM J. Control Optim., 41, , [6] S. Amstutz, The topological asymptotic for Navier-Stokes equations, ESAIM Control Optim. Calc. Var., 11, , [7] F. Chalot, M. Mallet and M. Ravachol, A comprehensive finite element Navier-Stokes solver for lowand high-speed aircraft design, AIAA, 814, [8] J. Chetboun, Design of aerodynamic shapes in presence of flow separations : control and optimization, Phd thesis, in preparation. [9] K. Waithe, Source term model for vortex generator vanes in a Navier-Stokes computer code, AIAA, 1236, [10] Q.V. Dinh, G. Rogé, C. Sevin and B. Stoufflet, Shape optimization in computational fluid dynamics, European Journal of Finite Elements, 5, , [11] L. Daumas, Q.V. Dinh, S. Kleinveld and G. Rogé, A CAD-based supersonic business jet optimization, ECCOMAS, [12] L. Hascoët, TAPENADE : a tool for Automatic Differentiation of programs, ECCOMAS,
8 Figure 6: Topological derivative for the addition of a second VG - Coarse mesh Figure 7: Topological derivative for the addition of a second VG - Fine mesh 8
9 COST COST ITERATION VORTEX! 1 POSITION-U, VORTEX! 1 POSITION-V VORTEX! 1 POSITION-U VORTEX! 1 POSITION-V ITERATION 5 VORTEX! 1 ANGLE 40 VORTEX! 1 ANGLE ITERATION VORTEX! 1 LENGTH, VORTEX! 1 HEIGHT VORTEX! 1 LENGTH VORTEX! 1 HEIGHT ITERATION Figure 8: Evolution of cost function and parameters during optimization process 9
Shape optimization for aerodynamic design: Dassault Aviation challenges and new trends
Direction Générale Technique Shape optimization for aerodynamic design: Dassault Aviation challenges and new trends Forum TERATEC 014 Atelier «Conception numérique optimale des systèmes complexes : état
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 informationAnalysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization
Analysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization Final Presentation Dylan Jude Graduate Research Assistant University of Maryland AMSC 663/664 May 4, 2017
More informationNUMERICAL SIMULATIONS OF FLOW THROUGH AN S-DUCT
NUMERICAL SIMULATIONS OF FLOW THROUGH AN S-DUCT 1 Pravin Peddiraju, 1 Arthur Papadopoulos, 2 Vangelis Skaperdas, 3 Linda Hedges 1 BETA CAE Systems USA, Inc., USA, 2 BETA CAE Systems SA, Greece, 3 CFD Consultant,
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 informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationCFD design tool for industrial applications
Sixth LACCEI International Latin American and Caribbean Conference for Engineering and Technology (LACCEI 2008) Partnering to Success: Engineering, Education, Research and Development June 4 June 6 2008,
More informationEstimating Vertical Drag on Helicopter Fuselage during Hovering
Estimating Vertical Drag on Helicopter Fuselage during Hovering A. A. Wahab * and M.Hafiz Ismail ** Aeronautical & Automotive Dept., Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More 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 informationDimensioning and Airflow Simulation of the Wing of an Ultralight Aircraft
Dimensioning and Airflow Simulation of the Wing of an Ultralight Aircraft Richárd Molnár 1 Gergely Dezső 2* Abstract: Increasing interest to ultralight aircrafts usually made of composite materials leads
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
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 informationNumerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models
Numerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models D. G. Jehad *,a, G. A. Hashim b, A. K. Zarzoor c and C. S. Nor Azwadi d Department of Thermo-Fluids, Faculty
More informationANSYS FLUENT. Airfoil Analysis and Tutorial
ANSYS FLUENT Airfoil Analysis and Tutorial ENGR083: Fluid Mechanics II Terry Yu 5/11/2017 Abstract The NACA 0012 airfoil was one of the earliest airfoils created. Its mathematically simple shape and age
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 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 informationCOUPLING THE LEVEL SET METHOD AND THE TOPOLOGICAL GRADIENT IN STRUCTURAL OPTIMIZATION
COUPLING THE LEVEL SET METHOD AND THE TOPOLOGICAL GRADIENT IN STRUCTURAL OPTIMIZATION Grégoire ALLAIRE Centre de Mathématiques Appliquées (UMR 7641) Ecole Polytechnique, 91128 Palaiseau, France gregoire.allaire@polytechnique.fr
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 informationKeywords: 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 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 informationApplication of Wray-Agarwal Turbulence Model for Accurate Numerical Simulation of Flow Past a Three-Dimensional Wing-body
Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 4-28-2016 Application
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationExtension and Validation of the CFX Cavitation Model for Sheet and Tip Vortex Cavitation on Hydrofoils
Extension and Validation of the CFX Cavitation Model for Sheet and Tip Vortex Cavitation on Hydrofoils C. Lifante, T. Frank, M. Kuntz ANSYS Germany, 83624 Otterfing Conxita.Lifante@ansys.com 2006 ANSYS,
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 informationModeling icing using cartesian grids, penalization & level sets
Introduction State of the art The proposed approach Proof of concept Modeling icing using cartesian grids, penalization & level sets Héloïse Beaugendre IMB Institut de Mathématiques de Bordeaux - INRIA
More informationEstimation of Flow Field & Drag for Aerofoil Wing
Estimation of Flow Field & Drag for Aerofoil Wing Mahantesh. HM 1, Prof. Anand. SN 2 P.G. Student, Dept. of Mechanical Engineering, East Point College of Engineering, Bangalore, Karnataka, India 1 Associate
More informationProgram: Advanced Certificate Program
Program: Advanced Certificate Program Course: CFD-Vehicle Aerodynamics Directorate of Training and Lifelong Learning #470-P, Peenya Industrial Area, 4th Phase Peenya, Bengaluru 560 058 www.msruas.ac.in
More informationSPC 307 Aerodynamics. Lecture 1. February 10, 2018
SPC 307 Aerodynamics Lecture 1 February 10, 2018 Sep. 18, 2016 1 Course Materials drahmednagib.com 2 COURSE OUTLINE Introduction to Aerodynamics Review on the Fundamentals of Fluid Mechanics Euler and
More informationAnalysis of Flow Dynamics of an Incompressible Viscous Fluid in a Channel
Analysis of Flow Dynamics of an Incompressible Viscous Fluid in a Channel Deepak Kumar Assistant Professor, Department of Mechanical Engineering, Amity University Gurgaon, India E-mail: deepak209476@gmail.com
More informationScuola Politecnica DIME
Scuola Politecnica DIME Ingegneria Meccanica - Energia e Aeronautica Anno scolastico 2017-2018 Fluidodinamica Avanzata Aircraft S-shaped duct geometry optimization Professor Jan Pralits Supervisor Joel
More informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More 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 informationLES Analysis on Shock-Vortex Ring Interaction
LES Analysis on Shock-Vortex Ring Interaction Yong Yang Jie Tang Chaoqun Liu Technical Report 2015-08 http://www.uta.edu/math/preprint/ LES Analysis on Shock-Vortex Ring Interaction Yong Yang 1, Jie Tang
More informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationGradient-based Optimization of Parameterized CAD Geometries
11 th World Congress on Structural and Multidisciplinary Optimization 7 th - 12 th, June 2015, Sydney Australia Gradient-based Optimization of Parameterized CAD Geometries Timothée Leblond 1,2, Pierre
More informationTurbulent 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 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 informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationMcNair Scholars Research Journal
McNair Scholars Research Journal Volume 2 Article 1 2015 Benchmarking of Computational Models against Experimental Data for Velocity Profile Effects on CFD Analysis of Adiabatic Film-Cooling Effectiveness
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY. Analyzing wind flow around the square plate using ADINA Project. Ankur Bajoria
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Analyzing wind flow around the square plate using ADINA 2.094 - Project Ankur Bajoria May 1, 2008 Acknowledgement I would like to thank ADINA R & D, Inc for the full
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationHigh-order mesh generation for CFD solvers
High-order mesh generation for CFD solvers M. Turner, D. Moxey, S. Sherwin, J. Peiró Department of Aeronautics, Imperial College London DiPaRT 2015 Annual Meeting, Bristol, UK 17 th November 2015 Overview
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 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 informationThe viscous forces on the cylinder are proportional to the gradient of the velocity field at the
Fluid Dynamics Models : Flow Past a Cylinder Flow Past a Cylinder Introduction The flow of fluid behind a blunt body such as an automobile is difficult to compute due to the unsteady flows. The wake behind
More informationInternational Power, Electronics and Materials Engineering Conference (IPEMEC 2015)
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Numerical Simulation of the Influence of Intake Grille Shape on the Aerodynamic Performance of a Passenger Car Longwei
More informationDesign Optimization of a Subsonic Diffuser. for a Supersonic Aircraft
Chapter 5 Design Optimization of a Subsonic Diffuser for a Supersonic Aircraft 5. Introduction The subsonic diffuser is part of the engine nacelle leading the subsonic flow from the intake to the turbo-fan
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 informationThree Dimensional Numerical Simulation of Turbulent Flow Over Spillways
Three Dimensional Numerical Simulation of Turbulent Flow Over Spillways Latif Bouhadji ASL-AQFlow Inc., Sidney, British Columbia, Canada Email: lbouhadji@aslenv.com ABSTRACT Turbulent flows over a spillway
More informationAirfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition
Airfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition.5.5.5.5.5.5.5..5.95.9.85.8.75.7 Patrick A. LeGresley and Juan J. Alonso Dept. of Aeronautics & Astronautics
More informationNumerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Numerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind Fan Zhao,
More informationTHE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS
March 18-20, 2013 THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS Authors: M.R. Chiarelli, M. Ciabattari, M. Cagnoni, G. Lombardi Speaker:
More informationApplication of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution
Application of A Priori Error Estimates for Navier-Stokes Equations to Accurate Finite Element Solution P. BURDA a,, J. NOVOTNÝ b,, J. ŠÍSTE a, a Department of Mathematics Czech University of Technology
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 informationFAR-Wake Workshop, Marseille, May 2008
Wake Vortices generated by an Aircraft Fuselage : Comparison of Wind Tunnel Measurements on the TAK Model with RANS and RANS-LES Simulations T. Louagie, L. Georges & P. Geuzaine Cenaero CFD-Multiphysics
More informationEffects of bell mouth geometries on the flow rate of centrifugal blowers
Journal of Mechanical Science and Technology 25 (9) (2011) 2267~2276 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-011-0609-3 Effects of bell mouth geometries on the flow rate of centrifugal
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 informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationcomputational Fluid Dynamics - Prof. V. Esfahanian
Three boards categories: Experimental Theoretical Computational Crucial to know all three: Each has their advantages and disadvantages. Require validation and verification. School of Mechanical Engineering
More informationRANS 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 informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More 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 informationAERODYNAMIC OPTIMIZATION OF A FORMULA STUDENT CAR
AERODYNAMIC OPTIMIZATION OF A FORMULA STUDENT CAR 1 Argyrios Apostolidis *, 2 Athanasios Mattas, 3 Aggelos Gaitanis, 4 Nikolaos Christodoulou 1 Aristotle Racing Team, Greece, 4 BETA CAE Systems S.A., Greece
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More informationNumerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert
Numerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert A. Kadyyrov 1 1 Research center for power engineering problems Federal government budgetary institution of science Kazan
More informationNumerical Simulation of Coastal Wave Processes with the Use of Smoothed Particle Hydrodynamics (SPH) Method
Aristotle University of Thessaloniki Faculty of Engineering Department of Civil Engineering Division of Hydraulics and Environmental Engineering Laboratory of Maritime Engineering Christos V. Makris Dipl.
More informationGEOMETRICAL CONSTRAINTS IN THE LEVEL SET METHOD FOR SHAPE AND TOPOLOGY OPTIMIZATION
1 GEOMETRICAL CONSTRAINTS IN THE LEVEL SET METHOD FOR SHAPE AND TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique Results obtained in collaboration with F. Jouve (LJLL, Paris 7), G. Michailidis
More informationTurbulence 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 informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
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 informationFlow and likely scour around three dimensional seabed structures evaluated using RANS CFD
Flow and likely scour around three dimensional seabed structures evaluated using RANS CFD By Guillaume de Hauteclocque, Justin Dix, David Lambkin, Stephen Turnock University of Southampton Ship science
More informationSHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY
SHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY A. N. Ryabinin Department of Hydroaeromechanics, Faculty of Mathematics and Mechanics, Saint-Petersburg State University, St. Petersburg, Russia E-Mail: a.ryabinin@spbu.ru
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 informationAERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET
25 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET Qingzhen YANG*,Yong ZHENG* & Thomas Streit** *Northwestern Polytechincal University, 772,Xi
More informationS-ducts and Nozzles: STAR-CCM+ at the Propulsion Aerodynamics Workshop. Peter Burns, CD-adapco
S-ducts and Nozzles: STAR-CCM+ at the Propulsion Aerodynamics Workshop Peter Burns, CD-adapco Background The Propulsion Aerodynamics Workshop (PAW) has been held twice PAW01: 2012 at the 48 th AIAA JPC
More informationALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS ALE METHODS FOR DETERMINING STATIONARY SOLUTIONS OF METAL FORMING PROCESSES
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationModeling Unsteady Compressible Flow
Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial
More informationAIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS
AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS Emre Alpman Graduate Research Assistant Aerospace Engineering Department Pennstate University University Park, PA, 6802 Abstract A new methodology
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationHigh-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder
High-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder Aerospace Application Areas Aerodynamics Subsonic through Hypersonic Aeroacoustics Store release & weapons bay analysis High lift devices
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 informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationFinal drive lubrication modeling
Final drive lubrication modeling E. Avdeev a,b 1, V. Ovchinnikov b a Samara University, b Laduga Automotive Engineering Abstract. In this paper we describe the method, which is the composition of finite
More informationSimulation of Turbulent Axisymmetric Waterjet Using Computational Fluid Dynamics (CFD)
Simulation of Turbulent Axisymmetric Waterjet Using Computational Fluid Dynamics (CFD) PhD. Eng. Nicolae MEDAN 1 1 Technical University Cluj-Napoca, North University Center Baia Mare, Nicolae.Medan@cunbm.utcluj.ro
More informationCAD-BASED WORKFLOWS. VSP Workshop 2017
CAD-BASED WORKFLOWS VSP Workshop 2017 RESEARCH IN FLIGHT COMPANY Established 2012 Primary functions are the development, marketing and support of FlightStream and the development of aerodynamic solutions
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 informationNumerical 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 informationCFD Analysis of a Fully Developed Turbulent Flow in a Pipe with a Constriction and an Obstacle
CFD Analysis of a Fully Developed Turbulent Flow in a Pipe with a Constriction and an Obstacle C, Diyoke Mechanical Engineering Department Enugu State University of Science & Tech. Enugu, Nigeria U, Ngwaka
More informationIntroduction to Aerodynamic Shape Optimization
Introduction to Aerodynamic Shape Optimization 1. Aircraft Process 2. Aircraft Methods a. Inverse Surface Methods b. Inverse Field Methods c. Numerical Optimization Methods Aircraft Process Conceptual
More informationNUMERICAL SIMULATION OF FLOW FIELD IN AN ANNULAR TURBINE STATOR WITH FILM COOLING
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL SIMULATION OF FLOW FIELD IN AN ANNULAR TURBINE STATOR WITH FILM COOLING Jun Zeng *, Bin Wang *, Yong Kang ** * China Gas Turbine Establishment,
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 informationAERODYNAMIC DESIGN AND OPTIMIZATION TOOLS ACCELERATED BY PARAMETRIC GEOMETRY PREPROCESSING
1 European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 ECCOMAS AERODYNAMIC DESIGN AND OPTIMIZATION TOOLS ACCELERATED BY PARAMETRIC
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
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 informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
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 informationAn Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm
An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More information