Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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,297 records · Page 72

Pairing Dynamics in Low-energy Nuclear Collisions

Superfluidity is a generic feature of various quantum systems at low temperatures and it is in particular important for the description of dynam- ics of low energy nuclear reactions. The Time-Dependent Density Func- tional Theory (TDDFT) is, to date, the only microscopic method which takes into account in a consistent way far-from-equilibrium dynamics of pairing field and single-particle degrees of freedom. Here, the local version of TDDFT, so-called TDSLDA, is particularly useful for the description of nuclear reactions and is well-suited for leadership class computers of hybrid (CPU+GPU) architecture. The preliminary results obtained for collisions involving both medium-mass and heavy nuclei at the energies around the Coulomb barrier are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Prokaryotic responses to a warm temperature anomaly in northeast subarctic Pacific waters

Recent studies on marine heat waves describe water temperature anomalies causing changes in food web structure, bloom dynamics, biodiversity loss, and increased plant and animal mortality. However, little information is available on how water temperature anomalies impact prokaryotes (bacteria and archaea) inhabiting ocean waters. This is a nontrivial omission given their integral roles in driving major biogeochemical fluxes that influence ocean productivity and the climate system. Here we present a time-resolved study on the impact of a large-scale warm water surface anomaly in the northeast subarctic Pacific Ocean, colloquially known as the Blob, on prokaryotic community compositions. Multivariate statistical analyses identified significant depth- and season-dependent trends that were accentuated during the Blob. Moreover, network and indicator analyses identified shifts in specific prokaryotic assemblages from typically particle-associated before the Blob to taxa considered free-living and chemoautotrophic during the Blob, with potential implications for primary production and organic carbon conversion and export.

54 ENVIRONMENTAL SCIENCES↗

Radiation damage effects in beryllium for next generation neutrino beam targetry (Final Technical Report)

Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.

36 MATERIALS SCIENCE↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite Sea Surface Salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (National Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF ( Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive) V4 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by the satellite at approximately 1 centimeter) and in situ measurements (typically measured by Argo floats at 3 meters). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from the GMAO S2S reanalyses that assimilates/withholds along-track (L2) SSS. In particular, we contrast forecasts of the overestimated 2014 El Nino, the big 2015 El Nino, and the minor 2016 La Nina. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrades the mixed layer depth leading to more accurate coupled air/sea interaction.

Hackert, E.↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO Seasonal Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite sea surface salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (North American Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF (Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive satellite) V4.1 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by satellite at approximately 1-centimeter depth) and in situ measurements (typically measured by Argo floats at 3-meters depth). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from GMAO S2S spring reanalyses that assimilate/withhold along-track (L2) SSS. In particular, we contrast forecasts of the big 2015 El Nino, the 2017 La Nina and the 2018 weak El Nino. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrade the mixed layer depth leading to more accurate coupled air/sea interaction. From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of gridded satellite sea surface salinity observations on dynamical ENSO forecasts for the big 2015 El Nino.

Hackert, E.↗

Direct imaging of plasma waves using ultrafast electron microscopy

A femtosecond plasma imaging modality based on a new development of ultrafast electron microscope is introduced. We investigated the laser-induced formation of high-temperature electron microplasmas and their subsequent non-equilibrium evolution. Based on a straightforward field imaging principle, we directly retrieve detailed information about the plasma dynamics, including plasma wave structures, particle densities, and temperatures. We discover that directly subjected to a strong magnetic field, the photo-generated microplasmas manifest in novel transient cyclotron echoes and form new wave states across a broad range of field strengths and different laser fluences. Intriguingly, the transient cyclotron waves morph into a higher frequency upper-hybrid wave mode with the dephasing of local cyclotron dynamics. The quantitative real-space characterizations of the non-equilibrium plasma systems demonstrate the feasibilities of a new microscope system in studying the plasma dynamics or transient electric fields with high spatiotemporal resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗

Calibrating Pressure Transducers at Cryogenic Temperatures

