Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “low pressure”

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 127 records · Page 7

Investigating the Influence of Hexanuclear Clusters in Isostructural Metal–Organic Frameworks on Toxic Gas Adsorption

The efficient capture of toxic gases, such as ammonia (NH 3 ) and sulfur dioxide (SO 2 ), can protect the general population and mitigate widespread air pollution. Metal–organic frameworks (MOFs) comprise a tunable class of adsorbents with high surface areas that can meet this challenge by selectively capturing these gases at low concentrations. In this work, we explored how modifying the metal ions in the node of an isostructural MOF series from a transition metal to a lanthanide or actinide influences the electronic environment of the node-based active site. Next, we investigated the adsorption properties of each MOF toward the relatively basic NH 3 and relatively acidic SO 2 gases. Within the NU-907 family of MOFs, we found that Zr 6 -NU-907 exhibits the best uptake toward NH 3 at low pressures, while Th 6 -NU-907 demonstrates the best low-pressure performance for SO 2 adsorption. Here, tracking the infrared (IR) stretching frequency of the node-based μ 3 -OH groups provides insights into the electronegativity of the metal ion and suggests that the most electronegative metal ion (Zr) affords the node with the best NH 3 uptake at low pressures. In contrast, the Th 6 node contains additional coordinated water groups relative to the other M 6 nodes, which appears to yield the MOF with the greatest affinity for SO 2 uptake that occurs predominately through reversible physisorption interactions. Finally, in situ NH3 IR spectroscopic studies indicate that both NH 4 + and Lewis-bound NH 3 species form during adsorption. Combined, these results suggest that tuning the electronic properties and structure of the node-based active site in an MOF presents a viable strategy to change the affinity of an MOF toward toxic gases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Activity coefficients of siderophile elements in Fe-Si liquids at high pressure

Metallic core formation in differentiated bodies in the inner solar system can take place between low pressures (near 1 bar) to much higher pressures (up to 100 GPa). Most thermodynamic models of metal-silicate equilibria utilise activity coefficients for metallic tracers in Fe liquids, nearly all of which have been carried out at low pressures. This study focuses on the effect of pressure on activity coefficients for Au, P, V, Mn, Ga, Zn, Cd, Sn, W, Pb, and Nb in liquid Fe-Si alloys. From a series of experiments at 10 GPa, 2373 K containing variable Si content in a metallic liquid we have derived epsilon interaction parameters in Fe-Si liquids ($ε^{Si}_{M}$). Comparison of 1 GPa and 10 GPa data shows no difference except for Nb. Epsilon parameters derived from low pressure experiments can thus be used to calculate activity coefficients for application to higher pressure processes (at least to 10 GPa).

58 GEOSCIENCES↗

Understanding the Mechanisms of Self-Organization in Atmospheric Pressure DC Glows

One of the most challenging areas in atmospheric pressure plasmas is controlling uniformity and contact area. Indeed it is this problem which is closely related to scale up that has impeded the widespread applications of these plasmas. The atmospheric pressure plasma glow is a type of atmospheric pressure discharge akin to a traditional low pressure glow particularly when compared with discharge morphology. But this is where the similarities end. Unlike low pressure glows, atmospheric pressure glows do not rely on electrode processes such as secondary electron emission to make the discharge self-sustaining, rather the discharge is sustained by typically processes such as the presence of a resistive layer or through favorable atomic collisional processes such as that which prevail in helium that stabilizes the discharge so that the thermal instability is mitigated. These discharges can feature either solid metal electrodes or geometries where at least one electrode is a liquid water electrolyte. Regardless of the implementation, DC 1 Atm glows are fascinating in that the plasma attachment under certain conditions can self-organize. The attachment surface area therefore changes as a result. Therefore, self-organization could be the key to optimizing plasma contact area and thereby open the door to applications in the real world. The plasma physical processes prevailing at the surface where the self-organization takes place remain largely unexplored. There is very little published research on the conditions that prevail at the surface and virtually no scholarship on measured plasma parameters of the self-organized attachments themselves or the surface electric field. This effort breaks new ground by elucidating these largely unexplored physical processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic ordering in terbium at high pressures and low temperatures

