Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance modeling”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗
DOE OSTI2020

Risk Assessment and Ventilation Modeling for Hydrogen Release in Vehicle Repair Garages

The availability of repair garage infrastructure for hydrogen fuel cell vehicles is becoming increasingly important for future industry growth. Ventilation requirements for hydrogen fuel cell vehicles can affect both retrofitted and purpose-built repair garages and the costs associated with these requirements can be significant. A hazard and operability study (HAZOP) was performed to identify key risk-significant scenarios related to hydrogen vehicles in a repair garage. Detailed simulations and modeling were performed using appropriate computational tools to estimate the location, behavior, and severity of hydrogen release based on key HAZOP scenarios. This work compares current fire code requirements to an alternate ventilation strategy to further reduce potentially hazardous conditions. Overall, the amount of flammable mass of hydrogen at any one time in the simulation is low compared to the total mass of hydrogen released, due to the low flow rate of a low pressure release. It is shown that position, direction, and velocity of ventilation have a significant impact on the amount of instantaneous flammable mass in the domain.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Scoping Study for Fast Flux Testing in the Advanced Test Reactor

The value of fast spectrum reactors remains prominent in the nuclear technology portfolio. The performance of these reactors can be maximized with advancements in nuclear fuel technologies, but development of these technologies is currently held back by lack of fast spectrum test reactors available to the United States. Spectral modification of experiment positions in the thermal spectrum Advanced Test Reactor (ATR) has long been used to support fast reactor fuel development, but these methods have not been progressed to their full potential. This study investigated the use of concentric rings of aluminum-clad fuel plates in ATR flux traps and thermal neutron absorbing filters to increase fast neutron flux on test specimens. This concept was termed the Boosted Energy Advanced Spectrum Test (BEAST). This approach will enable irradiation of advanced fuel designs in prototypic-length fuel pins and representative flux environment to support post irradiation exams, enable transient testing, and produce the type of data that will permit lead test assembly irradiations in true Sodium Fast Reactors (SFRs) when they become available. Neutronic predictions were performed to investigate BEAST design options and thermal hydraulic models were produced to ensure feasibility of BEAST. Two versions were considered based on the geometric limitations of ATR’s small and large flux traps. The small version was found to be preferable due to slightly higher fast flux and fast-to-thermal neutron ratio. Perhaps more influentially, the small flux trap option was also preferred to avoid conflict with ongoing very high temperature reactor fuel irradiation programs in ATR’s large northeast flux trap. The small flux trap option provided less than half the test volume of the large version, but still had adequate volume for seven SFR pins in cross section which could be stacked two-high in ATR’s 1.2m long core to accommodate up to 14 EBR-II size pins. The preference for the small flux trap configuration should be revisited if additional collaborative test programs emerge with the need to irradiate a significant volume of additional specimens. Calculations were performed regarding a lithium deuteride ring to convert thermal neutrons into 14 MeV fusion neutrons. At the time this report was written these calculations were partially complete and it remains to be seen whether the concept would be worth including in BEAST. Given the preference for the small flux trap option, which does not afford enough room for the 14 MeV ring, it was concluded to defer future work on the lithium deuteride ring. This decision could be revisited if fusion material research programs emerge for collaborative testing in BEAST. A cadmium-lined specimen holder design was found to be adequate in filtering thermal neutrons and preferred over other neutron absorbers based on past experience with cadmium baskets. It was acknowledged that cadmium-bearing hardware would become depleted and need to be replaced occasionally, which appeared feasible from a mechanical design perspective. Neutronic studies investigated different enrichment levels in the booster fuel using uranium-molybdenum alloy dispersion fuel which has performed well in past ATR irradiations. Both options were able to drive fuel pins to SFR-like fission heating rates. The high enriched booster fuel option outperformed the low enriched option by ~20% on key metrics including fast flux and fast-to-thermal ratio, but the low enriched option was favored in order to broaden options for potential fuel suppliers. The preferred BEAST design options including cadmium filter with low enriched booster fuel in the small flux trap configuration was predicted to achieve 6.2E14 n/cm2sec fast flux (>0.1 MeV) with a fast-to-thermal ratio of 44.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Heavy-Duty Low-Temperature and Diesel Combustion & Heavy-Duty Combustion Modeling (FY 2018 Annual Progress Report)

