A Real-Time Pressure Wave Model for Predicting Engine Knock

Author(s):  
Ruixue C. Li ◽  
Guoming G. Zhu

Abstract This paper proposes a control-oriented pressure wave model, utilizing outputs of a reaction-based two-zone engine combustion model developed earlier, to accurately predict the key knock characteristics. The model can be used for model-based knock prediction and control. An in-cylinder pressure wave model of oscillation magnitude decay is proposed and simplified to describe pressure oscillations due to knock combustion, and the boundary and initial conditions of the pressure wave model at knock onset are provided by the two-zone reaction-based combustion model. The proposed pressure wave model is calibrated using experimental data, and the chemical kinetic-based Arrhenius integral (ARI) and maximum amplitude of pressure oscillations (MAPO) are used as the evaluation criteria for predicting knock onset and intensity, and the knock frequency is studied with the fast Fourier transform (FFT). The calibrated model is validated for predicting knock onset timing, knock intensity and frequency. Simulation results are compared with the experimental ones to demonstrate the capability of predicting engine knock characteristics by the proposed model.

2019 ◽  
pp. 146808741986916 ◽  
Author(s):  
Ruixue C Li ◽  
Guoming G Zhu

This article develops a control-oriented real-time zero-dimensional knock pressure wave model for spark-ignited engines based on the reaction-based two-zone combustion model developed earlier. This knock pressure wave model is capable of predicting the in-cylinder pressure oscillations under knocking combustion in real-time and can be used for the model-based knock prediction and control. A pressure wave equation including the knock deadening behavior is proposed, simplified, and used to calculate the pressure perturbations generated by knock combustion. The boundary and initial conditions at knock onset are analyzed and the analytic solution of the pressure wave equation is obtained. The model is calibrated and validated over two engine operating conditions at knock limit. The chemical kinetic-based Arrhenius integral and the maximum amplitude of pressure oscillations are used as the evaluation methods for knock onset and intensity prediction, and the knock frequency is studied using a fast Fourier transform of the filtered in-cylinder pressure oscillations. The simulation results show that the proposed knock pressure wave model is able to predict the knock onset timing, frequency, and intensity accurately under two engine operating conditions. Furthermore, the knock characteristics associated with gas mixture properties at intake valve closing is analyzed based on the experimental data, and their effect to knock cycle-to-cycle variation is also studied for the proposed model.


Author(s):  
A. Ahmedi ◽  
O. Stenla˚a˚s ◽  
B. Sunde´n ◽  
R. Egnell ◽  
F. Mauss

Autoignition in SI engines is an abnormal combustion mode and may lead to engine knock in SI engines. Knock may cause damage and it is a source of noise in engines. It limits the compression ratio of the engine and a low compression ratio means low fuel conversion efficiency of the engine. In this paper a multi zone model based on an existing two zone model Hajireza et al., [1 and 12] and Stenla˚a˚s et al., [30] is developed and validated against the experimental results. The validation is done by using the same detailed chemical mechanism consisting of 141 species and about 1405 reactions under the same conditions. The model is a zero dimensional model capable of simulating a full engine cycle. The two zone combustion model consists of a burned and an unburned zone, separated by a thin adiabatic flame front. The multi zone model differs in the handling of the burned gas. In the multi zone case a number of burned zones are present. The number of zones is decided by the temperature difference between the flame front and the last generated burned zone. The detailed chemical mechanism is taken into account in each zone, while the propagating flame front is calculated from the Wiebe function. Each zone is assumed to be a homogeneous mixture with a uniform temperature, mole and mass fractions of species. The spatial variation of the pressure is neglected, i.e., it is assumed to be the same in the whole combustion chamber at every instant of time. Autoignition is handled by the chemical kinetic model. As the unburned zone is assumed homogeneous the effect of auto ignition is a single pressure peak. The model is not designed to predict the pressure oscillations seen in engine knock.


Author(s):  
A. Ahmedi ◽  
F. Mauss ◽  
B. Sunde´n

