Modeling the Transport of Fuel Mixture Perturbations and Entropy Waves in the Linearized Framework

Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner

Abstract In this paper a new method is introduced to model the transport of entropy waves and equivalence ratio fluctuations in turbulent flows. The model is based on the Navier-Stokes equations and includes a transport equation for a passive scalar, which may stand for entropy or equivalence ratio fluctuations. The equations are linearized around the mean turbulent fields. These serve as the input to the model in addition to a turbulent eddy viscosity, which accounts for turbulent diffusion of the perturbations. Based on these inputs, the framework is able to predict the linear response of the flow velocity and passive scalar to harmonic perturbations that are imposed at the boundaries of the computational domain. These, in this study, are fluctuations in the passive scalar and/or velocities at the inlet of a channel flow. The code is first validated against analytic results, showing very good agreement. Then the method is applied to predict the convection, mean flow dispersion and turbulent mixing of passive scalar fluctuations in a turbulent channel flow, which has been studied in previous work with Direct Numerical Simulations (DNS). Results show that our code reproduces the dynamics of coherent passive scalar transport in the DNS with very high accuracy and low numerical costs. Furthermore, we demonstrate that turbulent mixing has a significant effect on the transport of the passive scalar fluctuations. Finally, we apply the method to explain experimental observations of transport of equivalence ratio fluctuations in the mixing duct of a model burner.

Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner

Abstract In this paper a new method is introduced to model the transport of entropy waves and equivalence ratio fluctuations in turbulent flows. The model is based on the Navier-Stokes equations and includes a transport equation for a passive scalar, which may stand for entropy or equivalence ratio fluctuations. The equations are linearized around the mean turbulent fields, which serve as the input to the model in addition to a turbulent eddy viscosity, which accounts for turbulent diffusion of the perturbations. Based on these inputs, the framework is able to predict the linear response of the flow velocity and passive scalar to harmonic perturbations that are imposed at the boundaries of the computational domain. These in this study are fluctuations in the passive scalar and/or velocities at the inlet of a channel flow. The code is first validated against analytic results, showing very good agreement. Then the method is applied to predict the convection, mean flow dispersion and turbulent mixing of passive scalar fluctuations in a turbulent channel flow, which has been studied in previous work with Direct Numerical Simulations (DNS). Results show that our code reproduces the dynamics of coherent passive scalar transport in the DNS with very high accuracy and low numerical costs, when the DNS mean flow and Reynolds stresses are provided. Furthermore, we demonstrate that turbulent mixing has a significant effect on the transport of the passive scalar fluctuations. Finally, we apply the method to explain experimental observations of transport of equivalence ratio fluctuations in the mixing duct of a model burner.


2012 ◽  
Vol 712 ◽  
pp. 169-202 ◽  
Author(s):  
A. Busse ◽  
N. D. Sandham

AbstractThe effects of rough surfaces on turbulent channel flow are modelled by an extra force term in the Navier–Stokes equations. This force term contains two parameters, related to the density and the height of the roughness elements, and a shape function, which regulates the influence of the force term with respect to the distance from the channel wall. This permits a more flexible specification of a rough surface than a single parameter such as the equivalent sand grain roughness. The effects of the roughness force term on turbulent channel flow have been investigated for a large number of parameter combinations and several shape functions by direct numerical simulations. It is possible to cover the full spectrum of rough flows ranging from hydraulically smooth through transitionally rough to fully rough cases. By using different parameter combinations and shape functions, it is possible to match the effects of different types of rough surfaces. Mean flow and standard turbulence statistics have been used to compare the results to recent experimental and numerical studies and a good qualitative agreement has been found. Outer scaling is preserved for the streamwise velocity for both the mean profile as well as its mean square fluctuations in all but extremely rough cases. The structure of the turbulent flow shows a trend towards more isotropic turbulent states within the roughness layer. In extremely rough cases, spanwise structures emerge near the wall and the turbulent state resembles a mixing layer. A direct comparison with the study of Ashrafian, Andersson & Manhart (Intl J. Heat Fluid Flow, vol. 25, 2004, pp. 373–383) shows a good quantitative agreement of the mean flow and Reynolds stresses everywhere except in the immediate vicinity of the rough wall. The proposed roughness force term may be of benefit as a wall model for direct and large-eddy numerical simulations in cases where the exact details of the flow over a rough wall can be neglected.


1984 ◽  
Vol 148 ◽  
pp. 413-442 ◽  
Author(s):  
S. Biringen

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 1500. Three-dimensional incompressible Navier–Stokes equations are numerically integrated to obtain the time evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 × 51 × 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and characteristic of the various stages of transition that lead to final breakdown. In particular, evidence points to the formation of a A-shaped vortex and the subsequent system of horsehoe vortices inclined to the main flow direction as the primary elements of transition. Details of the resulting flow field after breakdown indicate the evolution of streaklike formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), the introduction of subgrid-scale terms seems necessary to obtain fully developed turbulence statistics.


2019 ◽  
Vol 865 ◽  
pp. 546-562 ◽  
Author(s):  
Giovanni Iacobello ◽  
Stefania Scarsoglio ◽  
J. G. M. Kuerten ◽  
Luca Ridolfi

