Analysis of Approaches, Attempts to Average, Upscale the Turbulent Models, Governing Equations for Lower Continuum Mechanics Scale Depiction
Development of the Upper Scale Turbulent Transport Modeling Governing Equations for Porous Media
NOMENCLATURE
REFERENCES
I would like to make more assessable the structure and logic of this subsection. Getting in this direction I can modify the material from previous published version (still shortened) with the closer look at the text's parts making more comments those seems trivial for professionals in turbulent transport, still are not familiar with the VAT scaling turbulent transport concepts and mathematics.
Here I would like to remind the main unsolved difficulties in scaling description of nonlinear general and turbulent transport in fluid mechanics, thermal physics, and chemical engineering at that time  beginning of 80th. They are unsolved for many fields up to now. Those were 
Most of these issues are not understood in many professional communities up to now.
The scope of these tasks is not for the web publishing, but we are trying, with the hope that the learning capabilities of internet are still to be found.
The turbulent transport processes in highly structured or porous media are of great importance due to the large variety problems of this kind in environment, earth sciences, also for the heat and mass exchange equipment used in modern technology. These include heterogeneous media for heat exchangers and grain layers, packed columns, reactors, etc. In all cases there occurs a jet or stalled flow of fluids in channels or around the obstacles. There are, however, few theoretical developments for the flow and heat exchange in channels of complex configuration or when flowing around nonhomogeneous bodies with randomly varied parameters. The advanced forms of laminar transport equations in porous media were developed in a paper by Crapiste et al. (1986). For turbulent transport in heterogeneous media, there are few modeling approaches and their theoretical basis and final modeling equations differ.
The lack of a sound theoretical basis impacts the development of mathematical models for turbulent transport in the complex geometrical environments found in nuclear reactors subchannels where rodbundle geometries are considered to be formed by subchannels. Processes in each subchannel are calculated separately (see Teyssedou et al., 1992). The equations used in these works have often been obtained from twophase transport modeling equations (Ishii, 1975) with heterogeneity of spacial phase distributions neglected in the bulk. Threedimensional twofluid flow equationswere obtained by Ishii (1975) using a statistical averaging method. In his development, he essentially neglected nonlinear phenomena and took the flux forms of the diffusive terms to avoid averaging of the second power differential operators. Ishii and Mishima (1984) averaged a twofluid momentum equation of the form where is the local void fraction, is the mean interfacial shear stress, is the turbulent stress for the phase, is the averaged viscous stress for the phase, is the mass generation and is the generalized interfacial drag. Using the area average in the second time averaging procedure, Ishii and Mishima (1984) introduced a distribution of parameters to take into consideration the nonlinearity of convective term averaging. This approach cannot strictly take into account the stochastic character of various kinds of spatial phase distributions. The equations used by Lahey and Lopez de Bertodano (1991) and Lopez de Bertodano et al. (1990) are very similar with the momentum equation being Here the index denotes interfacial phenomena and is the volumetric wall force on phase . Additional terms in equations (1) and (2) are usually based on separate micro modeling efforts and experimental data. One of the more detailed derivations of the twophase flow governing equations by Lahey and Drew (1988) is based on a volume averaging methodology. Among the problems was that the authors developed their own volume averaging technique without consideration of theoretical advancements developed by Whitaker and colleagues (1986a)(1997) and Gray et al. (1993) for laminar and halflinear transport equations. The most important weaknesses are the lack of nonlinear terms (apart from the convection terms) which naturally arise and the nonexistance of interphase fluctuations.
Zhang and Prosperetti (1994) derived averaged equations for the motion of equal sized rigid spheres suspended in a potential flow using an equation for the probability distribution. They used the small particle dilute limit approximation to "close" the momentum equations. After approximate resolution of the continuous phase fluctuation tensor and the vector the fluctuating particle volume flux tensor, , they recognized that (p. 199)  "Closure of the system requires an expression for the fluctuating particle volume flux tensor .... This missing information cannot be supplied internally by the theory without a specification of the initial conditions imposed on the particle probability distribution". They also considered the case of "finite volume fractions for the linear problem" where the problem equations were formulated for inviscid and unconvectional media. The development by Zhang and Prosperetti (1994) is a good example of the correct application of ensemble averaging. The equations they derive compare exactly with those derived from rigorous volume averaging theory (VAT) (1999b). Transport phenomena in tube bundles of nuclear reactors and heat exchangers can be modeled by treating them as porous media Khan et al. (1975). The twodimensional momentum equations for a constant porosity distribution usually have the following form Subbotin et al. (1979), where the physical quantities are written as averaged values and the solid phase effects are included in two coefficients of bulk resistence, and , and an effective eddy viscosity, , that is not equal to the turbulent eddy viscosity. These kinds of equations were not designed to deal with nonlinearities induced by the physics of the problem and the medium variable porosity or to take into account local inhomogeneities. Some of the more interesting applications of turbulent transport in heterogeneous media are to agrometeorology, urban planning and air pollution. The first significant works on momentum and pollutant diffusion in urban environment treated as a twophase medium were the papers by Popov (1974, 1975). In these investigations, an urban porosity function was defined based on statistical averaging of a characteristic function for the surface roughness that is equal to zero inside of buildings and other structures and equal to unity in an outdoor space. The turbulent diffusion equation for an urban roughness porous medium after ensemble averaging is where means porous volume ensemble averaging, and porosity. Closure of the two "morphological" terms, the first and the second terms on the right hand side, were obtained using a Boussinesque analogy A descriptive analysis of the deviation variables , and the effective diffusion coefficient was not given. In many studies of meteorology and agronomy, the only modeling of the increase in the volume drag resistance is by addition of a nonlinear term as done by Yamada (1982), where is fraction of the earth surface occupied by forest, is the area porosity due to a tree volume and is a Coriolis parameter. The averaging technique used by Raupach and Shaw (1982) to obtain a turbulent transport equation for a twophase medium of agro and forest cultures is a plain surface 2D averaging procedure where the averaged function is defined by with being the area within the volume occupied by air. Raupach et al. (1986) and Coppin et al. (1986) assumed that the dispersive covariances were unimportant, where is a fluctuation value within the canopy and . The contribution of these covariances was found by Raupach et al. (1986) to be small in the region just above the canopy from experiments with a regular rough morphology. This finding has been explained by Scherban et al. (1986), Primak et al. (1986) and Travkin and Catton (1992,1995) for regular porous (roughness) morphology in terms of VAT. Covariances are larger, however, as the result of irregular or random twophase media. When the surface averaging used by Raupach et al. (1986) is used instead of volume averaging, especially in the case of nonisotropic media, the neglect of one of the dimensions in the averaging process results in an incorrect value. This result should be called a 2D averaging procedure, particularly when 3D averaging procedures are replaced by 2D for nonisotropic urban rough layer (URL) when developing averaged transport equations. Raupach et al. (1986) later introduced a true volume averaging procedure within an air volume that yielded the following averaged equation for momentum conservation where is interfacial area. Development of this incorrect and incomplete equation (see more explanations below in the subsection 
of this website) was based on intrinsic averaged values of or , whereas averages of vector field variables over the entire REV are more correct (see, for example, Kheifets, Neimark, 1982, etc.). Raupach et al. (1986) next simplified all the closure requirements by developing a bulk overall drag coefficient. The second, third and fifth terms on the right hand side of equation (Raupach et al., 1986) are represented by a common drag resistance term. For a stationary fully developed boundary layer, they write
where is an element drag coefficient and is an element area density  frontal area per unit volume. A wide range of flow regimes is reported in papers by Fand et al. (1987) and Dybbs and Edwards (1982). The latter work revealed that there were four regimes for regular spherical packing, and that only when the Reynolds number based on pore diameter, , exceeded 350 could the flow regime be considered to be turbulent flow. The Fand et al. (1987) investigation of a randomly packed porous medium made up of single size spheres showed that the fully developed turbulent regime occurs when > 120 where is particle Reynolds number.
Some words need to be said on the paper by Antohe and Lage (1997), which we even commented in the letter to editor of "International Journal of Heat and Mass Transfer":
".... We believe the derived equations are deficient and will attempt to explain why in the following paragraphs by analysis of the paper.
The approach taken by the authors demonstrates a misunderstanding of the basic physical principles of hierarchical modeling to obtain a phenomenological description of a transport process in porous media. The authors choose a set of phenomenological equations that are themselves the result of assumptions and simplifications and treat them as if they have a pedigree. The development of a set of equations that are rigorous does not allow one to use correlation based models developed by others that are themselves based on approximate conceptions of what the physical processes are dependent on. These models or terms in the equations already include many observed effects. After all that was their purpose. It is inadmissable for one to include such correlations in the Navier Stokes equations as done by Antohe and Lage because this results in the effects being included in the governing equations twice. Doing so incorporates the Darcy term and other factors (porosity function, for example) in the NavierStokes equations with the goal of further developing  equations.
A number of serious deficiencies found in the paper are the following: 1) The authors initial set of equations are based on the assumption that the turbulent fluctuations and fluctuations caused by the porous medium are of the same nature. They are not and serious error can result if they are assumed to be the same. 2) The initial set of equations (1)(3) is incorrect Darcy, Forchheimer and Brinkman terms are arbitrarily included in this set of equations that are assumed to be valid at some undefined scale. Porosity is arbitrarily included in the energy equation as a multiplicative factor in the convective term, again assuming that this is the correct equation for scale level heat transport description. 3) Given these deficiencies, the derived equations, ((7), (8) as well as (11), (12), (17), (19), (32), (33), (42), (45) and others) are incorrect. 4) Use and manipulation of equations obtained from different levels of approximation and different scales, using equations (9) and (18) in (17) for example, leads to an equations set that has little rigor or generality. Given the above observations, the conclusions presented in the abstract of the paper that "Among them, this conclusion supports the hypothesis of having microscopic turbulence, known to exist at high speed flow, damped by the volume averaging process. Therefore, turbulence models derived directly from the general (macroscopic) equations will inevitably fail to characterize accurately turbulence induced by the porous matrix in a microscopic sense," are not correct. Before one can reach such conclusions, the derivations of the equations upon which it is based must be valid...."
Volume averaging procedures were used by Masuoka and Takatsu (1996) to derive their volume averaged turbulent transport equations. Like numerous other studies of multiphase transport, the major difficulties of averaging the terms on the right hand were overcome by using assumed artificial closure models for the stress components terms. As a result, the averaged turbulent momentum equation, for example, has conventional additional resistance terms like the averaged momentum equation developed by Vafai and Tien (1981) for laminar regime transport in porous medium. A major assumption is the linearity of the fluctuation terms obtained, for example, by neglect of additional terms in the momentum equation.
I would like to discuss here in more detail the main arguments following the
papers by Masuoka and Takatsu (1996) and Takatsu and Masuoka (1997), including
as much elements as it was done in those years.
The momentum equation in Masuoka, and Takatsu (1996) work is
where
where
is the turbulent kinetic energy. The averaging procedure for this vectorial
equation provided using the intrinsic phase average
Drawbacks of this kind of averaging for a vector variable equation was
addressed by Kheifets and Neimark (1982).
In further development the closure of fluctuations stress terms of the right
hand side momentum equation tems for both the momentum and onetemperature
energy equation were announced neglegable. The closure of right hand side very
complicated expressions was obtained like
and
"where
is the correction factor which is introduced to extend the concept of the
hydrodynamic conductance defined by Darcy's law to the turbulent flow." ?!
So, this is the clear intention to stay within the area of laminar phenomena
parameters.
Also, the biggest assumption in this work is that the linearity of equation's
terms related to fluctuations induced by porosity is presented as
Further, after that kind of closure the averaged momentum equation became
In their experimental work Takatsu and Masuoka, "Turbulent phenomena in flow through porous media", Journal of Porous Media, Vol. 1, No. 3, pp.243251, (1997), studied by visualization the flow in the flat vertical channel with very narrow width of the channel 1 cm and circular cylinders crossing the channel in staggered regular pattern as a bank of short tubes. In fact the experiments were set up in a 3D environment. They write in the abstract: "The present experiments support the fact that the Forchheimer flow resistance and dispersion have a close relevance to the turbulent mixing". To come to this conclusion they used the modeling equations, which are derived by voluntarilly methods and are closed to be as the Forchheimer pressure resistance model. So, while reading this conclusion one needs understand that it is obtained due to the fact that the model used in this work contains the Forchheimer type resisitance term which in turn points out to this conclusion. They used the following correlations for experimental data reduction  with the hydraulic diameter also they used the flow resistance factor which is and where a) the resistance factor byKnudsen and Katz (1958) for laminar regime across a bank of tubes this correlation after substitution of the and through the and gives and for turbulent flow across a bank of tubes b) and by Bird et al. (1960) for packed beds or which is very close (actually almost equal) to the Ergun's correlation obtained for in Travkin and Catton (1995)
where the production term should be like this tensorial relation They did accounting for the additional term in "averaged" kinetic energy equation then they did the convenient thing declaring that according to their definition this term should be as which is unexplainable.
Conclusion can be made that the procedures used in this work were constructed not on the consistent scientific approach ( justifiable or not )  but using heuristic judgements to bring the final form of governing equations close to conventional one phase, one scale fluid mechanics equations.
The paper by Wang and Takle (1995) is devoted to the same noble aim  to develop the turbulent filtration governing equation using averaging methods. Reading this paper: p. 75: "In this paper, we derive highwavenumberaveraged equations and highfrequency turbulence budget equations which hold in the entire space rather than only in the airspace of a porous medium". The latter statement is the correct note. Authors state incorrectly that: ''Time averaging followed by spatial averaging implies that obstacle elements interact only with timeaveraged flow and turbulence and do not interact with turbulent properties of the boundaylayer shear flow. The turbulence energycascade process is precluded under these assumptions.'' ''Therefore, schemes that first use the conventional timeaveraging method do not allow for detailed description of interactions between the wake turbulence and shear flow.''
Firstly: the assumptions mentioned in the text are not assumptions but merely applied averaging procedures.
Secondly: When writing these statements authors forget to admit that the whole idea of the VAT is to perform the modeling and mathematical description on both scales ( levels) of a medium. This implies that the corect modeling on the first level of hierarchy gives the correct physical picture  including also the boundarylayer features modeling within the fluid region of the REV and in vicinity to interface surface.
Thirdly: It is the situation when results and performance should be in clear correspondence with the admitted or expected physical and mathematical initial constructions. Meaning, that expectation of the results should not be evaluated on the basis of idealistical desires of improving the basics of turbulence theory using the VAT turbulent filtration governing equations and models, but strictly in agreement with undertaken or admitted just at the beginning the turbulence theory models themselves and their simplifications which are independently accepted in science many years ago for that or another purposes. The limitations, at this case, are directly inhereted from the initial turbulence assumptions  but not from the averaging procedures peculiarities.
Their averaged NS equation (15) is almost correct (excluding the second left hand term instead of correct ) by form when comparing with the VAT equation but actually this equation is incorrect because of the different definitions for their fluctuation variables . (In this above equation dropped two terms irrelavent to our analysis  influences of Coriolis and buoyancy forces). The fluctuations of velocity and other functions are not equal to intrinsic fluctuations nevertheless that both are the spacial fluctuations and according to their definition (10) but in VAT and as the consequence the averaging of this fluctuation gives which contradicts with their (11) making (11) incorrect ! Also their equation (17) is the equation which is defined as the kind of closure equation for the laminar fluctuations in the fluid region field and there is no trust to this equation  it has no integral terms, etc.? Assumption (18) which states  that the four terms in the equation (17) are the resistance terms  is just the pure arbitrary speculation. Authors used to represent the averaged velocity field as the sum of time averaged and instant averaged velocities and at the same time going from their previous definitions one can deduce that so or the full instant velocity is
Their mean motion equation (24) was developed on an easy basis because it obtains the form after simple time averaging of the previously averaged laminar regime equation WT1. This is really simple. The last term in this equation is the convenient resistance term commonly used in meteorology  which is the reflection of two last integral terms in the WT1 and had undertaken no change during the time averaging procedure  which is, of coarse, inappropriate for transient statements, but as a first approximation is acceptable.
Apart of the previous remark that this equation is incorrect due to incorrect treatment of fluctuation variable (plus incorrect writing of the second left hand term), still there are following questions:
1) The velocity field is not the turbulent velocity by definition, so what is it?
2) What is the relation of this velocity to the real turbulent velocity field happened to be in the porous media?
3) For the steady state regime the laminar velocity field that plays so important role in this equation development is the field without time dependent fluctuations so, where from the time fluctuations does appear in the further developed part of this work the turbulent form of equation ? But we know that there are the steady state mean turbulent flows exist ! And they have field's time fluctuations as the inseparable part of their physics.
4) The same kind of questions can be adressed to the issue of fluctuation variables , sense (meaning) ?
5) The impression is that in this development the porous medium itself is the major cause of the fluid turbulent motions which is not. Yes, it is the essential part of the turbulence characterized by porous medium, but not the major cause of the turbulent motion appearance.
6) The of the momentum equation (24) is not understandable at the moment  a) the fluctuation variable is the laminar regime fluctuation by definition (10); b) the fluctuation variable is not defined as usual instant point related velocity fluctuation, so it should be determined after laminar regime solution would be known.
When this equation compared with the VAT turbulent filtration (with the molecular viscosity terms neglected for simplicity) equation one would observe that in the latter equation for each variable exists the definition which is recognizable on each scale of hierarchy.
referring visitors to more observations which
can be found in  Fundamentals of Hierarchical Scaled Description in
Physics and Technology/ Are there any other Methods and Theories available?
http://travkinhspt.com/fundament/04.htm
Fluid flow in a porous layer or medium can be characterized by several modes. Let us single out from among them the three modes found in a highly porous media. The first is flow around isolated obstacle elements, or inside an isolated pore. The second is interaction of traces or a hyperturbulent mode. The third is fluid flow between obstacles or inside a blocked interconnected swarm of channels (filtration mode). The models developed by Scherban et al. (1986), Primak et al. (1986), and Travkin and Catton (1998) were primarily for laminar and turbulent nonlinear filtration and hyperturbulent modes transport. Specific features of flows in the channels of filtered media include the following: 1) increased drag due to microroughness on the channel boundary surfaces, 2) gravity effects for meteorological problems, 3) free convection effects, 4) the effects of secondary flows of the second kind and curved streamlines, 5) largescale vortex effects, and 6) the anisotropic nature of turbulent transfer and resulting anisotropy of turbulent viscosity.
It is wellknown that in spacial boundary flows, an important role is played by the gradients of normal Reynolds stresses and that this is the case for flows in porous medium channels as well. As a rule, flow symmetry is not observed in these channels. Therefore, in channel turbulence models, the shear components of the Reynolds stress tensors have a decisive effect on the flow characteristics. At present, however turbulence models that are less than second order can not be successfully employed for simulating such flows (Rodi, 1980; Lumley 1978; and Shvab and Bezprozvannykh, 1984). Derivation of the equations of turbulent flow and diffusion in a highly porous medium during the filtration mode is based on the theory of averaging of the turbulent transfer equation in the liquid phase and the transfer equations in the solid phase of a heterogeneous medium (Primak et al., (1986); and Shcherban et al., 1986) over a specified REV .
The initial turbulent transport equation set for the first level of the hierarchy, microelement or pore, was taken to be of the form ( see, for example Rodi, 1980; and Patel et al., 1985; among many others) Here and its fluctuation represent any scalar field that might be transported into either of the porous medium phases, and the last terms on the right hand side of (14) and (15) are source terms. Next we introduce free stream turbulence into the hierarchy. Let us represent the turbulent values as follows: where the index stands for the turbulent components independent of inhomogeneities of dimensions and properties of the multitude of porous media channels (pores), and stands for contributions due to the porous medium inhomogeneity. Being independent of the dimensions and properties of the inhomogeneities of the porous media configurations, sections and boundary surfaces does not mean that the distribution of values of , and are altogether independent of the distance to the wall, pressure distribution, etc. Thus, the values , or stand for the values generally accepted in the turbulence theory, i.e. when a plane surface is referred to, these values are those of a classical turbulent boundary layer. When a roundsection channel is involved, and even if the crosssection of this channel is not round, but without disturbing nonhomogeneities in the section, then the characteristics of this regular section (and flow ) may be considered to be those that could be marked with index . Hence, if a channel in a porous medium can be approximated by superposition of smooth regular (of regular shape ) channels, then it is possible to give such a flow its characteristics and designated them with the index , which stands for the basic (canonical) values of the turbulent quantities. Triple decomposition techniques have been used in papers by Brereton and Kodal (1992) and Bisset et al. (1991) among others. The latter utilized triple decomposition, conditional averaging and double averaging to analyze the structure of largescale organized motion over the rough plate.
It should be noted that there are problems where , and can be found from known theoretical or experimental expressions (correlations) where the definitions of , and are equivalent to the solution of an independent problem ( for example  turbulent flow in a curved channel). The same thing can be said about flow around a separate obstacle located on a plain surface. In this case one can write The term appears if the flow is through a nonuniform array of obstacles. If all the obstacles are the same and ordered, then can be taken equal to 0. Naturally, the term in this particular case doesn't equal the fluctuation vector over a smooth, plain surface. The following hypothes about the additive components is developed to correct the above deficiencies
It should be noted that solutions to the equations for the turbulent characteristics may be influenced by external parameters of the problem, namely, by the coefficients and boundary conditions which themselves can carry information about porous medium morphological features. The adoption of a hypothesis about the additive components of functions representing turbulent filtration facilitates the problem of averaging the equations for the Reynolds stresses and covariations of fluctuations (flows) in pores over the REV.
After averaging the basic initial set of turbulent transport equations over
the REV and using the averaging formalism developed in the works by Primak et
al. (1986), Shcherban et al. (1986), Primak and Travkin (1989), one obtains
the following equations for mass conservation
for turbulent filtration (with molecular viscosity terms neglected for simplicity)
and
for scalar diffusion (with molecular diffusivity terms neglected)
Many details and possible variants of the above equations with tensorial terms are found in Primak et al. (1986), Shcherban et al. (1986), Travkin and Catton (1995, 1998). Using an approximation to K  theory in an elementary channel (pore), the equation for turbulent diffusion of nth species takes the following more complex form after being averaged
In the more general case, the momentum flux integrals on the righthand sides of equations (16) through (17) do not equal zero, since there could be penetration through the phase transition boundary changing the boundary conditions in the microelement to allow for heat and mass exchange through the interface surface as the values of velocity, concentrations and temperature at do not equal zero (see also Crapiste et al., 1986). The first term on the righthand side of equation (18) is the divergence of the REV averaged product of velocity fluctuations and admixture concentration caused by random morphological properties of the medium being penetrated and is responsible for morphoconvectional dispersion of admixture in this particular porous medium. The third term on the righthand side of equation (18) can be associated with the notion of morphodiffusive dispersion of a substance or heat in a randomly nonhomogeneous medium. The term with may also reflect, specifically, the impact of microroughness from the previous level of the simulation hierarchy. The importance of accounting for this roughness has been demonstrated by many studies. The remaining step is to account for the microroughness characteristics of the previous level.
One dimensional mathematical statements will be used in what follows for simplicity. Admission of specific types of medium irregularity or randomness requires that complicated additional expressions be included in the generalized governing equations. Treatment of these additional terms becomes a crucial step once the governing averaged equations are written. An attempt to implement some basic departures from a porous medium with strictly regular morphology descriptions into a method for evaluation of some of the less tractable, additional terms is explained below.
The 1D momentum equation with terms representing a detailed description of the medium morphology is depicted as follows where is the turbulent eddy viscosity, is the square friction velocity at the upper boundary of surface roughness layer averaged over interface surface . General statements for energy transport in a porous medium require twotemperature treatments. Travkin et al. (1994), Gratton et al. (1995) showed that the proper form for the turbulent heat transfer equation in the fluid phase using oneequation Ktheory closure with primarily 1D convective heat transfer is while in the neighboring solid phase, the corresponding equation is or The generalized longitudinal 1D mass transport equation in the fluid phase, including description of potential morphofluctuation influences, for a medium morphology with only 1D fluctuations is written while the corresponding nonlinear equation for the solid phase is or
We recognize that the theory and the VAT mathematics those have been given in this subsection
are rather short for this formidable task, and hoping that the some additional materials shown in the following sections of this website 
and others, will fill some gaps in the exposition, while speaking on many more VAT turbulent and nonlinear governing equations those scattered throughout the website describing the specific problems.

