Search Results

In this work an integration between a 1D code (Gasdyn) with a CFD code (OpenFOAM®) has been applied to improve the performance of a Moto3 engine. The four-stroke, single cylinder S.I. engine was modeled, in order to predict the wave motion in the intake and exhaust systems and to study how it affects the cylinder gas exchange process. The engine considered was characterized by having an air induction system with integrated filter cartridge, air-box and intake runner, including two fuel injectors, resulting in a complex air-path from the intake mouth to the intake valves, which presents critical aspects when a 1D modeling is addressed. The exhaust and intake systems have been optimized form the point of view of the wave action. However, due to the high revolution speed reached by this type of engine, the interaction between the gas stream and the fuel spray becomes a key aspect to be addressed in order to achieve the best performance at the desired operating condition.

In July 1997, the Pacific Northwest and Alaska Regional Bioenergy Program, in cooperation with several industrial and institutional partners initiated a long-haul 322,000 km (200,000 mile) operational demonstration using a biodiesel and diesel fuel blend in a 324 kW (435 HP), Caterpillar 3406E Engine, and a Kenworth Class 8 heavy duty truck. This project was designed to: develop definitive biodiesel performance information, collect emissions data for both regulated and non-regulated compounds including mutagenic activity, and collect heavy-duty operational engine performance and durability information. To assess long-term engine durability and wear; including injector, valve and port deposit formations; the engine was dismantled for inspection and evaluation at the conclusion of the demonstration. The fuel used was a 50% blend of biodiesel produced from used cooking oil (hydrogenated soy ethyl ester) and 50% 2-D petroleum diesel.

High fidelity engine simulation requires realistic fuel models. Although typical automotive fuels consist of more than few hundreds of hydrocarbon species, researches show that the physical and chemical properties of the real fuels could be represented by appropriate surrogate fuel models. It is desirable to represent the fuel using the same set of physical and chemical surrogate components. However, when the reaction mechanisms for a certain physical surrogate component is not available, the chemistry of the unmatched physical component is described using that of a similar chemical surrogate component at the expense of accuracy. In order to reduce the prediction error while maintaining the computational efficiency, a method of on-the-fly reactivity adjustment (ReAd) of chemical reaction mechanism along with fuel re-distribution based on reactivity is presented and tested in this study.

Modeling plume interaction and collapse for direct-injection gasoline sprays is important because of its impact on fuel-air mixing and engine performance. Nevertheless, the aerodynamic interaction between plumes and the complicated two-phase coupling of the evaporating spray has shown to be notoriously difficult to predict. With the availability of high-speed (100 kHz) Particle Image Velocimetry (PIV) experimental data, we compare velocity field predictions between plumes to observe the full temporal evolution leading up to plume merging and complete spray collapse. The target “Spray G” operating conditions of the Engine Combustion Network (ECN) is the focus of the work, including parametric variations in ambient gas temperature. We apply both LES and RANS spray models in different CFD platforms, outlining features of the spray that are most critical to model in order to predict the correct aerodynamics and fuel-air mixing.

The paper illustrates and validates a novel predictive combustion model for the estimation of performances and pollutant production in CI engines. The numerical methodology was developed by the authors for near real-time applications, while aiming at an accurate description of the air mixing process by means of a multi-zone approach of the air-fuel mass. Charge stratification is estimated via a 2D representation of the fuel spray distribution that is numerically derived by an axial one-dimensional control-volume description of the direct injection. The radial coordinate of each control volume is reconstructed a posteriori by means of a local distribution function. Fuel mass clustered in each zone is further split in ‘liquid’, ‘unburnt’ and ‘burnt’ sub-zones, given the local properties of the fuel spray control volumes with respect to space-time location of modelled ignition delay, liquid length, and flame lift-off.

The dynamics and evolution of turbulent structures inside an engine-like geometry are investigated by means of Large Eddy Simulation. A simplified configuration consisting of a flat-top cylinder head with a fixed, axis-centered valve and low-speed piston has been simulated by the finite volume CFD code OpenFOAM®; the standard version of the software has been extended to include the compressible WALE subgrid-scale model, models for the generation of synthetic turbulence, some improvements to the mesh motion strategy and algorithms for LES data post-processing. In order to study both the initial transient and the quasi- steady operating conditions, ten complete engine cycles have been simulated. Phase and spatial averages have been performed over cycles three to ten in order to extract first and second moment of velocity; these quantities have then been used to validate the numerical procedure by comparison against experimental data.

