Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOILING”

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,009 records · Page 56

Listeria monocytogenes loss of cultivability on carrot is associated with the formation of mesosome-like structures

Raw carrot is known to have antimicrobial activity against Listeria monocytogenes, but the mechanism of action has not been fully elucidated. In this study, we examined carrot antilisterial activity against several strains of Listeria species (including L. grayi, L. innocua, L. seeligeri, and L. welshimeri) and L. monocytogenes. A representative strain of L. monocytogenes was subsequently used for further characterizing carrot antilisterial activity. Exposure to fresh-cut carrot for 15 min resulted in a similar loss of cultivability, ranging from 2.5 to 4.7 log units, across all Listeria strains evaluated. L. monocytogenes recovered from the fresh-cut surface of different raw carrots was 1.6 to 4.1 log lower than levels obtained from paired boiled carrot samples with abolished antilisterial activity. L. monocytogenes levels recovered from fresh-cut carrot were 2.8 to 3.1 log lower when enumerated by culture-dependent methods than by the culture-independent method of PMAxx-qPCR, a qPCR assay that is performed using DNA pre-treated to selectively sequester DNA from cells with injured membranes. These results suggested that L. monocytogenes loss of cultivability on fresh-cut carrot was not associated with a loss of L. monocytogenes cell membrane integrity and putative cell viability. Transmission electron microscopy imaging revealed that L. monocytogenes rapidly formed mesosome-like structures upon exposure to carrot fresh-cut surface but not upon exposure to boiled carrot surface, suggesting there may be an association between the formation of these mesosome-like structures and a loss of cultivability in L. monocytogenes. As a result, further research is necessary to conclude the causality of this association.

59 BASIC BIOLOGICAL SCIENCES↗

Computational examination of two-phase microchannel heat transfer correlations with conjugate heat spreading

Microchannel flow boiling is an attractive thermal management strategy for ever-growing volumetric heat dissipation demands associated with electronic systems. Due to difficulties related to measurement at the microscales, the majority of researchers have chosen to study relatively simple situations with uniform heat flux. However, many applications involve local hotspots which give rise to highly non-uniform heat flux and temperature gradients due to heat spreading. This necessitates the consideration of conjugate heat transfer for accurate analysis. The current work is aimed at investigating conjugate heat transfer in a two-phase microchannel array. Experimental data was collected on R134a flow boiling heat transfer for very small hydraulic diameter (<100 µm) silicon microchannels with very high heat fluxes (>1 kW cm -2 ) applied via platinum strip heaters with a footprint size of 1 cm × 1 mm. The collected experimental data was then combined with detailed computational modeling utilizing finite element modeling with COMSOL Multiphysics and MATLAB to examine the applicability of five published heat transfer correlations for use in determining local heat transfer coefficients. A two-phase correlation from Agostini and Bontemps, developed for markedly different test parameters, provided the best computational agreement with an RMS temperature difference from experiment of 3.3 °C and a predicted peak perimeter-averaged heat transfer coefficient of 116 MW m -2 K -1 . Modeling confirms the presence of highly non-uniform local heat flux and correspondingly non-uniform local heat transfer coefficient. Finally, the results of this study make clear the need for better micro-scale two-phase correlations developed to predict local heat transfer coefficients at these small scales and high local heat fluxes.

42 ENGINEERING↗

Spray process of multi-component gasoline surrogate fuel under ECN Spray G conditions

As modern gasoline direct injection (GDI) engines utilize sophisticated injection strategies, a detailed understanding of the air-fuel mixing process is crucial to further improvements in engine emission and fuel economy. In this study, a comprehensive evaluation of the spray process of single-component iso-octane (IC8) and multi-component gasoline surrogate E00 (36 % n-pentane, 46 % iso-octane, and 18 % n-undecane, by volume) fuels was conducted using an Engine Combustion Network (ECN) Spray G injector. High-speed extinction, schlieren, and microscopy imaging campaigns were carried out under engine-like ambient conditions in a spray vessel. Experimental results including liquid/vapor penetration, local liquid volume fraction, droplet size, and projected liquid film on the nozzle tip were compared under ECN G1 (573 K, 3.5 kg/m 3 ), G2 (333 K, 0.5 kg/m 3 ), and G3 (333 K, 1.01 kg/m 3 ) conditions. In addition to the experiments, preferential evaporation process of the E00 fuel was elucidated by Large–Eddy Simulations (LES). The three-dimensional liquid volume fraction measurement enabled by the computed tomographic reconstruction showed substantial plume collapse for E00 under the G2 and G3 conditions having wider plume growth and plume-to-plume interaction due to the fuel high vapor pressure. Here, the CFD simulation of E00 showed an inhomogeneity in the way fuel components vaporized, with more volatile components carried downstream in the spray after the end of injection. The high vapor pressure of E00 also results in ~4 μm smaller average droplet diameter than IC8, reflecting a higher rate of initial vaporization even though the final boiling point temperature is higher. Consistent with high vapor pressure, E00 had a wider plume cone angle and enhanced interaction with the wall to cover the entire surface of the nozzle tip in a film. However, the liquid fuel underwent faster evaporation, so the final projected tip wetting area was smaller than the IC8 under the flash-boiling condition.

