Conservation equations numerical method

Full numerical solution of the mass, energy, and momentum conservation equations (prognostic methods). These methods attempt to solve rigorously the full problem. Unfortunately, the full simulation of the system is an extremely difficult task and often results do not agree with available observations. Data assimilation techniques that force the model to approach observed values have been developed to overcome this problem. 

The major mechanism of a vapor cloud explosion, the feedback in the interaction of combustion, flow, and turbulence, can be readily found in this mathematical model. The combustion rate, which is primarily determined by the turbulence properties, is a source term in the conservation equation for the fuel-mass fraction. The attendant energy release results in a distribution of internal energy which is described by the equation for conservation of energy. This internal energy distribution is translated into a pressure field which drives the flow field through momentum equations. The flow field acts as source term in the turbulence model, which results in a turbulent-flow structure. Finally, the turbulence properties, together with the composition, determine the rate of combustion. This completes the circle, the feedback in the process of turbulent, premixed combustion in gas explosions. The set of equations has been solved with various numerical methods e.g., SIMPLE (Patankar 1980) SOLA-ICE (Cloutman et al. 1976). 

For most numerically solved models, a control-volume approach is used. This approach is based on dividing the modeling domain into a mesh. Between mesh points, there are finite elements or boxes. Using Taylor series expansions, the governing equations are cast in finite-difference form. Next, the equations for the two half-boxes on either side of a mesh point are set equal to each other hence, mass is rigorously conserved. This approach requires that all vectors be defined at half-mesh points, all scalars at full-mesh points, and all reaction rates at quarter-mesh points. The exact details of the numerical methods can be found elsewhere (for example, see ref 273) and are not the purview of this review article. The above approach is essentially the same as that used in CFD packages (e.g.. Fluent) or discussed in Appendix C of ref 139 and is related to other numerical methods applied to fuel-cell modeling. ... 

Field models estimate the fire environment in a space by numerically solving the conservation equations (i.e., momentum, mass, energy, diffusion, species, etc.) as a result of afire. This is usually accomplished by using a finite difference, finite element, or boundary element method. Such methods are not unique to fire protection they are used in aeronautics, mechanical engineering, structural mechanics, and environmental engineering. Field models divide a space into a large number of elements and solve the conservation equations within each element. The greater the number of elements, the more detailed the solution. The results are three-dimensional in nature and are very refined when compared to a zone-type model. 

The limitations encountered when obtaining an analytical solution to the conservation equations, as in the present work, differ from those encountered applying direct computational methods. For example, the cost of numerical computations is dependent on the grid and, especially, on the number of species for which conservation equations must be solved additional reactions do not add significantly to the computational effort. With RRA techniques, further limitations arise on the number of different reaction paths that can conveniently be included in the analysis. The analysis typically follows a sequence of reactions that make up the main path of oxidation, the most important reactions, while parallel sequences are treated as perturbations to the main solution and often are sufficiently unimportant to be neglected. The first step thus identifies a skeletal mechanism of 63 elementary steps by omitting the least important steps of the detailed mechanism [44]. 

Certain numerical methods benefit from writing the convective terms in a conservative form. For example, in the species conservation, show how the continuity equation can be used to write the substantial derivative as... 

For the systems that we have considered so far, the solutions behave smoothly in time and space. Often one can simply inspect the solution and decide if the mesh is sufficiently fine to represent it accurately. Refining mesh sizes and time steps is another simple method to assure oneself that a particular discretization was sufficient. Later we will be much more concerned about numerical accuracy and stability, especially when complex chemistry is considered. For now we take a somewhat cavalier approach, with the objective being mainly to explore the general numerical approaches to solving the conservation equations describing fluid flow. For relatively simple problems we can implement usable solutions with relative ease, for example, in a spreadsheet. 

The numerical methods for solving equations like (8.2.17), (8.2.22) and (8.2.23) are discussed in Section 5.1. In practice the conservative difference schemes are widely used for solving differential equations with the accuracy of the order 0(At + Ar2) [21, 26, 27] used as well 0(Af2 4- Ar2) [25], Unlike mathematically similar equations for the A + B —> 0 reaction (Section 5.1), where the correlation functions vary monotonously in time, the... 

