Convective carbon dioxide dissolution in a closed porous medium at high-pressure real-gas conditions
Not provided.
SEARCH · Search NASA
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.
Not provided.
This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.
Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.
Three models to compute the refractive index of gaseous mixtures at real gas conditions are presented with the purpose to improve the accuracy of state relationships between refractive index and thermodynamic properties. Models are compared with experimental data to determine one that is applicable to high pressure mixtures with non-ideal thermodynamic behavior near or above a fluids’ critical point. The optimal model is applied to analyze adiabatic thermal mixing of fuel and air at typical diesel engine conditions. Our results show that the ideal gas mixture law is appropriate in the vapor region, assuming it is known, for example, from experimental measurements. Finally, the model is applied for binary fuel-air mixing at supercritical conditions to demonstrate its full potential.
Not Available
Abstract not provided.
A novel continuously flowing liquid lithium limiter (FLiLi) which employs an in-vessel electro-magnetic pump to drive liquid Li flowing on the surface of limiter has been successfully designed and tested in the experimental advanced superconducting tokamak (EAST) device in 2014. In order to better control the surface temperature, an upgraded design and the real-time gas cooling were performed in 2016. Two kinds of cooling gas, helium (He) and argon (Ar) were tested on the FLiLi system prior to the experiment. It was found that the cooling rates of He and Ar at 2.5 MPa are 34.2 °C/min, 10.2 °C/min or 31.7 °C/kL and 28.1 °C/kL, respectively. The cooling performance of He is more effective than Ar. Experimental results show that the real-time He cooling could effectively contain the FLiLi surface temperature increase and prevent the strong passive evaporation of lithium during plasma discharges. Limiter heat flux up to 0.2 MW/m2 was removed during ohmic discharge with 2.5 MPa He gas. Finally, Because of the upgraded design and real-time gas cooling, the surface of FLiLi was not damaged by heat flux after entire plasma discharge.
Selenium is released into the environment via anthropogenic activities such as coal combustion. The U.S. Environmental Protection Agency’s Effluent Limitation Guideline (ELG) Rule limits the release of selenium species into U.S. waterways. To meet the ELG limit for selenium, coal-fired power plants must treat the wastewater produced by the flue gas desulfurization (FGD) systems. Selenium removal is of interest because small changes of the concentrations of selenium in the environment can cause detrimental ecological consequences. Biological treatment is used to remediate this wastewater, but little is known about the individual species in the microbial consortium that reduce toxic water-soluble selenium oxyanions to insoluble elemental selenium. To optimize current biological treatment capabilities of this water, the selenium-reducing capabilities must be leveraged. In order to do so, the selenium-reducing bacteria (SeRB) in biological treatment systems must be characterized. This study will be the first to characterize the microbial consortia from a real FGD WW biological treatment system and provide insights on how to better achieve selenium removal via biological treatment. This work will provide insights into the microorganisms predominantly responsible for the reduction of selenium oxyanions in FGD wastewater biological treatment systems, which can be used to better optimize more efficient biological treatment strategies of this industrial waste.
We developed a novel quantum platform with the complete capability to control individual atoms and their mutual interactions in real time. Ultracold atoms are employed since they have the lowest temperature and decoherence. This system is the ultimate dream of quantum physicists where the Hamiltonian of every single atom, including its kinetic energy, potential energy, and interactions, can be independently programmed in real time. We also initiated a new theoretical collaboration to explore its full potential.
Biological treatment is a recognized approach for removing toxic selenate and selenite oxyanions present in flue gas desulfurization (FGD) wastewater. However, the knowledge on the specific microbial species or communities responsible for reducing water-soluble selenium oxyanions to insoluble elemental selenium remains limited. In addition, the selenium oxyanion reduction genes and pathways have yet to be understood in these wastewaters. This study characterizes selenium oxyanion reducing bacteria (SeRB) native to FGD wastewater, and the resulting elemental selenium particles formed. By selecting for native SeRB microbes in a defined media, a novel resolution of these organisms has been achieved. This research identifies previously unrecognized selenium oxyanion reducing capabilities in Anaerosolibacter, alongside predominant SeRB from Mesobacillus and Tepidibacillus genera. This work encompasses both 16S and metagenomic techniques to recover novel metagenome assembled genomes (MAGs), distinct to this environment. The biogenic selenium produced by these organisms were predominantly of elemental selenium, in either amorphous or with a hexagonal structure. In addition, the elemental selenium particles formed where shown to increase in purity as the microbes were enriched. This study identifies the SeRB present in FGD wastewater and characterizes their selenium products, offering crucial insights to enhance the efficiency of biological treatment strategies and the potential of selenium recovery from this industrial waste.
The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.
Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.
Improving efficiency of hydrogen cooling in cryogenic conditions is important for the wider applications of hydrogen energy systems. The approach investigated in this study is based on a Ranque-Hilsch vortex tube (RHVT) that generates temperature separation in a working fluid. The simplicity of RHVT is also a valuable characteristic for cryogenic systems. In the present work, novel shapes of RHVT are computationally investigated with the goal to raise efficiency of the cooling process. Specifically, a smooth transition is arranged between a vortex chamber, where compressed gas is injected, and the main tube with two exit ports at the tube ends. Flow simulations have been carried out using STAR-CCM+ software with the real-gas Redlich-Kwong model for hydrogen at temperatures near 70 K. It is determined that a vortex tube with a smooth transition of moderate size manifests about 7% improvement of the cooling efficiency when compared vortex tubes that use traditional vortex chambers with stepped transitions and a no-chamber setup with direct gas injection.
Abstract Porous graphene and other atomically thin 2D materials are regarded as highly promising membrane materials for high‐performance gas separations due to their atomic thickness, large‐scale synthesizability, excellent mechanical strength, and chemical stability. When these atomically thin materials contain a high areal density of gas‐sieving nanoscale pores, they can exhibit both high gas permeances and high selectivities, which is beneficial for reducing the cost of gas‐separation processes. Here, recent modeling and experimental advances in nanoporous atomically thin membranes for gas separations is discussed. The major challenges involved, including controlling pore size distributions, scaling up the membrane area, and matching theory with experimental results, are also highlighted. Finally, important future directions are proposed for real gas‐separation applications of nanoporous atomically thin membranes.
Capabilities of the Plant Dynamics Code (PDC) have been extended to allow steady-state and transient simulation of graphite-cooled reactors with direct helium Brayton cycle. On the cycle side, the most significant code modification is the addition of helium properties, in the format required by the code’s equations. Since the code was already formulated to handle more complex real gas properties, adding helium as working fluid that behaves like ideal gas was fairly straightforward. A reactor module was added to PDC to simulate a reactor cooled by the working fluid of the Brayton cycle. Two options are supported: channel type, typical for graphite gas-cooled reactors, and pin type, typical for light-water and liquid metal-cooled reactors. The reactor module is an extension of the electrical heater model and simulates heat deposition in the fuel and transfer of this heat from the fuel to the coolant through the matrix and tube materials. The new reactor module becomes the third option in PDC for modeling heat addition to the cycle, besides previously modeled heat addition heat exchanger and electrical heater. In addition to those changes, other minor code modifications and improvements were introduced during the work of expanding PDC to modeling of gas-cooled reactors. These modifications are summarized in the last chapter of this report.
Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.
Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.
A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.