Formulation of the Mathematical Model

2 The time constants for the temperature and concentration profiles in the pellet to change are at least an order of magnitude faster than the time constants for the temperature and concentration profiles in the reactor bed. 

We also account for density, heat capacity, and molecular weight variations due to temperature, pressure, and mole changes, along with temperature-induced variations in equilibrium constants, reaction rate constants, and heats of reaction. Axial variations of the fluid velocity arising from axial temperature changes and the change in the number of moles due to the reaction are accounted for by using the overall mass conservation or continuity equation. 

The major assumptions underlying the original model are the following. 

3 Preliminary residence time distribution studies should be conducted on the reactor to test this assumption. Although in many cases it may be desirable to increase the radial aspect ratio (possibly by crushing the catalyst), this may be difficult with highly exothermic solid-catalyzed reactions that can lead to excessive temperature excursions near the center of the bed. Carberry (1976) recommends reducing the radial aspect ratio to minimize these temperature gradients. If the velocity profile in the reactor is significantly nonuniform, the mathematical model developed here allows predictive equations such as those by Fahien and Stankovic (1979) to be easily incorporated. 

A complete description of the reactor bed involves the six differential equations that describe the catalyst, gas, and thermal well temperatures, CO and C02 concentrations, and gas velocity. These are the continuity equation, three energy balances, and two component mass balances. The following equations are written in dimensional quantities and are general for packed bed analyses. Systems without a thermal well can be treated simply by letting hts, hlg, and R0 equal zero and by eliminating the thermal well energy equation. Adiabatic conditions are simulated by setting hm and hvg equal to zero. 

---

For the sake of brevity the reader is referred to Paper II for the details regarding the constitutive mathematical models of the method applied to measure the mass flow and stoichiometry of conversion gas as well as air factors for conversion and combustion system. Below is a condensed formulation of the mathematical models applied. Here a distinction is made between measurands and sought physical quantities of the method. 

Formulation of the mathematical model here adopts the usual assumptions of equimolar overflow, constant relative volatility, total condenser, and partial reboiler. Binary variables denote the existence of trays in the column, and their sum is the number of trays N. Continuous variables represent the liquid flow rates Li and compositions xj, vapor flow rates Vi and compositions yi, the reflux Ri and vapor boilup VBi, and the column diameter Di. The equations governing the model include material and component balances around each tray, thermodynamic relations between vapor and liquid phase compositions, and the column diameter calculation based on vapor flow rate. Additional logical constraints ensure that reflux and vapor boilup enter only on one tray and that the trays are arranged sequentially (so trays cannot be skipped). Also included are the product specifications. Under the assumptions made in this example, neither the temperature nor the pressure is an explicit variable, although they could easily be included if energy balances are required. A minimum and maximum number of trays can also be imposed on the problem. 

We further discuss how quantities typically measured in the experiment (such as a rate constant) can be computed with the new formalism. The computations are based on stochastic path integral formulation [6]. Two different sources for stochasticity are considered. The first (A) is randomness that is part of the mathematical modeling and is built into the differential equations of motion (e.g. the Langevin equation, or Brownian dynamics). The second (B) is the uncertainty in the approximate numerical solution of the exact equations of motion. 

The mathematical formulation comprises of a number of mass balances and scheduling constraints. Due to the nature of the processes involved, the time aspect is prevalent in all the constraints in some form or another. A superstructure is used in the derivation of the mathematical model, as discussed in the following section. A description of the sets, variables and parameters can be found in the nomenclature list. 

In Kiil et al. (2002a), the effects of seawater parameters and paint formulation parameters on the paint behaviour are simulated and discussed. As an example of the use of the mathematical model, the effect of CU2O particle size is reviewed here. All present self-polishing antifouling paints use the soluble... 

The fundamental basis for virtually all a prion mathematical models of air pollution is the statement of conservation of mass for each pollutant species. The formulation of a mathematical model of air poUution involves a number of basic steps, the first of which is a detaUed examination of the basis of the description of the diffusion of material released into the atmosphere. The second step requires that the form of interaction among the various physical and chemical processes be specified and tested against independent experiments. Once the appropriate mathematical descriptions have been formulated, it is necessary to implement suitable solution procedures. The final step is to assess the ability of the model to predict actual ambient concentration distributions. 