Detailed modelling, or numerical simulation, provides a method we can use to study complex reactive flow processes (1). Predictions about the behavior of a physical system are obtained by solving numerically the multi-fluid conservation equations for mass, momentum, and energy. Since the success of detailed modelling is coupled to one s ability to handle an abundance of theoretical and numerical detail, this field has matured in parallel with the increase in size and speed of computers and sophistication of numerical techniques. 

Computational fluid dynamics involves the analysis of fluid flow and related phenomena such as heat and/or mass transfer, mixing, and chemical reaction using numerical solution methods. Usually the domain of interest is divided into a large number of control volumes (or computational cells or elements) which have a relatively small size in comparison with the macroscopic volume of the domain of interest. For each control volume a discrete representation of the relevant conservation equations is made after which an iterative solution procedure is invoked to obtain the solution of the nonlinear equations. Due to the advent of high-speed digital computers and the availability of powerful numerical algorithms the CFD approach has become feasible. CFD can be seen as a hybrid branch of mechanics and mathematics. CFD is based on the conservation laws for mass, momentum, and (thermal) energy, which can be expressed as follows ... 

The velocity and concentration profiles are developed along the HFs by means of the mass conservation equation and the associated boundary conditions for the solute in the inner fluid. This analysis separates the effects of the operation variables, such as hydrodynamic conditions and the geometry of the system, from the mass transfer properties of the system, described by diffusion coefficients in the aqueous and organic phases and by membrane permeability. The solution of such equations usually involves numerical methods. Different applications can be found in the literature, for example, separation and concentration of phenol, Cr(VI), etc. [48-51]. 

Discussions in Chapter 2 may be referred to for explanations of the various symbols. It is straightforward to apply such conservation equations to single-phase flows. In the case of multiphase flows also, in principle, it is possible to use these equations with appropriate boundary conditions at the interface between different phases. In such cases, however, density, viscosity and all the other relevant properties will have to change abruptly at the location of the interface. These methods, which describe and track the time-dependent behavior of the interface itself, are called front tracking methods. Numerical solution of such a set of equations is extremely difficult and enormously computation intensive. The main difficulty arises from the interaction between the moving interface and the Eulerian grid employed to solve the flow field (more discussion about numerical solutions is given in Chapters 6 and 7). 

First principle mathematical models These models solve the basic conservation equations for mass and momentum in their form as partial differential equations (PDEs) along with some method of turbulence closure and appropriate initial and boundary conditions. Such models have become more common with the steady increase in computing power and sophistication of numerical algorithms. However, there are many potential problems that must be addressed. In the verification process, the PDEs being solved must adequately represent the physics of the dispersion process especially for processes such as ground-to-cloud heat transfer, phase changes for condensed phases, and chemical reactions. Also, turbulence closure methods (and associated boundary and initial conditions) must be appropriate for the dis-... 

In this chapter we will deal with steady-state and transient (or non steady-state) heat conduction in quiescent media, which occurs mostly in solid bodies. In the first section the basic differential equations for the temperature field will be derived, by combining the law of energy conservation with Fourier s law. The subsequent sections deal with steady-state and transient temperature fields with many practical applications as well as the numerical methods for solving heat conduction problems, which through the use of computers have been made easier to apply and more widespread. 

There are a variety of numerical methods available to solve the conservation equations. The most commonly used method in commercially available CFD software today is the finite volume method. Excellent descriptions of this method can be found in Refs. . With the finite volume method, an integral form of the conservation equations is solved by performing a mass and momentum balance over all faces of each computational cell. There are, however, many other methods available, such as the finite element method (where the equations are solved in differential form instead... 

Brenner (1980) has explored the subject of solute dispersion in spatially periodic porous media in considerable detail. Brenner s analysis makes use of the method of moments developed by Aris (1956) and later extended by Horn (1971). Carbonell and Whitaker (1983) and Koch et al. (1989) have addressed the same problem using the method of volume averaging, whereby mesoscopic transport coefficients are derived by averaging the basic conservation equations over a single unit cell. Numerical simulations of solute dispersion, based on lattice scale calculations of the Navier-Stokes velocity fields in spatially periodic structures, have also been performed (Eidsath et al., 1983 Edwards et al., 1991 Salles et al., 1993). These simulations are discussed in detail in the Emerging Areas section. 

---