Modelling equations

In this section, a model is presented for solid oxide electrolytes based upon two consecutive steps for oxygen permeation one for the surface exchange process at the oxide surface on both sides of the membrane, and another for the joint diffusion of oxygen ions and electron-holes through the solid. 

If it were supposed, for example, that the first of these reactions is at equilibrium and the second is rate-determining, the net flux of molecular oxygen through the interface at the high pressure side can be described by, 

Dou et al. [95] studied the isothermal oxygen permeation through calda-stabi- 

Dou et al. [95] performed their experiments on CSZ tubes with a homogeneous composition, having 10% pores by volume. The authors estimated that the 

---

The primary purpose for expressing experimental data through model equations is to obtain a representation that can be used confidently for systematic interpolations and extrapolations, especially to multicomponent systems. The confidence placed in the calculations depends on the confidence placed in the data and in the model. Therefore, the method of parameter estimation should also provide measures of reliability for the calculated results. This reliability depends on the uncertainties in the parameters, which, with the statistical method of data reduction used here, are estimated from the parameter variance-covariance matrix. This matrix is obtained as a last step in the iterative calculation of the parameters. 

The projection model (equations 1, 2) describes the mapping of 3D points = m,yw,2wY to 2D image points Pi = (ii,j/i) in a virtual, undistorted image plane. The... 

