Search Results

A method for automatically generating hex-dominant meshes for Computational Fluid Dynamics (CFD) applications is presented in this article. Two important regions of the mesh for any CFD simulation are the interior mesh and the boundary layer mesh. The interior mesh needs to be fine in the critical flow regions to ensure accurate solutions. The proposed method uses Bubble Mesh algorithm which packs bubbles inside the geometry to generate the mesh nodes. Algorithm was tested for sample flow problems and improvements were made to interior and boundary layer mesh generation methods. The interior mesh is generated using directionality and sizing control functions specified on the points of a 3D grid generated over the entire geometry. This offers a flexible control over mesh sizing and local mesh refinement. Boundary layer mesh is important to accurately model the physics of boundary layer near the geometry walls.

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.

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.

Cylinder liners exert a major influence on engine performance, reliability, durability and maintenance. Various combinations of nonmetallic reinforcements and coatings have been used to improve the tribological performance of sleeves or surfaces used in compressors and internal combustion engines in four stroke, two stroke and rotary configurations. In this paper we report the use of a hybrid composite containing silicon carbide and graphite in an aluminum alloy matrix to improve the performance of various small engines and compressors. Material properties of the base material, as well as comparative dynamometer testing, are presented.

A correlation was developed to predict the ignition delay of PRF blends at a wide range of engine-relevant operating conditions. Constant volume simulations were performed using Cantera coupled with a reduced reaction mechanism at a range of initial temperatures from 570-1860K, initial pressures from 10-100atm, oxygen mole percent from 12.6% to 21%, equivalence ratios from 0.30-1.5, and PRF blends from PRF0 to PRF100. In total, 6,480 independent ignition delay simulations were performed. The correlation utilizes the traditional Arrhenius formulation; with equivalence ratio (φ), pressure (p), and oxygen mole percentage (xo2) dependencies. The exponents α, β, and γ were fitted to a third order polynomial with respect to temperature with an exponential roll-off to a constant value at low temperatures to capture the behavior expressed by the reaction mechanism. The location and rate of the roll-off functions were modified by linear functions of PRF.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

The performance of a structural design significantly depends upon the assumptions made on input load. In order to estimate the input load, during the design and development stage of the suspension assembly of a BAJA car, designers and analysts invest immense amount of time and effort to formulate the mathematical model of the design. These theoretical formulations may include idealization errors which can affect the performance of the car as a final product. Due to the errors associated with the assumption of design load, several components might have more weight or may have less strength than needed. This discrepancy between the assumed input load (lab or theoretical studies) and the actual load from the environment can be eliminated by performing a real life testing process using load recovery methodology. Commercial load cells exist in industry to give engineers insight to understanding the complex real world loading of their structures.

Increasingly sophisticated vehicle automation can perform steering and speed control, allowing the driver to disengage from driving. However, vehicle automation may not be capable of handling all roadway situations and driver intervention may be required in such situations. The typical approach is to indicate vehicle capability through displays and warnings, but control algorithms can also signal capability. Psychophysical methods can be used to link perceptual experiences to physical stimuli. In this situation, trust is an important perceptual experience related to automation capability that is revealed by the physical stimuli produced by different control algorithms. For instance, precisely centering the vehicle in the lane may indicate a highly capable system, whereas simply keeping the vehicle within lane boundaries may signal diminished capability.

Distracted driving remains a serious risk to motorists in the US and worldwide. Over 3,000 people were killed in 2013 in the US because of distracted driving; and over 420,000 people were injured. A system that can accurately detect distracted driving would potentially be able to alert drivers, bringing their attention back to the primary driving task and potentially saving lives. This paper documents an effort to develop an algorithm that can detect visual distraction using vehicle-based sensor signals such as steering wheel inputs and lane position. Additionally, the vehicle-based algorithm is compared with a version that includes driving-based signals in the form of head tracking data. The algorithms were developed using machine learning techniques and combine a Random Forest model for instantaneous detection with a Hidden Markov model for time series predictions.

In Equation 10, there should be a plus sign between the first summation term and the quantity in brackets. The acceleration of the participant's vehicle was varied from −0.04 m/s2 to +0.06 m/s2 linearly as a function of the angle of the gas pedal with a 0.1-second delay.

The incorporation of detailed chemistry models in internal combustion engine simulations is becoming mandatory as local, globally lean, low-temperature combustion strategies are setting the path towards a more efficient and environmentally sustainable use of energy resources in transportation. In this paper, we assessed the computational efficiency of a recently developed sparse analytical Jacobian chemistry solver, namely ‘SpeedCHEM’, that features both direct and Krylov-subspace solution methods for maximum efficiency for both small and large mechanism sizes. The code was coupled with a high-dimensional clustering algorithm for grouping homogeneous reactors into clusters with similar states and reactivities, to speed-up the chemical kinetics solution in multi-dimensional combustion simulations.

