6. Well-Stirred Reactors III
|
|
- Ella Carroll
- 6 years ago
- Views:
Transcription
1 6. Well-Stirred Reactors III Reactors reaction rate or reaction velocity defined for a closed system of uniform pressure, temperature, and composition situation in a real reactor is usually quite different from this ideal situation types of reactors: mass exchange: batch reactor flow reactor heat exchange: adiabatic isothermal mechanical variables: constant pressure constant volume residence time: unique distribution (ideal: exponential) operation: transient stationary residence time = time spent in the reactor by each volume element of the reacting mass ideal reactor for the direct measurement of reaction rates is a flow, isothermal, constant-volume reactor Adv PChem 6.1
2 operating in the stationary state with such thorough mixing that the composition is the same everywhere in the reactor: continuous stirred tank reactor (CSTR) [continuous = continuous-flow] c i V T c i,f volumetric flow rate q; units L/s feed concentration of species i : c i,f ; units M reactor volume V ; units L reaction velocity v; units M/s reaction in the reactor: ν i J i = 0 mass-balance equation for each species (i = 1,2,..., N) d(v c i ) dt = qc i,f qc i + ν i V v(c 1,c 2,...,c N,T ) Adv PChem 6.2
3 divide by V : dc i dt = q V ( ci,f c i ) + νi v(c 1,c 2,...,c N,T ) define α q V flow rate units s 1 dc i dt = α( c i,f c i ) + νi v(c 1,c 2,...,c N,T ) stationary state 0 = dc i dt = α( c i,f c i ) + νi v(c 1,c 2,...,c N,T ) ν i v(c 1,c 2,...,c N,T ) = α ( c i c i,f ) v(c 1,c 2,...,c N,T ) = α ν i ( ci c i,f ) tracer experiment: CSTR initially filled with dye, then Adv PChem 6.3
4 flow distilled water dc dt = αc c(t) = c(0)exp( αt) distribution of residence times of particles in the reactor: exit age distribution E(t) E(t)dt is the probability that a particle exits the reactor between t and t + dt the number of tracer particles leaving the reactor between t and t + dt is dn (t) = N A qc(t)dt the number tracer particles initially in the reactor is N 0 = N A V c(0) which equals N 0 = 0 dn (t) = 0 N A qc(t)dt = N A q 0 c(t)dt Adv PChem 6.4
5 E(t)dt is the fraction of particles exiting between t and t + dt E(t)dt = or E(t)dt = So dn (t) N 0 = N Aqc(t)dt N A V c(0) = αc(t)dt c(0) dn (t) N 0 = N Aqc(t)dt N A q 0 c(t)dt = c(t)dt 0 c(t)dt E(t) = αc(t) c(0) = c(t) 0 c(t)dt Adv PChem 6.5
6 For a CSTR we find E(t) = αc(t) c(0) = αc(0)exp( αt) c(0) = αexp( αt) or c(t) E(t) = 0 c(t)dt c(0) exp( αt) = 0 c(0)exp( αt)dt exp( αt) = [ α 1 ] exp( αt) 0 = exp( αt) 1 α = αexp( αt) residence times in a CSTR are exponentially distribut- Adv PChem 6.6
7 ed τ = t = = 0 = α te(t)dt tαexp( αt)dt t exp( αt)dt x n exp( ax)dx = n! a n+1, a > 0 τ = α 1 α 2 τ = 1 α τ = residence time of CSTR [somewhat of an abuse of language] = inverse of flow rate Adv PChem 6.7
8 Gray-Scott model Gray, P. & Scott, S. K.: Sustained oscillations and other exotic patterns of behavior in isothermal reactions, J. Phys. Chem., 89, (1985), doi.org/ /j100247a009; Gray, P. & Scott, S. K.: Chemical Oscillations and Instabilities, Clarendon Press, (1990) U + 2V V k 1 k 1 3V (R1) k 2 k 2 W (R2) rate equations du dt = k 1u v 2 + k 1 v 3 + α ( u 0 u ) dv dt = k 1u v 2 k 1 v 3 k 2 v + k 2 w + α ( v 0 v ) dw dt = k 2 v k 2 w + α ( w 0 w ) Adv PChem 6.8
9 non-dimensionalize the rate equations u = u u, v = v 0 u, w = w 0 u 0 t = t k 1 u 2 0 k = k 2 k 1 u 2 0 α = α k 1 u 2 0 η 1 = k 1 k 1 η 2 = k 2 k 2 non-dimensionalized rate equations du dt = uv 2 + η 1 v 3 + α(1 u) dv dt = uv 2 η 1 v 3 kv + kη 2 w + α(v 0 v) dw dt = kv kη 2w + α(w 0 w) (1a) (1b) (1c) Adv PChem 6.9
10 adding the three rate equations we find d(u + v + w) dt = α ( 1 + v 0 + w 0 u v w ) define z(t) u(t) + v(t) + w(t), z v 0 + w 0 then dz dt = α( z 0 z ) thermodynamic equilibrium can only be attained if the net flow of matter through the reactor vanishes: α = 0 thermodynamic equilibrium state: α = 0, u = 0, v = 0, ẇ = 0 0 = u eq v 2 eq + η 1 v 3 eq 0 = u eq v 2 eq η 1 v 3 eq kv eq + kη 2 w eq Adv PChem 6.10
11 0 = kv eq kη 2 w eq either v eq = η 2 w eq u eq v 2 eq = η 1 v 3 eq v eq = 1 η 1 u eq [if v eq 0] w eq = 1 η 1 η 2 u eq u eq arbitrary or v eq = 0 w eq = 0 u eq arbitrary the second equilibrium state is a somewhat pathological feature of the reaction scheme: if initially no Adv PChem 6.11
12 V and W is present in the reactor, then no reaction will occur for α = 0, there is a further condition dz dt = 0 = z(t) = const = z(0) u eq + v eq + w eq = u(0) + v(0) + w(0) = z(0) therefore u eq + 1 η 1 u eq + 1 η 1 η 2 u eq = z(0) u eq = otherwise z(0) η η 1 η 2 if v(0) 0 or w(0) 0 u eq = u(0) the interesting case is the non-trivial first equilibrium state u eq = z(0) η η 1 η 2 Adv PChem 6.12
13 v eq = 1 η 1 u eq w eq = 1 η 1 η 2 u eq as α is increased from zero, the system is driven away from thermodynamic equilibrium determine nonequilibrium steady states and other asymptotic states of the reactor exploit the first integral of the CSTR dz dt = α( z 0 z ) z(t) = z(0)exp( αt) + z 0 [ 1 exp( αt) ] z(t) t z 0 in other words, the total concentration relaxes to the total concentration in the feed streams in experiments, one usually takes t to mean t = 4/α = 4τ; after that time, the initial condition has essentially been washed out of the reactor Adv PChem 6.13
14 therefore we will assume from now on that z(0) = z 0, which implies z(t) = z 0 eliminate w: w = z 0 u v du dt = uv 2 + η 1 v 3 + α(1 u) dv dt = uv 2 η 1 v 3 kv + kη 2 ( z0 u v ) + α(v 0 v) study the dynamical behavior of Gray-Scott model under the following simplifying assumptions: (i) the back reactions are negligible, irreversible Gray- Scott scheme η 1 1, η 2 1 Set η 1 = η 2 = 0 in the rate equations du dt = uv 2 + α(1 u) Adv PChem 6.14
15 dv dt = uv 2 kv + α(v 0 v) (ii) the autocatalytic intermediate V is not fed, v 0 = 0 du dt = uv 2 + α(1 u) dv dt = uv 2 kv αv find the steady states: 0 = u v 2 + α(1 u) 0 = u v 2 (k + α)v = v [ u v (k + α) ] consequently either u = 1, v = 0 [trivial steady state, no reaction] or u = k + α v Adv PChem 6.15
16 0 = k + α ( v v 2 + α 1 k + α ) v 0 = (k + α)v 2 + αv α(k + α) v 2 α k + α v + α = 0 v 2,3 = 1 α 4(k + α)2 1 ± 1 2 k + α α these two roots exist if 4(k + α)2 1 α 0 they coalesce into one double root at 4(k + α) 2 α = 1 4k 2 + 8kα + 4α 2 = α ( ) α 2 + α ± = 2k 1 4 ( k 1 8 α + k 2 = 0 ) ± k k k2 Adv PChem 6.16
17 α ± = ( k 1 8 ) ± k = the two non-trivial steady states (u, v) 2,3 appear and disappear at α and α + of course α ± must be a non-negative real number; = k 1 8 < 0 and k > 0 k < 1 8 and k < 1 16 consequently, if k < 1 16 then over the interval (α,α + ) of the flow rate the two non-trivial steady states exist note: α v 2,3 (α ) = k + α Adv PChem 6.17
18 α + v 2,3 (α + ) = k + α + = these two new branches of steady states do not appear via a bifurcation from the branch of trivial steady states, v = 0 the branches of non-trivial steady states form an isola for the non-trivial steady states u v (k + α) = 0 u v = k + α u = k + α v δ u 2,3 = u 2,3 = 4(k + α)2 α 1 2 k + α α 2(k + α)2 α [0 δ 1 see above] k + α [ 1 1 δ 1 ] 1 1 δ Adv PChem 6.18
19 u 2,3 = 1 4(k + α) α 1 1 δ u 2,3 = 1 2 δ δ u 2,3 = 1 ( 1 )( 1 δ 1 + ) 1 1 δ δ u 2,3 = 1 [ 1 ± ] 1 δ 2 steady states of the irreversible Gray-Scott model: (u 1, v 1 ) = (0,1) ( 1 (u 2, v 2 ) = 2 ( 1 (u 3, v 3 ) = 2 [ δ [ 1 ] 1 δ ], 1 α 2 k + α, 1 α 2 k + α [ 1 ] ) 1 δ [ 1 + ] ) 1 δ stability of the steady states Adv PChem 6.19
20 Jacobian ( ) v 2 α 2u v J = v 2 2u v (k + α) for the trivial steady state, (u 1, v 1 ) = (0,1): J = ( ) α 0 0 (k + α) the Jacobian is diagonal and the eigenvalues are λ 1 = α λ 2 = (k + α) the trivial steady state is stable for all flow rates this implies that a small initial amount of the autocatalyst V will decay to zero non-trivial steady states J = ( v 2 α 2(k + α) v 2 k + α ) Adv PChem 6.20
21 sign structure J = ( ) + + the Gray-Scott model is a cross activator-inhibitor scheme; V is the activator T = trj = v 2 + k = detj = ( v 2 + α ) (k + α) + 2(k + α) v 2 = (k + α) ( v 2 α + 2v 2) = (k + α) ( v 2 α ) ( for[ the lower branch of the isola, (u 2, v 2 ) = ] 1 δ, 1 2 k + α [ 1 ]) 1 δ : α = (k + α) ( 1 α 2 [ 4 (k + α) 2 1 α 2 [ 4 (k + α) 2 1 ] 2! 1 δ α < 0 1 ) ] 2 1 δ α Adv PChem 6.21
22 1 α 2 [ 4 (k + α) 2 1 ] 2! 1 δ < α [ 1 ] 2! 4(k + α)2 1 δ < = δ α δ + 1 δ! < δ 2 2δ! < 2 1 δ 1 δ! < 1 δ true for 0 < δ < 1 the determinant is negative for all steady states on the lower branch of the isola and becomes zero at α ± = these states are saddle points and are unstable for ( [ the upper branch of the isola, (u 3, v 3 ) = 12 1 ] 1 δ, 1 2 k + α [ 1 + ]) 1 δ : α ( 1 α 2 [ = (k + α) 4 (k + α) ) ] 2 1 δ α 1 α 2 [ 4 (k + α) ] 2! 1 δ α > 0 Adv PChem 6.22
23 1 α 2 [ 4 (k + α) ] 2! 1 δ > α [ 1 + ] 2! 4(k + α)2 1 δ > = δ α δ + 1 δ! > δ 2 1 δ! > 2δ 2 1 δ! > δ 1 true for 0 < δ < 1 the determinant is positive for all steady states on the upper branch of the isola and becomes zero at α ± = stability is determined by T for α near α +, T is negative; the steady states are stable for α near α, T is positive; the steady states are unstable The upper branch of the isola undergoes a Hopf bifur- Adv PChem 6.23
24 cation, T = 0: v 2 + k = 0 v 2 = k v = k 1 α [ 1 + ] 1 δ = k 2 k + α 1 α 4(k + α) = k 2 k + α α 4(k + α) = 2(k + α) k α α 4(k + α)2 1 = 2(k + α) k 1 α α [ 4(k + α)2 2(k + α) k 1 = 1 α α ] 2 4(k + α)2 1 α = 4(k + α)2 k α 2 4(k + α) k α + 1 Adv PChem 6.24
25 α(k + α) = (k + α)k α k ( α 2 + 2k ) k + k 2 = 0 α 1,2 = 1 ( 2k ) 1 ( k ± 2k ) 2 k k [ k 2k ± 4k 2 4k ] k + k 4k 2 α 1,2 = 1 2 [ k 2k ± α 1,2 = 1 2 k 4k k ] α 2 is the acceptable root; α 1 is a parasite root [the original equation is not zero at that solution]; so α H = 1 2 [ k 2k k ( 1 4 ) ] k as the flow rate α is decreased from α +, the upper branch of the isola undergoes a Hopf bifurcation at α H examples for two values of k: k = = 0.9 (1/16), α = , α + = , α H = Adv PChem 6.25
26 v α k = = 1/32, α = , α + = , α H = v α Adv PChem 6.26
27 the dynamical behavior of two-variables systems is restricted since trajectories cannot cross the possible asymptotic states are stationary states and simple oscillations (simple limit cycles) additional dynamical behavior and asymptotic states occur in systems with three or more variables consider the following reaction scheme which is related to the Gray-Scott model and known as the autocalator [Merkin, J. H.; Needham, D. J. & Scott, S. K.: Oscillatory Chemical Reactions in Closed Vessels, Proc. R. Soc. Lond. A, 406, (1986), royalsocietypublishing.org/content/406/1831/299. abstract] P U + 2V V U k 0 U (R3) k 1 3V (R4) k 2 W (R5) k u V (R6) (R7) Adv PChem 6.27
28 Taking into account that the precursor P decays with a first-order rate law, the nondimensionalized rate equations are given by du dt = µexp( ɛt) uv 2 u dv dt = uv 2 + u v this describes a closed system and provides a simple model to study transient phenomena with the pool chemical assumption, supply of the precursor, inexhaustible du dt = µ P uv 2 u dv dt = uv 2 + u v two-variable system, no new dynamical behavior Showalter and coworkers introduced a modification of the pool chemical autocatalator model which leads to a three-variable scheme [Peng, B.; Scott, S. K. Adv PChem 6.28
29 & Showalter, K.: Period doubling and chaos in a three-variable autocatalator, J. Phys. Chem., 94, (1990), /j100376a014] P U + 2V V U P + W k 0 U (R8) k 1 3V (R9) k 2 W (R10) k u V (R11) k c U + W (R12) W k 3 Z (R13) nondimensionalized rate equations du dt = κ + µw uv 2 u σ dv dt = uv 2 + u v Adv PChem 6.29
30 δ dw dt = v w this scheme dispalys simple and complex oscillations and a period-doubling sequence to chaos: aperiodic deterministic motion; strange attractor; sensitive dependence on initial conditions Lorenz model [Lorenz, E. N.: Simplified dynamic equations applied to the rotating-basin experiments, J. Atmos. Sci., 19, (1962), %3C0039%3ASDEATT%3E2.0.CO%3B2] simplified model of Rayleigh-Bénard convection; shallow fluid layer between two plates, upper plate has uniform temperature of T and lower plate a uniform temperature of T + T, T > 0; if T is sufficiently small, then no fluid motion and the temperature profile between the lower and upper plate is linear; if T > T c, convective motion in the form of rolls nondimensionalized evolution equations dx dt = σ(y x) Adv PChem 6.30
31 dy = x(r z) y dt dz dt = x y bz x velocity amplitude of convective motion, y difference in temperature between rising and descending currents, z amplitude of distortion of the temperature profile from linearity; σ, r and b are positive parameters steady states 0 = σ(y x) 0 = x(r z) y 0 = x y bz The first equation implies that y = x Substitute into second and third equation 0 = x(r z) x Adv PChem 6.31
32 0 = x 2 bz 0 = x(r 1 z) = x = 0 or z = r 1 Case 1: x = 0 = y = 0 = z = 0, no convection (trivial steady state) (x, y, z) 0 = (0,0,0) Case 2: z = r 1 = 0 = x 2 b(r 1) x = ± b(r 1) y = ± b(r 1) ( (x, y, z) 1,2 = ± b(r 1),± ) b(r 1),r 1 steady states 1 and 2 do not exist for r < 1 linear stability analysis Adv PChem 6.32
33 Jacobian σ σ 0 J = z + r 1 x y x b stability of the trivial steady state σ σ 0 J = r b eigenvalues det ( ) σ λ σ 0 J λi 3 = r 1 λ b λ = σ λ σ r 1 λ ( b λ) = [ λ 2 + (σ + 1)λ + σ(1 r ) ] ( b λ) = 0 Adv PChem 6.33
34 λ 0;1 = b [ λ 0;2,3 = 1 (σ + 1) ± 2 [ = 1 2 (σ + 1) ± ] (σ + 1) 2 4σ(1 r ) ] (σ 1) 2 + 4σr λ 0,1 and λ 0,3 are always real and negative λ 0,2 is always real and crosses 0 at r = r c : σ + 1 = (σ 1) 2 + 4σr c (σ + 1) 2 = (σ 1) 2 + 4σr c σ 2 + 2σ + 1 = σ 2 2σ σr c 4σ = 4σr c r c = 1 as expected at r = r c = 1 the quiescent state state becomes unstable and convective rolls appear Adv PChem 6.34
35 stability analysis of the convective branches (x, y, z) 1,2 : J = σ σ 0 r r 1 b(r 1) ± b(r 1) ± b(r 1) b σ σ 0 = 1 1 b(r 1) ± b(r 1) ± b(r 1) b eigenvalues det ( σ λ σ 0 ) J λi 3 = 1 1 λ b(r 1) ± b(r 1) ± b(r 1) b λ = 0 characteristic polynomial λ 3 + c 1 λ 2 + c 2 λ + c 3 = 0 c 1 = 1 + σ + b > 0 Adv PChem 6.35
36 c 2 = b(r + σ) > 0 c 3 = 2σb(r 1) > 0 linear stability analysis for systems with n 3 variables: characteristic polynomial ( 1) n det ( J λi n ) = 0 λ n + c 1 λ n 1 + c 2 λ n c n 1 λ + c n = 0 Routh-Hurwitz theorem: All roots of the characteristic polynomial have a negative real part, if and only if c 1 c 3 1 c 2 c 4 0 c l = 1 c c 2 > 0, l = 1,...,n 0 0 c l together with the condition c n > 0 Adv PChem 6.36
37 Set c l = 0 in l for l > n generic bifurcations of steady states in n-variable systems, n > 2: stationary bifurcation, where a real eigenvalues passes through zero, and the Hopf bifurcation, where a pair of complex conjugate eigenvalues crosses the imaginary axis stationary bifurcation: c n = 0 Hopf bifurcation: n 1 = 0 c n > 0 l > 0, l = 1,...,n 2 Lorenz model, the convective branches (x, y, z) 1,2 : 1 = c 1 = 1 + σ + b > 0 2 = c 1 c 3 1 c = c 1c 2 c 3 = b [r (1 + b σ) + σ(3 + b + σ)] 2 Adv PChem 6.37
38 c 1 c 3 c 5 c 1 c 3 0 c 3 = 1 c 2 c 4 = 1 c 2 0 = c 1 c c = c c 1 c 3 0 c 1 c 3 2 since c 3 = 2σb(r 1) > 0, the convective states (x, y, z) 1,2 cannot undergo a stationary bifurcation the convective states (x, y, z) 1,2 undergo a Hopf bifurcation, if σ > b + 1, at r H = σ(σ + b + 3) σ b 1 where 2 = 0 and 1 > 0 and c 3 > 0 it turns out that the Hopf bifurcation is subcritical all trajectories remain bounded above r H ; the following figures illustrate the behavior of the Lorenz model for σ = 10.0, b = 8/3 and r = 28.0 [the Hopf bifurcation occurs at r H = 24.74] initial condition: (x(0), y(0), z(0)) = (x, y, z) 1 + (0.001, 0.001,0.001) Adv PChem 6.38
39 time series x y t t z t Adv PChem 6.39
40 the asymptotic state is a strange attractor z x sensitive dependence on initial conditions: at t = 500 we have (x, y, z) = ( , , ); start a second trajectory nearby (x, y, z )(500) = ( , , ); figure shows the Euclidean distance [x(t) (t) = x (t) ] 2 [ + y(t) y (t) ] 2 [ + z(t) z (t) ] 2 Adv PChem 6.40
41 t main routes to chaos in chemical systems: (i) perioddoubling, (ii) mixed-mode oscillations, (iii) intermittency (regular oscillations are interrupted by chaotic bursts of random duration), (iv) quasiperiodicity (two incommensurate frequencies) Adv PChem 6.41
7. Well-Stirred Reactors IV
7. Well-Stirred Reactors IV The Belousov-Zhabotinsky Reaction: Models and Experiments Oregonator [based on the FKN mechanism; Field, R. J. & Noyes, R. M.: Oscillations in chemical systems. IV. Limit cycle
More informationSemi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical
ANZIAM J. 53 (EMAC211) pp.c511 C524, 212 C511 Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical M. R. Alharthi 1 T. R. Marchant 2 M. I.
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationInstabilities In A Reaction Diffusion Model: Spatially Homogeneous And Distributed Systems
Applied Mathematics E-Notes, 10(010), 136-146 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Instabilities In A Reaction Diffusion Model: Spatially Homogeneous And
More informationQuadratic autocatalysis in an extended continuousflow stirred tank reactor (ECSTR)
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 2016 Quadratic autocatalysis in an extended continuousflow
More information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationScenarios for the transition to chaos
Scenarios for the transition to chaos Alessandro Torcini alessandro.torcini@cnr.it Istituto dei Sistemi Complessi - CNR - Firenze Istituto Nazionale di Fisica Nucleare - Sezione di Firenze Centro interdipartimentale
More informationLecture 3 : Bifurcation Analysis
Lecture 3 : Bifurcation Analysis D. Sumpter & S.C. Nicolis October - December 2008 D. Sumpter & S.C. Nicolis General settings 4 basic bifurcations (as long as there is only one unstable mode!) steady state
More information6.2 Brief review of fundamental concepts about chaotic systems
6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationChaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad
Chaos Dr. Dylan McNamara people.uncw.edu/mcnamarad Discovery of chaos Discovered in early 1960 s by Edward N. Lorenz (in a 3-D continuous-time model) Popularized in 1976 by Sir Robert M. May as an example
More informationModelling and Mathematical Methods in Process and Chemical Engineering
Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series 3 1. Population dynamics: Gendercide The system admits two steady states The Jacobi matrix is ẋ = (1 p)xy k 1 x ẏ
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationTHE STATIC BIFURCATION DIAGRAM FOR THE GRAY SCOTT MODEL
International Journal of Bifurcation and Chaos, Vol. 11, No. 9 2001) 2483 2491 c World Scientific Publishing Company THE STATIC BIFURCATION DIAGRAM FOR THE GRAY SCOTT MODEL RODICA CURTU Department of Mathematical
More informationCubic autocatalytic reaction diffusion equations: semi-analytical solutions
T&T Proof 01PA003 2 January 2002 10.1098/rspa.2001.0899 Cubic autocatalytic reaction diffusion equations: semi-analytical solutions By T. R. Marchant School of Mathematics and Applied Statistics, The University
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationMATH 415, WEEK 12 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior
MATH 415, WEEK 1 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior 1 Higher-Dimensional Systems Consider the following system of differential equations: dx = x y dt dy dt = xy y dz dt
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationComplex Patterns in a Simple System
Complex Patterns in a Simple System arxiv:patt-sol/9304003v1 17 Apr 1993 John E. Pearson Center for Nonlinear Studies Los Alamos National Laboratory February 4, 2008 Abstract Numerical simulations of a
More informationPOPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL. If they co-exist in the same environment:
POPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL Next logical step: consider dynamics of more than one species. We start with models of 2 interacting species. We consider,
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
More informationChemical Reaction Engineering Lecture 5
Chemical Reaction Engineering g Lecture 5 The Scope The im of the Course: To learn how to describe a system where a (bio)chemical reaction takes place (further called reactor) Reactors Pharmacokinetics
More informationTHE SYNCHRONIZATION OF TWO CHAOTIC MODELS OF CHEMICAL REACTIONS
ROMAI J., v.10, no.1(2014), 137 145 THE SYNCHRONIZATION OF TWO CHAOTIC MODELS OF CHEMICAL REACTIONS Servilia Oancea 1, Andrei-Victor Oancea 2, Ioan Grosu 3 1 U.S.A.M.V., Iaşi, Romania 2 Erasmus Mundus
More informationLaurette TUCKERMAN Rayleigh-Bénard Convection and Lorenz Model
Laurette TUCKERMAN laurette@pmmh.espci.fr Rayleigh-Bénard Convection and Lorenz Model Rayleigh-Bénard Convection Rayleigh-Bénard Convection Boussinesq Approximation Calculation and subtraction of the basic
More informationReview: Nonideal Flow in a CSTR
L3- Review: Nonideal Flow in a CSTR Ideal CSTR: uniform reactant concentration throughout the vessel Real stirred tank Relatively high reactant concentration at the feed entrance Relatively low concentration
More informationComplex system approach to geospace and climate studies. Tatjana Živković
Complex system approach to geospace and climate studies Tatjana Živković 30.11.2011 Outline of a talk Importance of complex system approach Phase space reconstruction Recurrence plot analysis Test for
More informationThe Extended Malkus-Robbins dynamo as a perturbed Lorenz system
The Extended Malkus-Robbins dynamo as a perturbed Lorenz system Irene M. Moroz Mathematical Institute 24-29 St Giles, Oxford OX1 3LB UK March 26, 2004 For correspondence: moroz@maths.ox.ac.uk 1 Abstract
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Type of reactors Dr. Zifei Liu Ideal reactors A reactor is an apparatus in which chemical, biological, and physical processes (reactions) proceed intentionally,
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationIn this report, I discuss two different topics of thermodynamics and explain how they have
Thermodynamics of Far-from-Equilibrium Systems: A Shift in Perception of Nature Truong Pham January 31st, 11 AME 36099; Directed Readings Prepared for Professor J. M. Powers 1 Introduction In this report,
More informationChemical Reaction Engineering
Lecture 2 Chemical Reaction Engineering (CRE) is the field that studies the rates and mechanisms of chemical reactions and the design of the reactors in which they take place. 1 Lecture 2 Review of Lecture
More informationMixed quadratic-cubic autocatalytic reactiondiffusion equations: Semi-analytical solutions
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 214 Mixed quadratic-cubic autocatalytic reactiondiffusion
More informationDevelopment of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
More information1. Introductory Material
CHEE 321: Chemical Reaction Engineering 1. Introductory Material 1b. The General Mole Balance Equation (GMBE) and Ideal Reactors (Fogler Chapter 1) Recap: Module 1a System with Rxn: use mole balances Input
More informationDynamic Behavior. Chapter 5
1 Dynamic Behavior In analyzing process dynamic and process control systems, it is important to know how the process responds to changes in the process inputs. A number of standard types of input changes
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: April 06, 2018, at 14 00 18 00 Johanneberg Kristian
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationDiscussion of the Lorenz Equations
Discussion of the Lorenz Equations Leibniz Universität Hannover Proseminar Theoretische Physik SS/2015 July 22, 2015 (LUH) Lorenz Equations July 22, 2015 1 / 42 Outline 1 2 3 4 5 6 7 8 (LUH) Lorenz Equations
More informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationCHEM 515: Chemical Kinetics and Dynamics
Alejandro J. Garza S01163018 Department of Chemistry, Rice University, Houston, TX email: ajg7@rice.edu, ext. 2657 Submitted December 12, 2011 Abstract Spontaneous antispiral wave formation was observed
More informationChaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB
Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the
More informationNonlinear Dynamics of Neural Firing
Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials
More informationLecture 7, 8 and 9 : Thermodynamic process by: Asst. lect. Karrar Al-Mansoori CONTENTS. 7) Thermodynamic process, path and cycle 2
CONTENTS Topics pages 7) Thermodynamic process, path and cycle 8) Reversibility and irreversibility 4 9) Thermodynamic processes and calculation of work 5 9.: Constant pressure process or isobaric process
More informationThermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space
1/29 Thermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space Jungho Park Department of Mathematics New York Institute of Technology SIAM
More information16 Period doubling route to chaos
16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More informationarxiv:chao-dyn/ v1 20 May 1994
TNT 94-4 Stochastic Resonance in Deterministic Chaotic Systems A. Crisanti, M. Falcioni arxiv:chao-dyn/9405012v1 20 May 1994 Dipartimento di Fisica, Università La Sapienza, I-00185 Roma, Italy G. Paladin
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationCHEM-UA 652: Thermodynamics and Kinetics
CHEM-UA 65: Thermodynamics and Kinetics Notes for Lecture I. THE COMPLEXITY OF MULTI-STEP CHEMICAL REACTIONS It should be clear by now that chemical kinetics is governed by the mathematics of systems of
More informationEffect of Gravity Modulation and Internal Heat Generation on Rayleigh-Bénard Convection in Couple Stress Fluid with Maxwell-Cattaneo Law
Effect of Gravity Modulation and Internal Heat Generation on Rayleigh-Bénard Convection in Couple Stress Fluid with Maxwell-Cattaneo Law Maria Thomas Research scholar, Department of Mathematics, CHRIST
More information1 Background: The Gray-Scott Reaction-Diffusion System
MSE 350 Modeling 2-D Reaction-Diffusion Instructor: R.G. Erdmann In this project, worth 3 points, you will simulate the dynamics of the Gray-Scott reaction diffusion system. (see the following: http://www-swiss.ai.mit.edu/projects/amorphous/
More informationLagrangian chaos, Eulerian chaos, and mixing enhancement in converging diverging channel flow
Lagrangian chaos, Eulerian chaos, and mixing enhancement in converging diverging channel flow Cristina H. Amon, a) Amador M. Guzmán, a) and Benoit Morel b) Carnegie Mellon University, Pittsburgh, Pennsylvania
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationTransient chaos in a closed chemical system
Transient chaos in a closed chemical system Stephen K. Scott*) School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom Bo Peng Department of Chemistry, West Virginia University, Morgantown,
More informationLecture 8. Mole balance: calculations of microreactors, membrane reactors and unsteady state in tank reactors
Lecture 8 Mole balance: calculations of microreactors, membrane reactors and unsteady state in tank reactors Mole alance in terms of Concentration and Molar Flow Rates Working in terms of number of moles
More informationSPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE. Itishree Priyadarshini. Prof. Biplab Ganguli
SPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE By Itishree Priyadarshini Under the Guidance of Prof. Biplab Ganguli Department of Physics National Institute of Technology, Rourkela CERTIFICATE This is to
More informationLesson 4: Non-fading Memory Nonlinearities
Lesson 4: Non-fading Memory Nonlinearities Nonlinear Signal Processing SS 2017 Christian Knoll Signal Processing and Speech Communication Laboratory Graz University of Technology June 22, 2017 NLSP SS
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationSolutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du
Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and
More informationA First Course on Kinetics and Reaction Engineering. Class 20 on Unit 19
A First Course on Kinetics and Reaction Engineering Class 20 on Unit 19 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going Part III - Chemical Reaction Engineering A. Ideal
More informationChaos Control for the Lorenz System
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 4, 181-188 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.8413 Chaos Control for the Lorenz System Pedro Pablo Cárdenas Alzate
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationBy Nadha CHAOS THEORY
By Nadha CHAOS THEORY What is Chaos Theory? It is a field of study within applied mathematics It studies the behavior of dynamical systems that are highly sensitive to initial conditions It deals with
More informationEect of Process Design on the Open-loop Behavior of a. Rensselaer Polytechnic Institute, Howard P. Isermann Department of Chemical Engineering, Troy,
Eect of Process Design on the Open-loop Behavior of a Jacketed Exothermic CSTR Louis P. Russo and B. Wayne Bequette Rensselaer Polytechnic nstitute, Howard P. sermann Department of Chemical Engineering,
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable
More informationIDEAL REACTORS FOR HOMOGENOUS REACTION AND THEIR PERFORMANCE EQUATIONS
IDEAL REACTORS FOR HOMOGENOUS REACTION AND THEIR PERFORMANCE EQUATIONS At the end of this week s lecture, students should be able to: Differentiate between the three ideal reactors Develop and apply the
More informationSome Dynamical Behaviors In Lorenz Model
International Journal Of Computational Engineering Research (ijceronline.com) Vol. Issue. 7 Some Dynamical Behaviors In Lorenz Model Dr. Nabajyoti Das Assistant Professor, Department of Mathematics, Jawaharlal
More information1/r plots: a brief reminder
L10-1 1/r plots: a brief reminder 1/r X target X L10-2 1/r plots: a brief reminder 1/r X target X L10-3 1/r plots: a brief reminder 1/r X target X Special Case: utocatalytic Reactions Let s assume a reaction
More informationBifurcation Analysis of Non-linear Differential Equations
Bifurcation Analysis of Non-linear Differential Equations Caitlin McCann 0064570 Supervisor: Dr. Vasiev September 01 - May 013 Contents 1 Introduction 3 Definitions 4 3 Ordinary Differential Equations
More informationOn the Asymptotic Convergence. of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans. Research School Of Chemistry
1 On the Asymptotic Convergence of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans Research School Of Chemistry Australian National University Canberra, ACT 0200 Australia
More informationReport E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache
Potsdam, August 006 Report E-Project Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory
More informationCHAPTER FIVE REACTION ENGINEERING
1 CHAPTER FIVE REACTION ENGINEERING 5.1. Determination of Kinetic Parameters of the Saponification Reaction in a PFR 5.3. Experimental and Numerical Determination of Kinetic Parameters of the Saponification
More informationLecture 8. Mole balance: calculations of microreactors, membrane reactors and unsteady state in tank reactors
Lecture 8 Mole balance: calculations of microreactors, membrane reactors and unsteady state in tank reactors Mole alance in terms of oncentration and Molar low Rates Working in terms of number of moles
More informationCHAOTIC ANALYSIS OF MERIDIONAL CIRCULATION OVER SOUTHEAST ASIA BY A FORCED MERIDIONAL-VERTICAL SYSTEM IN THE CASE OF FX = 0, FY = RF AND FZ = 0
CHAOTIC ANALYSIS OF MERIDIONAL CIRCULATION OVER SOUTHEAST ASIA BY A FORCED MERIDIONAL-VERTICAL SYSTEM IN THE CASE OF FX = 0, FY = RF AND FZ = 0 NATTAPONG WAREEPRASIRT, 2 DUSADEE SUKAWAT,2 Department of
More informationTorus Doubling Cascade in Problems with Symmetries
Proceedings of Institute of Mathematics of NAS of Ukraine 4, Vol., Part 3, 11 17 Torus Doubling Cascade in Problems with Symmetries Faridon AMDJADI School of Computing and Mathematical Sciences, Glasgow
More informationLecture 4. Mole balance: calculation of membrane reactors and unsteady state in tank reactors. Analysis of rate data
Lecture 4 Mole balance: calculation of membrane reactors and unsteady state in tank reactors. nalysis of rate data Mole alance in terms of Concentration and Molar Flow Rates Working in terms of number
More informationMathematical Modeling of the Brusselator
Mathematical Modeling of the Brusselator prepared by: Matthew P. McDowell prepared for: Joseph M. Powers, AME 36099-01 Department of Aerospace and Mechanical Engineering University of Notre Dame Notre
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. #40 Problem solving: Reactor Design Friends, this is our last session
More informationNon-Linear Models. Non-Linear Models Cont d
Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method
More informationFROM EQUILIBRIUM TO CHAOS
FROM EQUILIBRIUM TO CHAOS Practica! Bifurcation and Stability Analysis RÜDIGER SEYDEL Institut für Angewandte Mathematik und Statistik University of Würzburg Würzburg, Federal Republic of Germany ELSEVIER
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationCHEMICAL ENGINEERING KINETICS/REACTOR DESIGN. Tony Feric, Kathir Nalluswami, Manesha Ramanathan, Sejal Vispute, Varun Wadhwa
CHEMICAL ENGINEERING KINETICS/REACTOR DESIGN Tony Feric, Kathir Nalluswami, Manesha Ramanathan, Sejal Vispute, Varun Wadhwa Presentation Overview Kinetics Reactor Design Non- Isothermal Design BASICS OF
More informationPhase Plane Analysis
Phase Plane Analysis Phase plane analysis is one of the most important techniques for studying the behavior of nonlinear systems, since there is usually no analytical solution for a nonlinear system. Background
More informationSlow Manifold of a Neuronal Bursting Model
Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,
More information5. Collection and Analysis of. Rate Data
5. Collection and nalysis of o Objectives Rate Data - Determine the reaction order and specific reaction rate from experimental data obtained from either batch or flow reactors - Describe how to analyze
More informationDynamic modelling J.P. CORRIOU. Reaction and Process Engineering Laboratory University of Lorraine-CNRS, Nancy (France) Zhejiang University 2016
Dynamic modelling J.P. CORRIOU Reaction and Process Engineering Laboratory University of Lorraine-CNRS, Nancy (France) Zhejiang University 216 J.P. Corriou (LRGP) Dynamic modelling Zhejiang University
More informationChemical reactors. H has thermal contribution, pressure contribution (often negligible) and reaction contribution ( source - like)
Chemical reactors - chemical transformation of reactants into products Classification: a) according to the type of equipment o batch stirred tanks small-scale production, mostly liquids o continuous stirred
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationA First Course on Kinetics and Reaction Engineering Unit 33. Axial Dispersion Model
Unit 33. Axial Dispersion Model Overview In the plug flow reactor model, concentration only varies in the axial direction, and the sole causes of that variation are convection and reaction. Unit 33 describes
More informationChE 6303 Advanced Process Control
ChE 6303 Advanced Process Control Teacher: Dr. M. A. A. Shoukat Choudhury, Email: shoukat@buet.ac.bd Syllabus: 1. SISO control systems: Review of the concepts of process dynamics and control, process models,
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More informationHardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field
Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity
More information