Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fractionalization”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase equilibria and microstructure evolution in the Ni-rich region of binary Ni-Ce and ternary Ni-Ce-X (X=Al, Nb, Cr, Ti) systems at 900 °C

The equilibrium compositions of all phases in the Ni-rich region of binary Ni-Ce and ternary Ni-Ce-X (X=Al, Nb, Cr, Ti) systems at 900 °C were determined experimentally through isothermal annealing for up to 500 h and compared with currently available computational thermodynamic descriptions. Additionally, the study examined how adding a third element to the binary system affects the microstructure by analyzing changes in the volume fraction of the intermetallic phase in the eutectic region. It was observed that Nb and Ti addition lowered the intermetallic fraction in the eutectic region while Al and Cr addition did the opposite. The fraction of intermetallic phase in the eutectic region is very important because it contributes to extreme brittleness of these alloys. Also, Al, Nb and Ti addition promoted different precipitates upon annealing when added in sufficient amounts. The collected data on composition and phase equilibria provides crucial information for improving the thermodynamic assessments of the binary and ternary alloy systems.

energy dispersive spectrometry↗
DOE OSTI2021

Insights into Defect-Mediated Nucleation of the Equilibrium B2 Phase in Face-Centered Cubic High-Entropy Alloys

A fine scale second phase distribution is important for strengthening metallic materials. High precipitate density is often difficult to achieve due to sluggish precipitation kinetics and limited nucleation sites. Based on the CALPHAD predictions, transition element-based complex concentrated alloys (CCAs)/ High Entropy Alloys (HEAs), have a highly stable face-centered cubic (FCC) phase and intermetallic B2 phase. The precipitation of the strengthening B2 phase is often limited to low fractions in these alloys due to a high nucleation barrier in the FCC matrix. Here, by illustrating in three different alloys, we present a defect mediated nucleation of a high fraction equilibrium B2 phase. The degree of cold work is used to manipulate the phase fraction, morphology, and distribution of B2 precipitates. Our approach synergistically improves the strength and ductility of these alloys and is widely applicable to many other complex alloys.

Sharma, Abhishek↗
DOE OSTI2024

Phase Identification in Synchrotron X-ray Diffraction Patterns of Ti–6Al–4V Using Computer Vision and Deep Learning

X-ray diffraction patterns contain information about the atomistic structure and microstructure (defect population) of materials, extracting detailed information from diffraction patterns is complex, demanding and relies on prior knowledge. Here, we hypothesize that deep-learning techniques can help to perform an effective and accurate analysis with high throughput rates. To demonstrate this concept, we applied a novel deep learning framework to determine the evolution of the β-phase volume fraction in a Ti–6Al–4V alloy during heat-treatment from video sequences of 2D diffraction patterns recorded in transmission and with highly monochromatic radiation in a synchrotron beamline. In particular, we studied the impact of network design on prediction reliability and computational performance. Networks of different architectures were trained using 3008 experimental 2D patterns. A well-tuned model was found to reproduce the phase fractions of another experimental data set, consisting of 1100 diffraction patterns, with a mean-square error as small as 2.6 x 10 -4 . The average prediction error of β-phase volume fraction was within 1.6 x 10 -2 (in each diffraction pattern) of the values obtained by conventional methods. Our work demonstrates that convolutional neural networks can evaluate high energy X-ray diffraction patterns with a remarkable level of reliability. Furthermore, it demonstrates the significance of network design on the reliability of predictions and computational performance. The most complex models do not necessarily result in highest accuracy and may even fail to learn from the data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pitting behavior of friction stir repair-welded 304L stainless steel in 3.5% NaCl solution at room temperature: role of grain and defect structures

Stainless steel (SS) canisters used for storing spent nuclear waste are prone to chloride-induced stress corrosion cracking. Friction stir welding (FSW), being a low heat input process, can be used for repair welding of these canisters. In this study, an artificial crack was introduced on the 304L SS coupons by an electric discharge machining process. The artificial crack was repair-welded by FSW by maintaining the tool temperature a constant at two levels (725 and 825 °C). Friction stirring significantly reduced the grain size of the stirred zone (SZ) from about 47 to 2–4 μm. The fraction of low-angle grain boundaries increased in the SZ from 2 to 37–43%. Here, on the other hand, the fraction of special grain boundaries (Σ3 and Σ9) that was ~ 50% in the base metal reduced to < 10% in the SZ. During friction stirring, the oxide layer of artificial crack was broken and aligned to form a spiral defect called lazy-S structure. All these microstructural changes affected the corrosion behavior of the FSW specimens when tested in 3.5% NaCl at room temperature using cyclic polarization, chronoamperometry, and electrochemical impedance spectroscopy. The FSW specimens showed lower polarization resistance and lower transpassive potential than those of the base metal specimens. However, the pitting protection potentials of the FSW specimens were higher than that of the base metal. The pitting behavior of the FSW specimens was influenced more by the preferential attack on the lazy-S region than by passive film breakdown. The flat band potentials of the passive film formed on the FSW specimens were more positive than that of base metal. The charge carrier density of passive film formed on the FSW specimens was higher than that of the base metal. The higher fraction of low-angle grain boundaries present in the FSW specimens could supply more number of misorientation dislocations at the metal/film interface which could form anion vacancies by a climb process leading to formation of oxide at these locations without stressing the substrate. Therefore, low-angle boundaries are considered helpful for formation of stress-free passive films that will be highly stable and enhance both pitting and stress corrosion cracking resistance.

