SEARCH · Engineering Papers

Results for “experimental conditions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Spectroscopy of the fractal Hofstadter energy spectrum

Hofstadter’s butterfly, the predicted energy spectrum for non-interacting electrons confined to a two-dimensional lattice in a magnetic field, is one of the most remarkable fractal structures in nature. At rational ratios of magnetic flux quanta per lattice unit cell, this spectrum shows self-similar distributions of energy levels that reflect its recursive construction. For most materials, Hofstadter’s butterfly is predicted under experimental conditions that are unachievable using laboratory-scale magnetic fields. More recently, electrical transport studies have provided evidence for Hofstadter’s butterfly in materials engineered to have artificially large lattice constants, such as those with moiré superlattices. Yet, so far, direct spectroscopy of the fractal energy spectrum predicted by Hofstadter nearly 50 years ago has remained out of reach. Here, in this work, we use high-resolution scanning tunnelling microscopy/spectroscopy (STM/STS) to investigate the flat electronic bands in twisted bilayer graphene (TBG) near the predicted second magic angle, an ideal setting for spectroscopic studies of Hofstadter’s spectrum. Our study shows the fractionalization of flat moiré bands into discrete Hofstadter subbands and discerns experimental signatures of self-similarity of this spectrum. Moreover, our measurements uncover a spectrum that evolves dynamically with electron density, showing phenomena beyond that of Hofstadter’s original model owing to the combined effects of strong correlations, Coulomb interactions and the quantum degeneracy of electrons in TBG.

electronic properties and materials↗

NASA NTRS2001

Determination of a Two-Phase Structure of Nanocrystals: GaN and SiC

The properties of nano-crystalline materials are critically dependent on the structure of the constituent grains. Experimental conditions necessary to perform structural analysis of nanocrystalline materials as a two-phase core-surface shell system are discussed. It is shown, that a standard X-ray diffraction measurements and analysis are insufficient and may lead to incorrect conclusions as to the real structure of the materials. A new method of evaluation of powder diffraction data based on the analysis of the shift of the Bragg reflections from their perfect-lattice positions was developed. "Apparent lattice parameters" quantity, alp, was introduced and calculated from the actual positions of each individual Bragg reflection. The alp values plotted versus diffraction vector (Q) show characteristic features that are used for evaluation of the experimental results. The study was based on modeling of nano-grains and simulations of theoretical intensity profiles using the Debye functions. The method was applied to the analysis of synchrotron X-ray diffraction data of GaN and SiC nanocrystals. A presence of strained surface shell and a considerable internal pressure (GaN) in the nanoparticles was concluded.

Palosz, W.↗

DOE OSTI2022

When physics-informed data analytics outperforms black-box machine learning: A case study in thickness control for additive manufacturing

Aerosol jet printing (AJP) has emerged as a promising noncontact additive manufacturing method for high-resolution printing for a wide range of material systems. A key challenge limiting the broader adoption of AJP in the material science community is the lack of methods to precisely control thickness. Herein, we develop a model-based design of experiment (MBDoE) framework that integrates physics-informed models, nonlinear regression, and information criteria to postulate, select and calibrate the best model to describe and optimize the AJP manufacturing process. Starting with already available data from system commissioning (e.g., prior single variable sensitivity analysis), four candidate physics-informed models are postulated and trained. MBDoE identifies a single additional optimal experiment to validate these predictive models with quantified uncertainties, which are then used to determine the best experimental conditions to control printed film thickness. As a comparative benchmark, the analysis is repeated using the same dataset with nonparametric Gaussian process regression (GPR) model that does not incorporate physical information. Using MBDoE principles, we find that only five experiments are necessary to calibrate the nonlinear physics-informed parametric model, and with said limited data, this model outperforms the black-box machine learning GPR model. This key result underscores an emerging trend in the data science community: incorporating physical information into predictive models often drastically reduces the data requirements. Leveraging MBDoE further increased the data efficiency. By design, the proposed data science framework is general in nature and can be easily extended to other experimental and additive manufacturing systems beyond AJP.

Aerosol jet printing↗

DOE OSTI2026

Sliding friction and wear behavior of nuclear graphite in high temperature inert environment: Influence of contact load, speed and temperature

