Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh compression”

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.

540 records · Page 30

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Experimental Study on Portable Air-Conditioning System with Enhanced PCM Condenser

A portable personal air-conditioning system with a vapor compression cycle (VCC), electrical battery, and thermal battery is an attractive option for localized cooling as compared with stationary air-conditioning systems. A phase change material (PCM) can be used for storing waste heat from the VCC condenser, thereby acting as a thermal battery. In the PCM condenser, copper tubes are submerged in the PCM container. Due to the low thermal conductivity of pure PCM, the overall heat transfer coefficient of the PCM condenser is limited resulting in high VCC condensing pressure and thus, low VCC coefficient of performance (COP). Therefore, it is necessary to improve VCC system performances by applying heat transfer augmentation in PCM. In this work, a portable air-conditioning system with heat transfer enhanced PCM that stores heat from condenser during VCC operation and is regenerated by a thermosyphon operation is proposed. The system also includes a receiver to address the charge difference between the thermosyphon and the VCC operation. The paper experimentally studies the effects of the following three types of enhanced PCM condensers on the system performance: (1) copper mesh enhanced PCM with helically coiled copper tube branches and large refrigerant distribution headers (CMPHX), (2) graphite matrix enhanced PCM with straight copper tube branches and small headers (GMPHX1), and (3) graphite matrix enhanced PCM with straight copper tube branches and large headers (GMPHX2). Results show that this system operates properly and effectively during both VCC and PCM regeneration operations. The system can provide four hours of continuous cooling, and the PCM can be fully regenerated after 7.4 hours by thermosyphon operation. The comparisons of three PCM condensers indicate that the system with GMPHX2 performs the best with the average VCC COP of 4.7. The VCC system with CMPHX has an average COP of 4.2, and GMPHX1 has an average COP of 3.1. GMPHX2 has the best heat transfer enhancement with an overall heat transfer coefficient of 632.1 W∙m-2∙K-1, which is more than twice of GMPHX1 and five times of CMPHX. Moreover, GMPHX2 can achieve more uniform refrigerant flows and PCM temperature distributions than GMPHX1. This leads to the higher utilization ratio of PCM latent heat. CMPHX is the second best but needs larger heat transfer surface area and inner volume as well as more refrigerant charge. In conclusion, the system design with GMPHX2 is recommended due to its small refrigerant charge, high VCC COP, and the use of good heat transfer enhancement material.

Qiao, Yiyuan↗

GAS ACCRETION IS DOMINATED BY WARM IONIZED GAS IN MILKY WAY MASS GALAXIES AT z {approx} 0

We perform high-resolution hydrodynamic simulations of a Milky Way mass galaxy in a fully cosmological setting using the adaptive mesh refinement code, Enzo, and study the kinematics of gas in the simulated galactic halo. We find that the gas inflow occurs mostly along filamentary structures in the halo. The warm-hot (10{sup 5} K 10{sup 6} K) ionized gases are found to dominate the overall mass accretion in the system (with M-dot = 3-5 M {sub Sun} yr{sup -1}) over a large range of distances, extending from the virial radius to the vicinity of the disk. Most of the inflowing gas (by mass) does not cool, and the small fraction that manages to cool does so primarily close to the galaxy (R {approx}< 100 kpc, with more pronounced cooling at smaller R), perhaps comprising the neutral gas that may be detectable as, e.g., high-velocity clouds. The neutral clouds are embedded within larger, accreting filamentary flows, and represent only a small fraction of the total mass inflow rate. The inflowing gas has relatively low metallicity (Z/Z {sub Sun} < 0.2). The outer layers of the filamentary inflows are heated due to compression as they approach the disk. In addition to the inflow, we find high-velocity, metal-enriched outflows of hot gas driven by supernova feedback. Our results are consistent with observations of halo gas at low z.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of internal and near-nozzle flow for a GDI fuel injector

