SEARCH · Engineering Papers

Results for “Constants”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2025

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗

DOE OSTI2021

Development of an anti-clogging perforated plate atomizer for a zero liquid discharge humidification-dehumidification desalination system

An anti-clogging perforated plate atomizer is designed for high temperature and salinity applications in a novel solar-powered zero liquid discharge humidification-dehumidification desalination system (US Patent Application US62882953). Herein this paper presents a detailed discussion on its design, and operation. Experiments are performed on seven atomizers to study the effect of design, and operating parameters on spray cone angle, and average droplet diameter. Spray cone angle remains constant with change in air mass flux, and increases with increasing orifice diameter, and manifold diameter. It also increases initially with water mass flow rate, and manifold height, before attaining a constant value. Average droplet diameter increases with increasing water mass flow rate, and decreasing air mass flux, orifice diameter, and manifold diameter. Further, anti-clogging performance of the atomizer is tested with hypersaline water of 100,000 ppm total dissolved solid (NaCl) at elevated temperatures such as air at 175 °C and saline water at 45 °C. Results show no clogging for 2–13 h of operation. Near-complete suppression of atomizer clogging makes our new additively-manufactured perforated plate atomizer an ideal fit for high salinity, and zero liquid discharge humidification-dehumidification desalination systems. Additionally, its open-surface design allows additional surface modifications to further reduce clogging, and enhance self-cleaning characteristics.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Kinetic study of hydrogen transport in graphite under molten fluoride salt environment

Here, in this work, a kinetic model, which describes the reactions occurring during the hydrogen charging process on a graphite surface, was optimized on Electrochemical Impedance Spectroscopic (EIS) to study the entry of hydrogen into graphite in high temperature molten FLiNaK salt. The results show that the rate constants of k 1 increases with decreasing applied potential (becoming more negative), and k 2 and k 3 are independent of the applied potential. At the same charging potential, all the three rate constants are independent of the water concentration in the FLiNaK salt. The surface coverage of absorbed hydrogen increases with decreasing the charging potential and increasing the moisture content in the salt. The adsorption efficiency of hydrogen increases with increasingly charging potential and decreases with increasing moisture content of the FLiNaK salt. These results demonstrate that EIS is a promising technique for measuring the kinetic parameters for hydrogen transport in graphite under molten fluoride salt environment.

EIS↗

DOE OSTI2021

Decoupling of particles and dissolved iron downstream of Greenlandic glacier outflows

Glaciers can be a significant and locally dominant source of iron (Fe), a biologically essential micronutrient, in high latitude coastal seas. The vast majority of this glacial Fe delivery is associated with particles, yet the speciation of the solid-phase Fe and specifically the relationships that govern exchange between particulate and dissolved Fe phases in these environments are poorly described. Here, in this work, we performed measurements of in situ dissolved Fe (dFe) along meltwater and particle plumes in three transects around Disko Bay and Ameralik Fjord (West Greenland). Measurements of dFe were combined with Fe K-edge X-ray absorption spectroscopy analysis of ~40 suspended sediment samples obtained from the same transects and from select depth profiles down to 300 m. We observed relatively constant dFe levels (4 to 10 nM for nearly all dFe measurements) across fjords with widely varying particulate Fe(II) contents (from 20 to 90% Fe(II)), indicating that dFe concentrations had little dependence on the oxidation state of Fe in the suspended sediment. Particulate Fe data were grouped by underlying bedrock geology, with suspended sediment consisting of 80-90% biotite-like Fe(II) in fjords with Precambrian shield geology and poorly-ordered Fe(III) particles (<20-30% Fe(II)) in one fjord with suspended sediments derived from tertiary basalts. Our characterization data indicated no significant change in the average Fe oxidation state and bonding environment of particles along the fjord transects, implying that Fe(II) in biotite-like coordination is not a readily labile Fe form on this spatial scale. Finally, our results suggest that dFe in these glacially-modified coastal waters is buffered at a relatively constant low nM concentration due to factors other than particle Fe mineralogy and that glacier-derived Fe phases are relatively inert on this spatial scale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗

DOE OSTI2021