Regulatory drivers and market demands for lower pollutant emissions, lower carbon dioxide emissions, and lower fuel consumption motivate the development of clean and fuel-efficient engine operating strategies. Most current production engines use a combination of both in-cylinder and exhaust emission control strategies to achieve these goals. The emissions and efficiency performance of in-cylinder strategies depend strongly on flow and mixing processes associated with fuel injection. Both heavy- and light-duty engine/vehicle manufacturers use multiple-injection strategies to reduce noise, emissions, and fuel consumption. For both conventional and low-temperature diesel combustion, the state of knowledge and modeling tools for multiple injections are far less advanced than for single-injection strategies. Engine efficiency is limited to some degree by tradeoffs that must be accepted to meet particulate matter (including soot) emissions limits. Recent work on this project has filled some knowledge gaps on soot oxidation with multiple injections, and the current work for Fiscal Year (FY) 2018 addresses knowledge gaps on soot formation for multiple injections. While total in-cylinder soot is readily measured, discerning formation from oxidation is difficult. The FY 2018 experiments are designed to create in-cylinder conditions at the threshold of soot formation, where processes that affect soot formation can be more readily discerned. Soot formation pathways under such conditions are fraught with uncertainties, and soot models significantly overpredict polyaromatic hydrocarbon (PAH) and soot, so experimental data at these conditions will provide much needed data for improvements to PAH and soot models.

02 PETROLEUM↗
DOE OSTI2024

In-Situ Species Concentration Measurements In Ammonia-Mix Flames Using Ftir Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH 3 /H 2 ) and ammonia-natural gas (NH 3 /CH 4 ) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H 2 O absorption.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In-Situ Species Concentration Measurements in Ammonia-Mix Flames Using FTIR Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH3/H2) and ammonia-natural gas (NH3/CH4) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H2O absorption.

Bedick, Clinton↗
DOE OSTI2020

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗
DOE OSTI2021

Ex Situ Catalytic Fast Pyrolysis of Lignocellulosic Biomass to Hydrocarbon Fuels: 2020 State of Technology

Biomass catalytic fast pyrolysis (CFP) research was funded by the Bioenergy Technologies Office (BETO) of the U.S. Department of Energy (DOE) to enable technical breakthroughs to help reduce the cost of conversion of biomass to liquid hydrocarbon fuels. The key goals were to achieve a mature plant modeled Minimum Fuel Selling Price of $3/GGE, with greater than 60% greenhouse gas (GHG) reduction over petroleum derived gasoline by 2022. Experimental results at NREL and PNNL, along with lifecycle analysis by ANL show that we achieved these goals based on performance at the bench scale. A Pt/TiO2 catalyst was used in a fixed bed reactor for upgrading pyrolysis vapors from an upstream fast pyrolysis reactor (ex-situ configuration) to produce CFP oil with 16.5% oxygen content. Subsequent co-hydrotreating of 10 vol% of CFP oil with 90 vol% straight run diesel (SRD) at the bench scale yielded greater than 95% carbon efficiency for the CFP oil portion of the feed. Light oxygenated coproducts were recovered from the off gases from the CFP process via adsorption and subsequent desorption. Aspen Plus modeling of purification of the desorbed stream using the UNIFAC property method showed that acetone and methyl-ethyl-ketone can be purified to chemicals grade and sold as coproducts to benefit the economics and reduce GHG emissions attributed to the fuel product. Some of the key risks and uncertainties associated with this process, and potential future work are mentioned in this report.

09 BIOMASS FUELS↗
NASA NTRS2016

CFD-Based Design of a Filming Injector for N+3 Combustors

An effort was undertaken to perform CFD analysis of fluid flow in Lean-Direct Injection (LDI) combustors with axial swirl-venturi elements coupled with a new fuel-filming injector design for next-generation N+3 combustors. The National Combustion Code (NCC) was used to perform non-reacting and two-phase reacting flow computations on a N+3 injector configuration, in a single-element and a five-element injector array. All computations were performed with a consistent approach towards mesh-generation, spray-, ignition- and kinetics-modeling with the NCC. Computational predictions of the aerodynamics of the injector were used to arrive at an optimal injector design that met effective area, aerodynamics, and fuel-air mixing criteria. LDI-3 emissions (EINOx, EICO and UHC) were compared with the previous generation LDI-2 combustor experimental data at representative engine cycle conditions.

