4.4 Buoyant plume from a steady heat source

Size: px
Start display at page:

Download "4.4 Buoyant plume from a steady heat source"

Transcription

1 1 Notes on 1.63 Advanced Envionmental Fluid Mechanics Instucto: C. C. Mei, Octobe 25, buoyplum.tex 4.4 Buoyant plume fom a steady heat souce [Refeence]: Gebhat, et. al. (Jalluia, Mahajan, Saammakia), Buoyancy-induced Flows and Tanspot, 1988, Hemisphee Publishing Copoation Let T = T T = tempeatue vaiation whee T i nfty is a constant (no ambient statification). Fo a stong enough heat souce, we expect bounday laye behavio, À x,uàv, p = The bounday laye equations ae (u) x + (v) u u x + v u = gβ(t T )+ ν u T x + v T = k The centeline = is an axis of symmety, = (4.4.1) T u (4.4.2) (4.4.3) Fa outside the plume v = u = T = (4.4.4) u andt T, ( T ) (4.4.5) Rewite (4.4.3) as (u T ) x + (v T ) = (u T ) x T (u) x + (v) + (v T ) = k T (4.4.6)

2 2 afte using continuity. Now integating the last equation fom =to = Z x = k2π 2πu Td+2π d T (v T ) d theefoe Z 2πu Td+2πv T x Using the bounday conditions, we get o =2πk T = = (4.4.7) 2πu Td = constant Note that 2πd uρc T = ate of buoyancy flux = ate of heat flux = Q(given ate of heat elease at x =) theefoe, Q = This is a bounday condition. Let the steam function ψ be defined by 2πdρuC T (4.4.8) u = ψ, v = ψ x (4.4.1) is automatically satisfied. Fom the momentum equation: 1 ψ 1 2 ψ x 1 ψ x 1 ψ = gβ T + ν " 1 # ψ (4.4.9) (4.4.1) Fom the enegy equation 1 ψ T x 1 ψ T x = k T (4.4.11) and fom the buoyancy flux condition Q =2πρC d 1 ψ T (4.4.12)

3 3 Ty a similaity solution with the one-paamete tansfomation Fom (4.4.1), fom (4.4.11) and fom (4.4.12) Fom these thee equations we get x λ a x, = λ b, ψ = λ c ψ, T = λ d T λ 2c 4b a = λ 2c 4b a = λ d = λ c 4b (4.4.13) λ c+d 2b a = λ d 2b (4.4.14) λ c+d = 1 (4.4.15) c a =1, b a = 1 2, d a = 1. We leave it as an execise to show that the similaity vaiable can be taken to be and the similaity solutions to be Afte much algeba, and noting we get fom (4.4.1) and fom (4.4.11). η = 1 x, 1/2 Befoe integating, let us nomalize : It follows fom (4.4.18) that η = x 1/2 (4.4.16) ψ = xf (η), and T = x 1 G(η) (4.4.17) η x = 1 2 x = 1 1 3/2 2 x 1/2 x = η 2x F νf + (F ν)+gβηg = (4.4.18) η νγ α 3 F + γ α 3 k(ηg ) +(FG) = (4.4.19) η = α η, F = γ F, G = σḡ. (4.4.2) F (γ F ν)+gβασ ηḡ = (4.4.21) η

4 4 whee pime denotes d/d η. Settingγ = ν and ν 2 α 3 = gβασ which elates σ and α, σ = ν2 (4.4.22) gβα 4 we get F F + ( η F 1) + ηḡ = (4.4.23) Simila nomalization of (4.4.19) gives. which can be simplified to.whee kασ α 2 ( ηḡ ) + γσ α ( F Ḡ) = (4.4.24) ( ηḡ ) + P( F Ḡ) = (4.4.25) P = ν k = Pandtl Numbe (4.4.26) Fo wate ν =1 2 cm 2 /s, k =1.42cm 2 /s, hence P =7. Foaiν =.145cm 2 /s, k =.22cm 2 /s, hence P =.75. We now integate (4.4.25)to give ηḡ + P F Ḡ =constant Since ψ(x, ) =, we must have F () = ; the constant above is zeo. Equation (4.4.27) can be witten ηḡ + P F Ḡ = (4.4.27) Ḡ Ḡ = P F η, d ln Ḡ F o = P d η η Z η ln Ḡ = P F d η +constant η Z η Ḡ( η) =Ḡ() exp F P η d η (4.4.28) Substituting Eqn. (4.4.28) into Eqn. (4.4.23), the esulting equation fo F must be integated numeically. Now let us find the bounday condtions fo F o F.