A numerical study of two-phase flow inside the nozzle holes and the issuing spray jets for a multi-hole direct injection gasoline injector has been presented in this work. The injector geometry is representative of the Spray G nozzle, an eight-hole counterbore injector, from the Engine Combustion Network (ECN). Simulations have been carried out for a fixed needle lift. Effects of turbulence, compressibility and non-condensable gases have been considered in this work. Standard k -ε turbulence model has been used to model the turbulence. Homogeneous Relaxation Model (HRM) coupled with Volume of Fluid (VOF) approach has been utilized to capture the phase change phenomena inside and outside the injector nozzle. Three different boundary conditions for the outlet domain have been imposed to examine non-flashing and evaporative, non-flashing and non-evaporative and flashing conditions. Noticeable hole-to-hole variations have been observed in terms of mass flow rates for all the holes under all the operating conditions considered in this study. Inside the nozzle holes mild cavitation-like and in the near-nozzle region flash boiling phenomena have been predicted when liquid fuel is subjected to superheated ambiance. Under favorable conditions considerable flashing has been observed in the near-nozzle regions. An enormous volume is occupied by the gasoline vapor, stantial computational cost. Volume-averaging instead of mass-averaging is observed to be more effective, especially for finer mesh resolutions.

42 ENGINEERING↗

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Self-consistent clustering analysis: an efficient multiscale scheme for inelastic heterogeneous materials

A predictive computational theory is shown for modeling complex, hierarchical materials ranging from metal alloys to polymer nanocomposites. The theory can capture complex mechanisms such as plasticity and failure that span across multiple length scales. This general multiscale material modeling theory relies on sound principles of mathematics and mechanics, and a cutting-edge reduced order modeling method named self-consistent clustering analysis (SCA) [Zeliang Liu, M.A. Bessa, Wing Kam Liu, “Self-consistent clustering analysis: An efficient multi-scale scheme for inelastic heterogeneous materials,” Comput. Methods Appl. Mech. Engrg. 306 (2016) 319–341]. SCA reduces by several orders of magnitude the computational cost of micromechanical and concurrent multiscale simulations, while retaining the microstructure information. This remarkable increase in efficiency is achieved with a data-driven clustering method. Computationally expensive operations are performed in the so-called offline stage, where degrees of freedom (DOFs) are agglomerated into clusters. The interaction tensor of these clusters is computed. In the online or predictive stage, the Lippmann-Schwinger integral equation is solved cluster-wise using a self-consistent scheme to ensure solution accuracy and avoid path dependence. To construct a concurrent multiscale model, this scheme is applied at each material point in a macroscale structure, replacing a conventional constitutive model with the average response computed from the microscale model using just the SCA online stage. A regularized damage theory is incorporated in the microscale that avoids the mesh and RVE size dependence that commonly plagues microscale damage calculations. The SCA method is illustrated with two cases: a carbon fiber reinforced polymer (CFRP) structure with the concurrent multiscale model and an application to fatigue prediction for additively manufactured metals. For the CFRP problem, a speed up estimated to be about 43,000 is achieved by using the SCA method, as opposed to FE2, enabling the solution of an otherwise computationally intractable problem. The second example uses a crystal plasticity constitutive law and computes the fatigue potency of extrinsic microscale features such as voids. This shows that local stress and strain are capture sufficiently well by SCA. This model has been incorporated in a process-structure-properties prediction framework for process design in additive manufacturing.

Liu, Z.↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Improved Thermoplastic/Iron-Particle Transformer Cores

