With stringent emission regulations, it has become more important for modern diesel engine manufacturers to accurately predict engine-out nitrogen oxide (NOx) emissions across a wide range of operating conditions. Thermal NO is the major source of engine-out NOx in modern diesel engines. For thermal NO formation, several earlier studies have recommended the forward and reverse reaction rate coefficients of the rate-limiting reaction (O + N2 ⇌ NO + N). However, due to deficiencies in sub-models and inadequacies of reduced chemical mechanisms to represent diesel combustion, these recommended values more often than not need to be adjusted in reduced order combustion models to accurately predict engine-out NOx. Hence, in this work a systematic and computationally efficient approach has been proposed to streamline the process of determining the optimum reaction rate coefficients. To develop the optimization approach, four different production diesel engines with different operating conditions in terms of speed, load, and exhaust-gas recirculation have been considered. Numerical simulations have been performed using a detailed zero-dimensional velocity-composition-frequency transported probability density function (0D-VCF-tPDF) model that uses hundreds of notional particles to capture in-cylinder stratification. Four different combinations of hydrocarbon and NOx chemical mechanisms were used to represent chemistry. It was found that for the rate-limiting reaction, the pre-exponent factors (Af1,Ar1) and activation energies (EA,f1,EA,r1) of the forward and reverse reaction rates follow a linear band in Af1−EA,f1 and Ar1−EA,r1 space where predicted engine-out NOx match the measured values closely. By encompassing such bands from different engines and considering constraints on activation energies, a reduced search domain of pre-exponent factors and activation energies was constructed that is expected to be applicable to any diesel engine. Eventually, computationally efficient three-line and one-line search approaches were proposed to determine the optimum values of the pre-exponent factors and activation energies that led to a minimum error between measured and predicted engine-out NOx. Finally, these three-line and one-line NOx optimization approaches were applied to a fifth production diesel engine for which the 0D-VCF-tPDF model showed a very good predictive performance in terms of predicting peak pressure, 50% burn rate, and engine-out NOx when compared to measured and 3D-CFD values.
Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.
Proper thermal management is key to the safe and efficient operation of electric motors, which relies on accurate prediction of electromagnetic loss. Electromagnetics is often coupled with thermal analysis to provide detailed loss information, which is referred as one-way coupled Multiphysics analysis. While one-way coupling improves model accuracy by incorporating detailed loss information, the electromagnetic analysis still relies on assumed temperature. To close the loop, a two-way coupling has been proposed to couple thermal information back to the electromagnetics analysis. The two-way coupling provides a complete modeling framework by considering the two physics via temperature dependent material properties, leading to improved electromagnetics accuracy in addition to thermal analysis. However, its application has been limited due to the complicated mapping and computing overhead. In this study, the two-way coupling has been investigated over full 3D geometry to examine benefits and pitfalls, focused on the energy conservation between the two physics.
Three stator cooling jacket concepts, including a new cooling jacket concept with parallel straight spiral channels, are evaluated for electric machine of mild hybrid electric car. Steady-state cooling jacket thermal flow analysis and transient-state full motor simulations with circumferential, axial, and spiral jacket concepts are carried out to predict their heat transfer performance and temperature distribution. A prototype motor has been made with the new jacket concept. Infrared thermography is performed to detect the transient temperature distribution of the motor stator windings and laminations for validation of the thermal flow analysis.
Fuel film that adhered on engine walls from spray impingement is considered a primary source of harmful combustion emissions. However, the physics of the wall film formation, propagation, and breakup is not fully understood yet because of its multiphase nature. Existing literature has revealed that the mass transportation within the fuel film takes a wave propagation form. This article aims to identify the dynamics of the wall film during spray impingement via high-speed laser diagnostics. In this work, a single-hole injector was used and the spray impinged onto a stage made of sapphire glass for diagnostics purposes. Iso-octane was used as the fuel and 10% ethanol was blended to dope rhodamine 6G as the fluorescent species for laser excitation. Simultaneous optical measurements, such as laser-induced fluorescence and Mie scattering, are performed to obtain the characteristics of the wall film quantitatively. Various aspects of the wall film, including the frequency of the wave, wave speed, and wave height, are inspected. The impact of fuel temperature and wall temperature under typical cold-start conditions are also studied to investigate the temperature dependence of the wall film dynamics. The research is also intended to provide quantitative experimental data for numerical models for impingement prediction and thus help the process of emission reduction and combustion optimization of internal combustion engines.
A computational approach to predict external spray characteristics for flashing and cavitating nozzles was presented and validated. It is developed as a fully Eulerian and compressible two-phase flow solver that simulates vaporization and condensation of the fuel using a Homogeneous Relaxation Model (HRM). The flow solver together with the interface area density model was applied to predict spray spreading angle. A method to identify spray plume boundary from the predicted flow field that offers a meaningful comparison to experimental definition was discussed in detail. Using the experimental data available in literature, a comparison between axi-symmetric and asymmetric nozzles was made to assess the nature of the influence of the nozzle geometry and the operating conditions on the ensuing spray. On the basis of this comparison, it was inferred that the spray plume angle of asymmetric nozzles is largely geometry dependent for a wide range of pressure ratios. (C) 2018 Elsevier Ltd. All rights reserved.
The collapse or merging of individual plumes of direct-injection gasoline injectors is of fundamental importance to engine performance because of its impact on fuel-air mixing. However, the mechanisms of spray collapse are not fully understood. The purpose of this work is to study the effects of injection duration and multiple injections on the interaction and/or collapse of multiplume gasoline direct injection sprays. High-speed (100 kHz) particle image velocimetry is applied along a plane between plumes to observe the full temporal evolution of plume interaction and potential collapse, resolved for individual injection events. Supporting information along a line of sight is obtained using diffused back illumination. Experiments are performed under simulated engine conditions using a symmetric 8-hole injector in a high-temperature, high-pressure vessel at the "Spray G" operating conditions of the Engine Combustion Network. Longer injection duration is found to promote plume collapse, while staging fuel delivery with multiple, shorter injections is resistant to plume collapse.
Gasoline direct injection (GDI) sprays are complex multiphase flows. When compared to multi-hole diesel sprays, the plumes are closely spaced, and the sprays are more likely to interact. The effects of multi-jet interaction on entrainment and spray targeting can be influenced by small variations in the mass fluxes from the holes, which in turn depend on transients in the needle movement and small-scale details of the internal geometry. In this paper, we present a comprehensive overview of a multi-institutional effort to experimentally characterize the internal geometry and near-nozzle flow of the Engine Combustion Network (ECN) Spray G gasoline injector. In order to develop a complete pictitre of the near-nozzle flow, a standardized setup was shared between facilities. A wide range of techniques were employed, including both X-ray and visible-light diagnostics. The novel aspects of this work include both new experimental measurements, and a comparison of the results across different techniques and facilities. The breadth and depth of the data reveal phenomena which were not apparent from analysis of the individual data sets. We show that plume-to-plume variations in the mass fluxes from the holes can cause large-scale asymmetries in the entrainment field and spray structure. Both internal flow transients and small-scale geometric features can have an effect on the external flow. The sharp turning angle of the flow into the holes also causes an inward vectoring of the plumes relative to the hole drill angle, which increases with time due to entrainment of gas into a low-pressure region between the plumes. These factors increase the likelihood of spray collapse with longer injection durations.
A computational study was performed to investigate the influence of transient needle motion on gasoline direct injection (GDI) internal nozzle flow and near-field sprays. Simulations were conducted with a compressible Eulerian flow solver modeling liquid, vapor, and non-condensable gas phases with a diffuse interface. Variable rate generation and condensation of fuel vapor were captured using the homogeneous relaxation model (HRM). The non-flashing (spray G) and flashing (spray G2) conditions specified by the Engine Combustion Network were modeled using the nominal spray G nozzle geometry. Transient needle lift and wobble were based upon ensemble averaged X-ray imaging preformed at Argonne National Lab. The minimum needle lift simulated was 5 pm and dynamic mesh motion was achieved with Laplacian smoothing. The results were qualitatively validated against experimental imaging and the experimental rate of injection profile was captured accurately using pressure boundary conditions and needle motion to actuate the injection. Low needle lift is shown to result in vapor generation near the injector seat. Finally, the internal injector flow is shown to be highly complex, containing many transient and interacting vortices which result in perturbations in the spray angle and fluctuations in the mass flux. This complex internal flow also results in intermittent string flash-boiling when a strong vortex is injected and the resulting swirling spray contains a thermal non-equilibrium vapor core. (C) 2016 Elsevier Ltd. All rights reserved.
Sprays are among the most intellectually challenging and practically important topics in fluid mechanics. This paper reviews needs, milestones, challenges, and a broad array of techniques for spray measurement. In addition, tabular summaries provide cross-referenced entry points to the vast literature by organizing over 300 citations according to key spray phenomena, physical parameters and measurement techniques for each of the principal spray regions (nozzle internal flow, near-field spray-formation region, far-field developed spray, and spray-wall interaction). The article closes with perspectives on some current issues in spray research, including the cost and complexity of apparatus for spray physics and spray engineering, the need for simultaneous diagnostic measurements under application-relevant conditions, and the effective comparison of spray measurements and numerical simulations.
Spark-ignition direct-injection engines operating in a stratified, lean-burn regime offer improved engine efficiency; however, seemingly random fluctuations in stratified combustion that result in partial-burn or misfire prevent widespread implementation. Eliminating these poor combustion events requires detailed understanding of engine flow, fuel delivery, and ignition, but knowing the dominant cause is difficult because they occur simultaneously in an engine. This study investigated the variability in fuel–air mixture linked to fuel injection hardware in a near-quiescent pressure vessel at high-temperature conditions representative of late, stratified-charge injection. An eight-hole spark-ignition direct-injection spray was interrogated using high-speed schlieren and Mie-scatter imaging from multiple, simultaneous views to acquire the vapor and liquid envelopes of the spray. The mixture fraction of vaporized sections of the spray was then quantified at a plane between plumes using Rayleigh scattering. P...
A high-speed imaging system capable of acquiring elastic scattering images and planar laser-induced fluorescence (PLIF) images in a near-simultaneous fashion has been developed. Acquiring both elastic scattering and PLIF images enables the liquid phase to be discriminated from the vapor phase. High-speed imaging allows the temporal evolution of flow structures to be evaluated. Images of sprays from a multi-hole diesel fuel injector operating under engine-like conditions were acquired. The vapor phase images reveal intricate fluid dynamic structures that exhibit a high degree of variability, indicative of a turbulent gas jet.
A novel high-speed imaging system capable of acquiring schlieren and elastic scattering images in a near-simultaneous fashion along the same line of sight has been developed. Liquid and vapor envelopes of sprays from a multi-hole gasoline fuel injector operating under engine-like conditions were systematically investigated. For the conditions evaluated, the maximum liquid penetration length was found to be inversely proportional to injection pressure. Penetration lengths derived from the elastic scattering envelops and schlieren envelopes were found to be quite similar up until the occurrence of the peak liquid penetration, implying that little vapor exists ahead of the penetrating liquid.