5 5 Eqn. (4.4.8 ) becomes Z Q 2πρC = d Theefoe, Let us choose so that 1 ψ G(η) x = d x1/2 F G x = d η F Ḡ = Q 2πρCνσ dη(f G)=νσ d η( F Ḡ) (4.4.29) (4.4.3) Q 2πρCνσ = 1 (4.4.31) d η F Ḡ = 1 (4.4.32) is the bounday condition fo F and Ḡ. Now (4.4.31) defines σ, thescaleofg. Note that lage Q implies lage σ and smalle α. Thus a stonge heat souce leads to a geate centeline tempeatue and a thinne plume. Also, u as hence u = 1 ψ = F η = ν F α 2 η, The adial velocity is, in geneal v = 1 ψ x = 1 F η F 2 as η η Since we must have, Clealy v as η, F () =. F ( η) = as η (4.4.33) The numeical esults by Mollendof & Gelhat, 1974, ae shown in Figs , fo vaious Pandtl numbes. A schlieian photogaph due to Gebhat (copied fom Van Dyke An Album of Fluid Motion) ishowninfiguefig:plumevd. Remak: u = 1 ψ = F µ = x η F 1 x 1/2 Along the centeline u(x, ) = ³ F = constant depending on P η. Why? Buoyancy acceleation is counteacted by entainment.

6 6 Remak: Let the adius of the plume be a which vaies as a x 1/2 This is consistent with the behavio that u x,and T x 1,since On the othe hand the mass flux ate is and the momentum flux ate is a 2 u T = Q ua 2 x u 2 a 2 x hence both appoach zeo at the souce. Thus a plume is the esult of enegy souce, not of mass o momentum.

7 Figue 4.4.1: Velocity and tempeatue pofiles in a themal plume 7

8 8 Figue 4.4.2: A 2D themal plume fom a line heat souce. Gebhat, Pea and Schoo 197, Fom Van Dyke, photo by

1 Similarity Analysis

1 Similarity Analysis ME43A/538A/538B Axisymmetic Tubulent Jet 9 Novembe 28 Similaity Analysis. Intoduction Conside the sketch of an axisymmetic, tubulent jet in Figue. Assume that measuements of the downsteam aveage axial

More information

In the previous section we considered problems where the

In the previous section we considered problems where the 5.4 Hydodynamically Fully Developed and Themally Developing Lamina Flow In the pevious section we consideed poblems whee the velocity and tempeatue pofile wee fully developed, so that the heat tansfe coefficient

More information

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Intoduction Pehaps the most citicism of Pulse Neuton Neuon (PNN) logging methods has been chage that PNN is to sensitive to the

More information

One-Dimensional, Steady-State. State Conduction with Thermal Energy Generation

One-Dimensional, Steady-State. State Conduction with Thermal Energy Generation One-Dimensional, Steady-State State Conduction with Themal Enegy Geneation Implications of Enegy Geneation Involves a local (volumetic) souce of themal enegy due to convesion fom anothe fom of enegy in

More information

Incompressible Viscous Flows

Incompressible Viscous Flows Incompessible Viscous Flows F an incompessible fluid, the continuity equation and the Navie-Stokes equation ae given as v = 0, () + v v = P + ν t v Using a vect identity, Equation () may be estated as

More information

Do not turn over until you are told to do so by the Invigilator.

Do not turn over until you are told to do so by the Invigilator. UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.

More information

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty

More information

An Exact Solution of Navier Stokes Equation

An Exact Solution of Navier Stokes Equation An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in

More information

Lecture 2 - Thermodynamics Overview

Lecture 2 - Thermodynamics Overview 2.625 - Electochemical Systems Fall 2013 Lectue 2 - Themodynamics Oveview D.Yang Shao-Hon Reading: Chapte 1 & 2 of Newman, Chapte 1 & 2 of Bad & Faulkne, Chaptes 9 & 10 of Physical Chemisty I. Lectue Topics:

More information

The Combined Effect of Chemical reaction, Radiation, MHD on Mixed Convection Heat and Mass Transfer Along a Vertical Moving Surface

The Combined Effect of Chemical reaction, Radiation, MHD on Mixed Convection Heat and Mass Transfer Along a Vertical Moving Surface Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-966 Vol. 5, Issue (Decembe ), pp. 53 53 (Peviously, Vol. 5, Issue, pp. 63 6) Applications and Applied Mathematics: An Intenational Jounal (AAM)

