Movement of Pollutants in Lakes

Size: px
Start display at page:

Download "Movement of Pollutants in Lakes"

Transcription

1 Movement of Pollutants in Lakes Simplest pattern of water movement in a lake: win-riven currents: Moule 2: Surface Waters, Lecture Chemical Fate an Transport in the Environment, 2 n eition. H.F. Hemon an E.J. Fechner-Levy. Acaemic Press. Lonon Principles of Surface Water Quality Moeling an Control. R.V. Thomann an J.A. Mueller. Harper & Row, New York Hemon an Fechner-Levy 2000 Fickian Mixing in Lakes A mass of tracer injecte into a lake will move by avection with the water currents, but will also sprea out into an ever-larger volume of water. Given enough time, it will ten to become completely mixe. This mixing is primarily ue to turbulence, carrying chemicals away from regions of higher concentrations to areas of lower concentrations. Hemon an Fechner-Levy

2 Concentrations for an instantaneous ischarge into a twoimensional boy of water (vertically mixe): Hemon an Fechner-Levy 2000 C( x, y, t) = 4πt M D D x y e 2 2 (( x Vxt) /(4Dx t) + ( y V yt ) / 4Dyt) kt e M is the mass of the chemical ischarge, per epth of water [M/L] x an y are the istances from the injection location [L] t is the time lapse since injection [T] V x an V y are the average velocity in the x an y irections [L/T] D x an D y are the Fickian transport coefficients in the x an y irections [L 2 /T] K is the first-orer ecay rate constant [T -1 ] The epth is the total epth for a vertically well-mixe lake, or the thickness of a layer in a stratifie lake. What woul lake temperatures look like in winter (very col area)? Hemon an Fechner-Levy

3 W = QeSe + QrSr + QT ST + PAsSp + SDV W= mass input [M/T] QeSe = waste effluent ischarge to lake QrSr = mass from main river QtSt = mass from tributary PAsSp = mass input from precipitation SV= seiment release Qe= effluent ischarge Qr = river flow Qt = tributary flow P = precipitation amount As = lake surface area V = lake volume Se = effluent concentration Sr = river concentration St = tributary concentration Sp = rain concentration S = seiment concentration

4 Concentration increases with time, since the start of the ischarge: Concentration Decrease after Discharge Stops Vs = W ( t) Qs kvs Change of mass with time = input mass (gain) mass outflow (loss) ecay (loss) K = Q+kV Assuming a constant Q an k over time Expaning the erivative: Vs = V s + s V V If V is temporarily constant: = 0 Then: Vs = V s Re-arranging results in: s W ( t) = V + Qs + kvs Concentration ecreases from initial conitions (at en of input) by a combination of flushing an ecay: Can simplify if efine: k ' = Q + kv Then: Resulting in: s W ( t) = V + k ' s 1 s = s 0 exp + k t t Where the lake etention time, t is efine as: V t = Q 4

5 Response ue to step loa = sum of initial conition (from previous loa), plus step loa effect: W Q Q s = + k t + s + k t Q + kv 1 exp V exp o V Total response an transitions can be etermine by calculating iniviual responses an summing the effect. Response ue to varying loa: (a) Step input an subsequent reuction to zero (b) Step input with reuction to new loa level Example Problem: A Lake (with initial s = 0), receives a loa of a slowing reacti ng pesticie (triallate) of 1080 /b/ay for 1.5 years an is then terminate. 1. Determine the equilibrium concentration 2. The maximum concentration. The time until a level of 100 µg/l is reache 5

6 (1) Determine the equilibrium concentration W W / Q s = = Q + kv 1+ kt 9 t V.15x10 ft ay year = = x x = 1. year Q 100 ft / sec 86,400sec 65ays 0 W / Q s = 1 + kt 1080lb / ay 100 ft / sec x86,400sec/ ay = = lb / ft ( 1.0year ) ay lb ft gal 454,000mg s = x x x = 1.6mg/ L = 160µ g / L ft 7.48gal.78L lb (2) Determine the maximum concentration (the max. concentration will occur at the en of the ischarge time, at t=1.5 years) s = s 1 exp 1 ( + kt ) t t yr s = 160 µ g / L 1 exp yr = 170 µ g / L yr 1.0 yr Never reaches the equilibrium concentration before it starts to ecrease. () Determine the time until a level of 100 µg/l is reache s = s 0 exp ( 1+ Kt ) t' t Where: t' = t 1. 5years 0.2 t' 100 µ g / L = 170µ g / Lexp yr yr 1.0yr t'= 2. 1 years 6

7 7

Transport of Chemicals, Kinetics, and Equilibrium

