Boundary momentum method

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. 

While the film and surface-renewal theories are based on a simplified physical model of the flow situation at the interface, the boundary layer methods couple the heat and mass transfer equation directly with the momentum balance. These theories thus result in anal3dical solutions that may be considered more accurate in comparison to the film or surface-renewal models. However, to be able to solve the governing equations analytically, only very idealized flow situations can be considered. Alternatively, more realistic functional forms of the local velocity, species concentration and temperature profiles can be postulated while the functions themselves are specified under certain constraints on integral conservation. Prom these integral relationships models for the shear stress (momentum transfer), the conductive heat flux (heat transfer) and the species diffusive flux (mass transfer) can be obtained. 

Laminar Free Convection. Sparrow and Gregg [33] were the first to use the boundary layer method to study laminar, gravity-driven film condensation on a vertical plate. They improved upon the approximate analysis of Nusselt by including fluid acceleration and energy convection terms in the momentum and energy equations, respectively. Their numerical results can be expressed as ... 

This transfomi also solves the boundary value problem, i.e. there is no need to find, for an initial position x and final position a ", tlie trajectory that coimects the two points. Instead, one simply picks the initial momentum and positionp, x and calculates the classical trajectories resulting from them at all times. Such methods are generally referred to as initial variable representations (IVR). 

One simplified method for determining stack height is a geometric method described in ASHRAE. The geometric method assumes an exhaust plume shape with a lower boundary having a 1 5 slope relative to the horizontal. The stack and plume are raised until the lower plume boundary is above rooftop penthouses, separation zones, and zones of high turbulence. ASHRAE provides equations for the sizes and locations of the separation and turbulence zones. A stack height reduction credit is provided to account for the vertical exhaust momentum. 

Thus, a velocity boundary layer and a thermal boundary layer may develop simultaneously. If the physical properties of the fluid do not change significantly over the temperature range to which the fluid is subjected, the velocity boundary layer will not be affected by die heat transfer process. If physical properties are altered, there will be an interactive effect between the momentum and heat transfer processes, leading to a comparatively complex situation in which numerical methods of solution will be necessary. 

Typically, the interface obtained with the versions of the VOF method described above is smeared over a few grid cells, which, on sufficiently fine grids, allows one to identify uniquely the simply connected volumes belonging to the different phases. Instead of regarding the dynamic conditions of Eqs. (132)-(134) as boundary conditions, surface tension can be implemented as a volume force in those cells where c lies between 0 and 1. In the method developed by Brackbill et al. [176], a momentum source term of the form... 

In this chapter we will have a closer look at the methods of the reconstruction of the momentum densities and the occupation number densities for the case of CuAl alloys. An analogous reconstruction was successfully performed for LiMg alloys by Stutz etal. in 1995 [3], It was found that the shape of the Fermi surface changed and its included volume grew with Mg concentration. Finally the Fermi surface came into contact with the boundary of the first Brillouin zone in the [110] direction. Similar changes of the shape and the included volume of the Fermi surface can be expected for CuAl [4], although the higher atomic number of Cu compared to that of Li leads to problems with the reconstruction, which will be examined. 

For turbulent flow on a rotating sphere or hemisphere, Sawatzki [53] and Chin [22] have analyzed the governing equations using the Karman-Pohlhausen momentum integral method. The turbulent boundary layer was assumed to originate at the pole of rotation, and the meridional and azimuthal velocity profiles were approximated with the one-seventh power law. Their results can be summarized by the... 

CFD may be loosely thought of as computational methods applied to the study of quantities that flow. This would include both methods that solve differential equations and finite automata methods that simulate the motion of fluid particles. We shall include both of these in our discussions of the applications of CFD to packed-tube simulation in Sections III and IV. For our purposes in the present section, we consider CFD to imply the numerical solution of the Navier-Stokes momentum equations and the energy and species balances. The differential forms of these balances are solved over a large number of control volumes. These small control volumes when properly combined form the entire flow geometry. The size and number of control volumes (mesh density) are user determined and together with the chosen discretization will influence the accuracy of the solutions. After boundary conditions have been implemented, the flow and energy balances are solved numerically an iteration process decreases the error in the solution until a satisfactory result has been reached. 

If these ideas are combined with the concept of a stable distribution, i.e., one which is but minimally altered by any slight changes at the boundary, one is led to seek the S for which these effects have operated to their fullest extent, i.e., to the one which minimizes /[J ] among all those that have the same total expected mass, momentum, and energy. This of course is the method of obtaining the canonical distribution noted in Section I. 

This momentum equation is a linear parabolic partial differential equation (for constant p) that can be solved by the method of separation of variables. In this approach the solution can be found to be a product of two functions as w(t, r) = f t)g r). The solution is represented as an infinite series that can be readily evaluated at any time or value of r. Such a solution is available for a variety of boundary conditions, including time-oscillating rotation rates. At this point, however, we choose to proceed with a numerical solution. 

Following the very brief introduction to the method of lines and differential-algebraic equations, we return to solving the boundary-layer problem for nonreacting flow in a channel (Section 7.4). From the DAE-form discretization illustrated in Fig. 7.4, there are several important things to note. The residual vector F is structured as a two-dimensional matrix (e.g., Fuj represents the residual of the momentum equation at mesh point j). This organizational structure helps with the eventual software implementation. In the Fuj residual note that there are two timelike derivatives, u and p (the prime indicates the timelike z derivative). As anticipated from the earlier discussion, all the boundary conditions are handled as constraints and one is implicit. That is, the Fpj residual does not involve p itself. 

For many metals, the "nearly free" electron description corresponds quite closely 10 the physical situation. The Fermi surface remains nearly spherical in shape. However, it may now he intersected by several Brillouin zone boundaries which break the surface into a number of separate sheets. It becomes useful to describe the Fermi surface in terms not only of zones or sheets filled with electrons, but also of zones or sheets of holes, that is. momentum space volumes which are empty of electrons. A conceptually simple method of constructing these successive sheets, often also referred lo as "first zone. "second zone." and so on was demonstrated by Harrison. An example of such construction is shown in Fig. 2. 

---