Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Size: px
Start display at page:

Download "Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi"

Transcription

1 Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec AML3; Indian Insiue of Technology Delhi Absrac Sudy of Viscoelasic subsance in an enclosed boundary volume. This erm paper aims a he developmen of a model for sudying and modelling a creep es in a viscoelasic maerial. The basic aim of our projec is narrowed down o he applicaion of recurring numerical mehods o he soluion of coefficiens of he Prony equaion/model of he represenaion of Viscoelasic maerials. We have used he creep es daa from NASA on a maerial and used he daa for he calculaion of he maerial properies namely he above coefficiens. Lieraure Survey Viscoelasic fluids (subsances o be exac) exhibi resemblance o an elasic solid under some condiions and a viscous liquid under oher. There are he following models amongs many for he sudy of Viscoelasic subsances: MODEL MAXWELL MODEL KELVIN-VOIGT MODEL MAIN ASSUMPTION This is represened by a purely viscous damper and a purely elasic spring conneced in series. This model employs a Newonian Damper and a Hookean elasic spring. CHARACTERISTIC EQUATION dε oal d = dε D d + dε S d = σ η + dσ E d σ() = Eε() + η dε() d SYMBOLIC MODEL STANDARD LINEAR SOLID MODEL I combines a Hookean Spring and he Maxwell Model in parallel E 2 dε ( η dσ η d = E 2 d + σ E ε) E + E 2 GENERALISED MAXWELL MODEL PRONY SERIES I considers ha he relaxaion does no occur a a single ime, bu wih a disribuion of ime inervals (i.e. a variable ime disribuion). I has muliple spring-dashpos in parallel as required for he accurae represenaion of he model Takes ino accoun he linear model The equaion for creep: only wih a ε() = σ combinaion of. J() + ( ξ) dσ(ξ) dashpo and spring.

2 P a g e 2 Linear viscoelasic consiuive models are represened by simple physical models composed of springs and dashpos. The spring is he linear-elasic componen, and is consiuive equaion is σ = E. ε. The dashpo is he viscous componen, and is consiuive equaion is σ = η. ε where η is he viscosiy consan. Mahemaical Modelling (Prony Series) The consiuive equaion of a viscoelasic maerial characerises dependence upon he boh he srain and he srain rae. The equaion for creep: ε() = σ. J() + ( ξ) dσ(ξ) The coupling of σ = E. ε and σ = η. ε gives σ() = Y(). ε where Y() = E. ( p i is he i h Prony Consan τ i is he i h Prony reardaion ime consan E is he insananeous modulus of he maerial n i= p i. ( e τ i In a creep es, a creep compliance funcion, J() is defined as ε() = J(). σ n Y() = E. ( p i i=. ( e τ i)) For ime =, Y() = E and for =, y( ) = E (-Σpi). The compliance funcion is hen deermined by procedures analogous o hose described To deermine he sress sae in a viscoelasic maerial a a given ime, he deformaion hisory mus be considered. For linear viscoelasic maerials, a superposiion of herediary inegrals describes he ime dependen response a. If a specimen is load free prior o he ime = O, a which a sress, σ + σ(), is applied he srain for ime > can be represened as follows. ε() = ε. J() + J( ξ) dσ(ξ) Where J() is he compliance funcion of he maerial and dσ(ξ) is he sress rae. A similar equaion can be used for he relaxaion model o obain he sress funcion inroduced by an arbirary srain funcion σ() = ε. Y() + Y( ξ) dε(ξ) where Y() is he relaxaion funcion (Equaion 4) and dε(ξ) _is he srain rae. An example of applying herediary inegrals for a muliple loading segmen process is shown nex. )) Herediary inegrals for a muliple loading process Herediary inegrals wih Prony series kernels can be applied o model a loading process such as he one shown in Figure 2. The process is divided ino four segmens for which srain and srain rae funcions are defined. The funcions are: Navnee Saini, Mayank Goyal, Vishal Bansal

