Derivation of the basic equations of fluid flows. No. Conservation of mass of a solute (applies to non-sinking particles at low concentration).

Size: px
Start display at page:

Download "Derivation of the basic equations of fluid flows. No. Conservation of mass of a solute (applies to non-sinking particles at low concentration)."

Transcription

1 Deriation of the basic eqations of flid flos. No article in the flid at this stage (net eek). Conseration of mass of the flid. Conseration of mass of a solte (alies to non-sinking articles at lo concentration). Conseration of momentm. Alication of these basic eqations to a trblent flid.

2 A fe concets before e get to the meat ensor (Stress), Vectors (e.g. osition, elocit) and scalars (e.g., S, CO ). We need to define a coordinate sstem, and an (infinitesimal) element of olme. We assme a continos flid, and that all the fields of interest are differentiable.

3 he Lagrangian frameork is the frameork in hich the las of classical mechanics are often stated. he coordinates of a oint (t) describe the trajector from ( t ). Densit,, can eole along the trajector. B the chain rle, along a arcel trajector: d d dt dt t d( ( t), t) dt D Dt tt t const const t t (, t) Conersion from Lagrangian to Elerian

4 Eamle: Let s assme that e are in a rier that feeds on glacial melt. he ater arms at a constant rate that is a fnction of distance from the sorce. If e drift don rier (A la Hckleberr Fin ), the temeratre increases ith time (D/Dt>). At one oint along the rier, hoeer, e ma see no change in temeratre ith time (/t), as the ater arriing there is alas at the same temeratre. he heat fl is adectie, (/>). In short, the conectie deriatie is: D Dt t t

5 Mass conseration (Elerian, differential aroach): Acconting for the change in mass inside a fied, constant- sie olme: (-Δ/,Δ/,Δ/) (Δ/,Δ/,Δ/) (-Δ/,-Δ/,Δ/) A ΔΔ AΔΔ AΔΔ AΔΔ Δ Δ AΔΔ (Δ/,Δ/,-Δ/) MassV A ΔΔ ( V ) ( ) ( ) A A A ( ) t 1 t Δ t (-Δ/,-Δ/,-Δ/) (Δ/,-Δ/,-Δ/) / / / / Δ Δ Δ Δ Δ / Δ / 1 ( ) ( ) 1 ( ) Δ Δ Δ / Δ / Δ / Δ / Δ / t ( ) ( ) ( ) ( ) Δ /

6 Mass conseration (Elerian, integral aroach): Acconting for the change in mass inside a fied, constant-olme olme (V ): d dt V dv t V t dv V V nds dv Where e sed the diergence theorem: It states that the olme total of all sinks and sorces, the olme integral of the diergence, is eqal to the net flo across the olmes bondar (WIKI).

7 Reiteration (no sinks/sorces): Mass conseration (Lagrangian, integral): D Dt dv V Mass conseration (Elerian, integral): d dt V dv nds V

8 Mass conseration: Note that: Can be ritten as: tt r ( ) 1D r DtD he nd term is the flid diergence (rate of otflo of olme er nit olme). his can be nonero onl for comressible flids. It is the rate of loss of densit de to comression/eansion. r For both ater and air e can assme that t in terms of their dnamics (e need comressibilit to ass sond ).

9 Mass balance for consered scalar: Adding moleclar diffsion: C dv t V r ( C K C ) S r nds Where V is the olme of the control olme and S its srface, and sing Fick s la. B the hel of the diergence theorem: V C C t r dv ( C K C) Since the olme is arbitrar, this can be tre if and onl if: C C t r ( C) ( K C)

10 Momentm balance (Naier-Stokes): Neton s nd la of motion states that the time rate of change of momentm of a article is eqal to the force acting on it. his la is Lagrangian, the time rate of change is ith resect to a reference sstem folloing the article. d dt dv gdv V ( t ) V ( t ) V ( t ) ds Where g is the bod force er nit mass (e.g. grait) and is the srface force er nit srface area bonding V. If the olme is small enogh, the integrands can be taken ot of the integral: d dt d dt d d dv dv dt dt V ( t) V ( t) d d( δ V ) δv δv δv dt dt ( ) ( δ V ) ( ) d d dt

11 Momentm balance (Naier-Stokes): he bod force is similarl il l treated: t gdv V ( t) gδv Defining a stress tensor (eanded on the net slide): Τ n And aling the diergence theorem: ds V ( t) V ( t) Τ dv ΤδV D Dt g Τ

12 Srface forcing: For an iniscid flid, the srface force eerted b the srronding flid is normal to the srface, i.e. n, and is called the ressre force. In general, iscos stress force, S, is also resent. For iscos flids: n S. B definition Τ n, and e no hae I, here S n and I is the identit tensor. For Netonian flids, Σ And the resltant Naier-Stokes eqations for incomressible flids are: D Dt g

13 Rotational smmetr: 1 ( ) 3 1

