Search NASA⌕ Search

SEARCH · Search NASA

Results for “EVAPORATION MODEL”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

THE QUASI-ROCHE LOBE OVERFLOW STATE IN THE EVOLUTION OF CLOSE BINARY SYSTEMS CONTAINING A RADIO PULSAR

We study the evolution of close binary systems formed by a normal (solar composition), intermediate-mass-donor star together with a neutron star. We consider models including irradiation feedback and evaporation. These nonstandard ingredients deeply modify the mass-transfer stages of these binaries. While models that neglect irradiation feedback undergo continuous, long-standing mass-transfer episodes, models including these effects suffer a number of cycles of mass transfer and detachment. During mass transfer, the systems should reveal themselves as low-mass X-ray binaries (LMXBs), whereas when they are detached they behave as binary radio pulsars. We show that at these stages irradiated models are in a Roche lobe overflow (RLOF) state or in a quasi-RLOF state. Quasi-RLOF stars have radii slightly smaller than their Roche lobes. Remarkably, these conditions are attained for an orbital period as well as donor mass values in the range corresponding to a family of binary radio pulsars known as ''redbacks''. Thus, redback companions should be quasi-RLOF stars. We show that the characteristics of the redback system PSR J1723-2837 are accounted for by these models. In each mass-transfer cycle these systems should switch from LMXB to binary radio pulsar states with a timescale of approximately one million years. However, there is recent and fast growing evidence of systems switching on far shorter, human timescales. This should be related to instabilities in the accretion disk surrounding the neutron star and/or radio ejection, still to be included in the model having the quasi-RLOF state as a general condition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Evaporation Induced Entrainment of Contaminants from Evaporating and Burning Liquids.

This interim report details model development, theory, and a literature review focusing on the evaporation induced entrainment (sub-boiling) of contaminated liquids. Entrainment from a variety of sources is the topic of DOE Handbook 3010, and this report deals more broadly with fire related airborne sources of contaminants in hazardous operations. Relatively few studies have examined sub-boiling behavior in the past, however, it can be a phenomenon that presents a fire related risk under hazardous operations. Molecular dynamics simulations are used to infer the gaseous evolution of coordinated complexes, and a model for a water/plutonium/nitrate system is deduced from the simulation results by evaluating the statistical trends of the results. Questions remain as to the chemical reactivity and longevity of entrained species. A generalized computer model capability and simple analytical model assumptions are developed for predicting the results of these and other (boiling and solid entrainment) scenarios. Verification related predictions using these models are illustrated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Accelerated Carbon and Water Cycles in the Amazon and Congo Basins Revealed From TRENDY Models and Remote Sensing Products

Tropical forests play a vital role in the global carbon cycle and land–atmosphere interactions. Estimating tropical forest carbon–water dynamics is challenging due to observational and modeling uncertainties. This study leverages the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project models and satellite observations to assess changes (2003–2021) in vegetation carbon, gross primary production (GPP), evapotranspiration (ET), and net biosphere production (NBP) in the Amazon and Congo. Atmospheric CO 2 , climate variability, and land use and land cover changes constrain these variables between 1700 and 2021 with the overall increasing trends of carbon stock and fluxes. The models overestimate vegetation carbon and GPP, while ET and NBP are consistent with observations. Fire-activated models predict lower values for vegetation carbon and GPP, ET, and NBP, aligning more closely with observations. The higher ET from fire-activated models may result from enhanced soil evaporation due to increased canopy openings. Fire-inactivated models could well estimate the magnitudes of NBP. The high vegetation carbon in nitrogen-enabled models points to simulation uncertainties and imbalance in model numbers regarding the nitrogen cycle. Although the nitrogen cycle enhances water use efficiency in both the Amazon and Congo, the models show a higher sensitivity to the nitrogen cycle in the Congo. This study highlights the challenges in accurately representing tropical biogeochemical cycles and the values of satellite products in model evaluations, underscoring the need for standard modeling protocols that address biogeochemical components (e.g., nutrient cycles) to better resolve process-based representations.

