Particle dynamics class, SMS 618, (Emmanuel Boss 11/19/2003) Van Rijn s TRANSPOR lab: computation of sediment transport in current and wave direction.

Similar documents
Skin Effect and Formation Damage

Photographing a time interval

VISUAL PHYSICS ONLINE EQUATION MINDMAPS

1 Fundamental Solutions to the Wave Equation

TRAVELING WAVES. Chapter Simple Wave Motion. Waves in which the disturbance is parallel to the direction of propagation are called the

Extra Examples for Chapter 1

User s Guide NBC 2005, Structural Commentaries (Part 4 of Division B)

Circular Motion Problem Solving

AVS fiziks. Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES

one primary direction in which heat transfers (generally the smallest dimension) simple model good representation for solving engineering problems

ANALYSIS AND TEMPERATURES CONTROL IN A TUBULAR CHEMICAL REACTOR

1 Fundamental Solutions to the Wave Equation

Mass Transfer (Stoffaustausch)

Physics 218, Spring March 2004

h=1 cm 1.03 g/cm 1.43 g/cm

On a proper definition of spin current

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

Mass Transfer. Dispersion. Lecture 13, , Dr. K. Wegner

Capacitance Extraction. Classification (orthogonal to 3D/2D)

Example 1. Centripetal Acceleration. Example 1 - Step 2 (Sum of Vector Components) Example 1 Step 1 (Free Body Diagram) Example

Non-Ideal Gas Behavior P.V.T Relationships for Liquid and Solid:

Chemical Reaction Engineering

Fall 2004/05 Solutions to Assignment 5: The Stationary Phase Method Provided by Mustafa Sabri Kilic. I(x) = e ixt e it5 /5 dt (1) Z J(λ) =

3. Perturbation of Kerr BH

(conservation of momentum)

Why do Golf Balls have Dimples on Their Surfaces?

8.022 (E&M) Lecture 13. What we learned about magnetism so far

Adaptive LQ Cascade Control of a Tubular Chemical Reactor

Time Dilation in Gravity Wells

