Search NASA⌕ Search

SEARCH · Search NASA

Results for “EVAPORATION RATE”

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 127 records · Page 7

Ex-Situ Recrystallization of CIGS via Metal Halides

Post Deposition Treatments (PDT) have been key to improving the efficiency of CIGS solar cells. While efficiency has continued to increase, for thermal evaporation processes, high substrate temperature and low deposition rates are a necessity. For many chemical processes, chemical vapor transport techniques can be used to both purify and grow both large and pure crystalline compounds. This is often done with a halide species and temperature. Here, by using the substituent metals of CIGS in confluence with a halide, we explored the possibility to grow a CIGS layer at low temperature and fast rate while still having sufficiently high efficiency.

14 SOLAR ENERGY↗

Control Oriented Model of Cabin-HVAC System in a Long-Haul Trucks for Energy Management Applications

Super Truck II is a 48V mild hybrid class 8 truck with an all auxiliary loads powered purely by the battery pack. Electric Heating Ventilation and Air Conditioning (HVAC) load is the most prominent battery load during the hotel period, when the truck driver is resting inside the sleeper. For the PACCAR Super Truck II (ST-II) project a 48 V battery system provides the required power during the hotel period. A cabin-HVAC model estimates the electric load on the 48V battery system, allowing the control system to implement an efficient energy management strategy that avoids engine idling during the hotel period. The thermal model accounts for the sun load due to the time of day and the geographic location of the truck during the hotel period. The cabin-HVAC model has two parts. First, a grey box model with two heat exchangers (Condenser and Evaporator) working in unison with refrigerant mass flow rate as an input and HVAC load as an output. Second, a two-node cabin model formulated to estimate the cabin temperature as a function of the Global Horizontal Irradiance (GHI), HVAC load and ambient temperature. The models are calibrated using experimental cabin-HVAC system data as for long-haul class 8 truck (e.g. ST-II). Here, the model simulations show that the overall Root Mean Square Error (RMSE) value of 0.4°C between the experimental and simulated cabin temperature.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Effect of Ice Number Concentration on the Evolution of Boundary Layer Clouds During Arctic Marine Cold‐Air Outbreaks

Abstract Marine cold‐air outbreaks (MCAOs) are crucial for Arctic Ocean heat loss, featuring convective cloud rolls that transition into convection cells downstream. Understanding factors controlling this transformation is the key for improving MCAO cloud representation in climate models. This study employs large‐eddy simulations to investigate how cloud ice number concentrations () affect cloud evolution using a case from the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) campaign. The simulations, performed in a Lagrangian framework following an air mass trajectory, are driven by ERA5 reanalysis data. Initially, all simulations produce similar cloud patterns, but higher leads to earlier breakup of cloud rolls. Between 4 and 10 hr, surface precipitation rates are similar across simulations, but precipitation initiates earlier, and the cloud‐base precipitation rates are higher when is higher. The stronger precipitation evaporation leads to increased stability of the boundary layer and reduced intensity of vertical mixing between the surface and cloud layer. An increased sink of cloud layer moisture via precipitation and decreased source through diminished vertical transport result in earlier cloud breakup in higher conditions. Simulations with different sea surface temperatures (SST) indicate that this cloud breakup mechanism remains valid for MCAOs of different strengths, although the cloud organization is more sensitive to SST changes in low environments. This work highlights the importance of accurate representations of ice processes in simulating MCAO clouds and suggests the need for observational constraints of ice nucleating particles and over the mixed‐phase cloud regimes.

54 ENVIRONMENTAL SCIENCES↗

Abstract for CRADA between Capture6 CORP and NETL

