Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution species”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
NASA NTRS1974

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗
DOE OSTI2020

Beyond Strain: Controlling the Surface Chemistry of CsPbI 3 Nanocrystal Films for Improved Stability against Ambient Reactive Oxygen Species

Colloidal halide perovskite nanocrystals (NCs) have the possibility of easy scale-up due to their batch synthesis and have demonstrated excellent optoelectronic properties. In particular, perovskite NCs have remarkably high photoluminescence quantum yields in solution and as thin films and impressive open circuit voltages in photovoltaic devices. Despite these promising results, little work has been done to understand the stability of CsPbI 3 NCs for optoelectronic device applications. It has been previously shown that the ligands impart tensile surface strain, which stabilizes the black three-dimensional (3D) perovskite phase against phase degradation, making CsPbI 3 NCs some of the most structurally robust inorganic halide perovskites to date. However, understanding exactly how CsPbI 3 NCs degrade under ambient conditions is critical. Additionally, we demonstrate that the degradation mechanism of NCs is unique from, and 2 orders of magnitude slower than, their polycrystalline thin-film counterparts. Under specific conditions, CsPbI 3 NC films show a compositional instability instead of the phase instability seen in large grain CsPbI 3 . This is mediated through reactions with superoxide and other reactive oxygen species, which are initiated from surface defect states, O 2 and light. We then use this mechanistic insight to identify multiple strategies to prolong the lifetimes of CsPbI 3 NC films, by going beyond surface strain to mitigate key surface chemistries. We demonstrate that (1) minimizing the number of surface defects (2) using an alkylammonium bromide ligand surface treatment and (3) encapsulation with an oxygen scavenging layer all increase NC film lifetimes by inhibiting various steps in the photo-oxidation degradation reaction.

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

Long-Term Stability of Ferri-/Ferrocyanide as an Electroactive Component for Redox Flow Battery Applications: On the Origin of Apparent Capacity Fade

We assess the suitability of potassium ferri-/ferrocyanide as an electroactive species for long-term utilization in aqueous organic redox flow batteries. A series of electrochemical and chemical characterization experiments was performed to distinguish between structural decomposition and apparent capacity fade of ferri-/ferrocyanide solutions used in the capacity-limiting side of a flow battery. Our results indicate that, in contrast with previous reports, no structural decomposition of ferri-/ferrocyanide occurs at tested pH values as high as 14 in the dark or in diffuse indoor light. Instead, an apparent capacity fade takes place due to a chemical reduction of ferricyanide to ferrocyanide, via chemical oxygen evolution reaction. We find that this parasitic process can be further exacerbated by carbon electrodes, with apparent capacity fade rates at pH 14 increasing with an increased ratio of carbon electrode surface area to ferricyanide in solution. Based on these results, we report a set of operating conditions that enables the long-duration cycling of alkaline ferri-/ferrocyanide electrolytes and demonstrate how apparent capacity fade rates can be engineered by the initial system setup. If protected from direct exposure to light, the structural stability of ferri-/ferrocyanide anions allows for their practical deployment as electroactive species in long duration energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗
DOE OSTI2023

High temperature gamma radiation-induced chromium redox chemistry via in situ spectroscopic measurements

