Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 775 records · Page 43

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗

SNPP VIIRS Reflective Solar Bands On-Orbit Calibration Using the Moon

The Visible Infrared Imaging Radiometer Suite (VIIRS) on board the Suomi National Polar-Orbiting Partnership (SNPP) satellite has been on orbit for more than eight years since its launch on October 28, 2011. The VIIRS has 22 spectral bands, among which 14 are reflective solar bands(RSBs) covering a spectral range from 0.41 to 2.25 m. The VIIRS RSBs are primarily calibrated on-orbit by an onboard solar diffuser (SD) panel and a solar diffuser stability monitor (SDSM). Besides the SD and SDSM calibration, the RSBs are scheduled to view the Moon approximately monthly through the instrument’s space view (SV). The lunar observations have also been used to calibrate the RSBs on-orbit since early mission. Due to the non uniformity of the SD’s degradation, the calibration coefficients that are derived from the SD/SDSM calibration have long-term biases, especially at short wavelengths. In addition, the SDSM has no capability to monitor the SD degradation beyond 0.935 m, resulting in long-term bias in the shortwave infrared bands, about 0.72% for band M8 (1.238 m). These biases induce significant errors and long-term drifts in the VIIRS sensor data records (SDR) and consequently in the environmental data records (EDR). Unlike the SD, the Moon is a known stable target and any temporal drifts observed while viewing the Moon can be attributed to the sensor’s degradation. Thus, the VIIRS lunar calibration is used to track the RSB on-orbit changes, especially to provide an accurate long-term baseline. Due to the non-uniformity of the lunar surface, the lunar irradiance, instead of the lunar radiance, is used to calibrate the RSBs. The lunar irradiance strongly depends on lunar view geometry and it is still a challenge to accurately characterize the geometric effects associated with the lunar measurements and any residual errors can induce seasonal oscillations in the derived calibration coefficients. The errors of the geometric dependence correction induce seasonal oscillations in the derived RSB lunar calibration coefficients. In this paper, the algorithms for the view geometric effect correction are significantly improved, resulting in a significant reduction in the seasonal oscillations observed in the calibration coefficient time-series. The lunar and SD/SDSM calibration results are properly incorporated to generate a set of hybrid calibration coefficients and implementation of these coefficients is shown to significantly improve the long-term stability of the VIIRS SDR. This is of fundamental importance in making accurate Earth observations from which reliable and high quality science products are generated. The consequent improvements in SNPP VIIRS RSB SDRs and EDRs are shown and discussed. The lunar calibration methodology can be directly applied to follow-on VIIRS instruments.

Junqiang Sun↗

Radiation-Induced Effects on the Extraction Properties of Hexa-n-octylnitrilo-triacetamide (HONTA) Complexes of Americium and Europium

