Search Results

For a series plug-in hybrid electric vehicle (PHEV), it is critical that batteries be sized to maximize vehicle performance variables, such as fuel efficiency, gasoline savings, and zero emission capability. The wide range of design choices and the cost of prototype vehicles calls for a development process to quickly and systematically determine the design characteristics of the battery pack, including its size, and vehicle-level control parameters that maximize the net present value (NPV) of a vehicle during the planning stage. Argonne National Laboratory has developed Autonomie, a modeling and simulation framework. With support from The MathWorks, Argonne has integrated an optimization algorithm and parallel computing tools to enable the aforementioned development process. This paper presents a study that utilized the development process, where the NPV is the present value of all the future expenses and savings associated with the vehicle.

Plug-in hybrid electric vehicle (PHEV) technologies have the potential for considerable petroleum consumption reductions, possibly at the expense of increased tailpipe emissions due to multiple “cold” start events and improper use of the engine for PHEV specific operation. PHEVs operate predominantly as electric vehicles (EVs) with intermittent assist from the engine during high power demands. As a consequence, the engine can be subjected to multiple cold start events. These cold start events may have a significant impact on the tailpipe emissions due to degraded catalyst performance and starting the engine under less than ideal conditions. On current hybrid electric vehicles (HEVs), the first cold start of the engine dictates whether or not the vehicle will pass federal emissions tests. PHEV operation compounds this problem due to infrequent, multiple engine cold starts.

In the development of an FAW SUV, one of the goals is to achieve a state of the art drag level. In order to achieve such an aggressive target, feedback from aerodynamics has to be included in the early stage of the design decision process. The aerodynamic performance evaluation and improvement is mostly based on CFD simulation in combination with some wind tunnel testing for verification of the simulation results. As a first step in this process, a fully detailed simulation model is built. The styling surface is combined with engine room and underbody detailed geometry from a similar size existing vehicle. From a detailed analysis of the flow field potential areas for improvement are identified and five design parameters for modifying overall shape features of the upper body are derived. In a second step, a response surface method involving design of experiments and adaptive sampling techniques are applied for characterizing the effects of the design changes.

Coatings have the potential to improve bearing tribological performance. However, every coating application process and material combination may create different residual stresses and coating microstructures, and their effect on bearing fatigue and wear performance is unclear. The aim of this work is to investigate coating induced residual stress effects on bearing failure indicators using a microstructural contact mechanics (MSCM) finite element (FE) model. The MSCM FE model consists of a two-dimensional FE model of a coated bearing surface under sliding contact where individual grains are represented by FE domains. Interactions between FE domains are represented using contact element pairs. Unique to this layered rolling contact FE model is the use of polycrystalline material models to represent realistic bearing and coating microstructural behavior. The MSCM FE model was compared to a second non-microstructural contact mechanics (non-MSCM) model.

A full understanding and characterization of the near-field of diesel sprays is daunting because the dense spray region inhibits most diagnostics. While x-ray diagnostics permit quantification of fuel mass along a line of sight, most laboratories necessarily use simple lighting to characterize the spray spreading angle, using it as an input for CFD modeling, for example. Questions arise as to what is meant by the “boundary” of the spray since liquid fuel concentration is not easily quantified in optical imaging. In this study we seek to establish a relationship between spray boundary obtained via optical diffused backlighting and the fuel concentration derived from tomographic reconstruction of x-ray radiography. Measurements are repeated in different facilities at the same specified operating conditions on the “Spray A” fuel injector of the Engine Combustion Network, which has a nozzle diameter of 90 μm.

Cavitation plays a significant role in high pressure diesel injectors. However, cavitation is difficult to measure under realistic conditions. X-ray phase contrast imaging has been used in the past to study the internal geometry of fuel injectors and the structure of diesel sprays. In this paper we extend the technique to make in-situ measurements of cavitation inside unmodified diesel injectors at pressures of up to 1200 bar through the steel nozzle wall. A cerium contrast agent was added to a diesel surrogate, and the changes in x-ray intensity caused by changes in the fluid density due to cavitation were measured. Without the need to modify the injector for optical access, realistic injection and ambient pressures can be obtained and the effects of realistic nozzle geometries can be investigated. A range of single and multi-hole injectors were studied, both sharp-edged and hydro-ground. Cavitation was observed to increase with higher rail pressures.