A method of fabricating improved transformer cores from composites of thermoplastic matrices and iron-particles has been invented. Relative to commercially available laminated-iron-alloy transformer cores, the cores fabricated by this method weigh less and are less expensive. Relative to prior polymer-matrix/ iron-particle composite-material transformer cores, the cores fabricated by this method can be made mechanically stronger and more magnetically permeable. In addition, whereas some prior cores have exhibited significant eddy-current losses, the cores fabricated by this method exhibit very small eddy-current losses. The cores made by this method can be expected to be attractive for use in diverse applications, including high-signal-to-noise transformers, stepping motors, and high-frequency ignition coils. The present method is a product of an experimental study of the relationships among fabrication conditions, final densities of iron particles, and mechanical and electromagnetic properties of fabricated cores. Among the fabrication conditions investigated were molding pressures (83, 104, and 131 MPa), and molding temperatures (250, 300, and 350 C). Each block of core material was made by uniaxial-compression molding, at the applicable pressure/temperature combination, of a mixture of 2 weight percent of LaRC (or equivalent high-temperature soluble thermoplastic adhesive) with 98 weight percent of approximately spherical iron particles having diameters in the micron range. Each molded block was cut into square cross-section rods that were used as core specimens in mechanical and electromagnetic tests. Some of the core specimens were annealed at 900 C and cooled slowly before testing. For comparison, a low-carbon-steel core was also tested. The results of the tests showed that density, hardness, and rupture strength generally increased with molding pressure and temperature, though the correlation was rather weak. The weakness of the correlation was attributed to the pores in the specimens. The maximum relative permeabilities of cores made without annealing ranged from 30 to 110, while those of cores made with annealing ranged from 900 to 1,400. However, the greater permeabilities of the annealed specimens were not associated with noticeably greater densities. The major practical result of the investigation was the discovery of an optimum distribution of iron-particle sizes: It was found that eddy-current losses in the molded cores were minimized by using 100 mesh (corresponding to particles with diameters less than or equal to 100 m) iron particles. The effect of optimization of particle sizes on eddy-current losses is depicted in the figure.

Wincheski, Russell A.↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

ASU’s DAC polymer-enhanced cyanobacterial bioproductivity (AUDACity)

ASU’s DAC polymer-enhanced cyanobacterial bioproductivity (AUDACity) project aims to demonstrate a novel, scalable method for removing carbon dioxide (CO 2 ) directly from ambient air and delivering it to cyanobacterial cultures to produce commodity biofuel, mid-value protein for supplements, and high value phycocyanin (PC), a natural blue colorant (Figure A). This approach uses low-cost, reusable anion exchange polymers embedded in modular mesh packets, which capture CO 2 during drying cycles when exposed to ambient air, and release concentrated CO 2 into aqueous cultivation systems. The project addresses a critical challenge in energy research needed for developing sustainable, economically viable methods of Direct Air Capture (DAC) that can be integrated with bio-based systems for fuel and chemical production. AUDACity contributes to scientific understanding by integrating materials chemistry, cyanobacterial biology, and system engineering to create a distributed CO 2 delivery platform. Key insights have emerged around the design of biocompatible sorbents, optimization of CO 2 capture-release cycles, and durability of packet-based delivery systems under outdoor conditions. Notably, the team has synthesized and tested a range of polymer sorbents, identified mechanisms of material degradation and fouling, and advanced both lab- and pilot-scale cultivation systems to evaluate performance. From a technical and economic standpoint, AUDACity shows promise for achieving cost-effective CO 2 capture and delivery into aqueous media and biofuel production. Preliminary techno-economic analysis (TEA) indicates that the DAC system based on current performance can reach $\$$680/tonne CO 2 delivered into aqueous solution; with reasonable improvements to sorbent lifetime, sorbent capacity, reducing water uptake the approach could reach $\$$66/tonne by avoiding the need for energy-intensive sorbent regeneration and CO 2 compression, making it more feasible for decentralized deployment. With these costs for CO 2 and by extracting and selling high-value PC ($\$$50/kg) and mid-value protein supplement ($\$$6/kg), the remaining biomass can be hydrothermally treated into biofuel for $\$$2.50/gallon, and would support a small first-of-a-kind biorefinery capable of producing 500 barrels per day of biofuel. The project offers meaningful public benefits by advancing carbon removal technologies that are low-energy, modular, and adaptable to non-arable land and brackish water use. It aligns with national goals to develop advanced biotechnology and supports future pathways for bio-based fuels and products. By enabling direct coupling of CO 2 transfer into aqueous medium and biological carbon utilization, AUDACity lays the groundwork for effective algae cultivation without wasteful CO 2 delivery and is a promising and innovative solution for low-carbon fuel and bioproduct generation contributing to a vigorous bioeconomy.

