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.

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.

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.

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

This work involves modeling internal and near-nozzle flows of a gasoline direct injection (GDI) nozzle. The Engine Combustion Network (ECN) Spray G condition has been considered for these simulations using the nominal geometry of the Spray G injector. First, best practices for numerical simulation of the two-phase flow evolution inside and the near-nozzle regions of the Spray G injector are presented for the peak needle lift. The mass flow rate prediction for peak needle lift was in reasonable agreement with experimental data available in the ECN database. Liquid plume targeting angle and liquid penetration estimates showed promising agreement with experimental observations. The capability to assess the influence of different thermodynamic conditions on the two-phase flow nature was established by predicting non-flashing and flashing phenomena.

In a Spark-Ignited engine, there will come a point, as load is increased, where the unburned air-fuel mixture undergoes auto-ignition (knock). The onset of knock represents the upper limit of engine output, and limits the extent of engine downsizing / boosting that can be implemented for a given application. Although effective at mitigating knock, requiring high octane fuel is not an option for most markets. Retarding spark timing can extend the high load limit incrementally, but is still bounded by limits for exhaust gas temperature, and spark retard results in a notable loss of efficiency. Likewise, enriching the air-fuel mixture also decreases efficiency, and has profound negative impacts on engine out emissions. In this current work, a Direct-Injected, Boosted, Spark-Ignited engine with Variable Valve Timing was tested under steady state high load operation. Comparisons were made among three fuels; an 87 AKI, a 91 AKI, and a 110 AKI off-road only race fuel.

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.

It is challenging to perform probabilistic analysis and design of large-scale structures because probabilistic analysis requires repeated finite element analyses of large models and each analysis is expensive. This paper presents a methodology for probabilistic analysis and reliability based design optimization of large scale structures that consists of two re-analysis methods; one for estimating the deterministic vibratory response and another for estimating the probability of the response exceeding a certain level. The deterministic re-analysis method can analyze efficiently large-scale finite element models consisting of tens or hundreds of thousand degrees of freedom and large numbers of design variables that vary in a wide range. The probabilistic re-analysis method calculates very efficiently the system reliability for many probability distributions of the design variables by performing a single Monte Carlo simulation.

An approach for Probabilistic Reanalysis (PRA) of a system is presented. PRA calculates very efficiently the system reliability or the average value of an attribute of a design for many probability distributions of the input variables, by performing a single Monte Carlo simulation. In addition, PRA calculates the sensitivity derivatives of the reliability to the parameters of the probability distributions. The approach is useful for analysis problems where reliability bounds need to be calculated because the probability distribution of the input variables is uncertain or for design problems where the design variables are random. The accuracy and efficiency of PRA is demonstrated on vibration analysis of a car and on system reliability-based optimization (RBDO) of an internal combustion engine.

In this study a dynamic simulation was conducted on a forged steel crankshaft from a single cylinder four stroke engine. Finite element analysis was performed to obtain the variation of the stress magnitude at critical locations. The dynamic analysis resulted in the development of the load spectrum applied to the crankpin bearing. This load was then applied to the FE model and boundary conditions were applied according to the engine mounting conditions. Results obtained from the aforementioned analysis were then used in optimization of the forged steel crankshaft. Geometry, material, and manufacturing processes were optimized using different geometric constraints, manufacturing feasibility, and cost. The first step in the optimization process was weight reduction of the component considering dynamic loading. This required the stress range under dynamic loading not to exceed the magnitude of the stress range in the original crankshaft.

During metal forming processes such as rolling and forging, deformable manganese sulfide (MnS) inclusions become elongated. Such elongated MnS inclusions can have considerable adverse effects on mechanical properties, if the inclusions are not aligned with the loading direction. The objectives of this study were to evaluate and compare fatigue, monotonic tensile and CVN impact behavior of SAE 4140 steel with high (0.077% S), low (0.012% S) and ultra low (0.004% S) sulfur contents at two hardness levels (40 HRC and 50 HRC). The longitudinally oriented samples at 40 HRC, where MnS inclusions were oriented along the loading direction, did not exhibit any significant sensitivity of tensile or fatigue properties to the sulfur content. For the transversely oriented MnS inclusions, however, the monotonic tensile test results indicate very low ductility of the high sulfur material at both hardness levels, where specimens failed shortly after yielding.

This paper discusses the effects of changes in specimen geometry, stress gradient, and residual stresses on fully-reversed constant amplitude uniaxial fatigue behavior of a medium carbon steel. Axial fatigue tests were performed on both flat and round specimens, while four-point rotating bending tests were performed only on round specimens. All the tests were performed using shot peened and unpeened flat and round samples, to investigate the effects of compressive residual stresses on fatigue behavior. The specimens in the rotating bending tests experienced longer life for a given stress amplitude than in the axial test. Shot-peening was found to be beneficial in the long life region, while in short life tests the shot-peened samples experienced a shorter life than the unpeened samples under both axial and bending test conditions.