Dry storage canister↗
DOE OSTI2022

Cavitation-resistant intergranular precipitates enhance creep performance of θ'-strengthened Al-Cu based alloys

Tensile and compressive creep properties of a quaternary Al-Cu-Mn-Zr (ACMZ) alloy and its commercial counterpart (Al-Cu-Mn-Zr with Ni, Co and Sb additions, RR350) are investigated at 300°C. At low stresses up to 30 MPa where diffusional creep dominates, creep resistance is the same in tension and compression and RR350 deforms more slowly than ACMZ, consistent with RR350 alloy's larger linear fraction of intergranular precipitates (Al7Cu2(NiFe) and Al9FeNi for RR350 vs. θ-Al2Cu for ACMZ) and a reduced fraction of precipitate-free zones near grain boundaries. At stresses between 30 and 80 MPa, dislocation creep with a stress exponent n ~ 3 becomes rate-limiting in compression, which is expected to be controlled by θ' precipitates within the grain bulk. By contrast, in tension, enhanced creep rate and higher apparent stress exponents are measured, consistent with cavitation at intergranular precipitates becoming increasingly dominant as the stress increases. In the dislocation creep regime, RR350 alloy is again more creep resistant than ACMZ alloy, which is related to three mechanisms (i) a reduced fraction of softer precipitate-free zones, (ii) more effective load transfer to intergranular precipitates, and (iii) reduced cavitation. A model for cavitation is applied to calculate tensile creep rates from compressive creep rates and the model successfully predicts the improved tensile creep resistance of the RR350 alloy. Overall, the present investigation underscores the importance of intergranular grain boundary precipitates, in addition to strengthening θ' precipitates, in enhancing the creep resistance of Al-Cu alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atmospheric humidity deficits tell us how soil moisture deficits down-regulate ecosystem evaporation

Knowing how actual evaporation is down-regulated from potential evaporation during periods with soil moisture deficits is one of the greatest challenges towards computing evaporation everywhere on a regular basis. We propose the hypothesis that vegetated landscapes transmit information on soil moisture deficits through its down-regulation of evaporation, which in turn affects the humidification and growth of the planetary boundary layer. To test this hypothesis, here, we examined how the evaporative fraction, defined as the ratio between actual and potential evaporation, corresponds with the soil moisture stress index, defined as the relative humidity with a power law dependence on vapor pressure deficit times a parameter, β. We tested the parameterized soil moisture stress index with direct eddy covariance measurements of actual evaporation and computations of potential evaporation based on meteorological conditions, averaged on monthly time steps. The analysis was conducted on a dataset obtained from 144 FLUXNET sites. These sites spanned much of the world's climates and biomes and contained over 5900 months of observations. Observations of the evaporative fraction fit the model of the soil moisture stress index best for semi-arid and arid ecosystems, and least for the humid and wet ecosystems. Consequently, a significant relationship between measurements of the evaporative fraction and soil moisture stress index held for about one-half of the population of sites. Under this condition, the median value of the parameter, β, was 1.41. To investigate the mechanism of this empirical soil moisture stress index, we diagnosed it with a coupled evaporation-planetary boundary layer model. The soil moisture stress index is strongly related to surface resistance, as defined by inverting Penman-Monteith equation. Consequently, this index provides an independent estimate of surface resistance based on easy to measure mean monthly weather conditions like relative humidity and temperature. Thereby, this soil moisture stress index has potential to be applied to weather, climate and biogeochemical models and the interpretation of satellite derived evaporation products, like the one provided by the ECOSTRESS mission.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation of two-phase gas transport into VERA for molten salt reactor analysis

