Search Results

Since transient vehicle HVAC computational fluids (CFD) simulations take too long to solve in a production environment, the goal of this project is to automatically create a lumped-parameter flow network from a steady-state CFD that solves nearly instantaneously. The data mining algorithm k-means is implemented to automatically discover flow features and form the network (a reduced order model). The lumped-parameter network is implemented in the commercial thermal solver MuSES to then run as a fully transient simulation. Using this network a “localized heat transfer coefficient” is shown to be an improvement over existing techniques. Also, it was found that the use of the clustering created a new flow visualization technique. Finally, fixing clusters near equipment newly demonstrates a capability to track localized temperatures near specific objects (such as equipment in vehicles).

In this study, the combustion system of a light-duty compression ignition engine running on a market gasoline fuel with Research Octane Number (RON) of 91 was optimized using computational fluid dynamics (CFD) and Machine Learning (ML). This work was focused on optimizing the piston bowl geometry at two compression ratios (CR) (17 and 18:1) and this exercise was carried out at full-load conditions (20 bar indicated mean effective pressure, IMEP). First, a limited manual piston design optimization was performed for CR 17:1, where a couple of pistons were designed and tested. Thereafter, a CFD design of experiments (DoE) optimization was performed where CAESES, a commercial software tool, was used to automatically perturb key bowl design parameters and CONVERGE software was utilized to perform the CFD simulations. At each compression ratio, 128 piston bowl designs were evaluated.

A Machine Learning-Genetic Algorithm (ML-GA) approach was developed to virtually discover optimum designs using training data generated from multi-dimensional simulations. Machine learning (ML) presents a pathway to transform complex physical processes that occur in a combustion engine into compact informational processes. In the present work, a total of over 2000 sector-mesh computational fluid dynamics (CFD) simulations of a heavy-duty engine were performed. These were run concurrently on a supercomputer to reduce overall turnaround time. The engine being optimized was run on a low-octane (RON70) gasoline fuel under partially premixed compression ignition (PPCI) mode. A total of nine input parameters were varied, and the CFD simulation cases were generated by randomly sampling points from this nine-dimensional input space. These input parameters included fuel injection strategy, injector design, and various in-cylinder flow and thermodynamic conditions at intake valve closure (IVC).

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.

In order to extend the operability limit of the gasoline compression ignition (GCI) engine, as an avenue for low temperature combustion (LTC) regime, the effects of parametric variations of engine operating conditions on the performance of six-stroke GCI (6S-GCI) engine cycle are numerically investigated, using an in-house 3D CFD code coupled with high-fidelity physical sub-models along with the Chemkin library. The combustion and emissions were calculated using a skeletal chemical kinetics mechanism for a 14-component gasoline surrogate fuel. Authors’ previous study highlighted the effects of the variation of injection timing and split ratio on the overall performance of 6S-GCI engine and the unique mixing-controlled burning mode of the charge mixtures during the two additional strokes. As a continuing effort, the present study details the parametric studies of initial gas temperature, boost pressure, fuel injection pressure, compression ratio, and EGR ratio.

Increasing regulatory demand for efficiency has led to development of novel combustion modes such as HCCI, GCI and RCCI for gasoline light duty engines. In order to realize HCCI as a compression ignition combustion mode system, in-cylinder compression temperatures must be elevated to reach the autoignition point of the premixed fuel/air mixture. This should be co-optimized with appropriate fuel formulations that can autoignite at such temperatures. CFD combustion modeling is used to model the auto ignition of gasoline fuel under compression ignition conditions. Using the fully detailed fuel mechanism consisting of thousands of components in the CFD simulations is computationally expensive. To overcome this challenge, the real fuel is represented by few major components of create a surrogate fuel mechanism. In this study, 9 variations of gasoline fuel sets were chosen as candidates to run in HCCI combustion mode.

Achieving robust ignitability for compression ignition of diesel engines at cold conditions is traditionally challenging due to insufficient fuel vaporization, heavy wall impingement, and thick wall films. Gasoline compression ignition (GCI) has shown the potential to offer an enhanced NOx-particulate matter tradeoff with diesel-like fuel efficiency, but it is unknown how the volatility and reactivity of the fuel will affect ignition under very cold conditions. Therefore, it is important to investigate the impact of fuel physical and chemical properties on ignition under pressures and temperatures relevant to practical engine operating conditions during cold weather. In this paper, 0-D and 3-D computational fluid dynamics (CFD) simulations of GCI combustion at cold conditions were performed.