We have investigated magnetic ordering in the heavy lanthanide terbium (Tb) for pressures up to 70 GPa and temperatures down to 20 K by synchrotron x-ray diffraction at the Advanced Photon Source, Argonne National Laboratory. The x-ray diffraction studies were complemented by neutron diffraction studies at the Spallation Neutron Source, Oak Ridge National Laboratory. We focused on the higher-pressure phases of Tb, such as the double hexagonal close-packed (dhcp, above 16 GPa), distorted face-centered cubic (hR24, above 30 GPa), and orthorhombic (oF16, above 51 GPa) phases. We observed the phenomenon of spontaneous magnetostriction in the dhcp and hR24 phases, where magnetic ordering at low temperatures gives rise to subtle splitting in x-ray diffraction peaks. The neutron diffraction study of the dhcp phase at 20 GPa revealed magnetic peaks that can be assigned a propagation vector k = ($\frac{1}{2}$, 0, $\frac{1}{2}$), and a magnetic ordering temperature that is consistent with the magnetostriction effect observed with x-ray diffraction. The high-pressure and low-temperature x-ray diffraction study at 70 GPa in the collapsed oF16 phase did not show any evidence of magnetostriction, indicating a significant decrease in the magnetic moment of the 4f-shell due to band broadening at ultra-high compression.

36 MATERIALS SCIENCE↗

Large Eddy Simulation of Rotationally Induced Ingress and Egress around an Axial Seal between Rotor and Stator Disks

In gas turbines, the hot gas exiting the combustor can have temperatures as high as 2000 °C, and some of this hot gas enter into the space between the stator and rotor disks (wheelspace). Since the entering hot gas could damage the disks, its ingestion must be minimized. This is carried out by rim seals and by introducing a cooler flow from the compressor (sealing flow) into the wheelspace. Ingress and egress into rim seals are driven by the stator vanes, the rotor and its rotation, and the rotor blades. This study focuses on the ingress and egress driven by the rotor and its rotation. This is carried out by performing wall-resolved large eddy simulation (LES) around an axial seal in a rotor–stator configuration without vanes and blades. Results obtained show the mechanisms by which the rotor and its rotation induce ingress, egress, and flow trajectories. Kelvin–Helmholtz instability was found to create a wavy shear layer and displacement thickness that produces alternating regions of high and low pressures around the rotor side of the seal. Vortex shedding on the backward-facing side of the seal and its impingement on the rotor side of the seal also produces alternating regions of high and low pressures. The locations of the alternating regions of high and low pressures were found to be statistically stationary and to cause ingress to start on the rotor side of the seal. Vortex shedding and recirculating flow in the seal clearance also cause ingress by entrainment. With the effects of the rotor and its rotation on ingress and egress isolated, this study enables the effects of stator vanes and rotor blades to be assessed.

42 ENGINEERING↗

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation-aided characterization of a versatile water-based condensation particle counter for atmospheric airborne research

Capturing the vertical profiles and horizontal variations of atmospheric aerosols often requires accurate airborne measurements. With the advantage of avoiding health and safety concerns related to the use of butanol or other chemicals, water-based condensation particle counters have emerged to provide measurements under various environments. However, airborne deployments are relatively rare due to the lack of instrument characterization under reduced pressure at flight altitudes. This study investigates the performance of a commercial “versatile” water-based condensation particle counter (vWCPC, model 3789, TSI, Shoreview, MN, USA) under various ambient pressure conditions (500–920 hPa) with a wide range of particle total number concentrations (1500–70 000 cm -3 ). The effect of conditioner temperature on vWCPC 3789 performance at low pressure is examined through numerical simulation and laboratory experiments. We show that the default instrument temperature setting of 30 °C for the conditioner is not suitable for airborne measurement and that the optimal conditioner temperature for low-pressure operation is 27°. Under the optimal conditioner temperature (27°), the 7 nm cut-off size is also maintained. Additionally, we show that insufficient droplet growth becomes more significant under the low-pressure operation. The counting efficiency of the vWCPC 3789 can vary up to 20 % for particles of different chemical compositions (e.g., ammonium sulfate and sucrose particles). However, such variation is independent of pressure.

54 ENVIRONMENTAL SCIENCES↗

Low-Cost Preform and Molding Processes