33 ADVANCED PROPULSION SYSTEMS↗

Modeling of hydrogen liquefaction using magnetocaloric cycles with permanent magnets

Hydrogen (H 2 ) is promising alternative to replace fossil fuels, but its transport and storage has been challenging. As H 2 fuel cell vehicles are gaining traction, the infrastructure for storing large amounts of liquid H 2 is needed. However, liquid H 2 would suffer from boil-off loss, and traditional vapor compression refrigeration systems would not be able to economically recover the lost H 2 due to the low efficiencies at cryogenic temperature. Magnetocaloric (MC) refrigeration systems could possess much higher coefficient of performance (COP) at cryogenic temperature compared to the vapor compression ones. Previous work on cryogenic MC systems, however, have only focused on large scale applications which use superconducting magnets to provide a large magnetic field but are prohibitively expensive to operate for small scale applications, such as that of a H 2 refilling station. In this work, we model the performance of a MC refrigeration cycle using 1-Tesla permanent magnets for H 2 liquefaction, with the objective of cooling H 2 from 80 K (using liquid nitrogen as the heat sink) to 20 K (boiling point of hydrogen). We evaluate main performance metrics including the total work input to the refrigeration system, COP, total MCM mass in the system, and total volume of the permanent magnets, etc. Our modeling results indicate that such a permanent magnet-based MC cooling system is feasible for small-scale H 2 liquefaction, with projected COP values significantly higher than those of vapor compression systems. In conclusion, this work provides design guidelines for future experimental efforts on permanent magnet MC cooling systems for cryogenic cooling.

08 HYDROGEN↗

Dynamics within geyser conduits, and sensitivity to environmental perturbations: Insights from a periodic geyser in the El Tatio geyser field, Atacama Desert, Chile

Despite more than 200. years of scientific study, the internal dynamics of geyser systems remain poorly characterized. As a consequence, there remain fundamental questions about what processes initiate and terminate eruptions, and where eruptions begin. Over a one-week period in October 2012, we collected down-hole measurements of pressure and temperature in the conduit of an exceptionally regular geyser (132. s/cycle) located in the Chilean desert. Here, we identified four stages in the geyser cycle: (1) recharge of water into the conduit after an eruption, driven by the pressure difference between water in the conduit and in a deeper reservoir; (2) a pre-eruptive stage that follows the recharge and is dominated by addition of steam from below; (3) the eruption, which occurs by rapid boiling of a large mass of water at the top of the water column, and decompression that propagates boiling conditions downward; and (4) a relaxation stage during which pressure and temperature decrease until conditions preceding the recharge stage are restored. Eruptions are triggered by the episodic addition of steam coming from depth, suggesting that the dynamics of the eruptions are dominated by geometrical and thermodynamic complexities in the conduit and reservoir. Further evidence favoring the dominance of internal processes in controlling periodicity is also provided by the absence of responses of the geyser to environmental perturbations (air pressure, temperature and probably also Earth tides).

Bubble trap↗

A conceptual geochemical model of the geothermal system at Surprise Valley, CA