The candidate An(III)/Ln(III) separation ligand hexa-n-octylnitrilo-triacetamide (HONTA) was irradiated under envisioned SELECT (Solvent Extraction from Liquid waste using Extractants of CHON-type for Transmutation) process conditions (n-dodecane/0.1 M HNO 3 ) using a solvent test loop in conjunction with cobalt-60 gamma irradiation. The extent of HONTA radiolysis and complimentary degradation product formation was quantified by HPLC-ESI-MS/MS. Further, the impact of HONTA radiolysis on process performance was evaluated by measuring the change in 243 Am and 154 Eu distribution ratios as a function of absorbed gamma dose. HONTA was found to decay exponentially with increasing dose, affording a dose coefficient of d = (4.48 ± 0.19) × 10 –3 kGy –1 . Multiple degradation products were detected by HPLC-ESI-MS/MS with dioctylamine being the dominant quantifiable species. Both 243 Am and 154 Eu distribution ratios exhibited an induction period of ~70 kGy for extraction (0.1 M HNO 3 ) and back-extraction (4.0 M HNO 3 ) conditions, after which both values decreased with absorbed dose. The decrease in distribution ratios was attributed to a combination of the destruction of HONTA and ingrowth of dioctylamine, which is capable of interfering in metal ion complexation. Furthermore, the loss of HONTA with absorbed gamma dose was predominantly attributed to its reaction with the n-dodecane radical cation (R˙ + ). These R˙ + reaction kinetics were measured for HONTA and its 241 Am and 154 Eu complexes using picosecond pulsed electron radiolysis techniques. All three second-order rate coefficients (k) were essentially diffusion limited in n-dodecane indicating a significant reaction pathway: k(HONTA + R˙ + ) = (7.6 ± 0.8) × 10 9 M –1 s –1 , k(Am(HONTA) 2 + R˙ + ) = (7.1 ± 0.7) × 10 10 M –1 s –1 , and k(Eu(HONTA) 2 + R˙ + ) = (9.5 ± 0.5) × 10 10 M –1 s –1 . HONTA-metal ion complexation afforded an order-of-magnitude increase in rate coefficient. Additional nanosecond time-resolved measurements showed that both direct and indirect HONTA radiolysis yielded the short-lived (< 100 ns) HONTA radical cation and a second long-lived (µs) species identified as the HONTA triplet excited state. The latter was confirmed by a series of oxygen quenching picosecond pulsed electron measurements, affording a quenching rate coefficient of k( 3 [HONTA]* + O 2 ) = 2.2 × 10 8 M –1 s –1 . Overall, both the HONTA radical cation and triplet excited state are important precursors to the suite of measured HONTA degradation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Water Flow Performance of a Superscale Model of the Fastrac Liquid Oxygen Pump

As part of the National Aeronautics and Space Administration's ongoing effort to lower the cost of access to space, the Marshall Space Flight Center has developed a rocket engine with 60,000 pounds of thrust for use on the Reusable Launch Vehicle technology demonstrator slated for launch in 2000. This gas generator cycle engine, known as the Fastrac engine, uses liquid oxygen and RP-1 for propellants and includes single stage liquid oxygen and RP-1 pumps and a single stage supersonic turbine on a common shaft. The turbopump design effort included the first use and application of new suction capability prediction codes and three-dimensional blade generation codes in an attempt to reduce the turbomachinery design and certification costs typically associated with rocket engine development. To verify the pump's predicted cavitation performance, a water flow test of a superscale model of the Fastrac liquid oxygen pump was conducted to experimentally evaluate the liquid oxygen pump's performance at and around the design point. The water flow test article replicated the flow path of the Fastrac liquid oxygen pump in a 1.582x scale model, including scaled seal clearances for correct leakage flow at a model operating speed of 5000 revolutions per minute. Flow entered the 3-blade axial-flow inducer, transitioned to a shrouded, 6- blade radial impeller, and discharged into a vaneless radial diffuser and collection volute. The test article included approximately 50 total and static pressure measurement locations as well as flush-mounted, high frequency pressure transducers for complete mapping of the pressure environment. The primary objectives of the water flow test were to measure the steady-state and dynamic pressure environment of the liquid oxygen pump versus flow coefficient, suction specific speed, and back face leakage flow rate. Initial results showed acceptable correlation between the predicted and experimentally measured pump head rise at low suction specific speeds. Likewise, only small circumferential variations in steady-state were observed from 80% to 120% of the design flow coefficient, matching the computational predictions and confirming that the integrated design approach has minimized any exit volute-induced distortions. The test article exhibited suction performance trends typically observed in inducer designs with virtually constant head rise with decreasing inlet pressure until complete pump head breakdown. Unfortunately, the net positive suction head at 3% head fall-off occurred far below that predicted at all tested flow coefficients, resulting in a negative net positive suction head margin at the design point in water. Additional testing to map the unsteady pressure environment was conducted and cavitation-induced flow disturbances at the inducer inlet were observed. Two distinct disturbances were identified, one rotating and one stationary relative to the fixed frame of reference, while the transition from one regime to the next produced significant effects on the steady state pump performance. The impact of the unsteady phenomena and the corresponding energy losses on the unexpectedly poor pump performance is also discussed.

Skelley, Stephen↗