Repeated dynamic interactions of graphitic components in pebble-bed gas-cooled nuclear reactors can cause abrasive wear-induced pebble surface damage, generate hazardous fine graphite debris, and alter fuel circulation dynamics due to changes in friction behavior. Comprehensive tribological characterization of nuclear graphitic materials in conditions relevant to reactor operation is needed to assess reactor long-term safety and performance. This work reports sliding friction and wear behavior of self-mated nuclear graphite ET-10 at various elevated temperatures (650 °C and 750 °C), sliding speeds (1 and 10 mm/s) and contact loads (20 and 40 N) in a controlled argon environment. The results revealed nonmonotonic frictional behavior with a higher running-in coefficient of friction (COF) followed by a lower steady-state COF, as a result of transition from two-body abrasion to three-body abrasion along with formation of a tribofilm. A key finding of this work is the sensitivity of the running-in COF to experimental conditions; maximum running-in values were lower at either elevated temperature (0.52–0.54) or reduced sliding speed (0.51–0.54). Conversely, the steady-state COF remained invariant at approximately 0.3 across all tested parameters. Transmission electron microscopy revealed a 0.5–2.0 μm thick nanocrystalline tribofilm that was thought to be formed by the compaction of the graphitic wear debris on the contact surface during the sliding process. The nanocrystalline nature of the tribofilm was further confirmed by Raman spectroscopy. As a result, the combination of tribological testing and morphological characterization provided a mechanistic understanding of the frictional behavior of nuclear graphite upon sliding.

Friction↗

DOE OSTI2021

Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery

One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.

03 NATURAL GAS↗

NASA NTRS1995

Spray combustion at normal and reduced gravity in counterflow and co-flow configurations

Liquid fuel dispersion in practical systems is typically achieved by spraying the fuel into a polydisperse distribution of droplets evaporating and burning in a turbulent gaseous environment In view of the nearly insurmountable difficulties of this two-phase flow, a systematic study of spray evaporation and burning in configurations of gradually increasing levels of complexity, starting from laminar sprays to fully turbulent ones, would be useful. A few years ago we proposed to use an electrostatic spray of charged droplets for this type of combustion experiments under well-defined conditions. In the simplest configuration, a liquid is fed into a small metal tube maintained at several kilovolts relative to a ground electrode few centimeters away. Under the action of the electric field, the liquid meniscus at the outlet of the capillary takes a conical shape, with a thin jet emerging from the cone tip (cone-jet mode). This jet breaks up farther downstream into a spray of charged droplets - the so-called ElectroSpray (ES). Several advantages distinguish the electrospray from alternative atomization techniques: (1) it can produce quasi-monodisperse droplets over a phenomenal size range; (2) the atomization, that is strictly electrostatic, is decoupled from gas flow processes, which provides some flexibility in the selection and control of the experimental conditions; (3) the Coulombic repulsion of homopolarly charged droplets induces spray self-dispersion and prevents droplet coalescence; (4) the ES provides the opportunity of studying regimes of slip between droplets and host gas without compromising the control of the spray properties; and (5) the compactness and potential controllability of this spray generation system makes it appealing for studies in reduced-gravity environments aimed at isolating the spray behavior from natural convection complications. With these premises, in March 1991 we initiated a series of experiments under NASA sponsorship (NAG3-1259 and 1688) in which the ES was used as a research tool to examine spray combustion in counter-flow and co-flow spray diffusion flames, as summarized below. The ultimate objective of this investigation is to examine the formation and burning of sprays of liquid fuels, at both normal and reduced gravity, first in laminar regimes and then in turbulent ones.

Gomez, Alessandro↗

NASA NTRS2007

Risk Reduction and Measures of Injury for EVA Associated Upper Extremity Medical Issues: Extended Vent Tube Study

During Neutral Buoyancy Lab (NBL) training sessions, a large amount of moisture accumulates in the EVA gloves. The glove design restricts the extension of the EVA suit s ventilation/cooling system to the hand. Subungual redness and fingernail pain develops for many astronauts following their NBL training sessions with subsequent oncholysis occurring over succeeding weeks. Various attempts have been made to reduce or avoid this problem. The causal role of moisture has yet to be defined. Methods: To determine the contribution that moisture plays in the injury to the fingers and fingernails during EVA training operations in NBL, the current Extravehicular Mobility Unit (EMU), with a Portable Life Support System (PLSS) was configured with a ventilation tube that extended down a single arm of the crewmember during the test and compared with the unventilated contralateral arm; with the ventilated hand serving as the experimental condition (E) and the opposite arm as the control (C). A cross-over design was used with opposite handedness for the vent tube on a subsequent NBL training run. Moisture content measures were conducted at six points on each hand with three types of moisture meters. A questionnaire was administered to determine subjective thermal hand discomfort, skin moisture perception, and hand and nail discomfort. Photographs and video were recorded. Measures were applied to six astronauts pre- and post-run in the NBL. Results: The consistent trends in relative hydration ratios at the dorsum, from 3.34 for C to 2.11 for E, and first ring finger joint locations, from 2.46 for C to 1.96 for E, indicated the extended vent tube promoted skin drying. The experimental treatment appeared to be more effective on the left hand versus the right hand, implying an interaction with hand anthropometry and glove fit. Video analyses differentiated fine and gross motor training tasks during runs and will be discussed. Conclusions: This potential countermeasure was effective in reducing the risks of hand and nail discomfort symptoms from moderate to low in two of six subjects. Improved design in the ventilation pattern of such a countermeasure is expected to improve the countermeasure s efficiency.