09 BIOMASS FUELS↗

A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

In this paper, the compressible Rayleigh–Taylor (RT) instability is studied by performing a suite of large eddy simulations (LES) using the Miranda and Ares codes. A grid convergence study is carried out for each of these computational methods, and the convergence properties of integral mixing diagnostics and late-time spectra are established. A comparison between the methods is made using the data from the highest resolution simulations in order to validate the Ares hydro scheme. We find that the integral mixing measures, which capture the global properties of the RT instability, show good agreement between the two codes at this resolution. The late-time turbulent kinetic energy and mass fraction spectra roughly follow a Kolmogorov spectrum, and drop off as k approaches the Nyquist wave number of each simulation. The spectra from the highest resolution Miranda simulation follow a Kolmogorov spectrum for longer than the corresponding spectra from the Ares simulation, and have a more abrupt drop off at high wave numbers. The growth rate is determined to be between around 0.03 and 0.05 at late times; however, it has not fully converged by the end of the simulation. Finally, we study the transition from direct numerical simulation (DNS) to LES. The highest resolution simulations become LES at around t/τ ≃ 1.5. Finally, to have a fully resolved DNS through the end of our simulations, the grid spacing must be 3.6 (3.1) times finer than our highest resolution mesh when using Miranda (Ares).

42 ENGINEERING↗

Trabecular Bone Mechanical Properties and Fractal Dimension

Countermeasures for reducing bone loss and muscle atrophy due to extended exposure to the microgravity environment of space are continuing to be developed and improved. An important component of this effort is finite element modeling of the lower extremity and spinal column. These models will permit analysis and evaluation specific to each individual and thereby provide more efficient and effective exercise protocols. Inflight countermeasures and post-flight rehabilitation can then be customized and targeted on a case-by-case basis. Recent Summer Faculty Fellowship participants have focused upon finite element mesh generation, muscle force estimation, and fractal calculations of trabecular bone microstructure. Methods have been developed for generating the three-dimensional geometry of the femur from serial section magnetic resonance images (MRI). The use of MRI as an imaging modality avoids excessive exposure to radiation associated with X-ray based methods. These images can also detect trabecular bone microstructure and architecture. The goal of the current research is to determine the degree to which the fractal dimension of trabecular architecture can be used to predict the mechanical properties of trabecular bone tissue. The elastic modulus and the ultimate strength (or strain) can then be estimated from non-invasive, non-radiating imaging and incorporated into the finite element models to more accurately represent the bone tissue of each individual of interest. Trabecular bone specimens from the proximal tibia are being studied in this first phase of the work. Detailed protocols and procedures have been developed for carrying test specimens through all of the steps of a multi-faceted test program. The test program begins with MRI and X-ray imaging of the whole bones before excising a smaller workpiece from the proximal tibia region. High resolution MRI scans are then made and the piece further cut into slabs (roughly 1 cm thick). The slabs are X-rayed again and also scanned using dual-energy X-ray absorptiometry (DEXA). Cube specimens are then cut from the slabs and tested mechanically in compression. Correlations between mechanical properties and fractal dimension will then be examined to assess and quantify the predictive capability of the fractal calculations.

Hogan, Harry A.↗

A Lagrangian discontinuous Galerkin hydrodynamic method

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for solving the two-dimensional gas dynamic equations on unstructured hybrid meshes. The physical conservation laws for the momentum and total energy are discretized using a DG method based on linear Taylor expansions. Three different approaches are investigated for calculating the density variation over the element. The first approach evolves a Taylor expansion of the specific volume field. The second approach follows certain finite element methods and uses the strong mass conservation to calculate the density field at a location inside the element or on the element surface. The third approach evolves a Taylor expansion of the density field. The nodal velocity, and the corresponding forces, are explicitly calculated by solving a multidirectional approximate Riemann problem. An effective limiting strategy is presented that ensures monotonicity of the primitive variables. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. Results from a suite of test problems are presented to demonstrate the robustness and expected second-order accuracy of this new method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗

Petrophysical and geomechanical properties of gas hydrate-bearing sediments recovered from Alaska North Slope 2018 Hydrate-01 Stratigraphic Test Well

Knowledge of petrophysical and geomechanical properties of gas hydrate-bearing sediments are essential for predicting reservoir response due to gas production from gas hydrate reservoirs. That information will be critical parameters for designing production well completion such as specification of depressurization pump, water storage tank, and mesh size of sand screen. In December 2018, Stratigraphic Test Well Hydrate-01 was drilled in the western part of the Prudhoe Bay Unit, Alaska North Slope as part of the technical planning effort for a future long-term production test by collaborative team of DOE/NETL, USGS, and MH21-S (Boswell et al., 2020, Collett et al., 2020, Okinaka et al., 2020). Data set of logging-while-drilling (LWD) were acquired (Haines et al., 2020, Suzuki et al., 2019) and core sampling depth was determined on-site.Side-wall pressure coring was conducted to recover gas hydrate-bearing sediments from two reservoir sections named Unit-B and Unit-D. A total of 34 cores were successfully recovered by 5 runs of a wire-line deployed pressure corer (CoreVault® System - Halliburton). Core analysis scheme of this project are shown in Figure 1. All cores were quenched in liquid nitrogen while at high pressure in the laboratory of Stratum Reservoir, LLC. at Anchorage (Figure 1, a)). And the cores were removed from the pressure corer autoclave with temperature support of dry ice and stored under liquid nitrogen at atmospheric pressure. 19 damaged cores were processed to index property measurements directly such as grain size, grain density. 4 of another 17 cores were depressurized and trimmed for making a plug to analyze petrophysical properties of host sediments. Unsteady-state permeability test was conducted to obtain relative water permeability to gas and core scale NMR T2 distribution measurement was performed for evaluating pore size distribution at Houston (Figure 1, b)). Remained high quality 13 cores were preserved with gas hydrate for advanced laboratory analysis. National Institute of Advanced Industrial Science and Technology, as a part of the Japanese National Hydrate Research Program (MH21-S, funded by Ministry of Economy, Trade and Industry), received the samples at Sapporo, Japan for advanced core analysis. High-resolution X-ray CT was used to analyze the quality of the samples, which showed undisturbed lithological layers. Cores were lathed into cylindrical shape and distributed for multi property measurements (Figure 1, c)).<p>As a result, sediment from Unit-D is categorized as silty sand at ~37% porosity with ~80% gas hydrate saturation. An average hydration number n = 6.16 was measured by Raman spectroscopy. An average intrinsic permeability of ~400 mD and in situ effective permeability (with hydrate) on the order of ~10 mD. The Unit B recovered cores consisted of well sorted sand at ~40% porosity with ~95% gas hydrate saturation. An average intrinsic permeability of ~1 Darcy and in situ effective permeability on the order of ~30 mD was measured for the Unit B cores. There was a small permeability reduction due to porosity loss with increasing effective stress that simulated consolidation behavior along with depressurization in the highly permeable sandy sediment. The apparent minimum change in porosity and permeability may be caused by the low compressibility of quartz sand grains in the recovered cores. XRD and thermal conductivity analysis also suggested high quartz content. Triaxial compression tests established internal friction angles based on the Mohr-Coulomb's failure criterion, which are 40° for hydrate-bearing sediment and 29.8° for hydrate free sediment.</p>

Yoneda, Jun↗

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Emissions mitigation technology for advanced water-lean solvent-based CO 2 capture processes

