Store Separation Simulation using Oct-tree Grid Based Solver
|
|
- Phillip Gilmore
- 6 years ago
- Views:
Transcription
1 SAROD Symposium on Applied Aerodynamics and Design of Aerospace Vehicles (SAROD 2009) December 10-12, 2009, Bengaluru, India Store Separation Simulation using Oct-tree Grid Based Solver Saurabh Pandey and Bharat R. Agrawal Zeus Numerix Private Limited ABSTRACT Simulating moving solid body in a static Cartesian Oct-tree grid is the topic of discussion. Previously many implementations of oct-tree based adaptive grid to static steady state problems have been demonstrated. This work enhances the scope of oct-tree grid implementation to unsteady problems. Initially a compressible Euler inviscid solver is validated over some standard test cases. Grid is then modified to take into account moving geometry with impermeable boundaries and maintaining the conservation because of the movement. Dynamics of solid body is simulated wherein the trajectory of the body is predicted using aerodynamic forces estimated from the CFD simulation. Solver s results are discussed at the end for the well known store separation problem. Keywords: Cartesian Grid, Compressible Solver, Oct Tree, Store Separation, Unsteady NOMENCLATURE CV = Control Volume CS = Boundary of the control volume U = Vector of conserved variables F = Vector of flux of conserved variables ρ = Density within a control volume u = x component of velocity of fluid in control volume v = y component of velocity of fluid in control volume w = z component of velocity of fluid in control volume e o = Total energy of the fluid in control volume v n = Velocity normal to the boundary of a control volume
2 Store Separation Simulation using Oct-tree Grid Based Solver 445 p = Pressure n x,n y,n z = x, y and z components of outward normal unit vector to the boundary h o = Total enthalpy e = Static internal energy h = Static enthalpy a = Speed of sound γ = Ratio of specific heat capacities C p = Specific heat capacity of a gas at constant pressure C v = Specific heat capacity of a gas at constant volume 1. INTRODUCTION CFD has started playing a crucial role in arriving at designs with required expectations. Its extensive use has still not seen the horizon because of the time constraints of advanced CFD analysis. Over and above if the analysis has to be done on large number of possible design configurations the problem of time constraints of simulation makes it impossible in the existing scenario. If the simulation time could be reduced to make such numerous analyses possible then non-optimum configurations could be easily discarded in an early stage. This would save a substantial amount of human and machine time as other analysis such as structural, etc. would be done on lesser number of possible configurations. The most time consuming step in a CFD analysis of a complex geometry is the grid generation process as the solver generally are already utilizing power of parallel computing. Also, this time required increases with the complexity of geometry and after a certain level of complexity grid generation becomes almost impossible. This creates a need for an automatic grid generation tool which can be readily utilized in creating grids over various configurations under consideration in the design process. Furthermore, this grid generation should be such that the time taken in creating grid is independent from geometric complexity. Current methodologies of grid generation are largely based on boundary fitted methods. These methods are highly dependent on geometry surface and hence limited by the complexity of geometry. In order to remove this dependence, the other method which is boundary immersed method is prescribed. These methods start with a domain which is filling the space where the analysis needs to be done and progressively creates discretized cells towards the geometry surface. This creates the flexible framework where a small change in geometry configuration does not require restart of whole grid generation process [4]. Oct-tree based grid generation is a popular boundary immersed method [1-4]. The hallmark of oct-tree grids is that they are recursive in nature. Hence cells are created by division of bigger cell only in a region where smaller cells are required to capture the gradients. Other significant features of oct-tree grids includes simplifying the advanced grid generation algorithms like adaptive mesh refinement (AMR), domain partitioning, changing cells near the geometry on the fly, etc [6].
3 446 Saurabh Pandey and Bharat R. Agrawal Having created the grid, a new solver is required to extract the maximum advantages from this grid. Finally, to completely solve the problem of store separation, a coupled system is needed which must comprise a highly dynamic grid generation module, a compressible solver for CFD and a 6DOF model for dynamics [8]. 2. GRID GENERATOR A grid generator module previously developed and tested with Euler incompressible solver is used in the current work [11]. Grid generation using oct-tree method over a complex geometry is constrained with three requirements from the solver. Firstly, the ratio of volume between two neighbouring cells should be at most eight or at least one-eight. This constraint is implemented to ease the gradient calculations. The second important constraint in the oct-tree grid is that the uncut volume of all neighbouring cells around a cut-cell should be same as uncut volume of the cut-cell itself. This is an essential requirement for quick merging of small cut-cells. Finally, volume and shape of cut-cell needs to be extracted with maximum possible accuracy so that the volume conservation is held near the geometry and the computations done at the boundary are accurate [5]. 2.1 Merging of cut-cells Oct-tree grid is infected with very well known Small Cell problem wherein as the geometry is inserted inside the grid it can intersect some cells in such a way that a very small cell is left in the computational domain. Such cells being very small in volume adversely affect the time stepping of the solver hence greatly increasing the simulation time [7]. Merging of such small cells is one of the well known remedies. It is done by choosing cells around the selected small cell and merging with it. The basis of selection is the area of the face between the cell selected for merging and the small cell. Neighbouring cell sharing the largest area with small cell is selected and merged. This process is continued as long as the merged cell s volume is above a threshold limit. 2.2 Maintenance of volume conservation Another important concern for simulating a moving solid boundary problem in a static Cartesian mesh is the conservation of volume. In order to get a judgment of change in volume of cell on its division volume of cell pre as well as post division is calculated and an error is evaluated. It is in general found that these operations don t produce significant change in the volume. Therefore cell division operations indeed conserve the volume. 3. COMPRESSIBLE SOLVER Finite Volume Method (FVM) based compressible solver is developed on these grids. The governing equation used is Euler inviscid which for a cell whose faces are not necessarily aligned along the x, y or z axis can be written as follows
4 Store Separation Simulation using Oct-tree Grid Based Solver 447 where,, and The velocity at the cell boundary is defined as the dot product of local velocity vector and unit outward normal to the boundary The total energy and enthalpy per unit mass in the above equations is composed of static and dynamic parts Static energy, enthalpy and pressure can be expressed as a function of local speed of sound, temperature of the fluid and the ratio of specific heats. where,,,, Flux schemes due to van Leer [9] and HLLC [10] are written and validated. The solver is made capable to account for AMR by maintaining conservation while such adaptive refinements are performed on the grid. Grid refinement is based on the difference in one of the primitive variables in neighbouring cells. If this difference in primitive is above a preset threshold value refinement is performed while maintain conservation. De-refinement on the other side is done by calculating the variance of a primitive variable in all child cells. If the variance is not large, again depending on a preset limit, then all the child cells are merged together to results in its parent cell. Needless to say that all the conserved distributed among children are summed up and given to the parent on derefinement. 4. DYNAMICS IMPLEMENTATION Final implementation in the setup to simulate store separation specifically or any unsteady flow phenomena involving motion of geometry with impermeable boundary is to incorporate dynamics. It involves calculation of aerodynamic forces on the solid body in consideration using CFD solver and thereafter moving the geometry with respect to the grid under the action of the aerodynamic and/or other body forces.
5 448 Saurabh Pandey and Bharat R. Agrawal Dynamics cycle is discussed here on store separation problem for simplicity and coherence but it is in general applicable for any relatively moving solid body simulation problem. The whole process is started by first simulating the aircraft with stores in place to its steady state. With the steady state results in hand the process of separation of store is started by detaching store geometry from the grid. This is achieved by removing cut information of store geometry from the grid. Steady state result obtained initially is used to calculate the aerodynamic forces on the store which further is used to move the store in 3D space. Grid is regenerated at this new moved location of the store and solver is now executed over this new configuration for few iterations. This in turn gives new forces on the store and it is again transported to a new location based on the dynamics. This cycle is repeated for a specified amount of time. Thus solving store separation problem. The main restrictive criterion here is the time for which the store is displaced is such that no part of the store geometry crosses more than one face in the computational domain. This is implemented to make the process of maintenance of conservation simple. It can also be viewed as a CFL type criterion for simulating dynamics. 4.1 Maintenance of mass, momentum and energy conservation Solid body movement relative to the grid brings new cells in the computational domain at the same time erases few cells from the fluid domain. Also adaptive refinement of cells as well as a cell s de-refinement creates and deletes cells from the domain. Mass, momentum and energy conservation has to be maintained in cells like these to correctly capture the physics of the problem. Cells which go inside the solid body and hence out of the fluid domain have to distribute the conserved quantities in the neighbouring cells which will remain in the domain even after the motion of the solid body. Hence conserved quantities from such cells are distributed among neighbouring cells in the ratio of interface area. Similarly for newly created cells, in order to maintain global conservation, are required to have ideally no conserved quantities. But zero initialization of conserved quantities will cause singularity in the solver as a result of which such newly born cells are initialized with infinitesimally small amount values. Refinement and de-refinement operations on cells are made conservative by interpolating the conserved data from old cell(s) to the new ones and then converting them to their primitive form. 5. DISCUSSION OF TEST RESULTS In this section 4 test cases are discussed. Starting with the validation of supersonic flow over a wedge and a forward step the discussion would move to the results obtained over a realistic geometry of an aircraft and the discussion ends with the dynamic simulation of store separating from the same aircraft.
6 5.1 Supersonic flow over a wedge of 10 angle Store Separation Simulation using Oct-tree Grid Based Solver 449 One of the basic validation cases for any compressible solver is how close the angle of a capture shock wave is. For this purpose, a flow of Mach number 2 is considered over a wedge of 10 angle with the flow direction. This problem is solved using flux splitting scheme due to van Leer [9]. The results are shown in Figures 1 and 2. Figure 1 shows the grid for solution. This case was initialized with a uniform Cartesian base grid and the clustering along the shocks appears completely due to the adaptive framework of the solver. Figure 2 shows the contour plot of the Mach number in the flow field. Comparison of obtained data with theoretical is as follows: Theoretical solution: shock angle = 39.31, pressure ratio across the shock wave = Computed solution: shock angle = ~40, pressure ratio across the shock wave = Figure1: Solution grid for wedge test case Figure2: Contour plot of Mach number for wedge test case 5.2 Supersonic flow over forward facing step In the next validation case, a flow of Mach number 3 inside a 2D duct with a forward facing step was considered. This test case was solved using HLLC scheme. Figure 3 shows the pressure contour plot from literature [9] and Figure 4 shows the pressure contour plot obtained using the present solver. In the obtained solution the maximum pressure is reached inside the Mach disc at the top boundary and the ratio of this pressure with upstream pressure is obtained as which match closely with the literature results.
7 450 Saurabh Pandey and Bharat R. Agrawal Figure3: Contour plot of pressure from literature for forward step problem Figure4: Contour plot of pressure for forward step problem 5.3 Supersonic flow over an aircraft To demonstrate the speed of the solver, supersonic flow at Mach number 2.0 is considered over aircraft geometry at angle of attack of 5. Figure 5 and 6 shows the pressure profile over the upper and lower surfaces of the wing respectively. It can be noted that red and blue represents highest and lowest values respectively. Thus, the differences in colors on these surfaces suggest a positive lift over this geometry. This is a converged solution and was obtained in less than 15 minutes on a standard modern day desktop machine. Figure5: Pressure over the upper surface of the aircraft
8 Store Separation Simulation using Oct-tree Grid Based Solver 451 Figure6: Pressure over the lower surface of the aircraft 5.4 Store separation Aircraft and generic store geometry was simulated at supersonic flow of Mach number 2.0. Figure 7 shows the pressure variation over both the geometries after the steady state has been reached by the solver. In this figure a blue patch is clearly seen on the lower surface of the aircraft with. It can be noted that most of the surface of aircraft is blue which is due to the fact that this region lies outside the actual computational domain and does not have any pressure values associated with it. Figure7: Pressure contours after reaching steady state Figure 8, 9 and 10 shows the pressure variations and the location of the store at 10 iterations of the dynamics solver. It can be seen that from Figure 7 to figure 8, the blue patch which is a low pressure area grows in size and then in subsequent figures it moves towards the trailing edge and fades.
9 452 Saurabh Pandey and Bharat R. Agrawal Figure8: Pressure contours after running 10 iterations of dynamics solver Figure9: Pressure contours after running 20 iterations of dynamics solver Figure10: Pressure contours after running 30 iterations of dynamics solver 6. CONCLUSION AND FUTURE SCOPE An oct-tree based adaptive compressible solver is validated against standard test cases of flow over a wedge and forward facing step problems. A technique for simulating moving solid body in a static Cartesian grid is demonstrated. The whole unsteady simulation is
10 Store Separation Simulation using Oct-tree Grid Based Solver 453 broken down into four main steps starting with the CFD simulation of the body and calculating the aerodynamic forces. Next is detaching the moving solid body from the grid. Third step involves movement of the body by a small amount based on the forces calculated in first step. Finally the body is re-immersed in the grid after the displacement. Most obvious step ahead from the given point is implementing Arbitrary Lagrangian Eulerian (ALE) formulation to account for the moving boundary of a solid body with respect to the computational domain. Further the flux calculation scheme has to be modified accordingly to account for moving boundaries. ACKNOWLEDGMENT We would like to thank Prof. G. R. Shevare of IIT Bombay for his valuable guidance and Zeus Numerix for providing support for the development. REFERENCES [1] Aftosmis, M.J., Melton, J.E., and Berger, M.J., Adaptive Cartesian Mesh Generation. Chapter 22 in Handbook of Grid Generation, Thompson, J, Weatherhill, N., and Soni, B. eds. CRC Press [2] Berger, M.J., and Aftosmis, M.J., Aspects (and aspect ratios) of Cartesian mesh methods. Proceedings of the 16 th International Conference on Numerical Methods in Fluid Mechanics, Arcachon, France, July [3] Aftosmis, M.J., Solution adaptive Cartesian grid methods for aerodynamic flows with complex geometries, von Karman Institute for Fluid Dynamics, Lecture Series , Rhode- Saint-Genèse, Belgium, Mar. 3-7, [4] Aftosmis, M.J., Berger, M.J., Melton, J.E., Robust and efficient Cartesian mesh generation for component based geometry, AIAA Paper , Jan [5] Berger, M.J., and Melton, J.E., An Accuracy Test of a Cartesian Grid Method for Steady flow in Complex Geometries, Proc. 8thrIntl. Conf. Hyp. Problems, Uppsala, Stonybrook, NY, Jun., also RIACS Report [6] De Zeeuw, D., and Powell, K., An Adaptively Refined Cartesian Mesh Solver for the Euler Equations, AIAA Paper , [7] Popinet, S., Gerris: A tree based adaptive solver for the incompressible euler equations in complex geometries. J. Comp. Phys. 190, , [8] Murman, S.M., Aftosmis, M.J., and Berger, M.J., Implicit Approaches for Moving Boundaries in a 3-D Cartesian Method, AIAA , [9] van Leer, B., Flux-vector splitting for the Euler equations, Lecture Notes in Physics.170. Springer, Berlin. pp , [10] Toro, E.F., Spruce, M., and Speares, W., Restoration of the contact surface in the HLL- Riemann solver, Shock Waves, Volume 4, pp , [11] Agrawal, B.R., and Pandey, S., Revisiting Projection Methods over Automatic Oct-tree Meshes, IISc Centenary International Conference and Exhibition on Aerospace Engineering, pp , 2009.
Estimation 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 informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
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 informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More 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 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 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 informationParallel Multigrid on Cartesian Meshes with Complex Geometry +
Parallel Multigrid on Cartesian Meshes with Complex Geometry + Marsha Berger a and Michael Aftosmis b and Gedas Adomavicius a a Courant Institute, New York University, 251 Mercer St., New York, NY 10012
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 informationValidation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons Stores Separation Introduction Flight Test Expensive, high-risk, sometimes catastrophic
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 informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationAvailable online at ScienceDirect. Procedia Engineering 99 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 99 (2015 ) 575 580 APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology, APISAT2014 A 3D Anisotropic
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationA STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.
More 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 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 informationSupersonic Flow Over a Wedge
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 2 Supersonic Flow Over a Wedge Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification A uniform supersonic stream encounters a wedge
More informationCFD Post-Processing of Rampressor Rotor Compressor
Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements
More informationComputation of Fictitious Gas Flow with Euler Equations
1 Computation of Fictitious Gas Flow with Euler Equations Pei Li and Helmut Sobieczky DLR Göttingen, Germany Abstract The Fictitious Gas Concept supports some computational design methods to construct
More informationAnalysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions
2010 NDIA GROUND VEHICLE SYSTEMS ENGINEERING AND TECHNOLOGY SYMPOSIUM MODELING & SIMULATION, TESTING AND VALIDATION (MSTV) MINI-SYMPOSIUM AUGUST 17-19 DEARBORN, MICHIGAN Analysis Comparison between CFD
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationCOMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY. Abstract. I. Introduction
COMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY John R. Cipolla 709 West Homeway Loop, Citrus Springs FL 34434 Abstract A series of computational fluid dynamic (CFD)
More informationRecent & Upcoming Features in STAR-CCM+ for Aerospace Applications Deryl Snyder, Ph.D.
Recent & Upcoming Features in STAR-CCM+ for Aerospace Applications Deryl Snyder, Ph.D. Outline Introduction Aerospace Applications Summary New Capabilities for Aerospace Continuity Convergence Accelerator
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More 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 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 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 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 informationCFD Analysis of conceptual Aircraft body
CFD Analysis of conceptual Aircraft body Manikantissar 1, Dr.Ankur geete 2 1 M. Tech scholar in Mechanical Engineering, SD Bansal college of technology, Indore, M.P, India 2 Associate professor in Mechanical
More informationExpress Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 -
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationA MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES
A MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES Khalil Bensassi, and Herman Deconinck Von Karman Institute for Fluid Dynamics Aeronautics & Aerospace
More informationGrid. Apr 09, 1998 FLUENT 5.0 (2d, segregated, lam) Grid. Jul 31, 1998 FLUENT 5.0 (2d, segregated, lam)
Tutorial 2. Around an Airfoil Transonic Turbulent Flow Introduction: The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a non-zero angle of attack. You will use the
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 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 informationCoupled Analysis of FSI
Coupled Analysis of FSI Qin Yin Fan Oct. 11, 2008 Important Key Words Fluid Structure Interface = FSI Computational Fluid Dynamics = CFD Pressure Displacement Analysis = PDA Thermal Stress Analysis = TSA
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 informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationAbstract. Introduction
EULER SOLUTIONS AS LIMIT OF INFINITE REYNOLDS NUMBER FOR SEPARATION FLOWS AND FLOWS WITH VORTICES Wolfgang Schmidt and Antony Jameson Dornier GmbH, D-7990 Friedrichshafen, FRG and Princeton University,
More informationA Hybrid Cartesian Grid and Gridless Method for Compressible Flows
rd AIAA Aerospace Sciences Meeting and Exhibit,, January 5, Reno, Nevada A Hybrid Cartesian Grid and Gridless Method for Compressible Flows Hong Luo and Joseph D. Baum Science Applications International
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationModeling & Simulation of Supersonic Flow Using McCormack s Technique
Modeling & Simulation of Supersonic Flow Using McCormack s Technique M. Saif Ullah Khalid*, Afzaal M. Malik** Abstract In this work, two-dimensional inviscid supersonic flow around a wedge has been investigated
More informationAn Adaptively-Refined Cartesian Mesh Solver for the Euler Equations
An Adaptively-Refined Cartesian Mesh Solver for the Euler Equations Darren De Zeeuw * Kenneth G. Powell t The University of Michigan Department of Aerospace Engineering Ann Arbor, MI 48109-2140 April,
More informationAIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN ABSTRACT
AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN Adeel Khalid *, Daniel P. Schrage + School of Aerospace Engineering, Georgia Institute of Technology
More informationThree dimensional meshless point generation technique for complex geometry
Three dimensional meshless point generation technique for complex geometry *Jae-Sang Rhee 1), Jinyoung Huh 2), Kyu Hong Kim 3), Suk Young Jung 4) 1),2) Department of Mechanical & Aerospace Engineering,
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 informationA cell by cell anisotropic adaptive mesh Arbitrary Lagrangian Eulerian method for the numerical solution of the Euler equations. J. M.
THE UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS A cell by cell anisotropic adaptive mesh Arbitrary Lagrangian Eulerian method for the numerical solution of the Euler equations. J. M. Morrell Thesis
More informationMissile External Aerodynamics Using Star-CCM+ Star European Conference 03/22-23/2011
Missile External Aerodynamics Using Star-CCM+ Star European Conference 03/22-23/2011 StarCCM_StarEurope_2011 4/6/11 1 Overview 2 Role of CFD in Aerodynamic Analyses Classical aerodynamics / Semi-Empirical
More informationAccurate and Efficient Turbomachinery Simulation. Chad Custer, PhD Turbomachinery Technical Specialist
Accurate and Efficient Turbomachinery Simulation Chad Custer, PhD Turbomachinery Technical Specialist Outline Turbomachinery simulation advantages Axial fan optimization Description of design objectives
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationAdaptive Mesh Refinement of Supersonic Channel Flows on Unstructured Meshes*
International Journal of Computational Fluid Dynamics, 2003 Vol. not known (not known), pp. 1 10 Barron Special Issue Adaptive Mesh Refinement of Supersonic Channel Flows on Unstructured Meshes* R.C. RIPLEY
More informationAirfoil shape optimization using adjoint method and automatic differentiation. Praveen. C
'th Annual AeSI CFD Symposium, -2 August 2009, Bangalore Airfoil shape optimization using adjoint method and automatic differentiation Praveen. C TIFR Center for Applicable Mathematics Post Bag No. 6503,
More informationGrid Dependence Study of Transonic/Supersonic Flow Past NACA Air-foil using CFD Hemanth Kotaru, B.Tech (Civil Engineering)
Grid Dependence Study of Transonic/Supersonic Flow Past NACA Air-foil using CFD Hemanth Kotaru, B.Tech (Civil Engineering) Abstract Computational fluid dynamics is a relatively young field in engineering.
More informationNumerical Simulations of Fluid-Structure Interaction Problems using MpCCI
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI François Thirifay and Philippe Geuzaine CENAERO, Avenue Jean Mermoz 30, B-6041 Gosselies, Belgium Abstract. This paper reports
More informationTVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II
TVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II EDISSON SÁVIO DE GÓES MACIEL IEA- Aeronautical Engineering Division ITA
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 informationHigh-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA )
High-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA9550-07-0195) Sachin Premasuthan, Kui Ou, Patrice Castonguay, Lala Li, Yves Allaneau,
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 informationAn efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment
The Aeronautical Journal November 2015 Volume 119 No 1221 1451 An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment J. A. Camberos
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 informationRAPID LARGE-SCALE CARTESIAN MESHING FOR AERODYNAMIC COMPUTATIONS
RAPID LARGE-SCALE CARTESIAN MESHING FOR AERODYNAMIC COMPUTATIONS Daisuke Sasaki*, Kazuhiro Nakahashi** *Department of Aeronautics, Kanazawa Institute of Technology, **JAXA Keywords: Meshing, Cartesian
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More 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 informationA high order moving boundary treatment for compressible inviscid flows 1. Abstract
A high order moving boundary treatment for compressible inviscid flows Sirui Tan and Chi-Wang Shu Abstract We develop a high order numerical boundary condition for compressible inviscid flows involving
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 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 informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationEXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS
EXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS Edisson Sávio de Góes Maciel, edissonsavio@yahoo.com.br Mechanical
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 informationTAU mesh deformation. Thomas Gerhold
TAU mesh deformation Thomas Gerhold The parallel mesh deformation of the DLR TAU-Code Introduction Mesh deformation method & Parallelization Results & Applications Conclusion & Outlook Introduction CFD
More 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 informationNumerical Methods in Aerodynamics. Fluid Structure Interaction. Lecture 4: Fluid Structure Interaction
Fluid Structure Interaction Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark
More informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
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 informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
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 informationResearch Article A Computational Investigation of Unsteady Aerodynamics of Insect-Inspired Fixed Wing Micro Aerial Vehicle s 2D Airfoil
Advances in Aerospace Engineering, Article ID 5449, 7 pages http://dx.doi.org/1.1155/214/5449 Research Article A Computational Investigation of Unsteady Aerodynamics of Insect-Inspired Fixed Wing Micro
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 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 informationExploiting a database to predict the in-flight stability of the F-16
Exploiting a database to predict the in-flight stability of the F-16 David Amsallem and Julien Cortial December 12, 2008 1 Introduction Among the critical phenomena that have to be taken into account when
More informationComputational shock and Mach waves visualization aiding the development of aerodynamic design techniques
Computational shock and Mach waves visualization aiding the development of aerodynamic design techniques H. Sobieczky, M. Hannemann Inst. of Fluid Mechanics, DLR German Aerospace Research Establishment,
More informationOzenCloud Case Studies
OzenCloud Case Studies Case Studies, April 20, 2015 ANSYS in the Cloud Case Studies: Aerodynamics & fluttering study on an aircraft wing using fluid structure interaction 1 Powered by UberCloud http://www.theubercloud.com
More informationBest Practices Workshop: Overset Meshing
Best Practices Workshop: Overset Meshing Overview Introduction to Overset Meshes Range of Application Workflow Demonstrations and Best Practices What are Overset Meshes? Overset meshes are also known as
More informationOptimization of Laminar Wings for Pro-Green Aircrafts
Optimization of Laminar Wings for Pro-Green Aircrafts André Rafael Ferreira Matos Abstract This work falls within the scope of aerodynamic design of pro-green aircraft, where the use of wings with higher
More informationRecent applications of overset mesh technology in SC/Tetra
Recent applications of overset mesh technology in SC/Tetra NIA CFD Seminar October 6, 2014 Tomohiro Irie Software Cradle Co., Ltd. 1 Contents Introduction Software Cradle SC/Tetra Background of Demands
More informationAn Object-Oriented Serial and Parallel DSMC Simulation Package
An Object-Oriented Serial and Parallel DSMC Simulation Package Hongli Liu and Chunpei Cai Department of Mechanical and Aerospace Engineering, New Mexico State University, Las Cruces, New Mexico, 88, USA
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 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 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 informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
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 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 informationLecture 1.1 Introduction to Fluid Dynamics
Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory
More 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 informationDNV GL s 16th Technology Week
OIL & GAS DNV GL s 16th Technology Week Advanced Simulation for Offshore Application: Application of CFD for Computing VIM of Floating Structures 1 SAFER, SMARTER, GREENER OUTLINE Introduction Elements
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More information