combustion chambers↗
DOE OSTI2024

In-situ determination of strain during transient burst testing and the temperature dependence of Zircaloy-4 claddings

Understanding fuel system behavior during postulated loss-of-coolant accidents is pertinent for continued safe and efficient operation of light water reactors, particularly as higher burnups are being pursued and safety margins re-evaluated. Conventional mechanical models for the incumbent Zr alloys typically rely on the assumption that steady-state creep is the dominant fuel cladding response during transient accident conditions. To investigate this assumption, simulated accident burst testing was performed on Zircaloy-4 claddings with balloon behavior measured in-situ. Here, two distinct loading conditions were utilized during burst testing: (1) constant-gas-inventory where pressure was allowed to increase with temperature and (2) constant pressure. In-situ strains and strain rates were measured via 2-dimensional digital image correlation techniques and synchronized with temperature to determine deformation dependencies. The temperature dependence of strain rate was characterized by a two segment Arrhenius relationship, with a distinct transition between the high and low temperature/strain regimes. The average activation energy of the lower temperature/strain regime was 328 ± 25 kJ/mol, in agreement with the ~320 kJ/mol used for conventional LOCA models. However, the higher temperature/strain segment, which encompassed most of ballooning, showed increased activation energies as well as a dependence on whether the burst region was in view. For tests that burst away from the camera view, the average high temperature/strain segment activation energy was 635 ± 150 kJ/mol. For samples where the rupture opening formed in view, the average activation energy was 1015 ± 179 kJ/mol. This observed shift in temperature dependence indicates a transition in deformation mechanism at the end of life, possibly to time independent failure mechanisms, which has not yet been visualized in the literature for Zr alloys. Parameters at the transition points were analyzed to determine thresholds for this change in behavior, which occurred at an average hoop strain of 6.9 ± 2.1 %.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sensitivity and Uncertainty Quantification of Transition Scenario Simulations

This report documents the first collective attempt at developing and applying capabilities to quantify uncertainties, assess parametric sensitivities, and optimize multiple parameters and metrics in fuel cycle simulations generated by the SA&I Campaign. To do this, external codes that were designed to perform sensitivity analysis and uncertainty quantification (SA&UQ) needed to be coupled to the SA&I Campaign’s nuclear fuel cycle simulators (NFCS). In FY20, two approaches were pursued: 1) coupling Cyclus to an ORNL-internal code called MOT (Metaheuristic Optimization Tool) and 2) coupling DYMOND to the opensource SA&UQ tool kit Dakota. The primary objective of having these NFCS/SA&UQ coupled capabilities is to better inform DOE-NE and other stakeholders on the results generated from the NFCS. For a given set of fuel cycle strategies, policies, and technology assumptions that make up a fuel cycle scenario, these NFCS have traditionally been used by the SA&I Campaign to provide quantitative answers in terms of year-by-year mass flows, infrastructure requirements, costs, etc. With these newly developed coupled capabilities, the SA&I Campaign can now efficiently simulate hundreds or thousands of these scenarios, sample large ranges of parameters and assumptions, and use the unique features of the SA&UQ tools to process the data. This enables providing answers with known and propagated uncertainties, determining the sensitivity of important metrics to different parameters and assumptions, quantifying how much fuel cycle and technology parameters impact each other, and producing optimized fuel cycle strategies for single and multiple variables. To demonstrate these new capabilities, the Cyclus/MOT was used to model several scenarios ranging from simple fleet retirements to transitions to advanced reactors. Specifically, for a transition scenario from LWRs to SFRs and advanced LWRs, uncertainty quantification, sensitivity analysis, and optimization studies were applied to cases involving single and multiple parameter (input) and single and multiple metric (output) variations. In addition, a similar transition scenario was modeled to demonstrate how to optimize the reprocessing capacity parameter to minimize two performance metrics while taking into account uncertainties from two other parameters. Lastly, a depletion module based on SCALE/ORIGEN was added in Cyclus to simulate the third scenario that was designed to quantify the impact of the modeling assumption that all LWR used nuclear fuel have the same burnup. The newly developed DYMOND/Dakota capability was also applied to a transition scenario from the existing fleet to small modular reactors and fast reactors. This particular scenario involves not only explicit isotopic depletion via ORIGEN-2, but also includes multirecycling and utilizing the criticality search feature to determine the fresh fuel composition of recycled fuel, a feature unique to the DYMOND NFCS. A large database of simulations were run with 4 main parameters that were sampled: start date of reprocessing, reprocessing capacity, energy demand growth rate, and advanced reactor share of the fleet. The 4 main metrics were uranium consumption, enrichment requirements, waste generation, and levelized cost of electricity using data from the Cost Basis Report. The demonstrated SA&UQ results include those that inform on how to choose parameters to avoid “failed” scenarios, Sobol’ indices that inform on the importance of various parameters individually and synergistically, and Analysis of Variance (ANOVA) studies that decompose parameter ranges into groups and informs on whether variations are statistically significant.