thermal diffusivity
[]
 mean drag resistance coefficient in the REV []
 mean skin friction coefficient over the turbulent area of
[]
 mean form resistance coefficient in the REV []
 drag resistance coefficient upon single sphere []
 mean skin friction coefficient over the laminar region inside of the REV
[]
 specific heat
[]
 constant coefficient in Kolmogorov turbulent exchange coefficient
correlation []
 character pore size in the cross section
[]
 diameter [m] of ith pore
[]
 particle diameter
[]
 interphase differential area in porous medium
[]
 molecular diffusion coefficient
],
also  tube or pore diameter
[]
 flat channel hydraulic diameter
[]
diffusion coefficient in
solid[]
 internal surface in the REV
[]
 averaged over
value
 intrinsic averaged variable
 value f, averaged over
in a REV  phase averaged variable
 value
,
averaged over
in a REV
 morphofluctuation value of
in a
 gravitational constant
[]
 width of the channel
[]
 averaged heat transfer coefficient over
[],
halfwidth of the channel
[]
 pore scale microroughness layer thickness
[]
 internal surface in the REV
[]
 fluid thermal conductivity
[]
 solid phase thermal conductivity
[]
 permeability
[]
 turbulent kinetic energy exchange coefficient
[]
 turbulent diffusion coefficient
