Search Results

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.

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.

Cycle-To-Cycle Variability (CCV) must be properly considered when modeling the ignition process in SI engines operating with ultra-lean mixtures. In this work, a strategy to model the impact of the ignition type on the CCV was developed using the RANS approach for turbulence modelling, performing multi-cycle simulations for the power-cycle only. The spark-discharge was modelled through a set of Lagrangian particles, introduced along the sparkgap and interacting with the surrounding Eulerian gas flow. Then, at each discharge event, the velocity of each particle was modified with a zero-divergence perturbation of the velocity field with respect to average conditions. Finally, the particles velocity was evolved according to the Simplified Langevin Model (SLM), which keeps memory of the initial perturbation and applies a Wiener process to simulate the stochastic interaction of each channel particle with the surrounding gas flow.

The Reynolds stress turbulence model (RSM) has been developed to go beyond the Boussinesq hypothesis and to improve turbulence modeling of flows with significant mean streamline curvature and secondary flow. In this paper the RSM in commercial CFD software CONVERGE is tested for its performance and robustness when applying to complex flows. Several validation cases including flow over flat plate, vortex combustor, diesel engine spray and combustion were selected to test the RSM. The swirling flow in vortex combustor, non-reacting but vaporizing ECN Spray A (free jet) and Sandia small bore diesel engine case are used to demonstrate the benefits of the RSM over the widely used RNG k-epsilon model without model tuning. The vortex combustor case shows the RSM can provide good prediction for strong swirling flow. ECN spray A case was used to demonstrate that the RSM can accurately predict the liquid and vapor penetration lengths of a free jet under diesel engine conditions.

Many studies have been carried out to optimize the aerodynamic performances of a single car or a single vehicle. In present days the traffic increases and sophisticated technologies are developing to guarantee the drivers safety, to minimize the fuel consumption and be more environmentally friendly. Within this research area a new technique that is being studied is Platooning: this means that different vehicles travel in a configuration that minimizes the aerodynamic drag and therefore the fuel consumption and the longitudinal space. In the present study platoons with different vehicles and configurations are taken into account, to analyze the influence of car shape and relative distance between the vehicles. The research has been carried out using CFD techniques to investigate the different flow fields around different platoons, while wind tunnel tests have been used to validate the results of the CFD simulations.

In the paper, a procedure to assess the quality of the shielding effect of wind barriers to protect large sided vehicles crossing the wake of a bridge pylon under cross wind conditions is proposed. The methodology is based on Multi-Body simulations of the response of the vehicle-driver system when it is subjected to the sudden change of the aerodynamic forces due to the wind-tower interaction. The aerodynamic forces that are instantaneously acting on the vehicle are computed according to a force distribution approach that relies on wind tunnel tests that may be performed on still scaled models. From the knowledge of the aerodynamic force distribution along the vehicle at different yaw angles and of the mean wind profile across the tower wake, the aerodynamic force, acting on the moving vehicle, is reconstructed at each time step taking into consideration the actual vehicle-driver dynamics.

The primary objective of the research discussed here was to compare the commercial computational fluid dynamics (CFD) software, CONVERGE, and a prevalent open-source code, OpenFOAM, with regard to their ability to predict spray and combustion characteristics. The high-fidelity data were obtained from the Engine Combustion Network (ECN) at Sandia National Laboratory in a constant-volume combustion vessel under well-defined, controlled conditions. The experiments and simulations were performed by using two diesel surrogate fuels (i.e., n-heptane and n-dodecane) under both non-reacting and reacting conditions. Specifically, ECN data on spray penetration, liquid length, vapor penetration, mixture fraction, ignition delay, and flame lift-off length (LOL) were used to validate both codes. Results indicate that both codes can predict the above experimental characteristics very well.

The OpenFOAM® CFD methodology is nowadays employed for simulation in internal combustion engines and a lot of work has been done for an appropriate description of all complex phenomena. At the moment in the RANS turbulence models available in the OpenFOAM® toolbox the turbulence modulation is not yet included, and the present work analyzes the predictive capabilities of the code in simulating high injection pressure fuel sprays after modeling the influence of the dispersed phase on the turbulence structure. Different experiments were employed for the validation. At first, non-evaporating diesel spray was considered in a constant volume and quiescent vessel. The validation was performed via the available experimental spray evolution in terms of penetrations and spatial/temporal fuel distributions. Then the Sandia combustion chamber was chosen for diesel spray simulation in non-reacting conditions.

