Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gas and solid dynamics”

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 55 records · Page 3

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

36 MATERIALS SCIENCE↗

Spectral Similarity Masks Structural Diversity at Hydrophobic Water Interfaces

The air-water and graphene-water interfaces represent quintessential examples of the liquid-gas and liquid-solid boundaries, respectively. While the sum-frequency generation (SFG) spectra of these interfaces show similarities, a consensus on their signals and interpretations has yet to be reached. Leveraging deep learning, we computed first-principles SFG spectra for both systems, addressing experimental discrepancies. Here, our findings reveal that similarities in SFG signals do not translate into comparable interfacial microscopic properties. Instead, graphene-water and air-water interfaces exhibit fundamental differences in SFG-active thicknesses, hydrogen-bonding networks, and surface dynamics. These distinctions underscore roughness suppression and electronic interactions present at the solid-liquid interface but absent at the gas-liquid interface.

Wang, Yong [Princeton Univ., NJ (United States)] (↗

Rapid Load Transition for Integrated Solid Oxide Fuel Cell - Gas Turbine (SOFC-GT) Energy Systems: A Demonstration of the Potential for Grid Response

Rapid load transition is an essential requirement for integrated energy systems to maintain grid resilience as more renewable resources are added to the grid. Integrated solid oxide fuel cell - gas turbine (SOFC-GT) systems can provide high efficiency and low carbon emissions over a broad range of turndown. These hybrids also have the potential to enable rapid grid response. The challenge has been to demonstrate effective control strategies to manage load transitions. In the present study, a load transition of ~50% was achieved in 10 s using a novel but simple strategy. Power demand on the SOFC and the GT were ramped down concurrently. During this transition, the SOFC anode fuel was manipulated to maintain SOFC fuel utilization while the cathode inlet air flow and temperature were also manipulated to thermally protect the SOFC. This study was conducted using the Hybrid Performance (Hyper) facility at the National Energy Technology Laboratory (NETL) in a co-simulation environment with the Idaho National Laboratory (INL)'s grid-simulation. The load ramping strategy was tested using a hardware-based cyber-physical simulation methodology. The results demonstrate a high-fidelity representation of SOFC-GT hybrid dynamics and validation of the control strategy. Thermal and electrochemical transients indicated that the SOFC was well protected during rapid load turndown without violating operability constraints. This demonstration revealed the non-linear nature of tightly coupled SOFC-GT system components, especially the non-linear response of SOFC cathode air flow and inlet temperature controls. These results highlight the needs and challenges in developing adaptive automatic controls for autonomous rapid load transitions. This work demonstrates that SOFC-GT hybrids are a viable option to provide the fast-ramping characteristics essential to accommodate high levels of variable renewable power while maintaining grid resilience, reliability, and environmental performance. The results also demonstrate the utility of co-simulation in advancing the tightly-coupled integrated energy systems needed to meet goals for zero-carbon power generation.

DIRECT ENERGY CONVERSION,POWER TRANSMISSION AND DI↗

In-situ monitoring of atomically dispersed Pt sites supported on OMS-2 during CO 2 activation

Atomically dispersed catalysts have drawn great interest lately, as they showcase a high density of active sites, selectivity, and high turnover frequencies in oxidation chemistry due to labile oxygen activation. In contrast, the applications of these catalysts have lagged in reduction reactions due to the ambiguity caused by the sintering and restructuring of active sites. To bridge this gap, the evolution of Pt 4+ isomorphically substituted into an octahedral molecular sieve structure (OMS-2) under reductive conditions was correlatively characterized using multiple in-situ analytical techniques such as ambient pressure X-ray photoelectron spectroscopy, environmental transmission electron microscopy, and solid-state nuclear magnetic resonance. The surface dynamics of the Pt single atoms were revealed during the Reverse Water Gas Shift (RWGS) reaction, where the active sites were identified as two-coordinated platinum single atoms. Under reaction, we show nonbinding atoms adjacent to the single atoms restructured the motif of the single atoms to Pt 2+ via ion mobility of potassium, increasing the activation energy by 25.6 kJ/mol. Here, this work also highlights the potential for increased stability of the single atom sites via isomorphic substitution of the metal oxide support, since the Pt-OMS-2 catalyst retained activity for about 33 h before deactivation, after which nanoparticles were observed in TEM images. This work offers a new perspective in single atom synthesis using the metal oxide as the host for the single atom site, instead of adatoms on the surface

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrological control of chemical weathering and rock-carbon fluxes: East River, Colorado

This data package is used in the manuscript entitled “Hydrological control of chemical weathering and rock-carbon fluxes”. The field study was conducted in a lower montane hillslope of the East River watershed, underlain by Mancos Shale, within the lower 140 m section of a transect that extends nearly 1 km to its local peak. The data in this package, in CSV file format, were collected from Fall 2016 to fall 2021, including depth-resolved dynamic water table depths, subsurface water fluxes, solid phase (soil to bedrock) chemical compositions, geochemical properties of porewater and pore-gas, including radio carbon concentrations. The detailed methods of field studies, laboratory chemical analyses, and calculations are described in the Methods section and in the dataset file: Wan_et_al_Methods.pdf. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.

54 ENVIRONMENTAL SCIENCES↗

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

Direct study of changes in catalyst structure-kinetic properties during redox transitions using a new time-resolved technique

Direct study of changes in catalyst structure-kinetic properties during redox transitions using a new time-resolved technique The Temporal Analysis of Products (TAP) pulse response methodology1 is a transient technique offering the time scale needed to deconvolve many reaction steps from the complex network typical to industrial catalytic processes. The traditional TAP measurement observes the gas phase dynamic response but, hereto now, has been devoid of any direct measurement of change in the catalyst itself. The development of an operando technique that couples gas phase transient kinetics to dynamic metal centers and surface species is presented. Significance Unification of the TAP methodology with time-resolved spectroscopic measurement can offer unprecedented insight into the complex kinetic phenomena regulated by the solid catalyst.

03 - NATURAL GAS↗

Probing Attosecond Bound Electron Dynamics Driven by Strong-Field Light Transients

The project aimed to drive strong-field excitation of bound electronic coherences in gas-phase atoms and molecules using a novel source of near-single-cycle optical field transients, based on industrial-grade Yb:solid-state laser technology. Key objectives included generating near-single cycle field transients and isolated attosecond pulses, reconstructing time-dependent wave packets in atoms and molecules, probing their decoherence using attosecond interferometry, and controlling the dynamic evolution of systems on attosecond to few-femtosecond timescales. A secondary goal of the project was to develop field-resolved measurement platforms as a novel approach for studying electron dynamics in atomic and molecular systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Let’s Stick Together: Interplay of Moisture, Particle Size, and Anatomical Fraction in the Flowability of Corn Stover Crumbles

Handling of agricultural waste biomass feedstocks (such as corn stover) is a persistent challenge in the production of biofuels and chemicals at integrated biorefineries. Inconsistent flowability and high feedstock variability creates equipment downtime and processing challenges. This work investigates the influence of anatomical fractionation and moisture content on the bulk solids handling behavior of corn stover crumbles. A combination of rotational shear testing, tribological measurements, and surface and sorption characterization techniques (inverse gas chromatography (IGC), dynamic vapor sorption (DVS), and electron microscopy) were used to evaluate the contributions of particle composition, size, and surface properties to flow behavior. Shear and frictional testing revealed that moisture content has the greatest effect on flowability, significantly increasing unconfined yield strength (ƒ c ) and reducing the flow function coefficient (F F C), regardless of anatomical fraction. While differences in surface energy were observed between fractions, particularly in polar contributions (γ AB ), these did not correspond to meaningful differences in flow behavior. Coarser 6 mm particles exhibited poorer flow performance than 4 mm crumbles, likely due to broader particle size distributions and increased particle interlocking. Taken together, these results suggest that anatomical fractionation provides minimal benefit from a bulk handling perspective. Moisture content and particle size heterogeneity dominate handling behavior, reinforcing the utility of unfractionated corn stover in processing environments. These insights can inform the design of more efficient and reliable feedstock processing systems for agricultural biomass in integrated biorefineries.

09 BIOMASS FUELS↗

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