Fuels in the gasoline auto-ignition range (Research Octane Number (RON) > 60) have been demonstrated to be effective alternatives to diesel fuel in compression ignition engines. Such fuels allow more time for mixing with oxygen before combustion starts, owing to longer ignition delay. Moreover, by controlling fuel injection timing, it can be ensured that the in-cylinder mixture is “premixed enough” before combustion occurs to prevent soot formation while remaining “sufficiently inhomogeneous” in order to avoid excessive heat release rates. Gasoline compression ignition (GCI) has the potential to offer diesel-like efficiency at a lower cost and can be achieved with fuels such as low-octane straight run gasoline which require significantly less processing in the refinery compared to today’s fuels.

Combustion systems with advanced injection strategies have been extensively studied, but there still exists a significant fundamental knowledge gap on fuel spray interactions with the piston surface and chamber walls. This paper is meant to provide detailed data on spray-wall impingement physics and support the spray-wall model development. The experimental work of spray-wall impingement with non-vaporizing spray characterization, was carried out in a high pressure-temperature constant-volume combustion vessel. The simultaneous Mie scattering of liquid spray and schlieren of liquid and vapor spray were carried out. Diesel fuel was injected at a pressure of 1500 bar into ambient gas at a density of 22.8 kg/m3 with isothermal conditions (fuel, ambient, and plate temperatures of 423 K). A Lagrangian-Eulerian modeling approach was employed to characterize the spray-gas and spray-wall interactions in the CONVERGETM framework by means of a Reynolds-Averaged Navier-Stokes (RANS) formulation.

Cyclic variability in internal combustion engines (ICEs) arises from multiple concurrent sources, many of which remain to be fully understood and controlled. This variability can, in turn, affect the behavior of the engine resulting in undesirable deviations from the expected operating conditions and performance. Shot-to-shot variation during the fuel injection process is strongly suspected of being a source of cyclic variability. This study focuses on the shot-to-shot variability of injector needle motion and its influence on the internal nozzle flow behavior using diesel fuel. High-speed x-ray imaging techniques have been used to extract high-resolution injector geometry images of the sac, orifices, and needle tip that allowed the true dynamics of the needle motion to emerge. These measurements showed high repeatability in the needle lift profile across multiple injection events, while the needle radial displacement was characterized by a much higher degree of randomness.

A 2-D CFD model was developed to describe the heat transfer, and reaction kinetics in a honeycomb structured ceramic diesel particulate trap. This model describes the steady state as well as the transient behavior of the flow and heat transfer during the trap regeneration processes. The trap temperature profile was determined by numerically solving the 2-D unsteady energy equation including the convective, heat conduction and viscous dissipation terms. The convective terms were based on a 2-D analytical flow field solution derived from the conservation of mass and momentum equations (Opris, 1997). The reaction kinetics were described using a discretized first order Arrhenius function. The 2-D term describing the reaction kinetics and particulate matter conservation of mass was added to the energy equation as a source term in order to represent the particulate matter oxidation. The filtration model describes the particulate matter accumulation in the trap.

A global optimization method has been developed for an engine simulation code and utilized in the search of optimal fuel injection strategies. This method uses a Lagrange interpolation function which interpolates engine output data generated at the vertices and the intermediate points of the input parameters. This interpolation function is then used to find a global minimum over the entire parameter set, which in turn becomes the starting point of a CFD-based optimization. The CFD optimization is based on a steepest descent method with an adaptive cost function, where the line searches are performed with a fast-converging backtracking algorithm. The adaptive cost function is based on the penalty method, where the penalty coefficient is increased after every line search. The parameter space is normalized and, thus, the optimization occurs over the unit cube in higher-dimensional space.

Transient engine operations are modeled and simulated with a 1-D code (GT Power) using heat release and emission data computed by a 3-D CFD code (Kiva3). During each iteration step of a transient engine simulation, the 1-D code utilizes the 3-D data to interpolate the values for heat release and emissions. The 3-D CFD computations were performed for the compression and combustion stroke of strategically chosen engine operating points considering engine speed, torque and excess air. The 3-D inlet conditions were obtained from the 1-D code, which utilized 3-D heat release data from the previous 1-D unsteady computations. In most cases, only two different sets of 3-D input data are needed to interpolate the transient phase between two engine operating points. This keeps the computation time at a reasonable level. The results are demonstrated on the load response of a generator which is driven by a medium-speed diesel engine.

The development of new high power diesel engines is continually going for increased mean effective pressures and consequently increased thermal loads on combustion chamber walls close to the limits of endurance. Therefore accurate CFD simulation of conjugate heat transfer on the walls becomes a very important part of the development. In this study the heat transfer and temperature on piston surface was studied using conjugate heat transfer model along with a variety of near wall treatments for turbulence. New wall functions that account for variable density were implemented and tested against standard wall functions and against the hybrid near wall treatment readily available in a CFD software Star-CD.