Effects of Ambient Inert Gas Properties on Cool-Flame Combustion of Normal-Alkane Droplets in Microgravity

Cool-flame-supported quasi-steady droplet combustion was first observed during droplet combustion experiments conducted onboard the International Space Station (ISS) in the Combustion Integrated Rack (CIR) with n-alkane fuels. These cool flames, not visible to the naked eye, result from the existence of the negative-temperature-coefficient (NTC) region during the combustion n-alkane fuels when the chemical-kinetic pathways experience a transition from the low-temperature chemistry to the normal, high-temperature-controlled chemistry. In this study, we present results for the cool-flame-supported droplet burning rates and extinction diameters for n-alkane fuels (n-heptane, n-octane, n-decane, and n-dodecane) in inert-diluent-substituted ambient environments. In these ISS experiments, some of the nitrogen in the ambient gas mixture is replaced by an inert gas, namely, carbon dioxide, helium, or xenon. Different diluents were observed to exert remarkably different influences on the cool-flame combustion of fuel droplets. The results show that the cool-flame extinction diameter is drastically increased with helium dilution compared to N2, CO2 or Xe, and Xe dilution leads to longer burning with a smaller extinction diameter. Unlike hot flames, the cool-flame burning rates are controlled by the cool-flame temperature established in the NTC region by the chemical kinetics, rather than by equilibrium thermodynamics, and the thermo-physical properties of the ambient gases, which control the heat feed back to the liquid droplet. The cool-flame extinction is triggered by the instability that occurs at the low-temperature end of the NTC region, below which cool diffusion flames become statically unstable. A quasi-steady theoretical model, presented earlier, which involves a six-step reduced-chemistry mechanism, valid around the lower end of the NTC region, is shown to explain the observed experimental trends qualitatively.

microgravity↗

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗

Polarimetric calibration of the Multi-Angle Imager for Aerosols (MAIA)

NASA’s Multi-Angle Imager for Aerosols (MAIA) mission, under development at the Jet Propulsion Laboratory, is designed to study the adverse health effects of different types of particulate air pollution. Planned for launch in 2023 for a 3-year mission, the MAIA satellite instrument will focus on a selected set of metropolitan target areas, where air quality monitors and health data are available. Aerosol concentration and speciation are inferred from multi-angle measurements of backscattered sunlight in 14 spectral bands from 350-2200 nm, with bands near 442, 645 and 1040 nm measuring the degree (DoLP) and angle of linear polarization (AoLP) in addition to radiance. The pushbroom camera has a ~240-km cross-track field of view with a nadir resolution of ~200 m, and is mounted onto a biaxial gimbal to provide along-track view angles within ±60°, to extend the field of regard to ±48°, and to view the instrument’s onboard calibrator (OBC) and dark target. The OBC consists of a sunlit transmissive diffuser, followed by 12 polarizers at different orientations. MAIA’s polarimetry is implemented using miniature wire grid polarizers on the focal plane array, and dual photoelastic modulators (PEMs) and achromatic quarter-wave plates to rapidly rotate the polarization. The resulting ~26-Hz intensity modulation encodes the linearly polarized and total radiance in each pixel, leaving the DoLP and AoLP insensitive to gain calibration. We report on the polarimetric calibration of the MAIA camera using a vacuum-compatible polarization state generator, consisting of a 1600W Xenon lamp, 12-inch integrating sphere, and rotating high-extinction polarizer. Mueller-matrix-based calibration coefficients for each detector pixel are derived from measurements at multiple polarizer angles, and are used to correct the measurements for instrumental polarization aberrations. Prior to flight, the calibrated MAIA camera is panned across the OBC to characterize its output, using uniform illumination with an irradiance similar to the Sun.

Werne, Thomas↗

The Influence of Temperature Inversion on the Vertical Distribution of Aerosols