Characterizing the geothermal system at Surprise Valley (SV), northeastern California, is important for determining the sustainability of the energy resource, and mitigating hazards associated with hydrothermal eruptions that last occurred in 1951. Previous geochemical studies of the area attempted to reconcile different hot spring compositions on the western and eastern sides of the valley using scenarios of dilution, equilibration at low temperatures, surface evaporation, and differences in rock type along flow paths. These models were primarily supported using classical geothermometry methods, and generally assumed that fluids in the Lake City mud volcano area on the western side of the valley best reflect the composition of a deep geothermal fluid. In this contribution, we address controls on hot spring compositions using a different suite of geochemical tools, including optimized multicomponent geochemistry (GeoT) models, hot spring fluid major and trace element measurements, mineralogical observations, and stable isotope measurements of hot spring fluids and precipitated carbonates. We synthesize the results into a conceptual geochemical model of the Surprise Valley geothermal system, and show that high-temperature (quartz, Na/K, Na/K/Ca) classical geothermometers fail to predict maximum subsurface temperatures because fluids re-equilibrated at progressively lower temperatures during outflow, including in the Lake City area. We propose a model where hot spring fluids originate as a mixture between a deep thermal brine and modern meteoric fluids, with a seasonally variable mixing ratio. The deep brine has deuterium values at least 3 to 4‰ lighter than any known groundwater or high-elevation snow previously measured in and adjacent to SV, suggesting it was recharged during the Pleistocene when meteoric fluids had lower deuterium values. The deuterium values and compositional characteristics of the deep brine have only been identified in thermal springs and groundwater samples collected in proximity to structures that transmit thermal fluids, suggesting the brine may be thermal in nature. On the western side of the valley at the Lake City mud volcano, the deep brine-meteoric water mixture subsequently boils in the shallow subsurface, precipitates calcite, and re-equilibrates at about 130 °C. On the eastern side of the valley, meteoric fluid mixes to a greater extent with the deep brine, cools conductively without boiling, and the composition is modified as dissolved elements are sequestered by secondary minerals that form along the cooling and outflow path at temperatures <130 °C. Re-equilibration of geothermal fluids at lower temperatures during outflow explains why subsurface temperature estimates based on classical geothermometry methods are highly variable, and fail to agree with temperature estimates based on dissolved sulfate-oxygen isotopes and results of classical and multicomponent geothermometry applied to reconstructed deep well fluids. Furthermore, the proposed model is compatible with the idea suggested by others that thermal fluids on the western and eastern side of the valley have a common source, and supports the hypothesis that low temperature re-equilibration during west to east flow is the major control on hot spring fluid compositions, rather than dilution, evaporation, or differences in rock type.

15 GEOTHERMAL ENERGY↗

Pyroclast cooling and saturation in water

In the submarine setting pyroclast cooling and saturation, by controlling the buoyancy of individual clasts and surrounding fluid, exert fundamental controls on the dispersal and fate of volcanic material. Here, we use laboratory experiments to understand and quantify how hot and air-filled pyroclasts cool and saturate in liquid water. By measuring internal temperature and clast submerged weight in over thirty experiments, we find that pumice cools and saturates in two stages. We observe that the first stage is marked by much higher rates of cooling and saturation than the second stage and that the transition between stages is sharp. We propose conceptual and quantitative models for each stage of pyroclast cooling and saturation. Specifically we argue that the first stage is defined by the presence of internal steam. That is, hot pyroclasts ingest liquid water and generate steam. Once a clast's internal pore space is filled with steam, heat loss to the surroundings causes condensation of internal vapor and clast saturation. Stage 1 ends when all internal steam has condensed. We test this idea with a Stefan Model with advection and conclude that heat loss, rather than permeability, controls the saturation of pyroclasts above the boiling temperature. We find that permeability matters for air-filled clasts with initial temperatures below the boiling point. Stage 2 begins after all internal vapor has condensed and we propose that stage 2 cooling and saturation are controlled by heat conduction and thermal contraction of trapped gas, respectively. From our understanding of the processes that govern the stage 1 to 2 transition, we derive an empirical average pyroclast cooling rate of q = 7.5 ± 0.5 W cm -2 . Finally, we develop quantitative models for cooling-controlled porous clast saturation and for buoyant clast rise to the ocean surface. We find that meter size clasts stay buoyant for more than 10 min - long enough to reach the surface from a depth of 1000 m. Lastly, the models developed here test our understanding of clast-scale processes and demonstrate that - in hot and porous clasts - the details of internal texture do not matter because clast saturation, and thus buoyancy and fate in the water column, are governed by heat transfer.

58 GEOSCIENCES↗

Optimizing the impacts of solid additives on the operational stability and processing reliability of organic solar cells

Previous reports have revealed that by leveraging solid additives, organic solar cells (OSCs) can surpass the device’s performance beyond the intrinsic limitations of host photoactive molecules, a remarkable advancement. However, the impacts of more complex interactions introduced by solid additives are not yet well understood. Herein, optimizing the fabrication process based on the traditional efficiency-guided approach fails to represent the ideal and most practical devices. In particular, achieving superior operational stability while minimizing the device performance scattering was found to require processing solvent evaporation to be synchronized with the volatility of the chosen solid additive. However, this may be challenging since most organic photoactive materials display excellent efficiencies only with selected solvents. Accordingly, this work also demonstrates the potential of dual and complementary solvents selection, consisting of low boiling point (primary) and high boiling point (secondary). This strategy allows for the suppression of any potential trade-offs in efficiency. Meanwhile, the operational stability and precision of device performance are substantially enhanced. Additionally, solid additives have demonstrated that the singlet exciton dissociation rate does not limit the free charge generation yield. Finally, these findings are expected to reformulate OSC device fabrication strategies towards more practical devices.

