Lecture # 12: Shock Waves and De Laval Nozzle

Similar documents
In this lecture... Subsonic and supersonic nozzles Working of these nozzles Performance parameters for nozzles

AerE 344: Undergraduate Aerodynamics and Propulsion Laboratory. Lab Instructions

Lecture # 12: Shock Waves and De Laval Nozzle

AerE 344: Undergraduate Aerodynamics and Propulsion Laboratory. Lab Instructions. Pressure Measurements in a de Laval Nozzle

A Propagating Wave Packet Group Velocity Dispersion

State-space behaviours 2 using eigenvalues

Construction of Mimetic Numerical Methods

San José State University Aerospace Engineering AE 138 Vector-Based Dynamics for Aerospace Applications, Fall 2016

4.2 Design of Sections for Flexure

ANALYTICAL PARAMETRIC CYCLE ANALYSIS OF CONTINUOUS ROTATING DETONATION EJECTOR-AUGMENTED ROCKET ENGINE HUAN V. CAO

5.80 Small-Molecule Spectroscopy and Dynamics

Laboratory work # 8 (14) EXPERIMENTAL ESTIMATION OF CRITICAL STRESSES IN STRINGER UNDER COMPRESSION

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

Outline. Thanks to Ian Blockland and Randy Sobie for these slides Lifetimes of Decaying Particles Scattering Cross Sections Fermi s Golden Rule

A General Thermal Equilibrium Discharge Flow Model

A Simplified Theory of Microwave Pulse Compression

Utilization of Universal Controls Analysis Tool and CFD in Parallel for a Cold Nozzle Flow

A New Approach to the Fatigue Life Prediction for Notched Components Under Multiaxial Cyclic Loading. Zhi-qiang TAO and De-guang SHANG *

Scattering States of l-wave Schrödinger Equation with Modified Rosen Morse Potential

as a derivative. 7. [3.3] On Earth, you can easily shoot a paper clip straight up into the air with a rubber band. In t sec

EXST Regression Techniques Page 1

Forces. Quantum ElectroDynamics. α = = We have now:

Title: Vibrational structure of electronic transition

Evaluating Reliability Systems by Using Weibull & New Weibull Extension Distributions Mushtak A.K. Shiker

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

The second condition says that a node α of the tree has exactly n children if the arity of its label is n.

September 23, Honors Chem Atomic structure.notebook. Atomic Structure

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

ECE507 - Plasma Physics and Applications

A coupled-channel calculation of the B 3 Σ u X 3 Σ g spectrum of S 2

University of Illinois at Chicago Department of Physics. Thermodynamics & Statistical Mechanics Qualifying Examination

EXPONENTIAL ENTROPY ON INTUITIONISTIC FUZZY SETS

Definition1: The ratio of the radiation intensity in a given direction from the antenna to the radiation intensity averaged over all directions.

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

1973 AP Calculus AB: Section I

IYPT 2000 Problem No. 3 PLASMA

Construction of asymmetric orthogonal arrays of strength three via a replacement method

(Upside-Down o Direct Rotation) β - Numbers