The National Energy Technology Laboratory (NETL) and Capture6 CORP (Participant) will collaborate to refine Capture6’s transformative air contactor design for optimal direct air capture (DAC) under an awarded project from the Department of Energy’s Direct Air Capture Pre-Commercialization Technology Prize. Current cooling tower designs require exponentially more resources as the concentration of CO 2 in the gas phase decreases, resulting in a low conversion rate of carbonates, a large physical footprint, and high evaporation of water, which consequently limits the energy and water benefits that a DAC process can create. Capture6’s innovative tower design AURA (Absorption Unit Reactor for Ambient air) seeks to improve the mass transfer characteristics of legacy cooling towers and reduces total energy demand by: (1) introducing laminar and turbulent flow to increase mass transfer, (2) generating a solvent gradient to improve CO 2 absorption, (3) leveraging venturi and packed towers to increase the surface for the mass-transfer process, (4) adding charged/cool plates to reduce overall water loss in the system, and (5) utilizing surfactant additives to increase mass transfer. The collaboration will facilitate the deployment of technologies that reduce the total energy required and carbon footprint of solvent-based carbon dioxide removal (CDR) technologies. The anticipated results support NETL’s goal of fostering rapid development and commercialization pathways for technologies that remove CO 2 from the atmosphere.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simple device modification simultaneously enhances indoor passive evaporative cooling power and duration

Passive evaporative cooling from a damp cloth or vessel can be used for localized cooling of perishables such as food and medicine, drinking water, and even people. Such applications have two key objectives, which are normally in tension with each other: minimizing the water consumption rate and maximizing the cooling (steady-state temperature swing between the cool cloth and the warm surroundings). Here, we present a simple device that delivers higher performance in both metrics simultaneously by adding a perforated aluminum foil sheet suspended over a stagnant air layer to thermally shield the evaporating surface while also facilitating mass transfer. Experimental results demonstrate a simultaneous 10%–25% increase in temperature swing and 2–3× reduction in water consumption rate as compared to conventional evaporative cooling. This improvement is achieved through careful modeling to optimize the coupled heat and mass transfer through the airgap (thermally insulating and vapor permeable) and perforated foil (infrared reflective and vapor permeable).

Chen, Sarah Mizuno↗

Thin hybrid capillary two-phase cooling system

A novel hybrid two-phase cooling system was developed that integrated a mechanically pumped two-phase loop with a capillary-driven two-phase cooling device. The latter cooling mechanism was based on evaporation/boiling from wick structures made by sintering copper particles on the interior surfaces of a copper cold plate. The cold plate provided cooling to two surfaces and each of them included four heaters in series. The novelty of the developed technology was preventing flooding of the evaporator wicks by isolating the evaporation surface from the pumped liquid flow that fed them. This arrangement allowed for a high liquid feed flow rate much greater than would be allowed by a capillary pumped system while maintaining a low thermal resistance at the evaporation surface. Using this approach, the cooling system removed over 850 W with a low pumping power below 1.0 W while using R245fa as the working fluid. The equivalent heat fluxes exceeded 970 W/cm 2 over areas less than 0.12 cm 2 . The measured thermal resistance was as low as 0.09 K-cm 2 /W. The presented thermal management solution enables an increase in the power of high heat flux electronic devices beyond the state-of-the-art.

42 ENGINEERING↗

Antifoam Development for Eliminating Flammability Hazards and Decreasing Cycle Time in the Defense Waste Processing Facility