3 P a g e 3 ε. ( ) < ε() = { ε ε. ( 3 2 ) < 2 ; 2 < 3 2 < 4 ε dε ( ) < d = { < 2 ; 2 < 4 ε ( 3 2 ) 2 < 3 where ε = ε() = and =. And we have used ε =.; ε 2 =; =; =5; 2 =55; 3 =6 For a maerial wih a relaxaion funcion in he form of a Prony series (Equaion 4), he sress funcions of he loading process can be derived as follows: Sep. ( < ) Subsiue Equaions 4 and firs srain and srain rae funcions of Equaion 8 ino Equaion 7 and obain: σ() = ε. Y() + Y( ξ) dε(ξ) ] = E. ε [ξ p i. ξ + p i. τ i. e ( ξ) τi = E. ε [ p i. + p i. τ i. e τ i ] where n is he number of erms in he Prony series. To simplify he expression, n will no be shown in following equaions. Pi and "ciare he consans in he i-h erm of he Prony series. Sep 2. ( < 2 ) Using he second srain rae funcion obain: σ() = ε. Y() + Y( ξ) dε(ξ) + Y( ξ) dε(ξ) = E. ε [ p i. + p i. τ i. e ( ) τi + p i. τ i. e τ i ] Sep 3. ( 2 < 3 ) The hird srain rae funcion yields: σ() = E. ε [ p i. + p i. τ i. e ( ) τi + p i. τ i. e τ i ] + E. ε ( 3 2 ) [ p i. + p i. τ i 2 + p i. 2 + p i. τ i. e ( 2) τi Noe, he firs porion of his equaion is equal o he equaion of sep 2. Sep 4. ( 3 < 4 ) Similarly, he equaion for he fourh sep can be wrien as follows: σ() = σ 2 () E. ε ( 3 2 ) [ p i. + p i. τ i 2 + p i. 2 + p i. τ i. e ( 2) τi ] Mehod o Find he Consans of Prony Series Weighed non-linear regression: The Prony series coefficiens and reardaion ime are deermined using regression analysis. Regression analysis is basically a mehod used o deermine relaion beween he wo variables using experimenal resuls. In his mehod one independen variable is varied while keeping oher consan and he variable wih which relaion is o be found ou is ploed, and hen from graph we use differen echniques o find he relaion. One mosly used echnique is error mehod. ] Navnee Saini, Mayank Goyal, Vishal Bansal

4 P a g e 4 In his error funcion (f) wih respec o unknown consans is defined as, n F(a) = { y i y(x i, a) } 2 σ i Where xi and yi are experimenal daa, funcion y (xi, a) is he model o be fied, and σi is he sandard deviaion of measuremen error of i h daa poin. A se of unknown consans (a) will be deermined ha minimize he error funcion approximaed by is Taylor series wih he quadraic form: a. D. a F (a) = c d. a + 2 Where c is a consan and d is he gradien of F wih respec o he parameers a, which will be zero where F is minimum. Marix D is he second parial derivaive of F wih respec o he parameers a. Since each Prony erm includes wo variables (pi and ci) and since he insananeous modulus (E) mus be deermined, he oal number of variables in he regression is 2n+l. Since viscoelasic effecs are mos significan a he beginning of he relaxaion period, he fiing error in his region is significan. Since he percenage of he number of daa poins a he beginning of he relaxaion period is less, he error funcion χ 2 is dominaed by a long uniform ail region of he relaxaion period. To reduce he error and improve he fi a beginning of he relaxaion period, he weigh funcion [w = ; < σ <] is used. σ Numerical Soluions Working Non Linear for Regression Mehod The specimen (2 in. x 2 in. x.768 in.) was loaded in 22 seconds o a maximum deflecion of.3 in. a he middle of he span. Then he deflecion was held for,7 seconds (Figure 3). Since no analyical mehod exiss o form he weigh funcion, a rial-and-error mehod based on maerial properies was used o obain a fi curve. The viscoelasic sress decays exponenially in he relaxaion es. The viscoelasic sress decays exponenially in he relaxaion es, herefore he viscoelasic effec is more significan during he loading period and a he beginning (< 3seconds) of he holding period. The number of daa poins in hese periods (<) is much less han he number of daa poins in ail region of he relaxaion period (>5). The error funcion (F) will be dominaed by along uniform ail region of he relaxaion period if a uniform weigh funcion is applied. Therefore, a piecewise weigh funcion was used o obain beer fis for hese periods and improve he accuracy of he regression. Figure4 shows he load relaxaion a he beginning of he process. The dos represen he es daa. Three regression resuls are shown. The dash-do curve is he resul of a regression analysis wihou he weigh funcion (w/o WF) for a wo erm Prony series. Navnee Saini, Mayank Goyal, Vishal Bansal

5 P a g e 5 The long-dash curve is he resul for a wo-erm Prony series wih weigh funcion number (WF) shown below: 5 < 22 w = { < 7 The shor-dash curve is he resul for a hree-erm Prony series wih weigh funcion 2 (WF2) as follows:. < 22 w = { 6 ( 2 )22 < 7 As is well known, he regression daa fis beer in a paricular region if he relaive weigh facor (weigh facor / sum of weigh facors) in ha region is greaer han he average value. In Figure 4, he curves wih weighed funcions were closer o he daa poins near he beginning of he relaxaion period. Since he weigh facor a = 22 seconds of funcion 2 (=266) is greaer han he value of he weigh funcion (=), he curve of funcion 2 is a beer mach o he daa han curve of funcion a he beginning of relaxaion. However, as shown in he Figure3, weigh funcion2 does no fi he daa as well as he oher curves afer seconds. I appears ha funcion2 was over weighed a he beginning of he relaxaion and he relaive weigh facor a he oher region was oo small. The Prony consans of he regression are shown in Table. The resuls showed ha a weigh funcion should be seleced based on he maerial properies and es daa disribuion. Properly selecing he weigh funcion in he mos significan region can improve accuracy of he regression. Navnee Saini, Mayank Goyal, Vishal Bansal

6 P a g e 6 Conclusions In his we came wih echnique o deermine he Prony coefficiens ha can represen viscoelasic sress and srain equaion. The herediary inegral mehod was employed o obain an analyical represenaion of maerial response when i is subjeced o rae dependen loading. The analyical represenaion was used in a nonlinear regression analysis, wih measured daa, o evaluae he Prony series consans. Several regression analyses were performed using differen weigh funcions. Noe, he daa of similar lengh in ime had been analyzed for loading and relaxaion regions. Oher weighing funcions may be needed for differen loading schedules. The mehod could be used o provide highly accurae represenaion of he maerial behavior in he rae dependen loading region. Since once afer consans are deermined, he response of a viscoelasic maerial for oher loading schedules can also be calculaed. I can also be used o deermine equaion for oher loading schedules where we can keep relaxaion ime from our choice. References Deermining a Prony Series for a Viscoelasic Maerial from Time Varying Srain Daa (26); Tzikang Chen; U.S. Army Research Laboraory, Vehicle Technology Direcorae, Langley Research Cener, Virginia; Naional Aeronauics and Space Adminisraion. Microfluidic Flows of Viscoelasic Fluids; Monica S. N. Oliveira, Manuel A. Alves, Fernando T. Pinho; Transpor and Mixing in Laminar Flows: From Microfluidics o Oceanic Currens. Engineering Viscoelasiciy (2); David Roylance; Deparmen of Maerials Science and Engineering, Massachuses Insiue of Technology. Aachmens. Sress.m 2. Srain.m 3. Presenaion.pp Navnee Saini, Mayank Goyal, Vishal Bansal

Polymer Engineering (MM3POE)

Polymer Engineering (MM3POE) Polymer Engineering (MM3POE) VISCOELASTICITY hp://www.noingham.ac.uk/~eazacl/mm3poe Viscoelasiciy 1 Conens Wha is viscoelasiciy? Fundamenals Creep & creep recovery Sress relaxaion Modelling viscoelasic

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

ES240 Solid Mechanics Fall 2007

ES240 Solid Mechanics Fall 2007 ES4 Solid Mechanics Fall 7 Viscoelasiciy References NG McCrum, CP Buckley, and CB Bucknall, Principles of Polymer Engineering, nd ediion, Oxford Universiy Press, 997 A good balance of heory and applicaion

More information

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/

More information

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification IOP Conference Series: Maerials Science and Engineering PAPE OPEN ACCESS New effecive moduli of isoropic viscoelasic composies. Par I. Theoreical jusificaion To cie his aricle: A A Sveashkov and A A akurov

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

MIT Amorphous Materials

MIT Amorphous Materials MIT 3.071 Amorphous Maerials 7: Viscoelasiciy and Relaxaion Juejun (JJ) Hu 1 Afer-class reading lis Fundamenals of Inorganic Glasses Ch. 8, Ch. 13, Appendix A Inroducion o Glass Science and Technology

More information

Time series Decomposition method

Time series Decomposition method Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Analytic nonlinear elasto-viscosity of two types of BN and PI rubbers at large deformations

Analytic nonlinear elasto-viscosity of two types of BN and PI rubbers at large deformations Bulgarian Chemical Communicaions, Volume 48, Special Issue E (pp. 59-64) 016 Analyic nonlinear elaso-viscosiy of wo ypes of BN and PI rubbers a large deformaions K. B. Hadjov, A. S. Aleksandrov, M. P.

More information

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3 The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped

More information

MECHANICS OF MATERIALS Poisson s Ratio

MECHANICS OF MATERIALS Poisson s Ratio Poisson s Raio For a slender bar subjeced o axial loading: ε x x y 0 The elongaion in he x-direcion i is accompanied by a conracion in he oher direcions. Assuming ha he maerial is isoropic (no direcional

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

INDEX. Transient analysis 1 Initial Conditions 1

INDEX. Transient analysis 1 Initial Conditions 1 INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Curling Stress Equation for Transverse Joint Edge of a Concrete Pavement Slab Based on Finite-Element Method Analysis

Curling Stress Equation for Transverse Joint Edge of a Concrete Pavement Slab Based on Finite-Element Method Analysis TRANSPORTATION RESEARCH RECORD 155 35 Curling Sress Equaion for Transverse Join Edge of a Concree Pavemen Slab Based on Finie-Elemen Mehod Analysis TATSUO NISHIZAWA, TADASHI FUKUDA, SABURO MATSUNO, AND

More information

CHAPTER 6: FIRST-ORDER CIRCUITS

CHAPTER 6: FIRST-ORDER CIRCUITS EEE5: CI CUI T THEOY CHAPTE 6: FIST-ODE CICUITS 6. Inroducion This chaper considers L and C circuis. Applying he Kirshoff s law o C and L circuis produces differenial equaions. The differenial equaions

More information

The equation to any straight line can be expressed in the form:

The equation to any straight line can be expressed in the form: Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he

More information

Role of Material Functions in Rheological Analysis. Role of Material Functions in Rheological Analysis

Role of Material Functions in Rheological Analysis. Role of Material Functions in Rheological Analysis Role of Maerial Funcions in Rheological Analysis QUALIY CONROL compare wih oher in-house daa on qualiaie basis unknown maerial QUALIAIVE ANALYSIS compare daa wih lieraure repors on arious fluids conclude

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

5.2. The Natural Logarithm. Solution

5.2. The Natural Logarithm. Solution 5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

Average Number of Lattice Points in a Disk

Average Number of Lattice Points in a Disk Average Number of Laice Poins in a Disk Sujay Jayakar Rober S. Sricharz Absrac The difference beween he number of laice poins in a disk of radius /π and he area of he disk /4π is equal o he error in he

More information

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

( ) = b n ( t) n  (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2. Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion

More information

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel, Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

1 Evaluating Chromatograms

1 Evaluating Chromatograms 3 1 Evaluaing Chromaograms Hans-Joachim Kuss and Daniel Sauffer Chromaography is, in principle, a diluion process. In HPLC analysis, on dissolving he subsances o be analyzed in an eluen and hen injecing

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Viscoelastic Catenary

Viscoelastic Catenary Viscoelasic Caenary Anshuman Roy 1 Inroducion This paper seeks o deermine he shape of a hin viscoelasic fluid filamen as i sags under is own weigh. The problem is an exension of he viscous caenary [1]

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

Modeling and Simulation of Viscoelastic Flows in Polymer Processing with Integral Constitutive Equations Evan Mitsoulis

Modeling and Simulation of Viscoelastic Flows in Polymer Processing with Integral Constitutive Equations Evan Mitsoulis Modeling and Simulaion of Viscoelasic Flows in Polymer Processing wih Inegral onsiuive Equaions Evan Misoulis AMP - R&D ompuer-aided Maerials Processing Rheology & Design School of Mining Engineering and

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

2. For a one-point fixed time method, a pseudo-first order reaction obeys the equation 0.309

2. For a one-point fixed time method, a pseudo-first order reaction obeys the equation 0.309 Chaper 3. To derive an appropriae equaion we firs noe he following general relaionship beween he concenraion of A a ime, [A], he iniial concenraion of A, [A], and he concenraion of P a ime, [P] [ P] Subsiuing

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Chapter 13 Homework Answers

Chapter 13 Homework Answers Chaper 3 Homework Answers 3.. The answer is c, doubling he [C] o while keeping he [A] o and [B] o consan. 3.2. a. Since he graph is no linear, here is no way o deermine he reacion order by inspecion. A

More information

ON LINEAR VISCOELASTICITY WITHIN GENERAL FRACTIONAL DERIVATIVES WITHOUT SINGULAR KERNEL

ON LINEAR VISCOELASTICITY WITHIN GENERAL FRACTIONAL DERIVATIVES WITHOUT SINGULAR KERNEL Gao F. e al.: On Linear Viscoelasiciy wihin General Fracional erivaives... THERMAL SCIENCE: Year 7 Vol. Suppl. pp. S335-S34 S335 ON LINEAR VISCOELASTICITY WITHIN GENERAL FRACTIONAL ERIVATIVES WITHOUT SINGULAR

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

ENV 6015 Solution to Mixing Problem Set

ENV 6015 Solution to Mixing Problem Set EN 65 Soluion o ixing Problem Se. A slug of dye ( ) is injeced ino a single ank wih coninuous mixing. The flow in and ou of he ank is.5 gpm. The ank volume is 5 gallons. When will he dye concenraion equal

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

04. Kinetics of a second order reaction

04. Kinetics of a second order reaction 4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius

More information

APPENDIX H FRACTURE PROPERTIES OF ASPHALT MIXTURES

APPENDIX H FRACTURE PROPERTIES OF ASPHALT MIXTURES APPENDIX H FRACTURE PROPERTIES OF ASPHALT MIXTURES LIST OF FIGURES Figure Page H- Loading ime under sress free emperaure... H-2 H-2 Mehod o evaluae he m mix of fracure properies for hermal case... H-3

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS 8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,

More information

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount Adv. Sudies Theor. Phys., Vol. 7, 3, no., 3-33 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/.988/asp.3.3999 Accurae RS Calculaions for Periodic Signals by Trapezoidal Rule wih he Leas Daa Amoun Sompop Poomjan,

More information

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he

More information

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Kragujevac J. Sci. 3 () 7-4. UDC 53.5:536. 4 THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Hazem A. Aia Dep. of Mahemaics, College of Science,King Saud Universiy

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

1.6. Slopes of Tangents and Instantaneous Rate of Change

1.6. Slopes of Tangents and Instantaneous Rate of Change 1.6 Slopes of Tangens and Insananeous Rae of Change When you hi or kick a ball, he heigh, h, in meres, of he ball can be modelled by he equaion h() 4.9 2 v c. In his equaion, is he ime, in seconds; c represens

More information

Finite Element Analysis of Structures

Finite Element Analysis of Structures KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Testing the Random Walk Model. i.i.d. ( ) r

Testing the Random Walk Model. i.i.d. ( ) r he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:

EE650R: Reliability Physics of Nanoelectronic Devices Lecture 9: EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model: Dynamic Economeric Models: A. Auoregressive Model: Y = + 0 X 1 Y -1 + 2 Y -2 + k Y -k + e (Wih lagged dependen variable(s) on he RHS) B. Disribued-lag Model: Y = + 0 X + 1 X -1 + 2 X -2 + + k X -k + e

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Chapter 2 Constitutive Equations for Isotropic and Anisotropic Linear Viscoelastic Materials

Chapter 2 Constitutive Equations for Isotropic and Anisotropic Linear Viscoelastic Materials Chaper 2 Consiuive Equaions for Isoropic and Anisoropic Linear Viscoelasic Maerials Jacek J. Skrzypek and Arur W. Ganczarski Absrac In case of isoropic maerial symmery, he elasic-viscoelasic correspondence

More information