14 otal stress tensor, Netonian flid: j i δ ij j i i j ij δ i, j{1,, 3}{,, } Stokes, 1845: 1 Σ linear fnction of elocit gradients,j {,, } {,, } 1. Σ ij linear fnction of elocit gradients.. Σ ij shold anish if there is no deformation of flid elements. 3. Relationshi beteen stress and shear shold be isotroic. ( ) i ˆ ( )

15 Naier-Stokes eqations:, g Dt D 1 t 1 t 1 g t Coriolis is added hen moing the frameork to an accelerating frameork. Hae to add bondar & initial conditions.

16 Naier-Stokes eqations (Bossinessq aroimation): Searate balance of flid at rest from moing flid. ( ) (,,, ) t (,,, t) First order balance (hdrostatic): 1 g nd order balance: D g, Dt

17 Eamle: stead flo nder grait don an inclined lane. g α From: Acheson, ν g sinα BCs: 1 : cos α g h :, a

18 Eamle: stead flo nder grait don an inclined lane. g α Soltion: a g g ν ( h ) cosα ( h ) sinα Q: hat ν shold e se?

19 Renolds decomosition of the N-S eqations Assme a trblent flo. At an gien oint in sace e searate the mean flo (mean can be in time, sace, or ensemble) and deiation from the mean sch that:,,, ensemble) and deiation from the mean sch that:,, etc Sbstitting into the continit eqation (linear): ;

20 1 Sbstitting into -momentm Naier-Stokes eqation: t 1 ( ) ( ) ( ) he eoltion of the mean is forced b correlations of flctating roerties. he correlation terms time densit are the Renolds stresses. Sbstitting into a scalar conseration eqation: Sbstitting into a scalar conseration eqation C C C K C C C t C ( ) ( ) ( ) C C C t

21 Note that the Renolds stress tensor is smmetric (as is the iscos stress tensor): the iscos stress tensor): he closre roblem: to deelo eqations for the eoltion of the Renolds stresses theseles, higher order correlation are needed (e.g. ) and so on. For this reason theories hae been deised to describe ij in terms of th fl the mean flo. For more, see: htt://.cfdonline.com/wiki/introdction to trblence/renolds aeraged eqations onl ne.com/w k /Introdct on_to_trblence/renolds_aeraged_eqat ons

22 One soltion to the closre roblem is to link the Renolds stress to mean-flo Qantities. For eamle: K K edd U edd his te of formlation is aealing becase it: a. Proide for don-gradient fl. b. Is reminiscent of moleclar l diffsion and iscosit. c. Proide closre to the eqations of the the mean fields. his te of formlation is roblematic becase: a. K edd is a roert of the flo and not the flid. b. K edd is likel to ar ith orientation, nlike moleclar rocesses. Ho is K edd related to the trblence?

23 Assme a gradient in a mean roert (momentm, heat, solte, etc. Remember: no mean gradient no fl). Assme a flctating elocit field: l is the distance a arcel traels before it loses its identit. he rate of ard ertical trblent transfer of <Ψ> is don the mean gradient: ψ l ψ l ψ K edd ψ

24 Ho is K edd related to the trblence? U Kedd γ ennekes and Lmle (197) aroach this roblem from dimensional i analsis based on assming a single length scale-l and a single elocit scale ω< > 1/. c ω ; c ~ O(1) he eddies inoled in momentm transfer hae orticities, ω/l; this orticit is maintained b the mean shear (l is the length scales of the eddies, e.g. the decorrelation scale). U U ω / l c ; c ~ O(1)

25 It follos that: Ke dd ~ l ω ~ l U In analog ith momentm fl, for heat e hae: c c γ It is most commonl assmed, and erified that γ K edd.

26 Edd-diffsion: ersectie from a de atch (figres from lectre notes of Bill Yong, UCSD) De atch << dominant scale of eddies. Dashed circle denotes initial osition of atch.

27 De atch dominant scale of eddies De atch >> dominant scale of eddies

28 Cheat sheet: 1 G di f l ( ) k j i φ φ φ φ ˆ ˆ ˆ 1. Gradient of a scalar (a ector): k j i φ. Diergence of a ector (a scalar): 3 Diergence of a ensor : j i ˆ ˆ 3. Diergence of a ensor, 1 (face)(direction) : k j ˆ k

29 Cheat sheet (cotined): 4 L l i f ( ) 4. Lalacian of a ector (a ector): ( ) î ( ) ˆ ˆ k j k j

Lecture 5. Differential Analysis of Fluid Flow Navier-Stockes equation

Lecture 5. Differential Analysis of Fluid Flow Navier-Stockes equation Lectre 5 Differential Analsis of Flid Flo Naier-Stockes eqation Differential analsis of Flid Flo The aim: to rodce differential eqation describing the motion of flid in detail Flid Element Kinematics An

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

Introducing Ideal Flow

Introducing Ideal Flow D f f f p D p D p D f T k p D e The Continit eqation The Naier Stokes eqations The iscos Flo Energ Eqation These form a closed set hen to thermodnamic relations are specified Introdcing Ideal Flo Getting

More information

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined

More information

