Search Results

An experimental study has been conducted to provide insight into heat transfer to the piston of a light-duty single-cylinder research engine under Conventional Diesel (CDC), Homogeneous Charge Compression Ignition (HCCI), and Reactivity Controlled Compression Ignition (RCCI) combustion regimes. Two fast-response surface thermocouples embedded in the piston top measured transient temperature. A commercial wireless telemetry system was used to transmit thermocouple signals from the moving piston. A detailed comparison was made between the different combustion regimes at a range of engine speed and load conditions. The closed-cycle integrated and peak heat transfer rates were found to be lower for HCCI and RCCI when compared to CDC. Under HCCI operation, the peak heat transfer rate showed sensitivity to the 50% burn location.

The focus of the present study was to characterize Reactivity Controlled Compression Ignition (RCCI) using a single-fuel approach of gasoline and gasoline mixed with a commercially available cetane improver on a multi-cylinder engine. RCCI was achieved by port-injecting a certification grade 96 research octane gasoline and direct-injecting the same gasoline mixed with various levels of a cetane improver, 2-ethylhexyl nitrate (EHN). The EHN volume percentages investigated in the direct-injected fuel were 10, 5, and 2.5%. The combustion phasing controllability and emissions of the different fueling combinations were characterized at 2300 rpm and 4.2 bar brake mean effective pressure over a variety of parametric investigations including direct injection timing, premixed gasoline percentage, and intake temperature. Comparisons were made to gasoline/diesel RCCI operation on the same engine platform at nominally the same operating condition.

In recent years, the increase of gasoline fuel injection pressure is a way to improve thermal efficiency and lower engine-out emissions in GDI homogenous combustion concept. The challenge of controlling particulate formation as well in mass and number concentrations imposed by emissions regulations can be pursued improving the mixture preparation process and avoiding mixture inhomogeneity with ultra-high injection pressure values up to 100 MPa. The increase of the fuel injection pressure in GDI homogeneous systems meets the demand for increased injector static flow, while simultaneously improves the spray atomization and mixing characteristics with consequent better combustion performance. Few studies quantify the effects of high injection pressure on transient gasoline spray evolution. The aim of this work was to simulate with OpenFOAM the spray morphology of a commercial gasoline injected in a constant volume vessel by a prototypal GDI injector.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

In this paper, we studied the accuracy of computational modeling of the ignition of a pilot injectionin the Sandia National Laboratories (SNL) light-duty optical engine facility, using the physical properties of a cetane/iso-cetane Diesel Primary Reference Fuel (DPRF) mixture and the reaction kinetics of a well-validated mechanism for primary reference fuels. Local fuel-air equivalence ratio measurements from fuel tracer based planar laser-induced fluorescence (PLIF) experiments were used to compare the mixture formation predictions with KIVA-ERC-based simulations. The effects of variations in injection mass from 1 mg to 4 mg, in-cylinder swirl ratio, and near-TDC temperatures on non-combusting mixture preparation were analyzed, to assess the accuracy of the model in capturing average jet behavior, despite its inability to model the non-negligible jet-by-jet variations seen in the experiments.

Dynamic tests have been performed on carbon fiber racing seats following the FIA regulations. The tests have shown, in rear impact tests, a relatively strong rebound leading to large forward bending of neck, and, in side impact tests, very large lateral displacement of the head, the latter protruding dangerously towards hard portions of the car structure. Stiffening the seat back by steel struts results in reducing strongly both the motion and the acceleration of the head. Simulations of the dynamics of the tests have been done with multi-body models, including the Hybrid III dummy and seat deflection, by means of the program VEDYAC. It has been found that computer simulation can predict very accurately the result of a test, provided the numerical models have been carefully calibrated to match the dummy tolerance bands. Once they have been calibrated and validated with a number of tests, the computer models can be very useful to extend the test results to different test conditions.

Partially premixed combustion (PPC) can be applied to decrease emissions and increase fuel efficiency in direct injection, compression ignition (DICI) combustion engines. PPC is strongly influenced by the mixing of fuel and oxidizer, which for a given fuel is controlled mainly by (a) the fuel injection, (b) the in-cylinder flow, and (c) the geometry and dynamics of the engine. As the injection timings can vary over a wide range in PPC combustion, detailed knowledge of the in-cylinder flow over the whole intake and compression strokes can improve our understanding of PPC combustion. In computational fluid dynamics (CFD) the in-cylinder flow is sometimes simplified and modeled as a solid-body rotation profile at some time prior to injection to produce a realistic flow field at the moment of injection. In real engines, the in-cylinder flow motion is governed by the intake manifold, the valve motion, and the engine geometry.

The current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

The CFD tools in aerodynamic design process have been commonly used in aerospace industry in last three decades. Although there are many CFD software algorithms developed for aerodynamic applications, the nature of a complex, three-dimensional geometry in incompressible highly separated, viscous flow made computational simulation of ground vehicle aerodynamics more difficult than aerospace applications. However, recent developments in computational hardware and software industry enabled many new engineering applications on computational environment. Traditional production process has largely influenced by computational design, analysis, manufacturing and visualization. Different aspects of linking advanced computational tools and aerodynamic vehicle design challenges are discussed in the present work. Key technologies like parallel computation, turbulence modeling and CFD/wind tunnel compatibility issues are presented.