Molten salt reactors (MSRs) are a class of next-generation nuclear reactors that have received recent industrial and research interest. A generalized species transport solver was implemented in the Virtual Environment for Reactor Applications (VERA) computing suite to extend this tool to analyze liquid-fueled MSRs. This core simulator has been extended to model the transport of fission product gases into a collection of circulating gas bubbles with the purpose of removing the gases. This paper presents the governing species transport equation, along with various nuclear source terms. Development of the source term for phase migration is discussed, along with a simplified interfacial area tracking method. Finally, a case study on a simplified MSR loop is presented in which modeling parameters were varied to assess their impact on gas removal. The steady state results show that parameters such as bubble diameter, gas injection rate and mass transfer coefficient have a low to moderate effect on the fraction of xenon in the core region. Removal efficiency has the greatest effect on the fraction in the core region. After the pump bowl, bubble diameter has a minor effect on the fraction of xenon in the gas void. Finally, these results point out that increasing parameters such as mass transfer coefficient, gas injection rate, and removal efficiency drives the xenon into the circulating gas void, while decreasing bubble diameter also drives xenon into the gas void by increasing interfacial area.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Simulation-based methodology to assess the lattice defects creation as energy storing process

In nuclear industry, damages caused to in-core structural materials by high-energy neutrons have been studied for years. Recently, the possibility of storing energy through the creation of defects in crystalline materials via heavy particle bombardment was investigated. To evaluate the potential of the Wigner effect as energy-storing process, a graphite-moderated Molten Salt Reactor was adopted as a test case. Neutronics simulations performed with OpenMC Monte Carlo code allowed estimating the amount of energy released during fission events, the fraction that is deposited in graphite, and the fraction that has the potential to cause radiation damages. Additionally, molecular dynamics simulations were then performed to study the evolution of cascades triggered by high-energy neutrons. As a major outcome, a more accurate estimate of the radiation damage with respect to the traditional approaches was obtained, and the fraction of the deposited neutron energy that can be stored through long-lived defects was quantified.

25 ENERGY STORAGE↗
DOE OSTI2022

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (V d ) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases V d for submicron particles and decreases it for larger particles by an average of 37% and –66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%–97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM 10 by 41% and positive fractional biases for PM 2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Utilization of guaiacol-based deep eutectic solvent for achieving a sustainable biorefinery

This study proposed a renewable deep eutectic solvent (DES) pretreatment using lignin-derived guaiacol as the hydrogen bond donor. The DES showed excellent biomass fractionation efficiency after the incorporation of trace AlCl 3 as the reinforcer, which removed 79.1 % lignin while preserving more than 90 % glucan. The pretreated bamboo exhibited 96.2 % glucan enzymatic hydrolysis yield at only 110 °C. The physicochemical properties of the pretreated solids were comprehensively investigated to explain how the DES fractionation overcame the biomass recalcitrance. The regenerated lignin from the DES pretreatment was also analyzed, which revealed that lignin β-O-4 bond was significantly cleaved. In conclusion, this guaiacol-based DES could greatly contribute to establish a closed-loop biorefinery sequence with high lignin fractionation efficiency and great solvent recyclability.

09 BIOMASS FUELS↗
DOE OSTI2024

Catalytic conversion of post-consumer recycled high-density polyethylene oil over Zn-impregnated ZSM-5 catalysts

The catalytic conversion of post-consumer recycled high-density polyethylene (PCR-HDPE) pyrolysis oil was investigated using modified ZSM-5 zeolites, prepared under different conditions. The PCR-HDPE oil, obtained through pyrolysis at 500 °C followed by distillation at 165 °C, contained a naphtha fraction hydrocarbons suitable for catalytic upgrading. H/ZSM-5 converts PCR-HDPE naphtha fraction oil primarily into BTX, along with light paraffins and olefins. Introducing mesoporosity led to increased propylene yield at the expense of BTX. Zn impregnation onto H/ZSM-5 and H-meso/ZSM-5 enhanced BTX yield by up to 13 %. However, Zn-meso/ZSM-5 exhibited a lower BTX yield than Zn/ZSM-5 with the same Zn loadings, especially after additional steam treatment at 700 °C. Spent catalysts revealed no noticeable increase in metal contaminants such as calcium and aluminum that were present in the PCR-HDPE. Here, the Zn/ZSM-5 catalyst was stable for more than 50 h time on stream and showed promising potential for upgrading PCR-HDPE naphtha fraction oils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The simultaneous removal of technetium and iodine from Hanford tank waste

