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).

The study of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of a single diesel droplet impinging on the wall with different inclined angles (α) is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process.

This paper investigates an optimal design of a diesel engine and after-treatment systems for a series hybrid electric forklift application. A holistic modeling approach is developed in GT-Suite® to establish a model-based hardware definition for a diesel engine and an after-treatment system to accurately predict engine performance and emissions. The used engine model is validated with the experimental data. The engine design parameters including compression ratio, boost level, air-fuel ratio (AFR), injection timing, and injection pressure are optimized at a single operating point for the series hybrid electric vehicle, together with the performance of the after-treatment components. The engine and after-treatment models are then coupled with a series hybrid electric powertrain to evaluate the performance of the forklift in the standard VDI 2198 drive cycle.

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.

Past experimental studies conducted by the current authors on a 13 liter 16.7:1 compression ratio heavy-duty diesel engine have shown that diesel-Compressed Natural Gas (CNG) Reactivity Controlled Compression Ignition (RCCI) combustion targeting low NOx emissions becomes progressively difficult to control as the engine load is increased. This is mainly due to difficulty in controlling reactivity levels at higher loads. For the current study, CFD investigations were conducted in CONVERGE using the SAGE combustion solver with the application of the Rahimi mechanism. Studies were conducted at a load of 5 bar BMEP to validate the simulation results against RCCI experimental data. In the low load study, it was found that the Rahimi mechanism was not able to predict the RCCI combustion behavior for diesel injection timings advanced beyond 30 degCA bTDC. This poor prediction was found at multiple engine speed and load points.

Gasoline compression ignition (GCI) technology shows the potential to obtain high thermal efficiencies while maintaining low soot and NOx emissions in light-duty engine applications. Recent experimental studies and numerical simulations have indicated that high reactivity gasoline-like fuels can further enable the benefits of GCI combustion. However, there is limited empirical data in the literature studying the gasoline compression ignition process at relevant in-cylinder conditions, which are required for further optimizing combustion system designs. This study investigates the temporal and spatial evolution of the compression ignition process of various high reactivity gasoline fuels with research octane numbers (RON) of 71, 74 and 82, as well as a conventional RON 97 E10 gasoline fuel. A ten-hole prototype gasoline injector specifically designed for GCI applications capable of injection pressures up to 450 bar was used.

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.

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.

This study focuses on how to predict hot spots of one of the cylinders of a V8 5.4 L FORD engine running at full load. The KIVA code with conjugate heat transfer capability to simulate the fast transient heat transfer process between the gas and the solid phases has been developed at the Michigan Technological University and will be used in this study. Liquid coolant flow was simulated using FLUENT and will be used as a boundary condition to account for the heat loss to the cooling fluid. In the first step of calculation, the coupling between the gas and the solid phases will be solved using the KIVA code. A 3D transient wall heat flux at the gas-solid interface is then compiled and used along with the heat loss information from the FLUENT data to obtain the temperature distribution for the engine metal components, such as cylinder wall, cylinder head, etc.

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.

One-dimensional single-zone and two-zone analyses have been exercised to calculate the mass fraction burned in an engine operating on ethanol/gasoline-blended fuels using the cylinder pressure and volume data. The analyses include heat transfer and crevice volume effects on the calculated mass fraction burned. A comparison between the two methods is performed starting from the derivation of conservation of energy and the method to solve the mass fraction burned rates through the results including detailed explanation of the observed differences and trends. The apparent heat release method is used as a point of reference in the comparison process. Both models are solved using the LU matrix factorization and first-order Euler integration.

Optimal fuel injection strategies are obtained with a micro-genetic algorithm and an adaptive gradient method for a nonroad, medium-speed DI diesel engine equipped with a multi-orifice, asynchronous fuel injection system. The gradient optimization utilizes a fast-converging backtracking algorithm and an adaptive cost function which is based on the penalty method, where the penalty coefficient is increased after every line search. The micro-genetic algorithm uses parameter combinations of the best two individuals in each generation until a local convergence is achieved, and then generates a random population to continue the global search. The optimizations have been performed for a two pulse fuel injection strategy where the optimization parameters are the injection timings and the nozzle orifice diameters.