The entry cost for prototyping a composite component for manufacture using automated, high rate processes is prohibitively expensive in many cases, especially for small business, where tooling costs may be several $100k. Discussions with industry also indicate that many small companies, tier 1 and 2 suppliers, have an interest to mold composite parts but do not want to deal with the capital cost, material handling issues, and labor associated with dry fiber preforming operations. While the molders may locate near the end user for logistics reasons, it may be more cost effective for the performer to remain regional and invest in capital equipment to support preform automation, thus keeping costs to a minimum. This project was designed to explore and demonstrate several options to meet these industry needs. Dry fiber preforming approaches were evaluated which allow for low pressure resin infusion, single sided tooling options such a vacuum assisted resin transfer molding (VARTM) or low pressure resin transfer molding (RTM-light). Unlike sheet molding compound, SMC compression molding where typical molding pressures of 1000 psi are required to push material into the desired location; positioning of a dry fiber preform into the desired location on the tool allows for low molding pressures of 10-50 psi. Lower molding pressures allow for use of low cost, additive fabrication of polymeric tooling. Polymeric tooling is suitable for rapid part prototyping and limited production. Dry fiber preforming approaches evaluated included use of commercial chopped strand mat, robotic chopper gun deposition, and continuous fiber preform augmentation using tailored fiber placement (TFP). Use of chopped strand mat does not require a robotic deposition method, however a cutting table is generally required and there is typically 20-30% scrap generation. While various fiber areal weights are available, the preform is not readily optimized for minimal fiber use or weight savings. In contrast, a robotic chopper gun approach allows for localized deposition where fiber is required to meet structural requirements. The robotic method is highly automated and minimizes fiber scrap, however the capital cost of the equipment and engineering labor for programming can result in higher preform cost compared to chopped strand mat in certain cases depending on preform complexity. Dry fiber preforming using the robotic chopper gun method allows for creation of three dimensional forms. This approach may be ideal for molding in-house, or if the preforms stack together densely to allow for efficient shipping. Applications evaluated for this program considered trade-off between fabrication of a fully 3D preform versus production of a flat preform which is designed to readily drape into the final desired shape. Such a preform design greatly simplifies robotic programming and requires no specialized tooling. The flat preforms are easily stacked and shipped to the final molding location. Flat preforms are much easier to augment with TFP continuous fiber to provide local reinforcement. The demonstration and evaluation of these preforming and tooling methods were completed on three component applications. The first application was a battery box cover for an electric vehicle which was highly three dimensional. The second demonstrator article was comprised of complex contours and was used to demonstrate the use of TFP and RTM-light molding process. The third demonstration article was the roof of an operator’s cab for large construction equipment. The roof is relatively flat however it is comprised of complex changes in thickness which clearly demonstrate the advantage of robotic chopper gun approach as compared to using numerous preform layers of chopped strand mat. The cost trades for the various preforming methods are summarized to help guide the reader as to preforming method considerations. Finally, these demonstrations all used glass fiber roving. A fourth, exploratory task was added to evaluate the ability to make preforms using Zoltek’s carbon fiber split tow roving. We were able to adapt the chopper gun to make flat preforms for laminate testing, but further development effort would be required to make suitable preforms.

36 MATERIALS SCIENCE↗

System Study: Residual Heat Removal 1998-2024

This report presents an unreliability evaluation of the residual heat removal (RHR) system in two modes of operation (low-pressure injection in response to a large loss-of-coolant accident and in response to post-trip shutdown cooling) at 93 U.S. commercial operating nuclear reactors. New Standardized Plant Analysis Risk (SPAR) models with the most recent SPAR parameter update results were used in this report. Demand, run hour, and failure data from 1998 to 2024 for selected components were obtained from the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS), formerly the INPO Consolidated Events Database. The unreliability results are trended for the most recent 10 year period and yearly estimates for system unreliability are provided for the entire active period. Statistically significant and highly statistically significant decreasing trends were observed for RHR low-pressure injection mode start-only unreliability and RHR low-pressure injection model 24-hour unreliability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗

Analytic fits to atom-in-jellium shear modulus predictions

