Search results

The purpose of this study is to present and demonstrate a numerical method for solving chemically reacting flows. These are important for energy conversion devices, which rely on chemical reactions as their operational mechanism, with heat generated from the combustion of the fuel, often gases, being converted to work.

The numerical study of such flows requires the set of Navier-Stokes equations to be extended to include multiple species and the chemical reactions between them. The numerical method implemented in this study also accounts for changes in the material properties because of temperature variations and the process to handle steep spatial fronts and stiff source terms without incurring any numerical instabilities. An all-speed numerical framework is used through simple low-dissipation advection upwind splitting (SLAU) convective scheme, and it has been extended in a multi-component species framework on the in-house density-based flow solver. The capability of solving turbulent combustion is also implemented using the Eddy Dissipation Concept (EDC) framework and the recent k-kl turbulence model.

The numerical implementation has been demonstrated for several stiff problems in laminar and turbulent combustion. The laminar combustion results are compared from the corresponding results from the Cantera library, and the turbulent combustion computations are found to be consistent with the experimental results.

This paper has extended the single gas density-based framework to handle multi-component gaseous mixtures. This paper has demonstrated the capability of the numerical framework for solving non-reacting/reacting laminar and turbulent flow problems. The all-speed SLAU convective scheme has been extended in the multi-component species framework, and the turbulent model k-kl is used for turbulent combustion, which has not been done previously. While the former method provides the capability of solving for low-speed flows using the density-based method, the later is a length-scale-based method that includes scale-adaptive simulation characteristics in the turbulence modeling. The SLAU scheme has proven to work well for unsteady flows while the k-kL model works well in non-stationary turbulent flows. As both these flow features are commonly found in industrially important reacting flows, the convection scheme and the turbulence model together will enhance the numerical predictions of such flows.

The purpose of this paper is to ascertain the harm to personnel and equipment caused by an external explosion during natural gas explosion venting. The external explosion characteristics induced by the indoor natural gas explosion are the focal points of the investigation.

Computational fluid dynamics technology was used to investigate the large-scale explosion venting process of natural gas in a 6 × 3 × 2.5 m room, and the characteristics of external explosion under different scaled vent size (K_v = A_v/V2/3, 0.05, 0.08, 0.13, 0.18) were numerically analyzed.

When K_v = 0.08, the length and duration of the explosion fireball are 13.39 and 450 ms, respectively, which significantly expands the degree and range of high-temperature hazards. The suitable flow-field structure causes the external explosion overpressure to be more than twice that indoors, i.e. the natural gas explosion venting overpressure may be considerably more hazardous in an outdoor environment than inside a room. A specific range for the K_v can promote the superposition of outdoor rupture waves and explosion shock waves, thereby creating a new overpressure hazard.

Little attention has been devoted to investigating systematically the external explosion hazards. Based on the numerical simulation and the analysis, the external explosion characteristics induced by the indoor large-scale gas explosion were obtained. The research results are theoretically significant for mitigating the effects of external gas explosions on personnel and equipment.

The Polytechnic Institute of Brooklyn celebrated in 1955 its centenary and as part of the celebrations a Conference on High Speed Aeronautics was organized. In this book are collected the lectures presented at the Conference. The list of lectures is impressive, whilst the lectures are all of great interest, many containing much that is new and fundamental in this ever‐expanding subject. As might be expected, most attention was focused on flight at very high Mach numbers and the aerodynamic and structural problems associated with kinetic heating effects, but that persistent source of problems, transonic flow, also received attention.

Aeronautical research scientists of the National Advisory Committee for Aeronautics have in the last six years indexed 7,070 non‐security classified NACA research reports under 458 different subject classification headings for a total of 18,619 subject entries. They further face the task of subject classifying over 400 similar research reports per year, plus the large number of security classified reports issued annually.

The complete system of equations for a theory of laminar flame equations is presented, taking into account both heat conduction and diffusion, and for the case of an arbitrary number of simultaneous reactions. The eigen value problem determining the flame velocity is formulated. Two examples are given in order to show that explicit analytical expressions for the flame velocity can be obtained, which are in good agreement with the results obtained by numerical integration of the equations. In the first example (hydrazine decomposition) one reaction is considered as global reaction. In the second example (ozone decomposition) a hypothesis is introduced for the concentration of the free radical O, which corresponds to the steady‐state approximation generally used in classical chemical kinetics. In both cases the measured flame velocities are between the flame velocities computed with no diffusion, and with a coefficient estimated by Professor Hirschfelder from the kinetic theory of gases. The approximate explicit formulae are obtained without drastic assumptions and using legitimate approximation methods. The assumption used for the ozone decomposition flame has a bearing on a better understanding of the mechanism of chain reactions in general. The method indicated in the paper gives hope that the more complicated chain reactions such as the combustion of hydrocarbons will also be made accessible to theoretical computation.