[]
 turbulent eddy viscosity
[]
 effective thermal conductivity of solid phase
[]
 turbulent eddy thermal conductivity
[]
 turbulence mixing length
[]
 scale
[]
 averaged porosity []
 surface porosity []
 number of pores []
 number of pores with diameter of type i []

=,
interface surface Nusselt number []
 pressure
[];
or pitch in regular porous 2D and 3D medium
[];
or phase function []

=,
Darcy velocity pore scale Peclet number []

=particle
radius Peclet number []

=
,
Prandtl number []
 outward heat flux
[]
 Reynolds number of pore hydraulic diameter []

=
Darcy velocity Reynolds number of pore hydraulic diameter []

=,
particle Reynolds number []

=
Reynolds number of general scale pore hydraulic diameter []
 total cross sectional area available to flow
[]
 specific surface of a porous medium
[]
 cross flow projected area of obstacles
[]
 temperature
 characteristic temperature for given temperature range
 solid phase temperature
 wall temperature
 reference temperature
 velocity in xdirection
 square friction velocity at the upper boundary of hr averaged over surface
[]
 velocity

=Darcy
velocity
[]
 velocity in zdirection

averaged heat transfer coefficient over

representative elementary volume (REV)
 pore volume in a REV
 solid phase volume in a REV