Feng, B.↗
DOE OSTI2022

Predicting Performance of Macroalgae Farms with Hydrodynamic and Biological Modeling

As part of the ARPA-E MARINER program a set of simulation tools and findings were developed for the hydrodynamic, biological, and economic modeling of large scale offshore macroalgae farms. Results suggest the utility of the tools in understanding the complex interplay of design choices and environmental conditions on the structural loading and farm performance which drive the costs for macroalgae production.

09 BIOMASS FUELS↗
DOE OSTI2021

Evaluating Safeguards Statistical Assumptions via Stochastic Simulation

Herein, the authors built and tested a stochastic simulation to estimate achieved detection probabilities (DPs) on a stratum basis, over a tailorable range of diverted amounts from 0 to 2 SQ, using typical IAEA inspection data: i.e., SQ in stratum, number of items, number of gross/partial/bias defect measurements conducted, and realistic relative standard deviation (RSD) values for typical IAEA verification measurements. For bulk strata, the model calculates achieved DP at 0.01 SQ diversion increments; for item strata, the model calculates DP using the smallest realistic diversion increment (e.g., a plate, pin, or coupon). After successfully benchmarking against IAEA deterministic models, the simulation was used to test the sensitivity of DP to certain standard assumptions and selected input parameters. First, the equal defect assumption was tested; the results suggest significant complexity in the effectiveness of partial defect measurements. Next, the authors explored the sensitivity of DP to the assumed RSD of attribute tests. Then, the authors compared non-normal models for instrument performance (e.g., logistic, step, or arbitrary functions) to the typical results from a normal distribution (characterized by RSD). This last comparison was supplemented with experimentally derived performance data for an HM-5. The HM-5 was used to make enrichment measurements on both LEU and HEU MTR fuel elements as plates were removed, and the results fit with logistic and step curves and applied in the simulation. These stochastic DP results were compared to DP estimates from a deterministic model assuming a normal curve and typical RSD, yielding insights that could improve effectiveness in the field. These early results illustrate the potential of stochastic models to better understand achieved DP and to improve safeguards effectiveness.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

A sensitivity analysis to predict the neutronics behavior of samples irradiated in the VTR rabbit system

We report a low-order neutronics model is developed to carry out hundreds of simulations efficiently and investigate the neutronics behavior of samples being irradiated in a test reactor setting under different geometrical constraints. The low-order model allowed for simulations that yield the expected neutronics behavior of any irradiated sample in any environment and allows for the calculation of highly accurate spatially averaged statistics and idealized spatial distributions in the neutron flux. Several benchmarks are performed to evaluate the performance and limitations of the low-order model revealing many important findings. The low-order model predicted the LHGR in the EBR-II driver fuel to within 2.34% by only simulating the fuel rod by itself, which served as a validation for the model. Sensitivity studies investigated 3% enriched UO 2 and U-10Zr being irradiated in the Versatile Test Reactor rabbit system. The analyses investigated a range of combinations of 15 radii and 5 heights for each sample in the rabbit system. Similar data sets are also provided for irradiations in the Advanced Test Reactor’s B-10 irradiation position, which is a thermal neutron spectrum environment. Generalized fits and fit coefficients are obtained for sample heating, reaction rate densities, and local multiplication rate characteristics, allowing the predictions of the neutronics behavior of the samples based on their geometrical constraints. The analyses and fits laid the groundwork for developing a user-end Multiphysics analysis framework to assist and accelerate irradiation experiment design and optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Effect of non-solvent from the phase inversion method on the morphology and performance of the anode supported microtubular solid oxide fuel cells