Jones, Jeffrey A.↗

NASA NTRS2006

Uptake and Dissolution of Gaseous Ethanol in Sulfuric Acid

The solubility of gas-phase ethanol (ethyl alcohol, CH3CH2OH, EtOH) in aqueous sulfuric acid solutions was measured in a Knudsen cell reactor over ranges of temperature (209-237 K) and acid composition (39-76 wt % H2SO4). Ethanol is very soluble under these conditions: effective Henry's law coefficients, H*, range from 4 x 10(exp 4) M/atm in the 227 K, 39 wt % acid to greater than 10(exp 7) M/atm in the 76 wt % acid. In 76 wt % sulfuric acid, ethanol solubility exceeds that which can be precisely determined using the Knudsen cell technique but falls in the range of 10(exp 7)-10(exp 10) M/atm. The equilibrium concentration of ethanol in upper tropospheric/lower stratospheric (UT/LS) sulfate particles is calculated from these measurements and compared to other small oxygenated organic compounds. Even if ethanol is a minor component in the gas phase, it may be a major constituent of the organic fraction in the particle phase. No evidence for the formation of ethyl hydrogen sulfate was found under our experimental conditions. While the protonation of ethanol does augment solubility at higher acidity, the primary reason H* increases with acidity is an increase in the solubility of molecular (i.e., neutral) ethanol.

Michelsen, Rebecca R.↗

DOE OSTI2024

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗

DOE OSTI2021

Shock tube investigation of high-temperature, extremely-rich oxidation of several co-optima biofuels for spark-ignition engines

To reduce the reliance on fossil fuels in the transportation sector and increase combustion efficiency, the Co-Optima initiative from US Department of Energy identified the top 10 biofuels for downsized, boosted, spark-ignition engines. Most of these biofuels have detailed reaction mechanisms available in literature developed based on studies at temperatures lower than 1700 K and an equivalence ratio of less than five. As such, the performance of these detailed mechanisms at high temperature and extremely rich conditions are unknown. It is important to validate kinetic mechanisms at these conditions because they are conducive to soot formation. Prediction of soot by chemical kinetic models relies on the prediction of underlying benchmark species like carbon monoxide and ethylene. In this work, we conduct high temperature (1700–2050 K) and high equivalence ratio(Φ=8.6) oxidation of these biofuels, namely 2,4,4-trimethyl-1-pentene (α-diisobutylene), ethanol, cyclopentanone, methyl acetate, and 2-methylfuran, blended in ethylene behind reflected shock waves at 4–4.7 atm pressure. Carbon monoxide and ethylene time histories are measured simultaneously using a continuous feedback quantum cascade laser near 4.9 µm and a tunable CO 2 gas laser at 10.532 nm, respectively. Results show that ethanol blend forms more carbon monoxide than other biofuels and consumes ethylene faster than the biofuel blends in the temperature range considered. Overall, the performance of different mechanisms in literature are evaluated against the experimental results. The novel reaction mechanism ‘the Co-Optima model’, which includes the sub mechanisms for all the biofuels in this study, was found to be the best mechanism for the experimental conditions studied.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Testing Instrument Extremes in the TREAT Facility

TREAT is a transient power-shaping reactor whose primary mission is to enhance safety performance by testing nuclear fuels and materials under thermodynamic and neutronic conditions ranging from off-normal to extreme. Transient testing of nuclear fuels is analogous to car crash testing, in that the dramatic changes commonly seen when comparing the test specimens' initial and final state points necessitate in-situ measuring in order to interpret and understand the evolution of the experiment. TREAT’s experiment design strategy utilizes a highly reconfigurable and accessible reactor core that affords flexibility when installing experiment devices, thanks to the lack of a primary coolant boundary and the inclusion of multiple access points in the reactor bioshield. Within the experiment devices, the test specimens (which are integrated with their test environments), desired instrumentation diagnostics, and safety considerations are all contained in a single engineered package. The experiments encompass various types of environments (e.g., water, sodium, and gas at temperature and pressure), as well as a range of control options (e.g., static or flowing coolants). The experimental conditions are uniquely representative of advanced reactors, providing challenges to sensor performance but also affording a unique opportunity for development and qualification of sensors suited to these environments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2003