Spray impingement is an important phenomenon that introduces turbulence into the spray that promotes fuel vaporization, air entrainment and flame propagation. However, liquid impingement on the surface leads to wall-wetting and film deposition. The film region is a fuel-rich zone and it has potentials to produce higher emission. Film deposition in a non-reacting spray was studied previously but not in a reacting spray. In the current study, the film deposition of a reacting diesel spray was studied through computational fluid dynamic (CFD) simulations under a variety of ambient temperatures, gas compositions and impinging distances. Characteristics of film mass, distribution of thickness, soot formation and temperature distributions were investigated. Simulation results showed that under the same impinging distance, higher ambient temperature reduced film mass but showed the same liquid film pattern.

Swirling flows are very dominant in applied technical problems, especially in IC engines, and their prediction requires rather sophisticated modeling. An adaptive low-pass filtering procedure for the modeled turbulent length and time scales is derived and applied to Menter' original k - ω SST turbulence model. The modeled length and time scales are compared to what can potentially be resolved by the computational grid and time step. If the modeled scales are larger than the resolvable scales, the resolvable scales will replace the modeled scales in the formulation of the eddy viscosity; therefore, the filtering technique helps the turbulence model to adapt in accordance with the mesh resolution and the scales to capture.

This paper presents a novel fully-automatic remeshing procedure, based on the level-set method and Delaunay triangulation, to model three-dimensional boundary problems and generate a new conformal body-fitted mesh. The proposed methodology is applied to long-term in-flight ice accretion, which is characterized by the formation of extremely irregular ice shapes. Since ice accretion is coupled with the aerodynamic flow field, a multi-step procedure is implemented. The total icing exposure time is subdivided into smaller time steps, and at each time step a three-dimensional body-fitted mesh, suitable for the computation of the aerodynamic flow field around the updated geometry, is generated automatically. The methodology proposed can effectively deal with front intersections, as shown with a manufactured example.

A computational methodology has been developed for loss prediction in intake regions of internal combustion engines. The methodology consists of a hierarchy of four major tasks: (1) proper computational modeling of flow physics; (2) exact geometry and high quality and generation; (3) discretization schemes for low numerical viscosity; and (4) higher order turbulence modeling. Only when these four tasks are dealt with properly will a computational simulation yield consistently accurate results. This methodology, which is has been successfully tested and validated against benchmark quality data for a wide variety of complex 2-D and 3-D laminar and turbulent flow situations, is applied here to a loss prediction problem from industry. Total pressure losses in the intake region (inlet duct, manifold, plenum, ports, valves, and cylinder) of a Caterpillar diesel engine are predicted computationally and compared to experimental data.

Very little is currently known of the aerodynamic interaction between neighboring cycloidal rotors. Such knowledge is, however, of crucial importance to tune the controller and rotor disposition of a cyclogiro aircraft. Thus, a three-dimensional computational fluid dynamics (CFD) model is developed, validated, and used to analyze the D-Dalus L1 four-rotor unmanned aircraft operating under several configurations. The model solves the Euler equations using the OpenFOAM toolbox in order to provide fast results on a desktop computer. Validation is performed against thrust forces and flow streamlines obtained during wind tunnel experiments at various flight velocities. Numerical results from CFD match the trends of the experimental data. Flow behavior matches the video footage of the wind tunnel tests. Although boundary layer effects are neglected, satisfactory results are obtained both qualitatively and quantitatively.

Passive regeneration (oxidation of particulate matter without using an external energy source) of particulate filters in combination with active regeneration is necessary for low load engine operating conditions. For low load conditions, the exhaust gas temperatures are less than 250°C and the PM oxidation rate due to passive regeneration is less than the PM accumulation rate. The objective of this research was to experimentally investigate active regeneration of a catalyzed particulate filter (CPF) using diesel fuel injection in the exhaust gas after the turbocharger and before a diesel oxidation catalyst (DOC) and to collect data for extending the MTU 1-D 2-layer model to include the simulation of active regeneration. The engine used in this study was a 2002 Cummins ISM turbo charged 10.8 L heavy duty diesel engine with cooled EGR. The exhaust after-treatment system consisted of a Johnson Matthey DOC and CPF (a CCRT®).

The flow field resulting from injecting a gas jet into a crossflow confined in a narrow square duct has been studied under steady regime using schlieren imaging and laser Doppler velocimetry (LDV). This transparent duct is intended to simulate the intake port of an internal combustion engine fueled by gaseous mixture, and the jet is issued from a round nozzle. The schlieren images show that the relative small size of the duct would confine the development of the transverse jet, and the interaction among jet and sidewalls strongly influences the mixing process between jet and crossflow. The mean velocity and turbulence fields have been studied in detail through LDV measurements, at both center plane and several cross sections. The well-known flow feature formed by a counter rotating vortex pair (CVP) has been observed, which starts to appear at the jet exit section and persists far downstream contributing to enhancing mixing process.

