Search NASA⌕ Search

SEARCH · Search NASA

Results for “RANS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ECFA Higgs, electroweak, and top Factory Study

The ECFA Higgs, electroweak, and top Factory Study ran between 2021 and 2025 as a broad effort across the experimental and theoretical particle physics communities, bringing together participants from many different proposed future collider projects. Activities across three main working groups advanced the joint development of tools and analysis techniques, fostered new considerations of detector design and optimisation, and led to a new set of studies resulting in improved projected sensitivities across a wide physics programme. This report demonstrates the significant expansion in the state-of-the-art understanding of the physics potential of future e+e- Higgs, electroweak, and top factories, and has been submitted as input to the 2025 European Strategy for Particle Physics Update.

Altmann, J. [Monash U.]↗
DOE OSTI2023

High Efficiency 2-Bladed Utility Wind Turbine Enhanced by CoFlow Jet Active Flow Control

The aerodynamic performance and power output of a high efficiency Co-Flow Jet (CFJ) 2-bladed Wind Turbine is studied in this paper. CFJ is a Zero-Net-Mass-Flux (ZNMF) active flow control method that dramatically increases airfoil lift coefficient and suppresses flow separation at a low energy expenditure. The 3D Reynolds Averaged Navier-Stokes (RANS) equations with one-equation Spalart-Allmaras (SA) turbulence model are solved to simulate the 3D flows of the wind turbines. The CFJ-Wind Turbine in this paper is modified from the NREL 3-bladed 3.4 MW reference wind turbine with the same RPM, blade length, distribution of chord, twist, and bending. The CFJ 2-bladed wind turbine utilizes the CFJ-NACA 64 series airfoils at most of the blade span from 31.88\% to the tip. The CFJ injection and suction slots are implemented along the full blade span to enhance the effectiveness and power output. The study shows that the 2-bladed CFJ turbine significantly increases power output compared with the 3-bladed baseline turbine at all wind speeds by 15.8\% at the wind speed of 9.8 m/s and by 25.9\% at 4 m/s respectively. The efficiency enhancement is substantially more at the low wind speed than at high wind speed. This would be very beneficial to increase the annual power output or capacity factor as wind turbines work at the speeds lower than the rated speed for most of the time. The results show that the flow field around the blade surface presents a significant enhanced suction peak near the blade leading edge, which results in a very high tangential force and thus increased wind turbine net power coefficient. Consistent with the previous study for the small NREL Phase VI research wind turbine, this paper demonstrates that CFJ 2-bladed wind turbine is an effective technology to extract more power from the wind at all speeds. This study is just an initial effort to demonstrate the advantage of the CFJ blade with no design optimization. It is believed that there is more potential to further enhance the efficiency improvement.

active flow control↗
DOE OSTI2024

Enriched Background Isotope Study (EBIS): Analysis of 14C-Enriched Carbon Cycle in Soils and Litter at Forested Oak Ridge and AmeriFlux Sites, 2001-2011