The Deflection Plate Analyzer: A Technique for Space Plasma Measurements Under Highly Disturbed Conditions

A technique has been developed to measure the characteristics of space plasmas under highly disturbed conditions; e.g., non-Maxwellian plasmas with strong drifting populations and plasmas contaminated by spacecraft outgassing. The present method is an extension of the capabilities of the Differential Ion Flux Probe (DIFP) to include a mass measurement that does not include either high voltage or contamination sensitive devices such as channeltron electron multipliers or microchannel plates. This reduces the complexity and expense of instrument fabrication, testing, and integration of flight hardware as compared to classical mass analyzers. The new instrument design is called the Deflection Plate Analyzer (DPA) and can deconvolve multiple ion streams and analyze each stream for ion flux intensity (density), velocity (including direction of motion), mass, and temperature (or energy distribution). The basic functionality of the DPA is discussed. The performance characteristics of a flight instrument as built for an electrodynamic tether mission, the Propulsive Small Expendable Deployer System (ProSEDS), and the instrument s role in measuring key experimental conditions are also discussed.

Wright, Kenneth H., Jr.↗

NASA NTRS2004

The Deflection Plate Analyzer: A Technique for Space Plasma Measurements Under Highly Disturbed Conditions

Wright, Kenneth H., Jr.↗

DOE OSTI2023

In-situ thermophysical measurement of flowing molten chloride salt using modulated photothermal radiometry

Molten salts are leading candidates for high-temperature heat transfer fluids (HTFs) for thermal energy storage and conversion systems in concentrated solar power (CSP) and nuclear energy power plants. The ability to probe molten salt thermal transport properties in both stationary and flowing status is important for the evaluation of their heat transfer performance under realistic operational conditions, including the temperature range and potential degradation due to corrosion and contamination. However, accurate thermal transport properties are usually challenging to obtain even for stagnant molten salts due to different sources of errors from convection, radiation, and corrosion, let alone at flowing status. To the best of authors’ knowledge, there is no available in-situ technique for measuring flowing molten salt thermal conductivity. Here, we report the first in-situ flowing molten salt thermal conductivity measurement using modulated photothermal radiometry (MPR). We could successfully perform the first in-situ thermal conductivity measurement of flowing molten NaCl-KCl-MgCl 2 in the typical operating temperature (520 and 580 °C) with flow velocities ranging from around 0.3 to 1.0 m s -1 . The relative change of the molten salt thermal conductivity was measured. Gnielinski’s correlation was also used to estimate the heat transfer coefficient h of the flowing NaCl-KCl-MgCl 2 in the given experimental condition. Furthermore, the work showed the potential of the MPR technique serving as an in-situ diagnostics tool to evaluate the heat transfer performance of flowing molten salts and other high-temperature HTFs.

14 SOLAR ENERGY↗

DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Throughput Optimization of Molybdenum Carbide Nanoparticle Catalysts in a Continuous Flow Reactor Using Design of Experiments

Transition metal carbides (TMCs) have attracted significant attention because of their applications toward a wide range of catalytic transformations. However, the practicality of their synthesis is still limited because of the harsh conditions in which most TMCs are prepared. Recently, a solution-phase synthesis of phase-pure a-MoC1-x nanoparticles was presented. While this synthetic route yielded nanoparticles with exceptional catalytic performance, the reaction parameter space was not explored, and catalyst throughput was not optimized for scale-up. Continuous flow platforms coupled with statistical design of experiments (DoE) can provide a powerful method for understanding the reaction parameter space for optimizations. Here, we demonstrate the use of statistical DoE in tandem with response surface methodology for a parametric screening analysis to optimize the throughput of a MoC1-x nanoparticle synthesis utilizing a millifluidic flow reactor. A full factorial design was implemented to evaluate four input variables (reaction temperature, flow rate, solvent fraction of oleylamine, and precursor concentration) that carry statistically significant effects on three responses (throughput, residence time, and isolated yield). A Doehlert matrix was implemented to investigate each significant variable at a higher number of levels to optimize throughput. Our results give a nonintuitive set of experimental conditions that resulted in an optimized throughput of 2.2 g h-1. This translates to a 50-fold increase in throughput compared to the previously reported batch method. The catalytic performance of the MoC1-x nanoparticles produced under optimized throughput was demonstrated in the CO2 hydrogenation reaction. This DoE screening analysis and throughput optimization of MoC1-x synthesis open the door to an increased feasibility for scale-up.

design of experiments↗