Search Results

Premixed Charge Compression Ignition (PCCI), a Low Temperature Combustion (LTC) strategy for diesel engines is of increasing interest due to its potential to simultaneously reduce soot and NOx emissions. However, the influence of mixture preparation on combustion phasing and heat release rate in LTC is not fully understood. In the present study, the influence of injection timing on mixture preparation, combustion and emissions in PCCI mode is investigated by experimental and computational methods. A sequential coupling approach of 3D CFD with a Stochastic Reactor Model (SRM) is used to simulate the PCCI engine. The SRM accounts for detailed chemical kinetics, convective heat transfer and turbulent micro-mixing. In this integrated approach, the temperature-equivalence ratio statistics obtained using KIVA 3V are mapped onto the stochastic particle ensemble used in the SRM.

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 recent years urea selective catalytic reduction (SCR) has become the principal method of NOx abatement within heavy duty (HD) diesel exhaust systems; however, with upcoming applications demanding NOx reduction efficiencies of above 96 % on engines producing upwards of 10 g·kWh−1 NOx, future diesel exhaust fluid (DEF) dosing systems will be required to operate stably at significantly increased dosing rates. Developing a dosing system capable of meeting the increased performance requirements demands an improved understanding of how DEF sprays interact with changing exhaust flows. This study has investigated four production systems representing a diverse range of dosing strategies in order to determine how performance is influenced by spray structure and identify promising strategies for further development. The construction of an optically accessible hot-air flow rig has enabled visualisation of DEF injection into flows representative of HD diesel exhaust conditions.

The drag coefficient at zero yaw angle is the single parameter usually used to define the aerodynamic drag characteristics of a passenger car. However, this is usually the minimum drag condition and will, for example, lead to an underestimate of the effect of aerodynamic drag on fuel consumption because the important influence of the natural wind has been excluded. An alternative measure of aerodynamic drag should take into account the effect of nonzero yaw angles and a variant of wind-averaged drag is suggested as the best option. A wind-averaged drag coefficient (CDW) is usually derived for a particular vehicle speed using a representative wind speed distribution. In the particular case where the road speed distribution is specified, as for a drive cycle to determine fuel economy, a relevant drag coefficient can be derived by using a weighted road speed.

In a real-world environment, a vehicle on the road is subjected to a range of flow yaw angles, the most severe of which can impact handling and stability. A fully coupled, six degrees-of-freedom CFD and vehicle handling simulation has modelled the complete closed loop system. Varying flow yaw angles are introduced via time dependent boundary conditions and aerodynamic loads predicted, whilst a handling model running simultaneously calculates the resulting vehicle response. Updates to the vehicle position and orientation within the CFD simulation are achieved using the overset grid method. Using this approach, a crosswind simulation that follows the parameters of ISO 12021:2010 (Sensitivity to lateral wind - Open-loop test method using wind generator input), was performed using the fastback variant of the DrivAer model. Fully coupled aerodynamic and vehicle response was compared to that obtained using the simplified quasi-steady and unsteady, one way coupled method.

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.

Fuel lean combustion and exhaust gas dilution are known to increase the thermal efficiency and reduce NOx emissions. In this study, experiments are performed to understand the effect of equivalence ratio on flame kernel formation and flame propagation around the spark plug for different low turbulent velocities. A series of experiments are carried out for propane-air mixtures to simulate engine-like conditions. For these experiments, equivalence ratios of 0.7 and 0.9 are tested with 20 percent mass-based exhaust gas recirculation (EGR). Turbulence is generated by a shrouded fan design in the vicinity of J-spark plug. A closed loop feedback control system is used for the fan to generate a consistent flow field. The flow profile is characterized by using Particle Image Velocimetry (PIV) technique. High-speed Schlieren visualization is used for the spark formation and flame propagation.

In the wind tunnel the effect of a wind input on the aerodynamic characteristics of any road vehicle is simulated by yawing the vehicle. This represents a wind input where the wind velocity is constant with height above the ground. In reality the natural wind is a boundary layer flow and is sheared so that the wind velocity will vary with height. A CFD simulation has been conducted to compare the aerodynamic characteristics of a DrivAer model, in fastback and squareback form, subject to a crosswind flow, with and without shear. The yaw simulation has been carried out at a yaw angle of 10° and with one shear flow exponent. It is shown that the car experiences almost identical forces and moments in the two cases when the mass flow in the crosswind over the height of the car is similar. Load distributions are presented for the two cases. The implications for wind averaged drag are discussed.

The motivation for this paper is to predict the flow of water over exterior surfaces of road vehicles. We present simulations of liquid flows on solid surfaces under the influence of gravity with and without the addition of aerodynamic forces on the liquid. This is done using an implementation of a Coupled Level Set Volume of Fluid method (CLSVOF) multiphase approach implemented in the open source OpenFOAM CFD code. This is a high fidelity interface-resolving method that solves for the velocity field in both phases without restrictions on the flow regime. In the current paper the suitability of the approach to Exterior Water Management (EWM) is demonstrated using the representative test cases of a continuous liquid rivulet flowing along an inclined surface with a channel located downstream perpendicular to the oncoming flow.

