Search Results

Obtaining and maintaining a stratified charge in a practical engine is a difficult problem. Consequently, many approaches have been proposed and reported in the scientific and patent literature. In attempting to assess the most profitable approach for future development work, it is important to group together similar approaches so that one can study their performance as a group. Making such a classification has the additional advantage of helping to standardize terminology used by different investigators. With this thought in mind, a literature study was made and a proposed classification chart prepared for the different engine combustion systems reported in the literature. For the sake of completeness, the finally proposed classification chart includes homogeneous combustion engines as well as heterogeneous combustion engines. Because of their similarity of combustion, rotary engines such as the Wankel engine are considered as “reciprocating” although gas turbines are not included.

Due to growing interest in hybrid and electric vehicles, li-ion battery modeling is receiving a lot of attention from designers and researchers. This paper presents a complete model for a li-ion battery pack. It starts from the Newman electrochemistry model to create the battery performance curves. Such information is then used for cell level battery equivalent circuit model (ECM) parameter identification. 28 cell ECMs are connected to create the module ECM. Four module ECMs are connected through a busbar model to create the pack ECM. The busbar model is a reduced order model (ROM) extracted from electromagnetic finite element analysis (FEA) results, taking into account the parasitic effects. Battery thermal performance is simulated first by computational fluid dynamics (CFD). Then, a thermal linear and time-invariant (LTI) ROM is created out of CFD solution. The thermal LTI ROM is then two-way coupled with the battery pack ECM to form a complete battery pack model.

A computer program which rapidly calculates the equilibrium mole fractions and the partial derivatives of the mole fractions with respect to temperature, pressure and equivalence ratio for the products of combustion of any hydrocarbon fuel and air is described. A subroutine is also given which calculates the gas constant, enthalpy, internal energy and the partial derivatives of these with respect to temperature, pressure and equivalence ratio. Some examples of the uses of the programs are also given.

Demands for higher power engines have led to higher pressures in fuel injectors. Internal nozzle flow plays a critical role in the near nozzle flow and subsequent spray pattern. The internal flow becomes more difficult to model when the injector pressure and internal shape make it more prone to cavitation. Two Bosch injectors, proposed for experimental and computational studies under the Engine Combustion Network (namely “Spray C” and “Spray D”) are modeled in the computational fluid dynamics code ANSYS Fluent. Both injectors operate with n-dodecane as fuel at 150 MPa inlet pressures. The computational model includes cavitation effects to characterize any cavitating regions. Including compressibility of both liquid and vapor is found to be critical. Also, due to high velocity gradients and stresses in the nozzle, turbulent viscous energy dissipation is considered along with pressure work resulting from significant pressure changes in the injector.

A linear parameter varying (LPV) reduced order model (ROM) is used to approximate the volume-averaged temperature of battery cells in one of the modules of the battery pack with varying mass flow rate of cooling fluid using uniform heat source as inputs. The ROM runs orders of magnitude faster than the original CFD model. To reduce the time it takes to generate training data, used in building LPV ROM, a divide-and-conquer approach is introduced. This is done by dividing the battery module into a series of mid-cell and end-cell units. A mid-cell unit is composed of a cooling channel sandwiched in between two half -cells. A half-cell has half as much heat capacity as a full-cell. An end-cell unit is composed of a cooling channel sandwiched in between full-cell and a half-cell. A mass flow rate distribution look-up-table is generated from a set of steady-state simulations obtained by running the full CFD model at different inlet manifold mass flow rate samples.

The results of two derivations relating to the fuel economy of hybrid-electric vehicles (vehicles which employ both a heat engine and electric drive system) are presented and their use is illustrated through the examples of the University of Wisconsin and TRW Systems Group hybrid-electric vehicles. The method of mileage estimation employs a specific fuel-consumption versus torque-speed map for the heat engine under study and knowledge of the hybrid-vehicle dynamics and road-load power. The method is easily extended to estimation of emission reductions through use of specific-emission-production versus torque-speed maps and is applicable to hybrid vehicles with other than electrical energy-storage systems.

Due to growing interest in hybrid and electric vehicles, the battery, being one of the critical components, is receiving a lot of attention from designers and researchers. Two battery-modeling approaches, though seemingly different, share the same mathematical challenge of robust non-linear curve-fitting. The two methods are battery equivalent circuit model and battery system level thermal modeling using the linear time-invariant (LTI) method. Both modeling approaches involve curve-fitting testing data or data from advanced models to identify four parameters in a circuit model consisting of two pairs of RC elements. Such curve-fitting is mathematically a non-linear least-squares (LS) problem. Standard methods like the Levenberg-Marquardt (LM) method can be used for non-linear curve-fitting, but the LM method is known to be sensitive to initial conditions.

This paper introduces a model-based systems and embedded software engineering, workflow for the design of control systems. The interdisciplinary approach that is presented relies on an integrated set of tools that addresses the needs of various engineering groups, including system architecture, design, and validation. For each of these groups, a set of best practices has been established and targeted tools are proposed and integrated in a unique platform, thus allowing efficient communication between the various groups. In the initial stages of system design, including functional and architectural design, a SysML-based approach is proposed. This solution is the basis to develop systems that have to obey both functional and certification standards such as ARINC 653 (IMA) and ARP 4754A. Detailed system design typically requires modeling and simulation of each individual physical component of the system by various engineering groups (mechanical, electrical, etc.).

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