This technical final report submitted to DOE/NETL presents all the research activities performed during the entirety of DE-FE0031660 project-Emissions Mitigation Technology for Advanced Water-Lean Solvent-Based CO 2 Capture Processes which spans from October 2018 through March 2022. RTI International has been conducting studies from fundamental and operational aspects to reduce the overall amine emissions from the advanced Water-Lean Solvent (WLS) systems, specifically RTI’s Non-Aqueous Solvent (NAS). This technical final report will highlight the key findings from project which align closely to the project objectives which are: Identify the contribution of vapor loss, entrainment, and aerosols to the overall emissions of water-lean systems; Determine the significance of CO 2 capture system operating parameters to the amine emissions; Develop an emissions model based on critical operating parameters; Evaluate the effectiveness of emissions mitigation devices to reduce the amine emissions to <1 ppm under flue coal-fired flue gas; and, Determine the contribution of the ECTs to the overall CO 2 capture cost. The following are the key findings based on numerous tests using both lab-scale setups and parametric testing performed at RTI’s Bench-scale Gas Absorption System (BsGAS). During the BP1, the aerosol generation system and monitoring equipment were installed at BsGAS to produce and determine the aerosol characteristics during the NAS CO 2 capture process. The aerosol produced by this setup produced aerosols with the peak diameter and concentration of 50 micron and 1.2E10 7 cm -3 , respectively. These particle sizes and concentrations are matched to those observed in the actual coal-fired power plant flue gases and expected to be found at the absorber inlet of the CO 2 capture system. Over 1,300 hours of parametric testing have been conducted to evaluate the impact of the aerosols and operating conditions during the CO 2 capture with NAS on the overall amine emissions in the treated flue gas. At the worse condition tested, the presence of the aerosols in the flue gas could increase the overall emissions by 10X compared to the baseline emissions from NAS’s vapor pressure. CO 2 capture rate was found to be a main factor impacting the overall emissions as well as aerosol size and concentrations in the absorber off-gas. The higher CO 2 capture rate, the higher amine emissions in the treated gas. The temperature difference between the temperature bulge seen in the absorber and the water wash temperature also impacts the particle growth where the larger the temperature difference, the more amine emissions from aerosols in the treated gas. The majority of the aerosols did not grow substantially in the system, and the particle concentrations remained nearly constant between the absorber inlet and wash outlet. Only a small portion of the particles were found to grow significantly. The high efficiency demister with mesh size of 5-10 micron can be installed to remove a portion of the aerosols from the gas stream leaving the water wash. Overall, these results from parametric testing have established the emission baseline and validate our assumption on the need of emission control technologies (ECT) in order to minimize the emissions from the baseline NAS CO 2 capture process. Over 2,000 of BsGAS operating hours was used to investigate a handful of process improvements which led to a selection of the vital few changes that effectively control the amine emissions. These process improvements are lime-coated-filters for absorber gas inlet, advanced demister at the top of the absorber, a second water wash with amine recovery unit were designed, installed, and tested at BsGAS at the end of BP1. The result showed that the NAS CO 2 capture process with these additional emission control devices could lower the amine emission in the treated gas to about 1 ppm using a simulated coal-fire flue gas stream. The main contributor in lowering the amine emission came from the second water wash with amine recovery unit where the amine concentration in the scrubbing water was kept below 2 wt% through a continuous amine removal via an adsorbent bed, resulting in a low amine vapor pressure. The adsorbent bed was regenerated via a direct steam regeneration and the recover amine was returned to the absorber to minimize wastewater and makeup amine. A flue gas generation system was designed and installed during the first half of BP2 to support the emission testing using a real coal-derived flue gas. The system is capable of generating both coal- and natural gas- derived flue gases with the composition of the gaseous species highly resemble to that of the power plant flue gases. The particulates detected in the coal-derived flue gas showed the mean diameter of 1 micron. The CO 2 capture operating was then proceed using the real coal-derived flue gas where the amine emission was controlled to be about 0-3 ppm for the total run time of about 200 hours. Similar testing was conducted with natural gas-derived flue gas and the result showed a highly amine emission of 30 ppm under the total run time of 200 hours. The Principal Component Analysis (PCA) and the Partial Least Squares Projection to Latent Structures (PLS) techniques were applied to the parametric testing data to derive a multivariate statistical model. The model was validated and trained with half of the data collected, and the predictive ability of the model was evaluated using the remaining half of the data. The resulting empirical model was capable of predicting the overall emissions from the NAS process without the ECTs with ±15% accuracy (average absolute deviation, AAD) in BP1. As more emission data were obtained under the real coal-flue gas in the BP2, the model incorporated these new set of data to reflect the final process configuration, operating parameters, and amine emission. This results in the updated empirical model predicting the amine emission from the NAS CO 2 capture process with 84% goodness-of-fit (R 2 ), 85% predictability (Q 2 ), and 15% AAD. The study evaluates the use of RTI’s Non-Aqueous Solvent technology for 90% CO 2 capture from a net 650 MWe pulverized coal power plant, downstream of the flue-gas desulfurization unit. The captured CO 2 has a purity of > 95% CO 2 , and is dried, compressed to 15.3 MPa (2,215 psia), ready for sequestration. The analysis uses Case B12B from the DOE Baseline study on Bituminous Coal, Revision 4 where the Cansolv CO 2 capture plant is replaced by the RTI CO 2 Capture plant. The CO 2 capture plant has been sized to capture >90% CO 2 from flue gas derived from a net 650 MWe supercritical pulverized coal power plant. The CO 2 capture plant is equipped with emission control technologies that limits the amine emissions to < 1 ppm. Two different cases were evaluated for the technoeconomic study. The key difference between the two cases is the regenerator pressure. In Case 1, the regenerator operates at 0.195 MPa (28.3 psia), whereas in Case 2, the regenerator pressure is 0.44 MPa (64 psia) thus removing the need for the first stage of compression of the eight-stage compression train. Results from the TEA are compared against the DOE reference cases for SCPC plant with and without CO 2 Capture (Case B12A and Case B12B of the DOE Baseline study, respectively). Case 2 with CO 2 regeneration at higher pressure results in the lower cost of CO 2 capture. The total capital cost of the capture process has been estimated using 2018 dollars in Aspen Process Economic Analyzer and was estimated to be $579 MM. The capture plant operation leads to a total parasitic power loss rate of 96 MWe, resulting in a decrease in pulverized coal power plant efficiency of 7.8% points. The resulting cost of electric power increases from 64.4 mills/kWh, for no capture, to 97.5 mills/kWh, with 90% capture, an increase of 51% in the COE. The cost of capturing 90% CO 2 was estimated to be $38.2/tonne-CO 2 , and meets the DOE target of $40/t-CO 2 . Emission control technologies (ECT) investigated in this project includes a second water wash with use of activated carbon beds for removal of amine from the wash water prior to recirculation in the water wash. These ECT allow operation of the CO 2 capture plant with < 1 ppm amine emissions with the treated flue gas and contributes to $2.4/t-CO 2 captured. Amine emissions derived from thermal and oxidative degradations were investigated under this project along with the emissions derived from aerosols for the NAS system. The thermally degraded of the lean NAS showed less than 4% decreased of the original total amine content in the NAS at 150 °C while the result obtained at 120 °C showed no drop in total amine content, suggesting that thermal degradation of the NAS is minimal. These results also suggested that the thermally degraded species are not likely formed and contributed to the emissions due to the low regeneration temperature of the NAS at 90-105 °C. The oxidative degradation, on the other hand, could become problematic as some of these oxidative degraded species were observed during the NAS-5 testing at National Carbon Capture Center (NCCC) and SINTEF in our previous project. The rapid screening of selected inhibitors suggested that oxidative degradation of NAS can be suppressed using thiol containing compounds in amounts of at least 1 mol%. The detailed mechanistic degradation pathway was conceived for a specific amine used in NAS formulation during BP2. he reduction of the nitrosamines caused by the NO x present in the flue gas was also examined. The study suggested that the thermo-chemical treatment of the NAS solvent would be a more effective and economically viable compared to removing NO x at the DCC.

01 COAL, LIGNITE, AND PEAT↗