LECTURE NOTES - VI. Prof. Dr. Atıl BULU

LECTURE NOTES - VI. Prof. Dr. Atıl BULU LECTURE NOTES - VI «FLUID MECHANICS» Istanbl Technical Uniersit College of Ciil Engineering Ciil Engineering Deartment Hdralics Diision CHAPTER 6 TWO-DIMENSIONAL IDEAL FLOW 6. INTRODUCTION An ideal flid

More information

Viscous Dissipation and Heat Absorption effect on Natural Convection Flow with Uniform Surface Temperature along a Vertical Wavy Surface

Viscous Dissipation and Heat Absorption effect on Natural Convection Flow with Uniform Surface Temperature along a Vertical Wavy Surface Aailable at htt://am.ed/aam Al. Al. Math. ISSN: 93-966 Alications and Alied Mathematics: An International Jornal (AAM) Secial Isse No. (Ma 6),. 8 8th International Mathematics Conference, March,, IUB Cams,

More information

Q2. The velocity field in a fluid flow is given by

Q2. The velocity field in a fluid flow is given by Kinematics of Flid Q. Choose the correct anser (i) streamline is a line (a) hich is along the path of a particle (b) dran normal to the elocit ector at an point (c) sch that the streamlines diide the passage

More information

Non-Linear Squeezing Flow of Casson Fluid. between Parallel Plates

Non-Linear Squeezing Flow of Casson Fluid. between Parallel Plates International Jornal o Matematical Analsis Vol. 9 5 no. 5 - HIKAI Ltd www.m-ikari.com tt://d.doi.org/.988/ijma.5.49 Non-Linear Sqeezing Flow o Casson Flid between Parallel Plates S. Ganes C. K. Kirbasankar

More information

Computer Animation. Rick Parent

Computer Animation. Rick Parent Algorithms and Techniqes Flids Sperficial models. Deep models comes p throghot graphics, bt particlarl releant here OR Directl model isible properties Water waes Wrinkles in skin and cloth Hi Hair Clods

More information

The Vorticity Equation

The Vorticity Equation The Vorticit Eqation Potential orticit Circlation theorem is reall good Circlation theorem imlies a consered qantit dp dt 0 P g 2 PV or barotroic lid General orm o Ertel s otential orticit: P g const Consider

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Kinetic derivation of a finite difference scheme for the incompressible Navier Stokes equation

Kinetic derivation of a finite difference scheme for the incompressible Navier Stokes equation deriation of a finite difference scheme for the incompressible Naier Stokes eqation Mapndi K. Banda Michael Jnk Axel Klar Abstract In the present paper the low Mach nmber limit of kinetic eqations is sed

More information

5 th Smosium on Integrating CF and Exeriments in Aerodnamics (Integration 0) 3-5 October 0 JAXA Chofu Aerosace Center, Toko, Jaan rag and Lift Predict

5 th Smosium on Integrating CF and Exeriments in Aerodnamics (Integration 0) 3-5 October 0 JAXA Chofu Aerosace Center, Toko, Jaan rag and Lift Predict 5th Smosium on Integrating CF and Exeriments in Aerodnamics (Integration 0) 45 5 th Smosium on Integrating CF and Exeriments in Aerodnamics (Integration 0) 3-5 October 0 JAXA Chofu Aerosace Center, Toko,

More information

dx dt V x V t V y a Dt Acceleration field z dz dt V dt V v y V u t a Dt

dx dt V x V t V y a Dt Acceleration field z dz dt V dt V v y V u t a Dt FUNDAMENTALS OF Chate 6 Flo Analsis FLUID MECHANICS Using Diffeential Methods MAIN TOPICS I. Flid Element Motion II. Conseation of Mass (Continit eqation) and Conseation of Linea Momentm (Naie-Stokes Eqation)

More information

SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS

SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS THERMAL SCIENCE, Year 011, Vol. 15, No. 1, pp. 45-55 45 SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS b Khdheer S. MUSHATET Mechanical Engineering Department,

More information

EFFECT OF DIMENSIONLESS NUMBER ON NATURAL CONVECTION IN AN OPEN SQUARE CAVITY HAVING PARTIALLY HEATED SQUARE CYLINDER

EFFECT OF DIMENSIONLESS NUMBER ON NATURAL CONVECTION IN AN OPEN SQUARE CAVITY HAVING PARTIALLY HEATED SQUARE CYLINDER Proceedings of the International Conference on echanical Engineering ICE) 8- December Dhaka Bangladesh ICE-TH- EFFECT OF DIEIOE BER O TR COECTIO I OPE QRE CIT HIG PRTI HETED QRE CIDER ohammed asir ddin

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

STATIC, STAGNATION, AND DYNAMIC PRESSURES

STATIC, STAGNATION, AND DYNAMIC PRESSURES STATIC, STAGNATION, AND DYNAMIC PRESSURES Bernolli eqation is g constant In this eqation is called static ressre, becase it is the ressre that wold be measred by an instrment that is static with resect

More information