This paper implements a coupled approach to integrate the internal nozzle flow and the ensuing fuel spray using a Volume-of-Fluid (VOF) method in the CONVERGE CFD software. A VOF method was used to model the internal nozzle two-phase flow with a cavitation description closed by the homogeneous relaxation model of Bilicki and Kestin [1]. An Eulerian single velocity field approach by Vallet et al. [2] was implemented for near-nozzle spray modeling. This Eulerian approach considers the liquid and gas phases as a complex mixture with a highly variable density to describe near nozzle dense sprays. The mean density is obtained from the Favreaveraged liquid mass fraction. The liquid mass fraction is transported with a model for the turbulent liquid diffusion flux into the gas.

Shock waves have been recently observed in high-pressure diesel sprays. In this paper, three-dimensional numerical simulations of supersonic diesel spray injection have been performed to investigate the underlying dynamics of the induced shock waves and their interactions with the spray. A Volume-of-Fluid based method in the CFD software (CONVERGE) is used to model this multiphase phenomena. An adaptive Mesh Refinement (AMR) scheme is employed to capture the front of the spray and the shock waves with high fidelity. Simulation results are compared to the available experimental observations to validate the numerical procedure. Parametric studies with different injection and ambient conditions are conducted to examine the effect of these factors on the generation of shock waves and their dynamics.

It is widely understood that cold ambient temperatures increase vehicle fuel consumption due to heat transfer losses, increased friction (increased viscosity lubricants), and enrichment strategies (accelerated catalyst heating). However, relatively little effort has been dedicated to thoroughly quantifying these impacts across a large set of real world drive cycle data and ambient conditions. This work leverages experimental dynamometer vehicle data collected under various drive cycles and ambient conditions to develop a simplified modeling framework for quantifying thermal effects on vehicle energy consumption. These models are applied over a wide array of real-world usage profiles and typical meteorological data to develop estimates of in-use fuel economy. The paper concludes with a discussion of how this integrated testing/modeling approach may be applied to quantify real-world, off-cycle fuel economy benefits of various technologies.

Using the total probability theorem, we propose a method to calculate the failure rate of a linear vibratory system with random parameters excited by stationary Gaussian processes. The response of such a system is non-stationary because of the randomness of the input parameters. A space-filling design, such as optimal symmetric Latin hypercube sampling or maximin, is first used to sample the input parameter space. For each design point, the output process is stationary and Gaussian. We present two approaches to calculate the corresponding conditional probability of failure. A Kriging metamodel is then created between the input parameters and the output conditional probabilities allowing us to estimate the conditional probabilities for any set of input parameters. The total probability theorem is finally applied to calculate the time-dependent probability of failure and the failure rate of the dynamic system. The proposed method is demonstrated using a vibratory system.

When utilizing large models containing numerous uncertain parameters, model calibration becomes a critical step in the analysis. Traditional methods of calibration involve adjusting uncertain parameters based on expert opinion or best estimates. While this traditional calibration may lead to better model predictions, it usually only yields better estimates for certain specific conditions. This drastically reduces the functionality of the model in question. Bayesian calibration is an alternative to traditional calibration methods which utilizes available information (simulation results and/or real world measured values) to iteratively refine uncertain parameters (either assumed or measured uncertainty) while considering not only parametric uncertainty, but also model, observational, and residual uncertainties at every step of the calibration process.

In order to ensure the safety of a structure, adequate strength for structural elements must be provided. Moreover, catastrophic deformations such as buckling must be prevented. Using the linear finite element method, deterministic buckling analysis is completed in two main steps. First, a static analysis is performed using an arbitrary ordinate applied loading pattern. Using the obtained element axial forces, the geometric stiffness of the structure is assembled. Second, an eigenvalue problem is performed between structure's elastic and geometric stiffness matrices, yielding the structure's critical buckling loads. However, these deterministic approaches do not consider uncertainty the structure's material and geometric properties. In this work, a new method for finite element based buckling analysis of a structure with uncertainty is developed. An imprecise probability formulation is used to quantify the uncertainty present in the mechanical characteristics of the structure.

