Advanced flamelet modelling of turbulent nonpremixed and partialy premixed combustion

Current work focuses on the development and performance evaluation of advanced flamelet models for turbulent non-premixed and partially premixed combustion in RANS and large eddy simulation (LES) based modelling. A RANS based combustion modelling strategy which has the ability to capture the detailed structure of turbulent non-premixed flames, including the pollutant NO, and account for the effects of radiation heat loss and transient evolution of NO, has been developed and incorporated into the in-house RANS code. The strategy employs an 'enthalpy-defect' based non-adiabatic flamelet model in conjunction with steady or unsteady nonadiabatic flamelets based NO submodels. The performanceo f the non-adiabaticm odel and its NO submodelsh asb eena ssessed against experimental measurements and steady flamelet model predictions for turbulent CH4/H2 bluff-body stabilized and CH4-air piloted jet flames. Appreciable improvements in the mean thermal structure predictions have been observed in the piloted jet flames by consideration of radiation heat loss through the non-adiabatic flamelet model. Since transient effects were weaker in the piloted jet flame, both unsteady and steady non-adiabatic NO submodels provided similar level of improvement in the pollutant NO predictions in comparison to their adiabatic counterpartsT. ransiente ffectsw ere, however,d ominanti n the bluff-body flame. The unsteady non-adiabatic NO submodel provided excellent agreement with measured NO distribution in comparison to the appreciably overpredicted distribution by its steadyc ounterpart.T he strategyo f non-adiabaticf lamelet model in conjunctionw ith unsteady non-adiabatic NO submodel seems to provide an accurate and robust alternative to the conventional strategy of steady flamelet model with steady NO submodel. While addressing the limitations of steady flamelet model in regard to radiation and slow chemistry of NO is one objective of this research, extending the applicability of the model to partially premixed combustion has been pursued as the second objective. Flamelet/progress variable (FPV) approach based combustion models, which have the potential to describe both non-premixed and partially premixed combustion, have been incorporated in the in-house RANS and LES codes. Based on the form of the PDF for reaction progress variable, two different formulations, FPV-8 function model and FPV-P function model, have been derived. (Continues...)

Advanced flamelet modelling of turbulent non-premixed and partially premixed combustion