The discrete droplet model is widely used to describe two-phase flows such as high-velocity dense sprays. The interaction between the liquid and the gas phase is modeled via appropriate source terms in the gas phase equations. This approach can lead to a strong dependence of the liquid-gas coupling on the spatial resolution of the gas phase. The liquid-gas coupling requires the computation of source terms using the gas phase properties, and, subsequently, these sources are then distributed onto the gas phase mesh. In this study, a Lagrange polynomial interpolation method has been developed to evaluate the source terms and also to distribute these source terms onto the gas mesh. The focus of this investigation has been on the momentum exchange between the two phases. The Lagrange polynomial interpolation and source term distribution methods are evaluated for non-evaporating sprays using KIVA3 as a modeling platform.

This study presents diesel spray breakup regimes and the wave model basic theory from literature. The RD wave model and the KH-RT wave model are explained. The implementation of the KH-RT wave model in a commercial CFD code is briefly presented. This study relies on experimental data from non-evaporating sprays that have earlier been measured at Helsinki University of Technology. The simulated fuel spray in a medium-speed diesel engine had a satisfactory match with the experimental data. The KH-RT wave model resulted in a much faster drop breakup than with the RD wave model. This resulted in a thin spray core with the KH-RT model. The fuel viscosity effect on drop sizes was well predicted by the KH-RT wave model.

Impingement of jet-to-jet has been found to give improved spray penetration characteristics and higher vaporization rates when compared to multi-hole outwardly injecting fuel injectors which are commonly used in the gasoline engine. The current work studies a non-reacting spray by using a 5-hole impinging-jet style direct-injection injector. The jet-to-jet collision induced by the inwardly opening nozzles of the multi-hole injector produces rapid and short jet breakup which is fundamentally different from how conventional fuel injectors operate. A non-reacting spray study is performed using a 5-hole impinging jet injector and a traditional 6-hole Bosch Hochdruck-Einspritzventil (HDEV)-5 gasoline direct-injection (GDI) injector with gasoline as a fuel injected at 172 bar pressure with ambient temperature of 653 K and 490 K and ambient pressure of 37.4 bar and 12.4 bar.

Dimethyl Ether (DME) is considered a clean alternative fuel to diesel due to its soot-free combustion characteristics and its capability to be produced from renewable energy sources rather than fossil fuels such as coal or petroleum. To mitigate the effect of strong wave dynamics on fuel supply lines caused due to the high compressibility of DME and to overcome its low lubricity, a hydraulically actuated electronic unit injector (HEUI) with pressure intensification was used. The study focuses on high pressure operation, up to 2000 bar, significantly higher than pressure ranges reported previously with DME. A one-dimensional HEUI injector model is built in MATLAB/SIMULINK graphical software environment, to predict the rate of injection (ROI) profile critical to spray and combustion characterization.

In this study a detailed 1-D engine system model coupled with 3-D computational fluid dynamics (CFD) analysis was used to investigate the air system design requirements for a heavy duty diesel engine operating with low reactivity gasoline-like fuel (RON70) under partially premixed combustion (PPC) conditions. The production engine used as the baseline has a geometric compression ratio (CR) of 17.3 and the air system hardware consists of a 1-stage variable geometry turbine (VGT) with a high pressure exhaust gas recirculation (HP-EGR) loop. The analysis was conducted at six engine operating points selected from the heavy-duty supplemental emissions test (SET) cycle, i.e., A75, A100, B25, B50, B75, and C100. The engine-out NOx target was set at 1 g/hp-hr (1.34 g/kWh) to address a future hypothetical tailpipe NOx limit of 0.02 g/hp-hr (0.027 g/kWh) while an engine-out particulate matter (PM) target of 0.01 g/hp-hr (0.013 g/kWh) was selected to comply with existing EPA 2010 regulations.

Bio-derived fuels are drawing more and more attention in the internal combustion engine (ICE) research field in recent years. Those interests in use of renewable biofuels in ICE applications derive from energy security issues and, more importantly, from environment pollutant emissions concerns. High fidelity numerical study of engine combustion requires advanced computational fluid dynamics (CFD) to be coupled with detailed chemical kinetic models. This task becomes extremely challenging if real fuels are taken into account, as they include a mixture of hundreds of different hydrocarbons, which prohibitively increases computational cost. Therefore, along with employing surrogate fuel models, reduction of detailed kinetic models for multidimensional engine applications is preferred. In the present work, a reduced mechanism was developed for primary reference fuel (PRF) using the directed relation graph (DRG) approach. The mechanism was generated from an existing detailed mechanism.