The Savannah River National Laboratory (SRNL) was requested to develop a new antifoam control method for the Defense Waste Processing Facility’s (DWPF) Chemical Process Cell (CPC). SRNL completed testing of both chemical and nonchemical foam controls. The nonchemical foam controls were either ineffective (or worse, created more foam) or impractical (a water spray can control foam, but excessive water is needed). As a result, the focus of this study was on finding a superwetter or commercial antifoam for controlling foam. Thirty potential antifoams were tested as part of this study. A series of tests were developed to help screen out ineffective alternatives including: 1. Spreading testing of superspreaders, 2. Foam column testing with physical simulants, 3. Boiling testing with physical and chemical simulants, 4. Days-only Sludge Receipt and Adjustment Tank (SRAT) process simulations with sludge (containing noble metals and mercury), Precipitate Reactor Feed Tank (PRFT), and Slurry Mix Evaporator Feed Tank (SEFT) simulants in the RC1 Reaction Calorimeter (purchased for antifoam testing), and 5. Around-the-clock SRAT and Slurry Mix Evaporator (SME) process simulations with sludge(containing noble metals and mercury), PRFT, and SEFT simulants in the RC1 Reaction Calorimeter. Evonik Surfynol® MD20, a commercially available defoamer, was relatively effective in controlling foam, while remaining chemically stable in SRAT and SME processing across the pH range of 4 to 13. No degradation products were detected in the offgas, in the condensate or in the SRAT and SME products. In nitric-glycolic acid flowsheet testing, 250 mg/kg Evonik Surfynol® MD20 was needed for foam control compared to 1,625 mg/kg for Antifoam 747, DWPF’s current antifoam. In nitric-formic acid flowsheet testing, 1,125 mg/kg of Evonik Surfynol® MD20 was needed to control foam throughout the SRAT and SME cycles. The commercially available superspreader Momentive™ Y-17112 was even more effective than Evonik Surfynol® MD20 as both a defoamer and an antifoam. Not only was the foam destroyed upon addition but also was less persistent between additions. It was the most effective antifoam in testing using both the nitric-glycolic acid flowsheet and the nitric-formic acid flowsheet. In nitric-glycolic acid flowsheet testing, only 100 mg/kg Momentive™ Y-17112 was needed to control foam throughout the SRAT and SME cycles. In nitric-formic acid flowsheet testing, 300 mg/kg Momentive™ Y-17112 was needed to control foam throughout the SRAT and SME cycles. Momentive™ Y-17112 is also resistant to hydrolysis as demonstrated by its chemical stability in SRAT and SME processing across the pH range of 4 to 13 and lack of degradation products in offgas or condensate. Both candidates were effective as potential replacements for Antifoam 747, with Y-17112 demonstrating superior foam control. During nitric-glycolic flowsheet testing 50% less antifoam was needed when using Momentive™ Y-17112 compared to MD20. During nitric-formic flowsheet testing 75% less antifoam was needed when using Momentive™ Y-17112 compared to MD20. Foam remediated with Momentive™ Y-17112 was less persistent throughout testing. In addition, no degradation products were detected in the offgas, in the condensate or in the SRAT and SME products. Based on this testing, Momentive™ Y-17112 is clearly superior to Evonik Surfynol® MD20 and Antifoam 747, especially for the nitric-formic acid flowsheet processing; it is recommended that Momentive™ Y-17112 replace Antifoam 747 in DWPF. An antifoam addition strategy is recommended for both the nitric-glycolic acid flowsheet and the nitric-formic acid flowsheet. Implementation of Momentive™ Y-17112 is expected to decrease SRAT and SME boiling times by up to 50%, eliminate the issues resulting from Antifoam 747 degradation products, and minimize foamovers. To validate the effectiveness of these defoaming agents, SRNL recommends irradiation of a SRAT or SME product simulant containing fresh antifoam. The goal of this testing is to determine whether the irradiation causes decomposition of the antifoam that would make it less effective or produce new species in the offgas or slurry. This testing began in April 2020. An evaluation should be completed to determine the thermolytic hydrogen and methane generation rate in downstream equipment, including the High-Level Waste evaporators.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Sensitivities of Large Eddy Simulations of Aerosol Plume Transport and Cloud Response

Abstract Cloud responses to surface‐based sources of aerosol perturbation partially depend on how turbulent transport of the aerosol to cloud base affects the spatial and temporal distribution of aerosol. Here, scenarios of plume injection below a marine stratocumulus cloud are modeled using large eddy simulations coupled to a prognostic bulk aerosol and cloud microphysics scheme. Both passive plumes, consisting of an inert tracer, and active plumes are investigated, where the latter are representative of saltwater droplet plumes such as have been proposed for marine cloud brightening. Passive plume scenarios show higher in‐plume cloud brightness (relative to out‐of‐plume) due to the predominant transport of the passive plume tracer from the near‐surface to the cloud layer within updrafts. These updrafts rise into brighter areas within the cloud deck, even in the absence of an aerosol perturbation associated with an active plume. Comparing albedo at in‐plume to out‐of‐plume locations associates the inert plume with the brightest cloud locations, without any causal effect of the plume on the cloud. Numerical sensitivities are first assessed to establish a suitable model configuration. Then sensitivity to particle injection rate is investigated. Trade‐offs are identified between the number of injected particles and the suppressive effect of droplet evaporation on plume loft and spread. Furthermore, as the near‐field in‐plume brightening effect does not depend significantly on injection rate given a suitable definition of perturbed versus unperturbed regions of the flow, plume area is a key controlling factor on the overall cloud brightening effect of an aerosol perturbation.