Temperature inversion plays an important role in the accumulation and diffusion of aerosols. In this study, the relationship between temperature inversion and the vertical distribution of aerosols is investigated based on Raman lidar observations taken from January 2010 to September 2015 at the Atmospheric Radiation Measurement site in the Southern Great Plains, USA. First, the diurnal and seasonal variations of the surface-based inversion (SBI) and elevated temperature inversion (EI) are investigated. The results indicate that the occurrence frequency of SBI and EI have different seasonal trends. SBI has the highest frequency in summer, while EI has the highest frequency in winter. The diurnal variation of SBI is obvious, with a higher frequency in nighttime and a lower frequency in daytime. The inversion intensity (ΔT) and inversion depth (ΔZ) of SBI and EI have consistent diurnal and seasonal trends. The effects of SBI and EI on the vertical distribution of aerosols are then analyzed. The mean aerosol optical depth (AOD) below the SBI height shows a clear seasonal variation, which is consistent with the seasonal trends of ΔT and ΔZ. This phenomenon also occurs on the AOD below EI top height. The sensitivity analysis shows that the mean AOD below SBI height or EI top height increases with an increase of the ΔT and ΔZ of SBI (EI). It indicates that ΔT and ΔZ are the key factors affecting the vertical distribution of aerosols. In addition, the variation of AOD below and above EI top height is opposite to that of AOD below and above EI bottom height under different ΔT and ΔZ conditions. The correlation coefficients between ΔT (ΔZ) of EI with AOD in EI were 0.62 (0.65). These results indicate that the space between EI bottom height and EI top height can store aerosols. The larger the ΔZ of EI, the more aerosols are stored. These findings contribute to our understanding of the effect of temperature inversion on the vertical distribution of aerosols.

54 ENVIRONMENTAL SCIENCES↗

Demonstration of a multi-channel fluidized bed particle–supercritical carbon dioxide heat exchanger for concentrating solar applications

High-temperature thermal energy storage in oxide particles at temperatures above 600°C can couple concentrated solar energy with high-efficiency thermal power cycles to provide dispatchable solar-driven electricity. Challenges remain in developing cost-effective primary heat exchangers, which require expensive alloys, to extract the high-temperature thermal energy from the particles to power cycle fluids, such as supercritical CO 2 (sCO 2 ) in recuperated Brayton cycles. To explore one pathway for cost-effective, high-temperature particle heat exchangers, the current study demonstrates a shell-and-plate, particle–sCO 2 heat exchanger with narrow- channel fluidized beds coupled with micro-channel sCO 2 flows in the heat exchanger walls. This study evaluates the feasibility of multiple parallel, narrow-channel fluidized beds in shell-and-plate particle–sCO 2 HXs, to achieve high bed-wall heat fluxes at elevated temperatures. A reduced-order model simulates the narrow- channel, fluidized-bed particle–sCO 2 heat exchanger to design the fluidized bed geometry, in terms of depth, height, and number of channels,for a nominal 40-kWth heat exchanger at particle and sCO 2 inlet temperatures up to 600 °C and 400 °C respectively. The resulting shell-and-plate heat exchanger design operates with bubbling fluidization of the downward-flowing oxide particles to enhance bed-wall heat transfer. The heat exchanger core is fabricated with etched sCO 2 micro-channels in thin wall plates that are diffusion bonded to spacer frames to form the shell-and-plate structure with 12 parallel, fluidized bed channels, 10.4 mm deep. The heat exchanger is tested at the National Solar Thermal Test Facility at Sandia National Laboratories with CARBOBEAD HSP particles at design particle flow rates of 0.20 kg s –1 and inlet temperatures up to 525 °C. Results show that fluidization across multiple parallel channel beds can maintain uniform particle inventory with a common freeboard zone above the heat exchanger core. Bubbling fluidization improves particle–wall heat transfer coefficients but also increases axial dispersion of particle thermal energy, which lowers the log- mean temperature difference such that total heat transfer remains relatively constant to within ±10% over a broad range of fluidization gas velocities. The axial dispersion required particle and sCO 2 flow rates to be increased by 25% over model-designed conditions to achieve the targeted 40 kWth, which indicates the importance of incorporating axial dispersion into heat exchanger design models and of deploying bed structures to suppress it. Furthermore, this study demonstrates the feasibility and preferred fluidizing gas conditions for particle heat exchangers for releasing high-temperature thermal energy storage systems.

