Kinetic initial value problem

Polynomial differential equations, kinetic differential equations, kinetic initial value problems... 

Presume that a problem is described as an ordinary-differential-equation initial-value problem, such as the mass-action kinetics or plug-flow problems discussed earlier. In the standard form, such a problem might be written as... 

Simultaneous to the graph creation, kinetic properties in each vRxn are used to create the appropriate reaction rate equations (ordinary differential equations, ODE). These properties include rate constants (e.g., Michaelis constant, Km, and maximum velocity, Vmax, for enzyme-catalyzed reactions, and k for nonenzymatic reactions), inhibitor constants, A) and modes of inhibition or allosterism. The total set of rate equations and specified initial conditions forms an initial value problem that is solved by a stiff ODE equation solver for the concentrations of all species as a function of time. The constituent transforms for the each virtual enzyme are compiled by carefully culling the literature for data on enzymes known to act on the chemicals and chemical metabolites of interest. 

Note to experienced users Sometimes FEMLAB has difficulty with plug flow reactor problems, when MATLAB did not, because FEMLAB is solving them as boundary value problems, whereas MATLAB solves them as initial value problems. They are actually initial value problems you fool FEMLAB by setting the diffusion coefficient to zero, and the thermal conductivity to zero. The other difficulty is that the kinetic expression makes huge changes when the temperature changes. Anytime you have to iterate to find... 

Validation of our proposed kinetic model is illustrated by the solid and dashed curves shown by the C vs time results shown in Figs. 6-8 and 6-9 for C0 of 50 and 100 mg L-1. Here all model parameters for both the multireaction and second-order models were based on adsorption data only. With the exception of p and 0, initial conditions for this initial-value problem were the only input required. Based on these predictions, we can conclude that both models predicted Cu desorption or release behavior satisfactorily. However, predictions of desorption isotherms were not considered adequate at the initial stages of desorption following adsorption. In addition, the model underpredicted amount sorbed that directly influences subsequent predictions for the desorption isotherms. Discrepancies between experimental and predicted are expected if the amounts of Cu in the various phases (C, Se, Sh and S2) at each desorption step were significantly different. These underpredictions also may be due to the inherent assumptions of the model. Specifically, the models may... 

Considering chemical application problems, a large number of them yields mathematical models that consist of initial-value problems (IVPs) for ordinary differential equations (ODEs) or of initial-boundary-value problems (IBVPs) for partial differential equations (PDEs). Special problems of this kind, which we have treated, are diffusion-reaction processes in chemical kinetics (various polymerizations), polyreactions in microgravity environment (photoinitiated polymerization with laser beams) and drying procedures of hygroscopic porous media. 

For an arbitrary kinetic model, the surface concentrations c are solved by Equation 8.130 and the flux thus obtained, N, is inserted into the differential Equation 8.129. This equation is then solved numerically as an initial value problem (Appendix 2). A simplified solution procedure is possible for a first-order reaction. 

The kinetic system of ODEs and its initial values together provide the following initial value problem ... 

This means that the number of equations which needs to be solved is much less than the original kinetic system as discussed in Sect. 7.7.3. The calculated a values can be converted to the full concentration vector at any time point using function h. The initial value problem in Eq. (7.88) contains only N - N variables, but the values of... 

At each time tk, S bs is the observed [S], and. S pred is the predicted value from (5.53). Here, there is no analytical expression for the cost function, as we must solve the initial value problem for [S] as a function of time numerically. fit enzyme batch im1. m uses ode45 to simulate the batch kinetics for input values of the rate law parameters in order to evaluate the cost function. Either f minsearch or fminunc is used to perform the optimization. Here, we rely upon the optimizer to estimate the gradient through finite difference approximations. The agreement between the fitted equation and the data is shown in Figure 5.10. 

We can also study the kinetics of the reaction in a batch reactor. We charge the reactor with known initial concentrations ca(0) and Cb(0), and measnre the concenlrations of the species as fhnetions of time, governed by the initial value problem ... 

The three-phase isothermal reactor model was applied to analyze and simulate the performance of a bench-scale reactor. The model solution for the experimental reactor is an initial-value problem as the concentrations of reactants and products are known at the reactor inlet. The model was solved with the kinetic parameters estimated from experiments as reported previously. 

Long before electronic computers were invented, it was realized that mathematical sophistication could be introduced into numerical integration in order to save computational elTort and improve accuracy. Textbooks of numerical analysis are full of ways to do this. The most popular of them, the Runge-Kutta and predictor-corrector algorithms, once were standard methods for numerical solution of the initial value problems of chemical kinetics. They have been replaced, however, by more suitable methods invented for the specific purpose of dealing with chemical kinetics problems. 

