Fluid Thermal Interaction of High Speed Compressible Viscous Flow Past Uncooled and Cooled Structures by Adaptive Mesh
|
|
- Miles Owens
- 5 years ago
- Views:
Transcription
1 The 20th Conference of echanical Engineering Network of Thailand October 2006, Nakhon Ratchasima, Thailand luid Thermal nteraction of High peed Compressible Viscous low Past Uncooled and Cooled tructures by Adaptive esh Wiroj Limtrakarn Department of echanical Engineering, aculty of Engineering, Thammasat University, Bangkok 12120, Thailand Tel: , Abstract nteraction behaviors of high speed compressible viscous flow and thermal structural response of uncooled and cooled structures are presented. The compressible viscous laminar flow behavior based on the Navier tokes equations is predicted by using an adaptive cellcentered finite-element method. The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite-element method. The mass transport element is used to predict temperature distribution of cooled structure. The finite-element formulation and computational procedure are described. The high-speed flow and structural heat phenomena of uncooled and cooled structures are studied by ach 10 flow past a flat plate to demonstrate their interaction. 1 ntroduction luid Thermal analysis methods have an important role in the design of high speed flight vehicles, such as hypersonic airbreathing vehicles [1], for predicting vehicles aerothermostructural performance. ignificant coupling occurs between high-speed flow phenomena, aerodynamic heating rates on structural surfaces, structural temperature and their gradients, as well as structural deformations and stresses, creating multidisciplinary interaction phenomena. These coupled effects indicate that the analysis of high-speed flowstructure interaction is an important consideration to highspeed vehicle design. uch coupled effects have been studied by a number of researchers recently. Computational fluid and structural dynamics commercial programs were combined together for predicting the flow and structure behaviors [2, 3]. An adaptive finite element method was proposed [4] to predict complex behavior of fluid thermal structure interaction. n the present paper, combination of mass transport and integrated flowthermal-structural analysis approach for predicting each disciplinary behavior and their interaction is presented. or high-speed compressible flows, the cell-centered finite-element method [5, 6] is combined with an adaptive meshing technique to solve the Navier-tokes equations. The Galerkin finite-element method is applied to solve the structural energy and mass transport equations for temperature distribution. The paper starts by explaining the theoretical formulation of high-speed compressible flow, structural heat transfer, and mass transport. The high speed flow, structural heat transfer, and mass transport behaviors are then studied by applying the present method to ach 10 flow past a uncooled and cooled flat plate to demonstrate their interdisciplinary coupling. 2. Theoretical ormulation and olution Procedure 2.1 Governing Equations The equations for the high-speed compressible flow, the structural heat transfer, mass transport, and the structural analysis in two dimensions are described below. High-speed compressible viscous flow The equations for high-speed compressible viscous laminar flow are represented by the conservation of mass, momentum, and energy. These equations are written in the conservation form [7] as { U } t x + { } + { } E y 0 (1) where the subscript denotes the fluid analysis. The vector { U } contains the fluid conservation variables defined by { U } ρu ρv ρε T ρ (2) where ρ is the fluid density, u and v are the velocity components in the x and y directions, respectively, and ε is the total energy. The vectors { E } and { } consist of the flux components in the x and y directions, respectively [6]. tructural heat transfer The thermal response of the structure is described by the energy equation in the conservation form as t U T + ET + T x y G T (3) where the subscript T denotes the structural heat transfer analysis. The vector U T contains the thermal conservation variable defined by U T ρ TcTTT (4)
2 where c T is the specific heat of structure. The heat flux components E T and T are and E T TT kt and T x G T is the heat source. TT kt (5) y ass Transport The thermal reponse of the cooled structure is presnted by adding the mass transport equation in the form as Γ e ({ G } + { G }) nˆ dγ (10) where the flux vectors { G } and { V} and viscous flux vectors of { E + } and { V V } V G are the inviscid E +, respectively, and nˆ is the unit vector normal to the element boundary, Γ e. Equation (10) is evaluated by summing the normal fluxes from all sides, Γ e, of the element. The fluxes normal to the element sides are then approximated by the numerical inviscid and viscous. By applying an explicit time marching algorithm [7], Eq. (10) becomes fluxes, { } G and { G V} t U + E x + x 0 (6) Ae t { Un+ 1 U } δ ({ } { GV} ) n G + (11) where the subscript denotes the mass tranport analysis. The vector U contains the mass transport conservation variable defined by U ρ ct (7) where c is the specific heat of structure. The heat flux components E T and T are E ρ cu T and 2.2 inite-element ormulation T k (8) y The cell-centered finite-element method is applied to the Navier-tokes equations to derive the finiteelement equations. The Galerkin finite element approach is applied to the structural heat transfer equation, mass transport equation, and the equilibrium equations to derive the corresponding finite-element equations. The derivation procedures are briefly described below. inite-element flow equations The method of weighted residuals [8] is applied to Eq. (1) over the element domain, Ω, by using a unit interpolation function as t { } U dω x E dω { } y dω { } The Gauss divergence theorem is then applied to the flux integral terms of Eq. (9) to yield, x { } E dω y dω + { } (9) where U n 1 + and U n are the conservation variables at the time steps n+1 and n, respectively; A e is the element area; δ is the length of the element side. t is the allowable time step following the CL and viscous stability requirement [6]. G V is the average viscous flux and G is the average inviscid flux which is given by G [ GL GR + A* (UL U )] 1 R 2 + (12) where the superscripts L and R denote the left and right elements, respectively. The last term in Eq. (12) may be viewed as an artificial diffusion needed for the solution stability. This diffusion is represented by the product of the Jacobian matrix A * and the difference between the left and right element conservation variables [4]. U L and U R By substituting Eq. (12) into Eq. (11), then Eq. (11) becomes Ae t 1 2 { Un+ 1 U n } δ [{ L} + { R} + ({ L} { R G G A* U U })] { V} δ G (13) inite-element structural heat transfer equations The method of weighted residuals is applied to Eq. (3), over the element domain, Ω, by assuming a linear distribution of the conservation variable, U T, and the flux components E T and T. The finite-element equations can then be derived in the form:
3 [ ]{ } n + U 1 T { R T} n 1 + { T} n 2 R (14) where [ ] is the mass matrix, and { U T} n + 1 { U T} n+ 1 { U T} n at time n+1. The { R T} n 1 and { R T} n 2 vectors are associated with the thermal fluxes within each element and across the element boundary, respectively, and are given by, N x { } n R T 1 N dω{ En T } Γ e N y N dω{ T n } { R T} n 2 { }( En + n Tnx T n y ) + (15) N dγ (16) shock line from the sharp leading edge. The accuracy of the shock resolution and the shock angle strongly depends on the finite element mesh near the sharp leading edge. A total of 15,090 triangular elements are generated in the inviscid region and 1,230 quadrilateral elements inside the boundary layer. Ten graded layers of quadrilateral elements are used inside the boundary layer to capture steep temperature gradients for the accurate aerodynamic heating rate prediction. The predicted flow solution is shown by the density contours in igure 2(b). 0.3 m inite-element mass transport equations The galerkin method is applied to Eq. (6) in the same fashion as in the structural heat transfer analysis. The finite element equations can then be deived the form: [ ]{ } n + C U 1 { R } n 1 + { } n 2 R (17) where [ C ] is the mass matrix, and { } n + U 1 { } n + U 1 { U } n at time n+1. The { R } n 1 and { R } n 2 vectors are associated with the thermal fluxes within each element and across the element boundary, respectively, and are given by, { R } n 1 { R } n 2 { } N N n dω T x k x Ω (18) N n N u dω U x Ω (19) 10 Computational Domain hock Wave y x Convective coolant igure 1. ach 10 flow past a flat plate. 8,865 nodes 16,320 elements 4. Application A ach number 10 flow past a flat plate is presented to evaluate performance of the adaptive cell-centered finite-element method for high-speed flow analysis. The problem statement of a mach 10 flow past a flat plate as shown in igure 1 had the flow entering through the left boundary of the computational fluid domain. The shock wave is created from the leading edge as highlighted in the figure. The inlet flow conditions consist of specifying ρ 4.303E-05 kg/m 3, u 1,418.7 m/s, v 0, ε J/kg, Re 5,580 with the wall temperature of K. The combined method of the cell-centered finite-element analysis and the adaptive meshing technique is applied to solve the problem. igure 2(a) shows the final adaptive mesh that consists of small elements clustered along the (a) inal adaptive mesh (b) Density distribution (x10-5 kg/m 3 ) igure 2. inal adaptive mesh and corresponding density contours for mach 10 viscous flow past a flat plate.
4 Aerodynamics heating rate from fluid analysis is applied along fluid structure interface surface. The bottom surface of flat plate is convectively cooled by hydrogen stream with the inlet temperature of K. The hydrogen flow rate is kg/s and the coolant film coefficient is kw/m 2. The finite element model consists of 2,360 triangular elements and 118 mass transport element representing the hydrogen coolant flow. aterials of flat plate is copper. Temperature, K x, m (a) Uncooled structure Temperature, K through the thickness. Copper and beryllium are the candidate materials. igure 3(b) compares temperature distribution of both materials along interface of cooled structure at time t 0.1 sec K and K are the maximum temperature at leading edge of copper and beryllium, respectively. Beryllium has a little higher temperature at leading edge and after x 0.02 m. beryllium has temperature 1.2 K less than copper. 5. Concluding Remarks The multidisciplinary interaction behaviors of highspeed compressible flow, structural heat transfer, mass transport, and structural response were presented using the adaptive finite-element method. The finite-element method based on the cell-centered algorithm was used to predict the high-speed compressible flow behavior. The method was combined with the adaptive meshing technique to improve the flow accuracy. The ach 10 flow past a flat plate was then used to study the flowstructure interaction and to evaluate the performance of the proposed analysis procedure. The example highlights the interaction behavior between the high-speed flow, mass transport, and the thermal-structural response of uncooled and cooled structure. The high-speed compressible viscous flow behavior is well predicted by adaptive finite element method and the fluid structure interaction behavior by the proposed finite-element coupling procedure. 26 Copper Beryllium x, m Acknowledgments This research is partially supported by the Thailand Research und (TR) for the enior cholar Professor Pramote Dechaumphai. (b) Cooled tructure igure 3. Temperature distribution along x distance at t 0.1 sec. igure 3 (a) shows the temperature distribution at time t 0.1 sec. along fluid structure interface based on uncooled structure with insulated boundary conditions along left, right and bottom surface. The maximum temperature is K on the leading edge. lat plate materials with high thermal conductivity are used to reduce the leading edge temperature, temperature gradient References [1] Glass, D. E., erski, N. R., Glass, C. E., Airframe Research and Technology for Hypersonic Airbreathing Vehicles, AAA Paper [2] Baum, J. D., Development of a Coupled CD/CD ethodology using an Embedded CD Approach, Conference on Computational Physics. [3] Lohner, R. Baum, J. D. estreau, E. harov, D. Charman, C. Pelessone, D., Adaptive Embedded Unstructured Grid ethods, AAA Paper [4] Limtrakarn, W., and Dechaumphai, P., Computations of High peed Compressible lows
5 with Adaptive Cell Centered inite Element ethod, J Chin nst Eng, Vol. 26, pp [5] Gnoffo, P. A., Application of Program LAURA to Three-dimensional AOTV lowfields, AAA Paper [6] Dechaumphai, P., and Limtrakarn, W., Adaptive Cell-Centered inite Element Technique for Compressible lows, Journal of Energy, Heat and ass Transfer, Vol. 21, pp [7] Hirsch, C., Numerical Computation of nternal and External lows. Vol. 1, New York: Wiley. [8] Zienkiewicz, O. C., and Taylor, R. L., The inite Element ethod. ifth Ed., Woburn: Butterworth-Heinemann.
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationParabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
More informationNumerical Investigation of Shock wave Turbulent Boundary Layer Interaction over a 2D Compression Ramp
Advances in Aerospace Science and Applications. ISSN 2277-3223 Volume 4, Number 1 (2014), pp. 25-32 Research India Publications http://www.ripublication.com/aasa.htm Numerical Investigation of Shock wave
More informationThermal Damage Improvement of Slider Fixture in Ion Mill Heatless Etching Process by FEM
Thermal Damage Improvement of Slider Fixture in Ion Mill Heatless Etching Process by FEM Wiroj Limtrakarn Department of Mechanical Engineering, Faculty of Engineering, Thammasat University, Pathum Thani
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationOutlines. simple relations of fluid dynamics Boundary layer analysis. Important for basic understanding of convection heat transfer
Forced Convection Outlines To examine the methods of calculating convection heat transfer (particularly, the ways of predicting the value of convection heat transfer coefficient, h) Convection heat transfer
More informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationCFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM
More informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More information2 Navier-Stokes Equations
1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationNumerical investigation on the effect of inlet conditions on the oblique shock system in a high-speed wind tunnel
American Journal of Physics and Applications 2013; 1(3): 91-98 Published online November 30, 2013 (http:// www.sciencepublishinggroup.com/j/ajpa) doi: 10.11648/j.ajpa.20130103.16 Numerical investigation
More informationMODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI
TASK QUARTERLY 12 No 3, 173 184 MODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI AND SŁAWOMIR DYKAS Institute of Power Engineering and
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationTHERMAL PROPERTIES OF THERMOBLOCK 9LB25-O
SCHOOL OF MINING AND METALLURGICAL ENG. LABORATORY OF METALLURGY Technical Report ΤΕ 07/2018 THERMAL PROPERTIES OF THERMOBLOCK 9LB25-O : MAY 2018 Table of Contents 1. INTRODUCTION...3 2. METHODOLOGY...3
More informationInvestigation of Transpiration Cooling Effectiveness for Air- Breathing Hypersonic Vehicles
17th AIAA International Space Planes and Hypersonic Systems and Technologies Conference 11-14 April 2011, San Francisco, California AIAA 2011-2253 Investigation of Transpiration Cooling Effectiveness for
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER
THERMAL SCIENCE: Year 2018, Vol. 22, No. 2, pp. 963-972 963 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER by Jitesh RANA, Anshuman SILORI, Rajesh MAITHANI *, and
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationDevelopment of Two-Dimensional Convergent-Divergent Nozzle Performance Rapid Analysis Project
International Forum on Energy, Environment Science and Materials (IFEESM 015) Development of Two-Dimensional Convergent-Divergent Nozzle Performance Rapid Analysis Project Yaxiong YANG 1,a *, Eri Qitai,b,
More informationTransport equation cavitation models in an unstructured flow solver. Kilian Claramunt, Charles Hirsch
Transport equation cavitation models in an unstructured flow solver Kilian Claramunt, Charles Hirsch SHF Conference on hydraulic machines and cavitation / air in water pipes June 5-6, 2013, Grenoble, France
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationSolving the Navier-Stokes Equations
FMIA F. Moukalled L. Mangani M. Darwish An Advanced Introduction with OpenFOAM and Matlab This textbook explores both the theoretical oundation o the Finite Volume Method (FVM) and its applications in
More information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationPHYSICAL MECHANISM OF CONVECTION
Tue 8:54:24 AM Slide Nr. 0 of 33 Slides PHYSICAL MECHANISM OF CONVECTION Heat transfer through a fluid is by convection in the presence of bulk fluid motion and by conduction in the absence of it. Chapter
More informationAnalysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data
45th AIAA Aerospace Science Meeting and Exhibit, January 8 11, 25/Reno, Nevada Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data M. Wu and M.P. Martin Mechanical
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 3B: Conservation of Mass C3B: Conservation of Mass 1 3.2 Governing Equations There are two basic types of governing equations that we will encounter in this course Differential
More informationCALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE
J. Comput. Fluids Eng. Vol.17, No.4, pp.69-74, 2012. 12 / 69 CALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE IN NONEQUILIBRIUM HYPERSONIC FLOW M. Ahn Furudate * Dept. of Mechatronics Engineering,
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationCFD ANALYSIS OF TRIANGULAR ABSORBER TUBE OF A SOLAR FLAT PLATE COLLECTOR
Int. J. Mech. Eng. & Rob. Res. 2013 Basavanna S and K S Shashishekar, 2013 Research Paper ISSN 2278 0149 www.imerr.com Vol. 2, No. 1, January 2013 2013 IJMERR. All Rights Reserved CFD ANALYSIS OF TRIANGULAR
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationA Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
More informationPROBLEM Node 5: ( ) ( ) ( ) ( )
PROBLEM 4.78 KNOWN: Nodal network and boundary conditions for a water-cooled cold plate. FIND: (a) Steady-state temperature distribution for prescribed conditions, (b) Means by which operation may be extended
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationChamber Dynamics and Clearing Code Development Effort
Chamber Dynamics and Clearing Code Development Effort A. R. Raffray, F. Najmabadi, Z. Dragojlovic University of California, San Diego A. Hassanein Argonne National Laboratory P. Sharpe INEEL HAPL Review
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationME 144: Heat Transfer Introduction to Convection. J. M. Meyers
ME 144: Heat Transfer Introduction to Convection Introductory Remarks Convection heat transfer differs from diffusion heat transfer in that a bulk fluid motion is present which augments the overall heat
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationVariable Normalization (nondimensionalization and scaling) for Navier-Stokes equations: a practical guide
TR-CFD-13-77 Variable Normalization (nondimensionalization and scaling) for Navier-Stokes equations: a practical guide Marc Montagnac Technical Report V. 1.0 October 013 Copyrighted by the author(s) Centre
More informationIncompressible Flow Solver by Means of Pseudo-Compressibility Method
Incompressible Flow Solver by Means of Pseudo-Compressibility Method Árpád Veress Department of Aircraft and Ships Budapest University of Technology and Economics H Budapest, Sztoczek u. 6. J. ép., Hungary
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationNumerical Simulation of Hypersonic Boundary-Layer Instability in a Real Gas with Two-Dimensional Surface Roughness
AIAA Aviation 22-26 June 2015, Dallas, TX 45th AIAA Fluid Dynamics Conference AIAA 2015-3077 Numerical Simulation of Hypersonic Boundary-Layer Instability in a Real Gas with Two-Dimensional Surface Roughness
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationCFD Simulation of Internal Flowfield of Dual-mode Scramjet
CFD Simulation of Internal Flowfield of Dual-mode Scramjet C. Butcher, K. Yu Department of Aerospace Engineering, University of Maryland, College Park, MD, USA Abstract: The internal flowfield of a hypersonic
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More information4.2 Concepts of the Boundary Layer Theory
Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationFLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.
FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.Bhosale 1 Research scholar, 2 Head of department & Asst professor Department
More informationComputational Analysis of Scramjet Inlet
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology Volume 3, Special Issue 3, March 2014 2014 IEEE International Conference
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationINTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA)
INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) Title 3 Column (full page) 2 Column What is Finite Element Analysis? 1 Column Half page The Finite Element Method The Finite Element Method (FEM)
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationHypersonic flow and flight
University of Stuttgart, Aerospace Engineering and Geodesy Dept. - Lecture - Hypersonic flow and flight Master Level, Specialization 4 lecture hours per week in WS, 3-6 LPs/ECTS Lecturer: Dr. Markus J.
More informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationMIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER
THERMAL SCIENCE, Year 2016, Vol. 20, No. 2, pp. 483-492 483 MIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER by Liping WEI, Youjun LU*, and Jinjia WEI State Key Laboratory of Multiphase
More informationTAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF
TAU Extensions for High Enthalpy Flows Sebastian Karl AS-RF Contents Motivation Extensions available in the current release: Numerical schemes for super- and hypersonic flow fields Models for gas mixtures,
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationDetached Eddy Simulation on Hypersonic Base Flow Structure of Reentry-F Vehicle
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (2014) 000 000 www.elsevier.com/locate/procedia APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology,
More informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
More informationSimplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis
1 Portál pre odborné publikovanie ISSN 1338-0087 Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis Jakubec Jakub Elektrotechnika 13.02.2013 This work deals with thermo-hydraulic processes
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationProblem 4.3. Problem 4.4
Problem 4.3 Problem 4.4 Problem 4.5 Problem 4.6 Problem 4.7 This is forced convection flow over a streamlined body. Viscous (velocity) boundary layer approximations can be made if the Reynolds number Re
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationExperimental and numerical investigations of heat transfer and thermal efficiency of an infrared gas stove
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Experimental and numerical investigations of heat transfer and thermal efficiency of an infrared gas stove To cite this article:
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationCOMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
Suranaree J. Sci. Technol. Vol. 20 No. 4; October - December 2013 257 COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationEffect of Periodic Variation of Sol-air Temperature on the Performance of Integrated Solar Collector Storage System
Engineering, 2010, 2, 832-840 doi:10.4236/eng.2010.210106 Published Online October 2010 (http://www.scirp.org/journal/eng) Effect of Periodic Variation of Sol-air Temperature on the Performance of Integrated
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationCavitation modeling using compressible Navier Stokes and Korteweg equations
Computational Methods in Multiphase Flow VIII 425 Cavitation modeling using compressible Navier Stokes and Korteweg equations T. Ait Ali 1, S. Khelladi 1, L. Ramirez 2 & X. Nogueira 2 1 DynFluid Laboratory,
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More information