Procdings of IC-IDC0 ( and (, ( ( and (, and (f ( and (, rspctivly. If two input signals ar compltly qual, phas spctra of two signals ar qual. That is

4 x 4, and. where x is Town Square

On the Hamiltonian of a Multi-Electron Atom

Properties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator

E hf. hf c. 2 2 h 2 2 m v f ' f 2f ' f cos c

Pipe flow friction, small vs. big pipes

ME311 Machine Design

Modern Physics. Unit 5: Schrödinger s Equation and the Hydrogen Atom Lecture 5.6: Energy Eigenvalues of Schrödinger s Equation for the Hydrogen Atom

COMPUTATIONAL NUCLEAR THERMAL HYDRAULICS

Content Skills Assessments Lessons. Identify, classify, and apply properties of negative and positive angles.

Numerical considerations regarding the simulation of an aircraft in the approaching phase for landing

1 Isoparametric Concept

General Notes About 2007 AP Physics Scoring Guidelines

Physics 2D Lecture Slides Lecture 14: Feb 1 st 2005

5.62 Physical Chemistry II Spring 2008

That is, we start with a general matrix: And end with a simpler matrix:

PHA 5127 Answers Homework 2 Fall 2001

15. Stress-Strain behavior of soils

On the irreducibility of some polynomials in two variables

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

Differential Equations

LINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM

Properties of Quarks ( ) Isospin. π = 1, 1

2008 AP Calculus BC Multiple Choice Exam

Numerical Investigation of Nozzle Shape Effect on Shock Wave in Natural Gas Processing Abstract

0 +1e Radionuclides - can spontaneously emit particles and radiation which can be expressed by a nuclear equation.

Statistical Thermodynamics: Sublimation of Solid Iodine

The Matrix Exponential

Effects of Electron Model on Three-Grid Ion Engine Analyses

Sec 2.3 Modeling with First Order Equations

ON RIGHT(LEFT) DUO PO-SEMIGROUPS. S. K. Lee and K. Y. Park

Finite Element Models for Steady Flows of Viscous Incompressible Fluids

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Square of Hamilton cycle in a random graph


CALIBRATION EXPERIMENTS OF A NEW ACTIVE FAST RESPONSE HEAT FLUX SENSOR TO MEASURE TOTAL TEMPERATURE FLUCTUATIONS

Determination of the cutting resistance model of soil. Jozef Karol Szymanski* and Saleh Balideh

The Matrix Exponential

Application of Vague Soft Sets in students evaluation

Finite Element Model of a Ferroelectric

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

What are those βs anyway? Understanding Design Matrix & Odds ratios

Problem Set 6 Solutions

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Diffractive Dijet Production with Leading Proton in ep Collisions at HERA

1997 AP Calculus AB: Section I, Part A

Development of In Situ Permeability Test Using Constant Discharge Method for Sandy Soils

Ultimate lateral load resistance of laterally loaded pile

VALUING SURRENDER OPTIONS IN KOREAN INTEREST INDEXED ANNUITIES

Clausius-Clapeyron Equation

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

+ f. e f. Ch. 8 Inflation, Interest Rates & FX Rates. Purchasing Power Parity. Purchasing Power Parity

Prod.C [A] t. rate = = =

MEMORIAL UNIVERSITY OF NEWFOUNDLAND

Mechanical Properties

Intro to Nuclear and Particle Physics (5110)

CE 530 Molecular Simulation

The University of Alabama in Huntsville Electrical and Computer Engineering Homework #4 Solution CPE Spring 2008

What are molecular simulations? Introduction. Numerical Tools : Fluid Mechanics

INC 693, 481 Dynamics System and Modelling: The Language of Bound Graphs Dr.-Ing. Sudchai Boonto Assistant Professor

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

Transcription:

ArE 311L & ArE343L Lctur Nots Lctur # 1: Shock Wavs and D Laval Nozzl Dr. Hui H Hu Dpartmnt of Arospac Enginring Iowa Stat Univrsity Ams, Iowa 50011, U.S.A

ArE311L Lab#3: rssur Masurmnts in a d Laval Nozzl Tank with comprssd air Tst sction Tap No. Distanc downstram of throat (inchs) Ara (Sq. inchs) 1-4.00 0.800-1.50 0.59 3-0.30 0.480 4-0.18 0.478 5 0.00 0.476 6 0.15 0.497 7 0.30 0.518 8 0.45 0.539 9 0.60 0.560 10 0.75 0.581 11 0.90 0.599 1 1.05 0.616 13 1.0 0.67 14 1.35 0.63 15 1.45 0.634

1st, nd and 3 rd critic conditions st critic condition 0 incrasing 1st critic condition 0 incrasing 3st critic condition

ArE311L Lab#4: rssur Masurmnts in a d Laval Nozzl 1. Undr-xpandd xpandd flow. 3rd critical 3. Ovr-xpandd xpandd flow with obliqu shocks 4. nd critical 5. Normal shock xisting insid th nozzl 6. 1st critical Rquird lots: lots of th masurd prssur (static and total prssur) as a function of distanc along th nozzl axis for th cass, 4, 5 and 6. lots of th thortically prdictd prssur (static and total prssur) as a function of distanc along th nozzl axis for th cass, 4, 5 and 6. lots with th masurd and prdictd valus of th wall prssur distribution on th sam graphs for th cass, 4, 5 and 6 for comparison. lots of th thortically prdictd and masurd Mach numbr as a function of distanc along th nozzl axis for th cass, 4, 5 and 6.

1st, nd and 3 rd critic conditions nd critical shock is at nozzl xit Flow clos to 3 rd critical Ovr-xpandd xpandd flow with shock btwn nozzl xit and throat Undr- xpandd flow Ovr- xpandd flow 1 st critical shock is almost at th nozzl throat.

rdiction of th rssur Distribution within a D Laval Nozzl by using Numrical Approach Throat, A* or A t T1 M<1 M>1 Shock A S T1 1 T M<1 A = atm = atm Using th ara ratio, th Mach numbr at any point up to th shock can b dtrmind: A * A = γ + 1 1 1 1 γ γ 1 + M M γ + 1 Aftr finding Mach numbr at front of shock, calculat Mach numbr aftr shock using: M γ 1 1+ M1 = γ 1 γ M1 Thn, calculat th A * γ 1 A = M As 1+ M γ + 1 ( ) γ + 1 γ 1 which allows us calculat th rmaining Mach numbr distribution A * A = γ + 1 1 1 1 γ γ 1 + M M γ + 1

rdiction of th rssur Distribution within a D Laval Nozzl by using Numrical Approach Throat, A* or A t T1 Shock A S 1 T1 T M<1 M>1 M<1 T A = atm = atm a.for 3rd Critical

rdiction of th rssur Distribution within a D Laval Nozzl by using Shock Tabl Mthod Throat, A* or A t T1 M<1 M>1 Shock A S T T1 1 T M<1 A = atm = atm mthod #1, by using quations: If th shockwav is locatd at position of tab#1: Mthod #: by using Isntropic Flow proprtis tabl (Appndix-A A of Andrson s s txtbook)

rdiction of th rssur Distribution within a D Laval Nozzl by using Shock Tabl Mthod By using th normal shock tabls with M1 = 1.64 w find that M = 0.686. (Appndix-B B of Andrson s s txtbook) Nxt, w find th sonic rfrnc ara bhind th shock using th ara-mach rlation. i.., M=0.686 (Appndix-A A of Andrson s s txtbook) Find sonic rfrnc bhind th shock using th ara-mach rlationship: i.., A*=0.557sq. Inchs If th shockwav is locatd at position of tab#1:

rdiction of th rssur Distribution within a D Laval Nozzl by using Shock Tabl Mthod With th xit prssur to b sa- lvl standard prssur. W now calculat th total prssur bhind th shock using this valu of xit prssur and th prssur ratio at th xit: t 1 t = = 14.7 = 19.53 0.758 Our last major task is to find th total prssur ahad of th shock, t1 = t1 1 t1 t 1 t

Som xampls of prvious studnt lab rports Thortical Data - Gaug rssur vs. osition 50 40 Gaug rssur, psi 30 0 10 0 35 30 Exprimntal 1-D Thory -10 0 1 3 4 5 osition along Nozzl Axis, inchs 3rd Critical nd Critical Normal Shock 1st Critical rssur (SI) 5 0 15 10 Gaug rssur, psi Exprimntal Rsults - rssur vs. Nozzl Location 90 80 70 60 50 40 30 0 10 0-10 0 1 3 4 5 Nozzl location, inchs 5 Total rssur = 34.336 SI 0-5 -4-3 - -1 0 1 Distanc from Troat (in) Undr-Expandd Flow Approximatly nd Critical Normal Shock