The influence of spray-wall interaction on air entrainment in an unsteady non-evaporating diesel spray was studied using laser Doppler anemometry. The spray was injected into confined quiescent air at ambient pressure and temperature and made to impact on a flat wall. The air velocity component normal to a cylindrical surface surrounding the spray was measured during the entire injection period, allowing to evaluate the time history of the entrained air mass flow rate. The influence of wall distance and spray impingement angle on air entrainment characteristics has been investigated and the results indicate that the presence of a wall increases the entrained mass flow rate in the region close to the surface, during the main injection period. Normal impingement appears to produce stronger effects than oblique incidence at 30 and 45 deg. A qualitative explanation of the results is also proposed, based on the drop-gas momentum exchange mechanism.

This work presents a comprehensive numerical model for ice accretion and Ice Protection System (IPS) simulation over a 2D component, such as an airfoil. The model is based on the Myers model for ice accretion and extended to include the possibility of a heated substratum. Six different icing conditions that can occur during in-flight ice accretion with an Electro-Thermal Ice Protection System (ETIPS) activated are identified. Each condition presents one or more layers with a different water phase. Depending on the heat fluxes, there could be only liquid water, ice, or a combination of both on the substratum. The possible layers are the ice layer on the substratum, the running liquid film over ice or substratum, and the static liquid film between ice and substratum caused by ice melting. The last layer, which is always present, is the substratum. The physical model that describes the evolution of these layers is based on the Stefan problem. For each layer, one heat equation is solved.

This paper presents a novel fully-automatic remeshing procedure, based on the level-set method and Delaunay triangulation, to model three-dimensional boundary problems and generate a new conformal body-fitted mesh. The proposed methodology is applied to long-term in-flight ice accretion, which is characterized by the formation of extremely irregular ice shapes. Since ice accretion is coupled with the aerodynamic flow field, a multi-step procedure is implemented. The total icing exposure time is subdivided into smaller time steps, and at each time step a three-dimensional body-fitted mesh, suitable for the computation of the aerodynamic flow field around the updated geometry, is generated automatically. The methodology proposed can effectively deal with front intersections, as shown with a manufactured example.

Diesel piston bowl geometry can affect turbulent mixing and therefore it impacts heat-release rates, thermal efficiency, and soot emissions. The focus of this work is on the effects of bowl geometry and injection timing on turbulent flow structure. This computational study compares engine behavior with two pistons representing competing approaches to combustion chamber design: a conventional, re-entrant piston bowl and a stepped-lip piston bowl. Three-dimensional computational fluid dynamics (CFD) simulations are performed for a part-load, conventional diesel combustion operating point with a pilot-main injection strategy under non-combusting conditions. Two injection timings are simulated based on experimental findings: an injection timing for which the stepped-lip piston enables significant efficiency and emissions benefits, and an injection timing with diminished benefits compared to the conventional, re-entrant piston.

Sector mesh modeling is the dominant computational approach for combustion system design optimization. The aim of this work is to quantify the errors descending from the sector mesh approach through three geometric modeling approaches to an optical diesel engine. A full engine geometry mesh is created, including valves and intake and exhaust ports and runners, and a full-cycle flow simulation is performed until fired TDC. Next, an axisymmetric sector cylinder mesh is initialized with homogeneous bulk in-cylinder initial conditions initialized from the full-cycle simulation. Finally, a 360-degree azimuthal mesh of the cylinder is initialized with flow and thermodynamics fields at IVC mapped from the full engine geometry using a conservative interpolation approach. A study of the in-cylinder flow features until TDC showed that the geometric features on the cylinder head (valve tilt and protrusion into the combustion chamber, valve recesses) have a large impact on flow complexity.

Turbulence modeling for fuel spray simulation plays a prominent role in the understanding of the flow behavior in Internal Combustion Engines (ICEs). Currently, a lot of research work is actively spent on Large Eddy Simulation (LES) turbulence modeling as a replacement option of standard Reynolds averaged approaches in the Eulerian-Lagrangian spray modeling framework, due to its capability to accurately describe flow-induced spray variability and to the lower dependence of the results on the specific turbulence model and/or modeling coefficients. The introduction of LES poses, however, additional questions related to the implementation/adaptation of spray-related turbulence sources and to the rise of conflicting numerics and grid requirements between the Lagrangian and Eulerian parts of the simulated flow.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.

Prediction of in-cylinder flows and fuel-air mixing are two fundamental pre-requisites for a successful simulation of direct-injection engines. Over the years, many efforts were carried out in order to improve available turbulence and spray models. However, enhancements in physical modeling can be drastically affected by how the mesh is structured. Grid quality can negatively influence the prediction of organized charge motion structures, turbulence generation and interaction between in-cylinder flows and injected sprays. This is even more relevant for modern direct injection engines, where multiple injections and control of charge motions are employed in a large portion of the operating map. Currently, two different approaches for mesh generation exist: manual and automatic. The first makes generally possible to generate high-quality meshes but, at the same time, it is very time consuming and not completely free from user errors.

The paper focuses on the development of a mesh moving method based on non-conformal topologically changing grids applied to the simulation of IC engines, where the prescribed motion of piston and valves is accomplished by rigidly translating the sub-domain representing the moving component. With respect to authors previous work, a more robust and efficient algorithm to handle the connectivity of non-conformal interfaces and a mesh-motion solver supporting multiple layer addition/removal of cells, to decouple the time-step constraints of the mesh motion and of the fluid dynamics, has been implemented as a C++ library to extend the already existing classes for dynamic mesh handling of the finite-volume, open-source CFD code OpenFOAM®. Other new features include automatic decomposition of large multiple region domains to preserve processors load balance with topological changes for parallel computations and additional tools for automatic preprocessing and case setup.