Search Results

Accurately predicting the evolution of soot mass and soot particle numbers under engine conditions is critical to advanced engine design. A detailed soot-chemistry model that can capture soot under gasoline and diesel conditions without tuning is necessary for such predictions. Building confidence in the predictive usage of the chemistry in engine simulations requires validating the soot kinetics over a wide range of operating conditions and fuels, using data from different experimental techniques, and using sources from laboratory flames to engines. This validation study focuses on a soot-chemistry model that considers multiple nucleation, growth, and oxidation reaction pathways. It involves 14 gas-phase precursors and considers the effect of different soot-particle surface sites.

The safety, performance, and operational life of power dense Lithium-ion batteries used in Hybrid and Electric Vehicles are dependent on the operating temperature. Modeling and simulation are essential tools used to accelerate the design process of optimal thermal management systems. However, high-fidelity 3D computational fluid dynamics (CFD) simulation of such systems is often difficult and computationally expensive. In this paper, we demonstrate a multi-part coupled system model for simulating the heating/cooling system of the traction battery at a module level. We have reduced computational time by employing reduced-order modeling (ROM) framework on separate 3D CFD models of the battery module and the cooling plate. The order of the thermal ROM has also been varied to study the trade-off between accuracy, fidelity, and complexity. The ROMs are bidirectionally coupled to an empirical battery model built from in-house test data.

This study focuses on Lagrangian spray models that are commonly used in engine CFD simulations. In modeling sprays, droplet collision is one of the physical phenomena that must be accounted for. There are two main parts of droplet collision models for sprays - detecting colliding pairs of droplets and predicting the outcomes of these collisions. For the first part, we focus on the efficiency of the algorithm. We present an implementation of the arbitrary adaptive collision mesh model of Hou and Schmidt [1], and examine its efficiency in dealing with large simulations. Through theoretical analysis and numerical tests, we show that the computational cost of this model scales pseudo-linearly with respect to the number of parcels in the sprays. Regarding the second part, we examine the variations in existing phenomenological models used for predicting binary droplet collision outcomes. A quantitative accuracy metric is used to evaluate the models with respect to the experimental data set.

Cylinder head design is one of the most involved disciplines in engine design. Whether designing an altogether new head or revamping an old one, several different coupled and inter-dependent technologies ranging from heat transfer, fluid flow, combustion, material non-linearity, structural and fatigue have to be accounted. Simultaneous designing of ports, jacket and combustion chamber leads to cylinder head design, which is then tested for its strength and durability. Traditionally a series of analytical, empirical, test-based and simulation based exercises are conducted to design cylinder heads. With increasing pressure on reducing cost and turnaround time, focus on moving towards a quasi-simulation based design and development of cylinder heads is gaining strength. This paper talks about how a simulation driven process for cylinder head design can be developed.

A commonly used physics based electrochemisty model for a lithium-ion battery cell was first proposed by professor Newman in 1993. The model consists of a tightly coupled set of partial differential equations. Due to the tight coupling between the equations and the 2d implementation due to the particle modeling, and thus called pseudo-2d in literature, numerically obtaining a solution turns out to be challenging even for a lot of commercial softwares. In this paper, the VHDL-AMS language is used to solve the set of equations. VHDL-AMS allows the user to focus on the physical modeling rather than numerically solving the governing equations. In using VHDL-AMS, the user only needs to specify the governing equations after spatial discretization. A simulation environment, which supports VHDL-AMS, can then be used to solve the governing equations and also provides both pre- and post- processing tools.

Accurate modeling of the characteristics of diesel-engine combustion leads to more efficient design. Accurate modeling in turn depends on correctly capturing spray dynamics, turbulence, and fuel chemistry. This work presents a computational fluid dynamics (CFD) investigation of a well characterized small-bore direct injection diesel engine at Sandia National Laboratories’ Combustion Research Facility. The engine has been studied for two piston-bowls geometries and various injection timings. Simulation of these conditions test the predictive capabilities of our approach to diesel engine modeling using Ansys Forte. An experimental database covering a wide range of operating conditions is provided by the Engine Combustion Network for this engine, which is used to validate our modeling approach. Automatic and solution-adaptive meshing is used, and the recommended settings are discussed.