1 Differential Equations for Solid Mechanics

1 Differential Equations for Solid Mechanics 1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to

More information

Fluid Physics 8.292J/12.330J

Fluid Physics 8.292J/12.330J Fluid Phsics 8.292J/12.0J Problem Set 4 Solutions 1. Consider the problem of a two-dimensional (infinitel long) airplane wing traeling in the negatie x direction at a speed c through an Euler fluid. In

More information

Turbulence and boundary layers

Turbulence and boundary layers Trblence and bondary layers Weather and trblence Big whorls hae little whorls which feed on the elocity; and little whorls hae lesser whorls and so on to iscosity Lewis Fry Richardson Momentm eqations

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

Numerical Simulation of Density Currents over a Slope under the Condition of Cooling Period in Lake Biwa

Numerical Simulation of Density Currents over a Slope under the Condition of Cooling Period in Lake Biwa Nmerical Simlation of Densit Crrents oer a Slope nder the Condition of Cooling Period in Lake Bia Takashi Hosoda Professor, Department of Urban Management, Koto Uniersit, C1-3-65, Kotodai-Katsra, Nishiko-k,

More information

FORCED CONVECTIVE HEAT TRANSFER ENHANCEMENT WITH PERFORATED PIN FINS SUBJECT TO AN IMPINGING FLOW ABSTRACT

FORCED CONVECTIVE HEAT TRANSFER ENHANCEMENT WITH PERFORATED PIN FINS SUBJECT TO AN IMPINGING FLOW ABSTRACT SEGi Reie ISSN 1985-567 Vol. 5, No. 1, Jl 01, 9-40 *Corresponding athor. E-mail: jjfoo@segi.ed.m FORCED CONVECIVE HEA RANSFER ENHANCEMEN WIH PERFORAED PIN FINS SUBJEC O AN IMPINGING FLOW *Ji-Jinn Foo 1,,

More information

MEG 741 Energy and Variational Methods in Mechanics I

MEG 741 Energy and Variational Methods in Mechanics I MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

Dynamics of the Atmosphere 11:670:324. Class Time: Tuesdays and Fridays 9:15-10:35

Dynamics of the Atmosphere 11:670:324. Class Time: Tuesdays and Fridays 9:15-10:35 Dnamics o the Atmosphere 11:67:34 Class Time: Tesdas and Fridas 9:15-1:35 Instrctors: Dr. Anthon J. Broccoli (ENR 9) broccoli@ensci.rtgers.ed 73-93-98 6 Dr. Benjamin Lintner (ENR 5) lintner@ensci.rtgers.ed

More information

Ch.1: Basics of Shallow Water Fluid

Ch.1: Basics of Shallow Water Fluid AOS611Chapter1,/16/16,Z.Li 1 Sec. 1.1: Basic Eqations 1. Shallow Water Eqations on a Sphere Ch.1: Basics of Shallow Water Flid We start with the shallow water flid of a homogeneos densit and focs on the

More information

Boundary layer develops in the flow direction, δ = δ (x) τ

Boundary layer develops in the flow direction, δ = δ (x) τ 58:68 Trblent Flos Handot: Bondar Laers Differences to Trblent Channel Flo Bondar laer develops in the flo direction, not knon a priori Oter part of the flo consists of interittent trblent/non-trblent

More information

Correlations for Nusselt Number in Free Convection from an Isothermal Inclined Square Plate by a Numerical Simulation

Correlations for Nusselt Number in Free Convection from an Isothermal Inclined Square Plate by a Numerical Simulation American Jornal of Mechanics and Applications 5; : 8-8 Pblished online Ma 4, 5 http://sciencepblishinggropcom/j/ajma doi: 648/jajma5 ISSN: 76-65 Print; ISSN: 76-6 Online Correlations for Nsselt Nmber in

More information

FLUID MECHANICS. 1. Division of Fluid Mechanics. Hydrostatics Aerostatics Hydrodynamics Gasdynamics. v velocity p pressure ρ density

FLUID MECHANICS. 1. Division of Fluid Mechanics. Hydrostatics Aerostatics Hydrodynamics Gasdynamics. v velocity p pressure ρ density FLUID MECHANICS. Diision of Fluid Mechanics elocit p pressure densit Hdrostatics Aerostatics Hdrodnamics asdnamics. Properties of fluids Comparison of solid substances and fluids solid fluid τ F A [Pa]

More information

BOUNDARY LAYER FLOW OF A VISCOELASTIC LIQUID NEAR A STAGNATION POINT WITH SLIP

BOUNDARY LAYER FLOW OF A VISCOELASTIC LIQUID NEAR A STAGNATION POINT WITH SLIP BOUNDARY LAYER FLOW OF A VISCOELASTIC LIQUID NEAR A STAGNATION POINT WITH SLIP ACSahoo, TBiswal Det of Mathematics, TemleCit Institte of Technolog & Engineering (TITE), F/, IID Center, Barnei, Khrda Odisha

More information

An example of Lagrangian for a non-holonomic system