These data provide a record of the multi-year, multi-institutional Enriched Background Isotope Studies (EBIS) projects that ran from 2000 through 2011. Elevated levels of 14C enriched CO2 in the air and soil atmosphere as well as leaf, stem, and root tissues were observed on the Oak Ridge Reservations (ORR) during the summer of 1999, and were attributed to local incinerator activities on and/or near the ORR (Trumbore et al. 2002). The isolated enrichment of the background levels of 14C in local forest ecosystem represented a unique opportunity to study unresolved carbon cycling processes such as the contribution of leaf versus root litter contributions to soil carbon accumulation, the rate of vertical transport of carbon into deep soil storage pools, and the differential contribution of physicochemical versus faunal driven processes to soil carbon cycling and sequestration. Leaf litter from the local enriched forest was transplanted to selected sites on the ORR and to selected AmeriFlux study sites to study soil C cycling across a range of soils and climates. The EBIS research projects provide data on C flux from litter sources to mineral soil sinks for United States eastern hardwood forests necessary for testing process hypotheses and judging efficacy of soil C cycling models. Experimental results from this study are being used to parameterize and refine existing carbon dynamics models, the quantification of the long-term fate of ecosystem carbon inputs and as a means to judge the potential for ecosystem carbon sequestration via enhance litter inputs to soil. EBIS observations support conclusions that intra- and inter-annual soil carbon cycling in hardwood forest soils should be characterized as a least a two-compartment system where surface leaf-litter and belowground root turnover represent primary carbon sources for organic-layer and mineral-soil carbon cycles, respectively. EBIS experiments were conducted to complete enriched litterfall maniplations in upland forests on Ultisol and Inceptisol soils of the Oak Ridge Reservation, Oak Ridge, Tennessee. We also collected additional14C-enriched materials for new experimental applications, and applied those materials to multiple AmeriFlux sites over a range of climatic, edaphic and biological conditions. The research provided data for addressing DOE's goal of understanding mechanisms controlling C flux, and for the improvement of models to be applied to policy discussions regarding the safe levels of greenhouse gases for the earth's system. There are 5 data files provided in comma separated (*.csv) format for vegetation, field litter, soil and air [C] and C isotope data from the EBIS studies and associated environmental data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
DOE OSTI2021

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗
DOE OSTI2023

Numerical and boundary condition effects on the prediction of detonation engine behavior using detailed numerical simulations

High-fidelity numerical simulations of an experimental rotating detonation engine with discrete fuel/air injection were conducted. A series of configurations with different feed-plenum pressures but with constant equivalence ratio were studied. Detailed chemical kinetics for the hydrogen/air system is used. A resolution study for the full rotating detonation engine (RDE) system simulation is also conducted. Two kinds of boundary conditions, a total pressure boundary and a constant mass flow rate boundary, are used to assess the effects of the inlet boundary. As mass flow rate is increased, the total pressure boundary causes more error in the axial pressure distribution while the constant mass flow rate gives a better solution for all cases ran. The simulations confirm experimental findings, and reproduce qualitative as well as some of the quantitative trends. These results demonstrate that a) fuel-air mixing is highly non-uniform within the detonation chamber, leading to variations in local equivalence ratio, b) the fuel and oxidizer injectors experience significant backflow as the detonation wave passes over, but recover at different rates which further augments the inefficiencies in mixing, and c) parasitic combustion in the mixing region makes the detonation wave weak by extending the reaction zone across the wave.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗
DOE OSTI2023

Simulations for Planning of Liquid Hydrogen Spill Test

In order to better understand the complex pooling and vaporization of a liquid hydrogen spill, Sandia National Laboratories is conducting a highly instrumented, controlled experiment inside their Shock Tube Facility. Simulations were run before the experiment to help with the planning of experimental conditions, including sensor placement and cross wind velocity. This paper describes the modeling used in this planning process and its main conclusions. Sierra Suite’s Fuego, an in-house computational fluid dynamics code, was used to simulate a RANS model of a liquid hydrogen spill with five crosswind velocities: 0.45, 0.89, 1.34, 1.79, and 2.24 m/s. Two pool sizes were considered: a diameter of 0.85 m and a diameter of 1.7. A grid resolution study was completed on the smaller pool size with a 1.34 m/s crosswind. A comparison of the length and height of the plume of flammable hydrogen vaporizing from the pool shows that the plume becomes longer and remains closer to the ground with increasing wind speed. The plume reaches the top of the facility only in the 0.45 m/s case. From these results, we concluded that it will be best for the spacing and location of the concentration sensors to be reconfigured for each wind speed during the experiment.

08 HYDROGEN↗
DOE OSTI2025

Assessing the Heat Transfer Modeling Capabilities of CFD Software for Involute-Shaped Plate Research Reactors