14 SOLAR ENERGY↗

A Global End-Member Approach to Derive aCDOM(440) from Near-Surface Optical Measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end members, i.e., a short and long wavelength pair. The global perspective is established by sampling "extremely" clear water plus a generalized extent in turbidity and optical properties that each span three decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the east and west coasts of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses, i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent, plus 15 subcategories of water masses that were not evolving conservatively. The dependence on subcategories was confirmed with an objective methodology based on cluster analysis techniques. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four different algorithms based on different validation quality end members were validated with accuracies of 1.–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing non-conservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1mm and a boundary constraint to mitigate wave focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Stanford B Hooker↗

Time-dependent modulation of galactic cosmic rays by merged interaction regions

Models that solve the one-dimensional, solar modulation equation have reproduced the 11-year galactic cosmic ray using functional representations of global merged interaction regions (MIRs). This study extends those results to the solution of the modulation equation with explicit time dependence. The magnetometers on Voyagers 1 and 2 provide local magnetic field intensities at regular intervals, from which one calculates the ratio of the field intensity to the average local field. These ratios in turn are inverted to form diffusion coefficients. Strung together in radius and time, these coefficents then fall and rise with the strength of the interplanetary magnetic field, becoming representations of MIRs. These diffusion coefficients, calculated locally, propagate unchanged from approx. 10 AU to the outer boundary (120 AU). Inside 10 AU, all parameters, including the diffusion coefficient are assumed constant in time and space. The model reproduces the time-intensity profiles of Voyager 2 and Pioneer 10. Radial gradient data from 1982-1990 between Pioneer 10 and Voyager 2 are about the same magnitude as those calculated in the model. It is also shows agreement in rough magnitude with the radial gradient between Pioneer 10 and 1 AU. When coupled with enhanced, time-dependent solar wind speed at the probe's high latitude, as measured by independent observers, the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source for the observed negative latitudinal gradients. The model exhibits the 11-year cyclical cosmic ray intensity behavior at all radii, including 1 AU, not just at the location of the spacecraft where the magnetic fields are measured. In addition, the model's point of cosmic ray maximum correctly travels at the solar wind speed, illustrating the well-known propagation of modulation. Finally, at least in the inner heliosphere this model accounts for the delay experienced by lower-rigidity protons in reaching their time-intensity peak. The actual delays in this model, however, are somewhat smaller than the data. In the outer heliosphere the models sees no delays, and the data are ambiguous as to their existence. It appears that strong magnetic field compression regions (merged interaction regions) that are 3-4 times the average field strength can, at least in a helioequatorial band, disrupt effects, such as drifts, that could dominate in quieter magnetic fields. The question remains: Is the heliosphere ever quiet enough to allow such effects to be unambiguously measured, at least in the midlatitudes?

Perko, J. S.↗

Understanding the ignition process and flame structure of conventional and oxygenated fuels under engine relevant conditions – An optical study

Renewably generated synthetic fuels such as poly-oxymethylene ethers (OME) have a significant potential to effectively break the soot-NOX trade-off in compression ignition engines by using exhaust gas recirculation (EGR) to maintain low nitrogen oxide (NOX) emissions while maintaining good efficiency and simultaneously contributing to circular carbon economy. However, owing to the fundamental differences in properties of OME when compared to fossil-based diesel fuels, it is critical to fully understand its ignition and combustion phenomenology to take advantage of this fuel to its utmost potential. In this context, this work outlines the results of a systematic experimental study performed in a heavy-duty, single-cylinder, optical engine probing the spatial and temporal progression of fuel decomposition and ignition behavior of OME when compared to n-dodecane, a diesel-fuel surrogate. Thermodynamic analysis and optical diagnostics techniques including simultaneous HCHO-PLIF and OH-PLIF complemented by high-speed OH* chemiluminescence were employed along with parametric sweeps of intake temperature and EGR dilution rates. OME does not exhibit any observable low temperature heat release irrespective of the ambient oxygen concentration. Differences in the observed diffusive flame structure such as longer flame lift-off length, less pronounced combustion recession, faster premixed burn at ignition (“volumetric” ignition), non-sooting behavior suggest that the inherent presence of fuel-bound oxygen in OME can skew the air-fuel ratio (AFR) distribution within the jet thereby reducing the reliance of combustion on mixing and air entrainment. This leads to rapid late-cycle oxidation leading to shorter combustion duration and favorable combustion phasing. Results also suggest that OME exhibits relatively weak negative temperature coefficient (NTC) behavior, however, the OME fuel-decomposition kinetic-pathways produce significant concentration of HCHO, which might be erroneously interpreted as a product of cool-flames.

