Non-smooth feedback control for Belousov- Zhabotinskii reaction-diffusion equations: semianalytical
|
|
- Allan Lloyd
- 5 years ago
- Views:
Transcription
1 University o Wollongong Research Online Faculty o Engineering and Inormation Sciences - Papers: Part A Faculty o Engineering and Inormation Sciences 6 Non-smooth eedback control or Belousov- Zhabotinskii reaction-diusion equations: semianalytical solutions Hassan Alii The University o Dammam, hyja973@uowmail.edu.au Timothy R. Marchant University o Wollongong, tim@uow.edu.au Mark Nelson University o Wollongong, mnelson@uow.edu.au Publication Details Alii, H. Y., Marchant, T. R. & Nelson, M. I. (6). Non-smooth eedback control or Belousov-Zhabotinskii reaction-diusion equations: semi-analytical solutions. Journal o Mathematical Chemistry, 54 (8), Research Online is the open access institutional repository or the University o Wollongong. For urther inormation contact the UOW Library: research-pubs@uow.edu.au
2 Non-smooth eedback control or Belousov-Zhabotinskii reactiondiusion equations: semi-analytical solutions Abstract The Belousov-Zhabotinskii reaction is considered in one and two-dimensional reaction-diusion cells. Feedback control is examined where the eedback mechanism involves varying the concentrations in the boundary reservoir, in response to the concentrations in the centre o the cell. Semi-analytical solutions are developed, via the Galerkin method, which assumes a spatial structure or the solution, and is used to approximate the governing delay partial dierential equations by a system o delay ordinary dierential equations. The orm o eedback control considered, whilst physically realistic, is non-smooth as it has discontinuous derivatives. A stability analysis o the sets o smooth delay ordinary dierential equations, which make up the ull non-smooth system, allows a band o Hop biurcation parameter space to be obtained. It is ound that Hop biurcations or the ull non-smooth system all within this band o parameter space. In the case o eedback with no delay a precise semi-analytical estimate or the stability o the ull nonsmooth system can be obtained, which corresponds well with numerical estimates. Examples o limit cycles and the transient evolution o solutions are also considered in detail. Disciplines Engineering Science and Technology Studies Publication Details Alii, H. Y., Marchant, T. R. & Nelson, M. I. (6). Non-smooth eedback control or Belousov-Zhabotinskii reaction-diusion equations: semi-analytical solutions. Journal o Mathematical Chemistry, 54 (8), This journal article is available at Research Online:
3 Non-smooth eedback control or Belousov-Zhabotinskii reaction-diusion equations: semi-analytical solutions H.Y. Alii T. R. Marchant M. I. Nelson August 5, 6 Abstract The Belousov-Zhabotinskii reaction is considered in one and two-dimensional reactiondiusion cells. Feedback control is examined where the eedback mechanism involves varying the concentrations in the boundary reservoir, in response to the concentrations in the centre o the cell. Semi-analytical solutions are developed, via the Galerkin method, which assumes a spatial structure or the solution, and is used to approximate the governing delay partial dierential equations by a system o delay ordinary dierential equations. The orm o eedback control considered, whilst physically realistic, is non-smooth as it has discontinuous derivatives. A stability analysis o the sets o smooth delay ordinary dierential equations, which make up the ull non-smooth system, allows a band o Hop biurcation parameter space to be obtained. It is ound that Hop biurcations or the ull non-smooth system all within this band o parameter space. In the case o eedback with no delay a precise semi-analytical estimate or the stability o the ull non-smooth system can be obtained, which corresponds well with numerical estimates. Examples o limit cycles and the transient evolution o solutions are also considered in detail. Keywords mathematical modelling, reaction-diusion-delay equations, Belousov-Zhabotinskii, Hop biurcations, non-smooth eedback control Mathematics Subject Classiication () 35,37,4 Introduction Oscillatory phenomena in chemical systems have been studied, by both theoreticians and experimentalists, or many decades. The Belousov-Zhabotinskii (BZ) reaction, Bray-Liebhasky and Briggs-Rauscher systems undergo periodic concentration variations and have the added interest that these oscillations can be visualized via colour changes, see []. The BZ reaction, discovered by Belousov [] in 95, is a classical one used or understanding periodic chemical and biological relaxation oscillations. The BZ reaction has a rich history o experimental, theoretical and numerical study. The range o phenomena or which the BZ reaction, and other chemical oscillator systems, prove a useul test-bed include multi-stability, chaos, bursting, reaction-diusion patterns and waves and eedback control, see Sagues and Epstein [3] or a comprehensive review o these phenomena in the context o chemical systems. H.Y. Alii School o Mathematics and Statistics, The University o Dammam, Dammam, Eastern, 344, Saudi Arabia T. R. Marchant M. I. Nelson School o Mathematics and Applied Statistics, The University o Wollongong, Wollongong, 5, N.S.W., Australia. tim marchant@uow.edu.au
4 H.Y. Alii et al. Field et al. [4] described the chemistry o the BZ reaction and presented experimental data illustrating the sustained oscillations, rate constants or the reactions and stated the ten component reactions undamental to the system. There have been many subsequent approaches used to simpliy the original system o BZ equations. Field and Noyes [5] developed the Oregonator model, which consists o three coupled ODEs, to model the ive most important BZ reactions, while [6,7] proposed the Oregonator model ɛ du = qv uv + u( u), dt δ dv dw = qv uv + w, dt where ɛ = 4, δ = 4, q = 8 4, dt = u w, () where u, v and w are the reactant concentrations. Marchant [8] considered the Gray & Scott cubic autocatalytic model in a reaction-diusion cell. The Galerkin method was used to obtain a lower-order ODE model, as an approximation to the governing PDE system. Singularity and biurcation theory theory was then used to obtain semianalytical steady-state solutions and biurcation diagrams, together with the region o parameter space, where Hop biurcations occur. The comparison between the semi-analytical and numerical solutions o the governing PDEs was ound to be excellent. Marchant [9] extended these ideas to the Gray Scott model with Michaelis Menten decay. The Fourier Galerkin series solution method has also been used to obtain numerical solutions to a steady-state diusive BZ equation by Forbes [, ]. He ound a stripy pattern, corresponding to a standing wave and that the spatial pattern is not necessarily unique. It was also shown that small amplitude patterns are not stable but that large-amplitude patterns may be quasi-stable. Experimental studies o spatial phenomena or oscillatory chemical systems have relied on the development o new types o reactors, which allow the inlux o resh reactants without stirring the reactor contents. This has been achieved by the use o gel illed reactors (which prevents advective motion) coupled to well-stirred reactant reservoirs at the boundaries. Early experiments with gel reactors were perormed by Noszticzius et al. [] and Tam et al. [3] who reported results or the BZ reaction while later studies by Bagyan et al. [4] and Lavrova et al. [5] considered glycolytic reactions. Feedback control can be applied to chemical systems to achieve the stabilization o limit cycles and unstable steady states and also to generate chaotic behaviour. Sriram [6] studied numerical and experimental simulations o electrical eedback or the BZ reaction in a CSTR reactor. The amplitude and period o the limit cycle oscillations were increased by the eedback and it was shown that the experimental observations were accurately modelled by the Oregonator model, with a eedback term added to one o the dynamic variables. Biurcation diagrams were drawn and the eect o positive eedback on the Hop biurcation parameter region was investigated numerically. Vanag et al. [7] perormed numerical simulations o the BZ reaction. They used a three-variable model o the BZ reaction, proposed by [8], and obtained good comparisons between their simulations and experimental data. The control parameter was the inlow rate or the CSTR. They showed that there are two simple ways to modiy the transition rom chaotic behaviour, by varying the strength o the eedback and the delay response. Lipták et al. [9] considered a general open CSTR system obeying the mass action law. They proposed a class o polynomial eedback that stabilizes the system, which can then be described by a generalized Hamiltonian orm. Vanag and Epstein [] reviewed the design and control o patterns in both batch oscillators and gel reactors. Some o the varied control methods discussed include the use o initial conditions, photochemical control, periodic orcing and temperature. di Bernardo et al. [] considered biurcation theory or non-smooth piecewise continuous ODE systems. Many important applications, such as control and switching problems, impact oscillators and riction systems, are governed by such systems. They reviewed biurcation theory or steadystate solutions, which lie on discontinuity boundaries, and described the new types o instabilities which can occur in the non-smooth system. Camlibel et al. [] considered the stability o a plane piecewise smooth linear system with two dependent variables and discontinuous derivatives at the steady-state solution. They derived the conditions or the overall stability o the nonsmooth system, which relate the complex eigenvalues o the two smooth systems. Csikja et al. [3]
5 Non-smooth eedback control or BZ reaction-diusion equations 3 considered a linear piecewise aine ODE model o hysteresis. They constructed piecewise smooth limit cycle solutions and considered their stability. In this paper, we study the Oregonator model () in a reaction-diusion cell with eedback control. The eedback consists o varying the concentrations in the boundary reservoirs in response to the concentrations in the centre o the cell. A delay in the response time is also modelled. Also, to be physically realistic, the concentrations at the boundary must be positive, which requires the introduction o a modulus term in the eedback response, leading to a non-smooth system. This study has a number o aims; to illustrate that the Galerkin averaging method is useul in accurately approximating the ull reaction-diusion BZ system, to explore the eect o eedback on the BZ system, and to illustrate how approximations can be ound that predict Hop biurcations or nonsmooth systems. Moreover the Oregonator model represents a classical prototype or analyzing chemical oscillations; analytical techniques ound to be useul here can potentially be applied to the vast range o other chemical oscillators known to exist (see, or example, igure 3 o Sagues and Epstein [3]) or which the number o chemical species is very large and the corresponding sets o model equations do not enjoy the simplicity o the classical Oregonator. In Section the semi-analytical model, consisting o a set o non-smooth delay ODEs, is derived by using the Galerkin method. In Section 3 steady-state solutions are ound. In Section 4 the prediction o regions o parameters space, in which Hop biurcations occur, is considered in detail. The non-smooth system consist o eight dierent sets o smooth delay ODEs. The stability o each smooth ODE system is ound and combined. This approach gives three regions, a region in which all smooth parts o the system are stable, a region in which all are unstable, and an intermediate region where some smooth systems are stable and some are unstable. Hop biurcations or the ull non-smooth system always occur in this intermediate band o parameter space. In the case o no delay a precise prediction o the Hop biurcation parameter space, o the ull non-smooth system, is ound by considering the dominant eigenvalues and the ideas o [, ]. A good comparison between the semi-analytical model and numerical results is also obtained, or steady-state solutions, the transient evolution to the steady-state and or limit cycles. The semi-analytical model. Model equations The BZ equations in a -D reaction-diusion cell have the orm du dt = k u + dv (qv uv + u( u)), ɛ dt = k v + (w qv uv), δ () dw dt = k w + u w, u = v = w = at x, y = ±, u = u a, v = v a, w = w a, t =, (3) The system () is the BZ model proposed by [6], with Dirichlet boundary conditions (3). In - D we consider the natural simpliication o () where the concentrations in the y-direction are uniorm. The equations are in non-dimensional orm with the scaled reactant concentrations, u, v and w. It is an open system; the reactor has a permeable boundary at x, y = ±, joined to a reservoir in which the reactants have zero concentrations, see [, 3] or experimental scenarios. The parameters k,, ɛ, δ, q, deined in (), are all positive where is termed the stoichiometric actor and k the diusion coeicient. Jahnke and Winree [6] report that the range [, 4] or chemical systems. We let (u s, v s, w s ) be the steady-state concentrations at the centre o the reactor x = y =. The initial concentrations are (u a, v a, w a ). We are interested in examining the eect o eedback on the reaction-diusion cell () so consider the ollowing eedback algorithm u = H u s u(,, t τ), v = H v s v(,, t τ) (4) w = H w s w(,, t τ), at x, y = ±,
6 4 H.Y. Alii et al. where the reservoir concentrations are altered, in response to the concentrations in the cell at the centre o the cell (located at x = y = ). Many studies have been undertaken o eedback control or CSTRs, where the low rate is altered in response to the concentrations in the reactor; (4) represents an analogous orm o eedback control or a reaction-diusion cell, see [3, 4, 5]. Experimentally this can be achieved using a diusive gel coupled to a CSTR, which represents the reservoir; [3] considered a BZ reaction while [4, 5] considered a glycolytic reaction. The inert gel medium prevents convective motion but allows diusion o the chemical species. By using high low rates the concentrations in the CSTR remain close to the input values and the chemical reactions in the CSTR can be neglected. The eedback is proportional to the dierence between the transient concentrations and the steady-state values at the centre o the cell, while H is the strength o the eedback and τ is the delay response. The eedback response is zero at the steady-state so (4) does not alter the steady-state solutions o () and we investigate the eect o this eedback control on the stability o the reaction-diusion cell. Note that the concentrations must always be positive in the reservoir, hence the modulus signs on the eedback terms. The modulus terms result in a continuous eedback system but with non-smooth derivatives. A Crank Nicolson inite-dierence scheme is used to ind the numerical solutions o the governing PDE system. This implicit scheme is unconditionally stable. A ourth-order Runge Kutta scheme is used to solve the semi-analytical ODE models. Other numerical methods, such as operator splitting, exist or this class o reaction diusions equations, see or example, Ropp and Shadid [7].. The Galerkin method The semi-analytical models or () in the -D and -D geometries are ound using a Galerkin method. This method assumes a spatial structure o the concentration proiles, see [8, 8]. The Galerkin method allows the governing delay PDE to be approximated by a system o delay ODEs. In -D we use the expansion ( π ) u(x, t) = (u (t) H u s + u s u d u d ) cos x ( ) 3π + u (t) cos x + H u s + u s u d u d, ( π ) v(x, t) = (v (t) H v s + v s v d v d ) cos x ( ) 3π + v (t) cos x + H v s + v s v d v d, (5) ( π ) w(x, t) = (w (t) H w s + w s w d w d ) cos x ( ) 3π + w (t) cos x + H w s + w s w d w d. The subscript d implies a delay, that is u d = u (t τ). The trial unctions are chosen so that u + u, v + v and w + w are the concentrations at the centre o the reaction-diusion cell and so the boundary conditions at x = ± are satisied. Note that the orm o (5) is not unique; a more symmetric orm or could be chosen or quadratic expressions or the spatial proiles could be used, but the level o accuracy o the method is usually independent o the orms o the basis unctions used. The PDEs () are not satisied exactly, but the ree parameters in this expansion are obtained by evaluating averaged versions o the governing equations, weighted by the basis unctions cos ( πx) and cos ( 3 πx). Then, the ODEs du dt = π [ kπ Hπ 4H + π 4 u ( 8 ɛπ 5 u v u u + 7 ) 35 u v
7 Non-smooth eedback control or BZ reaction-diusion equations 5 ( 8 ɛπ 5 u v u u + 8 ) ] 3 u v πqv πu + HM, dv dt = π [ kπ Hπ 4H + π 4 v ( 8 δπ 3 u v u v + 7 ) 35 u v 8 5δπ u v + δ w ] δ qv + HM, ] dw π = [ kπ dt Hπ 4H + π 4 w + u w + HM 3, (6) du 3π = [ 9kπ dt 4H + 3π 4 u ( 8 ɛπ 5 u 8 9 u u v + 7 ) 35 u v + ( 8 ɛπ 9 u v 7 35 u v 44 ) ] 35 u u + πqv + πu + HM 4, dv dt = 3π [ 9kπ 4H + 3π 4 v ( 8 δπ 5 u v u v + 7 ) 35 u v + 8 9δπ u v + δ w ] δ qv + HM 5, ] dw 3π = [ 9kπ dt 4H + 3π 4 w + u w + HM 6, are obtained, where the M i are given in the Appendix. The series in (5) has been truncated ater two terms. The number o terms that are used in the truncated series represents a trade-o between the accuracy and complexity o the semi-analytical solution. It is ound that a two-term method gives suicient accuracy without excessive expression swell. A one-term solution is ound by letting each o u, v and w equal to zero. As the M i includes modulus terms, the ODEs (6) represent a non-smooth system (as there are discontinuous derivatives). Moreover, as there are three dierent modulus terms the non-smooth system (6) is composed o eight dierent smooth ODE systems. For the -D geometry, the expansion ( ) ( ) u(x, y, t) = (u (t) H u s + u s u d u d ) cos πx cos πy ( ) ( ) ( ) ( ) u (t) cos πx cos πy + u (t) cos πx cos πy + H u s + u s u d u d, ( ) ( ) v(x, y, t) = (v (t) H v s + v s v d v d ) cos πx cos πy ( ) ( ) ( ) ( ) v (t) cos πx cos πy + v (t) cos πx cos πy + H v s + v s v d v d, (7) ( ) ( ) w(x, y, t) = (w (t) H w s + w s w d w d ) cos πx cos πy ( ) ( ) ( ) ( ) w (t) cos πx cos πy + w (t) cos πx cos πy + H w s + w s w d w d. is used, which also satisies the relevant -D geometry boundary conditions in (). Symmetry implies that two o the terms have the same coeicient. Averaging using the weights cos ( πx) cos ( πy) and cos ( πx) cos ( 3 πy), gives the ollowing ODE model du π = [ kπ dt Hπ 6H + π u ( 64 ɛπ 9 u v + 8 ) 45 u v π qv
8 6 H.Y. Alii et al. ( 876 ɛπ dv dt = π Hπ 6H + π dw dt du dt 575 u v u v u u + 64 ) ] 9 u π u + HN, [ kπ v ( 64 δπ 9 u v + 8 ) 45 u v + π qv δπ u v 8 45δπ u v + δ w + HN π = [ kπ Hπ 6H + π 3π = 6H + 3π ɛπ + qv ɛ dv dt = 3π 6H + 3π dw dt [ 5kπ u ɛπ ], ] w + u w + HN 3, (8) ( u v u v ) 575 u v ) ( u u u u v u π u ] + HN 4, ( 64 [ 5kπ v δπ 45 u v u v ) 575 u v δπ u v δ qv + ] δ w + HN 5, 3π ] = [ 5kπ 6H + 3π w + u w + HN 6. where N i are given in the Appendix. 3 Steady-state solutions In this section we study the steady-state solutions o the semi-analytical model or the -D and -D geometries. In order to ind steady-state solutions, we let u(t) = u s, v(t) = v s and w(t) = w s in the ODE models, which reduces them to sets o transcendental equations. At the steady state, the eedback terms, involving H, are all zero. The steady-state solutions or the -D and -D geometries are ound by solving the transcendental equations using a root-inding routine rom the Maple sotware package. For all igures in this section the diusion coeicient k = and ɛ, δ, q are given in (). Figure shows the steady-state reactant concentrations u, v and w versus stoichiometric actor, or the -D geometry. Shown are the one and two-term semi-analytical and numerical solutions at the centre o the domain, x =. There is a unique steady-state solution or the reactant concentrations. The igures show u and w decrease as increases, beore approaching a minimum at large. However, the curve or v increases as increases. Hence, or large, u and w are near zero while v increases linearly. There is an excellent comparison between the two-term semianalytical and numerical solutions, with less than.5% error or all values o stoichiometric actors up to =. The concentration versus response curve or the non-diusive BZ system (the classical Oregonator model ()), is qualitatively similar to igure with a unique steady-state solution (see Field and Noyes [5]). Figure shows the steady-state reactant concentrations u, v and w versus x, or the -D geometry. Shown are the one-term and two-term semi-analytical and numerical solutions o the governing PDE. The stoichiometric actor = 8. The solution or the reactant concentration u has two humps, or the two-term semi-analytical and numerical solutions, while the other reactants have a single central peak. The two-term solution can model the non-central peak accurately while the one-term solution cannot. The comparison between the two-term semi-analytical and numerical solutions is excellent, while the one-term solution is reasonably accurate at the centre o the domain, x =. These behaviours are qualitatively similar to the concentration proiles in cubic auto-catalytic reactions, see [8].
9 Non-smooth eedback control or BZ reaction-diusion equations 7 u.5 v w Fig. (color online) The steady-state reactant concentrations u, v and w versus stoichiometric actor, at x =, or the -D geometry. Shown are the one-term (dashed blue line) and two-term (black solid line) semi-analytical solutions and numerical (red dotted line) solutions o the governing PDEs. 3 u x -3 v 8 w x Fig. (color online) The steady-state reactant concentrations u, v and w versus x, or the -D geometry. The parameter is = 8. Shown are the one-term (dashed blue line) and two-term (black solid line) semi-analytical solutions and numerical (red dotted line) solutions o the governing PDEs. x
10 8 H.Y. Alii et al. u.5 v. w Fig. 3 (color online) The steady-state reactant concentrations u, v and w at x = y = versus stoichiometric actor, or the -D geometry. Shown are the one-term (dashed blue line) and two-term (black solid line) semi-analytical solutions and numerical (red dotted line) solutions o the governing PDEs (). u x - v x w x Fig. 4 (color online) The steady-state reactant concentrations u, v and w versus x, or the -D geometry. The parameter space is = 4. Shown are the one-term (dashed blue line) and two-term (black solid line) semi-analytical solutions and numerical (red dotted line) solutions o the governing PDEs. x
11 Non-smooth eedback control or BZ reaction-diusion equations 9 Figure 3 shows the steady-state reactant concentrations u, v and w versus the stoichiometric actor, at the centre o the domain, x = y =, or the -D geometry. As in the -D case, the one and two-term semi-analytical and numerical solutions are shown. The curves in this case are qualitatively similar to igure but the dierence between the semi-analytical and numerical solutions is slightly larger, or large values. There is a dierence o 8% between the two-term semi-analytical and numerical solutions, up to = 5. Semi-analytical solutions or a -D geometry generally have slightly larger errors than those in a -D geometry, see [8,9,8]. Figure 4 shows the steady-state reactant concentrations u, v and w versus x, or the -D case. Shown are the one-term and two-term semi-analytical solutions and numerical solutions o the governing PDE (). The stoichiometric actor = 4 and a slice or y = are shown. The numerical solutions or v show a latter concentration proile, which is more challenging to model with a series o trial unctions. The two-term semi-analytical solution is superior to the one-term solution, but extra terms in the trial unctions would be needed to urther improve the comparison. 4 Stability analysis and Hop biurcations In this section we discuss the stability o the BZ model in order to determine the parameter regions in which Hop biurcations points occur. Standard texts in biurcation theory and dynamic systems describe the theory o Hop biurcations or delay systems, see [3, 3]. Stability theory is well understood or systems o smooth ODEs but is less well developed or the non-smooth ODE system considered here, see [] or a review o current theories. We use two approaches or the analysis o our semi-analytical non-smooth ODE system, which is comprised o eight sets o smooth ODEs. The irst approach leads to the prediction o a band o parameter space, in which Hop biurcations or the ull non-smooth system occur. The second approach is to consider the dominant eigenvalues o the system and the ideas o [, ] who constructed a hybrid stability condition or a non-smooth system. This approach leads to a precise prediction o the region in which Hop biurcations occur, or the case o eedback with no delay. 4. Theoretical considerations The -D model (6) and -D model (8) consist o ODEs or u i, v i and w i. We write the modulus terms, o the orm p, as p sgn(p) in the ODEs, and separately consider the stability o each smooth part o the system. As there are three dierent sign unctions in the equations we get eight dierent sets o smooth ODEs. The smooth ODEs are expanded in a Taylor series about the steady-state solution. Let u i = u is + ɛce µt, v i = v is + ɛge µt, w i = w is + ɛme µt, i =, ɛ, (9) and substitute (9) into the systems (6) and (8), and then linearize around the steady state. The eigenvalues o the Jacobian matrix describes the growth o small perturbations in the system. This gives the characteristic equation F (µ) = m + im = or the decay rate µ = w + iw. Hop biurcation points may occur at points where µ is purely imaginary. Here, the Hop biurcation points or the -D and -D cases are ound by solving the system o equations i =, i =,,... 6 and w = m = m =, () where the i are the steady-state versions o the -D (6) and -D (8) models. For our system, this gives eight dierent regions in which Hop biurcations occur. Combining these results gives three regions: a region in which all sets o smooth ODEs o the system are stable; a region in which all smooth parts are unstable; and an intermediate region in which some smooth parts are stable and some are unstable. Hop biurcations or the ull non-smooth ODE system occur in this intermediate region o the parameter space.
12 H.Y. Alii et al. To help resolve the exact parameter space in which Hop biurcations occur the system considered in [] provides some useul insights. They considered a linear system ẋ = { A x i c T x, A + x i c T x, () where the eigenvalues o A ±, λ = ς ± ± iω ±, (ω ± > ), are complex, the steady state solution x = and A ± R, c = R. Then the discontinuous system () is stable i S < where S = ς+ ω + + ς ω. () An insight into this hybrid stability condition can be seen by considering the solutions x = de ς±t cos(ω ± t), (3) o the two ODEs ẋ = A ± x. The solution stays in its portion o the phase plane space or time t = π ω, beore the components o x change sign. Hence, the growth or decay during this time is ± ς ± π ω, which leads to the hybrid stability condition S <. ± For our system without delay, the one-term model consists o three ODEs and the two-term model consists o six ODEs. In the parameter region where the Hop biurcation occurs the long time behaviour is dominated by a complex conjugate pair o eigenvalues. The other eigenvalues have real parts that are more negative so can be ignored as they only aect the dynamics o the system at short times. Hence the dynamics o each smooth part o our system without delay is governed by a complex conjugate pair o eigenvalues, as in the linear system (). We now examine the perturbation (3) in more detail. I the components o the initial perturbation d are all positive (termed the ppp case), then the evolution o the perturbation (3) can only switch the system between the ppp and nnn (all negative components) cases. Hence in this case the perturbation to the non-smooth system will be stable i the dominant eigenvalues or the smooth ODE systems nnn and ppp satisy S <. There are three other pairs o ODE systems (the npn and pnp, ppn and nnp, pnn and npp cases) or which similar switching can occur, so or each pair, the eigenvalues must also satisy S <. The neutrally stable points are given by S = or each pair o eigenvalues and points in the (, k) plane which satisy S = are ound using a bisection method. However, or the case with delay the dynamics is no longer dominated by a single set o complex conjugate eigenvalues and a precise analytical description o the parameter space is not available. The reason or this is that the eigenvalue equation is not a cubic equation but a transcendental equation with an ininite number o solutions. 4. Hop biurcation regions and limit cycles or the -D geometry In this section, the Hop biurcation parameter regions are obtained or -D geometry in two cases: eedback with no delay and eedback with delay. In each case, the Hop biurcation curves o the eight dierent smooth systems are ound and curves plotted which divide the parameter space into three regions: one where no Hop points occur; one where all the smooth systems predict Hop points; and one region that has mixed stability. For the case with no delay precise semi-analytical estimates o the Hop biurcation parameter region are also obtained. Figure 5 represents the region o Hop biurcations in the k plane or the -D geometry with no delay. The parameters are H =. and τ =. Figure 5(a) shows the two-term semianalytical and numerical solutions. The semi-analytical predictions are o the stability o the sets o smooth ODEs, which divides the parameter space into three regions: one in which all smooth ODEs are stable, one in which they are all unstable; and one o mixed stability. The numerical Hop biurcation points lie inside the band o mixed stability. The Hop biurcation curves corresponding to the eight sets o smooth ODE systems do not intersect each other so the curves, which separate the parameter space into the three regions are smooth. Figure 5(b) shows the numerical Hop
13 Non-smooth eedback control or BZ reaction-diusion equations k (a) 4 3 k (b) Fig. 5 (color online) The regions o the -k plane in which Hop biurcations can occur or the -D geometry. The parameter values are H =. and τ =. Shown in (a) are: the two-term semi-analytical (blue and red line) and numerical solutions (black dots), Shown in (b) are: numerical points o Hop biurcation or PDEs (black dots) and the theoretical prediction o S = (red crosses). 7 6 Region C 5 Region B k 4 3 Region A Fig. 6 (color online) The regions o the -k plane in which Hop biurcations can occur or the -D geometry. Shown are the two-term semi-analytical (blue and red line) and numerical (black dots) solutions at parameter values H =. and τ =.
14 H.Y. Alii et al.. k.. 3 (a)- Feedback no delay (τ = ) k (b)- Feedback delay (τ = ) Fig. 7 (color online) The regions o the -k plane in which Hop biurcations can occur, or dierent values o H and small values o k. Shown are two-term semi-analytical solutions in the -D geometry or our dierent values o H: H = (black solid line), H = (red dotted line), H = (dashed blue line) and H = 4 (green dotted line). biurcation points and the theoretical predictions o S =, rom the two-term model. The neutral stability predictions, or each pair o ODEs, are the same, to graphical accuracy hence only one set o crosses are plotted on the igure. It can be seen that the theoretical estimates are very accurate with an error o less than 5% or all choices o. It is also worth noting that [, 4] or real chemical systems which corresponds to the lhs o the igures. For smaller the inner and outer stability regions lie close together, except or some variations at large k. Figure 6 represents the region o Hop biurcations in the k plane or the -D geometry with delay eedback control. The parameters are τ = and H =.. Shown are the two-term semi-analytical and numerical solutions. Again the numerical Hop biurcation points occur in the theoretical band o mixed stability. In the case with delay the smooth ODE Hop biurcation curves do intersect so the composite curves shown in Figure 6 are not smooth, unlike those in igure 5 or the no delay case. The stability ormula () does not generalise to our system with delay and a precise estimate o the Hop biurcation region is not available. However or smaller values o the diusion coeicient and the stoichiometric actor, k and 4 (which corresponds to chemically realistic values), the intermediate region o mixed stability is small and hence precise estimates o the region in which Hop biurcations occur, are available. For example, on the lhs o igure 6, at k =, the transition rom region C to region A occurs or =.65 and.7 or
15 Non-smooth eedback control or BZ reaction-diusion equations 3 u.5 v.5. w. 3 4 t Fig. 8 (color online) The reactant concentrations u, v and w at x = versus t or the -D geometry. The two-term semi-analytical solution (black solid line) and numerical solution (red dotted line) are shown. The parameters are τ =, H =., u a = v a = w a =., k = and =.6. u x - v x + w x Fig. 9 (color online) The reactant concentrations u, v and w at x = versus t or the -D geometry. The two-term semi-analytical solution (black solid line) and numerical solution (red dotted line) are shown. The parameters are τ =, H =., u a = v a = w a =., k = and = 3.
16 4 H.Y. Alii et al. v w..4 u Fig. (color online) The limit cycle when = 3, H =., k = and τ =. The two-term (black solid line) semi-analytical and numerical solutions (red dotted line) are shown. the numerical, and two-term semi-analytical solutions respectively. The two predictions are close, with only a 3% dierence. Figure 7 represents the Hop biurcation region in the k plane, or the -D geometry, or small values o the diusion coeicient. Shown is the two-term semi-analytical solution or H =,, and 4. (a) is the no delay τ = case and (b) is or τ =. For small values o the diusion coeicient the stability predictions rom the eight smooth systems are the same to graphical accuracy, so there is no band o mixed stability and a precise Hop prediction is obtained. We can see that appropriately chosen values o H can stabilize or destabilize regions o parameter space. For igure 7(a), the case with no delay, we can see that as the eedback parameter H increases, the region o instability is decreasing. Figure 7(b) shows that as H increases, the region o instability increases. The eect o increasing H is stabilizing or small τ and destabilizing or large τ; the critical value o τ at which the behaviour changes is τ.5. Figures 8 and 9 show the reactant concentrations u, v and w at x =, versus t or the -D geometry. The parameters are u a = v a = w a =., k =, H =. and τ =, with =.6 (in region C o igure 6) or igure 8 and = 3 (in region A o igure 6) or igure 9. The two-term semi-analytical and numerical solutions are shown. For igure 8 =.6 and the solution evolves to a steady state, with u s.47, v s.49 and w s.5 as the time becomes large, ater some initial relaxation oscillations. The comparisons between the two-term semi-analytical and numerical solutions is excellent with only a maximum % error at the steady state. For igure 9 = 3 and a periodic solution occurs. It can be seen that the maximum concentration o v is two orders o magnitude greater than that o u and w. This example represents a challenging test case or the semi-analytical solution method but the results prove to be highly accurate. The numerical amplitudes o the limit cycle or the reactant concentrations u, v and w are.6, 85.9 and., respectively. These values are very close to the two-term semi-analytical results o.6, 85.6 and., respectively. The errors in the two-term semi-analytical values are less than.5%. Figure is a view o the 3-D phase space. The two-term semi-analytical and the numerical solutions in the -D geometry are shown. The parameters are = 3, k = and τ =. The numerical period o the limit cycle or the reactants is 3., while the two-term semi-analytical period o the limit cycle is.97, a dierence o only %. The two-term semi-analytical approximation is airly close to the numerical solution over the whole parameter space and the semi-analytical limit cycle has many quantitative similarities to the numerical solution. 4.3 Hop biurcation regions and limit cycles or the -D geometry Figure represents the region o Hop biurcations in the k plane or the -D geometry with no delay. The parameters are H =. and τ =. Figure (a) shows the two-term semi-analytical and numerical solutions. The stability o the smooth ODE systems give a band o parameter space in which Hop biurcations occur. As in the -D, there are three regions and the numerical Hop biurcations occur in the band region which has mixed stability. As or the -D case the composite
17 Non-smooth eedback control or BZ reaction-diusion equations k (a).5.5 k (b) Fig. (color online) The region o the -k plane in which Hop biurcations can occur or the -D geometry. The parameters are H =. and τ =. Shown in (a) are the two-term semi-analytical (blue and red line) and numerical solutions (black dots), Shown in (b) are the numerical points o Hop biurcation (black dots) and the theoretical prediction S = (red crosses). curves in igure (a) are smooth. Figure (b) shows numerical Hop biurcations and theoretical predictions o S = or the our pairs o ODE systems rom the two-term model (which are all the same to graphical accuracy). Also, as or the case o -D geometry the prediction is excellent, with an error o less than 9% between them at all choices o. As or the -D geometry the inner and outer stability regions lie close together or small (which corresponds to chemically realistic values), except or some variations at large k. Figure represents the region o Hop biurcations in the k plane, or the -D geometry with delay. The parameters are τ = and H =.. Shown are the two-term semi-analytical and numerical solutions. As in the -D geometry, this igure shows three dierent regions: one stable, one unstable, and one o mixed stability. It can be seen that the numerical Hop biurcation points occur in the region o mixed stability and are close to the border with the unstable region. Hence nearly all o smooth ODE systems need to be unstable or the ull non-smooth system to destabilize. As or the -D geometry case the intermediate region is small or smaller values o the diusion coeicient k. One the let hand side o the igure at k =, or the transition rom region C to region A, the Hop points are given by =.9 and 3., or the numerical and two-term semi-analytical solutions respectively, which is a 7% dierence.
18 6 H.Y. Alii et al. 4 k Region C Region B Region A Fig. (color online) The region o the -k plane in which Hop biurcations can occur or the -D geometry. Shown are the two-term semi-analytical (blue and red line) and numerical (black dotted) solutions at parameters H =. and τ =. Figure 3 represents the Hop biurcation region in the k plane or the -D geometry, or small values o the diusion coeicient. Shown is the two-term semi-analytical solution or H =,.5, and. (a) is the no delay τ = case and (b) is or τ =. As in the -D case, it can see that appropriately chosen values o H can stabilize or destabilize regions o parameter space. For igure 3(a) shows the case with eedback with no delay and we can seen that as the eedback parameter H increases, the region o instability decreases. The igure 3(b) shows the case with delay eedback at τ =. In this case as H increases the region o instability grows. As in -D, the eect o increasing H or small τ is stabilizing and destabilizing or large τ; the critical τ.6. Figures 4 and 5 show the reactant concentrations u, v and w at x = y = versus t or the -D geometry. The parameters are u a = v a = w a =. and τ = with =.7 (rom region C in igure ) or igure 4 and = 3.5 (rom region A in igure ) or igure 5. The two-term semianalytical and numerical solutions are shown. For igure 4 =.7 and the solution evolves to a steady state, with u s.47, v s.57 and w s.9 as the time becomes large. The comparison between the numerical and the two-term semi-analytical solutions shows a 8% dierence in the steady state. For igure 5 = 3.5, so periodic solutions occur. The numerical amplitudes o the limit cycle or u, v and w are.5, 34.4 and.7 respectively. These values are airly close to the two-term values o.58, 4.9 and.8 respectively. As in the -D case, the two-term semi-analytical method is again accurate, with errors o less than 8%. Figure 6 is a view o the 3-D phase space. The two-term semi-analytical and the numerical solutions or the PDEs in the -D geometry are shown. The parameter choices are = 3.5, k = and τ =. The numerical and two-term semi-analytical periods o the limit cycle o the reactant at concentrations u, v and w are.55 and.49 respectively. The errors in the two-term semi-analytical values are less than 4%. As or the -D case, the semi-analytical solution or the ODEs model is close to the numerical solution o the PDEs, in this 3-D parameter space. 5 Conclusion This paper has presented semi-analytical solutions or the BZ model in a reaction-diusion cell with eedback control or both the -D and -D geometries. The Galerkin method was used to approximate the governing delay PDEs by a system o delay ODEs. A key eature o the problem is the non-smooth nature o the eedback control and the challenges this presents or analytical investigation. For the no delay case the consideration o the dominant eigenvalues together with a hybrid stability condition allows an accurate semi-analytical prediction o the Hop biurcation region to be ound. For cases with eedback delay a band o parameter space is ound in which the
19 Non-smooth eedback control or BZ reaction-diusion equations 7. k (a)- Feedback no delay (τ = ) k. 4 6 (b)- Feedback delay (τ = ) Fig. 3 (color online) The region o the -k plane in which Hop biurcations can occur, or dierent values o H and small values o k. Shown are two-term semi-analytical solutions in the -D geometry or three dierent examples o H: H = (black solid line), H =.5 (red dotted line) and H = (dashed blue line). Hop points occur. The eect o eedback is stabilizing or small delay and destabilizing or large delay. Examples o stable and unstable limit cycles were obtained with a good comparison between semi-analytical and numerical solutions. This work illustrates the useulness o the Galerkin averaging technique, or reaction diusion equations and also contributes to the understanding o stability or non-smooth systems with multiple delay terms. Future work could involve extending the method to other classes o oscillatory chemical systems, or which the model equations are much more complex than the Oregonator model. Also the results illustrate a range o interesting behaviours that can occur in the reaction-diusion cell with boundary eedback control. Hopeully this study will motivate new experimental work using this type o eedback scenario. Acknowledgement: The authors would like to thank an anonymous reeree or their useul comments
20 8 H.Y. Alii et al..6 u.3.8 v.4. w 3 4 t Fig. 4 (color online) The reactant concentrations u, v and w at x = y = versus t or the -D geometry. The twoterm semi-analytical solution (black solid line) and numerical solution (red dotted line) are shown. The parameters are τ =, H =., u a = v a = w a =., k = and =.7. 6 u x - v x + w x Fig. 5 (color online) The reactant concentrations u, v and w at x = y = versus t or the -D geometry. The twoterm semi-analytical solution (black solid line) and numerical solution (red dotted line) are shown. The parameters are τ =, H =., u a = v a = w a =., k = and = 3.5. t
21 Non-smooth eedback control or BZ reaction-diusion equations 9 v 4.4 w.8.3 u Fig. 6 (color online) The limit cycle when = 3.5, H =., k = and τ =. The two-term (black solid line) semi-analytical and numerical solutions (red dotted line) are shown. Appendix: Expressions or the semi-analytical ODEs. This appendix presents relevant expression or or the semi-analytical models. The M i or the -D model (6) are M = H ɛ u dv d + H ɛ u d H 3ɛπ u dv d ɛ u d + 6u 3ɛπ u d + 6u 5ɛπ u d + 8u 3ɛπ v d + kπ 4 u d q ɛ v d + 8v 5ɛπ u d v ɛ u d H 3ɛπ u d u ɛ u d u ɛ v d, M = H δ u dv d + 8u 3δπ v d H 3δπ u dv d + 8v 5δπ u d u δ v d 4q δπ v d + 8v 3δπ u d + kπ 4 v d v δ u d δ w d + 8u 5δπ v d + 4 δπ w d + q δ v d, M 3 = 4 π u d u d + w d + kπ 4 w d 4 π w d, M 4 = 8v 5ɛπ u d + 8u 5ɛπ v d ɛπ u d 4 3ɛπ u d + 4H 5ɛπ u dv d + 4H 5ɛπ u d + 7u 35ɛπ v d u ɛ u d + 6u 5ɛπ u d + 7v 35ɛπ u d v ɛ u d u ɛ v d + 4q 3ɛπ v d, M 5 = 4q 3δπ v d 4 3δπ w d + 7v 35δπ u d + 7u 35δπ v d + 4 5δπ u dv d + 8u 5δπ v d + 8v 5δπ u d u δ v d v δ u d, M 6 = 4 3π w d 4 3π u d. The N i or the -D model (8) are N = 64u 9ɛπ v d + 8v 45ɛπ u d + 6q ɛπ v d + 8u 45ɛπ u d + H ɛ u d + 64v 9ɛπ u d q ɛ v d + kπ u d + 56u 45ɛπ u d + 6 ɛπ u d 8H 9ɛπ u dv d + H ɛ u dv d 8H 9ɛπ u d u ɛ u d + 8u 9ɛπ u d v ɛ u d u ɛ v d ɛ u d, N = 8H 9δπ v du d + H δ v du d + 64v 9δπ u d + kπ v d + 8v 45δπ u d + 64u 9δπ v d v δ u d 6q δπ v d + q δπ v d + 6 δπ w d δ w d + 8u 45δπ v d u δ v d, N 3 = kπ w d 6 π w d + 6 π u d + w d u d, N 4 = 6q 3ɛπ v d u ɛ u d + 988u 575ɛπ v d v ɛ u d + 76H 45ɛπ u d + 76H 45ɛπ u dv d + 64v 45ɛπ u d + 988v 575ɛπ u d + 64u 45ɛπ v d + 876u 575ɛπ u d + 8u 45ɛπ u d u ɛ v d 6 3ɛπ u d, N 5 = 6q 3δπ v d v δ u d + 64v 45δπ u d + 988u 575δπ v d + 76H 45δπ u dv d + 64u 45δπ v d 6 3δπ w d + 988v 575δπ u d u δ v d, N 6 = 6 3π w d 6 3π u d. The delay terms are deined as u d = u s + u s u (t τ) u (t τ), w d = w s + w s w (t τ) w (t τ). v d = v s + v s v (t τ) v (t τ),
22 H.Y. Alii et al. Reerences. Corbel, J.M.L., Van Lingen, J.N.J., Zevenbergen, J.F., Gijzeman, O.L.J., Meijerink, A. Strobes: pyrotechnic compositions that show a curious oscillatory combustion. Angew. Chem. Int. Ed (3). Belousov, B.P. An oscillating reaction and its mechanism. Sborn. Reerat. Radiat. Med., (Medgiz, Moscow), 45 (959) 3. Sagues, F., Epstein, I.R. Nonlinear chemical dynamics. Dalton. Trans. 7 7 (3) 4. Field, R.J., Körös, E., Noyes, R. Oscillations in chemical systems. II. thorough analysis o temporal oscillation in the Bromate-Cerium-Malonic Acid system. J. Am. Chem. Soc (97) 5. Field, R.J., Noyes, R.M. Oscillations in chemical systems. IV. limit cycle behavior in a model o a real chemical reaction. J. Chem. Phys (974) 6. Tyson, J.J. The Belousov-Zhabotinskii Reaction. Springer-Verlag, New York (976) 7. Tyson, J.J. Oscillations, Bistability and echo waves in models o the Belousov-Zhabotinskii reaction. Ann. New York Acad. Sci (979) 8. Marchant, T.R. Cubic autocatalytic reaction-diusion equations: semi-analytical solutions. Proc. R. Soc. Lond. A () 9. Marchant, T.R. Cubic autocatalysis with Michaelis-Menten kinetics: semi-analytical solutions or the reactiondiusion cell. Chem. Engng. Sci (4). Forbes, L. Stationary patterns o chemical concentration in the Belousov-Zhabotinskii reaction. Physica D (99). Forbes, L. On stability and uniqueness o stationary one-dimensional patterns in the Belousov-Zhabotinsky reaction. Physica D (99). Noszticzius, Z., Horstemke, W., McCromick, W.D., Swinney, H.L., Tam, W., Sustained chemical waves in an annular gel reactor: a chemical pinwheel. Nature (987) 3. Tam, W., Horstemke, W., Noszticzius, Z., Swinney, H.L. Sustained sprial waves in a continuously ed unstrirred chemical reactor. J. Chem. Phys (988) 4. Bagyan, S., Mair, T., Dulos, E., Boissonade, J., DeKepper, P., Muller, S. Glycolytic oscillations and waves in an open spatial reactor: Impact o eedback regulation o phosphoructokinase. Biophys. Chem (5) 5. Lavrova, A., Bagyan, S., Mair, T., Hauser, M., Schimansky-Geier, L. Modeling o glycolytic wave propagation in an open spatial reactor with inhomogeneous substrate lux. Biosystems (5) 6. Sriram, K. Eects o positive electrical eedback in the oscillating Belousov-Zhabotinsky reaction: Experiments and simulations. Chaos Soliton Fract (6) 7. Zhu, R., Qian, L. Eliminating chaos in the Belousov-Zhabotinsky reaction by no-delay eedback and delayed eedback. Theor. Chem. Acc (3) 8. Györgyi, L., Field, R.J. A three-variable model o deterministic chaos in the Belousov-Zhabotinsky reaction. Nature (99) 9. Lipták, G., Szederkényi, G., Hangos, K. M. Hamiltonian eedback design or mass action law chemical reaction networks. IFAC Papers Online (5). Vanag, V., Epstein, I.R. Design and control o patterns in reaction-diusion systems. Chaos 8 67 (8). di Bernardo, M., Budd, C., Champneys, A.R., Kowalczyk, P. Piecewise-Smooth Dynamical Systems: Theory and Applications. London, Springer (8). Camlibel, M.K., Heemels, W., Schumacher, J.M. Stability and controllability o planar linear bimodal complementarity systems. In: Proceedings o the 4nd IEEE Conerence on decision and control, Hawai, USA (3) 3. Csikja, R., Garay, B.M., Tóth, J. Chaos via two-valued interval maps in a piecewise aine model example or hysteresis. In: Proceedings o the 9th International Symposium on Mathematical Theory o Networks and Systems, Hungary () 4. Takada, H., Shimizu, Y., Miyao, M. The number o autocatalytic reactions in systems o oscillating reactions. Forma (3) 5. Vanag, V., Zhabotinsky, A., Epstein, I. Pattern ormation in the Belousov-Zhabotinsky reaction with photochemical global eedback. J. Phys. Chem. A () 6. Jahne, W., Winree, A.T. A survey o spiral wave behaviours in the Oregonator model. Int. J. Biurc. Chaos 445 (99) 7. Ropp, D.L., Shadid, J.N. Stability o operator splitting methods or systems with indeinite operators: reactiondiusion systems. J. Comput. Phys (5) 8. Alii, H.Y., Marchant, T.R., Nelson, M.I. Generalised diusive delay logistic equations: semi-analytical solutions. Dynam. Cont. Dis. Ser. B () 9. Alii, H.Y., Marchant, T.R., Nelson, M.I. Semi-analytical solutions or the - and -D diusive Nicholson s blowlies equation. IMA J. Appl. Math (4) 3. Erneux, T. Applied Delay Dierential Equations. Springer, New York (9) 3. Hale, J. Theory o Functional Dierential Equations. Springer Verlag, New York (977)
Semi-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 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 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 informationThe diffusive Lotka-Volterra predator-prey system with delay
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 5 The diffusive Lotka-Volterra predator-prey system
More information7. 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 informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More information( 1) ( 2) ( 1) nan integer, since the potential is no longer simple harmonic.
. Anharmonic Oscillators Michael Fowler Landau (para 8) considers a simple harmonic oscillator with added small potential energy terms mα + mβ. We ll simpliy slightly by dropping the term, to give an equation
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationCircular dispersive shock waves in colloidal media
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part B Faculty of Engineering and Information Sciences 6 Circular dispersive shock waves in colloidal
More informationNONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION
NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics
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 informationComments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium
Comments on Magnetohydrodynamic Unsteady Flow o A Non- Newtonian Fluid Through A Porous Medium Mostaa A.A.Mahmoud Department o Mathematics, Faculty o Science, Benha University (358), Egypt Abstract The
More information6. Well-Stirred Reactors III
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
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationSliding Mode Control and Feedback Linearization for Non-regular Systems
Sliding Mode Control and Feedback Linearization or Non-regular Systems Fu Zhang Benito R. Fernández Pieter J. Mosterman Timothy Josserand The MathWorks, Inc. Natick, MA, 01760 University o Texas at Austin,
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationIntrinsic Small-Signal Equivalent Circuit of GaAs MESFET s
Intrinsic Small-Signal Equivalent Circuit o GaAs MESFET s M KAMECHE *, M FEHAM M MELIANI, N BENAHMED, S DALI * National Centre o Space Techniques, Algeria Telecom Laboratory, University o Tlemcen, Algeria
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationFeedback linearization control of systems with singularities: a ball-beam revisit
Chapter 1 Feedback linearization control o systems with singularities: a ball-beam revisit Fu Zhang The Mathworks, Inc zhang@mathworks.com Benito Fernndez-Rodriguez Department o Mechanical Engineering,
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationScattering of Solitons of Modified KdV Equation with Self-consistent Sources
Commun. Theor. Phys. Beijing, China 49 8 pp. 89 84 c Chinese Physical Society Vol. 49, No. 4, April 5, 8 Scattering o Solitons o Modiied KdV Equation with Sel-consistent Sources ZHANG Da-Jun and WU Hua
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationActive Control and Dynamical Analysis of two Coupled Parametrically Excited Van Der Pol Oscillators
International Reereed Journal o Engineering and Science (IRJES) ISSN (Online) 39-83X, (Print) 39-8 Volume 6, Issue 7 (July 07), PP.08-0 Active Control and Dynamical Analysis o two Coupled Parametrically
More informationBOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE WITH MAGNETIC FIELD IN A NANOFLUID
Proceedings o the International Conerence on Mechanical Engineering and Reneable Energy 7 (ICMERE7) 8 December, 7, Chittagong, Bangladesh ICMERE7-PI- BOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE
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 informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationCHAPTER 5 Reactor Dynamics. Table of Contents
1 CHAPTER 5 Reactor Dynamics prepared by Eleodor Nichita, UOIT and Benjamin Rouben, 1 & 1 Consulting, Adjunct Proessor, McMaster & UOIT Summary: This chapter addresses the time-dependent behaviour o nuclear
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More information(One Dimension) Problem: for a function f(x), find x 0 such that f(x 0 ) = 0. f(x)
Solving Nonlinear Equations & Optimization One Dimension Problem: or a unction, ind 0 such that 0 = 0. 0 One Root: The Bisection Method This one s guaranteed to converge at least to a singularity, i not
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationAnalysis of Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch. Abstract
The 22 International Congress and Exposition on Noise Control Engineering Dearborn, MI, USA. August 19-21, 22 Analysis o Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch P. Wickramarachi
More informationChaos in the Belousov-Zhabotinsky Reaction
Chaos in the Belousov-Zhabotinsky Reaction David Connolly Brad Nelson December 2, 2011 Abstract In this project, we investigate two different models of the Belousov-Zhabotinsky Reaction, the Oregonator
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationADAPTIVE CHAOS CONTROL AND SYNCHRONIZATION OF HYPERCHAOTIC LIU SYSTEM
International Journal o Computer Science, Engineering and Inormation Technology (IJCSEIT), Vol.1, No., June 011 ADAPTIVE CHAOS CONTROL AND SYNCHRONIZATION OF HYPERCHAOTIC LIU SYSTEM Sundarapandian Vaidyanathan
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationRATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions
RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.
More informationCHEMICAL OSCILLATIONS IN HOMOGENEOUS SYSTEMS 1. ESSENTIAL THERMODYNAMIC AND KINETIC CONDITIONS FOR THE OCCURRENCE OF OSCILLATIONS
CHEMICAL OSCILLATIONS IN HOMOGENEOUS SYSTEMS 1. ESSENTIAL THERMODYNAMIC AND KINETIC CONDITIONS FOR THE OCCURRENCE OF OSCILLATIONS Rodica Vîlcu and Daniela Bala abstract: This manuscript reviews the understanding
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 informationOscillatory Turing Patterns in a Simple Reaction-Diffusion System
Journal of the Korean Physical Society, Vol. 50, No. 1, January 2007, pp. 234 238 Oscillatory Turing Patterns in a Simple Reaction-Diffusion System Ruey-Tarng Liu and Sy-Sang Liaw Department of Physics,
More informationReceived: 30 July 2017; Accepted: 29 September 2017; Published: 8 October 2017
mathematics Article Least-Squares Solution o Linear Dierential Equations Daniele Mortari ID Aerospace Engineering, Texas A&M University, College Station, TX 77843, USA; mortari@tamu.edu; Tel.: +1-979-845-734
More informationCHAPTER 4 Reactor Statics. Table of Contents
CHAPTER 4 Reactor Statics Prepared by Dr. Benjamin Rouben, & Consulting, Adjunct Proessor, McMaster University & University o Ontario Institute o Technology (UOIT) and Dr. Eleodor Nichita, Associate Proessor,
More informationRadiation Effects on MHD Free Convective Heat and Mass Transfer Flow Past a Vertical Porous Flat Plate with Suction
International Journal o Science, Engineering and Technology Research (IJSETR), Volume 3, Issue 5, May 4 Radiation Eects on MHD Free Convective Heat and Mass Transer Flow Past a Vertical Porous Flat Plate
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationM. Eissa * and M. Sayed Department of Engineering Mathematics, Faculty of Electronic Engineering Menouf 32952, Egypt. *
Mathematical and Computational Applications, Vol., No., pp. 5-6, 006. Association or Scientiic Research A COMPARISON BETWEEN ACTIVE AND PASSIVE VIBRATION CONTROL OF NON-LINEAR SIMPLE PENDULUM PART II:
More informationA SIMPLE MODEL FOR BIOLOGICAL AGGREGATION WITH ASYMMETRIC SENSING
A SIMPLE MODEL FOR BIOLOGICAL AGGREGATION WITH ASYMMETRIC SENSING PAUL A. MILEWSKI, XU YANG Abstract. In this paper we introduce a simple continuum model or swarming o organisms in which there is a nonlocal
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationLecture 25: Heat and The 1st Law of Thermodynamics Prof. WAN, Xin
General Physics I Lecture 5: Heat and he 1st Law o hermodynamics Pro. WAN, Xin xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Latent Heat in Phase Changes Latent Heat he latent heat o vaporization or
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationOregonator model of the Belousov-Zhabotinsky reaction Richard K. Herz,
Oregonator model of the Belousov-Zhabotinsky reaction Richard K. Herz, rherz@ucsd.edu Boris Belousov in the 1950's discovered that a mixture of malonic acid, potassium bromate, and cerium sulfate an acidic
More informationApplication of Method of Lines on the Bifurcation Analysis of a Horizontal Rijke Tube
th ASPACC July 9, 5 Beijing, China Application o Method o Lines on the Biurcation Analysis o a Horizontal Rijke Tube Xiaochuan Yang, Ali Turan, Adel Nasser, Shenghui Lei School o Mechanical, Aerospace
More informationSWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS
J. Appl. Math. & Computing Vol. 23(2007), No. 1-2, pp. 243-256 Website: http://jamc.net SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS MIHAI POPESCU Abstract. This paper deals with determining
More informationNatural convection in a vertical strip immersed in a porous medium
European Journal o Mechanics B/Fluids 22 (2003) 545 553 Natural convection in a vertical strip immersed in a porous medium L. Martínez-Suástegui a,c.treviño b,,f.méndez a a Facultad de Ingeniería, UNAM,
More informationOn Picard value problem of some difference polynomials
Arab J Math 018 7:7 37 https://doiorg/101007/s40065-017-0189-x Arabian Journal o Mathematics Zinelâabidine Latreuch Benharrat Belaïdi On Picard value problem o some dierence polynomials Received: 4 April
More informationA REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING THICKNESS
VOL., NO. 8, APRIL 6 ISSN 89-668 ARPN Journal o Engineering and Applied Sciences 6-6 Asian Research Publishing Networ (ARPN). All rights reserved. A REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING
More informationPERSISTENCE AND EXTINCTION OF SINGLE POPULATION IN A POLLUTED ENVIRONMENT
Electronic Journal o Dierential Equations, Vol 2004(2004), No 108, pp 1 5 ISSN: 1072-6691 URL: http://ejdemathtxstateedu or http://ejdemathuntedu tp ejdemathtxstateedu (login: tp) PERSISTENCE AND EXTINCTION
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More information2015 American Journal of Engineering Research (AJER)
American Journal o Engineering Research (AJER) 2015 American Journal o Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-7, pp-33-40.ajer.org Research Paper Open Access The
More informationADAPTIVE STABILIZATION AND SYNCHRONIZATION OF HYPERCHAOTIC QI SYSTEM
ADAPTIVE STABILIZATION AND SYNCHRONIZATION OF HYPERCHAOTIC QI SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR Dr. SR Technical University Avadi, Chennai-600 062,
More informationIOSR Journal of Mathematics (IOSR-JM) e-issn: , p-issn: X.Volume12,Issue 1 Ver. III (Jan.-Feb.2016)PP
IOSR Journal o Mathematics (IOSR-JM) e-issn:78-578, p-issn: 39-765X.Volume,Issue Ver. III (Jan.-Feb.6)PP 88- www.iosrjournals.org Eect o Chemical Reaction on MHD Boundary Layer Flow o Williamson Nanoluid
More informationarxiv:chao-dyn/ v1 12 Feb 1996
Spiral Waves in Chaotic Systems Andrei Goryachev and Raymond Kapral Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON M5S 1A1, Canada arxiv:chao-dyn/96014v1 12
More informationUltra Fast Calculation of Temperature Profiles of VLSI ICs in Thermal Packages Considering Parameter Variations
Ultra Fast Calculation o Temperature Proiles o VLSI ICs in Thermal Packages Considering Parameter Variations Je-Hyoung Park, Virginia Martín Hériz, Ali Shakouri, and Sung-Mo Kang Dept. o Electrical Engineering,
More informationRobust Residual Selection for Fault Detection
Robust Residual Selection or Fault Detection Hamed Khorasgani*, Daniel E Jung**, Gautam Biswas*, Erik Frisk**, and Mattias Krysander** Abstract A number o residual generation methods have been developed
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More informationPower Spectral Analysis of Elementary Cellular Automata
Power Spectral Analysis o Elementary Cellular Automata Shigeru Ninagawa Division o Inormation and Computer Science, Kanazawa Institute o Technology, 7- Ohgigaoka, Nonoichi, Ishikawa 92-850, Japan Spectral
More informationEvolution of solitary waves for a perturbed nonlinear Schrödinger equation
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 2010 Evolution of solitary waves for a perturbed nonlinear Schrödinger
More informationEffect of the Darrieus-Landau instability on turbulent flame velocity. Abstract
1 Eect o the Darrieus-Landau instability on turbulent lame velocity Maxim Zaytsev 1, and Vitaliy Bychkov 1 1 Institute o Physics, Umeå University, S-901 87 Umeå, Sweden Moscow Institute o Physics and Technology,
More informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam Abstract. A reormulation o the planetary
More informationAdditional exercises in Stationary Stochastic Processes
Mathematical Statistics, Centre or Mathematical Sciences Lund University Additional exercises 8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More informationTime-periodic forcing of Turing patterns in the Brusselator model
Time-periodic forcing of Turing patterns in the Brusselator model B. Peña and C. Pérez García Instituto de Física. Universidad de Navarra, Irunlarrea, 1. 31008-Pamplona, Spain Abstract Experiments on temporal
More informationTransient Adsorption and Desorption in Micrometer Scale Features
Journal o The Electrochemical Society, 149 8 G461-G473 2002 0013-4651/2002/149 8 /G461/13/$7.00 The Electrochemical Society, Inc. Transient Adsorption and Desorption in Micrometer Scale Features Matthias
More informationBond strength model for interfaces between nearsurface mounted (NSM) CFRP strips and concrete
University o Wollongong Research Online Faculty o Engineering and Inormation Sciences - Papers: Part A Faculty o Engineering and Inormation Sciences 2014 Bond strength model or interaces between nearsurace
More informationMath Review and Lessons in Calculus
Math Review and Lessons in Calculus Agenda Rules o Eponents Functions Inverses Limits Calculus Rules o Eponents 0 Zero Eponent Rule a * b ab Product Rule * 3 5 a / b a-b Quotient Rule 5 / 3 -a / a Negative
More informationThu June 16 Lecture Notes: Lattice Exercises I
Thu June 6 ecture Notes: attice Exercises I T. Satogata: June USPAS Accelerator Physics Most o these notes ollow the treatment in the class text, Conte and MacKay, Chapter 6 on attice Exercises. The portions
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationAlgebra II Notes Inverse Functions Unit 1.2. Inverse of a Linear Function. Math Background
Algebra II Notes Inverse Functions Unit 1. Inverse o a Linear Function Math Background Previously, you Perormed operations with linear unctions Identiied the domain and range o linear unctions In this
More informationLimit Cycles in High-Resolution Quantized Feedback Systems
Limit Cycles in High-Resolution Quantized Feedback Systems Li Hong Idris Lim School of Engineering University of Glasgow Glasgow, United Kingdom LiHonIdris.Lim@glasgow.ac.uk Ai Poh Loh Department of Electrical
More informationAn Ensemble Kalman Smoother for Nonlinear Dynamics
1852 MONTHLY WEATHER REVIEW VOLUME 128 An Ensemble Kalman Smoother or Nonlinear Dynamics GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway PETER JAN VAN LEEUWEN Institute or Marine
More informationRoot Finding and Optimization
Root Finding and Optimization Ramses van Zon SciNet, University o Toronto Scientiic Computing Lecture 11 February 11, 2014 Root Finding It is not uncommon in scientiic computing to want solve an equation
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More information3. Several Random Variables
. Several Random Variables. Two Random Variables. Conditional Probabilit--Revisited. Statistical Independence.4 Correlation between Random Variables. Densit unction o the Sum o Two Random Variables. Probabilit
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationMathematical modelling for a C60 carbon nanotube oscillator
University o Wollongong Research Online Faculty o Inormatics - Papers (Archive) Faculty o Engineering and Inormation Sciences 006 Mathematical modelling or a C60 caron nanotue oscillator Barry J. Cox University
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationEstimation and detection of a periodic signal
Estimation and detection o a periodic signal Daniel Aronsson, Erik Björnemo, Mathias Johansson Signals and Systems Group, Uppsala University, Sweden, e-mail: Daniel.Aronsson,Erik.Bjornemo,Mathias.Johansson}@Angstrom.uu.se
More informationAnderson impurity in a semiconductor
PHYSICAL REVIEW B VOLUME 54, NUMBER 12 Anderson impurity in a semiconductor 15 SEPTEMBER 1996-II Clare C. Yu and M. Guerrero * Department o Physics and Astronomy, University o Caliornia, Irvine, Caliornia
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationMATHEMATICS: PAPER I TRIAL EXAMINATION 28 AUGUST 2015
MATHEMATICS: PAPER I TRIAL EXAMINATION 8 AUGUST 015 TIME: 3 HOURS TOTAL: 150 MARKS EXAMINATION NUMBER: PLEASE READ THE FOLLOWING INSTRUCTIONS CAREFULLY 1. Write your examination number on the paper.. This
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS
ROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS I. Thirunavukkarasu 1, V. I. George 1, G. Saravana Kumar 1 and A. Ramakalyan 2 1 Department o Instrumentation
More information