54 ENVIRONMENTAL SCIENCES↗

The evaluation of aluminum and iron metal oxide settling behaviors for Hanford insoluble solids waste preprocessing

The Hanford site is currently one of the largest and most expensive cleanup sites for hazardous, radioactive waste. Over 20% of the waste found at the Hanford site is in the form of a high activity sludge. The insoluble solids in the sludge will need to be concentrated prior to vitrification in the high-level waste (HLW) melter. This will minimize the amount of liquid that will be evaporated during the melting process and expedite the melter processing rate. One proposed option for concentrating the insoluble solids is gravity settling in the storage tanks. Metal oxide compounds containing aluminum and iron make up the majority of the insoluble solids in the sludge, therefore understanding the behavior of these compounds in various tank waste matrices can facilitate sludge pretreatment options. A study of non-radioactive slurry solutions containing Al(OH) 3 (gibbsite), AlO(OH) (boehmite), and Fe 2 O 3 (iron (III) oxide) was conducted to determine the time dependent interface behavior and settling rates of these compounds. Variations in solids loading and sodium concentration were evaluated to represent waste processing conditions and the results of these settling studies were compared with prior tank waste settling tests. Information gathered from these studies can be used to inform future decisions on sludge treatment processes of the insoluble solids processed at the Hanford site.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Condensation and evaporation of boson stars

Abstract Axion-like particles, including the QCD axion, are well-motivated dark matter candidates. Numerical simulations have revealed coherent soliton configurations, also known as boson stars, in the centers of axion halos. We study evolution of axion solitons immersed into a gas of axion waves with Maxwellian velocity distribution. Combining analytical approach with controlled numerical simulations we find that heavy solitons grow by condensation of axions from the gas, while light solitons evaporate. We deduce the parametric dependence of the soliton growth/evaporation rate and show that it is proportional to the rate of the kinetic relaxation in the gas. The proportionality coefficient is controlled by the product of the soliton radius and the typical gas momentum or, equivalently, the ratio of the gas and soliton virial temperatures. We discuss the asymptotics of the rate when this parameter is large or small.

Physics↗

Sludge Batch 10 Flowsheet Testing with Non-radioactive Simulants

Fourteen Chemical Processing Cell (CPC) simulations were performed with nonradioactive sludge simulants at the Aiken County Technology Laboratory in Aiken, SC. Four of these experiments were performed with Tank 51 sludge simulant. The remaining ten were performed with Tank 40 sludge simulant. The purpose of these experiments was to elucidate the chemistry and characteristics of Sludge Batch (SB) 10 as anticipated in the Defense Waste Processing Facility (DWPF). Experiments were performed at acid stoichiometries between 76% and 138% of the Koopman Minimum Acid requirement (85% - 144% of the Hsu acid requirement) and at REDuction/OXidation (REDOX) targets between 0.1 and 0.3. Testing examined the impact of coupled operations and sludge-only operations during Sludge Receipt and Adjustment Tank (SRAT) and Slurry Mix Evaporator (SME) processing at both design basis and nominal boilup rates. This report shares conclusions made as a result of this testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Improving tubular protonic ceramic fuel cell performance by compensating Ba evaporation via a Ba-excess optimized proton conducting electrolyte synthesis strategy