The ongoing efforts to convert High-Performance Research Reactors (HPRRs) using Highly Enriched Uranium (HEU) to Low-Enriched Uranium (LEU) fuel require reliable thermal–hydraulic assessments of modified core designs. The involute-shaped fuel plates used in several major HPRRs present unique modeling challenges due to their compact core geometries and high heat flux conditions. This study evaluates the capability of three commercial CFD tools, STAR-CCM+, COMSOL, and ANSYS CFX, to predict cladding-to-coolant heat transfer using Reynolds-Averaged Navier–Stokes (RANS) methods within the thermal–hydraulic regimes of involute-shaped plate reactors. Broad sensitivity analysis was conducted across a range of reactor-relevant parameters using two turbulence models (k−ϵ and k−ω SST) and different near-wall treatment strategies. The results were benchmarked against the Sieder–Tate correlation and experimental data from historic studies. The codes produced consistent results, showing good agreement with the empirical correlation of Sieder–Tate and the experimental measurements. The findings support the use of these commercial CFD codes as effective tools for assessing the thermal–hydraulic performance of involute-shaped plate HPRRs and guide future LEU core development.

CFD↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2024

Design and Characterization of Highly Diffusive Turbine Vanes Suitable for Transonic Rotating Detonation Combustors

In rotating detonation engines the turbine inlet conditions may be transonic with unprecedented unsteady fluctuations. To ensure an acceptable engine performance, the turbine passages must be suited to these conditions. This article focuses on designing and characterizing highly diffusive turbine vanes to operate at any inlet Mach number up to Mach 1. First, the effect of pressure loss on the starting limit is presented. Afterward, a multi-objective optimization with steady RANS simulations, including the endwall and 3D vane design is performed. Compared to previous research, significant reductions in pressure loss and stator-induced rotor forcing are obtained, with an extended operating range and preserving high flow turning. Finally, the influence of the inlet boundary layer thickness on the vane performance is evaluated, inducing remarkable increases in pressure loss and downstream pressure distortion. Employing an optimization with a thicker inlet boundary layer, specific endwall design recommendations are found, providing a notable improvement in both objective functions.

Grasa, Sergio↗
DOE OSTI2021

Utilizing Amino Acid Composition and Entropy of Potential Open Reading Frames to Identify Protein-Coding Genes

One of the main steps in gene-finding in prokaryotes is determining which open reading frames encode for a protein, and which occur by chance alone. There are many different methods to differentiate the two; the most prevalent approach is using shared homology with a database of known genes. This method presents many pitfalls, most notably the catch that you only find genes that you have seen before. The four most popular prokaryotic gene-prediction programs (GeneMark, Glimmer, Prodigal, Phanotate) all use a protein-coding training model to predict protein-coding genes, with the latter three allowing for the training model to be created ab initio from the input genome. Different methods are available for creating the training model, and to increase the accuracy of such tools, we present here GOODORFS, a method for identifying protein-coding genes within a set of all possible open reading frames (ORFS). Our workflow begins with taking the amino acid frequencies of each ORF, calculating an entropy density profile (EDP), using KMeans to cluster the EDPs, and then selecting the cluster with the lowest variation as the coding ORFs. To test the efficacy of our method, we ran GOODORFS on 14,179 annotated phage genomes, and compared our results to the initial training-set creation step of four other similar methods (Glimmer, MED2, PHANOTATE, Prodigal). We found that GOODORFS was the most accurate (0.94) and had the best F1-score (0.85), while Glimmer had the highest precision (0.92) and PHANOTATE had the highest recall (0.96).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Clarissa Family Age from the Yarkovsky Effect Chronology

The Clarissa family is a small collisional family composed of primitive C-type asteroids. It is located in a dynamically stable zone of the inner asteroid belt. In this work we determine the formation age of the Clarissa family by modeling planetary perturbations as well as thermal drift of family members due to the Yarkovsky effect. Simulations were carried out using the SWIFT-RMVS4 integrator modified to account for the Yarkovsky and Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effects. We ran multiple simulations starting with different ejection velocity fields of fragments, varying proportion of initially retrograde spins, and also tested different Yarkovsky/YORP models. Our goal was to match the observed orbital structure of the Clarissa family which is notably asymmetrical in the proper semimajor axis, a {sub p}. The best fits were obtained with the initial ejection velocities ≲20 m s{sup −1} of diameter D ≃ 2 km fragments, ∼4:1 preference for spin-up by YORP, and assuming that ≃80% of small family members initially had retrograde rotation. The age of the Clarissa family was found to be t {sub age} = 56 ± 6 Myr for the assumed asteroid density ρ = 1.5 g cm{sup −3}. Small variation of density to smaller or larger value would lead to slightly younger or older age estimates. This is the first case where the Yarkovsky effect chronology has been successfully applied to an asteroid family younger than 100 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stellar Activity Manifesting at a One-year Alias Explains Barnard b as a False Positive

