-
Bayesian inference of physics-based models of acoustically-forced laminar premixed conical flames
Authors:
Alessandro Giannotta,
Matthew Yoko,
Stefania Cherubini,
Pietro De Palma,
Matthew P. Juniper
Abstract:
We perform twenty experiments on an acoustically-forced laminar premixed Bunsen flame and assimilate high-speed footage of the natural emission into a physics-based model containing seven parameters. The experimental rig is a ducted Bunsen flame supplied by a mixture of methane and ethylene. A high-speed camera captures the natural emission of the flame, from which we extract the position of the f…
▽ More
We perform twenty experiments on an acoustically-forced laminar premixed Bunsen flame and assimilate high-speed footage of the natural emission into a physics-based model containing seven parameters. The experimental rig is a ducted Bunsen flame supplied by a mixture of methane and ethylene. A high-speed camera captures the natural emission of the flame, from which we extract the position of the flame front. We use Bayesian inference to combine this experimental data with our prior knowledge of this flame's behaviour. This prior knowledge is expressed through (i) a model of the kinematics of a flame front moving through a model of the perturbed velocity field, and (ii) a priori estimates of the parameters of the above model with quantified uncertainties. We find the most probable a posteriori model parameters using Bayesian parameter inference, and quantify their uncertainties using Laplace's method combined with first-order adjoint methods. This is substantially cheaper than other common Bayesian inference frameworks, such as Markov Chain Monte Carlo. This process results in a quantitatively-accurate physics-based reduced-order model of the acoustically forced Bunsen flame for injection velocities ranging from 1.75 m/s to 2.99 m/s and equivalence ratio values ranging from 1.26 to 1.47, using seven parameters. We use this model to evaluate the heat release rate between experimental snapshots, to extrapolate to different experimental conditions, and to calculate the flame transfer function and its uncertainty for all the flames. Since the proposed model relies on only seven parameters, it can be trained with little data and successfully extrapolates beyond the training dataset. Matlab code is provided so that the reader can apply it to assimilate further flame images into the model.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Linear and nonlinear optimal growth mechanisms for generating turbulent bands
Authors:
E. Parente,
J-Ch. Robinet,
P. De Palma,
S. Cherubini
Abstract:
Linear and nonlinear energy optimizations in a tilted domain are used to unveil the main mechanisms allowing the creation of a turbulent band in a channel flow. Linear optimization predicts an optimal growth for streamwse and spanwise wavenumbers $k_x = 1.2$, $k_z = - 1.75$, corresponding to the peak values of the premultiplied energy spectra of direct numerical simulations. At target time, the li…
▽ More
Linear and nonlinear energy optimizations in a tilted domain are used to unveil the main mechanisms allowing the creation of a turbulent band in a channel flow. Linear optimization predicts an optimal growth for streamwse and spanwise wavenumbers $k_x = 1.2$, $k_z = - 1.75$, corresponding to the peak values of the premultiplied energy spectra of direct numerical simulations. At target time, the linear optimal perturbation is composed by oblique streaks, which, for a sufficiently large initial energy, induce turbulence in the whole domain, due to the lack of spatial localization. When localization is achieved by adding nonlinear effects to the optimization, or by artificially confining the linear optimal to a localized region in the spanwise direction, a large-scale flow is created, which leads to the generation of a localised turbulent band. These results suggest that inducing transition towards turbulent bands in a tilted domain, two main elements are needed: a linear energy growth mechanism such as the lift-up for generating large-amplitude flow structures which produce inflection points; large-scale vortices ensuring spatial localisation. Remarkably, these two elements alone are able to generate an isolated turbulent band also in a large, non-tilted domain.
△ Less
Submitted 21 July, 2021;
originally announced July 2021.
-
Minimal seeds for turbulent bands
Authors:
E. Parente,
J-Ch. Robinet,
P. De Palma,
S. Cherubini
Abstract:
In this work, nonlinear variational optimization is used for obtaining minimal seeds for the formation of turbulent bands in channel flow. Using nonlinear optimization together with energy bisection, we have found that the minimal energy threshold for obtaining spatially-patterned turbulence scales with $Re^{-8.5}$ for $Re>1000$. The minimal seed is constituted by a localized spot-like structure s…
▽ More
In this work, nonlinear variational optimization is used for obtaining minimal seeds for the formation of turbulent bands in channel flow. Using nonlinear optimization together with energy bisection, we have found that the minimal energy threshold for obtaining spatially-patterned turbulence scales with $Re^{-8.5}$ for $Re>1000$. The minimal seed is constituted by a localized spot-like structure surrounded by a low-amplitude large-scale quadrupolar structure filling the whole domain. This minimal-energy perturbation of the laminar flow has dominant wavenumbers equal to $0.15$ and $4$ in the streamwise and spanwise directions, respectively, and is characterized by a more marked spatial localization when the Reynolds number increases. At $Re \lesssim 1200$, the minimal seed evolves in time creating an isolated oblique band. Whereas, for $Re\gtrsim 1200$, an almost spanwise-symmetric evolution is observed, giving rise to two distinct bands. A similar evolution is found also at low $Re$ for non-minimal optimal perturbations. This highlights two different mechanisms of formation of turbulent bands in channel flow, depending on the Reynolds number and initial energy of the perturbation. The selection of one of these two mechanisms appears to be affected by the probability of decay of the newly-created stripe, which increases with time, but decreases with the Reynolds number.
△ Less
Submitted 21 July, 2021;
originally announced July 2021.
-
Analysing the resilience of the European commodity production system with PyResPro, the Python Production Resilience package
Authors:
Matteo Zampieri,
Andrea Toreti,
Andrej Ceglar,
Pierluca De Palma,
Thomas Chatzopoulos
Abstract:
This paper presents a Python object-oriented software and code to compute the annual production resilience indicator. The annual production resilience indicator can be applied to different anthropic and natural systems such as agricultural production, natural vegetation and water resources. Here, we show an example of resilience analysis of the economic values of the agricultural production in Eur…
▽ More
This paper presents a Python object-oriented software and code to compute the annual production resilience indicator. The annual production resilience indicator can be applied to different anthropic and natural systems such as agricultural production, natural vegetation and water resources. Here, we show an example of resilience analysis of the economic values of the agricultural production in Europe. The analysis is conducted for individual time-series in order to estimate the resilience of a single commodity and to groups of time-series in order to estimate the overall resilience of diversified production systems composed of different crops and/or different countries. The proposed software is powerful and easy to use with publicly available datasets such as the one used in this study.
△ Less
Submitted 17 June, 2020; v1 submitted 16 June, 2020;
originally announced June 2020.
-
An efficient flamelet progress-variable method for modeling non-premixed flames in weak electric fields
Authors:
M. Di Renzo,
P. De Palma,
M. D. de Tullio,
G. Pascazio
Abstract:
Combustion stabilization and enhancement of the flammability limits are mandatory objectives to improve nowadays combustion chambers. At this purpose, the use of an electric field in the flame region provides a solution which is, at the same time, easy to implement and effective to modify the flame structure. The present work describes an efficient flamelet progress-variable approach developed to…
▽ More
Combustion stabilization and enhancement of the flammability limits are mandatory objectives to improve nowadays combustion chambers. At this purpose, the use of an electric field in the flame region provides a solution which is, at the same time, easy to implement and effective to modify the flame structure. The present work describes an efficient flamelet progress-variable approach developed to model the fluid dynamics of flames immersed in an electric field. The main feature of this model is that it can use complex ionization mechanisms without increasing the computational cost of the simulation. The model is based on the assumption that the combustion process is not directly influenced by the electric field and has been tested using two chemi-ionization mechanisms of different complexity in order to examine its behavior with and without the presence of heavy anions in the mixture. Using a one- and two-dimensional numerical test cases, the present approach has been able to reproduce all the major aspects encountered when a flame is subject to an imposed electric field and the main effects of the different chemical mechanisms. Moreover, the proposed model is shown to produce a large reduction in the computational cost, being able to shorten the time needed to perform a simulation up to 40 times.
△ Less
Submitted 11 August, 2017;
originally announced August 2017.
-
LES of the Sandia Flame D Using an FPV Combustion Model
Authors:
Mario Di Renzo,
Alessandro Coclite,
Marco Donato de Tullio,
Pietro De Palma,
Giuseppe Pascazio
Abstract:
The simulation of turbulent combustion phenomena is still an open problem in modern fluid dynamics. Considering the economical importance of hydrocarbon combustion in energy production processes, it is evident the need of an accurate tool with a relatively low computational cost for the prediction of this kind of reacting flows. In the present work, a comparative study is carried out among large e…
▽ More
The simulation of turbulent combustion phenomena is still an open problem in modern fluid dynamics. Considering the economical importance of hydrocarbon combustion in energy production processes, it is evident the need of an accurate tool with a relatively low computational cost for the prediction of this kind of reacting flows. In the present work, a comparative study is carried out among large eddy simulations, performed with various grid resolutions, a Reynolds averaged Navier-Stokes simulation, and experimental data concerning the well-known Sandia D flame test case. In all the simulations, a flamelet progress variable model has been employed using various hypotheses for the joint probability density function closure. The filtered approach proved to be more accurate than the averaged one, even for the coarser grid used in this work. In fact both approaches have shown poorly accurate predictions in the first part of the combustion chamber, but only by the large eddy simulation one is capable to recover the inlet discrepancies with respect to the experimental data going along the streamwise direction.
△ Less
Submitted 1 February, 2016;
originally announced February 2016.
-
Computing supersonic non-premixed turbulent combustion by an SMLD flamelet progress variable model
Authors:
A. Coclite,
L. Cutrone,
M. Gurtner,
P. De Palma,
O. J. Haidnd,
G. Pascazio
Abstract:
This paper describes the numerical simulation of the NASA Langley Research Center supersonic H2 -Air combustion chamber performed using two approaches to model the presumed probability density function (PDF) in the flamelet progress variable (FPV) framework. The first one is a standard FPV model, built presuming the functional shape of the PDFs of the mixture fraction, Z, and of the progress param…
▽ More
This paper describes the numerical simulation of the NASA Langley Research Center supersonic H2 -Air combustion chamber performed using two approaches to model the presumed probability density function (PDF) in the flamelet progress variable (FPV) framework. The first one is a standard FPV model, built presuming the functional shape of the PDFs of the mixture fraction, Z, and of the progress parameter, Λ. In order to enhance the prediction capabilities of such a model in high-speed reacting flows, a second approach is proposed employing the statistically most likely distribution (SMLD) techcnique to presume the joint PDF of Z and Λ, without any assumption about their behaviour. The standard and FPV-SMLD models have been developed using the low Mach number assumption. In both cases, the temperature is evaluated by solving the total-energy conservation equation, providing a more suitable approach for the simulation of supersonic combustion. By comparison with experimental data, the proposed SMLD model is shown to provide a clear improvement with respect to the standard FPV model, especially in the auto-ignition and stabilization regions of the flame.
△ Less
Submitted 23 September, 2015;
originally announced September 2015.
-
An SMLD joint PDF model for turbulent non-premixed combustion using the flamelet progress-variable approach
Authors:
Alessandro Coclite,
Giuseppe Pascazio,
Pietro De Palma,
Luigi Cutrone,
Matthias Ihme
Abstract:
This paper provides an improved flamelet/progress variable (FPV) model for the simulation of turbulent combustion, employing the statistically most likely distribution (SMLD) approach for the joint probability density function (PDF) of the mixture fraction, Z, and of the progress parameter, Λ . Steady-state FPV models are built presuming the func- tional shape of the joint PDF of Z and Λ in order…
▽ More
This paper provides an improved flamelet/progress variable (FPV) model for the simulation of turbulent combustion, employing the statistically most likely distribution (SMLD) approach for the joint probability density function (PDF) of the mixture fraction, Z, and of the progress parameter, Λ . Steady-state FPV models are built presuming the func- tional shape of the joint PDF of Z and Λ in order to evaluate Favre-averages of thermody- namic quantities. The mixture fraction is widely assumed to behave as a passive scalar with a mono-modal behaviour modelled by a \b{eta} -distribution. Moreover, under the hypothesis that Z and Λ are statistically independent, the joint PDF coincides with the product of the two marginal PDFs. In this work we discuss these two constitutive hypotheses. The proposed model evaluates the most probable joint distribution of Z and Λ , relaxing some crucial as- sumption on their statistical behaviour. This provides a more general model in the context of FPV approach and an effective tool to verify the adequateness of widely used hypotheses. The model is validated versus experimental data of well-known test cases, namely, the San- dia flames. The results are also compared with those obtained by the standard FPV approach, analysing the role of the PDF functional form on turbulent combustion simulations.
△ Less
Submitted 20 April, 2015;
originally announced April 2015.
-
Numerical investigation of high-pressure combustion in rocket engines using Flamelet/Progress-variable models
Authors:
A. Coclite,
L. Cutrone,
P. De Palma,
G. Pascazio
Abstract:
The present paper deals with the numerical study of high pressure LOx/H2 or LOx/hydrocarbon combustion for propulsion systems. The present research effort is driven by the continued interest in achieving low cost, reliable access to space and more recently, by the renewed interest in hypersonic transportation systems capable of reducing time-to-destination. Moreover, combustion at high pressure ha…
▽ More
The present paper deals with the numerical study of high pressure LOx/H2 or LOx/hydrocarbon combustion for propulsion systems. The present research effort is driven by the continued interest in achieving low cost, reliable access to space and more recently, by the renewed interest in hypersonic transportation systems capable of reducing time-to-destination. Moreover, combustion at high pressure has been assumed as a key issue to achieve better propulsive performance and lower environmental impact, as long as the replacement of hydrogen with a hydrocarbon, to reduce the costs related to ground operations and increase flexibility. The current work provides a model for the numerical simulation of high- pressure turbulent combustion employing detailed chemistry description, embedded in a RANS equations solver with a Low Reynolds number k-omega turbulence model. The model used to study such a combustion phenomenon is an extension of the standard flamelet-progress-variable (FPV) turbulent combustion model combined with a Reynolds Averaged Navier-Stokes equation Solver (RANS). In the FPV model, all of the thermo-chemical quantities are evaluated by evolving the mixture fraction Z and a progress variable C. When using a turbulence model in conjunction with FPV model, a probability density function (PDF) is required to evaluate statistical averages of chemical quantities. The choice of such PDF must be a compromise between computational costs and accuracy level. State- of-the-art FPV models are built presuming the functional shape of the joint PDF of Z and C in order to evaluate Favre-averages of thermodynamic quantities. The model here proposed evaluates the most probable joint distribution of Z and C without any assumption on their behavior.
△ Less
Submitted 2 December, 2014;
originally announced February 2015.
-
An extended SMLD approach for presumed probability density function in flamelet combustion model
Authors:
Alessandro Coclite,
Giuseppe Pascazio,
Pietro De Palma,
Luigi Cutrone
Abstract:
This paper provides an extension of the standard flamelet progress variable (FPV) approach for turbulent combustion, applying the statistically most likely distribution (SMLD) framework to the joint PDF of the mixture fraction, Z, and the progress variable, C. In this way one does not need to make any assumption about the statistical correlation between Z and C and about the behaviour of the mixtu…
▽ More
This paper provides an extension of the standard flamelet progress variable (FPV) approach for turbulent combustion, applying the statistically most likely distribution (SMLD) framework to the joint PDF of the mixture fraction, Z, and the progress variable, C. In this way one does not need to make any assumption about the statistical correlation between Z and C and about the behaviour of the mixture fraction, as required in previous FPV models. In fact, for state-of-the-art models, with the assumption of very-fast-chemistry,Z is widely accepted to behave as a passive scalar characterized by a $β$-distribution function. Instead, the model proposed here, evaluates the most probable joint distribution of Z and C without any assumption on their behaviour and provides an effective tool to verify the adequateness of widely used hypotheses, such as their statistical independence. The model is validated versus three well-known test cases, namely, the Sandia flames. The results are compared with those obtained by the standard FPV approach, analysing the role of the PDF functional form on turbulent combustion simulations.
△ Less
Submitted 20 June, 2013;
originally announced June 2013.
-
The role of presumed probability density function in the simulation of non premixed turbulent combustion
Authors:
Alessandro Coclite,
Giuseppe Pascazio,
Pietro De Palma,
Luigi Cutrone
Abstract:
Flamelet Progress Variable (FPV) combustion models allow the evaluation of all thermo chemical quantities in a reacting flow by computing only the mixture fraction Z and a progress variable C. When using such a method to predict a turbulent combustion in conjunction with a turbulence model, a probability density function (PDF) is required to evaluate statistical averages (e.g., Favre average) of c…
▽ More
Flamelet Progress Variable (FPV) combustion models allow the evaluation of all thermo chemical quantities in a reacting flow by computing only the mixture fraction Z and a progress variable C. When using such a method to predict a turbulent combustion in conjunction with a turbulence model, a probability density function (PDF) is required to evaluate statistical averages (e.g., Favre average) of chemical quantities. The choice of the PDF is a compromise between computational costs and accuracy level. The aim of this paper is to investigate the influence of the PDF choice and its modeling aspects in the simulation of non premixed turbulent combustion. Three different models are considered: the standard one, based on the choice of a beta distribution for Z and a Dirac distribution for C; a model employing a beta distribution for both Z and C; a third model obtained using a beta distribution for Z and the statistical most likely distribution (SMLD) for C. The standard model, although widely used, doesn't take into account the interaction between turbulence and chemical kinetics as well as the dependence of the progress variable not only on its mean but also on its variance. The SMLD approach establishes a systematic framework to incorporate informations from an arbitrary number of moments, thus providing an improvement over conventionally employed presumed PDF closure models. The rational behind the choice of the three PDFs is described in some details and the prediction capability of the corresponding models is tested versus well known test cases, namely, the Sandia flames, and a test case for supersonic combustion provided by Cheng.
△ Less
Submitted 21 January, 2014; v1 submitted 25 May, 2013;
originally announced May 2013.