More information

T. Raja Rani. Military Technological College, Muscat, Oman. Abstract

T. Raja Rani. Military Technological College, Muscat, Oman. Abstract ISSN: 78-8 Vol. Issue, Octobe - Fee Convection ove a Vaying all Vetical Cylinde embedded in a Poous medium ith effect of Radiation, Vaiable Fluid Popeties and Statification. T. Raja Rani Militay Technological

More information

7.2.1 Basic relations for Torsion of Circular Members

7.2.1 Basic relations for Torsion of Circular Members Section 7. 7. osion In this section, the geomety to be consideed is that of a long slende cicula ba and the load is one which twists the ba. Such poblems ae impotant in the analysis of twisting components,

More information

I. CONSTRUCTION OF THE GREEN S FUNCTION

I. CONSTRUCTION OF THE GREEN S FUNCTION I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,

More information

Superposition. Section 8.5.3

Superposition. Section 8.5.3 Supeposition Section 8.5.3 Simple Potential Flows Most complex potential (invicid, iotational) flows can be modeled using a combination of simple potential flows The simple flows used ae: Unifom flows

More information

MONTE CARLO SIMULATION OF FLUID FLOW

MONTE CARLO SIMULATION OF FLUID FLOW MONTE CARLO SIMULATION OF FLUID FLOW M. Ragheb 3/7/3 INTRODUCTION We conside the situation of Fee Molecula Collisionless and Reflective Flow. Collisionless flows occu in the field of aefied gas dynamics.

More information

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES Poceedings: Building Simulation 007 NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES A. Satipi, A. Laouadi, D. Naylo 3, R. Dhib 4 Depatment of Mechanical and Industial

More information

Analytical evaluation of 3D BEM integral representations using complex analysis

Analytical evaluation of 3D BEM integral representations using complex analysis BIR Wokshop 5w5052 Moden Applications of Complex Vaiables: Modeling, Theoy and Computation Analytical evaluation of 3D BEM integal epesentations using complex analysis onia Mogilevskaya Depatment of Civil,

More information

PROBLEM SET #3A. A = Ω 2r 2 2 Ω 1r 2 1 r2 2 r2 1

PROBLEM SET #3A. A = Ω 2r 2 2 Ω 1r 2 1 r2 2 r2 1 PROBLEM SET #3A AST242 Figue 1. Two concentic co-axial cylindes each otating at a diffeent angula otation ate. A viscous fluid lies between the two cylindes. 1. Couette Flow A viscous fluid lies in the

More information

1D2G - Numerical solution of the neutron diffusion equation

1D2G - Numerical solution of the neutron diffusion equation DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both

More information

f h = u, h g = v, we have u + v = f g. So, we wish