Barnard’s star is among the most studied stars given its proximity to the Sun. It is often considered the radial velocity (RV) standard for fully convective stars due to its RV stability and equatorial decl. Recently, an Msin i=3.3M{sub ⊕} super-Earth planet candidate with a 233 day orbital period was announced by Ribas et al. New observations from the near-infrared Habitable-zone Planet Finder (HPF) Doppler spectrometer do not show this planetary signal. We ran a suite of experiments on both the original data and a combined original + HPF data set. These experiments include model comparisons, periodogram analyses, and sampling sensitivity, all of which show the signal at the proposed period of 233 days is transitory in nature. The power in the signal is largely contained within 211 RVs that were taken within a 1000 day span of observing. Our preferred model of the system is one that features stellar activity without a planet. We propose that the candidate planetary signal is an alias of the 145 day rotation period. This result highlights the challenge of analyzing long-term, quasi-periodic activity signals over multiyear and multi-instrument observing campaigns.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Probability of Forming Gaps in the GD-1 Stream by Close Encounters of Globular Clusters

One of the most intriguing properties of the GD-1 stellar stream is the existence of three gaps. If these gaps were formed by close encounters with dark matter subhalos, the GD-1 stream opens an exciting window through which we can see the size, mass, and velocity distributions of the dark matter subhalos in the Milky Way. However, in order to use the GD-1 stream as a probe of the dark matter substructure, we need to disprove that these gaps are not due to the perturbations from baryonic components of the Milky Way. Here we ran a large number of test-particle simulations to investigate the probability that each of the known globular clusters (GCs) can form a GD-1-like gap, by using the kinematical data of the GD-1 stream and GCs from Gaia early data release 3 and by fully taking the observational uncertainty into account. We found that the probability that all of the three gaps were formed by GCs is as low as 1.7 × 10 -5 , and the expected number of gaps formed by GCs is only 0.057 in our fiducial model. Our result highly disfavors a scenario in which GCs form the gaps. Given that other baryonic perturbers (e.g., giant molecular clouds) are even less likely to form a gap in the retrograde-moving GD-1 stream, we conclude that at least one of the gaps in the GD-1 stream was formed by dark matter subhalos if the gaps were formed by flyby perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Detection of the Orbital Modulation of Fe Kα Fluorescence Emission in Centaurus X-3 Using the High-resolution Spectrometer Resolve on board XRISM

Abstract The Fe K α fluorescence line emission in X-ray spectra is a powerful diagnostic tool for various astrophysical objects to reveal the distribution of cold matter around photoionizing sources. The advent of the X-ray microcalorimeter on board the XRISM satellite will bring new constraints on the emission line. We present one of the first such results for the high-mass X-ray binary Centaurus X-3, which is composed of an O-type star and a neutron star (NS). We conducted a 155 ks observation covering an entire binary orbit. A weak Fe K α line was detected in all orbital phases at an equivalent width (EW) of 10–20 eV. We found for the first time that its radial velocity (RV) is sinusoidally modulated by the orbital phase. The RV amplitude is 248 ± 13 km s −1 , which is significantly smaller than the value (391 km s −1 ) expected if the emission is from the NS surface, but is consistent if the emission takes place at the O star surface. We discuss several possibilities of the line production site, including the NS surface, O star surface, O star wind, and accretion stream from the O star to the NS. We ran radiative transfer calculation for some of them assuming spherically symmetric density and velocity profiles and an isotropic distribution of X-ray emission from the NS. None of them explains the observed EW and velocity dispersion dependence on the orbital phase, suggesting that more elaborated modeling is needed. In other words, the present observational results have the capability to constrain deviations from these assumptions.

High energy astrophysics↗
DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