Transport of Chemicals, Kinetics, and Equilibrium Transport of Chemicals, Kinetics, and Equilibrium Module 1: Basic Concepts, Lecture Chemical Fate and Transport in the Environment, nd edition. H.F. Hemond and E.J. Fechner-Levy. Academic Press. London.

More information

ARCH 614 Note Set 5 S2012abn. Moments & Supports

ARCH 614 Note Set 5 S2012abn. Moments & Supports RCH 614 Note Set 5 S2012abn Moments & Supports Notation: = perpenicular istance to a force from a point = name for force vectors or magnitue of a force, as is P, Q, R x = force component in the x irection

More information

m (ft-lb/ft). Using the point-slope

m (ft-lb/ft). Using the point-slope ENGR 1990 Engineering athematics pplications of Derivatives E 560, E 570 Eample #1 Consier a long slener beam of length with a concentrate loa acting at istance a from the left en. Due to this loa, the

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b) LPC Physics A 00 Las Positas College, Physics Department Staff Purpose: To etermine that, for a boy in equilibrium, the following are true: The sum of the torques about any point is zero The sum of forces

More information

RFSS: Lecture 4 Alpha Decay

RFSS: Lecture 4 Alpha Decay RFSS: Lecture 4 Alpha Decay Reaings Nuclear an Raiochemistry: Chapter 3 Moern Nuclear Chemistry: Chapter 7 Energetics of Alpha Decay Geiger Nuttall base theory Theory of Alpha Decay Hinrance Factors Different

More information

2.5 SOME APPLICATIONS OF THE CHAIN RULE

2.5 SOME APPLICATIONS OF THE CHAIN RULE 2.5 SOME APPLICATIONS OF THE CHAIN RULE The Chain Rule will help us etermine the erivatives of logarithms an exponential functions a x. We will also use it to answer some applie questions an to fin slopes

More information

2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math

2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that

More information

( ) ( ) ( ) PAL Session Stewart 3.1 & 3.2 Spring 2010

( ) ( ) ( ) PAL Session Stewart 3.1 & 3.2 Spring 2010 PAL Session Stewart 3. & 3. Spring 00 3. Key Terms/Concepts: Derivative of a Constant Function Power Rule Constant Multiple Rule n Sum/Difference Rule ( ) Eercise #0 p. 8 Differentiate the function. f()

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

MATH 205 Practice Final Exam Name:

MATH 205 Practice Final Exam Name: MATH 205 Practice Final Eam Name:. (2 points) Consier the function g() = e. (a) (5 points) Ientify the zeroes, vertical asymptotes, an long-term behavior on both sies of this function. Clearly label which

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS

A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS Peifeng Ma 1 an Ole Secher Masen Accurate preiction of current velocity an bottom shear stress, which both can be significantly influence by

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Mathematics 1210 PRACTICE EXAM II Fall 2018 ANSWER KEY

Mathematics 1210 PRACTICE EXAM II Fall 2018 ANSWER KEY Mathematics 1210 PRACTICE EXAM II Fall 2018 ANSWER KEY 1. Calculate the following: a. 2 x, x(t) = A sin(ωt φ) t2 Solution: Using the chain rule, we have x (t) = A cos(ωt φ)ω = ωa cos(ωt φ) x (t) = ω 2

More information

To understand how scrubbers work, we must first define some terms.

To understand how scrubbers work, we must first define some terms. SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

Q(s, t) = S M = S M [ G 1 (t) G 2 (t) G 3 1(t) G 4 (t) ] T

Q(s, t) = S M = S M [ G 1 (t) G 2 (t) G 3 1(t) G 4 (t) ] T Curves an Surfaces: Parametric Bicubic Surfaces - Intro Surfaces are generalizations of curves Use s in place of t in parametric equation: Q(s) = S M G where S equivalent to T in Q(t) = T M G If G is parameterize

More information

ME338A CONTINUUM MECHANICS

ME338A CONTINUUM MECHANICS global vs local balance equations ME338A CONTINUUM MECHANICS lecture notes 11 tuesay, may 06, 2008 The balance equations of continuum mechanics serve as a basic set of equations require to solve an initial

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

Physics 121 for Majors

Physics 121 for Majors Physics 121 for Majors Scheule Do Post-Class Quiz #3 Do Pre-Class Quiz #4 HW #2 is ue Wenesay Quiz #1 is ue Saturay, Sept. 16 Lab #1 is set up now, ue Monay Some Department Resources Computers N-212 ESC

More information

Poroelasticity and tidal loading

