Wojciech Bartnik, Andrzej Strużyński
|
|
- Miles Riley
- 5 years ago
- Views:
Transcription
1 Wojcech Bartnk, Andrzej Strużyńsk THE INFLENCE OF THE HYDRALIC PARAMETERS ON THE BEGINNING OF BED LOAD TRANSPORT IN MONTAIN RIVERS OBTAINED BY MEANS OF THE NISA PROGRAM Agrcultural nversty of Cracow Department of Hydraulc Engneerng The flow characterstcs n mountan rvers and streams, defnng the begnnng of bed load transport n nonunform materal, requre further knowledge of the hydraulc parameters. The most accurate results are obtaned by means of the Gordon [92] method whch descrbes the shear velocty on the bass of velocty dstrbuton. The NISA program has been used to formulate the profle velocty dstrbuton. The numercal smulaton demonstrated qualtatve agreement wth the measurement of crtcal shear stresses τ and the dmensonless parameter f under natural condtons. NOTATION C v - specfc heat of water E - constant f - shear stresses for d fracton or body force J c - crtcal slope h c - crtcal depth of flow K - thermal conductvty n j (s) - outward unt vector normal to the boundary p - pressure q - total heat flux q j - heat flux q s - heat source s - angle of straght lne accordng to = s (log y) + b or boundary parameter SF - shape factor T - temperature t - boundary parameter u - velocty u k - velocty vector - mean velocty n x-drecton * - shear velocty c - crtcal velocty for d fracton est- estmated velocty y - dstance from the wall Y - depth of water w - weghtng functon κ - von Karman constant
2 λ - second vscosty parameter µ - dynamc vscosty (coeffcent) ν - knematc vscosty of water ρ - densty of water σ j - boundary tracton τ j` - stress tensor resultng from nteracton of fluctuaton τ - shear stress τ w - wall shear stress Γ - functon of boundary condton accordng to Drchlet or Neumann ϕ - shape functon for pressure Φ - dsspaton term Ψ - shape functon for velocty Superscrpts: - - mean quantty - fluctuaton quantty + - dmensonless wall functon quantty KEY WORDS turbulent flow, velocty dstrbuton, bed load transport, ncpent moton, NISA program INTRODCTION Although new research pertanng to flow characterstcs n mountan rvers and streams, defnng the begnnng of bed load transport n non-heterogeneous materals requres further knowledge of hydraulc parameters, the drect measurement of shear stresses n natural condtons s qute complcated. Therefore, these have been determned on the bass of the ndrect method from the calculaton of local shear velocty from the velocty dstrbuton. The turbulent velocty dstrbuton over a rough bad has been defned by Prandtl and Karman. Ther unversal equaton has been modfed by Semons and Sentürk [1977] : * 30 y = 5. 75log k The measurement of the bottom level by means of Prandtl and Karman unversal equaton s taken nto consderaton. Accurate results are obtaned on the bass of the Wlcock [1994] method. One may use the k value obtaned from k = 2,85 d m.. In some lterature another method may be found to descrbe the frcton velocty - namely Gordon and McMahon [1992]: * = s (2) (1)
3 Shear velocty s calculated on the bass of the Prandtl and Karman equaton. The slope lne estmated accordng to the velocty dstrbuton n the semlogarthmc graph s used to calculate shear velocty (eq. 2). It may be utlsed when the velocty profle n the semlogarthmc scale becomes close to a straght lne [Bartnk,Stru yñsk 1996]. Ths allows to calculate shear stresses 2 τ = ρ * (3) The slope of the lne s grows together wth the slope of the water level and bottom roughness. Accordng to Wlcock, the velocty profle s n agreement wth the logarthmc dstrbuton wthn 2d 90 <y<y/5. For the calculated profles, the range s at least equal to the above. In order to calculate the profle velocty dstrbuton, the NISA program has been used. NISA PROGRAM Theoretcal bass of Nsa program. NISA s a specal program for computer aded engneerng (CAE) snce t covers a wde spectrum of flud flow analyss [EMRC 1995]. NISA/3D-FLID computes lamnar and turbulent flows of Newtonan and Non-Newtonan fluds. The program s nterfaced wth a DISPLAY, whch allows to create a graphcal tool for the presentaton of results,.e. plots of pressure, velocty, temperature and densty number. The equatons of contnuty, momentum and energy become ρ u t + u k u x ρ + ( ρu k ) = 0 t x p λ u = x + x x k k k k µ u u + x + j x j x j + ρf (4) (5) ρc v T t + u k T p u K T k x = k x + x x k j j ρq s + Φ + (6) The equatons (4) - (6) are vald for both lamnar and turbulent flow, but a very fne grd s requred to predct turbulent moton n the flow feld. To overcame the dffculty, Reynolds averagng technque [Tennekes 1972] s used to modfy the Naver - Stokes equatons for turbulent flow. The physcal equatons charactersng the flow feld are wrtten as :
4 u = u + u` p = p + p` ρ = ρ + ρ` T = T + T ` (7) where, the quanttes appearng n the governng equatons can be consdered to consst of mean and fluctuatng parts. The equatons are then averaged over a fnte nterval of tme. Mathematcally, the governng equatons are a set of ellptc, second-order partal dfferental equatons. Therefore, the type of boundary condtons by Drchlet or Neumann are n fact close to the real state. These are wrtten as u = u (s,t) on Γ u σ j =σ (s,t) n j (s) on Γ σ (8) T = T (s,t) on Γ T q = q a (s,t) +q c (s) +q r (s) on Γ q where q a, q c and q r are the appled, convectve and radatve heat flux, respectvely. In turbulent flow, hgh gradents of velocty exst near a sold wall (fg.1). To resolve ths phenomenon, the wall functons for turbulent flow s used : y Ym + m = *, m =, j (9) v and as * = τ w ρ s the shear velocty. Next, a non-dmensonal velocty + m s defned + m m =, m =, j (10) τ Fg.1.Turbulent magnary plane
5 Based on the dstance Y m j and the prevous nteracton soluton, * can be calculated as follows or = Y for Y < (11) j j j ( ) j = EYj for Yj 11 6 κ ln. (12) Numercal problems are descrbed by usng Galerkn s method of Fnte Element Formulaton. The flow dscretzng s done by a number of smple shaped regons called elements. Wthn each element the dependent varable s approxmated by smple polynomal functons. The coeffcents of these polynomals are obtaned from nodal values of the dependent varable. Mathematcally, the velocty, pressure and temperature n an algnment can be wrtten as [Chung 1978] : T (, ) = ψ ( ) T ( ) ( ) T (, ) = θ T ( t) u x t u t p x, t = ϕ p t (13) T x t where the unknown u, p and T are column vectors of nodal ponts and Ψ, ϕ, θ are column vectors of shape functons. These are substtuted nto the governng equatons n the form of ( ψ ϕ θ ) f,,, u, p, T = R (14) where R s the resdual resultng from the use of the approxmatons of eq.13. The Galerkn form of the method of weghted resduals seeks to reduce the error R (resdual) to zero. Ths s done by achevng ortogonalty between the resdual and weghtng functons of the element whch s expressed as ( f w) dω ( R w) dω Ω E = 0 (15) NISA/3D-FLID contans a lbrary of element types whch can be used to solve a wde range of ssues. These element types can be dvded nto two man categores, namely two-dmensonal (or axsymmetrc) and three dmensonal. For two-dmensonal or axsymmetrc ssues, the followng elements may be used : a ) 4 node soparametrc quadrlateral, b ) 8 node soparametrc quadrlateral, c ) 3 node soparametrc trangle, d ) 6 node soparametrc trangle. Ω E
6 NISA/3D-FLID Element Lbrary Nsa/3-D consst of the followng executve commands for 2-D elements (NKTP=2) (Tab.1): Executve Commands Table1 Element type NKTP = 2 Analyss Types 1. FLID 2. HEAT 3. FLHT Degrees of Freedom 1. 2 per node :, V 2. 1 per node : T 3. 3 per node :, V, T NORDR - Quadrlateral : 4, 8 or 9 nodes Shape / number of nodes ( NORDR = 1, 2, 30 ) - Trangle : 3 or 6 nodes ( NORDR = 10, 11 ) Materal Propertes Newtonan Flud Natural Convecton Porous Meda 4 propertes : DEN, VISC, COND, SPEC 8 propertes : DEN, VISC, COND, SPEC,, GRAV,BETA, TEMB, ANGL 5 propertes : DEN, VISC, COND, SPEC, RIST Nodal Output - VELOCITY (, V ) - PRESSRE ( P ) - TEMPERATRE ( T ) - STREAM FNCTION ( S ) - Turbulent : Knetc energy ( K ), dsspaton rate ( E ), eddy vscosty ( EDD ) - Newtonan and Non-Newtonan flud : shear stress ( SXX, SYY, SXY ) - Vortcty ( VOR ) Calculaton of the velocty profle The followng executve commands were ON for the typcal flud analyss run n NISA/3D-FLID: analyss types: FLID degrees of freedom: 2 per mode:, V shape and No. of node: Quadrlateral:4, Trangle:3 materal propertes: DEN,VISC, nodal output: Velocty (,V) Turbulence
7 Fg.2. Mesh for sphercal grans. Fg.3. Mesh for ellpsodal grans. The velocty profles were constructed for unform grans. The roughness was defned by mesh for sphercal and for ellpsodal grans (fg.2,3). In ths numercal experment, the gran dameter for both cases was 0.08 m. The comparson of resultng for plot velocty dstrbuton shows greater nfluence of bed materal on sphercal grans than on ellpsodal grans. Accordng to ths phenomenon, the shear velocty for the known parameter s (fg.4,5) and the shear stresses (tab.3) were calculated. The numercal results of velocty dstrbuton, calculatons of frcton velocty and shear stresses were compared wth the experment results. The values of crtcal shear stresses for modelled grans are very close to Sheld s dagram [Kenedy 1995] for ncpent moton of unform gran (respectvely f = and 0.046). Recently, Bartnk and Mchalk [1994] presented a modfed Sheld's dagram for the ncpent moton of nonunform sedment. The correct determnaton of crtcal shear stresses s also strctly connected wth the shape factor of grans n ths dagram. Ths was acheved by
8 [m] 1 s = est [m s -1 ] Fg.4. Plot of velocty dstrbuton for sphercal grans. [m] 1 0,1 s = 1.40 est 0,01 0,70 1,20 1,70 2,20 [m s -1 ] Fg.5. Plot of velocty dstrbuton for ellpsodal grans. usng the radoactve tracer method and allowed to perform the drect measurement of crtcal parameters. The values of crtcal shear stresses for ellpsodal grans are very smlar to measurements under natural condtons ( f = and respectvely).
9 Results of velocty dstrbuton and shear stresses analyss Table 2. shape of grans sphercal ellpsodal accordng to eq.(2) accordng to eq.(1) SF = 1 SF = 0.6 y y y unform grans nonunform grans [m.] [ms -1 ] [ms -1 ] d [m] d m [m] SF SF s k *[ms -1 ] *[ms -1 ] τ [Nm -2 ] τ [Nm -2 ] f f Crtcal parameter of bed load moton under natural condton for the Raba rver Table 3. Year, number of run, rver Raba Raba Raba Raba d [m] > 0.08 h [m] I c τ c [Nm -2 ] f _ c [ms -1 ]
10 CONCLSIONS The paper presents a new methodology for calculatng crtcal condtons whch lead to the begnnng of bed load transport n mountan rvers and streams. Ths methodology dffers from the prevously employed methods based on measurements under natural condtons. The new methodology has the advantage of beng smple and flexble. Modfcatons can be easly made because more data becomes avalable. The tests also demonstrated a qualtatve agreement wth the measured crtcal shear stresses τ and the dmensonless parameter f for the natural condtons based on the Raba rver experments. Acknowledgements The materal presented s based on the NISA/DISPLAY program suppled by the INPOL COMPANY LTD of Cracow. REFERENCE Bartnk. W Fluval Hydraulcs of Mountan Rver wth Moble Bed. Begnnng of Bed Load Moton., Zeszyty Naukowe Agrcultural nversty of Cracow, no. 171, n polsh. Bartnk W., Mchalk Al Fluval Hydraulcs of Streams and Mountan Rvers wth Moble Bed., Hydraulcs Conference 1994, Buffalo, SA, Amercan Socety of Cvl Engneers, N.York, pp Bartnk W., Srużyńsk A Measurements of Basc Velocty Fluctuaton Characterstc n Rough Stream Bed., Zeszyty Naukowe Agrcultural nversty of Cracow., no. 16, pp.19-30, n polsh. Chung. T.J Fnte Element Analyss n flud dynamcs., McGraw Hll. Gordon. N. D., McMahon A Stream Hydrology., London. Jarrett. R. D Hydraulcs of Hgh-Gradent Streams., Journal of the Hydraulcs Dvson., Amercan Socety of Cvl Engneers, vol. 110, no.11, pp Kenedy J.F The Albert Shelds Story., Journal of Hydraulc Engneerng., Vol.121, No.11, pp Moore. M., Dplas. P Effects of Partcle Shape on Bed Load Transport., Hydraulc Engneerng 94, pp Semens. D. B., Sentürk. F Sedment Transport Technology., Ford Collns. Tennekes H., Lumley J.L A frst Course n Turbulence., The MIT Press, Cambrdge, MA. Wlcock. P. R., Barta. A. F., Shea. C. C Estmatng Local Bed Shear Stress In Large Gravel-Bed Rvers., Hydraulc Engneerng 94, pp Engneerng Mechancs Research Corporaton, NISA/3D-FLID (INCOMPRES- SIBLE) SERS MANAL., Centre for Engneerng and Computer Technology, Mchgan SA.
Module 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationIntroduction to Turbulence Modeling
Introducton to Turbulence Modelng Professor Ismal B. Celk West Vrgna nversty Ismal.Celk@mal.wvu.edu CFD Lab. - West Vrgna nversty I-1 Introducton to Turbulence CFD Lab. - West Vrgna nversty I-2 Introducton
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationLecture 12. Modeling of Turbulent Combustion
Lecture 12. Modelng of Turbulent Combuston X.S. Ba Modelng of TC Content drect numercal smulaton (DNS) Statstcal approach (RANS) Modelng of turbulent non-premxed flames Modelng of turbulent premxed flames
More informationPublication 2006/01. Transport Equations in Incompressible. Lars Davidson
Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport
More informationNUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD
THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and
More informationA Numerical Study of Heat Transfer and Fluid Flow past Single Tube
A Numercal Study of Heat ransfer and Flud Flow past Sngle ube ZEINAB SAYED ABDEL-REHIM Mechancal Engneerng Natonal Research Center El-Bohos Street, Dokk, Gza EGYP abdelrehmz@yahoo.com Abstract: - A numercal
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationTurbulent Flow. Turbulent Flow
http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and
More informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More information2) For a two-dimensional steady turbulent flow in Cartesian coordinates (x,y), with mean velocity components (U,V), write
058:68 Turbulent Flows 004 G. Constantnescu HOMEWORKS: Assgnment I - 01/6/04, Due 0/04/04 1) A cubcal box of volume L 3 s flled wth flud n turbulent moton. No source of energy s present, so that the turbulence
More informationA new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems
Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationSimulation of Turbulent Flow Using FEM
Internatonal Journal of Engneerng and Technology Volume 2 No. 8, August, 2012 Smulaton of Turbulent Flow Usng FEM Sabah Tamm College of Computng, AlGhurar Unversty, Duba, Unted Arab Emrates. ABSTRACT An
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More informationA Cartesian-grid integrated-rbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesan-grd ntegrated-rbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
More informationTurbulence and its Modelling
School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationUncertainty as the Overlap of Alternate Conditional Distributions
Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationCONTROLLED FLOW SIMULATION USING SPH METHOD
HERI COADA AIR FORCE ACADEMY ROMAIA ITERATIOAL COFERECE of SCIETIFIC PAPER AFASES 01 Brasov, 4-6 May 01 GEERAL M.R. STEFAIK ARMED FORCES ACADEMY SLOVAK REPUBLIC COTROLLED FLOW SIMULATIO USIG SPH METHOD
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
More informationSolutions to Exercises in Astrophysical Gas Dynamics
1 Solutons to Exercses n Astrophyscal Gas Dynamcs 1. (a). Snce u 1, v are vectors then, under an orthogonal transformaton, u = a j u j v = a k u k Therefore, u v = a j a k u j v k = δ jk u j v k = u j
More informationFlow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved
Smulaton of nose generaton and propagaton caused by the turbulent flow around bluff bodes Zamotn Krll e-mal: krart@gmal.com, cq: 958886 Summary Accurate predctons of nose generaton and spread n turbulent
More informationEXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017
EXAMPLES of THEORETICAL PROBLEMS n the COURSE MMV03 HEAT TRANSFER, verson 207 a) What s eant by sotropc ateral? b) What s eant by hoogeneous ateral? 2 Defne the theral dffusvty and gve the unts for the
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More information3-D Numerical computations of turbulence in a partially vegetated shallow channel
Rver Flow 010 - Dttrch, Koll, Aberle & Gesenhaner (eds) - 010 Bundesanstalt für Wasserbau ISBN 978-3-93930-00-7 3-D Numercal computatons of turbulence n a partally vegetated shallow channel D. Soulots
More informationTurbulent Heat Transfer Modelling in a Vessel Stirred by a Pitched Blade Turbine Impeller
Turbulent Heat Transfer Modellng n a Vessel Strred by a Ptched Blade Turbne Impeller BARBARA ZAKRZEWSKA, ZDZISŁAW JAWORSKI Chemcal Engneerng Faculty Szczecn Unversty of Technology Al. Pastó 42, 70-065
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationTURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March
TURBULENT FLOW A BEGINNER S APPROACH Tony Saad March 2004 http://tsaad.uts.edu - tsaad@uts.edu CONTENTS Introducton Random processes The energy cascade mechansm The Kolmogorov hypotheses The closure problem
More informationA solver for free-surface flow in heterogeneous porous media
A solver for free-surface flow n heterogeneous porous meda Olver Oxtoby Johan Heyns Aeronautc Systems, Councl for Scentfc and Industral Research Pretora, South Afrca Free-surface flow: Sloshng Smple small-ampltude
More informationFour Parameter Heat Transfer Turbulence Models for Heavy Liquid Metals
J. Energy Power Sources Vol., No., 015, pp. 63-74 Receved: January 9, 015, Publshed: February 8, 015 Journal of Energy and Power Sources www.ethanpublshng.com Four Parameter Heat Transfer Turbulence Models
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationThree-dimensional eddy current analysis by the boundary element method using vector potential
Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationNUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS
ISSN: 345-3109 RCEE Research n Cvl and Envronmental Engneerng www.rcee.com Research n Cvl and Envronmental Engneerng 014 (04) 190-198 NUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS Rasoul Daneshfaraz
More informationSupplemental Material: Causal Entropic Forces
Supplemental Materal: Causal Entropc Forces A. D. Wssner-Gross 1, 2, and C. E. Freer 3 1 Insttute for Appled Computatonal Scence, Harvard Unversty, Cambrdge, Massachusetts 02138, USA 2 The Meda Laboratory,
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More informationBasic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos
Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationSIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract
SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME PACS REFERENCE: 43.20.Mv Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Zeunerstr.
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationThe Analysis of Convection Experiment
Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 5) The Analyss of Convecton Experment Zlong Zhang School of North Chna Electrc Power Unversty, Baodng 7, Chna 469567@qq.com Keywords:
More informationResearch & Reviews: Journal of Engineering and Technology
Research & Revews: Journal of Engneerng and Technology Case Study to Smulate Convectve Flows and Heat Transfer n Arcondtoned Spaces Hussen JA 1 *, Mazlan AW 1 and Hasanen MH 2 1 Department of Mechancal
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationNUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST
Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationFeb 14: Spatial analysis of data fields
Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationThe Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially
Open Journal of Flud Dynamcs, 2015, 5, 183-187 Publshed Onlne June 2015 n ScRes. http://www.scrp.org/journal/ojfd http://dx.do.org/10.4236/ojfd.2015.52020 The Tangental Force Dstrbuton on Inner Cylnder
More informationHandout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models
Handout: Large Eddy mulaton I 058:68 Turbulent flows G. Constantnescu Introducton to ubgrd-cale (G) Models G tresses should depend on: Local large-scale feld or Past hstory of local flud (va PDE) Not all
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationSurvey of applications of discrete vortex method in civil engineering
Budownctwo Archtektura 5 (2009) 29-38 Survey of applcatons of dscrete vortex method n cvl engneerng Tomasz Nowck Lubln Unversty of Technology, Faculty of Cvl Engneerng and Archtecture, Department of Structural
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
More informationCFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL
CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL 1. Introducton Chrstophe Vallée and Thomas Höhne In dfferent scenaros of small break Loss of Coolant Accdent (SB-LOCA), stratfed twophase
More informationFORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO
ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO Mohammad Rahnama*, Seyed-Mad Hasheman*, Mousa Farhad**
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationNumerical analysis on the agglomeration behavior of fine particles in plane jet
Numercal analyss on the agglomeraton behavor of fne partcles n plane et Mn Guo 1, a, Ja L 2, Xn Su 1 and Gurong Yang 1 1 Tann Academy of Envronmental Scences, Tann, 300191, Chna; 2 The Fourth Research
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationIntroduction to Turbulence. by Håkan Gustavsson Division of Fluid Mechanics Luleå University of Technology
Introducton to Turbulence by Håkan Gustavsson Dvson of Flud Mechancs Luleå Unversty of Technology Foreword The ntenton wth these pages s to present the student wth the basc theoretcal concepts of turbulence
More informationGeoSteamNet: 2. STEAM FLOW SIMULATION IN A PIPELINE
PROCEEDINGS, Thrty-Ffth Workshop on Geothermal Reservor Engneerng Stanford Unversty, Stanford, Calforna, February 1-3, 010 SGP-TR-188 GeoSteamNet:. STEAM FLOW SIMULATION IN A PIPELINE Mahendra P. Verma
More informationPORE STRUCTURE AND THERMAL CONDUCTIVITY OF BURNT CLAY BRICKS INTRODUCTION
PORE STRUCTURE AND THERMAL CONDUCTIVITY OF BURNT CLAY BRICKS Olga Koronthalyova, Peter Matasovsky Insttute of Constructon and Archtecture, Slovak Academy of Scences, Dubravska 9, 845 43 Bratslava, Slovaka.
More informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationTurbulent Natural Convection of Heat with Localized Heating and Cooling on Adjacent Vertical Walls in an Enclosure
Internatonal Journal of Scence and esearch (IJS ISSN (Onlne: 319-7064 Turbulent Natural Convecton of Heat wth Localzed Heatng and Coolng on Adacent Vertcal Walls n an Enclosure 1 Chrstopher Onchaga, J.K.Sgey,
More informationAmplification and Relaxation of Electron Spin Polarization in Semiconductor Devices
Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More information