An example of Lagrangian for a non-holonomic system Uniersit of North Georgia Nighthaks Open Institutional Repositor Facult Publications Department of Mathematics 9-9-05 An eample of Lagrangian for a non-holonomic sstem Piotr W. Hebda Uniersit of North

More information

called the potential flow, and function φ is called the velocity potential.

called the potential flow, and function φ is called the velocity potential. J. Szantr Lectre No. 3 Potential flows 1 If the flid flow is irrotational, i.e. everwhere or almost everwhere in the field of flow there is rot 0 it means that there eists a scalar fnction ϕ,, z), sch

More information

ME 321: FLUID MECHANICS-I

ME 321: FLUID MECHANICS-I 8/7/18 ME 31: FLUID MECHANICS-I Dr. A.B.M. Toiqe Hasan Proessor Dearmen o Mechanical Engineering Bangladesh Uniersi o Engineering & Technolog BUET, Dhaka Lecre-13 8/7/18 Dierenial Analsis o Flid Moion

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

Homotopy Perturbation Method for Solving Linear Boundary Value Problems International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop

More information

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani

More information

SMS-618, Particle Dynamics, Fall 2003 (E. Boss, last updated: 10/8/2003) Conservation equations in fluids

SMS-618, Particle Dynamics, Fall 2003 (E. Boss, last updated: 10/8/2003) Conservation equations in fluids SMS-68 Parcle Dnamcs Fall 3 (E. Boss las daed: /8/3) onseraon eqaons n flds onces e need: ensor (Sress) ecors (e.g. oson eloc) and scalars (e.g. S O). Prode means o descrbe conseraon las h comac noaon

More information

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory

Derivation of 2D Power-Law Velocity Distribution Using Entropy Theory Entrop 3, 5, -3; doi:.339/e54 Article OPEN ACCESS entrop ISSN 99-43 www.mdpi.com/jornal/entrop Deriation of D Power-Law Velocit Distribtion Using Entrop Theor Vija P. Singh,, *, stao Marini 3 and Nicola

More information

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation Tye of Flows Continity Eqation Bernolli Eqation Steady Flow Energy Eqation Alications of Bernolli Eqation Flid Dynamics Streamlines Lines having the direction of the flid velocity Flids cannot cross a

More information

Comments on Vertical Vorticity Advection

Comments on Vertical Vorticity Advection Comments on Vertical Vorticity Advection It shold be fairly intitive that ositive maima in vertical vorticity are associated with cyclones, and ths ositive cyclonic vorticity advection might be a sefl

More information

Evaluation of the Lattice-Boltzmann Equation Solver PowerFLOW for Aerodynamic Applications

Evaluation of the Lattice-Boltzmann Equation Solver PowerFLOW for Aerodynamic Applications NASA/CR-2-255 ICASE Report No. 2-4 Ealation of the Lattice-Boltzmann Eqation Soler for Aerodnamic Applications Daid P. Lockard NASA Langle Research Center, Hampton, Virginia Li-Shi Lo ICASE, Hampton, Virginia

More information

The Faraday Induction Law and Field Transformations in Special Relativity

The Faraday Induction Law and Field Transformations in Special Relativity Apeiron, ol. 10, No., April 003 118 The Farada Indction Law and Field Transformations in Special Relatiit Aleander L. Kholmetskii Department of Phsics, elars State Uniersit, 4, F. Skorina Aene, 0080 Minsk

More information

Chapter 5 Darcy s Law and Applications

Chapter 5 Darcy s Law and Applications Chater 5 Darcy s La and Alications 4.1 Introdction Darcy' s la q reservoir Note : Time scale is added dn dt 4. Darcy s La; Flid Potential K h 1 l h K h l Different Sand Pac Different K The ressre at any

More information

The Numerical Simulation of Enhanced Heat Transfer Tubes