The distortion model (equations 3) describes the mapping of virtual, undistorted image points Pi = xi,yiY to distorted image points P = (i , y Y the real image using third order (2-dimensional) polynomals. 

The Washburn model is consistent with recent studies by Rye and co-workers of liquid flow in V-shaped grooves [49] however, the experiments are unable to distinguish between this and more sophisticated models. Equation XIII-8 is also used in studies of wicking. Wicking is the measurement of the rate of capillary rise in a porous medium to determine the average pore radius [50], surface area [51] or contact angle [52]. 

A fiirther refinement [143, 144] is to treat the atoms as being non-spherical by rewriting the repulsive part of the atom-atom exp-6 model, equation (A 1.5.71), as... 

To obtain equations describing the dusty gas model, equations (3.1) must be applied to a pseudo mixture of (n+1) species, in which the extra species, numbered n+1, represents the dust. We must also require... 

It is interesting to note that the dusty gas model equations also... 

These are equivalent to the dusty gas model equations, but are valid only for isobaric conditions, and this fact severely limits the capability of the model to represent Che behavior of systems with chemical reaction. To see this we need only remark that (8,7) and (3.8) together imply that ... 

Of course, these shortcomings of the Wakao-Smith flux relations induced by the use of equations (8.7) and (8.8) can be removed by replacing these with the corresponding dusty gas model equations, whose validity is not restricted to isobaric systems. However, since the influence of a strongly bidisperse pore size distribution can now be accounted for more simply within the class of smooth field models proposed by Feng and Stewart [49], it is hardly worthwhile pursuing this."... 

One of Che earliest examples of a properly conceived experimental investigation of the flux relations for a porous medium is provided by the work of Gunn and King [53] on the dusty gas model equations, and the following discussion is based largely on their work. Since all their experiments were performed on binary mixtures, the appropriate flux relations are (5.26) and (5,27). Writing... 

A third approach is suggested by Hugo s formulation of material balances at the limit of bulk diffusion control, described in Section 11.3. Hugo found expressions for the fluxes by combining the stoichiometric conditions and the Stefan-Maxvell relations, and this led to no inconsistencies since there are only n - 1 independent Stefan-Maxwell relations for the n fluxes. An analogous procedure can be followed when the diffusion is of intermediate type, using the dusty gas model equations in the form (5.10) and (5.11). Equations (5.11), which have the following scalar form ... 

Apart from the Knudsen limit equations (12.12), the only other reasonably compact solution of the dusty gas model equations is that given by equations (5.26) and (5.27), corresponding to binary mixtures. Consequently, if we are to study anything other than the Knudsen limit, attention will... 

The standard Galerkin technique provides a flexible and powerful method for the solution of problems in areas such as solid mechanics and heat conduction where the model equations arc of elliptic or parabolic type. It can also be used to develop robust schemes for the solution of the governing equations of... 

In the absence of body force, the dimensionless form of the governing model equations for two-dimensional steady-state incompressible creeping flow of a viscoelastic fluid are written as... 

It is evident that application of Green s theorem cannot eliminate second-order derivatives of the shape functions in the set of working equations of the least-sc[uares scheme. Therefore, direct application of these equations should, in general, be in conjunction with C continuous Hermite elements (Petera and Nassehi, 1993 Petera and Pittman, 1994). However, various techniques are available that make the use of elements in these schemes possible. For example, Bell and Surana (1994) developed a method in which the flow model equations are cast into a set of auxiliary first-order differentia] equations. They used this approach to construct a least-sciuares scheme for non-Newtonian flow equations based on equal-order C° continuous, p-version hierarchical elements. 

As described in Chapter 3, Section 5.1 the application of the VOF scheme in an Eulerian framework depends on the solution of the continuity equation for the free boundary (Equation (3.69)) with the model equations. The developed algorithm for the solution of the described model equations and updating of the free surface boundaries is as follows ... 

Here /, r and, v are unequal integers in the set 1, 2, 3. As already mentioned, in the thin-layer approach the fluid is assumed to be non-elastic and hence the stress tensor here is given in ternis of the rate of deforaiation tensor as r(p) = riD(ij), where, in the present analysis, viscosity p is defined using the power law equation. The model equations are non-dimensionalized using... 

After substitution of the leading terms of the expanded variables into the model equations and equating coefficients of equal powers of e from their sides, they are divided by common factors to obtain the following set ... 

Step 2 General structure of stiffness matrices derived for the model equations of Stokes flow in (x, 3O and (r, z) formulations (see Chapter 4) are compared. 

In this section we consider the boundary value problem for model equations of a thermoelastic plate with a vertical crack (see Khludnev, 1996d). The unknown functions in the mathematical model under consideration are such quantities as the temperature 9 and the horizontal and vertical displacements W = (w, w ), w of the mid-surface points of the plate. We use the so-called coupled model of thermoelasticity, which implies in particular that we need to solve simultaneously the equations that describe heat conduction and the deformation of the plate. The presence of the crack leads to the fact that the domain of a solution has a nonsmooth boundary. As before, the main feature of the problem as a whole is the existence of a constraint in the form of an inequality imposed on the crack faces. This constraint provides a mutual nonpenetration of the crack faces ... 

Here [ ] is the jump of a function across the crack faces and v is the normal to the surface describing the shape of the crack. Thus, we have to find a solution to the model equations of a thermoelastic plate in a domain with nonsmooth boundary and boundary conditions of the inequality type. 

Detailed Modeling Results. The results of a series of detailed calculations for an ideal isothermal plug-flow Langmuir system are summarized in Figure 15. The soHd lines show the form of the theoretical breakthrough curves for adsorption and desorption, calculated from the following set of model equations and expressed in terms of the dimensionless variables T, and P ... 

Solution of the model equations shows that, for a linear isothermal system and a pulse injection, the height equivalent to a theoretical plate (HETP) is given by... 

In general, the step count method of estimation can be apphed to any special situation to derive a model equation for that particular industry or group of processes. 

The nonequilibrium-model equations for the stage in Fig. 13-56 are as follows in residual form, where i = component (i = 1 to C), j = stage number (j = 1 to N), and V = a stage in another column that supplies an interlink. 

From the above list of rate-based model equations, it is seen that they total 5C -t- 6 for each tray, compared to 2C -t-1 or 2C -t- 3 (depending on whether mole fractious or component flow rates are used for composition variables) for each stage in the equihbrium-stage model. Therefore, more computer time is required to solve the rate-based model, which is generally converged by an SC approach of the Newton type. 

Equation (1) can be used in a general way to determine the variance resulting from the different dispersion processes that occur in an LC column. However, although the application of equation (1) to physical chemical processes may be simple, there is often a problem in identifying the average step and, sometimes, the total number of steps associated with the particular process being considered. To illustrate the use of the Random Walk model, equation (1) will be first applied to the problem of radial dispersion that occurs when a sample is placed on a packed LC column in the manner of Horne et al. [3]. 

Here, we shall examine a series of processes from the viewpoint of their kinetics and develop model reactions for the appropriate rate equations. The equations are used to andve at an expression that relates measurable parameters of the reactions to constants and to concentration terms. The rate constant or other parameters can then be determined by graphical or numerical solutions from this relationship. If the kinetics of a process are found to fit closely with the model equation that is derived, then the model can be used as a basis for the description of the process. Kinetics is concerned about the quantities of the reactants and the products and their rates of change. Since reactants disappear in reactions, their rate expressions are given a... 

It is often useful to write a model equation sueh as Equation 8-121 in terms of dimensionless variables. This introduees the Peelet number Npg = uL/Dg [, whieh represents the ratio of eharaeteristie dispersion time to eharaeteristie eonveetion time (residenee time), and the Damkdhler number,... 

The flow pattern is ealeulated from eonservation equations for mass and mometum, in eombination with the Algebraie Stress Model (ASM) for the turbulent Reynolds stresses, using the Fluent V3.03 solver. These equations ean be found in numerous textbooks and will not be reiterated here. Onee the flow pattern is known, the mixing and transport of ehemieal speeies ean be ealeulated from the following model equation ... 

---