Continuity equations general formulation

Besides the resuspension of particles, the perfect sink model also neglects the effect of deposited particles on incoming particles. To overcome these limitations, recent models [72, 97-99] assume that particles accumulate within a thin adsorption layer adjacent to the collector surface, and replace the perfect sink conditions with the boundary condition that particles cannot penetrate the collector. General continuity equations are formulated both for the mobile phase and for the immobilized particles in which the immobilization reaction term is decomposed in an accumulation and a removal term, respectively. Through such equations, one can keep track of the particles which arrive at the primary minimum distance and account for their normal and tangential motion. These equations were solved both approximately, and by numerical integration of the governing non-stationary transport equations. 

Rieckmann and Keil (1997) introduced a model of a 3D network of interconnected cylindrical pores with predefined distribution of pore radii and connectivity and with a volume fraction of pores equal to the porosity. The pore size distribution can be estimated from experimental characteristics obtained, e.g., from nitrogen sorption or mercury porosimetry measurements. Local heterogeneities, e.g., spatial variation in the mean pore size, or the non-uniform distribution of catalytic active centers may be taken into account in pore-network models. In each individual pore of a cylindrical or general shape, the spatially ID reaction-transport model is formulated, and the continuity equations are formulated at the nodes (i.e., connections of cylindrical capillaries) of the pore space. The transport in each individual pore is governed by the Max-well-Stefan multicomponent diffusion and convection model. Any common type of reaction kinetics taking place at the pore wall can be implemented. 

