An Exact Solution of Navier Stokes Equation
|
|
- Ezra Miller
- 6 years ago
- Views:
Transcription
1 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 solving the Navie Stokes equations (a set of nonlinea patial diffeential equations) aises fom the pesence of the nonlinea convective tem (V )V. Since thee ae no geneal analytical methods fo solving nonlinea patial diffeential equations exist, each poblem must be consideed individually. Fo most pactical flow poblems, convective acceleation of fluid paticles cannot be ignoed. Howeve, in geneal, exact solutions ae possible only when the nonlinea tems vanishes identically. Thee ae a few special cases fo which the convective acceleation vanish because of the natue of the geomety of the flow system. In these cases exact solutions ae usually possible and below we conside one of such poblems. Pipe Flow Induced by Movement of Wall Hee we conside an infinitely long hoizontal cicula pipe filled with a Newtonian fluid of density ρ and viscosity µ. Fo t < 0 both pipe and fluid ae at est. At time t = 0, the pipe impulsively stats moving in the axial diection with a unifom velocity U. As a consequence of this, the fluid movement is induced in the axial diection as sketched in figue. The momentum equation fo incompessible flow in cylindical coodinate system Figue : Fluid flow in a pipe induced by the motion of pipe walls.
2 is given by ( Du ρ Dt ( Duθ ρ Dt u2 θ + u u θ ) ρ Du z Dt = ρg p ( + µ ) = ρg θ p θ + µ = ρg z p z + µ 2 u z 2 u 2 u θ 2 θ u ) ( 2 2 u θ + 2 u 2 θ u ) θ 2 Fo the conditions of the pesent unsteady paallel flows in the absence of body foces, many tems disappea and the axial velocity u z is the only nonzeo velocity component. Theefoe the fist two equations become tivial and the last equation (axial momentum equation) educes to ρ u z t [ = µ ( u ) z + 2 ] u z 2 θ 2 Since the flow is otationally symmetic, equation () educes to u z = ν ( u ) z t The initial and bounday conditions ae u z (, 0) = 0 (initial condition) u z (0, t) is bounded (3) u z (R, t) = U (no-slip condition) Fo convenience, we intoduce the following dimensionless quantities: u = U u z U η = R τ = νt R 2 Thus the diffeential equation (2) can be witten in dimensionless quantities as follows: ν U u R 2 τ = ν U [ 2 u R 2 η 2 + ] u η η u τ = 2 u η 2 + u η η Coesponding dimensionless initial and bounday conditions ae u(η, 0) = (initial condition) u(0, τ) is bounded (5) u(, τ) = 0 (no-slip condition) () (2) (4) 2
3 Note that the new bounday conditions ae homogeneous. Although the poblem defined by (4) and (5) is time dependent, it is linea in u and confined to the bounded spatial domain, 0 η. Thus it can be solved by the method of sepaation of vaiables. In this method we fist find a set of eigensolutions that satisfy the diffeential equation (5) and the bounday condition at η = 0 and η = ; then we detemine the paticula sum of those eigensolutions that also satisfies the initial condition at τ = 0. The poblem (4) and (5) compises one example of the geneal class of so-called Stum Liouville poblems fo which an extensive theoy is available that ensues the existence and uniqueness of solutions constucted by means of eigenfunction expansions by the method of sepaation of vaiables. We begin with the basic hypothesis that a solution of (5) exists in the sepaable fom and choose the following ansatz: u(η, τ) = F(η)G(τ) (6) Substituting this ansatz into equation (4) to obtain F dg dτ = F Gd2 dη 2 + G df η dη (7) As G depends only on τ and F only on η, by sepaation of vaiables the following odinay diffeential equations fo G and F esult: dg G dτ = λ 2 (8) d 2 F F dη 2 + df F η dη = λ 2 (9) Hence the oiginal poblem govened by a PDE has now been tansfomed into two elated auxiliay poblems govened by ODEs. The choice of a negative constant λ 2 is due to the fact that the solution will decay to zeo as time inceases, i.e., u 0 as τ. The solution fo the diffeential equation (8) can be deived by integation: G = c 0 e λ 2 τ (0) whee c 0 is an integation constant to be detemined. In ode to detemine the solution of the diffeential equation fo F(η), equation (9) can be witten as follows: d 2 F dη 2 + df η dη + λ 2 F = 0 o η 2 d2 F dη 2 + η df dη + λ 2 η 2 F = 0 () 3
4 o, on intoducing a change of independent vaiables, y = λη we obtain y 2 d2 F dy 2 + ydf dy + y2 F = 0 (2) This is Bessel diffeential equation of ode 0. It has two linealy independent solutions, J 0 (y) and Y 0 (y) which ae known as Bessel functions of the fist and second kinds of ode 0. Thus the most geneal solution of equation (2) can be witten as F(y) = c J 0 (y) + c 2 Y 0 (y) (3) A plot showing the behavio of these two functions is shown in figue (2). Both oscillate Figue 2: The Bessel functions of the fist and second kinds of ode 0. back and foth acoss zeo, but Y 0 (y) as y 0. Substituting equations (0) and (3) into (6), we obtain the geneal solution of the fom u(y, τ) = e λ 2τ [c J 0 (y) + c 2 Y 0 (y)] (4) Note that the constant c 0 has been dopped since it is edundant hee. The solution is bounded at y = 0. This condition is satisfied only when the constant c 2 = 0. Hence, the geneal solution (4) (that is bounded at y = 0) takes a fom: u(y, τ) = c J 0 (y)e λ 2 τ (5) 4
5 Theefoe, the solution u, afte substituting y = λη, may be witten as u n (η, τ) = A n J 0 (λ n η)e λ 2 n τ (6) The subscipt n is added in anticipation of the fact that thee is an infinite, but discete, set of values possible fo λ such that the geneal solution (5), satisfies the bounday condition u = 0 at η =. This set of values of λ = λ n is known as the eigenvalues of the poblem, and the coesponding u n ae called the eigenfunctions. To detemine the eigenvalues λ n, we apply the bounday condition at η = to equation (6), that is, 0 = A n J 0 (λ n )e λ 2 n τ fo all Since setting A n = 0 would esults in a tivial solution, one must equie τ J 0 (λ n ) = 0 (7) fo the non-tivial solution. Theefoe, one obtains multiple values of λ n (eigenvalues) that satisfy the bounday conditions at the wall. Clealy these eigenvalues ae equal to the infinite set of zeoes of the Bessel function of ode zeo, J 0 (z). Refeing to figue 2, we have denoted those zeoes as s n, with the fist cossing fo the smallest value of z being s, namely, and thei values ae obtained as λ n = s n, n =, 2, 3,..., (8) s n = 2.405, 5.520, 8.654,.792, 4.93, 8.07, 2.22, , (9) Each of the solutions of λ n now constitutes an individual solution. Consideing the lineaity of the govening equation and bounday conditions (4) and (5), the complete solution fo u n (η, τ) is obtained by linea supeposition: u = u n (η, τ) = A n e λ n 2τ J 0 (λ n η) (20) whee A n ae abitay, constant coefficients. Equation (20) is called the Fouie Bessel seies. This solution satisfies the diffeential equation (4) and the bounday condition u = 0 at η = fo any choice of the constant coefficients A n. The final step is to choose the A n so that u(η, τ) satisfies the initial condition u(η, 0) =. The geneal Stum Liouville theoy guaantees that the eigenfunctions (6) fom a complete set of othogonal functions. Thus it is possible to expess the smooth initial condition by means of the Fouie Bessel seies (20) with τ = 0, that is, u(η, 0) = A n J 0 (λ n η) = A J 0 (λ η) + A 2 J 0 (λ 2 η) +... A n J 0 (λ n η) +... = (2) 5
6 To detemine the A n, we will take advantage of othogonality popeties of J 0 : J 0 (λ m η)j 0 (λ n η)η dη = 2 J 2(λ n) if n = m 0 0 if n m whee J is the Bessel function of fist kind of ode. Multiplying both sides of (2) by ηj 0 (λ m η) and integating ove η fom 0 to, we obtain A n 0 η J 0 (λ m η)j 0 (λ n η)dη = 0 η J 0 (λ m η)dη Due to the othogonality popety (22), the only nonzeo tem on the left hand side is that fo m = n hence, 0 η J 0 (λ n η)dη 0 A n = 0 η J0 2(λ nη)dη = η J 0 (λ n η)dη 2 J2 (λ (23) n) Fo evaluating the numeato of (23) we make use of the following popety of Bessel function: Theefoe (22) η p+ J p (λη)dη = λ η p+ J p+ (λη) (24) η J 0 (λ n η)dη = λ n ηj (λ n η) (25) and thus, the numeato of equation (23) becomes η J 0 (λ n η)dη = J (λ n ) (26) λ n 0 Substitution of equation (26) into (23) yields an expession fo A n : λ A n = n J (λ n ) 2 J2 (λ n) = 2 [J (λ n )] λ n Thus fo the velocity distibution accoding to (20), the following expession esults: u(η, τ) = 2 (27) λ n J (λ n ) e λ 2 n τ J 0 (λ n η) (28) The above Fouie Bessel seies has the popety of conveging vey quickly when the dimensionless time τ = νt/r 2 is lage. On the othe hand, the convegence is slow when τ is small. Reveting to dimensional vaiables, we can expess the solution of the full, oiginal poblem in tems of the axial velocity pofile: [ u z (, t) = U 2 ( λ n J (λ n ) e λ n 2 νt/r 2 ) ] J 0 λ n (29) R Obviously, as t, this solution evets to the steady-state unifom flow pofile. obtain othe details of this velocity pofile, it is necessay to evaluate the infinite seies numeically fo each value of t and. A typical numeical example of the esults is shown in figue 3, whee u z has been plotted vesus fo seveal values of t. 6 To
7 Figue 3: Tansient velocity pofiles in a pipe induced by the motion of the pipe walls. Shea stess distibution The shea stess distibution in the flowfield can be easily obtained fom the velocity distibution as follows In dimensionless fom whee dimensionless shea stess τ x = ( u τ x = µ x + u ) z = µ du z d τ x = du dη τ x µu/r. pofile, we make use of the following popety of Bessel function: Theefoe and thus the shea stess pofile is given by Fo evaluating the deivative of velocity d [ η p J p (λη) ] = λη p J p+ (λη) (30) dη d dη [J 0(λη)] = λj (λη) (3) τ x = du dη = 2 [J (λ n )] e λ n 2τ J (λ n η) (32) Fo the value of wall shea stess at the pipe wall, i.e. at η =, one obtains τ wall = τ x η= = du dη = 2 e λ n 2 τ η= (33) Reveting to dimensional vaiables, we can expess the wall shea stess in tems of oiginal vaiables: τ wall = du z d = 2µU =R R e λ n 2 νt/r 2 (34) 7
8 It can be seen that the wall shea stess has a finite value, even fo time t = 0. This is a supising esult when compaing to othe cases whee flow is induced by impulsively motion of the bounday. Note also that the wall shea stess τ wall 0 as t. This is obvious because the steady-state velocity pofile is unifom. 8
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 informationDo 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 information8 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 information15 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 informationStress, 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 information1 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 informationDymore 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-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.
The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationd 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c
Chapte 6 Geneal Relativity 6.1 Towads the Einstein equations Thee ae seveal ways of motivating the Einstein equations. The most natual is pehaps though consideations involving the Equivalence Pinciple.
More information7.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 informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationApplied Aerodynamics
Applied Aeodynamics Def: Mach Numbe (M), M a atio of flow velocity to the speed of sound Compessibility Effects Def: eynolds Numbe (e), e ρ c µ atio of inetial foces to viscous foces iscous Effects If
More information2. 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=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 informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More informationComputational 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 informationA 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 informationEM 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 informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationMathematical 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 informationASTR415: 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 informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationELASTIC 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 informationA Most Useful Device of Studying Electrode Processes: The Rotating Disk Electrode
A Most Useful Device of Studying Electode Pocesses: The Rotating Disk Electode the theoetical basis Soma Vesztegom Laboatoy of Electochemisty & Electoanalytical Chemisty Eötvös Loánd Univesity of Budapest
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationClassical Mechanics Homework set 7, due Nov 8th: Solutions
Classical Mechanics Homewok set 7, due Nov 8th: Solutions 1. Do deivation 8.. It has been asked what effect does a total deivative as a function of q i, t have on the Hamiltonian. Thus, lets us begin with
More informationI. 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, 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 informationyou of a spring. The potential energy for a spring is given by the parabola U( x)
Small oscillations The theoy of small oscillations is an extemely impotant topic in mechanics. Conside a system that has a potential enegy diagam as below: U B C A x Thee ae thee points of stable equilibium,
More informationReview: 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 informationThree 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 informationA 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 informationLiquid gas interface under hydrostatic pressure
Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,
More informationINTRODUCTION. 2. Vectors in Physics 1
INTRODUCTION Vectos ae used in physics to extend the study of motion fom one dimension to two dimensions Vectos ae indispensable when a physical quantity has a diection associated with it As an example,
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationNumerical 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 informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationStress Intensity Factor
S 47 Factue Mechanics http://imechanicaog/node/7448 Zhigang Suo Stess Intensity Facto We have modeled a body by using the linea elastic theoy We have modeled a cack in the body by a flat plane, and the
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More information3. Electromagnetic Waves II
Lectue 3 - Electomagnetic Waves II 9 3. Electomagnetic Waves II Last time, we discussed the following. 1. The popagation of an EM wave though a macoscopic media: We discussed how the wave inteacts with
More informationA new class of exact solutions of the Navier Stokes equations for swirling flows in porous and rotating pipes
Advances in Fluid Mechanics VIII 67 A new class of exact solutions of the Navie Stokes equations fo swiling flows in poous and otating pipes A. Fatsis1, J. Stathaas2, A. Panoutsopoulou3 & N. Vlachakis1
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationAE301 Aerodynamics I UNIT B: Theory of Aerodynamics
AE301 Aeodynamics I UNIT B: Theoy of Aeodynamics ROAD MAP... B-1: Mathematics fo Aeodynamics B-2: Flow Field Repesentations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis AE301
More informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationPHYS 110B - HW #7 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 0B - HW #7 Sping 2004, Solutions by David Pace Any efeenced euations ae fom Giffiths Poblem statements ae paaphased. Poblem 0.3 fom Giffiths A point chage,, moves in a loop of adius a. At time t 0
More informationLINEAR PLATE BENDING
LINEAR PLATE BENDING 1 Linea plate bending A plate is a body of which the mateial is located in a small egion aound a suface in the thee-dimensional space. A special suface is the mid-plane. Measued fom
More informationA method for solving dynamic problems for cylindrical domains
Tansactions of NAS of Azebaijan, Issue Mechanics, 35 (7), 68-75 (016). Seies of Physical-Technical and Mathematical Sciences. A method fo solving dynamic poblems fo cylindical domains N.B. Rassoulova G.R.
More informationMagnetic Field. Conference 6. Physics 102 General Physics II
Physics 102 Confeence 6 Magnetic Field Confeence 6 Physics 102 Geneal Physics II Monday, Mach 3d, 2014 6.1 Quiz Poblem 6.1 Think about the magnetic field associated with an infinite, cuent caying wie.
More informationPhysics 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 informationIncompressible 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 information3-7 FLUIDS IN RIGID-BODY MOTION
3-7 FLUIDS IN IGID-BODY MOTION S-1 3-7 FLUIDS IN IGID-BODY MOTION We ae almost eady to bein studyin fluids in motion (statin in Chapte 4), but fist thee is one cateoy of fluid motion that can be studied
More informationf(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 informationPreliminary Exam: Quantum Physics 1/14/2011, 9:00-3:00
Peliminay Exam: Quantum Physics /4/ 9:-: Answe a total of SIX questions of which at least TWO ae fom section A and at least THREE ae fom section B Fo you answes you can use eithe the blue books o individual
More information12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationis 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 information6 PROBABILITY GENERATING FUNCTIONS
6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to
More informationDonnishJournals
DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş
More informationTransformation of the Navier-Stokes Equations in Curvilinear Coordinate Systems with Maple
Global Jounal of Pue and Applied Mathematics. ISSN 0973-1768 Volume 12, Numbe 4 2016, pp. 3315 3325 Reseach India Publications http://www.ipublication.com/gjpam.htm Tansfomation of the Navie-Stokes Equations
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationBoundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.
Lectues 16 and 17 Bounday Layes and Singula Petubation A Regula Petubation In some physical poblems, the solution is dependent on a paamete K. When the paamete K is vey small, it is natual to expect that
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More informationMath 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 informationScattering 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 informationCBE 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 informationChapter 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 informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationPROBLEM 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 informationAveraging 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 informationNumerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method
Numeical Invesion of the Abel Integal Equation using Homotopy Petubation Method Sunil Kuma and Om P Singh Depatment of Applied Mathematics Institute of Technology Banaas Hindu Univesity Vaanasi -15 India
More informationDiffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.
Diffusion and Tanspot 10. Fiction and the Langevin Equation Now let s elate the phenomena of ownian motion and diffusion to the concept of fiction, i.e., the esistance to movement that the paticle in the
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationLecture 10. Vertical coordinates General vertical coordinate
Lectue 10 Vetical coodinates We have exclusively used height as the vetical coodinate but thee ae altenative vetical coodinates in use in ocean models, most notably the teainfollowing coodinate models
More informationQuantum Mechanics II
Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e
More informationChapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in
ME 01 DYNAMICS Chapte 1: Kinematics of a Paticle Chapte 1 Kinematics of a Paticle A. Bazone 1.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS Pola Coodinates Pola coodinates ae paticlaly sitable fo solving
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More informationChapter 6 Differential Analysis of Fluid Flow
1 Chapte 6 Diffeential Analysis of Fluid Flow Inviscid flow: Eule s equations of otion Flow fields in which the sheaing stesses ae zeo ae said to be inviscid, nonviscous, o fictionless. fo fluids in which
More informationFOURIER-BESSEL SERIES AND BOUNDARY VALUE PROBLEMS IN CYLINDRICAL COORDINATES
FOURIER-BESSE SERIES AND BOUNDARY VAUE PROBEMS IN CYINDRICA COORDINATES The paametic Bessel s equation appeas in connection with the aplace opeato in pola coodinates. The method of sepaation of vaiables
More information1) 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 informationHydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods
TEAM 2007, Sept. 10-13, 2007,Yokohama, Japan Hydoelastic Analysis of a 1900 TEU Containe Ship Using Finite Element and Bounday Element Methods Ahmet Egin 1)*, Levent Kaydıhan 2) and Bahadı Uğulu 3) 1)
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More informationB da = 0. Q E da = ε. E da = E dv
lectomagnetic Theo Pof Ruiz, UNC Asheville, doctophs on YouTube Chapte Notes The Maxwell quations in Diffeential Fom 1 The Maxwell quations in Diffeential Fom We will now tansfom the integal fom of the
More informationarxiv: v1 [physics.pop-ph] 3 Jun 2013
A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,
More informationPHYS 301 HOMEWORK #10 (Optional HW)
PHYS 301 HOMEWORK #10 (Optional HW) 1. Conside the Legende diffeential equation : 1 - x 2 y'' - 2xy' + m m + 1 y = 0 Make the substitution x = cos q and show the Legende equation tansfoms into d 2 y 2
More informationPhysics 505 Homework No. 9 Solutions S9-1
Physics 505 Homewok No 9 s S9-1 1 As pomised, hee is the tick fo summing the matix elements fo the Stak effect fo the gound state of the hydogen atom Recall, we need to calculate the coection to the gound
More informationarxiv: v1 [physics.flu-dyn] 21 Dec 2018
1 axiv:1812.921v1 [physics.flu-dyn] 21 Dec 218 The cicula capillay jump Rajesh K. Bhagat 1, and P. F. Linden 2, 1 Depatment of Chemical Engineeing and Biotechnology, Univesity of Cambidge, Philippa Fawcett
More informationThe 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 informationChapter 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 informationNewton 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 informationFREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS
Please cite this aticle as: Izabela Zamosa Fee tansvese vibations of non-unifom beams Scientific Reseach of the Institute of Mathematics and Compute Science Volume 9 Issue pages 3-9. The website: http://www.amcm.pcz.pl/
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More informationFinal Review of AerE 243 Class
Final Review of AeE 4 Class Content of Aeodynamics I I Chapte : Review of Multivaiable Calculus Chapte : Review of Vectos Chapte : Review of Fluid Mechanics Chapte 4: Consevation Equations Chapte 5: Simplifications
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.07: Electomagnetism II Septembe 5, 202 Pof. Alan Guth PROBLEM SET 2 DUE DATE: Monday, Septembe 24, 202. Eithe hand it in at the lectue,
More informationA Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements
Poceedings of the 213 Intenational Confeence on echanics, Fluids, Heat, Elasticity Electomagnetic Fields A Thee-Dimensional agnetic Foce Solution Between Axially-Polaied Pemanent-agnet Cylindes fo Diffeent
More informationFE FORMULATIONS FOR PLASTICITY
G These slides ae designed based on the book: Finite Elements in Plasticity Theoy and Pactice, D.R.J. Owen and E. Hinton, 970, Pineidge Pess Ltd., Swansea, UK. Couse Content: A INTRODUCTION AND OVERVIEW
More informationA 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 informationCompactly Supported Radial Basis Functions
Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically
More informationTHE 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