Poroelasticity and tidal loading Chapter 6: Poroelasticity an tial loaing 63 Chapter 6: Poroelasticity an tial loaing 6. Introuction It is establishe above, in Chapter 4, that there are many known examples of the tial moulation of seafloor

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

five moments Moments Moments Moments 8 a force acting at a different point causes a different moment:

five moments Moments Moments Moments 8 a force acting at a different point causes a different moment: ELEENTS O RCHITECTURL STRUCTURES: OR, EHVIOR, ND DESIGN DR. NNE NICHOLS SPRING 2014 forces have the tenency to make a boy rotate about an ais lecture five moments http://www.physics.um.eu same translation

More information

3-dimensional Evolution of an Emerging Flux Tube in the Sun. T. Magara

3-dimensional Evolution of an Emerging Flux Tube in the Sun. T. Magara 3-imensional Evolution of an Emerging Flux Tube in the Sun T. Magara (Montana State University) February 6, 2002 Introuction of the stuy Dynamical evolution of emerging fiel lines Physical process working

More information

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

Calculus Class Notes for the Combined Calculus and Physics Course Semester I Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Basic Differentiation Rules and Rates of Change. The Constant Rule

Basic Differentiation Rules and Rates of Change. The Constant Rule 460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of

More information

This section outlines the methodology used to calculate the wave load and wave wind load values.

This section outlines the methodology used to calculate the wave load and wave wind load values. COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview

More information

An unstructured-grid, finite-volume model of the Irish Sea and Celtic Sea: Application to tides, surges, and temperature evolution.

An unstructured-grid, finite-volume model of the Irish Sea and Celtic Sea: Application to tides, surges, and temperature evolution. An unstructure-gri, finite-volume moel of the Irish Sea an Celtic Sea: Application to ties, surges, an temperature evolution. Jiuing Xing, Alan avies, Phil Hall an Eric Jones National Oceanography Centre,

More information

CE2253- APPLIED HYDRAULIC ENGINEERING (FOR IV SEMESTER)

CE2253- APPLIED HYDRAULIC ENGINEERING (FOR IV SEMESTER) CE5-APPLIED HYDRAULIC ENGINEERING/UNIT-II/UNIFORM FLOW CE5- APPLIED HYDRAULIC ENGINEERING (FOR IV SEMESTER) UNIT II- UNIFORM FLOW CE5-APPLIED HYDRAULIC ENGINEERING/UNIT-II/UNIFORM FLOW CE5- APPLIED HYDRAULIC

More information

Math 308 Exam I Practice Problems

Math 308 Exam I Practice Problems Math 308 Exam I Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..

More information

PES 1120 Spring 2014, Spendier Lecture 36/Page 1

PES 1120 Spring 2014, Spendier Lecture 36/Page 1 PES 0 Spring 04, Spenier ecture 36/Page Toay: chapter 3 - R circuits: Dampe Oscillation - Driven series R circuit - HW 9 ue Wenesay - FQs Wenesay ast time you stuie the circuit (no resistance) The total

More information

Chapter 3. Modeling with First-Order Differential Equations

Chapter 3. Modeling with First-Order Differential Equations Chapter 3 Moeling with First-Orer Differential Equations i GROWTH AND DECAY: The initial-value problem x = kx, x(t 0 ) = x 0, (1) where k is a constant of proportionality, serves as a moel for iverse phenomena

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Chapter 2 Derivatives

Chapter 2 Derivatives Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

The thermal wind 1. v g

The thermal wind 1. v g The thermal win The thermal win Introuction The geostrohic win is etermine by the graient of the isobars (on a horizontal surface) or isohyses (on a ressure surface). On a ressure surface the graient of

More information

Math 190 Chapter 3 Lecture Notes. Professor Miguel Ornelas

Math 190 Chapter 3 Lecture Notes. Professor Miguel Ornelas Math 190 Chapter 3 Lecture Notes Professor Miguel Ornelas 1 M. Ornelas Math 190 Lecture Notes Section 3.1 Section 3.1 Derivatives of Polynomials an Exponential Functions Derivative of a Constant Function

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics. MATH A Test #2. June 25, 2014 SOLUTIONS

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics. MATH A Test #2. June 25, 2014 SOLUTIONS YORK UNIVERSITY Faculty of Science Department of Mathematics an Statistics MATH 505 6.00 A Test # June 5, 04 SOLUTIONS Family Name (print): Given Name: Stuent No: Signature: INSTRUCTIONS:. Please write

More information

Schöck Isokorb type KST

Schöck Isokorb type KST Schöck Isokorb type Schöck Isokorb type Contents Page Element arrangements/connection layouts 288-289 Views/Dimensions 290-293 Design an capacity table 294 Torsion spring strength/notes on calculations

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Lecture 10: Logistic growth models #2