This paper studies the level of confidence and applicability of CFD simulations using steady-state Reynolds-Averaged Navier-Stokes (RANS) in predicting aerodynamic performance losses on swept-wings due to contamination with ice accreted in-flight. The wing geometry selected for the study is the 65%-scale Common Research Model (CRM65) main wing, for which NASA Glenn Research Center’s Icing Research Tunnel has generated experimental ice shapes for the inboard, mid-span, and outboard sections. The reproductions at various levels of fidelity from detailed 3D scans of these ice shapes have been used in recent aerodynamic testing at the Office National d’Etudes et Recherches Aérospatiales (ONERA) and Wichita State University (WSU) wind tunnels. The ONERA tests were at higher Reynolds number range in the order of 10 million, while the WSU tests were in the order of 1 million.

A CFD simulation methodology is presented to calculate blockage due to ice crystal icing of the IGV passages of a gas turbine engine. The computational domain consists of six components and includes the nacelle, the full bypass and the air induction section up to the second stage of the low-pressure compressor. The model is of a geared turbofan with a fan that rotates at 4,100 rpm and a low-pressure stage that rotates at 8,000 rpm. The flight conditions are based on a cruising speed of Mach 0.67 in Appendix-D icing conditions with an ice crystal content is 4.24 g/m3. Crystal bouncing, and re-entrainment is considered in the calculations, along with variable relative humidity and crystal melting due to warmer temperatures within the engine core. Total time of icing is set to 20 seconds. The CFD airflow and ice crystal simulations are performed on the full 6-stage domain. The initial icing calculation determines which stage will be chosen for a more comprehensive analysis.

A 3D CFD methodology is presented to simulate ice build-up on propeller blades exposed to known icing conditions in flight, with automatic blade pitch variation at constant RPM to maintain the desired thrust. One blade of a six-blade propeller and a 70-passenger twin-engine turboprop are analyzed as stand-alone components in a multi-shot quasi-steady icing simulation. The thrust that must be generated by the propellers is obtained from the drag computed on the aircraft. The flight conditions are typical for a 70-passenger twin-engine turboprop in a holding pattern in Appendix C icing conditions: 190 kts at an altitude of 6,000 ft. The rotation rate remains constant at 850 rpm, a typical operating condition for this flight envelope.

A CFD simulation methodology for the inclusion of the post-impact trajectories of splashing/bouncing Supercooled Large Droplets (SLDs) and film detachment is introduced and validated. Several scenarios are tested to demonstrate how different parameters affect the simulations. Including re-injecting droplet flows due to splashing/bouncing and film detachment has a significant effect on the accuracy of the validations shown in the article. Validation results demonstrate very good agreement with the experimental data. This approach is then applied to a full-scale twin-engine turboprop to compute water impingement on the wings and the empennage.

The importance of the variation of relative humidity across turbomachinery engine components for in-flight icing is shown by numerical analysis. A species transport equation for vapor has been added to the existing CFD methodology for the simulation of ice growth and water flow on engine components that are subject to ice crystal icing. This entire system couples several partial differential equations that consider heat and mass transfer between droplets, crystals and air, adding the cooling of the air due to particle evaporation to the icing simulation, increasing the accuracy of the evaporative heat fluxes on wetted walls. Three validation cases are presented for the new methodology: the first one compares with the numerical results of droplets traveling inside an icing tunnel with an existing evaporation model proposed by the National Research Council of Canada (NRC).

With the emphasis of electrification in automotive industry, tremendous efforts are made to develop electric motors with high efficiency and power density, and reduce noise, vibration and harshness (NVH). A multiphysics simulation workflow is used to predict the eccentricity-induced noise for GM’s Bolt EV motor. Both static and dynamic eccentricities are investigated along with axial tilt. Analysis results show that these eccentricities play a critical role in the NVH behavior of the motor assembly. Transient electromagnetic (EM) analysis is performed first by extruding 2D stator and rotor sections to form 3D EM models. Sector model is duplicated to form full 360-degree model. Stator is split into three rotated sections to characterize stator skew, and the skew between two sections of rotor and magnets are also modelled. Sinusoidal current is applied and lumped-sum forces on each stator tooth are computed.

A full-automated CFD mesh generation technique has been developed and implemented for 3-D aircraft icing simulations to permit robust 45-minute ice accretion simulations in support of icing certification campaigns. The changes in the shape of the aircraft surfaces due to accreting ice and their effects on the air and droplet flow are accounted for in a quasi-steady manner by subdividing the total icing time into sequential steps of shorter duration, updating the computational grid at each step. This “multi-shot” ice accretion approach requires robust and accurate grid re-meshing for it to be embedded in engineering design and analysis workflows. ANSYS FENSAP-ICE has been coupled to Fluent Meshing to take advantage of generic and highly automated surface displacement and mesh wrapping tools. A wide spectrum of geometries is supported, ranging from full-size aircraft to air data probes, turbomachinery components, rotors and propellers.