Fine-wire resistance thermometers were used to measure compression gas temperatures in a motoring (nonfiring) cycle CFR engine. Temperature versus crankangle curves were obtained for the compression and expansion strokes by means of tungsten wires ranging in diameter from 0.15–1.00 mils and at speeds from 600–1800 rpm. The results were compared with the infrared pyrometer; the peak temperature and peak crankangle lags were determined as a function of the wire diameter and engine speed. Attempts to evaluate the instantaneous energy balance around the wire resulted in a negative heat transfer coefficient, for which no current satisfactory explanation is available, although other observers have reported similar phenomena. The tungsten resistance thermometer is simple to build, easy to install, and requires no modification of the engine block for use during motoring. Thus, it is suitable for comparing the compression temperatures of different design engines.

Battery thermal management for high power applications such as electrical/hybrid vehicles is crucial. Modeling is an indispensable tool to help engineers design better battery cooling systems. While Computational Fluid Dynamics (CFD) has been used quite successfully for battery thermal management, CFD models can be too large and too slow for repeated transient thermal analysis especially for a battery module or pack. An accurate but much smaller battery thermal model using a state space representation is proposed. The parameters in the state space model are extracted from CFD results. The state space model is then shown to provide identical results as those from CFD under transient power inputs. While a CFD model may take hours to run depending on the size of the problem, the corresponding state space model runs in seconds.

Cycle-to-cycle variation in combustion phasing and combustion rate cause knock to occur differently in every cycle. This is found to be true even if the end gas thermo-chemical time history is the same. Three cycles are shown that have matched combustion phasing, combustion rate, and time of knock onset, but have knock intensity that differs by a factor of six. Thus, the prediction of knock intensity must include a stochastic component. It is shown that there is a relationship between the maximum possible knock intensity and the unburned fuel energy at the time of knock onset. Further, for a small window of unburned energy at knock onset, the probability density function of knock intensity is self similar when scaled by the 95th percentile of the cumulative distribution, and log-normal in shape.

The development of a high speed, multichannel data acquisition system is described. A precision magnetic tape recorder is used to record analog data from highly transient phenomena. Analog-to-digital data conversion is performed on a hybrid computer and the digitized data is processed using large, high speed digital computers. A detailed example of the application of the system to the measurement of rates-of-injection, rates-of-heat release, and instantaneous rates-of-heat transfer from the cylinder gases to the cylinder walls in a high speed open-chamber diesel engine is presented.

This paper presents the current state of a three-layer surface icing model for ice crystal icing risk assessment in aircraft engines, being developed jointly by Ansys and Honeywell to account for possible heat transfer from inside an engine into the flow path where ice accretion occurs. The bottom layer of the proposed model represents a thin metal sheet as a substrate surface to conductively transfer heat from an engine-internal reservoir to the ice layer. The middle layer is accretion ice with a porous structure able to hold a certain amount of liquid water. A shallow water film layer on the top receives impinged ice crystals. A mass and energy balance calculation for the film determines ice accretion rate. Water wicking and recovery is introduced to transfer liquid water between film layer and porous ice accretion layer.

Modeling ignition kernel development in spark ignition engines is crucial to capturing the sources of cyclic variability, both with RANS and LES simulations. Appropriate kernel modeling must ensure that energy transfer from the electrodes to the gas phase has the correct timing, rate and locations, until the flame surface is large enough to be represented on the mesh by the G-Equation level-set method. However, in most kernel models, geometric details driving kernel growth are missing: either because it is described as Lagrangian particles, or because its development is simplified, i.e., down to multiple spherical flames. This paper covers the geometric aspects of kernel development, which makes up the core of a Triangulated Lagrangian Ignition Kernel model. One (or multiple, if it restrikes) spark channel is initialized as a one-dimensional Lagrangian particle thread.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

An experimental investigation of non-intrusive combustion sensing was performed using a tri-axial accelerometer mounted to the engine block of a small-bore high-speed 4-cylinder compression-ignition direct-injection (CIDI) engine. This study investigates potential techniques to extract combustion features from accelerometer signals to be used for cycle-to-cycle engine control. Selection of accelerometer location and vibration axis were performed by analyzing vibration signals for three different locations along the block for all three of the accelerometer axes. A magnitude squared coherence (MSC) statistical analysis was used to select the best location and axis. Based on previous work from the literature, the vibration signal filtering was optimized, and the filtered vibration signals were analyzed. It was found that the vibration signals correlate well with the second derivative of pressure during the initial stages of combustion.

A methodology has been implemented to calculate local turbulent flame speeds for spark ignition engines accurately and on-the-fly in 3-D CFD modeling. The approach dynamically captures fuel effects, based on detailed chemistry calculations of laminar flame speeds. Accurately modeling flame propagation is critical to predicting heat release rates and emissions. Fuels used in spark ignition engines are increasingly complex, which necessitates the use of multi-component fuels or fuel surrogates for predictive simulation. Flame speeds of the individual components in these multi-component fuels may vary substantially, making it difficult to define flame speed values, especially for stratified mixtures. In addition to fuel effects, a wide range of local conditions of temperature, pressure, equivalence ratio and EGR are expected in spark ignition engines.

Due to the increasingly stringent environmental regulations all around the world confronted by exhaust emission and energy consumption, improving fuel economy has been the top priority for most automotive manufacturers. In this context, the basic process for vehicle shape development has evolved into optimizing the design to achieve better aerodynamic characteristics, especially drag reduction. Of all the optimization approaches, the gradient-based adjoint method has currently received extensive attention for its high efficiency in calculating the objective sensitivity with respect to geometry parameters, which is the first and foremost step for subsequent shape modification. In this work, the main goal is to explore the adjoint method through optimizing the vehicle shape for a lower drag based on a production SUV. Firstly, the influence of different mesh schemes was discussed on sensitivity prediction of aerodynamic drag.