In Equation (4.12) the discretization of velocity and pressure is based on different shape functions (i.e. NjJ = l,n and Mil= l,m where, in general, mweight function used in the continuity equation is selected as -Mi to retain the symmetry of the discretized equations. After application of Green s theorem to the second-order velocity derivatives (to reduce inter-element continuity requirement) and the pressure terms (to maintain the consistency of the formulation) and algebraic manipulations the working equations of the U-V-P scheme are obtained as... 

Equation 10.1.1 represents a very general formulation of the first law of thermodynamics, which can be readily reduced to a variety of simple forms for specific applications under either steady-state or transient operating conditions. For steady-state applications the time derivative of the system energy is zero. This condition is that of greatest interest in the design of continuous flow reactors. Thus, at steady state,... 

This formulation results very insightful according to Equation 8.30, particles move under the action of an effective force — We , i.e., the nonlocal action of the quantum potential here is seen as the effect of a (nonlocal) quantum force. From a computational viewpoint, these formulation results are very interesting in connection to quantum hydrodynamics [21,27]. Thus, Equations 8.27 can be reexpressed in terms of a continuity equation and a generalized Euler equation. As happens with classical fluids, here also two important concepts that come into play the quantum pressure and the quantum vortices [28], which occur at nodal regions where the velocity field is rotational. 

As written. Equation (2) is the general formulation of a radiation balance for a single direction and a given amoimt of transported energy, but it is not very useful because we do not have sufficient information to evaluate the terms in its right-hand side. With the same approximation rigorously employed in continuous mechanics, the problem is solved with the introduction of a constitutive equation for each term, resorting to the... 

Both formulations of the constitutive equations for multicomponent diffusion, the Maxwell-Stefan equations and the generalized Fick s law, are most compactly written in matrix form. It might, therefore, be as well to begin by writing the continuity equations (Eq. 1.3.9) in n - 1 dimensional matrix form as well... 

The continuity and momentum equations that are common for these model versions are listed below. The model equations adopted for non-reactive mixtures can be deduced from the more general formulations (3.293) and (3.296), respectively. 

The second classification category distinguishes between approaches with continuous or discrete formulation for time and pipeline. In a continuous formulation, continuous coordinates are used to describe e.g. the position of a batch or a pipeline branch. The indicator u c, d indicates whether a continuous or discrete formulation is chosen, where the first component refers to the time aspect and the second refers to the pipeline. In general, time-discrete formulations are more intuitive, but to the costs of larger model dimensions in terms of the number of equations and variables. In contrast, time-continuous formulations are more compact but less intuitive. ... 

The lubrication approximation can be applied to the case when the Reynolds number is small and when the distances between the particle surfaces are much smaller than their radii of curvature (Reynolds 1886). When the role of surfactants is investigated an additional assumption is made - the Peclet number in the gap is small. Below the equations of lubrication approximation are formulated in a cylindrical coordinate system, Orz, where the droplet interface, 5, is defined as z = H(t,r)/2 and H is the local film thickness (see Fig. 1). In additional only axial symmetric flows are considered when all parameters do not depend on the meridian angle. The middle plane is z = 0 and the unit normal at the surface S pointed to the drop phase is n. The solution in the film (continuous phase) is assumed to be a mixture of nonionic and ionic surfactants and background electrolyte with relative dielectric permittivity 8f. The general formulation can be found in Kralchevsky et al. (2002). 

Stull [155] presents an excellent review of the procedures for formulating transport equations for the turbulent fluxes and variances applied to boundary layer meteorology. Wilcox [185], Pope [122] and Biswas and Eswaran [15] provide alternative texts intended for the engineering community, the latter textbook also considers experimental aspects of the field. Chou [23] was the first to derive and publish the generalized transport equation for the Reynolds stresses. The exact transport equation for the Reynolds stresses was established by use of the momentum equation, the continuity equation and a moderate amount of algebra. 

The next problem is the formulation of the transport fluxes. The precision with which it must be done for multicomponent mixtures is a matter on which there is diversity of opinion. This arises because the more precise formulations rapidly become more expensive in computer time but have uncertain returns in terms of accuracy because of shortcomings in the primary data. The author is of the opinion that a detailed transport flux model is a necessary reference representation on which approximations may be based so as to give the most satisfactory compromise for a specific problem. Fortunately, the various degrees of approximation result in formally similar continuity equations which can be treated by the same numerical techniques. Readers interested in the general aspects of modeling reactive flows are advised to omit this section on first reading and to accept the results given at the commencement of Section 4. 

Application of the weighted residual method to the solution of incompressible non-Newtonian equations of continuity and motion can be based on a variety of different schemes. Tn what follows general outlines and the formulation of the working equations of these schemes are explained. In these formulations Cauchy s equation of motion, which includes the extra stress derivatives (Equation (1.4)), is used to preseiwe the generality of the derivations. However, velocity and pressure are the only field unknowns which are obtainable from the solution of the equations of continuity and motion. The extra stress in Cauchy s equation of motion is either substituted in terms of velocity gradients or calculated via a viscoelastic constitutive equation in a separate step. 

Modelling is a process of continuous development, in which it is generally advisable to start off with the simplest conceptual representation of the process and to build in more and more complexities, as the model develops. Starting off with the process in its most complex form, often leads to confusion. A process of continuous validation is necessary, in which the model theory, data, equation formulation and model predictions must all be examined repeatedly. In formulating any model, it is therefore important to... 

As in Section II,A, a set of steady-state mass and energy balances are formulated so that the parameters that must be evaluated can be identified. The annular flow patterns are included in Regime II, and the general equations formulated in Section II,A,2,a, require a detailed knowledge of the hydrodynamics of both continuous phases and droplet interactions. Three simplified cases were formulated, and the discussion in this section is based on Case I. The steady-state mass balances are... 

In the complete Eulerian description of multiphase flows, the dispersed phase may well be conceived as a second continuous phase that interpenetrates the real continuous phase, the carrier phase this approach is often referred to as two-fluid formulation. The resulting simultaneous presence of two continua is taken into account by their respective volume fractions. All other variables such as velocities need to be averaged, in some way, in proportion to their presence various techniques have been proposed to that purpose leading, however, to different formulations of the continuum equations. The method of ensemble averaging (based on a statistical average of individual realizations) is now generally accepted as most appropriate. 

---