A procedure for fatigue life estimation of components and structures under variable amplitude multiaxial loadings based on simple and commonly available material properties is presented. Different aspects of the analysis consisting of load cycle counting method, plasticity model, fatigue damage parameter, and cumulative damage rule are presented. The only needed material properties for the proposed procedure are hardness and monotonic and axial cyclic deformation properties (HB, K, n, K′ and n′). Rainflow cycle counting method is used for identifying number of cycles. Non-proportional cyclic hardening is estimated from monotonic and axial cyclic deformation behaviors. A critical plane approach is used to quantify fatigue damage under variable amplitude multiaxial loading, where only material hardness is used to estimate the fatigue curve, and where the needed deformation response is estimated based on Tanaka's non-proportionality parameter.

Fuels in the gasoline auto-ignition range (Research Octane Number (RON) > 60) have been demonstrated to be effective alternatives to diesel fuel in compression ignition engines. Such fuels allow more time for mixing with oxygen before combustion starts, owing to longer ignition delay. Moreover, by controlling fuel injection timing, it can be ensured that the in-cylinder mixture is “premixed enough” before combustion occurs to prevent soot formation while remaining “sufficiently inhomogeneous” in order to avoid excessive heat release rates. Gasoline compression ignition (GCI) has the potential to offer diesel-like efficiency at a lower cost and can be achieved with fuels such as low-octane straight run gasoline which require significantly less processing in the refinery compared to today’s fuels.

The internal structure of Diesel fuel injectors is known to have a significant impact on the nozzle flow and the resulting spray emerging from each hole. In this paper the three-dimensional transient flow structures inside a Diesel injector is studied under nominal (in-axis) and realistic (including off-axis lateral motion) operating conditions of the needle. Numerical simulations are performed in the commercial CFD code CONVERGE, using a two-phase flow representation based on a mixture model with Volume of Fluid (VOF) method. Moving boundaries are easily handled in the code, which uses a cut-cell Cartesian method for grid generation at run time. First, a grid sensitivity study has been performed and mesh requirements are discussed. Then the results of moving needle calculations are discussed. Realistic radial perturbations (wobbles) of the needle motion have been applied to analyze their impact on the nozzle flow characteristics.

In this study a dynamic simulation was conducted on a crankshaft from a single cylinder four stroke engine. Finite element analysis was performed to obtain the variation of stress magnitude at critical locations. The pressure-volume diagram was used to calculate the load boundary condition in dynamic simulation model, and other simulation inputs were taken from the engine specification chart. The dynamic analysis was done analytically and was verified by simulation in ADAMS which resulted in the load spectrum applied to crank pin bearing. This load was applied to the FE model in ABAQUS, and boundary conditions were applied according to the engine mounting conditions. The analysis was done for different engine speeds and as a result critical engine speed and critical region on the crankshaft were obtained. Stress variation over the engine cycle and the effect of torsional load in the analysis were investigated.

In reliability design, often, there is scarce data for constructing probabilistic models. Probabilistic models whose parameters vary in known intervals could be more suitable than Bayesian models because the former models do not require making assumptions that are not supported by the available evidence. If we use models whose parameters vary in intervals we need to calculate upper and lower bounds of the failure probability (or reliability) of a system in order to make design decisions. Monte Carlo simulation can be used for this purpose, but it is too expensive for all but very simple systems. This paper proposes an efficient Monte-Carlo simulation approach for estimation of upper and lower probabilities. This approach is based on two ideas: a) use an efficient approach for reliability reanalysis of a system, which is introduced in this paper, and b) approximate the probability distribution of the minimum and maximum failure probabilities using extreme value statistics.

Probabilistic studies can be prohibitively expensive because they require repeated finite element analyses of large models. Re-analysis methods have been proposed with the premise to estimate accurately the dynamic response of a structure after a baseline design has been modified, without recalculating the new response. Although these methods increase computational efficiency, they are still not efficient enough for probabilistic analysis of large-scale dynamic systems with low failure probabilities (less or equal to 10-3). This paper presents a methodology that uses deterministic and probabilistic re-analysis methods to generate sample points of the response. Subsequently, tail modeling is used to estimate the right tail of the response PDF and the probability of failure a highly reliable system. The methodology is demonstrated on probabilistic vibration analysis of a realistic vehicle FE model.

In reliability design, often, there is scarce data for constructing probabilistic models. It is particularly challenging to model uncertainty in variables when the type of their probability distribution is unknown. Moreover, it is expensive to estimate the upper and lower bounds of the reliability of a system involving such variables. A method for modeling uncertainty by using Polynomial Chaos Expansion is presented. The method requires specifying bounds for statistical summaries such as the first four moments and credible intervals. A constrained optimization problem, in which decision variables are the coefficients of the Polynomial Chaos Expansion approximation, is formulated and solved in order to estimate the minimum and maximum values of a system’s reliability. This problem is solved efficiently by employing a probabilistic re-analysis approach to approximate the system reliability as a function of the moments of the random variables.

The AUTomotive Open System ARchitecture (AUTOSAR) Development Partnership has published early 2008 the specifications Release 3.0 [1], with a prime focus on the overall architecture, basic software, run time environment, communication stacks and methodology. Heavy developments have taken place in the OEM and supplier community to deliver AUTOSAR loaded cars on the streets starting 2008 [2]. The 2008 achievements have been: Improving the specifications in order to secure the exploitation for body, chassis and powertrain applications Adding major features: safety related functionalities, OBD II and Telematics application interfaces.