Direct-injection technology represents today a very interesting solution to the typical problems that are generally encountered in SI, gas-fueled engines such as reduced volumetric efficiency, backfire and knock. However, development of suitable injection systems and combustion chamber geometry is necessary to optimize the fuel-air mixing and combustion processes. To this end, CFD models are widely applied even if the influence of the mesh structure, numerical and turbulence models on the computed results are still matter of investigation. In this work, a numerical methodology for the simulation of the gas exchange and injection processes in gas-fueled engines was developed within the Lib-ICE framework, which is a set of libraries and applications for IC engine modeling developed using the OpenFOAM® technology. The gas exchange and fuel injection processes were simulated into a four-valve, pent-roof hydrogen-fueled engine with optical access.

A numerical study has been carried out to assess the effects of needle movement and internal nozzle flow on spray formation for a multi-hole Gasoline Direct Injection system. The coupling of nozzle flow and spray formation is dynamic in nature and simulations with pragmatic choice of spatial and temporal resolutions are needed to analyze the sprays in a GDI system. The dynamic coupling of nozzle flow and spray formation will be performed using an Eulerian-Lagrangian Spray Atomization (ELSA) approach. In this approach, the liquid fuel will remain in the Eulerian framework while exiting the nozzle, while, depending on local instantaneous liquid concentration in a given cell and amount of liquid in the neighboring cells, part of the liquid mass will be transferred to the Lagrangian framework in the form of Lagrangian parcels.

The Flamelet Generated Manifold (FGM) method is a promising technique in engine combustion modeling to include tabulated chemistry. Different methodologies can be used for the generation of the manifold. Two approaches, based on igniting counterflow diffusion flamelets (ICDF) and homogeneous reactors (HR) are implemented and compared with Engine Combustion Network (ECN) experimental database for the baseline n-heptane case. Before analyzing the combustion results, the spray model is optimized after performing a sensitivity study with respect to turbulence models, cell sizes and time steps. The standard High Reynolds (Re) k-ε model leads to the best match of all turbulence models with the experimental data. For the convergence of the mixture fraction field an appropriate cell size is found to be smaller than that for an adequate spray penetration length which appears to be less influenced by the cell size.

A set of reduced chemical mechanisms was developed for use in multi-dimensional engine simulations of premixed gasoline combustion. The detailed Primary Reference Fuel (PRF) mechanism (1034 species, 4236 reactions) from Lawrence Livermore National Laboratory (LLNL) was employed as the starting mechanism. The detailed mechanism, referred to here as LLNL-PRF, was reduced using a technique known as Parallel Direct Relation Graph with Error Propagation and Sensitivity Analysis. This technique allows for efficient mechanism reduction by parallelizing the ignition delay calculations used in the reduction process. The reduction was performed for a temperature range of 800 to 1500 K and equivalence ratios of 0.5 to 1.5. The pressure range of interest was 0.75 bar to 40 bar, as dictated by the wide range in spark timing cylinder pressures for the various cases. In order to keep the mechanisms relatively small, two reductions were performed.

The continuous pursuit of higher combustion efficiencies, as well as the possible usage of synthetic fuels with different properties than fossil-ones, require reliable and low-cost numerical approaches to support and speed-up engines industrial design. In this context, SI engines operated with homogeneous ultra-lean mixtures both characterized by a classical ignition configuration or equipped with an active prechamber represent the most promising solutions. In this work, for the classical ignition arrangement, a 3DCFD strategy to model the impact of the ignition system type on the CCV is developed using the RANS approach for turbulence modelling. The spark-discharge is modelled through a set of Lagrangian particles, whose velocity is modified with a zero-divergence perturbation at each discharge event, then evolved according to the Simplified Langevin Model (SLM) to simulate stochastic interactions with the surrounding gas flow.