electric permittivity [Fr/m]
 dynamic viscosity
[]
or
[]
 magnetic permeability [H/m]
 kinematic viscosity
[];
also
 frequency [Hz]
 density
[];
also
 electric charge density
[C/m]
 medium specific electric conductivity [A/V/m]
 electric scalar potential [V]
 particle intensity per unit energy (frequency)
 ensembleaveraged value of
 interface ensembleaveraged value of
with phase
being to the left
 angular frequency [rad/s]
 magnetic susceptibility []
=
 absorption coefficient [1/m]
=
 scattering coefficient [1/m]
Antohe, B.V. and Lage, J.L. (1997), "A General TwoEquation Macroscopic Turbulence Model for Incompressible Flow in Porous Media", Int. J. Heat Mass Transfer, Vol. 40, No. 13, pp. 30133024.
Bisset, D. K., Antonia, R.A., and Raupach, M.R., (1991). Topology and transport properties of largescale organized motion in a slightly heated rough wall boundary layer. Phys. Fluids A 3, 22202228.
Brereton, G. J. and Kodal, A. (1992). A frequencydomain filtering technique for triple decomposition of unsteady turbulent flow. J. Fluids Engin. 114, 4551.
Catton, I. and Travkin, V. S. (1996). Turbulent flow and heat transfer in high permeability porous media. In Proc. Intern. Conf. on Porous Media and Their Applic. Science, Engineer. and Ind., (K.Vafai and P.N.Shivakumar, eds.), Engin. Found. & Inst. Ind. Math. Sc., pp. 333391, New York.
Coppin, P. A., Raupach, M. R., and Legg, B. J. (1986). Experiments on scalar dispersion within a model plant canopy. Part II: an elevated plane source. BoundaryLayer Meteorology 35, 167191.
Crapiste, G.H., Rotstein, E., and Whitaker, S. (1986). A general closure scheme for the method of volume averaging. Chem. Engn. Sci. 41, 227235.
Dybbs, A., Edwards, R. V. (1982). A new look at porous media fluid mechanics Darcy to turbulent. In Proc. NATO Advanced Study Institution on Mechanics of Fluids in Porous Media, NATO ASI Series E 82, pp. 201256.
Fand, R.M., Kim, B.Y.K., Lam, A.C.C., and Phan, R.T. (1987). Resistance to the flow of fluids through simple and complex porous media whose matrices are composed of randomly packed spheres. J. Fluids Engin. 109, 268274.
Gratton, L., Travkin, V.S., and Catton, I. (1995). The impact of morphology irregularity on bulk flow and twotemperature heat transport in highly porous media. In Proc. A.S.M.E./J.S.M.E. Thermal Engn. Joint Conf. 3, pp. 339346.
Gratton, L., Travkin, V.S., and Catton, I. (1996). The influence of morphology upon twotemperature statements for convective transport in porous media. J. Enhanced Heat Transfer 3, 129145.
Kheifets, L.I. and Neimark, A.V., Multiphase Processes in Porous Media , Nadra, Moscow, 1982.
Howle, L., Behringer, R. P., and Georgiadis, J. G. (1992). Pattern formation near the onset of convection for fluid in a porous medium. Private communication.
Ishii (1975). Thermofluid dynamic theory of twophase flow. Eyrolles, Paris.
Ishii and Mishima, K. (1984). Twofluid model and hydrodynamic constitutive relations. Nuclear Engineering and Design 82, 107126.
Khan, E.U., Rohsenow, W. M., Sonin, A. A., and Todreas, N.E. (1975). A porous body model for predicting temperature distribution in wirewrapped rod assemblies operating in combined forced and free convection. Nuclear Engineering and Design 35, 199211.
Lahey, T. Jr. and Lopez de Bertodano, M. (1991). The prediction of phase distribution using twofluid models. In Proc. ASME/JSME Thermal Engineering Conf. 2, pp. 193200.
Lahey, R.T., Jr. and Drew, D.A., (1988). The threedimensional time and volume averaged conservation equations of twophase flow. In Advances in Nuclear Science and Technology, (Lewins and Becker, eds.), 20, pp. 169.
Legg, B. J., Raupach, M. R., and Coppin, P. A. (1986). Experiments on scalar dispersion within a model plant canopy. Part III: an elevated line source", BoundaryLayer Meteorology 35, 277302.
Lopez de Bertodano, M., Lee, SJ., Lahey, R. T. Jr., and Drew, D.A. (1990). The prediction of twophase turbulence and phase distribution phenomena using a Reynolds stress model. J. Fluids Engin. 112, 107113.
Lumley, J. L. (1978). Computational modelling of turbulent flows. 1Advances in Applied Mechanics 18, 123176.
Masuoka, T. and Takatsu, Y. (1996). Turbulence model for flow through porous media. Int. J. Heat Mass Transfer 39, 28032809.
Patel, V. C., Rodi, W., and Scheurer, G. (1985). Turbulence models for nearwall and low Reynolds number flows: a review. AIAA J. 23, 13081319.
Popov, A. M. (1974). On peculiarities of atmospheric diffusion over inhomogeneous surface. Izv. AN SSSR, AOPh. 10, 13091312, (in Russian).
Popov, A. M. (1975). Atmospheric boundary layer simulation within the roughness layer. Izv. AN SSSR, AOPh. 11, 574581, (in Russian).
Primak, A.V., Shcherban, A.N. and Travkin, V.S. (1986). Turbulent transfer in urban agglomerations on the basis of experimental statistical models of roughness layer morphological properties. In Transactions World Meteorological Organization Conference on Air Pollution Modelling and its Application, 2, pp. 259266, WMO, Geneva.
Primak, A.V. and Travkin, V.S. (1989). Simulation of turbulent transfer of meteoelements and pollutants under conditions of artificial anthropogenic action in urban roughness layer as in sorbing and biporous twophase medium. In Intern. Sym. Interconnec. Problems Hydrological Cycle and Atmospheric Proc. Under Conditions Anthropogenic Influences, Trans., Schopron.
Raupach. M. R. and Shaw, R. H. (1982). Averaging procedures for flow within vegetation canopies. BoundaryLayer Meteorology 22, 7990 .
Raupach, M. R., Coppin, P. A., and Legg, B. J. (1986). Experiments on scalar dispersion within a model plant canopy. Part I: the turbulence structure. BoundaryLayer Meteorology 35, 2152.
Rodi, W. (1980). Turbulence models for environmental problems. In Prediction Methods for Turbulent Flows, (W. Kollmann, ed.), pp. 259350. Hemisphere Publishing Corporation, New York
Scherban, A. N., Primak, A.V., and Travkin, V. S. (1986). Mathematical models of flow and mass transfer in urban roughness layer. Problemy Kontrolya i Zaschita Atmosfery ot Zagryazneniya 12, 310, (in Russian).
Shvab, V. A. and Bezprozvannykh, V.A. (1984). On turbulent flow simulation in rectilinear channels of noncircular crosssection. In Metody Aerodinamiki i Teplomassoobmena v Tekhnologicheskikh Protsessakh, pp. 325. Izdatelstvo TGU, Tomsk, (in Russian).
Subbotin, V. I., Kashcheev, V. M., Nomofilov, E. V., and Yur'ev, Yu.S. (1979). Computer Problem Solving in Nuclear Reactor Thermophysics. Atomizdat, Moscow, (in Russian).
Takatsu, Y., and Masuoka, T., "Turbulent phenomena in flow through porous media", Journal of Porous Media, Vol. 1, No. 3, pp.243251, (1997).
Teyssedou, A. Tapucu, and Camarero, R. (1992). Blocked flow subchannel simulation comparison with singlephase flow data. J. Fluids Engin. 114, 205213.
Travkin, V.S. and Catton, I. (1992). Models of turbulent thermal diffusivity and transfer coefficients for a regular packed bed of spheres. In Fundamentals of Heat Transfer in Porous Media (M. Kaviany, ed.), ASME HTD  193, pp.1523.
Travkin, V.S., Catton, I., and Gratton, L. (1993). Single phase turbulent transport in prescribed nonisotropic and stochastic porous media. In Heat Transfer in Porous Media, ASME HTD240, pp. 4348.
Travkin, V.S. and Catton, I. (1994). Turbulent transport of momentum, heat and mass in a two level highly porous media. In Heat Transfer 1994, Proc. Tenth Intern. Heat Transfer Conf. (G.F.Hewitt, ed.), 5, pp. 399404. Chameleon Press, London.
Travkin, V.S., Gratton, L., and Catton, I. (1994). A morphologicalapproach for twophase porous mediumtransport and optimum design applications in energy engineering. In Proc. 12th Symp. Energy Engin. Sciences, Argonne National Laboratory, Conf. 9404137, pp. 4855.
Travkin, V.S. and Catton, I. (1995). A twotemperature model for turbulent flow and heat transfer in a porous layer. J. Fluids Eng. 117, 181188.
Travkin, V. S., Catton, I. and Hu, K. (1998). Channel flow In porous media in the limit as porosity approaches unity. In Proc. ASME HTD3611, pp. 277284.
Travkin, V.S. and Catton, I. (1998). Porous media transport descriptions  nonlocal, linear and nonlinear against effective thermal/fluid properties. Advances in Colloid and Interface Science 7677, 389443.
Travkin, V. S., Hu, K., and Catton, I. (1999). Turbulent kinetic energy and dissipation rate equation models for momentum transport in porous media. In Proc. 3rd ASME/JSME Fluids Engineering Conf.  FEDSM997275, ASME, San Francisco.
Travkin, V.S. and Catton, I. (1999a). Nonlinear effects in multiple regime transport of momentum in longitudinal capillary porous medium morphology, J. Porous Media., 2, 277294.
Travkin, V.S. and Catton, I. (1999b). Critique of theoretical models of transport phenomena in heterogeneous media, (invited). Presentation at the 3rd ASME/JSME Fluids Engineering Conf.  FEDSM997922, July 1823, 1999, San Francisco.
Vafai, K., and Tien, C.L. (1981). Boundary and inertia effects on flow and heat transfer in porous media. Intern. J. Heat and Mass Transfer 24, 195203.
Wang, H. and Takle, E.S. (1995), "BoundaryLayer Flow and Turbulence Near Porous Obstacles", BoundaryLayer Meteorology, Vol. 74, pp. 7388.
Whitaker, S. (1986a). Flow in porous media I: A theoretical derivation of darcy's law. Transport in Porous Media 1, 325.
Whitaker, S. (1986b). Flow in porous media II: The governing equations for immiscible, twophase flow. Transport in Porous Media, 1, 105125.
Yamada, T. (1982). A numerical model study of turbulent airflow in and above a forest canopy. J. Meteorol. Soc. Jap. 60, 439454.
Zhang, D.Z. and Prosperetti, A. (1994). Averaged equations for inviscid disperse twophase flow. J. Fluid Mech. 267, 185219.