Air-fuel ratio↗

Understanding radiation effects in friction stir welded MA956 using ion irradiation and a rate theory model

We report an outstanding challenge in the manufacturing and joining of oxide dispersion strengthened steels is retaining the nanofeatures in the alloy throughout the fabrication and welding process. MA956 was friction stir welded with two different sets of welding parameters, resulting in a medium and high heat input. After welding, 5 MeV Fe ++ ion irradiations were performed at doses ranging from 50 to 200 dpa in the temperature range of 400 to 500°C. Post-irradiation characterization was performed with scanning transmission electron microscopy and energy-dispersive x-ray spectroscopy to investigate the Y-Al-O dispersoids, voids, and dislocations. After welding, the dispersoid microstructure coarsened, resulting in fewer and larger dispersoids regardless of heat input. After irradiation, the dispersoid behavior in the welded material was sensitive to temperature, exhibiting growth behavior attributed to Ostwald coarsening at 500 °C but a mixture of nucleation and more muted growth at 400 and 450 °C, attributed to competing mechanisms of radiation-enhanced diffusion and Ostwald coarsening. Void swelling correlated to heat input; being more prevalent in the welded conditions occurring at lower doses and in higher values relative to the base material. The low values of swelling despite microstructure coarsening caused by welding demonstrate the excellent swelling resistance of MA956, even after welding with the highest swelling values of 0.5% noted in the stir zone high heat input condition at 450 °C, 200 dpa. The dislocation behavior was inconsistent: the strongest trend was that network density was higher for welded versus base material, and an increase in loop diameter with temperature was observed. A rate theory model based on the observed microstructure suggests at high temperature interstitial loss to sinks was more likely to be dominant compared to mutual annihilation via point defect recombination, because of an increase of the radiation diffusion coefficient with temperature regardless of initial welded microstructure.

36 MATERIALS SCIENCE↗

A Modeling Study of Oceanic Response to Daily and Monthly Surface Forcing

The goal of this study is to investigate the effect of high-frequency surface forcing (wind stresses and heat fluxes) on upper-ocean response. We use the reduced-gravity quasi-isopycnal ocean model by Schopf and Loughe (1995) for this study. Two experiments are performed: one with daily and the other with monthly surface forcing. The two experiments are referred to as DD and MM, respectively. The daily surface wind stress is produced from the SSM/I wind data (Atlas et al. 1991) using the drag coefficient of Large and Pond (1982). The surface latent and sensible heat fluxes are estimated using the atmospheric mixed layer model by Seager et al. (1995) with the time-varying air temperature and specific humidity from the NCEP-NCAR reanalysis (Kalnay et al. 1996). The radiation is based on climatological shortwave radiation from the Earth Radiation Budget Experiment (ERBE) [Harrison et al. 1993] and the daily GEWEX SRB data. The ocean model domain is restricted to the Pacific Ocean with realistic land boundaries. At the southern boundary the model temperature and salinity are relaxed to the Levitus (1994) climatology. The time-mean SST distribution from MM is close to the observed SST climatology while the mean SST field from DD is about 1.5 C cooler. To identify the responsible processes, we examined the mean heat budgets and the heat balance during the first year (when the difference developed) in the two experiments. The analysis reveals that this is contributed by two factors. One is the difference in latent heat flux. The other is the difference in mixing processes. To further evaluate the responsible processes, we repeated the DD experiment by reducing the based vertical diffusion from 1e-4 to 0.5e-5. The resultant SST field becomes quite closer to the observed SST field. SST variability from the two experiments is generally similar, but the equatorial SST differences between the two experiments show interannual variations. We are investigating the possible mechanisms responsible for the different responses.