Safety improvement and pollutant reduction in many practical combustion systems and especially in aero-gas turbine engines require an adequate understanding of flame ignition and stabilization mechanisms. Improved software and hardware have opened up greater possibilities for translating basic knowledge and the results of experiments into better designs. The present study deals with the large eddy simulation (LES) of an ignition sequence in a conical shaped bluff-body stabilized burner involving a turbulent non-premixed flame. The purpose of this paper is to investigate the impact of spark location on ignition success. Particular attention is paid to the ease of handling of the numerical tool, the computational cost and the accuracy of the results.

The discrete particle ignition kernel (DPIK) model is used to capture the ignition kernel dynamics in its early stage of growth after the breakdown period. The ignition model is coupled with two combustion models based on the mixture fraction-progress variable formulation. An infinitely fast chemistry assumption is first done, and the turbulent fluctuations of the progress variable are captured with a bimodal probability density function (PDF) in the line of the Bray–Moss–Libby (BML) model. Thereafter, a finite rate chemistry assumption is considered through the flamelet-generated manifold (FGM) method. In these two assumptions, the classical beta-PDF is used to model the temporal fluctuations of the mixture fraction in the turbulent flow. To model subgrid scale stresses and residual scalars fluxes, the wall-adapting local eddy (WALE) and the eddy diffusivity models are, respectively, used under the low-Mach number assumption.

Numerical results of velocity and mixing fields, as well as the ignition sequences, are validated through a comparison with their experimental counterparts. It is found that by coupling the DPIK model with each of the two combustion models implemented in a LES-based solver, the ignition event is reasonably predicted with further improvements provided by the finite rate chemistry assumption. Finally, the spark locations most likely to lead to a complete ignition of the burner are found to be around the shear layer delimiting the central recirculation zone, owing to the presence of a mixture within flammability limits.

Some discrepancies are found in the radial profiles of the radial velocity and consequently in those of the mixture fraction, owing to a mismatch of the radial velocity at the inlet section of the computational domain. Also, unlike FGM methods, the BML model predicts the overall ignition earlier than suggested by the experiment; this may be related to the overestimation of the reaction rate, especially in the zones such as flame holder wakes which feature high strain rate due to fuel-air mixing.

This work is adding a contribution for ignition modeling, which is a crucial issue in various combustion systems and especially in aircraft engines. The exclusive use of a commercial computational fluid dynamics (CFD) code widely used by combustion system manufacturers allows a direct application of this simulation approach to other configurations while keeping computing costs at an affordable level.

This study provides a robust and simple way to address some ignition issues in various spark ignition-based engines, namely, the optimization of engines ignition with affordable computational costs. Based on the promising results obtained in the current work, it would be relevant to extend this simulation approach to spray combustion that is required for aircraft engines because of storage volume constraints. From this standpoint, the simulation approach formulated in the present work is useful to engineers interested in optimizing the engines ignition at the design stage.

The internal combustion engine operated on gaseous fuels shows great potential in terms of integration of the renewable and traditional sources for an effective solution for clean energy production challenge. Different fuel mixtures that can be used to power the engine are characterized by various combustion properties, which can affect its overall efficiency. The purpose of this paper is to provide reasonable answer, how the operation condition can change due to different fuel, without enormous cost of prototyping processes using physical models a digital model can be seen as promising technique.

Presented work discusses the application, and extensive description of two commercial codes Ansys Fluent and Forte for modeling stationary engine fueled by compressed natural gas (CNG) and biogas. To check the model accuracy, all carried out numerical results were compared against experimental data collected at in-house test rig of single cylinder four stroke engine. The impacts of tested gaseous fuel on the engine working conditions and emission levels were investigated.

Carried out simulations showed good agreement with experimental data for investigated cases. Application on numerical models give possibility to visualize flame front propagation and pollutant formation for tested fuels. The biogas fuel has shown the impaired early flame phase, which led to longer combustion, lower efficiency, power output, repeatability and in some cases higher HC and carbon monoxide (CO) emissions as a result of combustion during the exhaust stroke. Looking at the CO formation it was observed that it instantly accrue with flame front propagation as a result of methane oxidation, while for NOx formation revers effect was seen.

The application of new approach for modeling combustion process in stationary engines fueled by CNG and alternative biogas fuel has been discussed. The cons and pros of the Forte code in terms of its application for engine prototaping process has been discussed.