Shi, Mingjie [Pacific Northwest National Laborator↗
DOE OSTI2015

OPTIMAL ELECTRON ENERGIES FOR DRIVING CHROMOSPHERIC EVAPORATION IN SOLAR FLARES

In the standard model of solar flares, energy deposition by a beam of electrons drives strong chromospheric evaporation leading to a significantly denser corona and much brighter emission across the spectrum. Chromospheric evaporation was examined in great detail by Fisher et al., who described a distinction between two different regimes, termed explosive and gentle evaporation. In this work, we examine the importance of electron energy and stopping depths on the two regimes and on the atmospheric response. We find that with explosive evaporation, the atmospheric response does not depend strongly on electron energy. In the case of gentle evaporation, lower energy electrons are significantly more efficient at heating the atmosphere and driving up-flows sooner than higher energy electrons. We also find that the threshold between explosive and gentle evaporation is not fixed at a given beam energy flux, but also depends strongly on the electron energy and duration of heating. Further, at low electron energies, a much weaker beam flux is required to drive explosive evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Generalized Spallation Model

The Generalized Spallation Model (GSM) event generator is an advanced, robust, and scalable Fortran code developed from an extended and improved version of the earlier event generator codes CEM03.03 and LAQGSM03.03. The GSM implements its algorithms within the context of cascade exciton, particle equilibration, particle evaporation, and particle disintegration models sampled via Monte-Carlo techniques. Coalescence of particles is also considered. The GSM assumes the incident particles has sufficient energy to satisfy the assumptions utilized in the GSM. The GSM calculates reactions induced by hadrons and nuclei incident on nuclear targets at incident energies ranging from ~20 MeV up to several TeV per nucleon. Usage of GSM is recommended however for particles with incident energies exceeding ~100 MeV to ensure model-assumptions are satisfied. The GSM additionally provides its clients with sub-model libraries to allow them to model desired physics within their own physical models.

Juneau, Chase↗
DOE OSTI2024

A Semi‐Mechanistic Model for Partitioning Evapotranspiration Reveals Transpiration Dominates the Water Flux in Drylands

Abstract Popular evapotranspiration (ET) partitioning methods make assumptions that might not be well‐suited to dryland ecosystems, such as high sensitivity of plant water‐use efficiency (WUE) to vapor pressure deficit (VPD). Our objectives were to (a) create an ET partitioning model that can produce fine‐scale estimates of transpiration (T) in drylands, and (b) use this approach to evaluate how climate controls T and WUE across ecosystem types and timescales along a dryland aridity gradient. We developed a novel, semi‐mechanistic ET partitioning method using a Bayesian approach that constrains abiotic evaporation using process‐based models, and loosely constrains time‐varying WUE within an autoregressive framework. We used this method to estimate daily T and weekly WUE across seven dryland ecosystem types and found that T dominates ET across the aridity gradient. Then, we applied cross‐wavelet coherence analysis to evaluate the temporal coherence between focal response variables (WUE and T/ET) and environmental variables. At yearly scales, we found that WUE at less arid, higher elevation sites was primarily limited by atmospheric moisture demand, and WUE at more arid, lower elevation sites was primarily limited by moisture supply. At sub‐yearly timescales, WUE and VPD were sporadically correlated. Hence, ecosystem‐scale dryland WUE is not always sensitive to changes in VPD at short timescales, despite this being a common assumption in many ET partitioning models. This new ET partitioning method can be used in dryland ecosystems to better understand how climate influences physically and biologically driven water fluxes.

Environmental Sciences & Ecology↗
DOE OSTI2016

NSRD-11: Computational Capability to Substantiate DOE-HDBK-3010 Data.

Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE - HDBK - 3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this year, this research included improvements of the high-fidelity codes to model particle resuspension and multi-component evaporation for fire scenarios. We also began to model ceramic fragmentation experiments, and to reanalyze the liquid fire and powder release experiments that were done last year. The results show that the added physics better describes the fragmentation phenomena. Thus, this work provides a low-cost method to establish physics-justified safety bounds by taking into account specific geometries and conditions that may not have been previously measured and/or are too costly to perform.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Capturing the Page curve and entanglement dynamics of black holes in quantum computers

Quantum computers are emerging technologies expected to become important tools for exploring various aspects of fundamental physics in the future. Therefore, we pose the question of whether quantum computers can help us to study the Page curve and the black hole information dynamics, which has been a key focus in fundamental physics. In this regard, we rigorously examine the qubit transport model, a toy qubit model of black hole evaporation on IBM’s superconducting quantum computers, to shed light on this question. Specifically, we implement the quantum simulation of the scrambling dynamics in black holes using an efficient random unitary circuit. Furthermore, we employ the swap-based many-body interference protocol and the randomized measurement protocol to measure the entanglement entropy of Hawking radiation qubits in this model. Finally, by incorporating quantum error mitigation techniques into our challenging implementation of entanglement entropy measurement protocols on the IBM quantum hardware, we accurately determine the Rényi entropy in the qubit transport model, thus showcasing the utility of quantum computers for future investigations of complex quantum systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Three‐Moment Representation of Rain in a Bulk Microphysics Model

Abstract A bulk three‐moment representation for rain microphysics is developed and implemented in the Predicted Particle Properties (P3) microphysics scheme. In addition, a new parameterization for rain self‐collection and collisional breakup (RSCB) is presented using a lookup table approach, based on the Spectral‐Bin Model (SBM). To quantify the impacts of sedimentation, evaporation, and RSCB on drop size distributions (DSDs), a rain shaft model is applied to a wide range of atmospheric scenarios (i.e., initial conditions and regimes) and compared against results from the SBM. DSD shapes are mainly determined by both sedimentation and evaporation, except in heavy rain where the impact of RSCB on DSD shape becomes more important than evaporation. The new parameterization for RSCB has a considerable impact on the mean drop size, improving the agreement between P3 and SBM. Only 4% of the original two‐moment rainshaft simulations have mean drop sizes and rain rates within ±20% of the SBM results, but this increases to more than 95% agreement when the three‐moment rain representation is used together with the new parameterization for RSCB. Generally, the improvement is more significant for heavy rain than for light drizzle. Remaining differences between bin and bulk model are attributable to treatments of evaporation, and the restriction to gamma DSDs in P3.

precipitation↗
DOE OSTI2022

Grey-box Fault Models and Applications for Low Carbon Emission CO2 Refrigeration System

CO2 refrigeration systems have attracted more attentions in recent years because of their zero-carbon emission advantages compared with traditional refrigerants. These systems are widely used in certain commercial buildings such as supermarkets. However, very limited energy studies have been done for such refrigeration systems under faults. Furthermore, this study investigated five common faults for supermarket transcritical CO2 refrigeration systems with display cases. The following major tasks were implemented for this study. (1) Five fault models (gray-box format) regarding the power consumption and supply air temperature were proposed for low- and high-temperature display cases, respectively. The models cover: ice accumulation on evaporator coil, evaporator fan partial failure, expansion valve failure, display door open, and condenser blockage. (2) These fault models were calibrated with field tests, and the modeling accuracies were in good agreement with measurement datasets. (3) The proposed fault models demonstrated more than 90% of fault detection for Fault Detection and Diagnostics (FDD) purposes for both field measurement and simulated data.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Gray-Box Fault Models and Applications for Low-Carbon Emission CO 2 Refrigeration Systems

CO 2 refrigeration systems have attracted more attentions in recent years because of their zero-carbon emission advantages compared with traditional refrigerants. These systems are widely used in certain commercial buildings such as supermarkets. However, very limited energy studies have been done for such refrigeration systems under faults. This study investigated five common faults for supermarket transcritical CO 2 refrigeration systems with display cases. The following major tasks were implemented for this study. (1) Five fault models (gray-box format) regarding the power consumption and supply air temperature were proposed for low- and high-temperature display cases, respectively. The models cover: ice accumulation on evaporator coil, evaporator fan partial failure, expansion valve failure, display door open, and condenser blockage. (2) These fault models were calibrated with field tests, and the modeling accuracies were in good agreement with measurement datasets. (3) The proposed fault models demonstrated more than 90% of fault detection for Fault Detection and Diagnostics (FDD) purposes for both field measurement and simulated data.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Understanding the causes of satellite–model discrepancies in aerosol–cloud interactions using near-LES simulations of marine boundary layer clouds

Aerosol–cloud interactions (ACI) remain the largest source of uncertainty in model estimates of anthropogenic radiative forcing, primarily because of deficiencies in representing aerosol–cloud microphysical processes that lead to inconsistent cloud liquid water path (LWP) responses to aerosol perturbations between observations and models. To investigate this discrepancy, we conducted a series of large-eddy-scale simulations driven by realistic meteorology over the eastern North Atlantic, and evaluated LWP susceptibility, precipitation processes, and boundary layer thermodynamics using satellite and ground-based observations. Simulated LWP responses show a strong dependence on cloud state. Non-precipitating thin clouds exhibit a modest LWP decrease with increasing cloud droplet number concentration (N d ), consistent in sign but weaker in magnitude than satellite estimates, reflecting enhanced turbulent mixing and evaporation. The largest model-observation discrepancy occurs in non-precipitating thick clouds, where simulated LWP susceptibilities are strongly positive (+0.32) while observations indicate large negative values (−0.69). This discrepancy stems from excessive precipitation driven by underestimated entrainment, overly active accretion, and overly broad drop-size distributions in polluted conditions. While our high-resolution setup mitigates the excessive drizzling common in coarser models and captures key regime transitions, these biases persist – highlighting that improved parameterizations of cloud-top processes, precipitation, and aerosol effects are needed beyond simply increasing model resolution. Additionally, misrepresented moisture inversions in reanalysis introduce a moist bias in cloud-top relative humidity, further amplifying positive LWP susceptibility. Our results also suggest that large negative N d –LWP relationships in observations may reflect internal cloud processes rather than true ACI effects.

Aerosol-cloud interaction↗
DOE OSTI2021

Minimizing Cr-Evaporation from Balance Of Plant Components By Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al 2 O 3 -forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Minimizing Cr-Evaporation from Balance of Plant Components by Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al2O3-forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Absorption of ethanol, acetone, benzene and 1,2-dichloroethane through human skin in vitro: a test of diffusion model predictions

The overall goal of this research was to further develop and improve an existing skin diffusion model by experimentally confirming the predicted absorption rates of topically-applied volatile organic compounds (VOCs) based on their physicochemical properties, the skin surface temperature, and the wind velocity. In vitro human skin permeation of two hydrophilic solvents (acetone and ethanol) and two lipophilic solvents (benzene and 1,2-dichloroethane) was studied in Franz cells placed in a fume hood. Four doses of each {sup 14}C-radiolabed compound were tested — 5, 10, 20, and 40 μL cm{sup −2}, corresponding to specific doses ranging in mass from 5.0 to 63 mg cm{sup −2}. The maximum percentage of radiolabel absorbed into the receptor solutions for all test conditions was 0.3%. Although the absolute absorption of each solvent increased with dose, percentage absorption decreased. This decrease was consistent with the concept of a stratum corneum deposition region, which traps small amounts of solvent in the upper skin layers, decreasing the evaporation rate. The diffusion model satisfactorily described the cumulative absorption of ethanol; however, values for the other VOCs were underpredicted in a manner related to their ability to disrupt or solubilize skin lipids. In order to more closely describe the permeation data, significant increases in the stratum corneum/water partition coefficients, K{sub sc}, and modest changes to the diffusion coefficients, D{sub sc}, were required. The analysis provided strong evidence for both skin swelling and barrier disruption by VOCs, even by the minute amounts absorbed under these in vitro test conditions. - Highlights: • Human skin absorption of small doses of VOCs was measured in vitro in a fume hood. • The VOCs tested were ethanol, acetone, benzene and 1,2-dichloroethane. • Fraction of dose absorbed for all compounds at all doses tested was less than 0.3%. • The more aggressive VOCs absorbed at higher levels than diffusion model predictions. • We conclude that even small exposures to VOCs temporarily alter skin permeability.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

4-(Trimethylsilyl)morpholine: synthesis, characterization, and prospects of use in film deposition processes

4-(Trimethylsilyl)morpholine (TMSM) was synthesized and shown to be applicable as a precursor for plasma-enhanced chemical vapor deposition of films. The temperature depen dence of the saturated vapor pressure of TMSM was determined by tensimetry and the thermodynamic characteristics of its evaporation were calculated. Thermodynamic modeling (calculation of the 'equilibrium' compositions of the condensed and gas phases) was carried out for broad ranges of temperatures (300—1200 K), pressures (10{sup –2}—10 Torr), and input gas (inert gas, H{sub 2}, NH{sub 3}) to precursor flow ratios (0—50). The calculation results were used to predict the composition of deposited condensed products occurring in equilibrium with the gas phase depending on the composition of the initial gas mixture and reaction conditions.

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

Distillation-based Droplet Modeling of Non-Ideal Oxygenated Gasoline Blends: Investigating the Role of Droplet Evaporation on PM Emissions

In some studies, a relationship has been observed between increasing ethanol content in gasoline and increased particulate matter (PM) emissions from vehicles equipped with spark ignition engines. The fundamental cause of the PM increase seen for moderate ethanol concentrations is not well understood. Ethanol features a greater heat of vaporization (HOV) than gasoline and also influences vaporization by altering the liquid and vapor composition throughout the distillation process. A droplet vaporization model was developed to explore ethanol's effect on the evaporation of aromatic compounds known to be PM precursors. The evolving droplet composition is modeled as a distillation process, with non-ideal interactions between oxygenates and hydrocarbons accounted for using UNIFAC group contribution theory. Predicted composition and distillation curves were validated by experiments. Detailed hydrocarbon analysis was applied to fuel samples and to distillate fractions, and used as input for the initial droplet composition. With composition calculated throughout the distillation, the changing HOV and other physical properties can be found using reference data. The droplet can thus be modeled in terms of energy transfer, which in turn provides the transient mass transfer, droplet temperature, and droplet diameter. Model predictions suggest that non-ideal vapor-liquid equilibrium along with an increase in HOV can alter the droplet composition evolution. Results predict that the presence of ethanol causes enrichment of the higher boiling fractions (T90+) in the aromatic components as well as lengthens the droplet lifetime. A simulation of the evaporation process in a transient environment as experienced within an engine cylinder predicts a decrease in mixing time of the heaviest fractions of the fuel prior to spark initiation, possibly explaining observations linking ethanol to PM.

02 PETROLEUM↗