f h = u, h g = v, we have u + v = f g. So, we wish Answes to Homewok 4, Math 4111 (1) Pove that the following examples fom class ae indeed metic spaces. You only need to veify the tiangle inequality. (a) Let C be the set of continuous functions fom [0,

More information

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in Supplementay Figue 1. Cicula paallel lamellae gain size as a function of annealing time at 50 C. Eo bas epesent the σ uncetainty in the measued adii based on image pixilation and analysis uncetainty contibutions

More information

Chapter 5 Page 5.1 CHAPTER 5. r Force times distance has units of energy. Therefore, fxr=u, or f / is dimensionless.

Chapter 5 Page 5.1 CHAPTER 5. r Force times distance has units of energy. Therefore, fxr=u, or f / is dimensionless. Chapte 5 Page 5.1 CHAPTER 5 Poblem 5.1: 1 (a) u () 4 0.90.93 3.0 (b) Foce times distance has units of enegy. Theefoe, fx=u, o f/ is dimensionless. d f = d u 1 d f 4ε 1 = f = 4 ε1 d 13 f = 4 ε 1 f ε = 4

More information

Physics 121 Hour Exam #5 Solution

Physics 121 Hour Exam #5 Solution Physics 2 Hou xam # Solution This exam consists of a five poblems on five pages. Point values ae given with each poblem. They add up to 99 points; you will get fee point to make a total of. In any given

More information

EFFECT OF A TEMPERATURE-DEPENDENT THERMAL CONDUCTIVITY ON A FIXED UNBOUNDED SOLID WITH A CYLINDRICAL CAVITY

EFFECT OF A TEMPERATURE-DEPENDENT THERMAL CONDUCTIVITY ON A FIXED UNBOUNDED SOLID WITH A CYLINDRICAL CAVITY U.P.B. Sci. Bull., Seies A, Vol. 78, Iss. 4, 016 ISSN 13-707 EFFECT OF A TEMPERATURE-DEPENDENT THERMAL CONDUCTIVITY ON A FIXED UNBOUNDED SOLID WITH A CYLINDRICAL CAVITY Ashaf M. ZENKOUR 1 This aticle investigates

More information

8 Separation of Variables in Other Coordinate Systems

8 Separation of Variables in Other Coordinate Systems 8 Sepaation of Vaiables in Othe Coodinate Systems Fo the method of sepaation of vaiables to succeed you need to be able to expess the poblem at hand in a coodinate system in which the physical boundaies

More information

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba

More information

is the instantaneous position vector of any grid point or fluid

is the instantaneous position vector of any grid point or fluid Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in

More information

Chapter 3: Theory of Modular Arithmetic 38

Chapter 3: Theory of Modular Arithmetic 38 Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences

More information

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2 THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace

More information

Dymore User s Manual Two- and three dimensional dynamic inflow models

Dymore User s Manual Two- and three dimensional dynamic inflow models Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1

More information

A Cross Section surface tension viscosity σ U 2 10 inertia gravity gd

A Cross Section surface tension viscosity σ U 2 10 inertia gravity gd .5 Final Exam 005 SOLUTION Question U A Coss Section Photo emoved fo copyight easons. Souce: Figue 7l in Clanet, C. "Dynamics and stability of wate bells." J. Fluid Mech 40 (00): -47. R d Tooidal im U

More information

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C Umeå Univesitet, Fysik 1 Vitaly Bychkov Pov i teknisk fysik, Fluid Dynamics (Stömningsläa), 2013-05-31, kl 9.00-15.00 jälpmedel: Students may use any book including the textbook Lectues on Fluid Dynamics.

More information

ASTR415: Problem Set #6

ASTR415: Problem Set #6 ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal

More information

Sensor and Simulation Notes. Note 525. Oct Lens Design for a Prolate-Spheroidal Impulse radiating Antenna (IRA)

Sensor and Simulation Notes. Note 525. Oct Lens Design for a Prolate-Spheroidal Impulse radiating Antenna (IRA) Senso and Simulation Notes Note 55 Oct 7 Lens Design fo a Polate-Spheoidal Impulse adiating Antenna (IRA) Sehat Altunc, Cal E. Baum, Chistos G. Chistodoulou and Edl Schamiloglu Univesity of New Mexico

More information

Computational Methods of Solid Mechanics. Project report

Computational Methods of Solid Mechanics. Project report Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate

More information

Chemical Reaction and Thermal Diffusion Effects on Mass Transfer Flow through an Inclined Plate

Chemical Reaction and Thermal Diffusion Effects on Mass Transfer Flow through an Inclined Plate Open Jounal of Fluid Dynamics, 6, 6, 6-74 Published Online Mach 6 in SciRes. http://www.scip.og/jounal/ojfd http://dx.doi.og/.436/ojfd.6.66 Chemical Reaction and Themal Diffusion Effects on Mass Tansfe

More information

EM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1)

EM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1) EM- Coulomb s law, electic field, potential field, supeposition q ' Electic field of a point chage ( ') E( ) kq, whee k / 4 () ' Foce of q on a test chage e at position is ee( ) Electic potential O kq

More information

Chemical Engineering 412

Chemical Engineering 412 Chemical Engineeing 41 Intoductoy Nuclea Engineeing Lectue 16 Nuclea eacto Theoy III Neuton Tanspot 1 One-goup eacto Equation Mono-enegetic neutons (Neuton Balance) DD φφ aa φφ + ss 1 vv vv is neuton speed

More information

2. Plane Elasticity Problems

2. Plane Elasticity Problems S0 Solid Mechanics Fall 009. Plane lasticity Poblems Main Refeence: Theoy of lasticity by S.P. Timoshenko and J.N. Goodie McGaw-Hill New Yok. Chaptes 3..1 The plane-stess poblem A thin sheet of an isotopic

More information

MECHANICAL PULPING REFINER MECHANICAL PULPS

MECHANICAL PULPING REFINER MECHANICAL PULPS MECHANICAL PULPING REFINER MECHANICAL PULPS Histoy of efine mechanical pulping Fo many yeas all mechanical pulp was made fom stone goundwood (SGW). This equied whole logs. Stating in the 950s, but eally

More information

Phys-272 Lecture 18. Mutual Inductance Self-Inductance R-L Circuits