The advent of jet and rocket propulsion in our present age has brought the problems of combustion right into the foreground. While much is known, still more theoretical and experimental research is necessary to extend our knowledge of combustion phenomena in combustion chambers and rockets. This impelled the Advisory Group for Aeronautical Research and Development (AGARD) of N.A.T.O. to collect all the latest information and results of research in the field of combustion science by forming a Combustion Panel and by inviting competent experts from Europe and the United States to discuss some of the most acute problems, and thus inspire interest and suggest directions for further research. The first Colloquium was held in December 1952 in Rome, and a further and more systematic one in December 1953 in Cambridge.

This study aims to deal with the large eddy simulation (LES) of an ignition sequence and the resulting steady combustion in a swirl-stabilized liquid-fueled combustor. Particular attention is paid to the ease of handling the numerical tool, the accuracy of the results and the reasonable computational cost involved. The primary aim of the study is to appraise the ability of the newly developed computational fluid dynamics (CFD) methodology to retrieve the spark-based flame kernel initiation, its propagation until the full ignition of the combustion chamber, the flame stabilization and the combustion processes governing the steady combustion regime.

The CFD model consists of an LES-based spray module coupled to a subgrid-scale ignition model to capture the flame kernel initiation and the early stage of the flame kernel growth, and a combustion model based on the mixture fraction-progress variable formulation in the line of the flamelet generated manifold (FGM) method to retrieve the subsequent flame propagation and combustion properties. The LES-spray module is based on an Eulerian-Lagrangian approach and includes a fully two-way coupling at each time step to account for the interactions between the liquid and the gaseous phases. The Wall-Adapting Local Eddy-viscosity (WALE) model is used for the flow field while the eddy diffusivity model is used for the scalar fluxes. The fuel is liquid kerosene, injected in the form of a polydisperse spray of droplets. The spray dynamics are tracked using the Lagrangian procedure, and the phase transition of droplets is calculated using a non-equilibrium evaporation model. The oxidation mechanism of the Jet A-1 surrogate is described through a reduced reaction mechanism derived from a detailed mechanism using a species sensitivity method.

By comparing the numerical results with a set of published data for a swirl-stabilized spray flame, the proposed CFD methodology is found capable of capturing the whole spark-based ignition sequence in a liquid-fueled combustion chamber and the main flame characteristics in the steady combustion regime with reasonable computing costs.

The proposed CFD methodology simulates the whole ignition sequence, namely, the flame kernel initiation, its propagation to fully ignite the combustion chamber, and the global flame stabilization. Due to the lack of experimental ignition data on this liquid-fueled configuration, the ability of the proposed CFD methodology to accurately predict ignition timing was not quantitatively assessed. It would, therefore, be interesting to apply this CFD methodology to other configurations that have experimental ignition data, to quantitatively assess its ability to predict the ignition timing and the flame characteristics during the ignition sequence. Such further investigations will not only provide further validation of the proposed methodology but also will potentially identify its shortfalls for better improvement.

This CFD methodology is developed by customizing a commercial CFD code widely used in the industry. It is, therefore, directly applicable to practical configurations, and provides not only a relatively straightforward approach to predict an ignition sequence in liquid-fueled combustion chambers but also a robust way to predict the flame characteristics in the steady combustion regime as significant improvements are noticed on the prediction of slow species.

The incorporation of the subgrid ignition model paired with a combustion model based on tabulated chemistry allows reducing computational costs involved in the simulation of the ignition phase. The incorporation of the FGM-based tabulated chemistry provides a drastic reduction of computing resources with reasonable accuracy. The CFD methodology is developed using the platform of a commercial CFD code widely used in the industry for relatively straightforward applicability.

This paper aims to introduce a two-dimensional smoothed particle hydrodynamics (SPH) framework for simulating the evaporation and combustion process of fuel droplets.

To solve the gas–liquid two-phase flow problem, a multiphase SPH method capable of handling high density-ratio problems is established. Based on the Fourier heat conduction equation and Fick’s law of diffusion, the SPH discrete equations are derived. To effectively characterize the phase transition problem, inspired by volume of fluid method, the concept of liquid phase mass fraction of the SPH particles is proposed. The one-step global reaction model of n-hexane is used for the vapor combustion.

The evaporation and combustion process of single droplet conforms to the law. The framework works out well when the evaporation of multiple droplets involves coalescence process. Three different kinds of flames are observed in succession in the combustion process of a single droplet at different inflow velocity, which agree well with the results of the experiment.

To the best of the authors’ knowledge, this is the first computational framework that has the capability to simulate evaporation and combustion with SPH method. Based on the particle nature of SPH method, the framework has natural advantages in interface tracking.