The other state variables are the fugacity of dissolved methane in the bulk of the liquid water phase (fb) and the zero, first and second moment of the particle size distribution (p0, Pi, l )- The initial value for the fugacity, fb° is equal to the three phase equilibrium fugacity feq. The initial number of particles, p , or nuclei initially formed was calculated from a mass balance of the amount of gas consumed at the turbidity point. The explanation of the other variables and parameters as well as the initial conditions are described in detail in the reference. The equations are given to illustrate the nature of this parameter estimation problem with five ODEs, one kinetic parameter (K ) and only one measured state variable. 

For a non-premixed homogeneous flow, the initial conditions for (5.299) will usually be trivial Q(C 0 = 0. Given the chemical kinetics and the conditional scalar dissipation rate, (5.299) can thus be solved to find ((pip 0- The unconditional means (y>rp) are then found by averaging with respect to the mixture-fraction PDF. All applications reported to date have dealt with the simplest case where the mixture-fraction vector has only one component. For this case, (5.299) reduces to a simple boundary-value problem that can be easily solved using standard numerical routines. However, as discussed next, even for this simple case care must be taken in choosing the conditional scalar dissipation rate. 

This equation is similar to the Newton s equation (remember that Q is weighted by the mass) except that the component of the force in the direction perpendicular to the path is projected out and the mass is replaced by twice the kinetic energy. Prom the differential equation a finite difference formula for Q as a function of s can be obtained similarly to the initial difference formula we have for X as a function of the time t. This equation is not so popular with initial value solvers since the term E — U can go to zero, or becomes (numerically) even negative causing significant implementation problems. 

According to the relevant power and momentum balance, Eqs. (38) and (39), the electron kinetics in steady-state plasmas is characterized by tbe conditions that at any instant the power and the momentum input from the electric field are dissipated by elastic and inelastic electron collisions into the translational and internal energy of the gas particles. This instantaneous complete compensation of the respective gain from the field and the loss in collisions usually does not occur in time-dependent plasmas, and often the collisional dissipation follows with a more or less large delay—for example, the temporally varying action of a time-dependent field. Thus, the temporal response of the electrons to certain disturbances in the initial value of their velocity distribution or to rapid changes of the electric field becomes more complicated, and the study of kinetic problems related to time-dependent plasmas naturally becomes more complex and sophisticated. Despite this extended interplay between the action of the binary electron collisions and the action of the electric field, the electron kinetics in time-... 

The system of partial differential equations of first order, Eqs. (44), usually has to be treated as an initial-boundary-value problem on an appropriate energy region 0 < U < U°° and for times t > 0, where the time represents the evolution direction of the kinetic problem. Initial values for each of the distributions fo(U, i) and/,( /, t), suitable for the problem under consideration, have to be fixed, for example at t = 0. Appropriate boundary conditions for the system are given by the requirements /o([7 > U°°, t) = 0 and /,(0, t) = 0. 

As detailed below, the parabolic equation, Eq. (54), describes the evolution of the isotropic distribution and has to be solved as an initial-boundary-value problem on a nonrectangular solution region whose boundaries are partly determined by the spatial course of the electric field and thus by the specific kinetic problem considered. The parabolic problem has to be completed by appropriate initial and boundary conditions, which are briefly described below. 

Modeling of process equipment often requires a careful inspection of start-up problems that lead to differential equations of initial value type, for example, modeling of chemical kinetics in a batch reactor and modeling of a plug flow reactor. For the first example with two reactions in series... 

The popular problems of kinetics theory is the derivation of hydrodynamic equations, in certain conditions, solution of f (r, v,t) transport equation is similar the form that can relate directly to continuous or hydrodynamic description. In certain conditions the transport process is like hydrodynamic limit. In 1911 David Hilbert was who ptropwsed the existence Boltzmann equations solutions (named normal solutions), and these are determinate by initial values of hydrodynamic variables it return to collision invariant (mass, momentum and kinetics energy), Sydney Chapman and David Enskog in 1917 were whose imroUed a systematic process for derivate the hydrodynamic equations (and their corrections of superior order) for these variables. 

The existence and uniqueness of the solution of mass-action-type kinetic differential equations (or, more precisely, initial value or Cauchy problems for this type of differential equations) are ensured by general theorems, such as the Picard-Lindelof theorem (see the textbooks cited above). 

Solution of the optimization problem is simplified markedly if the initial value of guessed at once. Usually this becomes possible if one succeeds in determining the physieal, kinetic meaning of the conjugate function for the dynamic system under consideration and for the specified problem. This question will be elucidated in Section 4.3. 

The initializing of the kinetic parameter values is another problem frequently foimd in NLRA that may converge to local minimum of the objective hinction usually given by the sum of square errors between experimental and predicted product yields. 

Let us not forget, that entry conditions in kinetic problem are the initial values of concentrations of the reagents in the moment of time 1 = 0, and time-dependences of the same reagents concentrations are the desired functions. 

---