Phys-272 Lecture 18. Mutual Inductance Self-Inductance R-L Circuits Phys-7 ectue 8 Mutual nductance Self-nductance - Cicuits Mutual nductance f we have a constant cuent i in coil, a constant magnetic field is ceated and this poduces a constant magnetic flux in coil. Since

More information

Mobility of atoms and diffusion. Einstein relation.

Mobility of atoms and diffusion. Einstein relation. Mobility of atoms and diffusion. Einstein elation. In M simulation we can descibe the mobility of atoms though the mean squae displacement that can be calculated as N 1 MS ( t ( i ( t i ( 0 N The MS contains

More information

Absolute Specifications: A typical absolute specification of a lowpass filter is shown in figure 1 where:

Absolute Specifications: A typical absolute specification of a lowpass filter is shown in figure 1 where: FIR FILTER DESIGN The design of an digital filte is caied out in thee steps: ) Specification: Befoe we can design a filte we must have some specifications. These ae detemined by the application. ) Appoximations

More information

Solution to HW 3, Ma 1a Fall 2016

Solution to HW 3, Ma 1a Fall 2016 Solution to HW 3, Ma a Fall 206 Section 2. Execise 2: Let C be a subset of the eal numbes consisting of those eal numbes x having the popety that evey digit in the decimal expansion of x is, 3, 5, o 7.

More information

Physics 181. Assignment 4

Physics 181. Assignment 4 Physics 181 Assignment 4 Solutions 1. A sphee has within it a gavitational field given by g = g, whee g is constant and is the position vecto of the field point elative to the cente of the sphee. This

More information

CBE Transport Phenomena I Final Exam. December 19, 2013

CBE Transport Phenomena I Final Exam. December 19, 2013 CBE 30355 Tanspot Phenomena I Final Exam Decembe 9, 203 Closed Books and Notes Poblem. (20 points) Scaling analysis of bounday laye flows. A popula method fo measuing instantaneous wall shea stesses in

More information

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below. Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings

More information

The second law of thermodynamics - II.

The second law of thermodynamics - II. Januay 21, 2013 The second law of themodynamics - II. Asaf Pe e 1 1. The Schottky defect At absolute zeo tempeatue, the atoms of a solid ae odeed completely egulaly on a cystal lattice. As the tempeatue

More information

Scaling laws for dynamos in rotating spherical shells and applications to planetary magnetic fields

Scaling laws for dynamos in rotating spherical shells and applications to planetary magnetic fields Scaling laws fo dynamos in otating spheical shells and applications to planetay magnetic fields Ulich Chistensen Max-Planck-Institute fo Sola System Reseach, Katlenbug-Lindau, Gemany In collaboation with

More information

rt () is constant. We know how to find the length of the radius vector by r( t) r( t) r( t)

rt () is constant. We know how to find the length of the radius vector by r( t) r( t) r( t) Cicula Motion Fom ancient times cicula tajectoies hae occupied a special place in ou model of the Uniese. Although these obits hae been eplaced by the moe geneal elliptical geomety, cicula motion is still

More information

Motithang Higher Secondary School Thimphu Thromde Mid Term Examination 2016 Subject: Mathematics Full Marks: 100

Motithang Higher Secondary School Thimphu Thromde Mid Term Examination 2016 Subject: Mathematics Full Marks: 100 Motithang Highe Seconday School Thimphu Thomde Mid Tem Examination 016 Subject: Mathematics Full Maks: 100 Class: IX Witing Time: 3 Hous Read the following instuctions caefully In this pape, thee ae thee

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position

More information

A Relativistic Electron in a Coulomb Potential

A Relativistic Electron in a Coulomb Potential A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify

More information

EM Boundary Value Problems

EM Boundary Value Problems EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do

More information

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,

More information

Three dimensional flow analysis in Axial Flow Compressors

Three dimensional flow analysis in Axial Flow Compressors 1 Thee dimensional flow analysis in Axial Flow Compessos 2 The ealie assumption on blade flow theoies that the flow inside the axial flow compesso annulus is two dimensional means that adial movement of

More information

A thermodynamic degree of freedom solution to the galaxy cluster problem of MOND. Abstract

A thermodynamic degree of freedom solution to the galaxy cluster problem of MOND. Abstract A themodynamic degee of feedom solution to the galaxy cluste poblem of MOND E.P.J. de Haas (Paul) Nijmegen, The Nethelands (Dated: Octobe 23, 2015) Abstact In this pape I discus the degee of feedom paamete

More information

2 E. on each of these two surfaces. r r r r. Q E E ε. 2 2 Qencl encl right left 0

2 E. on each of these two surfaces. r r r r. Q E E ε. 2 2 Qencl encl right left 0 Ch : 4, 9,, 9,,, 4, 9,, 4, 8 4 (a) Fom the diagam in the textbook, we see that the flux outwad though the hemispheical suface is the same as the flux inwad though the cicula suface base of the hemisphee

More information

International Journal of Mathematical Archive-9(5), 2018, Available online through ISSN

International Journal of Mathematical Archive-9(5), 2018, Available online through   ISSN Intenational Jounal of Mathematical Achive-9(5), 18, -3 Available online though www.ijma.info ISSN 9 546 MHD FLOW OVER A MOVING INFINITE VERTICAL POROUS PLATE WITH UNIFORM HEAT FLUX IN THE PRESENCE OF

More information

Many Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued

Many Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued Many Electon Atoms The many body poblem cannot be solved analytically. We content ouselves with developing appoximate methods that can yield quite accuate esults (but usually equie a compute). The electons

More information

Numerical Calculation of Circulation Flow Rate in the Degassing Rheinstahl-Heraeus Process

Numerical Calculation of Circulation Flow Rate in the Degassing Rheinstahl-Heraeus Process ISIJ Intenational, Vol. 40 (2000), No. 8, pp. 749 755 Numeical Calculation of Ciculation Flow Rate in the Degassing Rheinstahl-Heaeus Pocess Young-Geun PARK, Won-Chul DOO, Kyung-Woo YI and Sang-Bok AN

More information

THERMODYNAMIC OPTIMIZATION OF TUBULAR HEAT EXCHANGERS BASED ON MINIMUM IRREVERSIBILITY CRITERIA

THERMODYNAMIC OPTIMIZATION OF TUBULAR HEAT EXCHANGERS BASED ON MINIMUM IRREVERSIBILITY CRITERIA THERMODYNAMIC OPTIMIZATION OF TUBULAR HEAT EXCHANGER BAED ON MINIMUM IRREVERIBILITY CRITERIA As. dd. ing. Adina GHEORGHIAN, Pof. d. ing. Alexandu DOBROVICECU, As. dd. ing. Andeea MARIN,.l. d. ing. Claudia

More information

Lecture 7: Angular Momentum, Hydrogen Atom

Lecture 7: Angular Momentum, Hydrogen Atom Lectue 7: Angula Momentum, Hydogen Atom Vecto Quantization of Angula Momentum and Nomalization of 3D Rigid Roto wavefunctions Conside l, so L 2 2 2. Thus, we have L 2. Thee ae thee possibilities fo L z

More information

1) Consider an object of a parabolic shape with rotational symmetry z

1) Consider an object of a parabolic shape with rotational symmetry z Umeå Univesitet, Fysik 1 Vitaly Bychkov Pov i teknisk fysik, Fluid Mechanics (Stömningsläa), 01-06-01, kl 9.00-15.00 jälpmedel: Students may use any book including the tetbook Lectues on Fluid Dynamics.

