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

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

3 2D Elastostatic Problems in Cartesian Coordinates

1 Differential Equations for Solid Mechanics

STATIC, STAGNATION, AND DYNAMIC PRESSURES

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

Second-Order Wave Equation

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

Chapter 1: Differential Form of Basic Equations

Elements of Coordinate System Transformations

Method of Images

4 Exact laminar boundary layer solutions

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

AMS 212B Perturbation Methods Lecture 05 Copyright by Hongyun Wang, UCSC

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

5. The Bernoulli Equation

PhysicsAndMathsTutor.com

Axial Compressor Design Parameters

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

4 Primitive Equations

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

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

Formal Methods for Deriving Element Equations

THE EFFECTS OF RADIATION ON UNSTEADY MHD CONVECTIVE HEAT TRANSFER PAST A SEMI-INFINITE VERTICAL POROUS MOVING SURFACE WITH VARIABLE SUCTION

Partial Differential Equations with Applications

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

NATURAL CONVECTION No mechanical force to push the fluid pump, fan etc. No predefined fluid flowrate and velocity can t prescribe Reynolds

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS

08.06 Shooting Method for Ordinary Differential Equations

Offshore Hydromechanics Module 1

UNIT V BOUNDARY LAYER INTRODUCTION

Theoretical Fluid Mechanics Turbulent Flow Velocity Profile By James C.Y. Guo, Professor and P.E. Civil Engineering, U. of Colorado at Denver

Elastico-Viscous MHD Free Convective Flow Past an Inclined Permeable Plate with Dufour Effects in Presence of Chemical Reaction

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

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Conceptual Questions. Problems. 852 CHAPTER 29 Magnetic Fields

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

OUTLINE FOR Chapter 3

Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of

STUDY OF THE NON-DIMENSIONAL SOLUTION OF DYNAMIC EQUATION OF MOVEMENT ON THE PLANE PLAQUE WITH CONSIDERATION OF TWO-ORDER SLIDING PHENOMENON

Section 7.4: Integration of Rational Functions by Partial Fractions

Determining of temperature field in a L-shaped domain

FEA Solution Procedure

Chapter 6. Inverse Circular Functions and Trigonometric Equations. Section 6.1 Inverse Circular Functions y = 0

ECON3120/4120 Mathematics 2, spring 2009

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

EXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.

3.4-Miscellaneous Equations

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

Inertial and gravitational mass in relativistic mechanics

Mean Value Formulae for Laplace and Heat Equation

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

L = 2 λ 2 = λ (1) In other words, the wavelength of the wave in question equals to the string length,

Turbomachinery Lecture Notes

THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.

Motion in an Undulator

PHY 113 C General Physics I 11 AM 12:15 PM TR Olin 101


Calculations involving a single random variable (SRV)

MATH2715: Statistical Methods


Digital Image Processing. Lecture 8 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep.

Boundary Layer Flow and Heat Transfer over a. Continuous Surface in the Presence of. Hydromagnetic Field

Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction

Wall treatment in Large Eddy Simulation

Introducing Ideal Flow

Fluidmechanical Damping Analysis of Resonant Micromirrors with Out-of-plane Comb Drive

Prandl established a universal velocity profile for flow parallel to the bed given by

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63

Dimensional analysis. Scaling and similitude

Material Transport with Air Jet

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

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

CFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion

Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature

Chapter 2 Difficulties associated with corners

Created by T. Madas VECTOR OPERATORS. Created by T. Madas

TRANSIENT FREE CONVECTION MHD FLOW BETWEEN TWO LONG VERTICAL PARALLEL PLATES WITH VARIABLE TEMPERATURE AND UNIFORM MASS DIFFUSION IN A POROUS MEDIUM

Problems set # 2 Physics 169 February 11, 2015

Fluid Mechanics II. Newton s second law applied to a control volume

Spring Semester 2011 April 5, 2011

Curves - Foundation of Free-form Surfaces

Polymer confined between two surfaces

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

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

Finite Element Analysis of Heat and Mass Transfer of a MHD / Micropolar fluid over a Vertical Channel

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions

INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1

Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

2. Find the coordinates of the point where the line tangent to the parabola 2

3. Several Random Variables

arxiv: v1 [physics.flu-dyn] 11 Mar 2011