Sui, Chung-Hsiung↗

The Effects of Gravity on the Crystallization Behavior of Heavy Metal Fluoride Glasses

Heavy metal fluoride glasses are used in such applications as fiber lasers and laser amplifiers. ZrF4-BaF2-LaF3-AlF3-NaF (ZBLAN) is one of the more commonly used heavy metal fluoride glasses. ZBLAN is an infrared transmitter and has a theoretical attenuation coefficient of 0.002 db/km. However, due to impurities and small crystallites this attenuation coefficient has not been achieved to date. ZBLAN is a fragile glass which can lead to rapid crystallization, if the glass is not cooled rapidly to below the glass transition temperature or if the glass is reheated near the crystallization temperature for any period of time. Studies carried on at Marshall Space Flight Center and the University of Alabama in Huntsville since 1993 have shown that heating ZBLAN glass at the crystallization temperature in reduced gravity results in a suppression of crystallization when compared to ZBLAN processed in unit gravity. These studies utilized NASA's KC-135 aircraft and the Conquest sounding rocket. In the first series of experiments, short lengths of ZBLAN fiber were heated to the crystallization temperature in reduced gravity on board the KC- 135 and the Conquest sounding rocket and compared with fibers heated in unit gravity. The fibers processed in reduced gravity showed no evidence of crystallization when studied with x-ray diffraction and scanning electron microscopy. However, the fibers processed in unit gravity were completely crystallized. Subsequent experiments included heating small pieces of ZBLAN glass at the crystallization temperature while viewing with a video camera to follow the crystallization phenomenon. In this experiment crystallization was observed in reduced gravity, however, it was suppressed when compared to heating in unit gravity. In the most recent experiment on board the KC-135, rapid thermal analysis of ZBLAN was performed. A mechanism to explain the observations has been proposed. This mechanism is based on shear thinning whereby, the glass undergoes a reduction in viscosity in the 10(exp 5) - 10(exp 7) poise range, allowing more rapid diffusion and thus crystallization. It is proposed that this mechanism is suppressed in reduced gravity. An experiment is presently being conducted to test this theory. With increased knowledge of ZBLAN behavior in reduced gravity, three low earth orbit tiber drawing facilities have been designed. One would be suitable for use on the International Space Station, another while aboard the Space Shuttle and the third system is a fully automated facility which would operate independently of the ISS or Shuttle in a free float mode. The primary benefits of free floating a facility in LEO includes a higher quality of microgravity and reduced safety concerns since it is not in a manned environment.

Tucker, Dennis S.↗

Ablation of ATJ Graphite at High Temperatures

Results of an extensive experimental program on the ablation of ATJ graphite in air at surface pressures of 0.3 to 4.4 atm and surface temperatures of 2570 to 4030 K are presented. The measured mass-loss rates are correlated with surface temperature, pressure, and effective nose radius. The results are compared with several equilibrium thermochemical ablation theories which differ basically in the chemical species that are considered and in the thermodynamic properties assigned to these species. All the theories predict about the same mass-loss rate in the diffusion-controlled oxidation regime and are in good agreement with the experimental results. At higher temperatures, however, the experimental and theoretical results do not agree. At temperatures above 3700 K, the experimental mass-loss rate becomes independent of pressure and an exponential function of temperature. A comparison of the high-temperature results with the theories indicates that all the theories, except the one based on the thermodynamic properties from the JANAF tables, are invalid because they overpredict the mass-loss rate. On the other hand, the theory based on the JANAF properties underpredicts the measured mass-loss rate by a factor of three at 4000 K. At least a portion of this difference between experiment and theory is attributed to particulate mass loss which is visually observed.

Thermochemical ablation↗

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