Lecture 10: Logistic growth models #2 Lecture 1: Logistic growth moels #2 Fugo Takasu Dept. Information an Computer Sciences Nara Women s University takasu@ics.nara-wu.ac.jp 6 July 29 1 Analysis of the stochastic process of logistic growth

More information

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 2 Solutions

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 2 Solutions MA 214 Calculus IV (Spring 2016) Section 2 Homework Assignment 2 Solutions 1 Boyce and DiPrima, p 60, Problem 2 Solution: Let M(t) be the mass (in grams) of salt in the tank after t minutes The initial-value

More information

Define each term or concept.

Define each term or concept. Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan

More information

FLUID MECHANICS UNIVERSITY OF LEEDS. May/June Examination for the degree of. BEng/ MEng Civil Engineering. Time allowed: 2 hours

FLUID MECHANICS UNIVERSITY OF LEEDS. May/June Examination for the degree of. BEng/ MEng Civil Engineering. Time allowed: 2 hours This question paper consists of printe pages, each of which is ientifie by the Coe Number CIVE 4 UNIVERSITY OF LEEDS May/June Examination for the egree of BEng/ MEng Civil Engineering FLUID MECANICS Time

More information

Numerical Hydraulics

Numerical Hydraulics ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction

More information

Math 308 Exam I Practice Problems

Math 308 Exam I Practice Problems Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

water adding dye partial mixing homogenization time

water adding dye partial mixing homogenization time iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,

More information

AOS 104 Fundamentals of Air and Water Pollution

AOS 104 Fundamentals of Air and Water Pollution AOS 104 Fundamentals of Air and Water Pollution Dr. Jeffrey Lew lew@atmos.ucla.edu AIM: jklew888 TuTh 1:00-3:05 MS 7124A MS1961 310-825-3023 1 Grades Homework 160 pts 2 Midterms 290 pts Final Exam Total

More information

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring

More information

R is the radius of the sphere and v is the sphere s secular velocity. The

R is the radius of the sphere and v is the sphere s secular velocity. The Chapter. Thermal energy: a minnow, an E. Coli an ubiquinone a) Consier a minnow using its fins to swim aroun in water. The minnow must o work against the viscosity of the water in orer to make progress.

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes Proceeings of the 3th WSEAS nternational Conference on SYSTEMS New Simple Controller Tuning Rules for ntegrating an Stable or Unstable First Orer plus Dea-Time Processes.G.ARVANTS Department of Natural

More information

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late

More information

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function

More information

Theoretical Estimation of Characteristics of Thermoelectric Materials Made of Nanopowders

Theoretical Estimation of Characteristics of Thermoelectric Materials Made of Nanopowders Theoretical Estimation of Characteristics of Thermoelectric Materials Mae of Nanopowers A.I. Holopkin V.N. Abrutin V.A. Romanko P.V. Miroevskiy I.A. Zotov AV-Engineering Co. Lt. Office 6 Builing 5 B. Tolmachevskiy

More information

Exam 2 Review Solutions

Exam 2 Review Solutions Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon

More information

Torque OBJECTIVE INTRODUCTION APPARATUS THEORY

Torque OBJECTIVE INTRODUCTION APPARATUS THEORY Torque OBJECTIVE To verify the rotational an translational conitions for equilibrium. To etermine the center of ravity of a rii boy (meter stick). To apply the torque concept to the etermination of an

More information

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the

More information

Estimation amount of snow deposits on the road

Estimation amount of snow deposits on the road Estimation amount of snow eposits on the roa Olga. Glaysheva oronezh State University of Architecture an Civil Engineering, Russia Email: glaov@ox.vsi.ru ABSTRACT The article uner consieration gives the

More information

Multi-robot Formation Control Using Reinforcement Learning Method

Multi-robot Formation Control Using Reinforcement Learning Method Multi-robot Formation Control Using Reinforcement Learning Metho Guoyu Zuo, Jiatong Han, an Guansheng Han School of Electronic Information & Control Engineering, Beijing University of Technology, Beijing

More information

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - III RAPIDLY VARIED FLOW

CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - III RAPIDLY VARIED FLOW CE 6 APPLIED HYDRAULIC ENGINEERING UNIT - III RAPIDLY VARIED FLOW Application of the energy equation for RVF - Critical epth an velocity - Critical, Sub-critical an Super-critical flow - Application of

More information

Mass transport in an artificial heterogeneous aquifer: Experiments and numerical modelling