Chromium ions can make their way into the primary coolant of nuclear power reactors from the corrosion of stainless-steel reactor components, decreasing the material's corrosion resistance and resulting in increased transport of further corrosion products. Despite these potential effects, the radiation-induced redox speciation of chromium ions in aqueous solution is not well understood, especially at the elevated temperatures experienced by reactor coolants. Here, in the present work, we report new experimental results demonstrating that in aerated aqueous solution, the radiolytic oxidation of Cr(III) to Cr(VI) occurs at pH 4, while the reduction of Cr(VI) to Cr(III) occurs at pH 2. The oxidation of Cr(III) is primarily attributed to the reaction of the hydroxyl radical (˙OH) with the Cr(OH) 2+ species, while the reduction of Cr(VI) is attributed to reactions involving the hydrated electron (e aq - ) and hydrogen atom (H˙). Additionally, the steady-state equilibrium yield of Cr(VI) from the gamma irradiation of pH 4 Cr(III) solutions decreased with increasing temperature (over a range of 37–195 °C). This observation indicates that the activation energy of the Cr(VI) reduction reactions is higher than that for the Cr(III) oxidation reactions, such that it becomes relatively more favorable at higher temperatures. Overall, these data are important for the development of complementary multiscale models for the prediction of metal ion speciation in high temperature radiation environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2026

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
DOE OSTI2021

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H + (aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H + ·(H 2 O) n cluster ions isolated in the gas phase shed light on the chemical nature of H + (aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared–infrared (IR–IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H + (aq).

Chemistry↗
DOE OSTI2020

Intensified Flue Gas Desulfurization Water Treatment for Reuse, Solidification, and Discharge

The goal of this final project report is to summarize the work conducted on project DE-FE0031555. In accordance with the Statement of Project Objectives (SOPO), the University of Kentucky Center for Applied Energy Research (UK CAER) (Recipient) developed a hybrid electrochemical-membrane technology for treating wet flue gas desulfurization (WFGD) blowdown from coal-fired power generation plants that (1) removes government-regulated species and total suspended solids (TSS) via electrocoagulation with solid/liquid separation to avoid further treatments by a biological method specifically for selenium, and (2) reduces total dissolved solids (TDS) in a nanofiltration unit to promote water recycle and reduce freshwater intake. Additional technological configurations were pursued, including (1) nanofiltration combined with zeolite dewatering for a near-zero liquid discharge (ZLD) scenario, and (2) polishing discharge from zeolite dewatering using capacitive deionization (CDI) with 90+% salt rejection to produce water with less than 100ppm TDS for more beneficial use. Project results verified that UK CAER’s hybrid approach addressed the complexities of WFGD blowdown with simplified unit operation at relatively low capital investment by demonstrating: (1) an effective in-situ generated solid sorbent material, green-rust (GR), that specifically removes regulated dissolved species in WFGD blowdown, including selenium, arsenic, and nitrates; (2) a low-cost, high flux nanofiltration treatment that removes both monovalent and divalent species at >80% rejection to adequately provide water recycle and extend the life of FGD equipment by mitigating chloride aided pitting corrosion mechanisms; and (3) a solidification and stabilization process using pozzolanic agents that retain regulated species and satisfy the Resource Conservation and Recovery Act (RCRA) leachate requirements. UK CAER proved electrocoagulation with GR to be an intensified solution for effluent limit guidelines (ELG), and based on preliminary techno-economic analysis (TEA), in comparison to DOE’s Case 1 for Biological Wastewater Treatment that had a purchased equipment cost (PEC) of $\$$5.26 million and a Total plant cost (TPC) of $\$$30.0 million, the EC process had a PEC of $\$$2.67 million and the TPC of $\$$17.0 million. While the electricity consumption for Case 1 was 53 kW, it was 110 kW for the EC process due to dissolving iron for a three-stage reaction design. However, electricity consumption for the EC process could be significantly reduced by adopting a counter-current reactor design for at least a 6-stage chemical reaction operation.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1976

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

Penny, M. M.↗
NASA NTRS1981

Electromagnetic studies of redox systems for energy storage

Both chromium and iron couples were studied on various electrode surfaces in acidic perchlorate solution by using rotating ring-disk techniques. It was found that chloride which forms inner sphere coordination complexes with the redox species enhances the electrode kinetics dramatically. The effects of lead underpotential deposition and surface alloy formation on the kinetics of the chromium couple on gold were studied using both linear sweep voltammetry and potential step techniques. The lad underpotential deposition was found to slow down the kinetics of the reduction of the Cr species on gold surfaces although increase the hydrogen overvoltage. The effect on the chromium kinetics can be explained in terms of principally a double layer effect. The underpotential deposition lead species with its positive charge results in a decrease in the concentration of the Cr species at the electrode surface. Similar phenomena were also observed with bismuth underpotential deposition on gold for the iron couple.

Wu, C. D.↗
NASA NTRS1984

Stability of a laminar premixed supersonic free shear layer with chemical reactions

The stability of a two-dimensional compressible supersonic flow in the wake of a flat plate is discussed. The fluid is a multi-species mixture which is undergoing finite rate chemical reactions. The spatial stability of an infinitesimal disturbance in the fluid is considered. Numerical solutions of the eigenvalue stability equations for both reactive and nonreactive supersonic flows are presented and discussed. The chemical reactions have significant influence on the stability behavior. For instance, a neutral eigenvalue is observed near the freestream Mach number of 2.375 for the nonreactive case, but disappears when the reaction is turned on. For reactive flows, the eigenvalues are not very dependent on the free stream Mach number.

Menon, S.↗
NASA NTRS1986

Nonlinear upper hybrid drift waves for a longitudinal electric field perpendicular to a uniform magnetic field in the Vlasov-Maxwell approximation

Upper hybrid drift waves are found as a special solution to a Vlasov-Maxwell plasma which has a longitudinal electric field and a perpendicular uniform magnetic field. A single-species plasma with a constant-density mobile neutralizing background supports spatially varying disturbances that oscillate at the upper hybrid frequency. The general functional dependences of the electric field, the plasma number density, and the one-particle distribution function for the special case are found from more general Vlasov-Maxwell equations invariant under a Lie group point transformation. The one-particle distribution function for the plasma is a function of the Liouville invariant, which is the energy in the generalized Bernstein-Greene-Kruskal (BGK) reference frame, and the momentum in the drift direction.

Abraham-Shrauner, B.↗
NASA NTRS1990

Enhanced thermochemical nonequilibrium computations of flow around the aeroassist flight experiment vehicle

A three-dimensional explicit, finite-rate, shock-capturing numerical algorithm is used to calculate thermochemical nonequilibrium flowfields about the Aeroassist Flight Experiment vehicle at one of its flight trajectory points. The full Navier-Stokes equations and an eleven species chemical model with the latest reaction rates are incorporated into the code. Results are compared against experimental data and other numerical solutions. The effects of changes to the physical model on the stagnation line and base region flow and surface quantities on the forebody are investigated.

Palmer, Grant↗
NASA NTRS1991

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗
NASA NTRS1991

HANA - A three-dimensional Navier-Stokes code for chemically reacting flows

A recently developed computer program for analyzing chemically reacting flows is described. The program uses a diagonal implicit algorithm to solve the equations of viscous flow for a gas in thermochemical nonequilibrium. The diagonal scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The program uses multiple zones of grids patched together and includes a k-epsilon turbulence model. Solutions are presented for hypersonic flows of air, hypersonic flows of Martian atmospheric gases, combusting flows of hydrogen/air mixtures, and combusting flows of methane/air mixtures.

Imlay, Scott T.↗