36 MATERIALS SCIENCE↗

Multivariate analysis of gamma spectra to characterize used nuclear fuel

Abstract—The Multi-Isotope Process (MIP) Monitor provides an efficient approach to monitoring the process conditions in used nuclear fuel reprocessing facilities to support process verification and validation. The MIP Monitor applies multivariate analysis to gamma spectroscopy of reprocessing streams in order to detect small changes in the gamma spectrum, which may indicate changes in process conditions. This research extends the MIP Monitor by characterizing a used fuel sample after initial dissolution according to the type of reactor of origin (pressurized or boiling water reactor), initial enrichment, burn up, and cooling time. Simulated gamma spectra were used to develop and test three fuel characterization algorithms. The classification and estimation models employed are based on the partial least squares regression (PLS) algorithm. A PLS discriminate analysis model was developed which perfectly classified reactor type. Locally weighted PLS models were fitted on-the-fly to estimate continuous fuel characteristics. Burn up was predicted within 0.1% root mean squared percent error (RMSPE) and both cooling time and initial enrichment within approximately 2% RMSPE. This automated fuel characterization can be used to independently verify operator declarations of used fuel characteristics and inform the MIP Monitor anomaly detection routines at later stages of the fuel reprocessing stream to improve sensitivity to changes in operational parameters and material diversions.

Multivariate Analysis for Characterization of Used↗

Application of Jacobian-free Newton–Krylov method in implicitly solving two-fluid six-equation two-phase flow problems: Implementation, validation and benchmark

This work represents a first-of-its-kind successful application to employ advanced numerical methods in solving realistic two-phase flow problems with two-fluid six-equation two-phase flow model. These advanced numerical methods include high-resolution spatial discretization scheme with staggered grids (high-order) fully implicit time integration schemes, and Jacobian-free Newton–Krylov (JFNK) method as the nonlinear solver. The computer code developed in this work has been extensively validated with existing experimental flow boiling data in vertical pipes and rod bundles, which cover wide ranges of experimental conditions, such as pressure, inlet mass flux, wall heat flux and exit void fraction. Additional code-to-code benchmark with the RELAP5-3D code further verifies the correct code implementation. The combined methods employed in this work exhibit strong robustness in solving two-phase flow problems even when phase appearance (boiling) and realistic discrete flow regimes are considered. Transitional flow regimes used in existing system analysis codes, normally introduced to overcome numerical difficulty, were completely removed in this work. As a result, this in turn provides the possibility to utilize more sophisticated flow regime maps in the future to further improve simulation accuracy.

42 ENGINEERING↗

Heat up and failure of BWR upper internals during a severe accident

In boiling water reactors, the shroud dome, separators, and dryers above the core are made of approximately 100,000 kg of stainless steel. During a severe accident in which the coolant boils away and exothermic oxidation of zirconium occurs, gases (steam and hydrogen) are superheated in the core region and pass through the upper internals. In this scenario, the upper internals can also be heated by thermal radiation from the hot degrading core. Historically, models of the upper internals have been relatively simple in severe accident codes. The upper internals are typically modeled in MELCOR as two lumped volumes with simplified heat transfer characteristics and no structural integrity considerations, and with limited ability to oxidize, melt, and relocate. The potential for and the subsequent impact of the upper internals to heat up, oxidize, fail, and relocate during a severe accident was investigated. A higher fidelity representation of the shroud dome, steam separators, and steam driers was developed in MELCOR v1.8.6 by extending the core region upwards. The MELCOR modeling effort entailed adding 45 additional core cells and control volumes, 98 flow paths, and numerous control functions. The model accounts for the mechanical loading and structural integrity, oxidation, melting, flow area blockage, and relocation of the various components. Consistent with a previous study, the results indicate that the upper internals can reach high temperatures during a severe accident sufficient to lose their structural integrity and relocate. Finally, the additional 100 metric tons of stainless steel debris influences the subsequent in-vessel and ex-vessel accident progression.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental studies on metallic fuel relocation in a single-pin core structure of a sodium-cooled fast reactor