In order to achieve the purpose of saving energy and reducing emission, the improvement of aerodynamic performance plays an increasingly crucial role for car manufacturers. Previous studies have confirmed the validity of gradient-based adjoint algorithm for its high efficiency in shape optimization. In this paper, two important aspects of adjoint approach were explored. One is vehicle aerodynamic optimization with multiple objectives, and the other is using time-averaged flow results as the primal solution, both are issues of high interest in recent applications. First, adjoint shape optimization with steady-state and time-averaged flow simulations were respectively calculated and comparatively discussed based on a production SUV. The shape modifications of the two cases indicated that the impact of primal solution on design change could not be neglected, due to the different intrinsic codes of steady and transient turbulence models.

Brake squeal is an instability issue with many parameters. This study attempts to assess the effect of thermal load on brake squeal behavior through finite element computation. The research can be divided into two parts. The first step is to analyze the thermal conditions of a brake assembly based on ANSYS Fluent. Modeling of transient temperature and thermal-structural analysis are then used in coupled thermal-mechanical analysis using complex eigenvalue methods in ANSYS Mechanical to determine the deformation and the stress established in both the disk and the pad. Thus, the influence of thermal load may be observed when using finite element methods for prediction of brake squeal propensity. A detailed finite element model of a commercial brake disc was developed and verified by experimental modal analysis and structure free-free modal analysis.

In this case study, the thermoforming of an automotive instrument panel is considered. The effect of different oven settings on the final material distribution is studied using structural FEA simulation. The variable thickness distribution of the thermoformed part is mapped onto a structural model using a new simple mapping algorithm, and a structural FEA simulation is carried out to examine the final warpage of the instrument panel. The simulation predicts that the minimum thickness of the formed part can be increased by 10% by optimizing the oven settings. Although the optimized process uses oven settings that are less uniform than the baseline settings, the model indicates that warpage experienced by the optimized part will be less than that of the baseline case.

Aromatics constitute a significant portion of refinery fuels. Characterizing the impact of various aromatic components on combustion and emissions facilitates formulation of surrogate fuels for engine simulations. The impact of blending aromatics in fuel surrogates is usually nonlinear for ignition characteristics responsible for knocking in spark engines and for combustion phasing in diesel engines. In this work, we have characterized the behavior of nine aromatics components under engine-relevant conditions. A self-consistent and validated detailed kinetics mechanism has been developed for gasoline and diesel surrogates that contains toluene, ethylbenzene, n-propylbenzene, n-butylbenzene, isomers of xylene, 1,2,4-trimethylbenzene, and 1-methylnaphthalene. Numerical experiments using 0-D and 1-D models have been performed to study the relative behavior of these aromatics for different reacting conditions.

This paper introduces the Lagrangian particle tracking technology readily available in Ansys Fluent in the in-flight icing simulation workflow, which normally uses the Eulerian approach for droplet flows. The Lagrangian solver is incorporated in the Fluent Icing workspace which is to become the next-gen in-flight icing simulation tool provided by Ansys. Lagrangian tracking will eventually be used for SLD and ice crystal rebound and re-impingement calculations in the Ansys workflow. Here we introduce some preliminary results with the current state of its implementation as of Fluent Icing release 2023R2. Example cases include several selections from the 1st Ice prediction workshop with experimental comparisons as well as results obtained earlier with the Eulerian droplet solution strategy. Collection efficiency comparisons on clean geometries show good agreement between Eulerian and Lagrangian methods when the particle seeds are in the millions range.

In-flight icing significantly influences the design of large passenger aircraft. Relevant aspects include sizing of the main aerodynamic surfaces, provision of anti-icing systems, and setting of operational restrictions. Empennages of large passenger aircraft are particularly affected due to the small leading edge radius, and the requirement to generate considerable lift for round out and flare, following an extended period of descent often in icing conditions. This paper describes a CFD-based investigation of the effects of sweep on the aerodynamic performance of a novel forward-swept horizontal stabilizer concept in icing conditions. The concept features an unconventional forward sweep, combined with a high lift leading edge extension (LEX) located within a fuselage induced droplet shadow zone, providing passive protection from icing.

Airflow over the Ahmed body is simulated by means of the Embedded LES technique. This is a zonal approach which allows the LES turbulence model to be used in a sub-domain while the rest of the domain is solved using a RANS turbulence model, i.e., the kω-SST model. This allows the LES zone to be restricted to the rear end of the vehicle where the unsteadiness of the flow must be accurately predicted. A comparison with full RANS modeling and experiments is reported.