A presumed joint pdf model for turbulent combustion with varying equivalence ratio

Modeling of the joint probability density function of the mixture fraction and progress variable with a given covariance value is studied. This modeling is validated using experimental and direct numerical simulation (DNS) data. A very good agreement with experimental data of turbulent stratified flames and DNS data of a lifted hydrogen jet flame is obtained. The effect of using this joint pdf modeling to calculate the mean reaction rate with a flamelet closure in Reynolds averaged Navier-Stokes (RANS) calculation of stratified flames is studied. The covariance effect is observed to be large within the flame brush. The results obtained from RANS calculations using this modeling for stratified jet- and rod-stabilized V-flames are discussed and compared to the measurements as a posteriori validation for the joint probability density function model with the flamelet closure. The agreement between the computed and measured values of flame and turbulence quantities is found to be good. © 2012 Copyright Taylor and Francis Group, LLC.

Bayesian non-parametrics and the probabilistic approach to modelling.

Modelling is fundamental to many fields of science and engineering. A model can be thought of as a representation of possible data one could predict from a system. The probabilistic approach to modelling uses probability theory to express all aspects of uncertainty in the model. The probabilistic approach is synonymous with Bayesian modelling, which simply uses the rules of probability theory in order to make predictions, compare alternative models, and learn model parameters and structure from data. This simple and elegant framework is most powerful when coupled with flexible probabilistic models. Flexibility is achieved through the use of Bayesian non-parametrics. This article provides an overview of probabilistic modelling and an accessible survey of some of the main tools in Bayesian non-parametrics. The survey covers the use of Bayesian non-parametrics for modelling unknown functions, density estimation, clustering, time-series modelling, and representing sparsity, hierarchies, and covariance structure. More specifically, it gives brief non-technical overviews of Gaussian processes, Dirichlet processes, infinite hidden Markov models, Indian buffet processes, Kingman's coalescent, Dirichlet diffusion trees and Wishart processes.

Zonal RANS-LES modeling for turbines in aeroengines

The cost of large-eddy simulation (LES) modeling in various zones of gas turbine aeroengines is outlined. This high cost clearly demonstrates the need to perform hybrid Reynolds-averaged Navier-Stokes-LES (RANS-LES) over the majority of engine zones because the Reynolds number is too high for pure LES. The RANS layer is used to cover over the fine streaks found in the inner part of the boundary layer. The hybrid strategy is applied to various engine zones, which is shown to typically give much greater predictive accuracy than pure RANS simulations. However, the cost estimates show that the RANS layer should be disposed within the low-pressure turbine zone. Also, the nature of the flow physics in this zone makes LES most sensible. © 2012 by Begell House, Inc.

The influence of geometry on jet plume development

Our recent efforts of using large-eddy simulation (LES) type methods to study complex and realistic geometry single stream and co-flow nozzle jets and acoustics are summarized in this paper. For the LES, since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving a numerical type LES (NLES). To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier-Stokes) model is smoothly blended in the LES making a hybrid RANS-NLES approach. Several complex nozzle geometries including the serrated (chevron) nozzle, realistic co-axial nozzles with eccentricity, pylon and wing-flap are discussed. The hybrid RANS-NLES simulations show encouraging predictions for the chevron jets. The chevrons are known to increase the high frequency noise at high polar angles, but decrease the low frequency noise at lower angles. The deflection effect of the potential core has an important mechanism of noise reduction. As for co-axial nozzles, the eccentricity, the pylon and the deployed wing-flap are shown to influence the flow development, especially the former to the length of potential core and the latter two having a significant impact on peak turbulence levels and spreading rates. The studies suggest that complex and real geometry effects are influential and should be taken into count when moving towards real engine simulations. © 2012 Elsevier Ltd. All rights reserved.

LES of heat transfer in electronics

Numerical methods based on the Reynolds Averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) equations are applied to the thermal prediction of flows representative of those found in and around electronics systems and components. Low Reynolds number flows through a heated ribbed channel, around a heated cube and within a complex electronics system case are investigated using linear and nonlinear LES models, hybrid RANS-LES and RANS-Numerical-LES (RANS-NLES) methods. Flow and heat transfer predictions using these techniques are in good agreement with each other and experimental data for a range of grid resolutions. Using second order central differences, the RANS-NLES method performs well for all simulations. © 2011 Elsevier Inc.

Numerical simulation of single-stream jets from a serrated nozzle

Hybrid large-eddy type simulations for cold jet flows from a serrated nozzle are performed at an acoustic Mach number Ma ac = 0.9 and Re = 1.03×10 6. Since the solver being used tends towards having dissipative qualities, the subgrid scale (SGS) model is omitted, giving a numerical type LES (NLES) or implicit LES (ILES) reminiscent procedure. To overcome near wall streak resolution problems a near wall RANS (Reynolds averaged Navier-Stokes) model is smoothly blended to the LES making a hybrid RANS-ILES. The geometric complexity of the serrated nozzle is fully considered without simplification or emulation. An improved but still modest hexahedral multi-block grid with circa 20 million grid points (with respect to 12.5 million in Xia et al.; Int J Heat Fluid Flow 30:1067-1079, 2009) is used. Despite the modest grid size, encouraging and improved results are obtained. Directly resolved mean and second-order fluctuating quantities along the jet centerline and in the jet shear layer compare favorably with measurements. The radiated far-field sound predicted using the Ffowcs Williams and Hawkings (FW-H) surface integral method shows good agreement with the measurements in directivity and sound spectra. © 2011 Springer Science+Business Media B.V.