The Numerical Simulation of Enhanced Heat Transfer Tubes Aailable online at.sciencedirect.com Phsics Procedia 4 (01 70 79 01 International Conference on Applied Phsics and Indstrial Engineering The Nmerical Simlation of Enhanced Heat Transfer Tbes Li Xiaoan,

More information

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS

DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS Forth International Conference on CFD in the Oil and Gas, Metallrgical & Process Indstries SINTEF / NTNU Trondheim, Noray 6-8 Jne 005 DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS John MORUD 1 1 SINTEF

More information

Wall treatment in Large Eddy Simulation

Wall treatment in Large Eddy Simulation Wall treatment in arge Edd Simlation David Monfort Sofiane Benhamadoche (ED R&D) Pierre Sagat (Université Pierre et Marie Crie) 9 novembre 007 Code_Satrne User Meeting Wall treatment in arge Edd Simlation

More information

ECE Notes 4 Functions of a Complex Variable as Mappings. Fall 2017 David R. Jackson. Notes are adapted from D. R. Wilton, Dept.

ECE Notes 4 Functions of a Complex Variable as Mappings. Fall 2017 David R. Jackson. Notes are adapted from D. R. Wilton, Dept. ECE 638 Fall 017 Daid R. Jackson Notes 4 Fnctions of a Comple Variable as Mappings Notes are adapted from D. R. Wilton, Dept. of ECE 1 A Fnction of a Comple Variable as a Mapping A fnction of a comple

More information

Numerical investigation of natural convection of air in vertical divergent channels

Numerical investigation of natural convection of air in vertical divergent channels Adanced Comptational Methods in Heat ransfer X 13 Nmerical inestigation of natral conection of air in ertical diergent channels O. Manca, S. Nardini, D. Ricci & S. ambrrino Dipartimento di Ingegneria Aerospaziale

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

On the importance of horizontal turbulent transport in high resolution mesoscale simulations over cities. A. Martilli (CIEMAT, Spain), B. R.

On the importance of horizontal turbulent transport in high resolution mesoscale simulations over cities. A. Martilli (CIEMAT, Spain), B. R. On the importance of horizontal trblent transport in high resoltion mesoscale simlations over cities. A. Martilli (CIEMAT, Spain), B. R. Rotnno, P. Sllivan, E. G. Patton, M. LeMone (NCAR, USA) In an rban

More information

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element

More information

3. Several Random Variables

3. Several Random Variables . Several Random Variables. To Random Variables. Conditional Probabilit--Revisited. Statistical Independence.4 Correlation beteen Random Variables Standardied (or ero mean normalied) random variables.5

More information

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow

More information

The New (2+1)-Dimensional Integrable Coupling of the KdV Equation: Auto-Bäcklund Transformation and Non-Travelling Wave Profiles

The New (2+1)-Dimensional Integrable Coupling of the KdV Equation: Auto-Bäcklund Transformation and Non-Travelling Wave Profiles MM Research Preprints, 36 3 MMRC, AMSS, Academia Sinica No. 3, December The New (+)-Dimensional Integrable Copling of the KdV Eqation: Ato-Bäcklnd Transformation and Non-Traelling Wae Profiles Zhena Yan

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

Circular Motion Act. Centripetal Acceleration and. SPH3UW: Circular Motion, Pg 1 -> SPH3UW: Circular Motion, Pg 2. Page 1. Uniform Circular Motion

Circular Motion Act. Centripetal Acceleration and. SPH3UW: Circular Motion, Pg 1 -> SPH3UW: Circular Motion, Pg 2. Page 1. Uniform Circular Motion SPH3UW Centripetal Acceleration and Circular Motion Uniform Circular Motion What does it mean? How do we describe it? What can we learn about it? SPH3UW: Circular Motion, Pg 1 SPH3UW: Circular Motion,

More information

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0 Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.

More information

Introduction to Convection

Introduction to Convection Chater 6 Introdction to Conection he Conection Bondar aer Velocit Bondar aer A coneqence o ico eect aociated with relatie motion between a lid and a race A region between the race and the ree tream whoe

More information

3 2D Elastostatic Problems in Cartesian Coordinates

3 2D Elastostatic Problems in Cartesian Coordinates D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates

More information

Lecture 9: 3.4 The Geometry of Linear Systems

Lecture 9: 3.4 The Geometry of Linear Systems Lectre 9: 3.4 The Geometry of Linear Systems Wei-Ta Ch 200/0/5 Dot Prodct Form of a Linear System Recall that a linear eqation has the form a x +a 2 x 2 + +a n x n = b (a,a 2,, a n not all zero) The corresponding

More information

Rossby waves (waves in vorticity)

Rossby waves (waves in vorticity) Rossb waes (waes in orticit) Stationar (toograhicall orced) waes NCEP Reanalsis Z500 Janar mean 2 3 Vorticit eqation z w z w z w t 2 1 Change in relatie (ertical comonent o) orticit at a oint, Adection

More information

Effect of Applied Magnetic Field on Pulsatile Flow of Blood in a Porous Channel

Effect of Applied Magnetic Field on Pulsatile Flow of Blood in a Porous Channel Dsmanta Kmar St et al, Int. J. Comp. Tec. Appl., ol 6, 779-785 ISSN:9-69 Effect of Applied Magnetic Field on Plsatile Flow of Blood in a Poros Cannel Sarfraz Amed Dsmanta Kmar St Dept. of Matematics, Jorat

More information

FRONT TRACKING FOR A MODEL OF IMMISCIBLE GAS FLOW WITH LARGE DATA

FRONT TRACKING FOR A MODEL OF IMMISCIBLE GAS FLOW WITH LARGE DATA FONT TACKING FO A MODEL OF IMMISCIBLE GAS FLOW WITH LAGE DATA HELGE HOLDEN, NILS HENIK ISEBO, AND HILDE SANDE Abstract. In this aer we stdy front tracking for a model of one dimensional, immiscible flow

More information

1. THE MOMENTUM EQUATIONS FOR SYNOPTIC-SCALE FLOW IN THE ROTATING COORDINATE SYSTEM

1. THE MOMENTUM EQUATIONS FOR SYNOPTIC-SCALE FLOW IN THE ROTATING COORDINATE SYSTEM NOTES FO THE THEOY OF WKD 35. THE MOMENTUM EQUATIONS FO SYNOPTIC-SCALE FLOW IN THE OTATING COODINATE SYSTEM Scalin o the momentm eqations or snotic scale circlation (>000km dimension) reslted in the elimination

More information

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile

More information

Baroclinic Buoyancy-Inertia Joint Stability Parameter

Baroclinic Buoyancy-Inertia Joint Stability Parameter Jornal of Oceanography, Vol. 6, pp. 35 to 46, 5 Baroclinic Boyancy-Inertia Joint Stability Parameter HIDEO KAWAI* 3-8 Shibagahara, Kse, Joyo, Kyoto Pref. 6-, Japan (Receied 9 September 3; in reised form

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

MODELING OF POLLUTION DISPERSION AROUND A CUBIC OBSTACLE

MODELING OF POLLUTION DISPERSION AROUND A CUBIC OBSTACLE BBAA VI International Colloqim on: Blff Bodies Aerodnamics & Applications Milano, Ital, Jl, 0-4 008 MODELING OF POLLUTION DISPERSION AROUND A CUBIC OBSTACLE Department of Chemical and Petrolem Engineering,

More information

Lecture 17 Errors in Matlab s Turbulence PSD and Shaping Filter Expressions

Lecture 17 Errors in Matlab s Turbulence PSD and Shaping Filter Expressions Lectre 7 Errors in Matlab s Trblence PSD and Shaping Filter Expressions b Peter J Sherman /7/7 [prepared for AERE 355 class] In this brief note we will show that the trblence power spectral densities (psds)

More information

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7 754 Chapter 6 Additional Topics in Trigonometry 115. Yo ant to fly yor small plane de north, bt there is a 75-kilometer ind bloing from est to east. a. Find the direction angle for here yo shold head the

More information

Inertial and gravitational mass in relativistic mechanics

Inertial and gravitational mass in relativistic mechanics Träge nd schere Masse in der Relativitätsmechanik Ann. Phs. (Leipig) (4) 4 (1913) 856-878. Inertial and gravitational mass in relativistic mechanics B Gnnar Nordstrøm Translated b D. H. Delphenich In several

More information

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently Relatiity II I. Henri Poincare's Relatiity Principle In the late 1800's, Henri Poincare proposed that the principle of Galilean relatiity be expanded to inclde all physical phenomena and not jst mechanics.

More information

Module 6. Lecture 2: Navier-Stokes and Saint Venant equations

Module 6. Lecture 2: Navier-Stokes and Saint Venant equations Modle 6 Lectre : Naer-Stokes and Sant Venant eqatons Modle 6 Naer-Stokes Eqatons Clade-Los Naer Sr George Gabrel Stokes St.Venant eqatons are dered from Naer-Stokes Eqatons for shallo ater flo condtons.

More information

2. Fluid-Flow Equations

2. Fluid-Flow Equations . Fli-Flo Eqation Governing Eqation Conervation eqation for: ma momentm energy (other contitent) Alternative form: integral (control-volme) eqation ifferential eqation Integral (Control-olme) Approach

More information

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n 12.4 The Cross Prodct 873 12.4 The Cross Prodct In stdying lines in the plane, when we needed to describe how a line was tilting, we sed the notions of slope and angle of inclination. In space, we want

More information

Parametrization of turbulent fluxes in the outer layer

Parametrization of turbulent fluxes in the outer layer Parametriation of trblent flxes in the oter layer Irina Sand Oerie of models Blk models Local -closre -profile closre 0 order 1 st order ED/F closre non-local TE closre 1.5 th order Crrent closre in the

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

GEF2500 GEOPHYSICAL FLUID MECHANICS

GEF2500 GEOPHYSICAL FLUID MECHANICS GEF5 GEOPHYSICAL FLUID MECHANICS Jan Erik H. Weber Deartment of Geosciences Section for Meteorolog and Oceanograh Uniersit of Oslo. E-mail: j.e.weber@geo.uio.no Blindern Januar 3 Contents. FLUID MECHANICS

More information

Lesson 81: The Cross Product of Vectors

Lesson 81: The Cross Product of Vectors Lesson 8: The Cross Prodct of Vectors IBHL - SANTOWSKI In this lesson yo will learn how to find the cross prodct of two ectors how to find an orthogonal ector to a plane defined by two ectors how to find

More information

Change of Variables. (f T) JT. f = U

Change of Variables. (f T) JT. f = U Change of Variables 4-5-8 The change of ariables formla for mltiple integrals is like -sbstittion for single-ariable integrals. I ll gie the general change of ariables formla first, and consider specific

More information

Spring Semester 2011 April 5, 2011

Spring Semester 2011 April 5, 2011 METR 130: Lectre 4 - Reynolds Averaged Conservation Eqations - Trblent Flxes (Definition and typical ABL profiles, CBL and SBL) - Trblence Closre Problem & Parameterization Spring Semester 011 April 5,

More information

Momentum and Energy. Relativity and Astrophysics Lecture 24 Terry Herter. Energy and Momentum Conservation of energy and momentum

Momentum and Energy. Relativity and Astrophysics Lecture 24 Terry Herter. Energy and Momentum Conservation of energy and momentum Momentum and Energy Relatiity and Astrohysics Lecture 4 Terry Herter Outline Newtonian Physics Energy and Momentum Conseration of energy and momentum Reading Sacetime Physics: Chater 7 Homework: (due Wed.

More information

CHAPTER 8 CONVECTION IN EXTERNAL TURBULENT FLOW

CHAPTER 8 CONVECTION IN EXTERNAL TURBULENT FLOW CHAPTER 8 CONVECTION IN EXTERNAL TURBULENT FLOW 8.1 Introdction Common phsical phenomenon, bt comple Still relies on empirical data and rdimentar conceptal drawings Tremendos growth in research over last

More information

6.4 VECTORS AND DOT PRODUCTS

6.4 VECTORS AND DOT PRODUCTS 458 Chapter 6 Additional Topics in Trigonometry 6.4 VECTORS AND DOT PRODUCTS What yo shold learn ind the dot prodct of two ectors and se the properties of the dot prodct. ind the angle between two ectors

More information

Boundary Layer Theory:

Boundary Layer Theory: Mass Transfer Bondar Laer Theor: Mass and Heat/Momentm Transfer Letre,..7, Dr. K. Wegner 9. Basi Theories for Mass Transfer Coeffiients 9. Flid-Flid Interfaes (letre of 5..7) Flid-flid interfaes are tpiall

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

MATH2715: Statistical Methods

MATH2715: Statistical Methods MATH275: Statistical Methods Exercises III (based on lectres 5-6, work week 4, hand in lectre Mon 23 Oct) ALL qestions cont towards the continos assessment for this modle. Q. If X has a niform distribtion

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

UNIT V BOUNDARY LAYER INTRODUCTION

UNIT V BOUNDARY LAYER INTRODUCTION UNIT V BOUNDARY LAYER INTRODUCTION The variation of velocity from zero to free-stream velocity in the direction normal to the bondary takes place in a narrow region in the vicinity of solid bondary. This

More information

INVESTIGATION OF THERMAL PERFORMANCE OF AIR TO WATER HEAT EXCHANGER USING NANO- FLUIDS

INVESTIGATION OF THERMAL PERFORMANCE OF AIR TO WATER HEAT EXCHANGER USING NANO- FLUIDS IIUM Engineering Jornal, Special Isse, Mechanical Engineering, 11 INVESIGAION OF HERMAL PERFORMANCE OF AIR O WAER HEA EXCHANGER USING NANO- FLUIDS NAWAF H. SAEID AND AN HENG CHIA Department of Mechanical,

More information

Propagation of error for multivariable function

Propagation of error for multivariable function Proagation o error or mltiariable nction No consider a mltiariable nction (,,, ). I measrements o,,,. All hae ncertaint,,,., ho ill this aect the ncertaint o the nction? L tet) o (Eqation (3.8) ± L ),...,,

More information

Dimensional analysis. Scaling and similitude

Dimensional analysis. Scaling and similitude Dimensional analsis Scaling - a powerfl idea Similitde Bckingham Pi theorem Eamples of the power of dimensional analsis Usefl dimensionless qantities and their interpretation Scaling and similitde Scaling

More information

Conservation of Energy Thermodynamic Energy Equation

Conservation of Energy Thermodynamic Energy Equation Conseration of Energy Thermodynamic Energy Equation The reious two sections dealt with conseration of momentum (equations of motion) and the conseration of mass (continuity equation). This section addresses

More information

2 Faculty of Mechanics and Mathematics, Moscow State University.

2 Faculty of Mechanics and Mathematics, Moscow State University. th World IMACS / MODSIM Congress, Cairns, Astralia 3-7 Jl 9 http://mssanz.org.a/modsim9 Nmerical eamination of competitie and predator behaior for the Lotka-Volterra eqations with diffsion based on the

More information

λ. It is usually positive; if it is zero then the constraint is not binding.

λ. It is usually positive; if it is zero then the constraint is not binding. hater 4 Utilit Maimization and hoice rational consmer chooses his most referred bndle of commodities from the set of feasible choices to consme The rocess of obtaining this otimal bndle is called tilit-maimization

More information

An alternative approach to evaluate the average Nusselt number for mixed boundary layer conditions in parallel flow over an isothermal flat plate

An alternative approach to evaluate the average Nusselt number for mixed boundary layer conditions in parallel flow over an isothermal flat plate An alternative approach to evalate the average Nsselt nber for ied bondary layer conditions in parallel flo over an isotheral flat plate Viacheslav Stetsyk, Krzysztof J. Kbiak, ande i and John C Chai Abstract

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas.

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas. Two identical flat sqare plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHAE areas. F > F A. A B F > F B. B A C. FA = FB. It depends on whether the bondary

More information