To calibrate transducers immersed in cryogenic environment, sealed reciprocating bellows generates sinusoidal dynamic pressure. Transducer manifold assembly and bellows immersed in cryogenic environment to minimize gradients in controllable cryogenic temperature.

Gibbens, B. V.↗

Corrections for Convective Heat Flux Gauges Subjected to a Surface Temperature Discontinuity

A two-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been carried out in an effort to determine the convective heat transfer corrections for circular heat flux gauges subjected to a surface temperature discontinuity. Solutions were obtained at a Reynolds number Of 1 x 10(exp 6) and a Mach number of 4. The CFD results are compared with the existing correlations for the correction factors. In general, the CFD corrections exceed those provided by the correlations. The discrepancy increases with increasing upstream surface temperature, thus indicating the role of property variations, which are not accounted for in the correlations. A quasi-two-dimensional analysis is also performed to treat the cylindrical geometry of the heat flux gauges by area-averaging the computed two-dimensional results from CFD.

Kandula, M.↗

Anatomy of Local Structural Disorder of Ni(II) Species in MgCl 2 –KCl Molten Salts

Understanding the speciation of metal ions dissolved in molten salts (MS) is critical for enabling a broad range of high-temperature energy applications, including MS nuclear reactors and concentrated solar power plants. However, due to the inherent dynamicity of metal species in the MS environment and the strong temperature dependencies of their multiple coexisting forms, they are difficult to resolve structurally. Herein, we show that combining in situ X-ray absorption spectroscopy (XAS) with ab initio molecular dynamics (AIMD) simulations is necessary to uncover and quantify the coexisting coordination states of Ni(II) in molten MgCl 2 –KCl mixtures and explain how the temperature and salt composition control their relative populations. Furthermore, from the interionic angle and distance distributions of nickel in different coordination states obtained from AIMD simulations, it is evident that for each coordination state, the width and skewness of their bonding distributions increase with increasing coordination number. In conclusion, the combination of XAS with first-principles modeling to resolve metastable metal species in MS is critical for understanding their behavior over a wide range of temperatures and chemical environments in nuclear and solar applications.

36 MATERIALS SCIENCE↗

Preliminary testing of a SNAP-19 TAGS RTG in support of the Pioneer F and G missions.

This paper presents the results of a preliminary performance evaluation of the SNAP-19 TAGS thermoelectric generator which was selected as the power source for the Pioneer Space Flight Mission. The performance of the SNAP-19 generator was evaluated at various input powers and fin-root temperatures. Results of this evaluation are analyzed and discussed. Tests were also performed to provide the Pioneer Flight Project Management with applicable data such as generator thermal dynamic behavior during launch and prior to deployment. These special tests measure generator performance as a function of fin-root temperature and the dynamic frequency response of the generator.

Stapfer, G.↗

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗

Flux flow and flux dynamics in high-T(sub c) superconductors

Because high temperature superconductors, including BYCO and BSCCO, are type 2 superconductors with relatively low H(sub c 1) values and high H(sub c 2) values, they will be in a critical state for many of their applications. In the critical state, with the applied field between H(sub c 1) and H(sub c 2), flux lines have penetrated the material and can form a flux lattice and can be pinned by structural defects, chemical inhomogeneities, and impurities. A detailed knowledge of how flux penetrates the material and its behavior under the influence of applied fields and current flow, and the effect of material processing on these properties, is required in order to apply, and to improve the properties of, these superconductors. When the applied field is changed rapidly, the time dependence of flux change can be divided into three regions, an initial region which occurs very rapidly, a second region in which the magnetization has a 1n(t) behavior, and a saturation region at very long times. A critical field is defined for depinning, H(sub c,p) as that field at which the hysteresis loop changes from irreversible to reversible. As a function of temperature it is found that H(sub c,p) is well described by a power law with an exponent between 1.5 and 2.5. The behavior of H(sub c,p) for various materials and its relationship to flux flow and flux dynamics are discussed.

Bennett, L. H.↗