Validation of the mathematical model and then the optimum formulation were then made. 

A mathematical model is expressed as an equation or a system of equations. The variables involved correspond to measurable or calculable characteristics of a real process or operation. They are functionally related to each other so that they behave much as the variables in the real system that they simulate. Some of the variables may assume, at least over a limited range, any values at the discretion of the operations planner. Others may be determined by external factors. Still others are the dependent variables, or responses, of the system. The convenient prediction of these responses for chosen values of the controllable variables is the principal reason for the formulation of a mathematical model. 

There are many different types of search routines used to locate optimum operating conditions. One approach is to make a large number of runs at different combinations of temperature, reaction time, hydrogen partial pressure, and catalyst amount, and then run a multivariable computer search routine (like the Hooke-Jeeves method or Powell method). A second approach is to formulate a mathematical model from the experimental results and then use an analytical search method to locate the optimum. The formulation of a mathematical model is not an easy task, and in many cases, this is the most critical step. Sometimes it is impossible to formulate a mathematical model for the system, as in the case of the system studied here, and an experimental search must be performed. 

The crucial step in model building is model formulation, since the mathematical modeling is intended to represent a large network of multiple biochemical reactions, controlled by complex regulatory processes that... 

These two publications indicated the importance and feasibility of mathematical approaches that describe the relations between molecular systems. They stimulated the formulation of further mathematical models and theories of chemistry that are of interest in their own right but are particularly useful for computer-assistance in chemistry. 

Some authors consider diffusion (a), (b) as consecutive processes, and assume the existence of colliding pairs [7-9]. Other models stress the importance of segmental diffusion of the active ends in a common volume of the two colliding macro molecules [10-12]. A common drawback of the mathematical models is the lack of a generally formulated expression for the effective diffusion coefficient of the active end in a coiling chain. Most models try to solve this difficulty by introducing suitable parameters with some physical meaning. 

In the various formulations of the mathematical theory of linear viscoelasticity, one should differentiate clearly the measurable and non-measurable fimctions, especially when it comes to modelling apart from the material functions quoted above, one may also define non measurable viscoelastic functions which Eu-e pure mathematical objects, such as the distribution of relaxation times, the distribution of retardation times, and tiie memory function. These mathematical tools may prove to be useful in some situations for example, a discrete distribution of relaxation times is easy to handle numerically when working with constitutive equations of the difierential type, but one has to keep in mind that the relaxation times derived numerically by optimization methods have no direct physical meaning. Furthermore, the use of the distribution of relaxation times is useless and costs precision when one wishes simply to go back and forth from the time domain to the frequency domain. This warning is important, given the large use (and sometimes overuse) of these distribution functions. 

The ideal model of chromatography, which has great importance in nonlinear chromatography, has little interest in linear chromatography. Along an infinitely efficient column, with a linear isotherm, the injection profile travels unaltered and the elution profile is the same as the injection profile. We also note here that, because of the profound difference in the formulation of the two models, the solutions of the mass balance equation of chromatography for the ideal, nonlinear model and the nonideal, linear model rely on entirely different mathematical techniques. 

The mathematical formulation of the PBPK model is dependent on several factors routes of intake of a chemical or sites of drug administration, target tissues of interest, physiological components to be explicitly modeled (kinetically important tissues and organs and the linkages between them), transport processes of the chemical (flow, diffusion, disposition, clearance, etc.), and metabolic processes involved. 

The mathematical formulation of the PBPK model is derived by applying mass balance rules across multiple compartments. The general form of mass balance equations is the same for the fat, slowly perfused, rapidly perfused, and kidney compartments, whereas the mass balance equations for the liver, gut, and lung compartments are unique. 