Experiments dropping molten uranium into test sections of single fuel pin geometry filled with sodium were conducted to investigate relocation behavior of metallic fuel in the core structures of sodium-cooled fast reactors during a hypothetical core disruptive accident. Metallic uranium was used as a fuel material and HT-9M was used as a fuel cladding material in the experiment in order to accurately mock-up the thermo-physical behavior of the relocation. The fuel cladding failed due to eutectic formation between the uranium and HT-9M for all experiments. The extent of the eutectic formation increased with increasing molten uranium temperature. Voids in the relocated fuel were observed for all experiments and were likely formed by sodium boiling in contact with the fuel. In one experiment, numerous fragments of the relocated fuel were found. In conclusion, it could be concluded that the injected metallic uranium fuel was fragmented and dispersed in the narrow coolant channel by sodium boiling

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal-hydraulic analyses of MIT reactor LEU transition cycles

The Massachusetts Institute of Technology Reactor (MITR) is a research reactor in Cambridge, Massachusetts designed primarily for experiments using neutron beam and in-core irradiation facilities. In the framework of non-proliferation policy, research and test reactors have started a program to convert highly enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. A new type of LEU fuel based on a high-density alloy of uranium and molybdenum (U-10Mo) is expected to allow conversion. A transition core plan, from 22 fresh LEU fuel elements to 24 in an equilibrium state, is proposed. This study evaluates the thermal-hydraulic safety margins of each transition cycle and state (i.e. beginning, middle and end). The STAT7 and RELAP5 codes are used in the study. STAT7 provides an integrated platform for statistical propagation of uncertainties in determining steady-state thermal-hydraulic margins. A RELAP5 model was created to verify MITR's STAT7 model, and provide transient analyses. The steady-state analyses determined that, during the entire transition from a fresh to equilibrium core, the Limiting System Safety Settings (LSSS) has a 10.2% margin to the minimum power at which onset of nucleate boiling (ONB) is precluded at 3σ confidence. The loss of flow (LOF) transient, which is the most limiting anticipated transient for the MITR, was simulated. With a conservative initial state of the reactor (operating at μ + 3σ power level and LSSS mass flow rate), the maximum calculated centerline fuel temperature was 115 degrees C, which is a significantly lower value than the blistering limit of 350 degrees C. The cladding wall temperatures did not exceed the ONB temperature throughout the transient. Therefore, during LOF, no nucleate boiling is expected, precluding any critical heat flux trigger. Overall, this study indicates that the proposed LEU transition core specifications have significant margin to thermal-hydraulic limits, during steady-state operation and the LOF transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fuel Properties of Oxymethylene Ethers with Terminating Groups from Methyl to Butyl

Oxymethylene ethers (OMEs) have been studied as possible additives or replacements for diesel fuels. Typically, studies have considered only methyl-terminated OMEs. Recent structure-property relationship models suggest that extended-alkyl OMEs may provide improvements to many of the properties of methyl-terminated OMEs that make them less suitable as diesel fuel blendstocks. In this work, we describe the synthesis and characterization of 16 different OMEs with methyl, ethyl, propyl, butyl, isopropyl, and isobutyl terminating alkyl groups with varying oxymethylene chain length. Indicated Cetane Number, Lower Heating Value, Flash Point, Density, Viscosity, Vapor Pressure, and Oxidative Stability are tested via ASTM standard methods. Additionally, Water Solubility, Boiling Point, seal material compatibility, and sooting propensity (via the Yield Sooting Index) are measured for these fuels. For diesel compatibility, all tested OMEs except smaller methyl and ethyl OMEs, and the branched isopropyl OME, meet cetane number requirements. Further, extending the alkyl end group increases the heating value, but all OMEs, due to their oxygen content, have heating values less than diesel; despite this, all OMEs show significant reductions in soot production per unit heating value. Only the heaviest OMEs meet diesel viscosity requirements, and most are higher density than diesel. OMEs with larger alkyl groups show the highest stability under accelerated auto-oxidation conditions. Increases in alkyl group length cause order of magnitude reduction in water solubility, from hundreds of g/L for methyl terminated OMEs to hundreds of mg/L for butyl terminated OMEs. Limited seal material testing indicates that PEEK polymers are unaffected by OMEs; while extended alkyl groups may improve compatibility with FKM (Viton), other common elastomers (NBR, silicone) remain incompatible with all tested OMEs. Overall, it is found that methyl-terminated OMEs exhibit the most potential for soot reduction, but OMEs with larger propyl and butyl terminating alkyl groups show improved compatibility with existing diesel systems.

09 BIOMASS FUELS↗