Nonlinear one-variable systems

Thermodynamics Far from Equilibrium Linear and Nonlinear One-Variable Systems... 

Before discussing further properties of this state function, we can proceed to nonlinear one-variable systems, which also have only one intermediate. 

The simplest nonlinear phenomenon is bistability of stationary states. It arises already in one-variable systems. In N-NDR systems this variable is the double layer potential, DL, and bistability is encountered under galvanostatic control and under potentiostatic control if the electrolyte resistance exceeds a critical value. 

The nonlinearity of the system of partial differential equations (51) and (52) poses a serious obstacle to finding an analytical solution. A reported analytical solution for the nonlinear problem of diffusion coupled with complexation kinetics was erroneous [12]. Thus, techniques such as the finite element method [53-55] or appropriate change of variables (applicable in some cases of planar diffusion) [56] should be used to find the numerical solution. One particular case of the nonlinear problem where an analytical solution can be given is the steady-state for fully labile complexes (see Section 3.3). However, there is a reasonable assumption for many relevant cases (e.g. for trace elements such as... 

The fundamental question in transport theory is Can one describe processes in nonequilibrium systems with the help of (local) thermodynamic functions of state (thermodynamic variables) This question can only be checked experimentally. On an atomic level, statistical mechanics is the appropriate theory. Since the entropy, 5, is the characteristic function for the formulation of equilibria (in a closed system), the deviation, SS, from the equilibrium value, S0, is the function which we need to use for the description of non-equilibria. Since we are interested in processes (i.e., changes in a system over time), the entropy production rate a = SS is the relevant function in irreversible thermodynamics. Irreversible processes involve linear reactions (rates 55) as well as nonlinear ones. We will be mainly concerned with processes that occur near equilibrium and so we can linearize the kinetic equations. The early development of this theory was mainly due to the Norwegian Lars Onsager. Let us regard the entropy S(a,/3,. ..) as a function of the (extensive) state variables a,/ ,. .. .which are either constant (fi,.. .) or can be controlled and measured (a). In terms of the entropy production rate, we have (9a/0f=a)... 

This set of equations is a nonlinear eigenvalue time delay differential equation. Such equations, even for one variable, often have periodic or chaotic solutions and, from the physics of the problem are also certain of having pulse-like solutions in some systems. 

The properties characteristic for electrochemical nonlinear phenomena are determined by the electrical properties of electrochemical systems, most importantly the potential drop across the electrochemical double layer at the working electrode (WE). Compared to the characteristic length scales of the patterns that develop, the extension of the double layer perpendicular to the electrode can be ignored.2 The potential drop across the double layer can therefore be lumped into one variable, DL, and the temporal evolution law of DL at every position r along the (in general two-dimensional) electrode electrolyte interface is the central equation of any electrochemical model describing pattern formation.3 It results from a local charge bal-... 

The nonlinear one-zone models can be generalized by the introduction of explicitly stochastic terms for any of the rates. The easiest, and physically most interesting, term to introduce is that of time-variable infall. In this case, Ferrini et al. have shown that there are analytic solutions possible for the simple two-component gas model which agree well with both the equilibrium behavior predicted by the linearized models and the deterministic systems. 

Whatever the approach selected, in order to optimize a nonlinear, preparative chromatographic system, one needs to consider four different factors, the objective function selected, the experimental parameters that can be optimized, the decision variable, and the constraints that must be satisfied [1]. The objective functions are discussed later, in Section 18.2.2). They include the production rate, its cost, the recovery )deld, the specific solvent consumption, or some combination of the above. Many experimental parameters of a preparative chromatography separation cannot be changed during the optimization process. Typically, such parameters are the nature of the feed components, their relative compositions, and the nature or even the brand of packing material. Sometimes, the coliunn diameter is also fixed by prior investments. The decisions variables are those parameters that can be changed during the optimization process, in order to maximize the ob-... 

First, we will study the linearization of a nonlinear equation with one variable and then we will extend it to multivariable systems. 

In previous sections we developed the linearized approximation of a nonlinear dynamic system that had only one variable. Let us now extend that approach to systems with more than one variable. 

A final comment is in order. In the previous and present sections we considered the presence of state variables only in the nonlinear functions. Thus for systems with one variable, we had only the state x, and for systems with two variables we had only states x, and x2. The formulation above should not be perceived as restrictive, but it is easily expanded to include the presence of input variables, like the manipulated variables and the disturbances. The following example demonstrates this point. 

If the Jacobian is evaluated numerically, it is not convenient to increment one variable at a time and to perform a call to the nonlinear system. 