Fundamental understanding of the sources of fuel-derived Unburned Hydrocarbon (UHC) emissions in heavy duty diesel engines is a key piece of knowledge that impacts engine combustion system development. Current emissions regulations for hydrocarbons can be difficult to meet in-cylinder and thus after treatment technologies such as oxidation catalysts are typically used, which can be costly. In this work, Computational Fluid Dynamics (CFD) simulations are combined with engine experiments in an effort to build an understanding of hydrocarbon sources. In the experiments, the combustion system design was varied through injector style, injector rate shape, combustion chamber geometry, and calibration, to study the impact on UHC emissions from mixing-controlled diesel combustion.

The motivation for this paper is to consider the effect of rear end geometry on rear soiling using a representative generic SUV body. In particular the effect of varying the top slant angle is considered using both experiment and Computational Fluid Dynamics (CFD). Previous work has shown that slant angle has a significant effect on wake shape and drag and the work here extends this to investigate the effect on rear soiling. It is hoped that this work can provide an insight into the likely effect of such geometry changes on the soiling of similarly shaped road vehicles. To increase the generality of results, and to allow comparison with previously obtained aerodynamic data, a 25% scale generic SUV model is used in the Loughborough University Large Wind Tunnel. UV doped water is sprayed from a position located at the bottom of the left rear tyre to simulate the creation of spray from this tyre.

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.

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.

Powertrain thermal encapsulation has the potential to improve fuel consumption and CO2 via heat retention. Heat retained within the powertrain after a period of engine-off, can increase the temperature of the next engine start hours after key-off. This in turn reduces inefficiencies associated with sub-optimal temperatures such as friction. The Ambient Temperature Correction Test was adopted in the current work which contains two World-wide harmonised Light duty Test Procedure (WLTP) cycles separated by a 9-hour soak period. A coupled 1D - 3D computational approach was used to capture heat retention characteristics and subsequent warm-up effects. A 1-D powertrain warm-up model was developed in GT-Suite to capture the thermal warm-up characteristics of the powertrain. The model included a temperature dependent friction model, the thermal-hydraulic characteristics of the cooling and lubrication circuits as well as parasitic losses associated with pumps.

Significant aerodynamic drag reduction is obtained on a bluff body by tapering the rear body. In the 1930’s it was found that a practical low drag car body could be achieved by cutting off the tail of a streamlined shape. The rear end of a car with a truncated tail is commonly referred to as a Kamm back. It has often been interpreted as implying that the drag of this type of body is almost the same as that for a fully streamlined shape. From a review of the limited research into truncated streamlined tails it is shown in this paper that, while true for some near axisymmetric bodies, it is not the case for many more car-like shapes. For these shapes the drag reduction from an elongated tail varies almost linearly with the reduction in cross section area. A CFD simulation to determine the drag reduction from a truncated streamlined tail of variable length on the simple Windsor Body is shown by way of confirmation.

It is estimated that operating continuously on a B20 fuel containing the current allowable ASTM specification limits for metal impurities in biodiesel could result in a doubling of ash exposure relative to lube-oil-derived ash. The purpose of this study was to determine if a fuel containing metals at the ASTM limits could cause adverse impacts on the performance and durability of diesel emission control systems. An accelerated durability test method was developed to determine the potential impact of these biodiesel impurities. The test program included engine testing with multiple DPF substrate types as well as DOC and SCR catalysts. The results showed no significant degradation in the thermo-mechanical properties of cordierite, aluminum titanate, or silicon carbide DPFs after exposure to 150,000 mile equivalent biodiesel ash and thermal aging. However, exposure of a cordierite DPF to 435,000 mile equivalent aging resulted in a 69% decrease in the thermal shock resistance parameter.

Single-cylinder engine experiments and computational fluid dynamics (CFD) modeling were used in this study to conduct a comprehensive evaluation of the accuracy of the modeling approach, with a focus on soot emissions. A semi-empirical soot model, the classic two-step Hiroyasu model with Nagle and Strickland-Constable oxidation, was used. A broad range of direct-injected (DI) combustion systems were investigated to assess the predictive accuracy of the soot model as a design tool for modern DI diesel engines. Experiments were conducted on a 2.5 liter single-cylinder engine. Combustion system combinations included three unique piston bowl shapes and seven variants of a common rail fuel injector. The pistons included a baseline “Mexican hat” piston, a reentrant piston, and a non-axisymmetric piston similar to the Volvo WAVE design. The injectors featured six or seven holes and systematically varied included angles from 120 to 150 degrees and hole sizes from 170 to 273 μm.

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.