Sequential differentiation method

As a second example let us consider the fed-batch bioreactor used by Ka-logerakis and Luus (1984) to illustrate sequential experimental design methods for dynamic systems. The governing differential equations are (Lim et al., 1977) ... 

The classic methods use an ODE solver in combination with an optimization algorithm and solve the problem sequentially. This solution strategy is referred to as a sequential solution and optimization approach, since for each iteration the optimization variables are set and then the differential equation constraints are integrated. Though straightforward, this approach is generally inefficient because it requires the accurate solution of the model equations at each iteration within the optimization, even when iterates are far from the final optimal solution. 

The set of partial differential equations developed for the simultaneous transfer of moisture, hear, and reactive chemicals under saturated/unsaturated soil conditions has been solved by the Galerkin finite element method. The chemical transport equations are formulated in terms of the total analytical concentration of each component species, and can be solved sequentially (Wu and Chieng, 1995). 

Since then, a considerable amount of structural and mechanistic information has been collected and yeast enolase is probably the best understood sequential enzyme to date. It is a homodimer and requires two Mg + ions per active site for catalytic activity under physiological conditions, although magnesium can be replaced with a variety of divalent metal ions in vitro. During a catalytic turnover, the metal ions bind to the active site in a kinetically ordered, sequential manner with differential binding affinities. The mode of action of yeast enolase is illustrated in Figure 26 and is unusually well understood since several solid-state structures for each intermediate identified with kinetic methods have been determined. 

One advantage in the sequential approach is that only the parameters that are used to discretize the control variable profile are considered as the decision variables. The optimization formulated by this approach is a small scale NLP that makes it attractive to apply for solving the optimal control with large dimensional systems that are modeled by a large number of differential equations. In addition, this approach can take the advantage of available IVP solvers. However, the limitation of the sequential method is a difficulty to handle a constraint on state variables (path constraint). This is because the state variables are not directly included in NLP. 

The text reviews the methodology of kinetic analysis for simple as well as complex reactions. Attention is focused on the differential and integral methods of kinetic modelling. The statistical testing of the model and the parameter estimates required by the stochastic character of experimental data is described in detail and illustrated by several practical examples. Sequential experimental design procedures for discrimination between rival models and for obtaining parameter estimates with the greatest attainable precision are developed and applied to real cases. 

Since the suggestion of the sequential QM/MM hybrid method, Canuto, Coutinho and co-authors have applied this method with success in the study of several systems and properties shift of the electronic absorption spectrum of benzene [42], pyrimidine [51] and (3-carotene [47] in several solvents shift of the ortho-betaine in water [52] shift of the electronic absorption and emission spectrum of formaldehyde in water [53] and acetone in water [54] hydrogen interaction energy of pyridine [46] and guanine-cytosine in water [55] differential solvation of phenol and phenoxy radical in different solvents [56,57] hydrated electron [58] dipole polarizability of F in water [59] tautomeric equilibrium of 2-mercaptopyridine in water [60] NMR chemical shifts in liquid water [61] electron affinity and ionization potential of liquid water [62] and liquid ammonia [35] dipole polarizability of atomic liquids [63] etc. 

PAH is an organic anion that has been used extensively for the quantitation of renal plasma flow. PAH is approximately 17% bound to plasma proteins and is eliminated extensively by active tubular secretion. Because PAH elimination is active, saturation of the transport processes have historically been anticipated, at concentrations of PAH in plasma above 10 to 20 mg/L. Recently, Dowling and associates used a sequential infusion technique and only observed concentration-dependent renal clearance of PAH at concentrations above 100 mg/dL. Furthermore, PAH is also metabolized, possibly within the kidney, to A-acetyl-PAH, and the analytical method must be able to differentiate the parent compound from the metabohte if one desires to obtain an accurate assessment of RPF. Prescott and coworkers noted that the renal clearance of PAH alone decreases at low plasma concentrations, while the clearance of the acetyl metabolite increases. Further studies are necessary to evaluate the mechanisms and significance of these findings. The extraction ratio (ER) for PAH is 70% to 90% at plasma concentrations of 10 to 20 mg/L, hence the term effective renal plasma flow (ERPF) has been used when the clearance of PAH is not corrected for the extraction ratio or if it is assumed to be 1. Normal values for ERPF are about 650 160 mL/min for men and 600 150 mL/min for women. Children will reach normahzed adult values by 3 years of age, and ERPE will begin to decline as a function of age after 30 years, reaching... 

Equations (9), (20), and (21), and the boundary conditions define a nonlinear and coupled system of partial differential equations, solved by an FVM. The equations were linearized around a guessed value. The guessed values were updated iteratively to convergence before executing the next time step. Since the electroneutrality constraint tightly couples the potential and concentration fields, the discretized sets of algebraic equations at each node point were solved simultaneously. Attempts were made to employ a sequential solver in which the electrical field was assumed for determination of the concentration of each species. In this way, the concentration fields appear decoupled and could be determined easily with a commercial, convection-diffusion solver. A robust method for converging upon the correct electrical field was, however, not found. 

---