Analysis of an RNG based turbulence model for separated flows

A two-equation turbulence model of the K-epsilon type was recently derived by using Renormalization Group (RNG) methods. It was later reported that this RNG based model yields substantially better predictions than the standard K-epsilon model for turbulent flow over a backward facing step - a standard test case used to benchmark the performance of turbulence models in separated flows. The improvements obtained from the RNG K-epsilon model were attributed to the better treatment of near wall turbulence effects. In contrast to these earlier claims, it is shown in this paper that the original version of the RNG K-epsilon model substantially underpredicts the reattachment point in the backstep problem. This is a deficiency that is traced to the modeling of the production of dissipation term. However, with the most recent improvements in the RNG K-epsilon model, excellent results for the backstep problem are now obtained

Turbulent separated flow past a backward-facing step: A critical evaluation of two-equation turbulence models

The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem

Bounded energy states in homogeneous turbulent shear flow: An alternative view

The equilibrium structure of homogeneous turbulent shear flow is investigated from a theoretical standpoint. Existing turbulence models, in apparent agreement with physical and numerical experiments, predict an unbounded exponential time growth of the turbulent kinetic energy and dissipation rate; only the anisotropy tensor and turbulent time scale reach a structural equilibrium. It is shown that if vortex stretching is accounted for in the dissipation rate transport equation, then there can exist equilibrium solutions, with bounded energy states, where the turbulence production is balanced by its dissipation. Illustrative calculations are present for a k-epsilon model modified to account for vortex stretching. The calculations indicate an initial exponential time growth of the turbulent kinetic energy and dissipation rate for elapsed times that are as large as those considered in any of the previously conducted physical or numerical experiments on homogeneous shear flow. However, vortex stretching eventually takes over and forces a production-equals-dissipation equilibrium with bounded energy states. The validity of this result is further supported by an independent theoretical argument. It is concluded that the generally accepted structural equilibrium for homogeneous shear flow with unbounded component energies is in need of re-examination

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion

On the geometry of loop quantum gravity on a graph

We discuss the meaning of geometrical constructions associated to loop quantum gravity states on a graph. In particular, we discuss the "twisted geometries" and derive a simple relation between these and Regge geometries.Comment: 6 pages, 1 figure. v2: some typos corrected, references update

A semiclassical tetrahedron

We construct a macroscopic semiclassical state state for a quantum tetrahedron. The expectation values of the geometrical operators representing the volume, areas and dihedral angles are peaked around assigned classical values, with vanishing relative uncertainties.Comment: 10 pages; v2 revised versio

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to two severe test cases. The turbulent flows considered include the incompressible flat plate boundary layer with the adverse pressure gradients and incompressible flow past a backward facing step. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data

Twisted geometries: A geometric parametrisation of SU(2) phase space

A cornerstone of the loop quantum gravity program is the fact that the phase space of general relativity on a fixed graph can be described by a product of SU(2) cotangent bundles per edge. In this paper we show how to parametrize this phase space in terms of quantities describing the intrinsic and extrinsic geometry of the triangulation dual to the graph. These are defined by the assignment to each triangle of its area, the two unit normals as seen from the two polyhedra sharing it, and an additional angle related to the extrinsic curvature. These quantities do not define a Regge geometry, since they include extrinsic data, but a looser notion of discrete geometry which is twisted in the sense that it is locally well-defined, but the local patches lack a consistent gluing among each other. We give the Poisson brackets among the new variables, and exhibit a symplectomorphism which maps them into the Poisson brackets of loop gravity. The new parametrization has the advantage of a simple description of the gauge-invariant reduced phase space, which is given by a product of phase spaces associated to edges and vertices, and it also provides an abelianisation of the SU(2) connection. The results are relevant for the construction of coherent states, and as a byproduct, contribute to clarify the connection between loop gravity and its subset corresponding to Regge geometries.Comment: 28 pages. v2 and v3 minor change

Physical boundary state for the quantum tetrahedron

We consider stability under evolution as a criterion to select a physical boundary state for the spinfoam formalism. As an example, we apply it to the simplest spinfoam defined by a single quantum tetrahedron and solve the associated eigenvalue problem at leading order in the large spin limit. We show that this fixes uniquely the free parameters entering the boundary state. Remarkably, the state obtained this way gives a correlation between edges which runs at leading order with the inverse distance between the edges, in agreement with the linearized continuum theory. Finally, we give an argument why this correlator represents the propagation of a pure gauge, consistently with the absence of physical degrees of freedom in 3d general relativity.Comment: 20 pages, 6 figure

Grasping rules and semiclassical limit of the geometry in the Ponzano-Regge model

We show how the expectation values of geometrical quantities in 3d quantum gravity can be explicitly computed using grasping rules. We compute the volume of a labelled tetrahedron using the triple grasping. We show that the large spin expansion of this value is dominated by the classical expression, and we study the next to leading order quantum corrections.Comment: 18 pages, 1 figur