More information

15 Solving the Laplace equation by Fourier method

15 Solving the Laplace equation by Fourier method 5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the

More information

Heat transfer has direction as well as magnitude. The rate of heat conduction

Heat transfer has direction as well as magnitude. The rate of heat conduction cen58933_ch2.qd 9/1/22 8:46 AM Page 61 HEAT CONDUCTION EQUATION CHAPTER 2 Heat tansfe has diection as well as magnitude. The ate of heat conduction in a specified diection is popotional to the tempeatue

More information

Newton s Laws, Kepler s Laws, and Planetary Orbits

Newton s Laws, Kepler s Laws, and Planetary Orbits Newton s Laws, Keple s Laws, and Planetay Obits PROBLEM SET 4 DUE TUESDAY AT START OF LECTURE 28 Septembe 2017 ASTRONOMY 111 FALL 2017 1 Newton s & Keple s laws and planetay obits Unifom cicula motion

More information

T x. T k x. is a constant of integration. We integrate a second time to obtain an expression for the temperature distribution:

T x. T k x. is a constant of integration. We integrate a second time to obtain an expression for the temperature distribution: ME 336 Fall 8 HW solution Poblem - The geneal fom of the heat diffusion equation is: T cp = ( T) + eg t - one-dimensional conduction (along the x - diection only): = ˆi and T = T( x) x - steady state conditions:

More information

Scattering in Three Dimensions