Fits to AJ shear modulus calculations from ρ0 to the limit of the corresponding AJ EOS table (usually 1000 ρ0), adjusted to reproduce STP shear modulus where possible. Fitting minimized the fractional difference, so low-pressure points would be reproduced accurately despite the wide-ranging tabulation. AJ calculations were at T=0. Isochoric variation G(T) has been small for all cases spot-tested so far, and may have either sign. The fitted equation usually matches the AJ data to within a few percent. Between numerical noise in the AJ calculation and probably-physical structure not captured by the equation, the deviation could be up to 20% in some places in most models, and 30% in a few. AJ does not capture structural phase transitions or phases stabilized by directional bonds, where the shear modulus may vary by a greater amount. Where possible, STP ρ,G were used as parameters, and low pressure AJ points were de-weighted or removed if necessary. Otherwise, such as where AJ fails to capture solid phases with a significantly different shear modulus, the AJ data were fitted as far down in pressure as possible, G0 was also fitted if necessary, and ρ0 was also adjusted if needed to keep G0>0. The resulting models are not intended for use at low pressure, though some are probably adequate for practical purposes

36 MATERIALS SCIENCE↗

Measurement of Convective Heat Transfer Coefficients with Supercritical CO2 in Novel Additively Manufactured Helically Patterned Pin Fin Tubes Using the Wilson Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO2 flowing in pin-fin patterned pipes in the Heat Exchange and Experimental Testing (HEET) facility at the US DoE’s National Energy Technology Laboratory (NETL) in Morgantown, WV. The measurement procedures in the HEET rig were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over typical Reynolds number range in sCO2, the measured Nu and friction factors were within 7% of classical correlations for smooth tube flow.The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section. Pin length to diameter aspect ratios were 1.33 and 2, while the pin diameter to tube diameter ratio was 0.188 and 0.125. Tests were performed for ReD varying from 6.9×104 to 2.2×105 and at conditions equivalent to the low pressure outlet (8.69 MPa, 361 K) of the low temperature recuperator (LTR) in an indirect sCO2 power cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients.For the better performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.188), the local heat transfer coefficient increased by 112% relative to the Dittus-Boelter correlation at the LTR low pressure outlet. This corresponded to a 282% increase in the product of the heat transfer coefficient and the surface area. Large pressure drops across the test articles were observed.

Searle, Matthew↗

Decontamination of urban surfaces contaminated with radioactive materials and consequent onsite recycling of the waste water

Enhancing rapid remediation strategies is paramount for recovery after a large-scale nuclear contamination event in an urban environment. Some current strategies recommend use of readily available equipment, materials, and facilities to expedite recovery. For example, applying pressurized water to contaminated surfaces may effectively remove radioactive contamination. In this study, a commercial power washer removes soluble forms of 152 Eu 3 + , 85 Sr 2 + , and 137 Cs + contamination from common porous building materials, and computer simulations characterize the recycling of the resultant contaminated wash water. Pressure washing the porous building materials under spray conditions typical with do-it-yourself units improved decontamination factors (DFs) for 152 Eu compared to low-pressure application of tap water (majority of two-tailed t-test p-values < 0.1), but pressure did not improve DFs for 137 Cs or 85 Sr. For both pressurized and low-pressure applications, adding potassium ions (K + ) to promote ion exchange reactions produced significantly higher DFs for tested radionuclides on asphalt, brick, and concrete. The resultant contaminated wash water can be processed through self-prepared chemical filtration beds of clay and sand. Modeled in a prior study, the beds yielded linear trends (R 2 > 0.98) in sensitivity analyses between most bed configuration variables and bed performance variables, permitting flexible ad-hoc bed design. The experimental and simulation results led to estimates of the remediation rate and waste generated after cleaning 250 m 2 of cesium-contaminated concrete from the combined deployment of a power washer and two different mobile treatment beds. Furthermore, the first treatment bed was designed to reduce treatment time and processed 1900 L of wash solution in 70 min using 880 kg of clay/sand infill material. Designed to reduce the solid waste generated, the second bed processed the same solution volume in 1040 min (17 h) using 170 kg of clay/sand infill material. The results of this analysis warrant further investigation of power washing with recycled salt solution as an effective rapid decontamination method with manageable waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ab initio prediction of an order-disorder transition in Mg 2 Ge O 4 : Implication for the nature of super-Earth's mantles