Constant volume combustion is studied, using a zero-dimensional model, which is a wide-ranging chemical kinetic simulation that allows a closed system of gases to be described on the basis of a set of initial conditions. The model provides an engine- or reactor-like environment in which the engine simulations allow for a variable system volume and heat transfer both to and from the system. The combustion chamber is divided into two zones as burned and unburned ones, which are separated by a thin adiabatic flame front in the combustion model used in this work. A detailed chemical mechanism is applied in each zone to calculate the temperature and pressure history. Equilibrium assumptions have been adopted for the modeling of the thermal ionization, in which Saha’s equation was derived for singly ionized molecules. The investigation is focused on the thermal ionization and electron attachment of 13 chemical species by solving a set of 6 chemical reactions dynamically, the equilibrium calculation using Saha’s equation is performed in a post process, using the temperature and pressure history from the previous model. The experiments that were used for the validation of this model were performed in constant-volume bomb. The outputs generated by the model are temperature profiles, species concentration profiles, ionization degree and an electron density for each zone. The model also predicts the pressure cycle and the ion current. The results from the simulation show good agreement with the experimental measurements and literature data.


Aerospace ◽  
2021 ◽  
Vol 8 (11) ◽  
pp. 341
Author(s):  
Evgeny Strokach ◽  
Victor Zhukov ◽  
Igor Borovik ◽  
Andrej Sternin ◽  
Oscar J. Haidn

In this study, a single injector methane-oxygen rocket combustor is numerically studied. The simulations included in this study are based on the hardware and experimental data from the Technical University of Munich. The focus is on the recently developed generalized k–ω turbulence model (GEKO) and the effect of its adjustable coefficients on the pressure and on wall heat flux profiles, which are compared with the experimental data. It was found that the coefficients of ‘jet’, ‘near-wall’, and ‘mixing’ have a major impact, whereas the opposite can be deduced about the ‘separation’ parameter Csep, which highly influences the pressure and wall heat flux distributions due to the changes in the eddy-viscosity field. The simulation results are compared with the standard k–ε model, displaying a qualitatively and quantitatively similar behavior to the GEKO model at a Csep equal to unity. The default GEKO model shows a stable performance for three oxidizer-to-fuel ratios, enhancing the reliability of its use. The simulations are conducted using two chemical kinetic mechanisms: Zhukov and Kong and the more detailed RAMEC. The influence of the combustion model is of the same order as the influence of the turbulence model. In general, the numerical results present a good or satisfactory agreement with the experiment, and both GEKO at Csep = 1 or the standard k–ε model can be recommended for usage in the CFD simulations of rocket combustion chambers, as well as the Zhukov–Kong mechanism in conjunction with the flamelet approach.


Author(s):  
Peter Düben ◽  
Nils Wedi ◽  
Sami Saarinen ◽  
Christian Zeman

<p>Global simulations with 1.45 km grid-spacing are presented that were performed with the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). Simulations are uncoupled (without ocean, sea-ice or wave model), using 62 or 137 vertical levels and the full complexity of weather forecast simulations including recent date initial conditions, real-world topography, and state-of-the-art physical parametrizations and diabatic forcing including shallow convection, turbulent diffusion, radiation and five categories for the water substance (vapour, liquid, ice, rain, snow). Simulations are evaluated with regard to computational efficiency and model fidelity. Scaling results are presented that were performed on the fastest supercomputer in Europe - Piz Daint (Top 500, Nov 2018). Important choices for the model configuration at this unprecedented resolution for the IFS are discussed such as the use of hydrostatic and non-hydrostatic equations or the time resolution of physical phenomena which is defined by the length of the time step. </p><p>Our simulations indicate that the IFS model — based on spectral transforms with a semi-implicit, semi-Lagrangian time-stepping scheme in contrast to more local discretization techniques — can provide a meaningful baseline reference for O(1) km global simulations.</p>


2020 ◽  
Author(s):  
Shrabanti Roy ◽  
Omid Askari