Protonic ceramic fuel cells (PCFCs) are emerging as a promising technology for reduced temperature ceramic energy conversion devices. The BaCe 0.4 Zr 0.4 Y 0.1 Yb 0.1 O 3–δ (BCZYYb4411) electrolyte is notable for its high proton conductivity. However, the tendency of barium to volatilize in BCZYYb4411 during high-temperature sintering compromises its chemical stability and performance. This study investigates the effects of intentionally incorporating excess barium into BCZYYb4411, formulated as Ba 1+x Ce 0.4 Zr 0.4 Y0.1Yb 0.1 O 3–δ (where x = 0, 0.1, 0.2, and 0.3), with the aim of compensating barium evaporation and enhancing the physical and chemical properties. We find that excess barium results in a greater shrinkage rate, facilitating a denser electrolyte structure. This barium-enriched electrolyte demonstrates improved electrochemical performance by effectively counteracting the deleterious effects of barium evaporation. Applying this strategy to tubular PCFCs, we achieved a peak power density of 480 mW•cm –2 at 600 °C. This unique approach provides a simple, tunable, and easy-to-implement processing modification to achieve high-performance tubular PCFC.

25 ENERGY STORAGE↗

Evaluation of Thermolytic Hydrogen Generation Rate Models at High-Temperature/High-Hydroxide Regimes

This report describes the results of testing performed to extend the applicable ranges of temperature and hydroxide concentration for use within the Glycolate and Global Total Organic Carbon (TOC) Hydrogen Generation Rate (HGR) expressions. Seven experimental conditions (six simulants of the 242-25H Evaporator system chosen as a D-optimal set of experiments and a single test conducted at an elevated boiling point of 170 °C) were investigated in the presence of sodium glycolate and Xiameter TM AFE-1010. Glycolate was employed to study the extension of the Glycolate Thermolytic HGR expression while Xiameter TM AFE-1010 was employed to study the extension of the Global TOC Thermolytic HGR expression. The following conclusions were derived from this testing: The Glycolate Thermolytic HGR expression may be confidently used to predict thermolytic HGRs from glycolate at temperatures as high as 170 °C and hydroxide concentrations as high as 23 M.; The hydroxide and temperature-dependence predicted by the Global TOC Thermolytic HGR expression has been confirmed at temperatures as high as 170 °C and hydroxide concentrations as high as 23 M, suggesting that the Global TOC Thermolytic HGR expression may be used at these ranges.; Methane was observed from tests with Xiameter TM AFE-1010 at production rates higher than those observed for hydrogen. These rates were observed at temperatures higher than 100 °C.; Preliminary models suggest that increasing hydroxide/temperature causes an increase in Methane Generation Rate (MGR) from Xiameter TM AFE-1010. The following recommendations are based on this testing: The existing equations for thermolytic HGR from glycolate and non-glycolate organics should be used at Concentration, Storage, and Transfer Facilities (CSTF) storage and evaporation conditions, including temperatures and hydroxide concentrations exhibited in the 242-25H Evaporator.; Further investigation should be made into the influence of methylsilanes on CSTF flammability. This investigation should include: determination of the types of methylsilanes historically added to the CSTF, determination of methane formation rates from each type of methylsilane, and determination of the extent of degradation of methylsilanes in CSTF waste.; Characterization techniques should be developed by Savannah River National Laboratory (SRNL) to assist in the speciation of methylsilane-containing waste in the CSTF.; Additional testing with radioactive waste should be performed to determine the MGRs possible in radioactive waste and better inform model predictions made from testing with simulants.

08 HYDROGEN↗

Determination of the 60 Zn level density from neutron evaporation spectra