Here we present an ab initio prediction of an order-disorder transition (ODT) from a I¯42d-type to Th 3 P 4 -type phase in the cation sublattices of Mg 2 GeO 4 , a post-post-perovskite phase. This uncommon type of prediction is achieved by carrying out a high-throughput sampling of atomic configurations in a 56-atom supercell followed by a Boltzmann ensemble statistics calculation. Mg 2 GeO 4 is a low-pressure analog of I¯42d-type Mg 2 SiO 4 , a predicted major planet-forming phase of super-Earths' mantles. Therefore, a similar ODT is anticipated in I¯42d-type Mg 2 SiO 4 as well, which should impact the internal structure and dynamics of these planets. Furthermore, the prediction of this Th 3 P 4 -type phase in Mg 2 GeO 4 further enhances the relationship between the crystal structures of Earth/planet-forming silicates and oxides at extreme pressures and those of rare-earth sesquisulfides at low pressures.

36 MATERIALS SCIENCE↗

Cooling plate assembly for plasma windows positioned in a beam accelerator system

A beam accelerator system operable to produce a medical isotope, including an ion accelerator that generates an ion beam; a low-pressure chamber; an anode adjacent and fluidly connected to the low-pressure chamber; a plasma window adjacent and fluidly connected to the anode; and a cathode housing adjacent and fluidly connected to the plasma window. The plasma window has a plurality of plates, each plate having an aperture that is aligned with an aperture in one or more adjacent plates to form a plasma channel. One or more plates in the plurality of plates includes a unitary plate having an aperture therein, and one or more cooling channels entering the unitary plate at a first side of the unitary plate and exiting the unitary plate at a second side of the unitary plate. The one or more cooling channels run through a thickness of the unitary plate.

Gribb, Tye↗

The charge reduction rate for multiply charged polymer ions via ion–ion recombination at atmospheric pressure

The charge reduction of multiply charged macromolecular ions via recombination with small ions in the gas phase is commonly employed to modulate the charge on macromolecules prior to mass spectrometric and mobility analyses. We employ a recently developed continuum-Molecular Dynamics (MD) calculation approach to determine the recombination rate coefficient of multiply charged (1 to 7 excess positive charged) polyethylene glycol ions (mass of 4600 Da) with smaller singly charged anions, modeled as NO 2 – ions. The continuum-MD approach accounts explicitly for the influence of the background gas on the recombination process, accounts explicitly for ion translational, vibrational, and rotational motion, and enables recombination rate coefficient calculation in nitrogen near atmospheric pressure, wherein neither low pressure nor high pressure recombination theories are strictly applicable. Continuum-MD simulations yield recombination rate coefficients near 3.9 × 10 –14 m 3 s –1 for singly charged ions, increasing to 3.0 × 10 –11 m 3 s –1 for the +7 ion. Pre-existing collision rate coefficient expressions for rigid ions are found to be within a factor of 2–5 of calculations for all charge states, but their use requires knowledge of an appropriate collision distance, which is not well-defined for flexible polymer ions. Continuum-MD-inferred rate coefficients are incorporated into a model of charge reduction, and the charge state distribution versus anion concentration determined with it is compared to charge reduction measurements made via atmospheric pressure differential mobility analysis. Furthermore, good agreement is observed between simulations and experiments; although results are extremely sensitive to the recombination rate coefficients, experimental results are bound by models utilizing rates within a factor of 2 (0.5–2.0×) of the continuum-MD rates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Electrolyzer Integrated Modular Nano-Array Monolithic Catalytic Reactors

The project objective is to design and develop a new class of electrolyzer-integrated, modular, nanostructure array (nano-array) monolithic catalytic reactors for high-flux, robust and selective methanol synthesis. The modular reactor features metal oxide nano-array grown on an additively-manufactured channeled honeycomb monolith made of ceramic or metals; this low-pressure carbon dioxide hydrogenation reactor is fueled with hydrogen from a low-temperature water electrolyzer having reduced platinum group metal (PGM) content. Through these nanostructured monolithic reactors, methanol synthesis is enabled with cost-effective hydrogenation reactions under low-temperature (<200°C) and low-pressure (<10 atm) conditions, thus significantly reducing the energy demand of the entire process.

08 HYDROGEN↗

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