Iron, magnesium, and titanium isotopic fractionations between garnet, ilmenite, f ayalite, biotite, and tourmaline: Results from NRIXS, ab initio , and study of mineral separates from the Moosilauke metapelite

Interpreting isotopic signatures documented in natural rocks requires knowledge of equilibrium isotopic fractionation factors. Here, we determine equilibrium Fe isotope fractionation factors between several common rock-forming minerals using a comparative approach involving three independent methods: (i) isotopic analyses of natural minerals from a metapelite from Mt. Moosilauke, New Hampshire, for which equilibration temperature and pressure are well constrained to be near the aluminosilicate triple point (T ≃ 500 °C, P ≃ 4 kbar), (ii) Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements of Fe force constants of minerals, and (iii) Density Functional Theory (DFT) ab initio calculations of Fe force constants of minerals. The minerals studied for Fe isotopes include, in increasing order of their β-factors: garnet < ilmenite ≈ fayalite < biotite < tourmaline < muscovite ≈ plagioclase. Some of this ordering is affected by the presence of Fe 3+ in the minerals, which tends to form stiffer bonds and be associated with heavy Fe isotope enrichments relative to Fe 2+ . We are, however, able to assess the magnitude of the effect of the ratio Fe 3+ /ΣFe on equilibrium fractionation factors, notably on the ilmenite-hematite solid solution. Equilibrium Fe isotopic fractionation factors between garnet, ilmenite, biotite, tourmaline and fayalite are determined. We also report Mg and Ti isotopic compositions of selected Moosilauke minerals that allow us to better constrain the equilibrium fractionation factors for garnet-biotite-tourmaline (Mg isotopes) and biotite-ilmenite (Ti isotopes). We show how the newly determined equilibrium fractionation factors can be used to address diverse problems in Earth and planetary sciences, notably (i) Fe and Mg isotopic fractionation during anatexis, (ii) Fe isotopic fractionation in lunar ilmenite, and (iii) Ti isotopic fractionation during fluvial transport of minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗

DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗

DOE OSTI2021

Droplet Evaporation on Hot Micro-Structured Superhydrophobic Surfaces: Analysis of Evaporation from Droplet Cap and Base Surfaces

In this study, evaporation of sessile water droplets on hot micro-structured superhydrophobic surfaces is experimentally and theoretically investigated. Water droplets of 4 µL are placed on micro-pillared silicon substrates with the substrate temperature heated up to 120°C. A comprehensive thermal circuit model is developed to analyze the effects of substrate roughness and substrate temperature on the sessile droplet evaporation. For the first time, two components of heat and mass transfer, i.e., one from the droplet cap surface and the other from the droplet base surface, during droplet evaporation are distinguished and systematically studied. As such, the evaporation heat transfer rates from both the droplet cap surface and the interstitial liquid-vapor interface between micropillars at the droplet base are calculated in various conditions. For droplet evaporation on the heated substrates in the range of 40°C – 80°C, the predicted droplet cap temperature matches well with the experimental results. Furthermore, during the constant contact radius mode of droplet evaporation, the decrease of evaporation rate from the droplet base contributes most to the continuously decreasing overall evaporation heat transfer rate, whereas the decrease of evaporation rate from the droplet cap surface is dominant in the constant contact angle mode. The influence of internal fluid flow is considered for droplet evaporation on substrates heated above 100°C, and an effective thermal conductivity is adopted as a correction factor to account for the effect of convection heat transfer inside the droplet. Temperature differences between the droplet base and the substrate base are estimated to be about 2°C, 5°C, 8°C, 13°C and 18°C for droplet evaporation on substrates heated at 40°C, 60°C, 80°C, 100°C, and 120°C, respectively, elucidating the delayed or depressed boiling of water droplets on a heated rough surface due to evaporative cooling.

42 ENGINEERING↗

DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

DOE OSTI2020

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A thermodynamic model of integrated liquid-to-liquid thermoelectric heat pump systems