This paper describes a variable magnetomotive force interior permanent magnet (IPM) machine for use as a traction motor on automobiles in order to reduce total energy consumption during duty cycles and cut costs by using Dy-free magnets. First, the principle of a variable magnetomotive force flux-intensifying IPM (VFI-IPM) machine is explained. A theoretical operating point analysis of the magnets using a simplified model with nonlinear B-H characteristics is presented and the results are confirmed by nonlinear finite element analysis. Four types of magnet layouts were investigated for the magnetic circuit design. It was found that a radial magnetization direction with a single magnet is suitable for the VFI-IPM machine. Magnetization controllability was investigated with respect to the magnet thickness, width and coercive force for the prototype design. The estimated variable motor speed and torque characteristics are presented.

Knock in a Rotax-914 engine was modeled and investigated using an improved version of the KIVA-3V code with a G-equation combustion model, together with a reduced chemical kinetics model. The ERC-PRF mechanism with 47 species and 132 reactions [1] was adopted to model the end gas auto-ignition in front of the flame front. The model was validated by a Caterpillar SI engine and a Rotax-914 engine in different operating conditions. The simulation results agree well with available experimental results. A new engineering quantified knock criterion based on chemical mechanism was then proposed. Hydroperoxyl radical (HO2) shows obvious accumulation before auto-ignition and a sudden decrease after auto-ignition. These properties are considered to be a good capability for HO2 to investigate engine knock problems.

The CFD tools in aerodynamic design process have been commonly used in aerospace industry in last three decades. Although there are many CFD software algorithms developed for aerodynamic applications, the nature of a complex, three-dimensional geometry in incompressible highly separated, viscous flow made computational simulation of ground vehicle aerodynamics more difficult than aerospace applications. However, recent developments in computational hardware and software industry enabled many new engineering applications on computational environment. Traditional production process has largely influenced by computational design, analysis, manufacturing and visualization. Different aspects of linking advanced computational tools and aerodynamic vehicle design challenges are discussed in the present work. Key technologies like parallel computation, turbulence modeling and CFD/wind tunnel compatibility issues are presented.

To model sprays from pressure-swirl atomizers, the connection between the injector and the downstream spray must be considered. A new model for pressure-swirl atomizers is presented which assumes little knowledge of the internal details of the injector, but instead uses available observations of external spray characteristics. First, a correlation for the exit velocity at the injector exit is used to define the liquid film thickness. Next, the film must be modeled as it becomes a thin, liquid sheet and breaks up, forming ligaments and droplets. A linearized instability analysis of the breakup of a viscous, liquid sheet is used as part of the spray boundary condition. The spray angle is estimated from spray photographs and patternator data. A mass averaged spray angle is calculated from the patternator data and used in some of the calculations.

The three-dimensional CFD codes KIVA-II and KIVA-3 have been used together to study the effects of intake generated in-cylinder flow structure on fuel-air mixing and combustion in a direct injected (DI) Diesel engine. In order to more accurately account for the effect of intake flow on in-cylinder processes, the KIVA-II code has been modified to allow for the use of data from other CFD codes as initial conditions. Simulation of the intake and compression strokes in a heavy-duty four-stroke DI Diesel engine has been carried out using KIVA-3. Flow quantities and thermodynamic field information were then mapped into a computational grid in KIVA-II for use in the study of mixing and combustion. A laminar and turbulent timescale combustion model, as well as advanced spray models, including wave breakup atomization, dynamic drop drag, and spray-wall interaction has been used in KIVA-II.

An overview of a design methodology based on finite element and fatigue analyses is described for the design of welded structures. Fatigue life is a primary design consideration for agricultural equipment structural members. The re-design of an agricultural implement frame and hitch assembly is used to illustrate the methodology for optimal design and improved fatigue life. A description of the finite element mode of the agricultural implement and the simulated loading history are described. Results from finite element analysis are used to optimize the frame member sizes and calculate the fatigue strength characteristics. This design analysis strategy provides improved fatigue life characteristics for the implement frame or welded structure.

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.

A detailed mathematical model of the thermodynamic events of a crankcase scavenged, two-stroke, SI engine is described. The engine is divided into three thermodynamic systems: the cylinder gases, the crankcase gases, and the inlet system gases. Energy balances, mass continuity equations, the ideal gas law, and thermodynamic property relationships are combined to give a set of coupled ordinary differential equations which describe the thermodynamic states encountered by the systems of the engine during one cycle of operation. A computer program is used to integrate the equations, starting with estimated initial thermodynamic conditions and estimated metal surface temperatures. The program iterates the cycle, adjusting the initial estimates, until the final conditions agree with the beginning conditions, that is, until a cycle results.

The infrared radiation method of compression and end-gas temperature measurement was applied to the problem of measuring gas temperatures up to the time of knock. Pressure data were taken for each run on a CFR engine with mixtures of isooctane and n-heptane under both knocking and nonknocking conditions. Main engine parameters studied were the intake pressure, intake temperature, and engine speed. The rate and extent of chemical energy release were calculated from the temperature and pressure histories using an energy balance. The computed rates of chemical energy release were correlated to a chain-type kinetic model