Scattering in Three Dimensions Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.

More information

Review: Electrostatics and Magnetostatics

Review: Electrostatics and Magnetostatics Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion

More information

Energy Levels Of Hydrogen Atom Using Ladder Operators. Ava Khamseh Supervisor: Dr. Brian Pendleton The University of Edinburgh August 2011

Energy Levels Of Hydrogen Atom Using Ladder Operators. Ava Khamseh Supervisor: Dr. Brian Pendleton The University of Edinburgh August 2011 Enegy Levels Of Hydogen Atom Using Ladde Opeatos Ava Khamseh Supeviso: D. Bian Pendleton The Univesity of Edinbugh August 11 1 Abstact The aim of this pape is to fist use the Schödinge wavefunction methods

More information

Numerical Integration

Numerical Integration MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.

More information

c n ψ n (r)e ient/ h (2) where E n = 1 mc 2 α 2 Z 2 ψ(r) = c n ψ n (r) = c n = ψn(r)ψ(r)d 3 x e 2r/a0 1 πa e 3r/a0 r 2 dr c 1 2 = 2 9 /3 6 = 0.

c n ψ n (r)e ient/ h (2) where E n = 1 mc 2 α 2 Z 2 ψ(r) = c n ψ n (r) = c n = ψn(r)ψ(r)d 3 x e 2r/a0 1 πa e 3r/a0 r 2 dr c 1 2 = 2 9 /3 6 = 0. Poblem {a} Fo t : Ψ(, t ψ(e iet/ h ( whee E mc α (α /7 ψ( e /a πa Hee we have used the gound state wavefunction fo Z. Fo t, Ψ(, t can be witten as a supeposition of Z hydogenic wavefunctions ψ n (: Ψ(,

More information

Averaging method for nonlinear laminar Ekman layers

Averaging method for nonlinear laminar Ekman layers Unde consideation fo publication in J. Fluid Mech. 1 Aveaging method fo nonlinea lamina Ekman layes By A. A N D E R S E N 1, 2, B. L A U T R U P 3 A N D T. B O H R 1 1 The Technical Univesity of Denmak,

More information

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4 Fouie tansfom pais: f(x) 1 f(k) e p 2 (k) p e iax 2 (k a) 2 e a x a a 2 + k 2 e a2 x 1 2, a > 0 a p k2 /4a2 e 2 1 H(x) ik p 2 + 2 (k) The fist few Y m Y 0 0 = Y 0 1 = Y ±1 1 = l : 1 Y2 0 = 4 3 ±1 cos Y

More information

Modeling Fermi Level Effects in Atomistic Simulations

Modeling Fermi Level Effects in Atomistic Simulations Mat. Res. Soc. Symp. Poc. Vol. 717 Mateials Reseach Society Modeling Femi Level Effects in Atomistic Simulations Zudian Qin and Scott T. Dunham Depatment of Electical Engineeing, Univesity of Washington,

More information

1 Fundamental Solutions to the Wave Equation

1 Fundamental Solutions to the Wave Equation 1 Fundamental Solutions to the Wave Equation Physical insight in the sound geneation mechanism can be gained by consideing simple analytical solutions to the wave equation. One example is to conside acoustic

More information

On the Sun s Electric-Field

On the Sun s Electric-Field On the Sun s Electic-Field D. E. Scott, Ph.D. (EE) Intoduction Most investigatos who ae sympathetic to the Electic Sun Model have come to agee that the Sun is a body that acts much like a esisto with a

More information

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc. Numeical solution of diffusion mass tansfe model in adsoption systems Pof., D.Sc. Agenda Mass Tansfe Mechanisms Diffusion Mass Tansfe Models Solving Diffusion Mass Tansfe Models Paamete Estimation 2 Mass

More information

Math 124B February 02, 2012

Math 124B February 02, 2012 Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial

More information

Water flows through the voids in a soil which are interconnected. This flow may be called seepage, since the velocities are very small.

Water flows through the voids in a soil which are interconnected. This flow may be called seepage, since the velocities are very small. Wate movement Wate flows though the voids in a soil which ae inteconnected. This flow may be called seepage, since the velocities ae vey small. Wate flows fom a highe enegy to a lowe enegy and behaves

More information

Chapter 3 Optical Systems with Annular Pupils

Chapter 3 Optical Systems with Annular Pupils Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The

More information

Germán Sierra Instituto de Física Teórica CSIC-UAM, Madrid, Spain

Germán Sierra Instituto de Física Teórica CSIC-UAM, Madrid, Spain Gemán Siea Instituto de Física Teóica CSIC-UAM, Madid, Spain Wo in pogess done in collaboation with J. Lins and S. Y. Zhao (Univ. Queensland, Austalia) and M. Ibañez (IFT, Madid) Taled pesented at the

More information

Nuclear size corrections to the energy levels of single-electron atoms

Nuclear size corrections to the energy levels of single-electron atoms Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea

More information

Basic Bridge Circuits

Basic Bridge Circuits AN7 Datafoth Copoation Page of 6 DID YOU KNOW? Samuel Hunte Chistie (784-865) was bon in London the son of James Chistie, who founded Chistie's Fine At Auctionees. Samuel studied mathematics at Tinity

More information

Stellar Structure and Evolution

Stellar Structure and Evolution Stella Stuctue and Evolution Theoetical Stella odels Conside each spheically symmetic shell of adius and thickness d. Basic equations of stella stuctue ae: 1 Hydostatic equilibium π dp dp d G π = G =.

More information

=0, (x, y) Ω (10.1) Depending on the nature of these boundary conditions, forced, natural or mixed type, the elliptic problems are classified as

=0, (x, y) Ω (10.1) Depending on the nature of these boundary conditions, forced, natural or mixed type, the elliptic problems are classified as Chapte 1 Elliptic Equations 1.1 Intoduction The mathematical modeling of steady state o equilibium phenomena geneally esult in to elliptic equations. The best example is the steady diffusion of heat in

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lectue 5 Cental Foce Poblem (Chapte 3) What We Did Last Time Intoduced Hamilton s Pinciple Action integal is stationay fo the actual path Deived Lagange s Equations Used calculus

More information

1) Emits radiation at the maximum intensity possible for every wavelength. 2) Completely absorbs all incident radiation (hence the term black ).

1) Emits radiation at the maximum intensity possible for every wavelength. 2) Completely absorbs all incident radiation (hence the term black ). Radiation laws Blackbody adiation Planck s Law Any substance (solid, liquid o gas) emits adiation accoding to its absolute tempeatue, measued in units of Kelvin (K = o C + 73.5). The efficiency at which

More information

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing

More information

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS LSN 10-: MOTION IN A GRAVITATIONAL FIELD Questions Fom Reading Activity? Gavity Waves? Essential Idea: Simila appoaches can be taken in analyzing electical

More information

Chapter 13 Gravitation

Chapter 13 Gravitation Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects

More information

Goodness-of-fit for composite hypotheses.

Goodness-of-fit for composite hypotheses. Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test

More information

CHAPTER 25 ELECTRIC POTENTIAL

CHAPTER 25 ELECTRIC POTENTIAL CHPTE 5 ELECTIC POTENTIL Potential Diffeence and Electic Potential Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic foce on the paticle given by F=E. When

More information

Pulse Neutron Neutron (PNN) tool logging for porosity

Pulse Neutron Neutron (PNN) tool logging for porosity Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Hotwell Handelsges.m.b.H Oedenbuge Stasse 6 7013 Klingenbach, AUSTRIA Tel.: +43 (0) 687-48058 Fax: +43 (0) 687 48059 office@hotwell.at

More information

/6/4 5 Stuctue-Induced Sediment Scou ef: Eosion and Sedimentation, P.Y. Julien, 998 The Mechanics of Scou in the Maine Envionment, B.M. Sume and J. Fedsoe, Evaluating Scou at Bidges (HEC-8), E.. ichadson

More information

Physics 235 Chapter 5. Chapter 5 Gravitation

Physics 235 Chapter 5. Chapter 5 Gravitation Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus

More information

4. Some Applications of first order linear differential

4. Some Applications of first order linear differential August 30, 2011 4-1 4. Some Applications of fist ode linea diffeential Equations The modeling poblem Thee ae seveal steps equied fo modeling scientific phenomena 1. Data collection (expeimentation) Given

More information

Stress, Cauchy s equation and the Navier-Stokes equations

Stress, Cauchy s equation and the Navier-Stokes equations Chapte 3 Stess, Cauchy s equation and the Navie-Stokes equations 3. The concept of taction/stess Conside the volume of fluid shown in the left half of Fig. 3.. The volume of fluid is subjected to distibuted

More information

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions EN10: Continuum Mechanics Homewok 7: Fluid Mechanics Solutions School of Engineeing Bown Univesity 1. An ideal fluid with mass density ρ flows with velocity v 0 though a cylindical tube with cosssectional

More information