Bioimmittance is frequency dependent. In dielectric or electrolytic models there is a choice between a step (relaxational) and sinusoidal (single-frequency) waveform excitation. As long as the step response waveform is exponential and linear conditions prevail, the information gathered is the same. At high voltage and current levels, the system is nonlinear, and models and parameters must be chosen with care. Results obtained with one variable cannot necessarily be recalculated to other forms. In some cases, one single pulse may be the best waveform because it limits heat and sample destruction. 

X he most commonly encountered mathematical models in engineering and science are in the form of differential equations. The dynamics of physical systems that have one independent variable can be modeled by ordinary differential equations, whereas systems with two, or more, independent variables require the use of partial differential equations. Several types of ordinary differential equations, and a few partial differential equations, render themselves to analytical (closed-form) solutions. These methods have been developed thoroughly in differential calculus. However, the great majority of differential equations, especially the nonlinear ones and those that involve large sets of... 

For the interaction between a nonlinear molecule and an atom, one can place the coordinate system at the centre of mass of the molecule so that the PES is a fiinction of tlie three spherical polar coordinates needed to specify the location of the atom. If the molecule is linear, V does not depend on <() and the PES is a fiinction of only two variables. In the general case of two nonlinear molecules, the interaction energy depends on the distance between the centres of mass, and five of the six Euler angles needed to specify the relative orientation of the molecular axes with respect to the global or space-fixed coordinate axes. 

Hamiltonian, but in practice one often begins with a phenomenological set of equations. The set of macrovariables are chosen to include the order parameter and all otlier slow variables to which it couples. Such slow variables are typically obtained from the consideration of the conservation laws and broken synnnetries of the system. The remaining degrees of freedom are assumed to vary on a much faster timescale and enter the phenomenological description as random themial noise. The resulting coupled nonlinear stochastic differential equations for such a chosen relevant set of macrovariables are collectively referred to as the Langevin field theory description. 

Considering the similarity between Figs. 1 and 2, the electrode potential E and the anodic dissolution current J in Fig. 2 correspond to the control parameter ft and the physical variable x in Fig. 1, respectively. Then it can be said that the equilibrium solution of J changes the value from J - 0 to J > 0 at the critical pitting potential pit. Therefore the critical pitting potential corresponds to the bifurcation point. From these points of view, corrosion should be classified as one of the nonequilibrium and nonlinear phenomena in complex systems, similar to other phenomena such as chaos. 

Let us consider the genera class of systems described by a system of n nonlinear parabolic or hyperbolic partial differential equations. For simplicity vve assume that we have only one spatial independent variable, z. 

An alternative method of solving the equations is to solve them as simultaneous equations. In that case, one can specify the design variables and the desired specifications and let the computer figure out the process parameters that will achieve those objectives. It is possible to overspecify the system or to give impossible conditions. However, the biggest drawback to this method of simulation is that large sets (tens of thousands) of nonlinear algebraic equations must be solved simultaneously. As computers become faster, this is less of an impediment, provided efficient software is available. 

At this point we can see the advantage of working with the reduced problem. Most published algorithms carry a nonlinear variable for each chemical component plus one for each mineral in the system. The number of nonlinear variables in the method presented here, on the other hand, is the number of components minus the number of minerals. Depending on the size of the problem, the savings in computing effort in evaluating Equation 4.31 can be dramatic (Fig. 4.4). 

If a more complex mathematical model is employed to represent the evaporation process, you must shift from analytic to numerical methods. The material and enthalpy balances become complicated functions of temperature (and pressure). Usually all of the system parameters are specified except for the heat transfer areas in each effect (n unknown variables) and the vapor temperatures in each effect excluding the last one (n — 1 unknown variables). The model introduces n independent equations that serve as constraints, many of which are nonlinear, plus nonlinear relations among the temperatures, concentrations, and physical properties such as the enthalpy and the heat transfer coefficient. 

The system thus obtained involves N — n + 1 variables, including t//, related by the same number of equations. Since N — n of these are nonlinear equations because of the aq term, an iteration procedure is needed. One starts from a set of q values obtained for a = 0. The equations then become linear and the Gauss elimination method may thus be used to obtain these starting q values. In a second round, these values are used in the aq term and a new set of q values are obtained by... 

The notion of chaos is interwoven with the discussion of time evolution, which we do not pursue in this volume. It is worthwhile, however, to note that it is, by now, well understood that a quantum-mechanical system with a finite Hamiltonian matrix cannot satisfy many of the purely mathematical characterizations of chaos. Equally, however, over long periods of time such systems can manifest many of the qualitative features that one associates with classically chaotic systems. It is not our intention to follow this most interesting theme. Instead we seek a more modest aim, namely, to forge a link between the elementary notions of classical nonlinear dynamics and the algebraic approach. This turns out to be possible using the action-angle variables of classical mechanics. In this section we consider only the nonlinear dynamics aspects. We complete the bridge in Chapter 7. 

---