Thermoelectric (TE) heat pumps (TEHPs) are advantageous for heating and cooling in various applications because of their modularity and simple design. A TEHP system includes the TE modules with p- and n-type materials bonded to substrates, plus heat exchangers, thermal interfaces to the heat exchangers, and heat transfer fluids. Although modeling an individual TE module has been extensively studied, limited studies have reported performance at the larger system-level. Furthermore, no prior study has addressed the impact of temperature-dependent TE material properties (e.g., electric resistivity, thermal conductivity, and Seebeck coefficient) on overall heat-pump-system-level performance. Here, this work presents a mathematical model for TEHP system performance based on Goldsmid's approach for TE material performance, “effective” TE material properties, Gnielinski's correlation for convective heat transfer, and thermal balance theory for a heat exchange network. This combined approach provides an accurate model of the liquid-to-liquid TEHP system. Three different approaches—one empirical, one based on the manufacturer's specifications, and one drawn from the literature—were then used to determine values for TE material properties. The first two methods treated properties as constants, while the last approach treated properties as surface-temperature-based functions. Finally, experimental TEHP data was used to validate the models, all with relative absolute deviations of approximately 10% when predicting heating capacity and 10%–25% when forecasting cooling capacity up to a 30 K surface temperature lift. The results demonstrated that, at the TEHP system level, the TE material properties could be treated as constants, avoiding solver iterations and reducing the performance uncertainty by up to 95%.

42 ENGINEERING↗

DOE OSTI2026

A single-zone zero-dimensional study of HCCI combustion of methanol dehydration products to enable ignition of direct-injected methanol

Methanol is an alternative fuel gaining traction in the maritime sector. Its direct adoption, however, is accompanied by a unique set of technical challenges, such as low cetane number and high latent heat of vaporization. An approach to overcome these challenges is being developed at the US Department of Energy’s Oak Ridge National Laboratory, where onboard generation of dimethyl ether (DME) via catalytic dehydration of methanol can be used to assist in the mixing controlled combustion of direct-injected (DI) methanol. The generated mixture from this dehydration process can be premixed with intake air to condition the cylinder via. homogeneous charge compression ignition (HCCI) for subsequent DI methanol. In this preliminary work, various catalyst or reactor conversion efficiencies were simulated (using bottles) at constant DME and water flow at low load on a single-cylinder marine-variant of a CAT® C18 18 L engine with a 145 mm bore. To substantiate the experimental findings, a zero-dimensional engine model was developed in Cantera using a DME mechanism with 79 species and 658 reactions. Results presented include experimental and simulation heat release rate comparisons, species evolution information, and constant volume ignition delay (ID) for DI methanol with and without background species from HCCI of the premixed products from different reactor efficiencies. The results suggest that thermal effects dominate the DI methanol ignition process, and this work provides a chemical kinetic foundation or guideline for developing future control schemes.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗

DOE OSTI2021

The Delplot kinetic method applied to systems with adsorbates: Hydrodeoxygenation of benzofuran on a bimetallic CoPd phosphide catalyst supported on KUSY

The Delplot method provides a means of analyzing conversion and selectivity data to derive a reaction sequence and has been applied so far to gas-phase or liquid-phase species. This study analyses the applicability of the method for catalytic reactions and considers for the first time adsorbed intermediates. The method is applied to a contact time study of the hydrodeoxygenation (HDO) of benzofuran at 0.5 MPa and 350 °C over a catalyst consisting of bimetallic CoPd phosphide supported on a potassium ion-exchanged ultra-stable Y (KUSY) zeolite. Both simple and detailed reaction networks were derived. The simple networks considered only gas-phase species and were modeled by a first-order sequence of steps, whereas the detailed network took into account adsorbed intermediates and was simulated using a rake mechanism. The networks were compared using F-statistics, which accounted for the differences in the number of fitting parameters. The detailed network gave a better fit to the experimental data because it represented a more realistic description of the transformation. A suggested sequence from the Delplot method was consistent with the simple network but not with the detailed network. The lack of applicability of the Delplot analysis to the network with adsorbates was linked to overly large equilibrium constants, a determination supported by Delplot fitting for a model sequence. Further, this study indicated the inapplicability of the Delplot method in determining reaction sequences that involve adsorbed species with large equilibrium constants for formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Identification and quantification of distinct active sites in Hf-Beta zeolites for transfer hydrogenation catalysis