Nuclear reactions of interest for astrophysics and applications often rely on statistical model calculations for nuclear reaction rates, particularly for nuclei far from β stability. However, statistical model parameters are often poorly constrained, where experimental constraints are particularly sparse for exotic nuclides. For example, our understanding of the breakout from the NiCu cycle in the astrophysical rp-process is currently limited by uncertainties in the statistical properties of the proton-rich nucleus 60 Zn. We have determined the nuclear level density of 60 Zn using neutron evaporation spectra from 58 Ni( 3 He,n) measured at the Edwards Accelerator Laboratory. We compare our results to a number of theoretical predictions, including phenomenological, microscopic, and shell-model-based approaches. Notably, we find the 60 Zn level density is somewhat lower than expected for excitation energies populated in the 59 Cu(p,γ) 60 Zn reaction under rp-process conditions. This includes a level density plateau from roughly 5 to 6 MeV excitation energy, which is counter to the usual expectation of exponential growth and all theoretical predictions that we explore. Here, a determination of the spin distribution at the relevant excitation energies in 60 Zn is needed to confirm that the Hauser-Feshbach formalism is appropriate for the 59 Cu(p,γ) 60 Zn reaction rate at x-ray burst temperatures

59 ≤ A ≤ 89↗

Effect of Metal Halides Treatment on High Throughput Low Temperature CIGS Solar Cells

Copper indium gallium diselenide (CIGS) semiconductor thin films were deposited at high rate and low temperature using single-stage thermal co-evaporation process on molybdenum back contact. A post deposition treatment was done by flashing AgBr at 350 ºC to induce recrystallization. Changes in morphology were confirmed by SEM, with an observed increase in grain size, as well as by XRD measurements, with a decrease in FWHM. Device results show an improvement of the performance after the AgBr vapor treatment, as all the photovoltaic parameters enhanced. Altogether, AgBr seems to be a suitable transport agent and beneficial for device fabrication.

14 SOLAR ENERGY↗

Refueling Infrastructure Scoping and Feasibility Assessment for Hydrogen Rail Applications

The feasibility and component cost of hydrogen rail refueling infrastructure is examined. Example reference stations can inform future studies on components and systems specifically for hydrogen rail refueling facilities. All of the 5 designs considered assumed the bulk storage of liquid hydrogen on-site, from which either gaseous or liquid hydrogen would be dispensed. The first design was estimated to refuel 10 multiple unit trains per day, each train containing 260 kg of gaseous hydrogen at 350 bar on-board. The second base design targeted the refueling of 50 passenger locomotives, each with 400 kg of gaseous hydrogen on-board at 350 bar. Variations from this basic design were made to consider the effect of two different filling times, two different hydrogen compression methods, and two different station design approaches. For each design variation, components were sized, approximate costs were estimated for major components, and physical layouts were created. For both gaseous hydrogen-dispensing base designs, the design of direct-fill using a cryopump design was the lowest cost due to the high cost of the cascade storage system and gas compressor. The last three base designs all assumed that liquid hydrogen was dispensed into tender cars for freight locomotives that required 7,500 kg of liquid hydrogen, and the three different designs assumed that 5, 50, or 200 tender cars were refueled every day. The total component costs are very different for each design, because each design has a very different dispensing capacity. The total component cost for these three designs are driven by the cost of the liquid hydrogen tank; additionally, delivering that much liquid hydrogen to the refueling facility may not be practical. Many of the designs needed the use of multiple evaporators, compressors, and cryopumps operating in parallel to meet required flow rates. In the future, the components identified here can be improved and scaled-up to better fit the needs of heavy-duty refueling facilities. This study provides basic feasibility and first-order design guidance for hydrogen refueling facilities serving emerging rail applications.

08 HYDROGEN↗

Reduction of the neutron imaginary potential off the stability line and its possible impact on neutron capture rates

The effect of the isovector imaginary optical potential has been studied experimentally by using the particle-evaporation technique for the 11 B + 48 Ca reaction with a 21.8 MeV 11 B beam. Spectra of neutron, proton, and α particles emitted from the neutron-rich compound nucleus 59 Mn have been measured and analyzed with traditional optical model potentials with their original parametrizations as well as with adjusted isovector imaginary components. It is shown that the isovector component of the imaginary potential is indispensable in the reproduction of proton and α-particle yields from this reaction and even needs to be enhanced compared with the suggestions of the original model parametrizations. Furthermore, this can lead to important consequences for astrophysical reaction-rate calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