This study examines the impact of slip in aero-thermal conditions of supercooled large droplets (SLD) produced in an Icing Wind Tunnel (IWT) on the ice accretion characteristics. The study identifies potential biases in the SLD model development based on IWT data and numerical predictions that assume the SLD to be in aerothermal equilibrium with the IWT airflow. To obtain realistic temperature and velocity data for each droplet size class in the test section of the Braunschweig Icing Wind Tunnel (BIWT), a Lagrangian droplet tracking solver was used within a Monte Carlo framework. Results showed that SLDs experience considerable slips in velocity and temperature due to their higher inertia and short residence time in the Braunschweig IWT. Large droplets were found to be warmer and slower than the flow in the test section, with larger droplets experiencing larger aerothermal slips.

In the recent years, the interest in heavy-duty engines fueled with Compressed Natural Gas (CNG) is increasing due to the necessity to comply with the stringent CO2 limitation imposed by national and international regulations. Indeed, the reduced number of carbon atoms of the NG molecule allows to reduce the CO2 emissions compared to a conventional fuel. The possibility to produce synthetic methane from renewable energy sources, or bio-methane from agricultural biomass and/or animal waste, contributes to support the switch from conventional liquid fuels to CNG. To drive the engine development and reduce the time-to-market, the employment of numerical analysis is mandatory. This requires a continuous improvement of the simulation models toward real predictive analyses able to reduce the experimental R&D efforts. In this framework, 1D numerical codes are fundamental tools for system design, energy management optimization, and so on.

Ducted Fuel Injection (DFI) has the potential to reduce soot emissions in Diesel engines thanks to the enhanced mixing rate resulting from the liquid fuel flow through a small cylindrical pipe located at a certain distance from the nozzle injector hole. A consolidated set of experiments in constant-volume vessel and engine allowed to understand the effects of ambient conditions, duct geometry and shape on fuel-air mixing, combustion and soot formation. However, implementation of this promising technology in compression-ignition engines requires predictive numerical models that can properly support the design of combustion systems in a wide range of operating conditions. This work presents a computational methodology to predict fuel-air mixing and combustion with ducted fuel injection. Attention is mainly focused on turbulence and combustion modelling.

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.

Turbulence modeling for fuel spray simulation plays a prominent role in the understanding of the flow behavior in Internal Combustion Engines (ICEs). Currently, a lot of research work is actively spent on Large Eddy Simulation (LES) turbulence modeling as a replacement option of standard Reynolds averaged approaches in the Eulerian-Lagrangian spray modeling framework, due to its capability to accurately describe flow-induced spray variability and to the lower dependence of the results on the specific turbulence model and/or modeling coefficients. The introduction of LES poses, however, additional questions related to the implementation/adaptation of spray-related turbulence sources and to the rise of conflicting numerics and grid requirements between the Lagrangian and Eulerian parts of the simulated flow.

Homogeneous Charge Compression Ignition (HCCI) combustion technology has demonstrated a profound potential to decrease both emissions and fuel consumption. In this way, the significance of the 2-stroke HCCI engine has been underestimated as it can provide more power stroke in comparison to a 4-stroke engine. Moreover, the mass of trapped residual gases is much larger in a 2-stroke engine, causing higher initial charge temperatures, which leads to easier auto-ignition. For controlling 2-stroke HCCI engines, it is vital to find optimized simulation approaches of HCCI combustion with a focus on ignition timing. In this study, a Computational Fluid Dynamic (CFD) model for a 2-stroke gasoline engine was developed coupled to a semi-detailed chemical mechanism of iso-octane to investigate the simulation capability of the considered chemical mechanism and the effects of different simulation parameters such as the turbulence model, grid density and time step size.

The paper focuses on the development of a mesh moving method based on non-conformal topologically changing grids applied to the simulation of IC engines, where the prescribed motion of piston and valves is accomplished by rigidly translating the sub-domain representing the moving component. With respect to authors previous work, a more robust and efficient algorithm to handle the connectivity of non-conformal interfaces and a mesh-motion solver supporting multiple layer addition/removal of cells, to decouple the time-step constraints of the mesh motion and of the fluid dynamics, has been implemented as a C++ library to extend the already existing classes for dynamic mesh handling of the finite-volume, open-source CFD code OpenFOAM®. Other new features include automatic decomposition of large multiple region domains to preserve processors load balance with topological changes for parallel computations and additional tools for automatic preprocessing and case setup.