The effect of flow direction towards the spark plug electrodes on ignition parameters is analyzed using an innovative spark aerodynamics fixture that enables adjustment of the spark plug gap orientation and plug axis tilt angle with respect to the incoming flow. The ignition was supplied by a long discharge high energy 110 mJ coil. The flow was supplied by compressed air and the spark was discharged into the flow at varying positions relative to the flow. The secondary ignition voltage and current were measured using a high speed (10MHz) data acquisition system, and the ignition-related metrics were calculated accordingly. Six different electrode designs were tested. These designs feature different positions of the electrode gap with respect to the flow and different shapes of the ground electrodes. The resulting ignition metrics were compared with respect to the spark plug ground strap orientation and plug axis tilt angle about the flow direction.

Fuel spray interactions with piston surfaces and cylinder walls in internal combustion engines have been extensively studied in the past decades. However, there still exists an imperative knowledge gap on the fundamental understanding of dynamic droplet-wall interactions. Particularly, the impinging angle of droplet has been barely investigated as it renders asymmetrical droplet behaviors. This paper aims to provide detailed data of droplet-inclined surface impingement physics which could further support spray-wall model development. The experimental work of single diesel droplet impinging on an inclined dry surface was conducted under isothermal (25°C) conditions. Various droplet impact angle (φ) was achieved by adjusting surface tilting angle which was set from 0° to 45° in current study. A single diesel droplet impinged onto the inclined surface with different Weber number (around 20 ~ 800).

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.

Modeling the combustion process of a diesel-biodiesel fuel spray in a 3-dimensional (3D) computational fluid dynamics (CFD) domain remains challenging and time-consuming despite the recent advancement in computing technologies. Accurate representation of the in-cylinder processes is essential for CFD studies to provide invaluable insights into these events, which are typically limited when using conventional experimental measurement techniques. This is especially true for emerging new fuels such as biodiesels since fundamental understanding of these fuels under combusting environment is still largely unknown. The reported work here is dedicated to evaluating the Adaptive Local Mesh Refinement (ALMR) approach in OpenFOAM® for improved simulation of reacting biodiesel fuel spray. An in-house model for thermo-physical and transport properties is integrated to the code, along with a chemical mechanism comprising 113 species and 399 reactions.

Wheel and wheelhouses contribute up to 20-30% of the aerodynamic drag of passenger cars. Simulating the flow field around wheels is challenging due to the complexity of the flow structures generated by tires and rims, wheel rotation, tire deformation and contact with the ground. High accuracy is usually obtained with transient simulations that treat rim rotation with the Sliding Mesh (SM) approach, which is also computationally expensive. Previous studies have confirmed that the application of a tangential velocity component to the rim surface is unphysical for open rims, while a Moving Reference Frame (MRF) is lacking accuracy and the averaged results depend on the initial spokes position. These methods do not consider the dynamic nature of the problem. This work proposes the use of the Actuator Line (AL) and Rotor Disk (RD) approaches as alternatives for simulating open rims with much lower computational cost.

A KIVA-based CFD tool has been utilized to simulate the effect of a Common-Rail injection system applied to a large, uniflow-scavenged, two-stroke diesel engine. In particular, predictions for variations of injection pressure and injection duration have been validated with experimental data. The computational models have been evaluated according to their predictive capabilities of the combustion behavior reflected by the pressure and heat release rate history and the effects on nitric oxide formation and wall temperature trends. In general, the predicted trends are in good agreement with the experimental observations, thus demonstrating the potential of CFD as a design tool for the development of large diesel engines equipped with Common-Rail injection. Existing deficiencies are identified and can be explained in terms of model limitations, specifically with respect to the description of turbulence and combustion chemistry.

The definition of a robust methodology to perform a full-cycle CFD simulation of IC engines requires as first step the availability of a reliable grid generation tool, which does not only have to guarantee a high quality mesh but also has to prove to be efficient in terms of required time. In this work the authors discuss a novel approach entirely based on the OpenFOAM technology, in which the available 3D grid generator was employed to automatically create meshes containing hexahedra and split-hexahedra from triangulated surface geometries in Stereolithography (STL) format. The possibility to introduce local refinements and boundary layers makes this tool suitable for IC engine simulations. Grids are sequentially generated at target crank angles which are automatically determined depending on user specified settings such as maximum mesh validity interval and quality parameters like non-orthogonality, skewness and aspect ratio.