A quasi-steady 1-dimensional computer model of a catalyzed particulate filter (CPF) capable of simulating active regeneration of the CPF via diesel fuel injection upstream of a diesel oxidation catalyst (DOC) or other means to increase the exhaust gas temperature has been developed. This model is capable of predicting gaseous species concentrations (HC's, CO, NO and NO2) and exhaust gas temperatures within and after the CPF, for given input values of gaseous species and PM concentrations before the CPF and other inlet variables such as time-varying temperature of the exhaust gas at the inlet of the CPF and volumetric flow rate of exhaust gas.

As the incentive to produce cleaner and more efficient engines increases, diesel engines will become a primary, worldwide solution. Producing diesel engines with higher efficiency and lower emissions requires a fundamental understanding of the interaction of the injected fuel with air as well as with the surfaces inside the combustion chamber. One aspect of this interaction is spray impingement on the piston surface. Impingement on the piston can lead to decreased combustion efficiency, higher emissions, and piston damage due to thermal loading. Modern high-speed diesel engines utilize high pressure common-rail direct-injection systems to primarily improve efficiency and reduce emissions. However, the high injection pressures of these systems increase the likelihood that the injected fuel will impinge on the surface of the piston.

A 2-dimensional numerical model (MTU-FILTER) for a single channel of a honeycomb ceramic diesel particulate trap has been developed. The mathematical modeling of the filtration, flow, heat transfer and regeneration behavior of the particulate trap is described. Numerical results for the pressure drop and particulate mass were compared with existing experimental results. Parametric studies of the diesel particulate trap were carried out. The effects of trap size and inlet temperature on the trap performance are studied using the trap model. An approximate 2-dimensional analytical solution to the simplified Navier-Stokes equations was used to calculate the velocity field of the exhaust flow in the inlet and outlet channels. Assuming a similarity velocity profile in the channels, the 2-dimensional Navier-Stokes equations are approximated by 1-dimenisonal conservation equations, which is similar to those first developed by Bissett.

An analysis of vehicle engine cooling airflow by means of a one-dimensional, transient, compressible flow model was carried out and revealed that similarity theory could be applied to investigate the variation of the airflow with ambient and operating conditions. It was recognized that for a given vehicle engine cooling system, the cooling airflow behavior could be explained using several dimensionless parameters that involve the vehicle speed, fan speed, heat transfer rate through the radiator, ambient temperature and pressure, and the system characteristic dimension. Using the flow resistance and fan characteristics measured from a prototype cooling system and the computer simulation for the one-dimensional compressible flow model, a quantitative correlation of non-dimensional mass flow rate to three dimensionless parameters for a prototype heavy-duty truck was established. The results are presented in charts, tables, and formulas.

A one-dimensional incompressible flow model covering the mechanisms involved in the airflow through an automotive radiator-shroud-fan system with no heat transfer was developed. An analytical expression to approximate the experimentally determined fan performance characteristics was used in conjunction with an analytical approach for this simplified cooling airflow model, and the solution is discussed with illustrations. A major result of this model is a closed form equation relating the transient velocity of the air to the vehicle speed, pressure rise characteristics and speed of the fan, as well as the dimensions and resistance of the radiator. This provides a basis for calculating cooling airflow rate under various conditions. The results of the incompressible flow analysis were further compared with the computational results obtained with a previously developed one-dimensional, transient, compressible flow model.

An accurate estimation of cycle-by-cycle in-cylinder mass and the composition of the cylinder charge is required for spark-ignition engine transient control strategies to obtain required torque, Air-Fuel-Ratio (AFR) and meet engine pollution regulations. Mass Air Flow (MAF) and Manifold Absolute Pressure (MAP) sensors have been utilized in different control strategies to achieve these targets; however, these sensors have response delay in transients. As an alternative to air flow metering, in-cylinder pressure sensors can be utilized to directly measure cylinder pressure, based on which, the amount of air charge can be estimated without the requirement to model the dynamics of the manifold.

It is beneficial but challenging to operate spark-ignition engines under highly lean and dilute conditions. The unstable ignition behavior can result in downgraded combustion performance in engine cylinders. Numerical approach is serving as a promising tool to identify the ignition requirements by providing insight into the complex physical/chemical phenomena. An effort to simulate the early stage of flame kernel initiation in lean and dilute fuel/air mixture has been made and discussed in this paper. The simulations are set to validate against laboratory results of spark ignition behavior in a constant volume combustion vessel. In order to present a practical as well as comprehensive ignition model, the simulations are performed by taking into consideration the discharge circuit analysis, the detailed reaction mechanism, and local heat transfer between the flame kernel and spark plug.

This paper aims to extend the existing Volume of Fluid (VOF) model by implementing an evaporation sub-model in an open source Computational Fluid Dynamics (CFD) software, OpenFOAM. The paper applies the new model to numerically study the evaporation of spherical n-heptane droplets impinging on a hot wall at atmospheric pressure and a temperature above the Leidenfrost temperature. Volume of Fluid (VOF) method is chosen to track the liquid gas interface and the capability of VOF method implemented in interDyMFoam solver of OpenFOAM to simulate hydrodynamics during droplet-droplet interaction and droplet-film interaction is explored. Firstly, the in-built solver is used to simulate problems in isothermal conditions and the simulation results are compared qualitatively with the published results to validate the solver. A numerical method for modeling heat and mass transfer during evaporation is implemented in conjunction with the VOF.