A twist beam rear suspension system is modeled, analyzed and optimized in this paper. An ADAMS model is established based on the REC (Rigid-Elastic Coupling) Theory, which is verified by FEM (Finite Element Method) approach, the effects of the geometric parameters on the twist beam suspension performance are investigated. In order to increase the calculation efficiency and improve the simulation accuracy, a neural network model and NSGA II (Non-dominated Sorting Genetic Algorithm II) are adopted to conduct a multi-objective optimization on a twist beam rear suspension system.

Multi-mode combustion strategies may provide a promising pathway to improve thermal efficiency in light-duty spark ignition (SI) engines by enabling switchable combustion modes, wherein an engine may operate under advanced compression ignition (ACI) at low load and spark-assisted ignition at high load. The extension from the SI mode to the ACI mode requires accurate control of intake charge conditions, e.g., pressure, temperature and equivalence ratio, in order to achieve stable combustion phasing and rapid mode-switches. This study presents results from computational fluid dynamics (CFD) analysis to gain insights into mixture charge formation and combustion dynamics pertaining to auto-ignition processes. The computational study begins with a discussion of thermal wall boundary condition that significantly impacts the combustion phasing.

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.

This paper reports investigations on diesel jet transients, accounting for internal nozzle flow and needle motion. The calculations are performed with Large Eddy Simulation (LES) turbulence model by coupling the internal and external multiphase flows simultaneously. Short and multiple injection strategies are commonly used in internal combustion engines. Their features are significantly different from those generally found in steady state conditions, which have been extensively studied in the past, however, these conditions are seldom reached in modern engines. Recent researches have shown that residual gas can be ingested in the injector sac after the end-of-injection (EOI) and undesired dribbles can be produced. Moreover, a new injection event behaves differently at the start-of-injection (SOI) depending on the sac initial condition, and the initial spray development can be affected for the first few tens of μs.

Quantitative measurements of direct injection fuel spray density and mixing are difficult to achieve using optical diagnostics, due to the substantial scattering of light and high optical density of the droplet field. For multi-hole sprays, the problem is even more challenging, as it is difficult to isolate a single spray plume along a single line of sight. Time resolved x-ray radiography diagnostics developed at Argonne's Advanced Photon Source have been used for some time to study diesel fuel sprays, as x-rays have high penetrating power in sprays and scatter only weakly. Traditionally, radiography measurements have been conducted along any single line of sight, and have been applied to single-hole and group-hole nozzles with few plumes. In this new work, we extend the technique to multi-hole gasoline direct injection sprays.

A computational fluid dynamics (CFD) guided combustion system optimization was conducted for a heavy-duty compression-ignition engine with a gasoline-like fuel that has an anti-knock index (AKI) of 58. The primary goal was to design an optimized combustion system utilizing the high volatility and low sooting tendency of the fuel for improved fuel efficiency with minimal hardware modifications to the engine. The CFD model predictions were first validated against experimental results generated using the stock engine hardware. A comprehensive design of experiments (DoE) study was performed at different operating conditions on a world-leading supercomputer, MIRA at Argonne National Laboratory, to accelerate the development of an optimized fuel-efficiency focused design while maintaining the engine-out NOx and soot emissions levels of the baseline production engine.

The interaction of fuel sprays and in-cylinder flow in direct-injection engines is expected to alter kinetic energy and integral length scales at least during some portions of the engine cycle. High-speed particle image velocimetry was implemented in an optical four-valve, pent-roof spark-ignition direct-injection single-cylinder engine to quantify this effect. Non-firing motored engine tests were performed at 1300 RPM with and without fuel injection. Two fuel injection timings were investigated: injection in early intake stroke represents quasi-homogenous engine condition; and injection in mid compression stroke mimics the stratified combustion strategy. Two-dimensional crank angle resolved velocity fields were measured to examine the kinetic energy and integral length scale through critical portions of the engine cycle. Reynolds decomposition was applied on the obtained engine flow fields to extract the fluctuations as an indicator for the turbulent flow.

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