Despite the significant progress made in characterizing different framework heteroatom sites that exist in Lewis acidic zeolites with probe molecule adsorption and spectroscopy, methods to reliably quantify site counts remain indefinite and have been primarily limited to Sn and Ti Lewis acid sites. Here, methods to quantify framework Lewis acidic Hf 4+ sites in zeolite Beta (Hf-Beta) with two Lewis base titrants (pyridine, deuterated acetonitrile) were developed using infrared (IR) spectroscopy. Lewis acid site counts for Hf-Beta zeolites were validated by measuring integrated molar extinction coefficients (IMECs; ε, cm μmol –1 ) on Sn-Beta zeolites using identical Lewis base titrants to benchmark site counts with established literature procedures to quantify Lewis acid sites in Beta zeolites from IR spectra. Highlighting the importance of benchmarking active site counts against well-established experimental protocols, IMECs of CD 3 CN bound to open (ε(Sn; 2316 cm –1 ): 1.80 ± 0.25) and closed (ε(Sn; 2308 cm –1 ): 3.76 ± 0.33) Sn sites were ~1.8x larger than those previously reported while total Lewis acid site counts agreed with those measured by pyridine (ε(Sn; 1451 cm –1 ): 1.58 ± 0.16) on six different Sn-Beta zeolites. IMECs measured for IR peaks reflecting pyridine bound to Lewis acidic Hf sites (ε(Hf; 1448 cm –1 ): 1.54 ± 0.21) and CD 3 CN bound to open (ε(Hf; 2313 cm –1 ): 2.40 ± 0.22) and closed (ε(Hf; 2307 cm –1 ): 3.55 ± 0.41) Hf sites, gave similar counts for the total number of Lewis acidic sites across six Hf-Beta zeolites (Si/Hf = 100–413). Consistent with previous reports with Sn-Beta catalysts where open Sn sites are responsible for catalytic turnover, apparent first and zero-order MPVO rate constants (0.01–1 M cyclohexanone in 2-butanol; per total Hf, 373 K) correlated with the total number of open Hf sites, per total Hf, but not with the total number of closed Hf sites or total Lewis acid site counts. Measured initial MPVO rates (0.1 M cyclohexanone in 2-butanol, per open Hf, 373 K) were ~25x higher on hydrophobic Hf-Beta-F than on hydrophilic Hf-Beta-OH zeolites. Overall, the apparent first-order MPVO rate constants (2-butanol solvent, per open Hf, 373 K) were ~6x higher on Hf-Beta-F than on Hf-Beta-OH zeolites. The characterization methods reported here enable normalization of MPVO turnover rates on Sn- and Hf-Beta zeolites by their number of open sites. Finally, this enables performing quantitative rate comparisons across Lewis acid zeolites of varying active site identity, solvation, and pore topology used in liquid-phase catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Impact of inland boundary conditions on seawater intrusion modeling in stratified coastal aquifers under sea-level rise

This study systematically examines the impact of inland boundary conditions on seawater intrusion modeling in general stratified coastal aquifers under sea-level rise (SLR). We derived analytical solutions for seawater wedge toe positions in stratified aquifers with general head boundary (GHB) conditions and compiled solutions for constant head boundary (CHB), constant flux boundary (CFB), and GHB conditions in both homogeneous and stratified aquifers. Additionally, we introduced the concept of equivalent boundary conditions and developed equations to ensure consistent steady-state toe positions across different boundary types. Unlike previous studies focusing on specific layering configurations, we conducted a comprehensive sensitivity analysis of seawater intrusion in general stratified aquifers using transmissivity centroid elevation (TCE). Our results show that SLR generally exacerbates seawater intrusion, with unconfined aquifers exhibiting higher sensitivity across all boundary conditions. Among boundary types, CHB conditions were the most sensitive to SLR, CFB the least, and GHB showed an intermediate response depending on its distance from the model boundary. Furthermore, seawater intrusion increases with the TCE of the stratified aquifer under SLR, suggesting that high-permeability or preferential flow layers near the surface amplify the impact of SLR. These findings emphasize the importance of selecting appropriate boundary conditions when evaluating SLR-driven seawater intrusion in stratified coastal aquifers.

Boundary condition↗