Three-Dimensional Thermal Boundary Layer Corrections for Circular Heat Flux Gauges Mounted in a Flat Plate with a Surface Temperature Discontinuity

Three-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been performed in an effort to determine thermal boundary layer correction factors for circular convective heat flux gauges (such as Schmidt-Boelter and plug type)mounted flush in a flat plate subjected to a stepwise surface temperature discontinuity. Turbulent flow solutions with temperature-dependent properties are obtained for a free stream Reynolds number of 1E6, and freestream Mach numbers of 2 and 4. The effect of gauge diameter and the plate surface temperature have been investigated. The 3-D CFD results for the heat flux correction factors are compared to quasi-21) results deduced from constant property integral solutions and also 2-D CFD analysis with both constant and variable properties. The role of three-dimensionality and of property variations on the heat flux correction factors has been demonstrated.

Kandula, M.↗

Resolving Anodic Current and Temperature Distributions in a Polymer Electrolyte Membrane Water Electrolysis Cell Using a Pseudo-Two-Phase Computational Fluid Dynamics Model

Expanding upon our prior experimental work, we constructed a three-dimensional model of a polymer electrolyte membrane water electrolyzer using computational fluid dynamics. We applied the assumption of pseudo-two-phase flow, the flow of two phases with equal velocity. Experimental data were used to obtain parameters and to determine the conditions under which this model was valid. Anodic distributions of current density, temperature, liquid saturation, and relative humidity were obtained at various flow rates. The overall current density and temperature difference from inlet to outlet at the anode agreed strongly with experimental measurements under most circumstances. This verification allowed us to further examine the apparent gas coverage calculated from experimental and model temperature data. Results suggested a low liquid saturation and low relative humidity at the anode due to the consumption of liquid water and water vapor. However, we questioned the accuracy of the pseudo-two-phase assumption at low water feed rates. We concluded that the model was applicable to systems with liquid water feed rates greater than 0.6 ml min -1 cm -2 . Therefore, it is a fair screening method that can advise which operating conditions lead to excessive temperatures or drying at the anode, thereby promoting the longevity of the membrane and catalyst.

25 ENERGY STORAGE↗

Impact of Oxygen Non-Stoichiometry on Near-Ambient Temperature Ionic Mobility in Polaronic Mixed-Ionic-Electronic Conducting Thin Films

Enhanced ionic mobility in mixed ionic and electronic conducting solids contributes to improved performance of memristive memory, energy storage and conversion, and catalytic devices. Ionic mobility can be significantly depressed at reduced temperatures, for example, due to defect association and therefore needs to be monitored. Measurements of ionic transport in mixed conductors, however, proves to be difficult due to dominant electronic conductivity. This study examines the impact of different levels of quenched-in oxygen deficiency on the oxygen vacancy mobility near room temperature as measured by a novel dynamic current-voltage analysis. A Pr 0.1 Ce 0.9 O 2-δ film was grown by pulsed laser deposition and subsequently annealed, from 400-600° C, in various oxygen partial pressures to modify its oxygen vacancy concentration while minimizing microstructural growth and cation segregation. To monitor changes in film non-stoichiometry, we leverage the existence of an optical absorption center, related to the oxidation state of Pr ions in Pr0.1Ce0.9O2-δ. The oxygen vacancy migration enthalpy was found to exhibit a small increase from 0.73± 0.04 to 0.79 ± 0.02 eV with increasing oxygen deficiency, while the pre-factor was found to increase by a factor of 135. Here, a nearly 13-fold increase in ionic mobility at 60 °C for increases in oxygen non-stoichiometry from 0.032 ± 0.001 to 0.042± 0.001 was thereby detected. Raman spectroscopy was employed to rule out changes in strain as the primary cause for the significant change in mobility. Several factors potentially contributing to the large pre-factor changes are examined and discussed. Insights into how ionic defect concentration can markedly impact ionic mobility should help in elucidating the origins of variations seen in nanoionic devices.

36 MATERIALS SCIENCE↗