Abstract Biofuels are considered as an alternative source of energy which can decrease the growing consumption of fossil fuel, hence decreasing pollution. Anisole (methoxybenzene) is a potential source of biofuel produced from cellulose base compounds. It is mostly available as a surrogate of phenolic rich compound. Because of the attractive properties of this fuel in combustion, it is important to do detail kinetic study on oxidation of anisole. In this study a detail chemical mechanism is developed to capture the chemical kinetics of anisole oxidation. The mechanism is developed using an automatic reaction mechanism generator (RMG). To generate the mechanism, RMG uses some known set of species and initial conditions such as temperature, pressure, and mole fractions. Proper thermodynamic and reaction library is used to capture the aromaticity of anisole. The generated mechanism has 340 species and 2532 reactions. Laminar burning speed (LBS) calculated through constant volume combustion chamber (CVCC) at temperature ranges from 460–550 K, pressure of 2–3 atm and equivalence ratio of 0.8–1.4 is used to validate the generated mechanism. Some deviation with experimental result is observed with the newly generated mechanism. Important reaction responsible for LBS calculation, is selected through sensitivity analysis. Rate coefficient of sensitive reactions are collected from literature to modify and improve the mechanism with experimental result. The generated mechanism is further validated with available ignition delay time (IDT) results ranging from 10–20 atm pressure, 0.5–1 equivalence ratio and 870–1600 K temperature. A good agreement of results is observed at different operating ranges. Oxidation of anisole at stoichiometric condition and atmospheric pressure in jet stirred reactor is also used to compare the species concentration of the mechanism. This newly generated mechanism is considered as a good addition for further study of anisole kinetics.


2018 ◽  
Vol 141 (1) ◽  
Author(s):  
Mingzhang Pan ◽  
Haiqiao Wei ◽  
Dengquan Feng

Exhaust gas recirculation (EGR) has gained prominence as a significant method to control port fuel injection engine knock caused by high compression ratio and high intake pressure (IP). In this paper, the effect of EGR on knock intensity was investigated under various conditions which included different compression ratios (9:1, 10:1, 11:1), IPs (1.0 bar, 1.2 bar, 1.4 bar) and intake temperatures (ITs, 20 °C, 40 °C, 60 °C). The torque output being a crucial variant was also considered. The results showed that EGR effectively reduced the maximum amplitude of pressure oscillations (MAPO) and knock intensity factor (KI20). The effect of EGR on knock resistance was more significant at higher compression ratio, IP, and IT. The output torque of the engine reached a peak value with a suitable EGR ratio which also controlled the intensity of knock under different conditions.


2017 ◽  
Vol 821 ◽  
pp. 85-116 ◽  
Author(s):  
Hong-Yueh Lo ◽  
Philip L.-F. Liu

This paper presents a suite of analytical solutions, for both the free-surface elevation and the flow velocity, for landslide-generated water waves. The one-dimensional (horizontal, 1DH) analytical solutions for water waves generated by a solid landslide moving at a constant speed in constant water depth were obtained for the linear and weakly dispersive wave model as well as the linear and fully dispersive wave model. The area enclosed by the landslide was shown to have stronger lasting effects on the generated water waves than the exact landslide shape. In addition, the resonance solution based on the fully dispersive wave model was examined, and the growth rate was derived. For the 1DH linear shallow water equations (LSWEs) on a constant slope, a closed-form analytical solution, which could serve as a useful benchmark for numerical models, was found for a special landslide forcing function. For the two-dimensional (horizontal, 2DH) LSWEs on a plane beach, we rederived the solutions using the quiescent water initial conditions. The difference between the initial conditions used in the new solutions and those used in previous studies was found to have a permanent effect on the generated waves. We further noted that convergence rate of the 2DH LSWE analytical solutions varies greatly, and advised that case-by-case convergence tests be conducted whenever the modal analytical solutions are numerically evaluated using a finite number of modes.


2008 ◽  
Vol 21 (3) ◽  
pp. 576-583 ◽  
Author(s):  
David Ferreira ◽  
Claude Frankignoul

Abstract The transient atmospheric response to interactive SST anomalies in the midlatitudes is investigated using a three-layer QG model coupled in perpetual winter conditions to a slab oceanic mixed layer in the North Atlantic. The SST anomalies are diagnosed from a coupled run and prescribed as initial conditions, but are free to evolve. The initial evolution of the atmospheric response is similar to that obtained with a prescribed SST anomaly, starting as a quasi-linear baroclinic and then quickly evolving into a growing equivalent barotropic one. Because of the heat flux damping, the SST anomaly amplitude slowly decreases, albeit with little change in pattern. Correspondingly, the atmospheric response only increases until it reaches a maximum amplitude after about 1–3.5 months, depending on the SST anomaly considered. The response is similar to that at equilibrium in the fixed SST case, but it is 1.5–2 times smaller, and then slowly decays away.


Sign in / Sign up

Export Citation Format

Share Document