Current work focuses on the development and performance evaluation of advanced flamelet models for turbulent non-premixed and partially premixed combustion in RANS and large eddy simulation (LES) based modelling. A RANS-based combustion modelling strategy which has the ability to capture the detailed structure of turbulent non-premixed flames, including the pollutant NO, and account for the effects of radiation heat loss and transient evolution of NO, has been developed and incorporated into the in-house RANS code. The strategy employs an 'enthalpy defect'-based non-adiabatic flamelet model in conjunction with steady or unsteady nonadiabatic flamelets based NO submodels. [Continues.

A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations

summary:A new error correction method for the stationary Navier-Stokes equations based on two local Gauss integrations is presented. Applying the orthogonal projection technique, we introduce two local Gauss integrations as a stabilizing term in the error correction method, and derive a new error correction method. In both the coarse solution computation step and the error computation step, a locally stabilizing term based on two local Gauss integrations is introduced. The stability and convergence of the new error correction algorithm are established. Numerical examples are also presented to verify the theoretical analysis and demonstrate the efficiency of the proposed method

LES modelling of non-premixed and partially premixed turbulent flames

A large eddy simulation (LES) model has been developed and validated for turbulent non-premixed and partially premixed combustion systems. LES based combustion modelling strategy has the ability to capture the detailed structure of turbulent flames and account for the effects of radiation heat loss. Effects of radiation heat loss is modelled by employing an enthalpy-defect based non-adiabatic flamelet model (NAFM) in conjunction with a steady non-adiabatic flamelet approach. The steady laminar flamelet model (SLFM) is used with multiple flamelet solutions through the development of pre-integrated look up tables. The performance of the non-adiabatic model is assessed against experimental measurements of turbulent CH4/H2 bluff-body stabilized and swirl stabilized jet flames carried out by the University of Sydney combustion group. Significant enhancements in the predictions of mean thermal structure have been observed with both bluff body and swirl stabilized flames by the consideration of radiation heat loss through the non-adiabatic flamelet model. In particular, mass fractions of product species like CO2 and H2O have been improved with the consideration of radiation heat loss. From the Sydney University data the HM3e flame was also investigated with SLFM using multiple flamelet strategy and reasonably fair amount of success has been achieved. In this work, unsteady flamelet/progress variable (UFPV) approach based combustion model which has the potential to describe both non-premixed and partially premixed combustion, has been developed and incorporated in an in-house LES code. The probability density function (PDF) for reaction progress variable and scalar dissipation rate is assumed to follow a delta distribution while mixture fraction takes the shape of a beta PDF. The performance of the developed model in predicting the thermal structure of a partially premixed lifted turbulent jet flame in vitiated co-flow has been evaluated. The UFPV model has been found to successfully predict the flame lift-off, in contrast SLFM results in a false attached flame. The mean lift-off height is however over-predicted by UFPV-δ function model by ~20% for methane based flame and under-predicted by ~50% for hydrogen based flame. The form of the PDF for the reaction progress variable and inclusion of a scalar dissipation rate thus seems to have a strong influence on the predictions of gross characteristics of the flame. Inclusion of scalar dissipation rate in the calculations appears to be successful in predicting the flame extinction and re-ignition phenomena. The beta PDF distribution for the reaction progress variable would be a true prospect for extending the current simulation to predict the flame characteristics to a higher degree.EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Advanced flamelet modelling of turbulent nonpremixed and partialy premixed combustion

Current work focuses on the development and performance evaluation of advanced flamelet models for turbulent non-premixed and partially premixed combustion in RANS and large eddy simulation (LES) based modelling. A RANS based combustion modelling strategy which has the ability to capture the detailed structure of turbulent non-premixed flames, including the pollutant NO, and account for the effects of radiation heat loss and transient evolution of NO, has been developed and incorporated into the in-house RANS code. The strategy employs an 'enthalpy-defect' based non-adiabatic flamelet model in conjunction with steady or unsteady nonadiabatic flamelets based NO submodels. The performanceo f the non-adiabaticm odel and its NO submodelsh asb eena ssessed against experimental measurements and steady flamelet model predictions for turbulent CH4/H2 bluff-body stabilized and CH4-air piloted jet flames. Appreciable improvements in the mean thermal structure predictions have been observed in the piloted jet flames by consideration of radiation heat loss through the non-adiabatic flamelet model. Since transient effects were weaker in the piloted jet flame, both unsteady and steady non-adiabatic NO submodels provided similar level of improvement in the pollutant NO predictions in comparison to their adiabatic counterpartsT. ransiente ffectsw ere, however,d ominanti n the bluff-body flame. The unsteady non-adiabatic NO submodel provided excellent agreement with measured NO distribution in comparison to the appreciably overpredicted distribution by its steadyc ounterpart.T he strategyo f non-adiabaticf lamelet model in conjunctionw ith unsteady non-adiabatic NO submodel seems to provide an accurate and robust alternative to the conventional strategy of steady flamelet model with steady NO submodel. While addressing the limitations of steady flamelet model in regard to radiation and slow chemistry of NO is one objective of this research, extending the applicability of the model to partially premixed combustion has been pursued as the second objective. Flamelet/progress variable (FPV) approach based combustion models, which have the potential to describe both non-premixed and partially premixed combustion, have been incorporated in the in-house RANS and LES codes. Based on the form of the PDF for reaction progress variable, two different formulations, FPV-8 function model and FPV-P function model, have been derived. (Continues...).EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Mathematical Architecture for Models of Fluid Flow Phenomena

This thesis is a study of several high accuracy numerical methods for fluid flow problems and turbulence modeling.First we consider a stabilized finite element method for the Navier-Stokes equations which has second order temporal accuracy. The method requires only the solution of one linear system (arising from an Oseen problem) per time step. We proceed by introducing a family of defect correction methods for the time dependent Navier-Stokes equations, aiming at higher Reynolds' number. The method presented is unconditionally stable, computationally cheap and gives an accurate approximation to the quantities sought. Next, we present a defect correction method with increased time accuracy. The method is applied to the evolutionary transport problem, it is proven to be unconditionally stable, and the desired time accuracy is attained with no extra computational cost. We then turn to the turbulence modeling in coupled Navier-Stokes systems - namely, MagnetoHydroDynamics. Magnetically conducting fluids arise in important applications including plasma physics, geophysics and astronomy. In many of these, turbulent MHD (magnetohydrodynamic) flows are typical. The difficulties of accurately modeling and simulating turbulent flows are magnified many times over in the MHD case. We consider the mathematical properties of a model for the simulation of the large eddies in turbulent viscous, incompressible, electrically conducting flows. We prove existence, uniqueness and convergence of solutions for the simplest closed MHD model. Furthermore, we show that the model preserves the properties of the 3D MHD equations. Lastly, we consider the family of approximate deconvolution models (ADM) for turbulent MHD flows. We prove existence, uniqueness and convergence of solutions, and derive a bound on the modeling error. We verify the physical properties of the models and provide the results of the computational tests

Large-eddy simulation of turbulent wall-pressure fluctuations using the finite element method

​In the present dissertation, turbulent wall-pressure fluctuations are characterized. To capture the turbulent characteristics of the flow, large-eddy simulation is used to resolve the large scale motions of the flow directly. A wall-adapting local eddy-viscosity model is selected to account for the effect of small scale motions. The streamwise/upwind Petrov-Galerkin method is chosen to discretize the computational domain and a second-order backward difference formula is applied for the time integration. Maintaining turbulent flow throughout the simulation domain to properly characterize turbulence is critical in investigating wall-pressure fluctuations. In order to reduce the size of the simulation domain an inflow generation method, a variant of the recycling and rescaling method, is used. In this method, the turbulent velocity profile from a specific plane within the computational domain is recycled and rescaled propriately, and re-introduced at the inlet of the domain at every time step iteration. In the proposed method, the mean velocity profile is fixed at the inlet while the velocity fluctuations are recycled and rescaled to obtain the desired turbulent characteristics. This method is simple and effective and maintains the turbulent flow throughout the simulation domain. The non-reflecting boundary conditions with a sponge layer are applied at the top and exit of the computational domain to remove unwanted reflections from the boundary. In order to examine the present inflow generation method and the ability to capture the wall-pressure fluctuations, numerical results are verified on a flat plate with a zero pressure gradient. The mean velocity profile, the RMS velocity fluctuations, and the friction velocity over time are investigated to show the effectiveness of the present inflow turbulent generation method. Computed wall-pressure fluctuations are evaluated using the time-averaged statistics and the spectra, to show that they are characterized well using the present method

Nodally exact Ritz discretizations of 1D diffusion–absorption and Helmholtz equations by variational FIC and modified equation methods

The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-005-0011-zThis article presents the first application of the Finite Calculus (FIC) in a Ritz-FEM variational framework. FIC provides a steplength parametrization of mesh dimensions, which is used to modify the shape functions. This approach is applied to the FEM discretization of the steady-state, one-dimensional, diffusion–absorption and Helmholtz equations. Parametrized linear shape functions are directly inserted into a FIC functional. The resulting Ritz-FIC equations are symmetric and carry a element-level free parameter coming from the function modification process. Both constant- and variable-coefficient cases are studied. It is shown that the parameter can be used to produce nodally exact solutions for the constant coefficient case. The optimal value is found by matching the finite-order modified differential equation (FOMoDE) of the Ritz-FIC equations with the original field equation. The inclusion of the Ritz-FIC models in the context of templates is examined. This inclusion shows that there is an infinite number of nodally exact models for the constant coefficient case. The ingredients of these methods (FIC, Ritz, MoDE and templates) can be extended to multiple dimensions.Peer ReviewedPostprint (author's final draft