Comparison of les to LDA and PIV measurements of a small scale high speed coflowing jet

Measurements and predictions are made of a short cowl co-flowing jet with a bypass ratio of 8:1. The Reynolds number for computations and measurements are matched at 300,000 and the Mach numbers representative of realistic jet conditions with core and co flow velocities of 240m/s and 216m/s respectively. The low Reynolds number of the measurements makes the case well suited to the assessment of large eddy resolving computational strategies. Also, the nozzle concentricity was carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both two dimensional coincident LDA and PIV systems. The computational simulations are completed on a modest 12×106 mesh. The simulation is now being run on a 50×106 mesh using hybrid RANSNLES (Numerical Large Eddy Simulation). Close to the nozzle wall a k-l RANS model is used. For an axisymmetric jet, comparison is made between simulations which use NLES, RANSNLES and also a simple imposed velocity profile where the nozzle is not modeled. The use of a near wall RANS model is shown to be beneficial. When compared with the measurements the NLES results are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Strategies for modeling turbulent flows in electronics

Hybrid methods based on the Reynolds Averaged Navier Stokes (RANS) equations and the Large Eddy Simulation (LES) formulation are investigated to try and improve the accuracy of heat transfer and surface temperature predictions for electronics systems and components. Two relatively low Reynolds number flows are studied using hybrid RANS-LES, RANS-Implicit-LES (RANS-ILES) and non-linear LES models. Predictions using these methods are in good agreement with each other, even using different grid resolutions. © 2008 IEEE.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Predictive techniques for turbulent oscillatory flows

The prediction of turbulent oscillatory flow at around transitional Reynolds numbers is considered for an idealized electronics system. To assess the accuracy of turbulence models, comparison is made with measurements. A stochastic procedure is used to recover instantaneous velocity time traces from predictions. This procedure enables more direct comparison with turbulence intensity measurements which have not been filtered to remove the oscillatory flow component. Normal wall distances, required in some turbulence models, are evaluated using a modified Poisson equation based technique. A range of zero, one and two equation turbulence models are tested, including zonal and a non-linear eddy viscosity models. The non-linear and zonal models showed potential for accuracy improvements.

High-order DES and zonal LES of a labyrinth seal

Hybrid numerical large eddy simulation (NLES) and detached eddy simulation (DES) methods are assessed on a labyrinth seal geometry. A high sixth order discretization scheme is used and is validated using a test case of a two dimensional vortex. The hybrid approach adopts a new blending function and along with DES is initially validated using a simple cavity flow. The NLES method is also validated outside of RANS zones. It is found that there is very little resolved turbulence in the cavity for the DES simulation. For the labyrinth seal calculations the DES approach is problematic giving virtually no resolved turbulence content. It is seen that over the tooth tips the extent of the LES region is small and is likely to be a strong contributor to excessive flow damping in these regions. On the other hand the zonal Hamilton-Jacobi approach did not suffer from this trait. In both cases the meshes used are considered to be hybrid RANS-LES adequate. Fortunately (or perhaps unfortunately) the DES profiles are in agreement with the time mean experimental measurements. It is concluded that for an inexperienced CFD practitioner this could have wider implications particularly if transient results such as unsteady loading are desired. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Functional requirements for reward-modulated spike-timing-dependent plasticity.

Recent experiments have shown that spike-timing-dependent plasticity is influenced by neuromodulation. We derive theoretical conditions for successful learning of reward-related behavior for a large class of learning rules where Hebbian synaptic plasticity is conditioned on a global modulatory factor signaling reward. We show that all learning rules in this class can be separated into a term that captures the covariance of neuronal firing and reward and a second term that presents the influence of unsupervised learning. The unsupervised term, which is, in general, detrimental for reward-based learning, can be suppressed if the neuromodulatory signal encodes the difference between the reward and the expected reward-but only if the expected reward is calculated for each task and stimulus separately. If several tasks are to be learned simultaneously, the nervous system needs an internal critic that is able to predict the expected reward for arbitrary stimuli. We show that, with a critic, reward-modulated spike-timing-dependent plasticity is capable of learning motor trajectories with a temporal resolution of tens of milliseconds. The relation to temporal difference learning, the relevance of block-based learning paradigms, and the limitations of learning with a critic are discussed.

Bayesian updating of a unified soil compression model

Some amount of differential settlement occurs even in the most uniform soil deposit, but it is extremely difficult to estimate because of the natural heterogeneity of the soil. The compression response of the soil and its variability must be characterised in order to estimate the probability of the differential settlement exceeding a certain threshold value. The work presented in this paper introduces a probabilistic framework to address this issue in a rigorous manner, while preserving the format of a typical geotechnical settlement analysis. In order to avoid dealing with different approaches for each category of soil, a simplified unified compression model is used to characterise the nonlinear compression behavior of soils of varying gradation through a single constitutive law. The Bayesian updating rule is used to incorporate information from three different laboratory datasets in the computation of the statistics (estimates of the means and covariance matrix) of the compression model parameters, as well as of the uncertainty inherent in the model.

Fluid-structure interaction simulations of a tension-cone inflatable aerodynamic decelerator

A series of fluid-structure interaction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluid-structure interaction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure. Copyright © 2013 by Christopher Porter, R. Mark Rennie, Eric J. Jumper.