Refine
Year of publication
Document Type
Language
- English (28)
Keywords
- Charcoal bed (2)
- Combustion and combustion processes (2)
- Discrete element method (2)
- Euler-Lagrange (2)
- Fluidized bed (2)
- Fuel injection (2)
- Gasification (2)
- Linear eddy model (2)
- ODT (2)
- Pyrolysis (2)
Institute
This paper presents an investigation of primary breakup of planar turbulent liquid jets and breakup properties at the surface of turbulent jets in still air at standard conditions. Numerical simulations are carried out for jet exit Reynolds number 23000 and Weber numbers in the range [102–107].
Due to the limitation of direct numerical simulation (DNS) to moderate Reynolds numbers, a stochastic 1D ansatz based on the one-dimensional turbulence (ODT) model is used to simulate a planar liquid jet with a high lateral resolution. ODT permits an affordable high resolution of interface and single-phase property gradients which are key for understanding the local behavior. ODT is a stochastic model simulating turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. The most relevant mechanisms that influence the primary breakup of liquid jets are found to be represented accurately based on comparisons to experiments and correlations reported in the literature. Building on this finding, future work will focus on the statistics of droplets generated by primary breakup, both to investigate their relationship to breakup mechanisms and to provide input to models of secondary breakup and subsequent spray evolution.
We present a new modelling approach for turbulence modulation by droplets on the subgrid-scale (SGS) level of Large-Eddy-Simulations (LES). Many SGS models exist for the effect of gas phase SGS on the droplet phase, but very few for the mechanisms vice versa on the turbulent intensity of the gas phase. The reasons are a lack of physical understanding and limited computational resources for extensive DNS studies. To address both problems a dimension-reduced and consequently less costly model, namely One-Dimensional- Turbulence (ODT), is used to gather information about this specific flow phenomena. ODT is a stochastic tool simulating turbulent flows along a notional 1D line of sights. For modeling the turbulent advection instantaneous maps are applied to the line which represent the effect of individual eddies on property fields and the dispersed phase. After validating the general test case of a droplet-laden shear flow against DNS data, a concept is presented on how to gather turbulence modulation for several parameter ranges in a data base and how to make them accessible on the flight for LES. The three most significant parameters, the unladen flow Reynolds number, the droplet loading and the particle momentum number, are chosen to construct an efficient data base.
Increasing the injection pressure in gasoline direct injection engines has a substantial potential to reduce emissions while maintaining a high efficiency in spark ignition engines. Present gasoline injectors are operating in the range of 20 MPa to 25 MPa. Now there is an interest in higher fuel injection pressures, for instance, around 40 MPa, 60 MPa and even higher pressures, because of its potential for further emission reduction and fuel efficiency improvements. In order to fully utilize the high-pressure fuel injection technology, a fundamental understanding of gasoline spray characteristics is vital to gain insight into spray behavior under such high injection pressures. The understanding achieved may also be beneficial to improve further model development and facilitate the integration of such advanced injection systems into future gasoline engines. In the present study, a gasoline fuel spray has been investigated over a range of fuel injection pressures from 40 to 150 MPa through a numerical simulation study. The numerical calculations have been performed in a constant volume chamber under non-vaporizing conditions to best match the experimental setup. The numerical model utilized a large-eddy simulation (LES) approach for the gas flow and a standard Lagrangian spray model for the liquid phase. The spray atomization has been modeled using the Kelvin Helmholtz - Rayleigh Taylor (KH-RT) atomization model with a droplet size distribution from the injector assumed to follow a Rosin-Rammler distribution function. Simulation results for the spray liquid penetration length are validated with experimental findings under different fuel injection pressures. Afterwards, an arithmetic mean droplet diameter (D10) and a Sauter mean droplet diameter (D32) as a function of pressure are compared against the measured droplet diameters. Simulated drop size distributions are presented and compared with measured droplet sizes. The results indicate that a high fuel injection pressure increases the liquid penetration length and significantly reduces droplet sizes. The results also exhibit that the SMD decreases from 13.4 μm to 7.5 μm, when injection pressure changes from 40 MPa to 150 MPa and that probability of finding the 5-9 μm droplet diameter decreases from 72% to 40% for the injection pressure drops from 150 MPa to 40 MPa.
The effect of two thermal barrier coatings and their surface roughness on heat transfer, combustion, and emissions has been investigated in a single-cylinder light-duty diesel engine. The evaluated thermal barrier coating materials were plasma-sprayed yttria-stabilized zirconia and hard anodized aluminum, which were applied on the piston top surface. The main tool for the investigation was cylinder pressure analysis of the high-pressure cycle, from which the apparent rate of heat release, indicated efficiency, and heat losses were derived. For verification of the calculated wall heat transfer, the heat flow to the piston cooling oil was measured as well. Application of thermal barrier coatings can influence engine operating conditions like charge temperature and ignition delay. Therefore, extra attention was paid to choosing stable and repeatable engine operating points. The experimental data were modeled using multiple linear regression to isolate the effects of the coatings and of the surface roughness. The results from this study show that high surface roughness leads to increased wall heat losses and a delayed combustion. However, these effects are less pronounced at lower engine loads and in the presence of soot deposits. Both thermal barrier coatings show a reduction of cycle-averaged wall heat losses, but no improvement in indicated efficiency. The surface roughness and thermal barrier coatings had a significant impact on the hydrocarbon emissions, especially for low-load engine operation, while their effect on the other exhaust emissions was relatively small.
The interface/turbulence interaction between two fluids in a turbulent environment has an important role in many technical processes, e.g. primary liquid atomization in combustion devices. Primary atomization has a significant role in spray formation and its characteristics. The resulting dynamics typically span 4–6 orders of magnitude in length scales, making detailed numerical simulations exceedingly expensive. This motivates the need for modeling approaches based on spatial filtering such as large-eddy simulation (LES). In this paper, a new approach based on One-Dimensional turbulence (ODT) is presented to describe the subgrid interface dynamics. ODT is a stochastic model simulating turbulent flow evolution along a notional one-dimensional line of sight by applying instantaneous maps that represent the effects of individual turbulent eddies on property fields. It provides affordable high resolution of interface creation and property gradients within each phase, which are key for capturing the local behavior as well as overall trends. ODT has previously been shown to reproduce the main features of an experimentally determined regime diagram for primary jet breakup. Here a new approach called VODT is presented which produces a size-conditioned as well as a total time rate of generation of droplets for given flow conditions at an interface. At the LES level, the total droplet generation from VODT is interpreted as a rate of mass conversion of LES-resolved liquid into unresolved droplets. Preliminary results of applying VODT to a cell with a planar-shear-layer are discussed at the end of the paper.
The evolution of polydisperse systems is governed by population balance equations. A group of efficient solution approaches are the moment methods, which do not solve for the number density function (NDF) directly but rather for a set of its moments. While this is computationally efficient, a specific challenge arises when describing the fluxes across a boundary in phase space for the disappearance of elements, the so-called zero-flux. The main difficulty is the missing NDF-information at the boundary, which most moment methods cannot provide. Relevant physical examples are evaporating droplets, soot oxidation or particle dissolution.
In general, this issue can be solved by reconstructing the NDF close to the boundary. However, this was previously only achieved with univariate approaches, i.e. considering only a single internal variable. Many physical problems are insufficiently described by univariate population balance equations, e.g. droplets in sprays often require the temperature or the velocity to be internal coordinates in addition to the size.
In this paper, we propose an algorithm, which provides an efficient multivariate approach to calculate the zero-fluxes. The algorithm employs the Extended Quadrature Method of Moments (EQMOM) with Beta and Gamma kernel density functions for the marginal NDF reconstruction and a polynomial or spline for the other conditional dimensions. This combination allows to reconstruct the entire multivariate NDF and based on this, expressions for the disappearance flux are derived. An algorithm is proposed for the whole moment inversion and reconstruction process. It is validated against a suite of test cases with increasing complexity. The influence of the number of kernel density functions and the configuration of the polynomials and splines on the accuracy is discussed. Finally, the associated computational costs are evaluated.
Large-eddy simulation study of combustion cyclic variation in a lean-burn spark ignition engine
(2019)
Multi-cycle large-eddy simulation (LES) was performed to investigate combustion cyclic variability (CCV) in a single cylinder spark ignition engine with a homogeneous lean ( = 1.25) isooctane-air mixture. The aim was to obtain physical insights into the early stage of combustion and its influence on CCV. Propagation of the flame was modeled by a transport equation for the filtered flame surface density within the LES framework. The ignition process was represented by the imposed stretch spark ignition model (ISSIM-LES). Ten consecutive cold flow LES cycles followed by two initialization cycles (12 cycles in total) were used to perform the reactive simulations concurrently. The simulation results were compared with experimental data. Although the number of computed cycles was fairly low, the LES was able to reproduce the cyclic variability observed in experiments both quantitatively and qualitatively. Firstly, validation of the simulation was done by comparing measured pressure traces. Secondly, correlations between the timing of the 10% fuel burnt mass fraction with early flame kernel growth and initial-to-turbulent transition period (in which there was an asymmetric flame kernel that persisted through the early development periods) were determined. Thirdly, calculated results of the flame propagation were analyzed at two cross-sections (in swirl and tumble planes) of the combustion chamber, which highlighted differences in instantaneous flame structures and propagation characteristics between the fastest and slowest cycles. Good overall agreement was obtained between the measurements and simulation data. The results revealed that the instantaneous velocity and fluctuation of flows around the spark vicinity affect growth of the early flame kernel and cause combustion cyclic variability.