A temporal complex network-based approach is proposed as a novel formulation to investigate turbulent mixing from a Lagrangian viewpoint. By exploiting a spatial proximity criterion, the dynamics of a set of fluid particles is geometrized into a time-varying weighted network. Specifically, a numerically solved turbulent channel flow is employed as an exemplifying case. We show that the time-varying network is able to clearly describe the particle swarm dynamics, in a parametrically robust and computationally inexpensive way. The network formalism enables us to straightforwardly identify transient and long-term flow regimes, the interplay between turbulent mixing and mean flow advection and the occurrence of proximity events among particles. Thanks to their versatility and ability to highlight significant flow features, complex networks represent a suitable tool for Lagrangian investigations of turbulent mixing. The present application of complex networks offers a powerful resource for Lagrangian analysis of turbulent flows, thus providing a further step in building bridges between turbulence research and network science.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


1971 ◽  
Vol 47 (2) ◽  
pp. 321-335 ◽  
Author(s):  
Jon Lee

The truncated Burgers models have a unique equilibrium state which is defined continuously for all the Reynolds numbers and attainable from a realizable class of initial disturbances. Hence, they represent a sequence of convergent approximations to the original (untruncated) Burgers problem. We have pointed out that consideration of certain degenerate equilibrium states can lead to the successive turbulence-turbulence transitions and finite-jump transitions that were suggested by Case & Chiu. As a prototype of the Navier–Stokes equations, Burgers model can simulate the initial-value type of numerical integration of the Fourier amplitude equations for a turbulent channel flow. Thus, the Burgers model dynamics display certain idiosyncrasies of the actual channel flow problem described by a truncated set of Fourier amplitude equations, which includes only a modest number of modes due to the limited capability of the computer at hand.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


2014 ◽  
Vol 16 (5) ◽  
pp. 901-918 ◽  

<div> <p>Three-dimensional calculations were performed to simulate the flow around a cylindrical vegetation element using the Scale Adaptive Simulation (SAS) model; commonly, this is the first step of the modeling of the flow through multiple vegetation elements. SAS solves the Reynolds Averaged Navier-Stokes equations in stable flow regions, while in regions with unstable flow it goes unsteady producing a resolved turbulent spectrum after reducing eddy viscosity according to the locally resolved vortex size represented by the von Karman length scale. A finite volume numerical code was used for the spatial discretisation of the rectangular computational domain with stream-wise, cross-flow and vertical dimensions equal to 30D, 11D and 1D, respectively, which was resolved with unstructured grids. Calculations were compared with experiments and Large Eddy Simulations (LES). Predicted overall flow parameters and mean flow velocities exhibited a very satisfactory agreement with experiments and LES, while the agreement of predicted turbulent stresses was satisfactory. Calculations showed that SAS is an efficient and relatively fast turbulence modeling approach, especially in relevant practical problems, in which the very high accuracy that can be achieved by LES at the expense of large computational times is not required.</p> </div> <p>&nbsp;</p>


2013 ◽  
Vol 721 ◽  
pp. 541-577 ◽  
Author(s):  
Amin Rasam ◽  
Geert Brethouwer ◽  
Arne V. Johansson

AbstractIn Marstorpet al. (J. Fluid Mech., vol. 639, 2009, pp. 403–432), an explicit algebraic subgrid stress model (EASSM) for large-eddy simulation (LES) was proposed, which was shown to considerably improve LES predictions of rotating and non-rotating turbulent channel flow. In this paper, we extend that work and present a new explicit algebraic subgrid scalar flux model (EASSFM) for LES, based on the modelled transport equation of the subgrid-scale (SGS) scalar flux. The new model is derived using the same kind of methodology that leads to the explicit algebraic scalar flux model of Wikströmet al. (Phys. Fluids, vol. 12, 2000, pp. 688–702). The algebraic form is based on a weak equilibrium assumption and leads to a model that depends on the resolved strain-rate and rotation-rate tensors, the resolved scalar-gradient vector and, importantly, the SGS stress tensor. An accurate prediction of the SGS scalar flux is consequently strongly dependent on an accurate description of the SGS stresses. The new EASSFM is therefore primarily used in connection with the EASSM, since this model can accurately predict SGS stresses. The resulting SGS scalar flux is not necessarily aligned with the resolved scalar gradient, and the inherent dependence on the resolved rotation-rate tensor makes the model suitable for LES of rotating flow applications. The new EASSFM (together with the EASSM) is validated for the case of passive scalar transport in a fully developed turbulent channel flow with and without system rotation. LES results with the new model show good agreement with direct numerical simulation data for both cases. The new model predictions are also compared to those of the dynamic eddy diffusivity model (DEDM) and improvements are observed in the prediction of the resolved and SGS scalar quantities. In the non-rotating case, the model performance is studied at all relevant resolutions, showing that its predictions of the Nusselt number are much less dependent on the grid resolution and are more accurate. In channel flow with wall-normal rotation, where all the SGS stresses and fluxes are non-zero, the new model shows significant improvements over the DEDM predictions of the resolved and SGS quantities.


Sign in / Sign up

Export Citation Format

Share Document