Transcription:

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 that gradϕ. Sch a flow is called the potential flow, and fnction φ is called the velocit potential. ϕ ϕ We have: z ϕ z In case of potential flow of an incompressible flid the mass conservation eqation is transformed into the Laplace eqation: ρ t div ϕ ϕ ϕ z ρ ) divgradϕ ϕ 0 0 The Laplace eqation is linear, that means that the sm of its soltions is also a soltion. In practise, ver complicated potential fnctions, describing comple flows, ma be composed of the simple fnctions describing so called elementar flows.

The potential flows are particlarl sited for mathematical modelling of the flid motion in the regions otside the bondar laers and wakes, where the inflence of the flid viscosit is negligibl small. The method of formation of the complicated potential flows will be demonstrated sing the eample of twodimensional flows. In this case we have: where: ϕ, ) ϕ, ) ψ ψ, ), ) ϕ ψ - the velocit potential - the eqipotential lines - the stream fnction - the stream lines ϕ ψ

The elementar potential flows 1. The niform flow b ϕ a ϕ ) b a, ϕ The eqipotential lines: The velocit potential: b a The stream fnction: The stream lines: a b ) b a, ψ

The constant vales of the potential φ correspond to the constant vales of the radis r, i.e. the eqipotential lines are circles having a common centre.. The sorce positive or negative) ±πr r r r The sorce is a singlar point in the field of flow, in which an otflow of flid of a certain volmetric intensit takes place. This otflow is niform in all directions. In case of a negative sorce or a sink), the flid flows into the sink and disappears. Ths we have: ± r πr or: where: - the radial velocit ϕ r ± ϕ ± πr π ln r

In the artesian sstem of co-ordinates we have: r cosθ π and frther we have: ) the complete differential of the potential: dϕ where: θ r sinθ arctg π d d π ) ) π ) the complete differential of the stream fnction: dψ d d π ) π ) and frther, after integration we have: ϕ 1 π ln ) ψ π arctg Stream lines are straight lines passing throgh the sorce

Eample: Sperposition of an niform flow and a sorce The potential and stream fnction of the resltant flow are the sms of potentials and stream fnctions We assme that the niform flow is parallel to the ais. The potential: The stream fnction: ϕ ln 4π ψ arctg π ) The zero-th stream line: ψ 0 arctg 0

Soltion for the zero-th stream line: ctg π If we replace the zero-th stream line with a solid wall, the pictre of flow will not change and we obtain the flow arond a solid halfbod. The velocit components: π π ) ) The location of the stagnation point on the ais: 0 π The pressre p in an arbitrar point of the flow ma be calclated according to the Bernolli eqation: p ρ 1 1 ρ p p p p ρ

After sbsittion we get: p 1 π 4 ) π In the stagnation point we get: p π π π 1 Bt on the zero-th stream line at the point 0 we have: ± π or: 4 184 π 1,

3. The doblet the dipole) The dipole reslts from sperposition of a positive and negative sorce of the same absolte intensit. The measre of the dipole intensit is the so called moment Ma. ontrar to a sorce, a dipole has directional characteristics, becase it emits flid in a given direction and scks in flid from the opposite direction. Hence its location in space is important. For a dipole at 0, 0, directed in the positive direction, we have: The dipole potential: ϕ M π The dipole stream fnction: ψ M π

4 The stream lines for a dipole: The stream lines are circles of radii /, the centres of which lie on the ais in the points /. 4 The eqipotential lines for a dipole: The eqipotential lines are circles of radii /, the centres of which lie on the ais in the points /.

Eample: Irrotational flow arond a circlar clinder Sperposition: the niform flow the dipole

In the sstem Orθ we have: M πa where: a clinder radis a a 1 r r 4 cos θ where: θ arctg ρ p p 1 for ra i.e. on the clinder srface we have: p θ ρ p 1 4sin θ ) The components of the force acting on the clinder ma be calclated: P P a a π pθ cosθdθ 0 0 π pθ sinθdθ 0 0 This is so called d Alembert parado

The d Alembert parado means, that in a potential flow the forces eerted b the flowing flid on the immersed bodies are eqal zero, what does not agree with a common eperience. This is a direct conseqence of the smmetr of the calclated pressre field, which is smmetrical both with respect to the ais and ais. In realit, the pressre field on the clinder is asmmetrical with respect to the ais, what is shown in the diagram on the previos slide and in the sketch and photograph below, showing the calclated and eperimentall visalised flows. The differences are first of all visible on the trailing side of the clinder.