TRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE.
|
|
- Marybeth Atkinson
- 6 years ago
- Views:
Transcription
1 TRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE Chiara Manletti Space Launcher Systems Analysis (SART), DLR, Cologne, Germany, chiara.manletti@dlr.de ABSTRACT It is renowned that transient simulation is a powerul tool which reveals undamental characteristics o any liquid engine. Using the program TLR-Sim (Transient Liquid Rocket Engine Simulator), whose development has recently begun, both LOX/Kerosene and LOX/Methane engine components have been closer examined whilst seeking or a clear indication o advantages associated with the use o one or the other in a number o potential engine applications. Results o the perormed simulations are presented and commented. The program TLR-Sim and its development are shortly presented. A short summary o the program structure and main characteristics precedes the detailed engine analysis and comparison. A brie conclusion will ollow, outlining both the main results o the inspection and a uture roadmap or a subsequent program evolution. NOMENCLATURE a Speed o sound m/s d (hydraulic) diameter m E Total energy per unit volume J/m 3 F Cross-sectional Area m h Surace roughness m h Surace roughness coeicient - l length m m Mass low kg/s n Factor o margin - p pressure bar κ Thermal conductivity W/mK Re Reynolds number - T temperature K u velocity m/s V volume m³ Finite dierence n.a. µ Viscosity Pa s ρ density kg/m 3 ϖ Frequency Hz ξ Friction actor /m 4. INTRODUCTION In the rame o an internal study at DLR, the program TLR-Sim is under development with the goal o perorming transient simulations o liquid rocket engine components. In the attempt o inding the best alternative to the cryogenic propellant couple LOX/LH, studies have been concentrated on comparisons between numerous propellant couples. Amongst the most targeted are LOX/Kerosene, or the latter's storability characteristics, and LOX/Methane, or possible lower acquisition costs, potentially higher achievable speciic impulse, lower pressure drop in cooling channels, superior cooling properties and higher coking limits []. Engine nominal operational point comparisons have already been examined and presented by the DLR Space Launcher Systems Analysis (SART) department. These have led to the conclusion that although the increased motor mass, booster size and hence drag, a reusable methane engine could provide a competitive solution to the extensively studied kerosene engines, with the premise o urther studies. The purpose o the study underway is to ill one o the gaps existing in previous studies, namely the examination and comparison o the transient behaviours o the propellants themselves under speciic conditions.. PROGRAM STRUCTURE The program TLR-Sim possesses a highly modular structure in order to assure as high a degree o lexibility as possible. The general program structure is illustrated in the scheme in Figure.
2 major techniques in the category o more detailed methods are the ollowing: - Method o characteristics (MOC) - Riemann Solvers (RS) The Method o Characteristics has a wide range o applications, thus an equally vast class o possible constructions. One can choose, or example, between explicit and implicit solving methods. A literature survey has revealed that implicit methods, when applied to sharp transients, can lead to initial errors which are carried through the entire calculation and may lead to questionable results. The choice thus alls on the explicit solving routines. Further decisions, coherent with the problem at hand, must be made in the uture. Fig.. TLR-Sim General Flow Diagram. TLR-Sim has been programmed to autonomously generate the components order o computation subsequently to having processed the user-deined input. The simulation is then perormed in three clearly deined steps. The irst part o this step-wise approach consists in calculating components that have been deined as elements. These include pumps, valves, combustion devices and the like. Subsequently, ducts, pipes, and lexibles are computed with the newly acquired boundary conditions stemming rom their neighbouring elements. Only ater this second step, does a complete engine run-through take place. 3. DUCT SIMULATION METHODS Within TLR-Sim, numerical routines or the simulation o a number o dierent components have already been implemented. Being, however, the ocus o this paper luid low through pipe lines, this section concentrates only on the simulation methods existing within TLR- Sim or the computation o duct lows. 3. Ducts All components alling under this category, that is, rigid and lexible pipes can be simulated by applying a number o dierent techniques. These techniques may be classiied according to the detail to which they examine the low. General ODEs (ordinary dierential equations) can be used to compute pressure changes, mass low variations etc. (see Section 5) More detailed schemes allow, not only the computation o the low properties, but also the evaluation o its structure. Two Riemann Solvers solve or the Riemann problem (RP) which is none other than a special case o the Initial Value Problem (IVP). Riemann solvers can be adopted to solve a local RP at any position along a duct. The solution to the RP is then used to solve or the Euler equations with or without additional source terms. Although being somewhat more complex, the use o Riemann solvers can be implemented in the solving o problems where the MOC is not suitable. As will be described in Section 6, the study here presented has made use o the Riemann solvers to solve the Navier- Stokes Equations. 4. KEROSENE METHANE COMPARISON 4. Analysis Background In previous studies conduced in the Launcher s Systems Analysis department o the German Aerospace Centre, DLR, two engines in the RD-80 class have been hypothesised to compare the perormance o methane and kerosene under nominal engine operation. The paper at hand has taken these engines as its starting point []. Table summarises the tank conditions or each o the methane and kerosene RD-80 class engines. Table.Propellant Tank Parameters Methane Kerosene RP Tank Pressure bar 3 Tank Temperature K 0 90 Density kg/m
3 4. Analysis Breakdown The study here presented is based on a two-old analysis. The irst is based on non-linear ODEs (ordinary dierential equations) describing the evolution o pressure and mass low along a duct and is perormed with an ODE integrating tool. The second aspect o the analysis is a more detailed one consisting in the implementation o the Navier-Stokes equations to evaluate the potential dierences between kerosene and methane. Viscosity [ kg/ms ] 4,0e-03 3,0e-03,0e-03,0e-03 Jet A Viscosity vs. Temp. (along cnst. p curves) Jet A 0. MPa Jet A MPa Jet A MPa Jet A 3 MPa JetA 4 MPa In the ollowing sections each o the above-mentioned assessments is introduced, the underlying model is highlighted and the out-coming major results are presented. 4.3 Propellant Properties It should be noted that Jet A data has been used due to its easy availability. Dierence between RP and Jet A due exist but or present purposes are o minor importance [4], [6]. Figures and 3 illustrate the viscosity o methane and kerosene as unctions o temperature and pressure. The plots show that, within the pressure and temperature range o interest, methane viscosity values can be an order o magnitude smaller than those o kerosene. Similar curves depicting methane and kerosene density show in this case the ratio is only two-old. Viscosity [ Pa*s ],5e-04,0e-04,5e-04,0e-04 5,0e-05 Methane Viscosity vs. Temp. (along cnst. p curves) 0.0 MPa 0. MPa 0.5 MPa 0.5 MPa MPa 5 MPa 0 MPa 0 MPa 0,0e Temperature [ K ] Fig.. Methane Viscosity Plot as a ( T, p ) Fn. 0,0e Temperature [ K ] Fig. 3. Kerosene Viscosity Plot as a ( T, p ) 5. ODE ANALYSIS 5. Introduction Fn. The ODE analysis has been perormed with a program widely used at DLR or the simulation control laws such as those existing in aviation light control systems. The tool has the capability o integrating ordinary dierential equations and allows or eedback loops. 5. ODE Equations The luid parameters, pressure and mass low, characterising the luid low though a duct can be expressed according to the ollowing ODEs [3]: dm dp dt in dt F = l a = V p in p out [ m in m out ] ρ ξ m out where the subscript stands or { M } kerosene and methane respectively. () () K, or The riction actor ξ is a unction o the Reynolds number, and or turbulent lows, i.e. Re 300 is deined as: 4 68 l ξg ξ = 0. h + + (3) Re d F F h where h =. d
4 5.3 ODE Model The model used in the ODE analysis is aimed at reproducing, in more simple terms, the start sequence o a rocket engine employing a burst disk and includes a tank directly connected with a duct. Bends and areavariations have not been taken into account. Acceleration eects stemming rom the launcher system have also been neglected. In the model the burst disk is assumed to be located just beore the irst duct section. The duct is broken down into sections each having an inlet pressure, an outlet pressure and inlet and outlet mass lows. Figure 4 illustrates the model used. Each duct section is o equal length and is user deined according to: πa l max (4) ϖ n max where ϖ max is the maximum requency that is to be examined. Under the spectrum o the ODE analysis, three models have been used. Each model diers in the length o its duct sections. Table summarises the main details o each model coniguration. Table. Model Conigurations and Corresponding Frequencies Model Section Length [m] ϖ [Hz] max Model Model Model ODE Results Fig. 4. ODE Model The graphs shown below have been obtained by implementing the three dierent models described above. Although each model diers rom the others in the length o its duct sections the total duct length examined has been kept. Each graph depicts normalised pressure, i.e. the pressure in the duct normalised to the pressure in the respective tanks (3 bar or methane and bar or kerosene), against time in seconds. In addition to centre o the duct where the evolution o pressure is shown, two urther locations have been plotted or each model. These are the entrance and exit locations o the duct. The initial pressure within the duct has been assumed to be one bar lower than pressure existing within the tank (corresponding to a normalised pressure o 0.5 or kerosene and approximately 0.67 or methane)..0e+00 Kerosene & Methane Normalised Pressure vs. Time - Model.00E E E E-0.00E E Fig. 5. ODE Model Test Results Methane Duct Entrance Methane Duct Exit Kerosene Duct Entrance Kerosene Duct Exit
5 .0E+00 Kerosene & Methane Normalised Pressure vs. Time - Model.00E E E-0 Methane Duct Entrance 4.00E-0 Methane Duct Exit Kerosene Duct Entrance.00E-0 Kerosene Duct Exit 0.00E Fig. 6. ODE Model Test Results.0E+00 Kerosene & Methane Normalised Pressure vs. Time - Model 3.00E E E E-0.00E E+00 Methane Duct Entrance Methane Duct Exit Kerosene Duct Entrance Kerosene Duct Exit Fig. 7. ODE Model 3 Test Results 6. -DIMENSIONAL NAVIER-STOKES ANALYSIS 6. Introduction The one-dimensional Navier-Stokes analysis has been built on the Euler equations with the addition o a series o sources terms deemed as important or luid low in ducts. These source terms include area variations as well as viscous and heat contributions. 6. Navier-Stokes Equations and Riemann Solvers The Navier-Stokes equations in conservative orm are as ollows: U t ρ ρu E t + + u F( U ) ρu x ( ρu + p ) ( E + p ) x = S + S = S + S (5) where S and S are source terms accounting or areavariation and viscous and heat eects and which upgrade the one-dimensional Euler equations to the one-dimensional Navier-Stokes equations with areavariations.
6 By adopting the WAF (Weight Averaged Flux) method proposed by Toro [5], and using an exact Riemann solver and an appropriate lux limiter, the Riemann problem can be solved or each U ( i, i + ), along the duct. The result to the local Riemann problem in conjunction with a conventional TVD (Total Variation Diminishing) scheme can then be implemented to obtain the result to the Euler equations or all positions, i along the duct or a given time, n. Subsequently the source term eects are added to obtain the solution to the Navier-Stokes equations. 6.3 Results Figure 8, as Figures 5 through 7 or the ODE analysis, depicts normalised pressure vs. time. The results correspond to a shock-tube problem set up to have the same starting conditions as the ODE analysis. The purpose o this evaluation is to provide an indirect comparison to the ODE results or a somewhat dierent problem and thus validate the general characteristics shown by kerosene and methane. 9.00E E-0 Normalised Pressure vs. Time (Navier Stokes) 8.00E E E E E E E Fig. 8. Navier-Stokes Results 7. CONCLUSIONS AND OUTLOOK The two-old analysis here presented has highlighted three undamental aspects. Firstly: the dynamic behaviours o kerosene and methane have been shown to possess substantial dierences, thus underlying the importance o a transient analysis. Secondly: the ODE and Navier-Stokes analyses have led to coherent results though applied to dierent test cases. This shows that the transient analysis can be perormed with tools having dierent ocuses, thus allowing the choice o the best tool to suit the problem at hand and thus holding a considerable advantage. Thirdly: a more advanced transient simulation which includes additional engine components and examines their interactions or a range o operating conditions represents an advantage in engine analysis and should be pursued. Thus, urther examinations and program developments are underway not only or dierent ambient conditions but or more complex models, comprehending other LRE components such as valves and combustion devices. 8. REFERENCES [] Anderson, J. D.: Fundamentals o Aerodynamics, Mc-Graw Hill, Singapore, 99. [] Burkhardt, H., et al.: Kerosene vs Methane: A Propellant Tradeo or Reusable Liquid Booster Stages, to appear in Journal o Spacecrat and Rockets. [3] Belyaev, E.N., et al.: Mathematical Simulation o the Operating Processes o Liquid Rocket Engines, MAI, Moscow, 999. [4] Rachner, M.: Die Stoeigenschaten von Kerosin Jet A-, Deutsches Zentrum uer Lut- und Raumahrt e.v., Koeln, March 998. [5] Toro, E.F.: Riemann Solvers and Numerical Methods or Fluid Dynamics, A Practical Introduction, Springer-Verlag, Berlin Heidelberg, 997. [6] N.N.: NIST Chemistry WebBook, [Online March 004].
When water (fluid) flows in a pipe, for example from point A to point B, pressure drop will occur due to the energy losses (major and minor losses).
PRESSURE DROP AND OSSES IN PIPE When water (luid) lows in a pipe, or example rom point A to point B, pressure drop will occur due to the energy losses (major and minor losses). A B Bernoulli equation:
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationFiltration. Praktikum Mechanical Engineering. Spring semester ML F16 Tel.:
Praktikum Mechanical Engineering Spring semester 2018 Filtration Supervisor: Davide Stucchi ML F16 stucchid@ptl.mavt.ethz.ch Tel.: 044 632 25 05 1 1 Table o Contents 1 TABLE OF CONTENTS... 2 2 INTRODUCTION...
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationFLOW CHARACTERISTICS OF HFC-134a IN AN ADIABATIC HELICAL CAPILLARY TUBE
E HEFAT7 5 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics Sun City, South Arica Paper number: KM1 FLOW CHARACTERISTICS OF HFC-1a IN AN ADIABATIC HELICAL CAPILLARY TUBE Khan
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationKuldeep Rawat*, Ayushman Srivastav* *Assistant Professor, Shivalik College of Engineering, Dehradun.
International Journal o Scientiic & Engineering search, Volume 7, Issue 12, December-16 348 ISSN 2229-18 NUMERICAL INVESTIGATION OF HEAT TRANSFER ENHANCEMENT OVER RECTANGULAR PERFORATED FIN Abstract Kuldeep
More informationChapter 3 Water Flow in Pipes
The Islamic University o Gaza Faculty o Engineering Civil Engineering Department Hydraulics - ECI 33 Chapter 3 Water Flow in Pipes 3. Description o A Pipe Flow Water pipes in our homes and the distribution
More informationAvailable online at ScienceDirect. Energy Procedia 83 (2015 ) Václav Dvo ák a *, Tomáš Vít a
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 83 (205 ) 34 349 7th International Conerence on Sustainability in Energy and Buildings Numerical investigation o counter low plate
More informationScienceDirect. Heat transfer and fluid transport of supercritical CO 2 in enhanced geothermal system with local thermal non-equilibrium model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 7644 7650 GHGT-12 Heat transer and luid transport o supercritical CO 2 in enhanced geothermal system with local thermal
More informationLocal Heat Transfer Coefficient Measurements, Using a Transient Imaging Method With an Inverse Scheme
Local Heat Transer Coeicient Measurements, Using a Transient Imaging Method With an Inverse Scheme A. EL ABBADI, D. BOUGEARD, B.BAUDOIN Ecole des Mines de Douai, Département Energétique Industrielle, 941,
More informationAn Improved Expression for a Classical Type of Explicit Approximation of the Colebrook White Equation with Only One Internal Iteration
International Journal o Hydraulic Engineering 06, 5(): 9-3 DOI: 0.593/j.ijhe.06050.03 An Improved Expression or a Classical Type o Explicit Approximation o the Colebrook White Equation with Only One Internal
More informationANALYSIS OF ELECTRO-RHEOLOGICAL FLUID IN HYDRAULIC SYSTEM
ANALYSIS OF ELECTRO-RHEOLOGICAL FLUID IN HYDRAULIC SYSTEM Földi László, Jánosi László Institute or Mechanical Engineering Technology, Faculty o Mechanical Engineering, Szent István University, Gödöllő,
More informationFLUID MECHANICS. Lecture 7 Exact solutions
FLID MECHANICS Lecture 7 Eact solutions 1 Scope o Lecture To present solutions or a ew representative laminar boundary layers where the boundary conditions enable eact analytical solutions to be obtained.
More informationTransport Properties: Momentum Transport, Viscosity
Transport Properties: Momentum Transport, Viscosity 13th February 2011 1 Introduction Much as mass(material) is transported within luids (gases and liquids), linear momentum is also associated with transport,
More informationChapter 8 Flow in Pipes. Piping Systems and Pump Selection
Piping Systems and Pump Selection 8-6C For a piping system that involves two pipes o dierent diameters (but o identical length, material, and roughness connected in series, (a the low rate through both
More informationCOMPUTATIONAL STUDY OF CHEMICALLY REACTING HYPERSONIC FLOW
COMPUTATIONAL STUDY OF CHEMICALLY REACTING HYPERSONIC FLOW Yoshiuru Funahashi Department o Aeronautics and Astronautics, Graduate School o Engineering, The University o Tokyo, Tokyo, JAPAN Keywords: Hypersonic
More informationFLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 1
FLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 1 1. A pipe 100 mm bore diameter carries oil of density 900 kg/m3 at a rate of 4 kg/s. The pipe reduces
More informationPiping Systems and Flow Analysis (Chapter 3)
Piping Systems and Flow Analysis (Chapter 3) 2 Learning Outcomes (Chapter 3) Losses in Piping Systems Major losses Minor losses Pipe Networks Pipes in series Pipes in parallel Manifolds and Distribution
More informationKeywords Perforated pinned heat sinks, Conjugate heat transfer, Electronic component cooling.
Eect o Dierent Perorations Shapes on the Thermal-hydraulic Perormance o Perorated Pinned Heat Sinks Amer Al-Damook 1,, J.L. Summers 1, N. Kapur 1, H. Thompson 1 mnajs@leeds.ac.uk, j.l.summers@leeds.ac.uk,
More informationNUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID
NUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID Saleh Etaig*, Etaig.Mahmoud@Northumbria.ac.uk Reaz Hasan, Reaz.Hasan@Northumria.ac.uk
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More informationEXPERIMENTAL STUDY AND NUMERICAL SIMULATION OF PREFORM INFRARED RADIATIVE HEATING
EXPERIMENTAL STUDY AND NUMERICAL SIMULATION OF PREFORM INFRARED RADIATIVE HEATING Serge Monteix (*), F. Schmidt (*), Y. LeMaoult (*), G. Denis (**), M. Vigny (**) (*) Ecole des Mines d Albi-Carmaux Campus
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationDesign And Analysis Of Thrust Chamber Of A Cryogenic Rocket Engine S. Senthilkumar 1, Dr. P. Maniiarasan 2,Christy Oomman Jacob 2, T.
Design And Analysis Of Thrust Chamber Of A Cryogenic Rocket Engine S. Senthilkumar 1, Dr. P. Maniiarasan 2,Christy Oomman Jacob 2, T. Vinitha 2 1 Research Scholar, Department of Mechanical Engineering,
More informationFin efficiency of the newly developed Compartmented Coil of a Single Coil Twin Fan System
Fin eiciency o the newly developed Compartmented Coil o a Single Coil Twin Fan System ABSTRACT In predicting the perormance o any cooling coil, HVAC designers ace multiold challenges in designing the system
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
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 informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More informationNumerical Simulation of Supersonic Expansion in Conical and Contour Nozzle
Numerical Simulation of Supersonic Expansion in Conical and Contour Nozzle Madhu B P (1), Vijaya Raghu B (2) 1 M.Tech Scholars, Mechanical Engineering, Maharaja Institute of Technology, Mysore 2 Professor,
More informationHydraulic validation of the LHC cold mass heat exchanger tube.
Hydraulic validation o te LHC cold mass eat excanger tube. LHC Project Note 155 1998-07-22 (pilippe.provenaz@cern.c) Pilippe PROVENAZ / LHC-ACR Division Summary Te knowledge o te elium mass low vs. te
More informationDEVELOPMENT OF A COMPRESSED CARBON DIOXIDE PROPULSION UNIT FOR NEAR-TERM MARS SURFACE APPLICATIONS
DEVELOPMENT OF A COMPRESSED CARBON DIOXIDE PROPULSION UNIT FOR NEAR-TERM MARS SURFACE APPLICATIONS Erin Blass Old Dominion University Advisor: Dr. Robert Ash Abstract This work has focused on the development
More informationHeat-fluid Coupling Simulation of Wet Friction Clutch
3rd International Conerence on Mechatronics, Robotics and Automation (ICMRA 2015) Heat-luid Coupling Simulation o Wet Friction Clutch Tengjiao Lin 1,a *, Qing Wang 1, b, Quancheng Peng 1,c and Yan Xie
More informationCRITICAL MASS FLOW RATE THROUGH CAPILLARY TUBES
Proceedings o the ASME 010 rd Joint US-European Fluids Engineering Summer Meeting and 8th International Conerence FESM-ICNMM010 August 1-5, 010, Montreal, Canada Proceedings o ASME 010 rd Joint US-European
More informationBuoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center
July 4-6 2012 London U.K. Buoyancy Driven Heat Transer o Water-Based CuO Nanoluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center Ahmet Cihan Kamil Kahveci and Çiğdem Susantez
More informationCFD MODELLING AND VALIDATION OF HEAD LOSSES IN PIPE BIFURCATIONS
CFD MODELLING AND VALIDATION OF HEAD LOSSES IN PIPE BIFURCATIONS Kasturi Sukhapure* a, Alan Burns a, Tariq Mahmud a, Jake Spooner b. a School of Chemical and Process Engineering, University of Leeds, Leeds
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More information9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions
9. Graphing Functions by Plotting Points, The Domain and Range o Functions Now that we have a basic idea o what unctions are and how to deal with them, we would like to start talking about the graph o
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationCOMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION
HEFAT014 10 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 14 6 July 014 Orlando, Florida COMPARISON OF THERMA CHARACTERISTICS BETWEEN THE PATE-FIN AND PIN-FIN HEAT SINKS
More informationTFY4102 Exam Fall 2015
FY40 Eam Fall 05 Short answer (4 points each) ) Bernoulli's equation relating luid low and pressure is based on a) conservation o momentum b) conservation o energy c) conservation o mass along the low
More informationChannel Structure Influence on the Thermal-Hydraulic Performance of. Zigzag PCHE
The 6th International Supercritical CO2 Power Cycles Symposium March 27-29, 218, Pittsburgh, Pennsylvania Channel Structure Inluence on the Thermal-Hydraulic Perormance o Zigzag PCHE Yichao Gao Wenkai
More informationGas-side mass transfer coefficient of a laboratory column equipped with one sieve tray
Gas-side mass transer coeicient o a laoratory column equipped with one sieve tray Zhivko Ivanov, Zhelcho Steanov, Bogdan Bogdanov Astract: The inluence o plate geometry on the characteristics o luid low
More informationRobust Residual Selection for Fault Detection
Robust Residual Selection or Fault Detection Hamed Khorasgani*, Daniel E Jung**, Gautam Biswas*, Erik Frisk**, and Mattias Krysander** Abstract A number o residual generation methods have been developed
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationA NUMERICAL STUDY OF SINGLE-PHASE FORCED CONVECTIVE HEAT TRANSFER WITH FLOW FRICTION IN ROUND TUBE HEAT EXCHANGERS
www.arpapress.com/volumes/vol6issue4/ijrras_6_4_05.pd A NUMERICAL STUDY OF SINGLE-PHASE FORCED CONVECTIVE HEAT TRANSFER WITH FLOW FRICTION IN ROUND TUBE HEAT EXCHANGERS Pedram Mohajeri Khameneh 1,*, Iraj
More informationInvestigation of Initial Fouling Rates of Calcium Sulfate Solutions under Non-Boiling Conditions (Work-in-Progress)
eereed Proceedings Heat Exchanger Fouling and Cleaning: Fundamentals and Applications Engineering Conerences International Year 23 Investigation o Initial Fouling ates o Calcium Sulate Solutions under
More informationDescription of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator
This is a 1D model o an active magnetic regenerative rerigerator (AMRR) that was developed in MATLAB. The model uses cycle inputs such as the luid mass low and magnetic ield proiles, luid and regenerator
More informationROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)
AE301 Aerodynamics I UNIT D: Applied Aerodynamics ROAD MAP... D-1: Aerodynamics o 3-D Wings D-2: Boundary Layer and Viscous Eects D-3: XFLR (Aerodynamics Analysis Tool) AE301 Aerodynamics I : List o Subjects
More informationFigure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m
1. For the manometer shown in figure 1, if the absolute pressure at point A is 1.013 10 5 Pa, the absolute pressure at point B is (ρ water =10 3 kg/m 3, ρ Hg =13.56 10 3 kg/m 3, ρ oil = 800kg/m 3 ): (a)
More informationA Semi-Analytical Solution for a Porous Channel Flow of a Non-Newtonian Fluid
Journal o Applied Fluid Mechanics, Vol. 9, No. 6, pp. 77-76, 6. Available online at www.jamonline.net, ISSN 735-357, EISSN 735-3645. A Semi-Analytical Solution or a Porous Channel Flow o a Non-Newtonian
More informationCHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION
CHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION OBJECTIVE The objective of the experiment is to compare the heat transfer characteristics of free and forced convection.
More informationLEAKLESS COOLING SYSTEM V.2 PRESSURE DROP CALCULATIONS AND ASSUMPTIONS
CH-1211 Geneva 23 Switzerland EDMS No. ST/CV - Cooling of Electronics & Detectors GUIDE LEAKLESS COOLING SYSTEM V.2 PRESSURE DROP CALCULATIONS AND ASSUMPTIONS Objectives Guide to Leakless Cooling System
More information3D Numerical Modelling of Convective Heat Transfer through Two-sided Vertical Channel Symmetrically Filled with Metal Foams
P Periodica Polytechnica Mechanical Engineering P 60(4), pp. 193-202, 2016 DOI: 10.3311/PPme.8511 Creative Commons Attribution b 3D Numerical Modelling o Convective Heat Transer through Two-sided Vertical
More informationFundamentals of Fluid Mechanics
Sixth Edition Fundamentals of Fluid Mechanics International Student Version BRUCE R. MUNSON DONALD F. YOUNG Department of Aerospace Engineering and Engineering Mechanics THEODORE H. OKIISHI Department
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More information12d Model. Civil and Surveying Software. Version 7. Drainage Analysis Module Hydraulics. Owen Thornton BE (Mech), 12d Model Programmer
1d Model Civil and Surveying Sotware Version 7 Drainage Analysis Module Hydraulics Owen Thornton BE (Mech), 1d Model Programmer owen.thornton@1d.com 9 December 005 Revised: 10 January 006 8 February 007
More informationOPTIMIZATION AND DESIGN GUIDELINES FOR HIGH FLUX MICRO-CHANNEL HEAT SINKS FOR LIQUID AND GASEOUS SINGLE-PHASE FLOW
OPTIMIZATION AND DESIGN GIDELINES FOR HIGH FLX MICRO-CHANNEL HEAT SINKS FOR LIID AND GASEOS SINGLE-PHASE FLOW Norbert Müller, Luc G. Fréchette Mechanical Engineering Columbia niversity in the City o New
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationEXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER
EXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER 1.1 AIM: To determine the co-efficient of discharge of the orifice meter 1.2 EQUIPMENTS REQUIRED: Orifice meter test rig, Stopwatch 1.3 PREPARATION 1.3.1
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationComments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium
Comments on Magnetohydrodynamic Unsteady Flow o A Non- Newtonian Fluid Through A Porous Medium Mostaa A.A.Mahmoud Department o Mathematics, Faculty o Science, Benha University (358), Egypt Abstract The
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Stokes' law and Reynold number Dr. Zifei Liu The motion of a particle in a fluid environment, such as air or water m dv =F(t) - F dt d - 1 4 2 3 πr3
More informationApplication of Method of Lines on the Bifurcation Analysis of a Horizontal Rijke Tube
th ASPACC July 9, 5 Beijing, China Application o Method o Lines on the Biurcation Analysis o a Horizontal Rijke Tube Xiaochuan Yang, Ali Turan, Adel Nasser, Shenghui Lei School o Mechanical, Aerospace
More informationME 328 Machine Design Vibration handout (vibrations is not covered in text)
ME 38 Machine Design Vibration handout (vibrations is not covered in text) The ollowing are two good textbooks or vibrations (any edition). There are numerous other texts o equal quality. M. L. James,
More informationFlow of aerosol in 3D alveolated bifurcations: experimental measurements by Particle Image Velocimetry and Particle Tracking Velocimetry
Flow o aerosol in 3D alveolated biurcations: experimental measurements by Particle Image Velocimetry and Particle Tracking Velocimetry Vincent Ruwet, Patricia Corieri, Ra Theunissen, Baoshun Ma, Michel
More informationPIPE FLOW. General Characteristic of Pipe Flow. Some of the basic components of a typical pipe system are shown in Figure 1.
PIPE FLOW General Characteristic of Pipe Flow Figure 1 Some of the basic components of a typical pipe system are shown in Figure 1. They include the pipes, the various fitting used to connect the individual
More informationSolving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah
Solving Partial Dierential Equations Numerically Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Overview What are partial dierential equations? How do we solve them? (Example)
More informationChem 406 Biophysical Chemistry Lecture 1 Transport Processes, Sedimentation & Diffusion
Chem 406 Biophysical Chemistry Lecture 1 Transport Processes, Sedimentation & Diusion I. Introduction A. There are a group o biophysical techniques that are based on transport processes. 1. Transport processes
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationAEROSPACE ENGINEERING DEPARTMENT. Second Year - Second Term ( ) Fluid Mechanics & Gas Dynamics
AEROSPACE ENGINEERING DEPARTMENT Second Year - Second Term (2008-2009) Fluid Mechanics & Gas Dynamics Similitude,Dimensional Analysis &Modeling (1) [7.2R*] Some common variables in fluid mechanics include:
More informationPart A: 1 pts each, 10 pts total, no partial credit.
Part A: 1 pts each, 10 pts total, no partial credit. 1) (Correct: 1 pt/ Wrong: -3 pts). The sum of static, dynamic, and hydrostatic pressures is constant when flow is steady, irrotational, incompressible,
More informationITB J. Eng. Sci. Vol. 41, No. 1, 2009,
ITB J. Eng. Sci. Vol. 41, No. 1, 2009, 65-76 65 ressure Drop Correlation Covering Dilute to Dense Regimes o Solid article-gas Flow in a Vertical Conveying ipe Y. Bindar, N.A. Sutrisniningrum & D. Santiani
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationHydraulic resistance at sudden pipe expansion-the influence of cavitation
Hydraulic resistance at sudden pipe expansion-the influence of cavitation I Department of Hydraulic Structures and Water Resources Management, Technical University Graz 2 Department of Hydraulic Machinery,
More informationV/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0
UNIT III FLOW THROUGH PIPES 1. List the types of fluid flow. Steady and unsteady flow Uniform and non-uniform flow Laminar and Turbulent flow Compressible and incompressible flow Rotational and ir-rotational
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationPressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers
Pressure Losses for Fluid Flow Through Abrupt Area Contraction in Compact Heat Exchangers Undergraduate Research Spring 004 By Bryan J. Johnson Under Direction of Rehnberg Professor of Ch.E. Bruce A. Finlayson
More informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationChapter (3) Water Flow in Pipes
Chapter (3) Water Flow in Pipes Water Flow in Pipes Bernoulli Equation Recall fluid mechanics course, the Bernoulli equation is: P 1 ρg + v 1 g + z 1 = P ρg + v g + z h P + h T + h L Here, we want to study
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationLetters to the editor
Letters to the editor Dear Editor, I read with considerable interest the article by Bruno Putzeys in Linear udio Vol 1, The F- word - or, why there is no such thing as too much eedback. In his excellent
More informationThe Effect of Internal Obstructions in Naturally Ventilated Greenhouse Applications
HEFAT27 5 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 1-4 July 27, Sun City, South Arica Paper number: KS2 The Eect o Internal Obstructions in Naturally Ventilated Greenhouse
More informationPerformance. 5. More Aerodynamic Considerations
Performance 5. More Aerodynamic Considerations There is an alternative way of looking at aerodynamic flow problems that is useful for understanding certain phenomena. Rather than tracking a particle of
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationMulti-Fidelity Computational Flow Assurance for Design and Development of Subsea Systems and Equipment Simon Lo
Multi-Fidelity Computational Flow Assurance for Design and Development of Subsea Systems and Equipment Simon Lo CD-adapco, Trident House, Basil Hill Road, Didcot, OX11 7HJ, UK Multi-Fidelity Computational
More informationmywbut.com Hydraulic Turbines
Hydraulic Turbines Hydro-electric power accounts for up to 0% of the world s electrical generation. Hydraulic turbines come in a variety of shapes determined by the available head and a number of sizes
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationLarge Scale Simulations of Turbulent Flows for Industrial Applications. Lakhdar Remaki BCAM- Basque Centre for Applied Mathematics
Large Scale Simulations o Turbulent Flows or Industrial Alications Lakhdar Remaki BCAM- Basque Centre or Alied Mathematics Outline Flow Motion Simulation Physical Model Finite-Volume Numerical Method CFD
More information1-Reynold s Experiment
Lect.No.8 2 nd Semester Flow Dynamics in Closed Conduit (Pipe Flow) 1 of 21 The flow in closed conduit ( flow in pipe ) is differ from this occur in open channel where the flow in pipe is at a pressure
More informationME 309 Fluid Mechanics Fall 2010 Exam 2 1A. 1B.
Fall 010 Exam 1A. 1B. Fall 010 Exam 1C. Water is flowing through a 180º bend. The inner and outer radii of the bend are 0.75 and 1.5 m, respectively. The velocity profile is approximated as C/r where C
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More information