In summary, the microscale description provides two important pieces of information needed for the development of mesoscale models. First, the mathematical formulation of the microscale model, which includes all of the relevant physics needed to completely describe a disperse multiphase flow, provides valuable insights into what mesoscale variables are needed and how these variables interact with each other at the mesoscale. These insights are used to formulate a mesoscale model. Second, the detailed numerical solutions from the microscale model are directly used for validation of a proposed mesoscale model. When significant deviations between the mesoscale model predictions and the microscale simulations are observed, these differences lead to a reformulation of the mesoscale model in order to improve the physical description. Note that it is important to remember that this validation step should be done by comparing exact solutions to the mesoscale model with the microscale results, not approximate solutions that result... 

At the time that the basic formulation and testing of the mathematical models of quantitative structure-activity correlations were being made, another type of approach, the linear free-energy related model, was introduced (2). Using the basic Hammett equation (22, 36) for the chemical reactions of benzoic acid derivatives (Equation 12), several investigators attempted quantitative correlations between physicochemical properties... 

The chemisorbed molecules, whether on the external surface for non-porous pellets or the internal surface for porous catalyst pellets, undergo surface reaction producing chemisorbed product molecules. This surface reaction is the truly intrinsic reaction step. However, in chemical reaction engineering it is usual practice to consider that intrinsic kinetics include this surface reaction step coupled with the chemisorption steps. This is due to the difficulty of separating these steps experimentally and the ease by which they are combined mathematically in the formulation of the kinetic model. 

Thus, if the 5 formulations of the experimental design are being tested in rats, and 10 rats can be treated at a time, then the block consists of the duplicated design. At least one, and perhaps two, further blocks will be necessary. The form of the mathematical model is identical to that of the preceding case, except that a block effect, with mean value zero over all blocks is introduced in the model for the block i. The error term e includes all other variation, including the variation between dijferent animals. With two blocks there will be 12 degrees of freedom as before. This would normally be enough provided the experimental precision is adequate. However, since the repeatability includes intersubject variation, the... 

Mathematical models are very valuable because they permit the use of empirical data for calculation of other useful quantities and prediction of complex variables. Mathematical models usually explain reasons for observed behavior by giving the relationships and data used in development and validation of mathematical models. Accumulation of knowledge and data is a usual prerequisite to formulation of a mathematical model. In this sense, existence of a mathematical model usually indicates that sufficient experimental work was conducted to interpret data in a fundamental way. Below, some of these existing relationships, which help to use data on plasticizers, are discussed. 

For a given position of the interfaces 5, the solution of the mathematical model (Eqs. 15, 16, and 17) for the velocity u(xo) at a given point Xo can be obtained by means of the boundary integral formulation. The equations are non-dimensionalized and, for simplicity, the same notations for the dimensionless velocity, pressure, space, and time variables are used as those for the corresponding dimensional variables from the previous equations ... 

The statistical theory of open systems is not yet developed enough to be applied to physico-chemical problems. Both catastrophe and dissipative structure theories are of more general philosophic rather than practical value. So, only the classic Poincare-Andronov s bifurcatirMi theory gives real tools for the formulation and investigation of the mathematical models of the processes developing in physical and chemical systems far away from equilibrium. Some examples are presented in Chap. 5 where these tools were successfully applied to electrochemical systems. Main principles of such applications are given below. 

The formulation of a mathematical model becomes mode dependent. In one mode it may be an explicit state space model. If storage elements become dependent the model turns into a set of differential-algebraic equations (DAEs). 

At this stage in your career, we will show you, in some detail, two of the ten stages (marked with ) of OR application the formulation of a mathematical model and an application algorithm. In earlier chapters, you were continuously tested to formulate mathematical expressions to solve a wide variety of problems, so you are fully prepared to learn the right procedure to formulate and, with the help of an appropriate tool, solve these interesting and engaging problems. 

Chiu and co-workers [44] measured the cylindrical orthotropic thermal conductivity of spiral woven fabric composites using a mathematical model that they had devised previously. A parameter estimation technique was used to evaluate the thermal properties of spiral woven fabric composites to verify the predictability of the mathematical model. Good agreement was found between the temperatures measured in a transient heat conduction experiment and those calculated using the prediction equations formulated by the estimated parameters. 

---