Mass transport in an artificial heterogeneous aquifer: Experiments and numerical modelling Contaminant Transport in Grounwater, Kobus & Kinzelbach (es) 1989 Balkema, Rotteram SBN 9 6191879 Mass transport in an artificial heterogeneous aquifer: xperiments an numerical moelling GSchafer & HKobus

More information

Study on aero-acoustic structural interactions in fan-ducted system

Study on aero-acoustic structural interactions in fan-ducted system Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,

More information

Physics 2212 K Quiz #2 Solutions Summer 2016

Physics 2212 K Quiz #2 Solutions Summer 2016 Physics 1 K Quiz # Solutions Summer 016 I. (18 points) A positron has the same mass as an electron, but has opposite charge. Consier a positron an an electron at rest, separate by a istance = 1.0 nm. What

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

Interpreting Petterssen s Extratropical Cyclone* Self-Development Eqn. (based in part on notes by James T. Moore)

Interpreting Petterssen s Extratropical Cyclone* Self-Development Eqn. (based in part on notes by James T. Moore) Interpreting Petterssen s Extratropical Cyclone* Self-Development Eqn. (base in part on notes by James T. Moore) P.S. Market ATMS-97 Avance Synoptic Meteorology University of Missouri *Works for anticyclones

More information

Math 308, Sections 301, 302, Summer 2008 Lecture 5. 06/6/2008

Math 308, Sections 301, 302, Summer 2008 Lecture 5. 06/6/2008 Math 308, Sections 301, 302, Summer 2008 Lecture 5. 06/6/2008 Chapter 3. Mathematical methods and numerical methods involving first order equations. Section 3.3 Heating and cooling of buildings. Our goal

More information

arxiv:quant-ph/ v1 29 Jun 2001

arxiv:quant-ph/ v1 29 Jun 2001 Atomic wave packet basis for quantum information Ashok Muthukrishnan an C. R. Strou, Jr. The Institute of Optics, University of Rochester, Rochester, New York 14627 (March 15, 2001) arxiv:quant-ph/0106165

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

2 The governing equations. 3 Statistical description of turbulence. 4 Turbulence modeling. 5 Turbulent wall bounded flows

2 The governing equations. 3 Statistical description of turbulence. 4 Turbulence modeling. 5 Turbulent wall bounded flows 1 The turbulence fact : Definition, observations an universal features of turbulence 2 The governing equations PART VII Homogeneous Shear Flows 3 Statistical escription of turbulence 4 Turbulence moeling

More information

POINT SOURCES OF POLLUTION: LOCAL EFFECTS AND IT S CONTROL Vol. II - Contaminant Fate and Transport Process - Xi Yang and Gang Yu

POINT SOURCES OF POLLUTION: LOCAL EFFECTS AND IT S CONTROL Vol. II - Contaminant Fate and Transport Process - Xi Yang and Gang Yu CONTAMINANT FATE AND TRANSPORT PROCESS Xi Yang and Gang Yu Department of Environmental Sciences and Engineering, Tsinghua University, Beijing, P. R. China Keywords: chemical fate, physical transport, chemical

More information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

Contaminant Fate and Transport in the Environment EOH 2122; Lecture 6, Transport in Surface Waters

Contaminant Fate and Transport in the Environment EOH 2122; Lecture 6, Transport in Surface Waters Contaminant Fate and Transport in the Environment EOH 2122; Lecture 6, Transport in Surface Waters Conrad (Dan) Volz, DrPH, MPH Bridgeside Point 100 Technology Drive Suite 564, BRIDG Pittsburgh, PA 15219-3130

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

a) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use

a) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use EAMPLE PROBLEM for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

CHAPTER 1. Chapter 1 Page 1.1. Problem 1.1: (a) Because the system is conservative, ΔE = 0 and ΔK = ΔU. G m 2. M kg.

CHAPTER 1. Chapter 1 Page 1.1. Problem 1.1: (a) Because the system is conservative, ΔE = 0 and ΔK = ΔU. G m 2. M kg. Chapter Page. CHAPTER Problem.: (a) Because the system is conservative, ΔE = 0 an ΔK = ΔU M 5.970 4 kg G 6.670 m newton R 6.370 6 m kg ΔK= v e = MmG = Mm G r R r=r so, v e M G v e.84 km R sec (b) A circular

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

Australian Journal of Basic and Applied Sciences

Australian Journal of Basic and Applied Sciences Australian Journal of Basic an Applie Sciences, 9(1) January 015, Pages: 38-45 AENSI Journals Australian Journal of Basic an Applie Sciences ISSN:1991-8178 Journal home page: www.ajbasweb.com Velocity

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information