The microstructure of the anode in anode-supported solid oxide fuel cells has significant influence on the cell performance. Here, microtubular Ni-yttria stabilized zircona (Zr 0.8 Y 0.2 O 2 , YSZ) anode support has been prepared by the phase inversion method. Different compositions of non-solvent have been used for the fabrication of the Ni-YSZ anode support, and the correlation between non-solvent composition and characteristics of the microstructure of the anode support has been investigated. The presence of ethanol or isopropanol in the non-solvent can inhibit the growth of the finger-like pores in the anode support. With the increase of the concentration of ethanol or isopropanol in the non-solvent, a thin dense layer can be observed on the top of the prepared tubular anode support. In addition, the mechanism of pore formation is explained based on the inter-diffusivity between the solvent and the non-solvent. Lastly, the prepared microtubular solid oxide fuel cells (MT-SOFCs) have been tested, and the influence of the anode microstructure on the cell electrochemical performance is analyzed based on a polarization model.

36 MATERIALS SCIENCE↗
NASA NTRS1974

JT8D-100 turbofan engine, phase 1

The JT8D turbofan engine, widely used in short and medium range transport aircraft, contributes substantially to airport community noise. The jet noise is predominant in the JT8D engine and may be reduced in a modified engine, without loss of thrust, by increasing the airflow to reduce jet velocity. A configuration study evaluated the effects of fan airflow, fan pressure ratio, and bypass ratio on noise, thrust, and fuel comsumption. The cycle selected for the modified engine was based upon an increased diameter, single-stage fan and two additional core engine compressor stages, which replace the existing two-stage fan. Modifications were also made to the low pressure turbine to provide the increased torque required by the larger diameter fan. The resultant JT8D-100 engine models have the following characteristics at take-off thrust, compared to the current JT8D engine: Airflow and bypass ratio are increased, and fan pressure ratio and engine speed are reduced. The resultant engine is also longer, larger in diameter, and heavier than the JT8D base model, but these latter changes are compensated by the increased thrust and decreased fuel comsumption of the modified engine, thus providing the capability for maintaining the performance of the current JT8D-powered aircraft.

Source record↗
NASA NTRS1983

Determination of the static performance of a Cuyuna cc model UL-430RR engine

The results of static performance tests were carried out on an ultralight engine. A Cuyuna 430 cc, 2 stroke, 2 cylinder model was used to measure thrust, torque, propeller RPM, fuel flow, cylinder head temperature, and change in pressure through the propeller. Measurement of each of the above parameters was taken at specific values of RPM. The propeller's RPM's ranged from idle at approximately 750 RPM to a maximum value of 2810 RPM. The test results were then manipulated to obtain thrust coefficient, power coefficient, shaft horse power, and shaft specific fuel consumption.

Roskam, J.↗
NASA NTRS1990

Calculations of gaseous H2/O2 thruster

Calculations were made for the gaseous H2/O2 thruster for the Space Station by using the RPLUS code which employs an implicit finite volume, LUSSOR scheme to solve the Navier-Stokes equations and the species equations. The combustion processes of hydrogen and oxygen are modeled by a 9 species and 18 step reaction mechanism and the turbulence is simulated by the Baldwin-Lomax turbulence model for the thruster wall boundary layer and the modified Prandtl's mixing length model for the reacting shear layer. Results are presented for different mixture ratios and fuel film cooling percents and compared with the experimental data. The calculated performance predictions for the thruster agree well with the experimental data and the results demonstrate that the RPLUS code can be used for design and analysis of thrusters and rockets.

Kim, S. C.↗