The simultaneous removal of radionuclides technetium-99 and iodine-129 from an actual decontaminated Hanford tank waste sample (a mixture of decontaminated waste from tanks 241-AP-105 and 241-AP-107) was demonstrated for the first time in this work. A series of commercially available ion exchange resins were evaluated in batch contact tests in the tank waste, and all showed removal of both Tc and I. The highest Tc removal was observed for Purolite A530e while the highest iodine removal was observed for ResinTech SIR-110-MP. Batch tests in simulated tank waste with these two resins showed that the SIR-110-HP-MP had consistently higher K d for both pertechnetate and iodide and much higher K d than previous works on Tc removal from Hanford waste. As such, the SIR-110-MP was evaluated in a dual -column (lead/lag) test processing 5.2L of the tank waste mixture showing 60% breakthrough of Tc on the lead column and no significant breakthrough on the lag after 625 bed volumes (BV, 6 mL size) while significant iodine breakthrough (>50%) occurred after 28 BV. The limited iodine uptake was attributed to the column conditions generating mass transfer limitations. A fraction of the Tc and I was not captured by the resin (<10%) in either the batch tests or column tests. The iodine fraction was identified to be an iodide, likely organo-iodide. The fraction of the Tc was identified as a non-pertechnetate species, which is the first time non-pertechnetate has been identified in AP-105 and AP-107 tanks, although the exact species is still unknown.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗
DOE OSTI2022

Industrial scale mesofluidic particle separation

Here we explore the potential of a mesofluidic separator that fits within industrial piping and is designed to remove large aspherical particles out of an abrasive slurry at fast flow rates. Mesofluidic separators have recently demonstrated separation of spherical particles at flow rates of ~1 gallon/min (gpm) (6.3.10 -5 m 3 /s). Yet, the performance of these separators at flow rates achievable within industrial piping have not been explored. Here we evaluate the permeate flowrate fraction, large particle mass fraction, and pressure drop from approximately 35-100 gpm (2.2.10-3-6.3.10-3 m 3 /s) in a three inch tube using an abrasive and modestly non-spherical slurry. Faster flow rates may be achieved with larger pumps and larger diameters. We find the permeate fraction and pressure drop both increase with flow rate, though the pressure drop remains modest. Only two feet of inserts removed 97.6±2.1% of large particles for conditions tested. Separation was demonstrated in both vertical and horizontal orientations with the horizontal orientation benefiting marginally from gravitational settling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Separation of rare earth element radioisotopes by reverse-phase high-speed counter-current chromatography

Analytical scale purification of rare earth element (REE) radioisotopes is typically accomplished using cation-exchange resins (e.g. AG 50W-X8) and high-performance liquid chromatography (HPLC). Despite the variety of improvements made since the development of this separation process in the 1950s, nearest neighbor separations remain a challenge, as does the issue of irreversible sample adsorption. Herein, we report a study that evaluates the potential of high-speed counter-current chromatography (HSCCC) as an alternative method for purifying REE elements, with specific reference to separations of fission product REE of interest to nuclear forensics. Complementary HSCCC REE separation experiments, one spiked with radiotracer and REE fission product activity, allowed for in depth analysis of resulting fractions from both an elemental (inductively coupled plasma atomic emission spectroscopy, ICP-AES) and radiological (gamma-ray spectrometry, beta counting) purity perspective. The highly reproducible nature of separation profiles generated from HSCCC instruments was leveraged to simplify work-up of samples containing radioisotopes. Subsequent radioanalytical evaluation revealed minimal carryover of Eu into neighboring Sm and Tb fractions (as indicated by presence of 150Eu), and trace contamination of the Tb fraction with Y (as indicated by presence of 91Y). Subtle differences in stationary phase retention across the two columns were reflected in significant variations in decontamination factors of duplicate parallel separations. Furthermore, these differences paired with obtained distribution of radioisotopes provided valuable insights into future improvements. Collectively, this study represents a significant step forward in development of HSCCC technology for task specific REE radioisotope purification.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

A two-equation soot-in-flamelet modeling approach applied under Spray A conditions

Soot production (including formation and oxidation) is studied in the transient, high-pressure and turbulent n-dodecane Spray A flames from the Engine Combustion Network (ECN) using computational fluid dynamics (CFD) simulations. A two-equation soot-in-flamelet modeling approach is applied within the framework of the Unsteady Flamelet Progress Variable (UFPV) model and results are validated against experimental data. Equations for soot mass fraction and soot number density derived in the mixture fraction space are solved in the context of detailed flamelet calculations. Source terms for the different steps in the soot chemistry are tabulated and incorporated in the flamelet manifold. For the reference condition, the modeling approach based on the tabulated flamelet manifold reduces the computational cost of a CFD calculation by approximately 40 times compared to a non-tabulated well-mixed (WM) modeling approach. The soot-in-flamelet approach is then extended to study the effect of ambient oxygen concentration, ambient mixture composition and ambient temperature on soot production. Results show that the modeling approach is able to capture the experimental trends for the soot volume fraction (SVF) with good quantitative agreement, especially in the soot ramp-up region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