% ionic character = {1 exp 0.25 X 2

Finite Difference Solution of Mixed Convective Heat Transfer Transient Flow along a Continuously Moving Cooled Plate

Neoclassical transport processes in weakly collisional plasmas with fractured velocity distribution functions

Lecture 2 Date:

Thermal-Fluids I. Chapter 17 Steady heat conduction. Dr. Primal Fernando Ph: (850)

Generalized Vapor Pressure Prediction Consistent with Cubic Equations of State

Lecture No. 6 (Waves) The Doppler Effect

Lecture 2 - Thermodynamics Overview

EGN 3353C Fluid Mechanics

8 Separation of Variables in Other Coordinate Systems

Forging Analysis - 2. ver. 1. Prof. Ramesh Singh, Notes by Dr. Singh/ Dr. Colton


Determining the Best Linear Unbiased Predictor of PSU Means with the Data. included with the Random Variables. Ed Stanek

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

Chapter 2 Lecture 5 Longitudinal stick fixed static stability and control 2 Topics

Suppose the medium is not homogeneous (gravity waves impinging on a beach,

A Relativistic BCS Theory of S-Wave Superconductivity

2 Parallel-Plate Transmission Line (Geometric Model) = c Assume it s a plane wave propagate in the z with polarization in y direction. d dz ~ ˆ.

But for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,

In electrostatics, the electric field E and its sources (charges) are related by Gauss s law: Surface

A) At each point along the pipe, the volume of fluid passing by is given by dv dt = Av, thus, the two velocities are: v n. + ρgy 1

Water Tunnel Experiment MAE 171A/175A. Objective:

PROPAGATION OF PHOTON IN RESTING AND MOVING MEDIUM. J. Zaleśny. Institute of Physics, Technical University of Szczecin, A b s t r a c t

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

E(r,t) = e 3. r 3. (b) Show that the transverse current, J t,is 3n(n e 3 ) e 3

FI 2201 Electromagnetism

T h e C S E T I P r o j e c t

TOWARD IMPROVING PREDICTION OF SEDIMENT TRANSPORT OVER WAVE-INDUCED RIPPLES

not to be republished NCERT ELECTROMAGNETIC WAVES Chapter Eight MCQ I

FE FORMULATIONS FOR PLASTICITY

4.4 Buoyant plume from a steady heat source

Inference for A One Way Factorial Experiment. By Ed Stanek and Elaine Puleo

Chemical Reaction Engineering

CBE Transport Phenomena I Final Exam. December 19, 2013

COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE

Question 1: The dipole

PROCESSING QUALITY OF NEW POTATO CULTIVARS FOLLOWING PROLONGED STORAGE

PHZ 3113 Fall 2017 Homework #5, Due Friday, October 13

Solution to Problem First, the firm minimizes the cost of the inputs: min wl + rk + sf

Physics Courseware Electromagnetism

Three dimensional flow analysis in Axial Flow Compressors

Ultrasonic Measurement Models for Imaging with Phased Arrays

1 Fundamental Solutions to the Wave Equation

Part 2: CM3110 Transport Processes and Unit Operations I. Professor Faith Morrison. CM2110/CM Review. Concerned now with rates of heat transfer

16.1 Permanent magnets

POISSON S EQUATION 2 V 0

LEARNING FROM MISTAKES

Correspondence Analysis & Related Methods

V7: Diffusional association of proteins and Brownian dynamics simulations

6.4 Period and Frequency for Uniform Circular Motion

Revised Newtonian Formula of Gravity and Equation of Cosmology in Flat Space-Time Transformed from Schwarzschild Solution

Simulation thus allows us to check the efficacy of our screening tools.

10. Groundwater in geotechnical problems

Stability of Rock Berm under Wave and Current Loading

Ion-sound waves (electrostatic low frequency waves)

NONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor

Chapter 1 Functions and Graphs

Eddy Currents and Magnetic Calibrations in LDX using a Copper Plasma. D.P. Boyle, PPPL M.E. Mauel, D.T. Garnier, Columbia J.

Gravity. David Barwacz 7778 Thornapple Bayou SE, Grand Rapids, MI David Barwacz 12/03/2003

Special Relativity in Acoustic and Electromagnetic Waves Without Phase Invariance and Lorentz Transformations 1. Introduction n k.

Hopefully Helpful Hints for Gauss s Law

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

2.25 Advanced Fluid Mechanics

A method for solving dynamic problems for cylindrical domains

Current, Resistance and

PROBLEM SET 5. SOLUTIONS March 16, 2004

Correlation between System (Lagrangian) concept Control-volume (Eulerian) concept for comprehensive understanding of fluid motion?

GCSE: Volumes and Surface Area

cos kd kd 2 cosθ = π 2 ± nπ d λ cosθ = 1 2 ± n N db

Retrieval of three-dimensional distribution of rainfall parameters for rain attenuation correction using multi-parameter radar

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

Transcription:

Patile dynami la SMS 618 (Emmanuel Bo 11/19/003 Van Rijn TRANSPOR lab: omputation of ediment tanpot in uent and ave dietion Handout: Appendix A in van Rijn 1993 Piniple of ediment tanpot in ive etuaie and oatal ea Aqua Publiation Note: in ti ite-up log i log 10 ile ln i log e Input (ee Fig A1 in appendix fo dietion: -ate dept [m] v - dept-aveaged veloity in main uent dietion [m/] u - time-aveaged and dept-aveaged etun veloity belo ave toug ompenating te ma tanpot beteen ave et and toug [m/] u b - time-aveaged nea-bed veloity due to ave ind o lateal denity gadient [m/] H- ignifiant ave eigt (ignifiant ave eigt H i appoximately equal to te aveage of te iget one-tid of te ave H i alulated uing: H = 40 * qt(m 0 ee m 0 i te vaiane of te ave diplaement time eie aquied duing te ave aquiition peiod Vaiane an alo be alulated uing te nondietional ave petum aoding to te folloing elationip: m 0 = um(s(f*df ee te ummation of petal denity S(f i ove all fequeny band f of te nondietional ave petum and df i te bandidt of ea band Wave analyi ytem typially um ove te ange fom 003 to 040 Hz it fequeny bandidt of 001 Hz (ee ttp://ndbnoaagov/avealtml [m] Tp- Dominant o peak ave peiod (te peiod oeponding to te fequeny band it te maximum value of petal denity in te nondietional ave petum It i te eipoal of te peak fequeny fp: Tp = 1/fp Dominant peiod i epeentative of te ige ave enounteed duing te ave ampling peiod [] φ-angle beteen te dietion of ave popagation and main uent dietion [ ] d 50 - median diamete of bed mateial [m] d 90-90 t peentile diamete of bed mateial [m] d -epeentative diamete of upended mateial (ill be in te ange of (06 1 d 50 [m] k - uent-elated bed ougne eigt (minimum 001m (likely ange 001 1m [m] k - ave-elated bed ougne eigt (minimum 001m (likely ange 001 1m fo eet flo 001m [m] Te- tempeatue of fluid [ C] SA- alinity of fluid [ppt] Aumption: ρ =650kg/m 3 ediment denity κ=04 Von Kaman ontant Geneal paamete omputed: Cloinity: CL=(SA-003/1805 Fluid denity: ρ=1000+1455cl-00065(te-4+04cl [Kg/m 3 ] Kinemati vioity: ν=(4/(0+te10-5 [m /]

Sinking veloity: ( 1 gd 1< d < 100µ m 18ν 3 10ν 001( 1 gd 100 < d < 1000µ m 1+ 1 ; = ρ d ν 1000µ m< d 11 ( 1 gd ρ Sediment aateiti omputed: = ρ /ρ elative denity D * =d 50 [(-1g/ν ] 1/3 - Patile paamete 1 04D* 1< D* < 4 064 014D* 4< D* < 10 01 Sield paamete: θ 004D* 10 < D* < 0 09 0013D 0 < D * * < 150 0055 150 < D* =( ρ -ρgd 50 θ itial ea te u = 575 ( 1 gd θ log 4 d Citial dept-aveaged veloity U ( 50 90 ( 50 p ( 109( 1 gd 50 Tp /3 01( 1 g d T d < 00005m -- Citial peak obital veloity 50 0571 075 00005m< d50 Wave paamete omputed L Dopple ifted avelengt (ange in avelengt due to peene of mean flo: L' gl' π find L te olution of: vr oφ = tan T p π L' T p - Dopple ifted ave peiod: T p Tp ' = 1 o / L' ( vt R p φ Nea-bed peak obital exuion: A = in( π L' Nea-bed peak obital veloity: π H U = T' in π L' p H ( Wave-bounday laye tikne: = 007A ( A k 05

Nea-bed peak obital veloity in foad dietion: 3π H U + 4 001gTp < 4 Tp ' L' ( in( π L' U f H 1 03 U 001gTp > + Nea-bed peak obital veloity in bakad dietion: 3π H U 4 001gTp < 4 Tp ' L' ( in( π L' U b H 1 03 U 001gTp > 015 gh Retun veloity ma tanpot: u = ( 095 035H U f U b Nea bed ave indued mean veloity: u 005 05 b = U U f U + b Te lat to ae input paamete To ave tem omputed by te pogam inet 9 at te input line Bed paamete omputed: Appaent bed ougne: U ka 10k k exp γ < π ka ( v =08+ -03 R + u γ β β β = φ R 360 oteie 10k Effetive time-aveaged bed-ea tee omputed: log ( 13d90 Effiieny fato fo uent: µ log ( 1k exp 6 + 5 A 3d Effiieny fato fo ave: µ = exp 6 + 5 A k Wave-uent inteation oeffiient: α ( 90 ( 019 019 ( k o ( k ln 90 1+ ln 30 a = 1 ln ( 90 k 1+ ln( 30ka

1 04 Bed ea te due to uent: = ρ ( vr + ur 8 ( log ( 1k 1 exp 6 5 = ρ A k + U 8 Te total ea i te um of te ea: = + Bed ea te due to ave: ( 019 Effetive bed ea veloity uent: = ( u α µ ρ '* Bed ea te paamete omputation: ( α µ + µ Dimenionle bed-ea te fo bed load tanpot: T = Dimenionle bed-ea te fo efeene onentation at z=a: ( α µ + ( 06 / D* Ta = Bot ae zeo if negative Cuent veloity pofile omputation: ( urz vrz ( ur vr ln[ 30z ka] 1+ ln[ 30ka] ( ur vr ln90 [ ka] ln30 [ z k ] 1 ln30 [ k] ln90 [ k ] z 3 z < 3 + a Sediment eddy diffuion pofile omputation: Eddy diffuion due to mean uent: κβu* z( 1 z z < 05 ε ; 05κβu* z 05 ( ( u g u * β + v 1+ < 15 ee: β ; u* = 15 oteie 18log ( 4d90 Eddy diffuion due to ave: εbed = 0004D* U z < ε εmax = 0035 H Tp z 05 ; z εbed + ε max εbed 05 < z< 05 =03 H but =0m & =005m max min Combined uent and ave eddy diffuion: ε = [ ε + ε ] 05

Sediment volume onentation omputation: Refeene level: a=maximum(k k 15 d50 Ta Bed onentation (fo z a: Ca = 0015 03 a D* Fo z a te onentation gadient i give by: 5 ( 1 ( ( dc = dz ε 08 04 ( 1 + / 065 / 065 Obtain C by integation to te ufae fom a (ee onentation i knon to Time aveaged upended load tanpot omputation: ρ vr zcdz uent dietion a q ρ ave dietion ur zcdz a Time aveaged and intantaneou bed load tanpot ate omputation: vrln[ 30 ka] vr 1+ ln[ 30ka] Veloitie ae omputed at =max(3 k : = u Rln[ 30 v ka] ur v 1 + ln[ 30 ka] U = U o ϕ i te pae of te ave Intantaneou ave veloity: ( ϕ Intantaneou veloity along uent: U x U oφ vr ( ub ur Intantaneou veloity ao uent: U y = U inφ + ( ub + ur inφ Intantaneou veloity: U = U + U R x y Intantaneou fition oeffiient: α = v R = + + + oφ v R + U Intantaneou bed ea te: 1 ln( 30k + 04 019 ' 05 05 ( 1 exp 6 5 ( ( b = ρ α α A 3d + + U ln 30 k ( log ( 4d 90 90 R Intantaneou bed-load tanpot:

03 ' b ' b b ' b > b 05γρd50 D* H S qb ρ γ = max 031 b 0 oteie Time aveaged value ae obtained by aveaging ove a ave peiod (eg ove it pae ϕ Bed fom alulation: Ti i beyond te ope of ou la Cla aignment: Ho ave affet te total tanpot Ue te TRANSPOR model it te folloing input (maybe appopiate to te Dok in font of te Daling Maine ente: =1m d 50 =150µm d 90 =400µm d =08 d 50 Te=1 C SA=30pu u = 0 v = 05 m/ k = k = 00 Ho do you expet te tanpot to ange it ave paamete (H Tp and φ ompaed to aving no ave? Wite don you pedition Invetigate o te total tanpot vaie it ange in te ignifiant ave eigt (0 m ave peiod (1 8 and angle (0 90 beteen ave and uent Ho do te eult ompae it pedition? Ho do tey (qualitatively ompae it van-rijn alulation attaed to te andout? Fo one ae vay d fom 06d 50 to d 50 Ho ould you expet te total tanpot to ange? Ho mu ange do you obeve in te total tanpot?