SEARCH · Engineering Papers

Results for “Homogenization theory”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Pd–Methyl Bond Energy─Property Correlations, Noncorrelations, Machine Learning Models, and Application to Polymerization Catalysis

Metal–carbon bonds are a key intermediate in a variety of homogeneous organometallic transformations and often determine the critical thermodynamics and kinetics of catalytic processes. Surprisingly, the influence of different ligands on metal–carbon bond strengths has been largely overlooked. Here, in this study, we evaluated nearly 700 experimental Pd–methyl complexes by calculating their bond dissociation energies using density functional theory (DFT) and compared these bond strengths to several fundamental molecular properties, and this revealed several surprising correlations and noncorrelations. Most surprising was that several fundamental properties, such as the bond length, bond force constant, and bond electron density, have no correlation with bond strength, despite these correlations often holding for main-group compounds. We were indeed able to identify key ligand-dependent chemical features/descriptors that provided a highly accurate machine learning model and provided insight into the general factors that control the Pd–carbon bond strength, such as radical delocalization and nucleophilicity. Insights gained from the Pd–Me bond energy analysis were then applied to CO migratory insertion steps that are part of copolymerization reactions.

binding energy↗

DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗

DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗

DOE OSTI2020

Theoretical Correlation of Elemental Distribution of Nd and Pr in Ce-Fe-B Microstructure With Hard Magnetic Properties

Relatively resource-rich but property inferior-Ce-Fe-B magnet can be improved by partial replacement of Ce by Nd and/or Pr. In addition to the amount of Nd/Pr, their distribution profile in microstructure plays an important role. From our first principles density functional theory (DFT) calculation, the substitution energy of Ce by Pr/Nd is negative in Ce 2 Fe 14 B (2:14:1) while that for laves phase, CeFe 2 is positive, implying that Nd/Pr stabilize 2:14:1 and suppress the formation of the CeFe 2 phase. Further, micromagnetic simulation indicates that homogenized distribution of Nd/Pr improves squareness of demagnetization curve, while core (Ce-rich)-shell (Nd/Pr-rich) 2:14:1 grain structure enhances coercivity. Magnetic properties of Ce-Fe-B can be optimized by manipulating distribution profile of chemical element in microstructure based on their subtle difference in thermodynamic property, which is an effective pathway to design optimized chemical composition and processing route for high-performance magnet.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Effects of surface species and homogeneous reactions on rates and selectivity in ethane oxidation on oxide catalysts

Selective alkane oxidations on metal oxide catalysts involve complex mechanisms with multiple reactions in series and parallel, different types of reduced and oxidized surface species, and potential contributions from gas-phase reactions. We report kinetics and thermodynamics of elementary steps involved in C 2 H 6 -O 2 reactions on SiO 2 -supported small vanadium oxide domains are determined using density functional theory. These surface reactions together with gas-phase mechanisms are incorporated in kinetic simulations to determine how surface and gaseous reactions interact and contribute to rates and selectivity. The results show that gas-phase reactions within pore volumes in contact with the catalyst contribute significantly to C 2 H 6 activation rates, even at conditions where gas-phase reactions in empty volumes without catalyst are negligible. The majority of C 2 H 6 activations occur on the surface, via H abstraction by vanadium oxo species present at terminal lattice oxygens. The gas-phase activations via H-abstraction by OH radicals also exhibit significant contributions. The reduced centers formed by reactions at vanadium oxo species are re-oxidized rapidly and, therefore, are present in very small concentrations at reaction conditions. The re-oxidation steps lead to the formation of HO 2 radicals and surface peroxo species that are also rapidly consumed and are present in small concentrations. The peroxo species preferentially convert C 2 H 4 to its epoxide product and influence selectivity even at low concentrations. The gas-phase reactions decrease the concentrations of peroxo species and improve selectivity slightly. The effects of reaction conditions and catalyst site density provide further insights into how factors beyond conversions at lattice oxygens influence rates and selectivity in alkane oxidation reactions of significant industrial importance.

42 ENGINEERING↗

DOE OSTI2023

Logarithmic scaling of higher-order temperature moments in the atmospheric surface layer

A generalized logarithmic law for high-order moments of homogeneous passive scalars is proposed for turbulent boundary layers. This law is analogous to the generalized log law that has been proposed for high-order moments of the turbulent longitudinal velocity and is derived by combining the random sweeping decorrelation hypothesis with a spectral model informed by the attached eddy hypothesis. The proposed theory predicts that the high-order moments of passive scalar fluctuations within the inertial sublayer will vary logarithmically with wall-normal distance (), and is evaluated using high frequency time-series measurements of temperature and streamwise velocity fluctuations obtained in the first meter of the atmospheric surface layer (ASL) under near-neutral thermal stratification. The logarithmic dependence with within the inertial sublayer is observed in both the air temperature and velocity moments, with good agreement to the predictions from the proposed theory. Surprisingly, the proposed theory appears to be as, if not more, valid for transported passive scalars than for the longitudinal velocity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The role of turbulent fluctuations in aerosol activation and cloud formation

Aerosol indirect effects are one of the leading contributors to cloud radiative properties relevant to climate. Aerosol particles become cloud droplets when the ambient relative humidity (saturation ratio) exceeds a critical value, which depends on the particle size and chemical composition. In the traditional formulation of this problem, only average, uniform saturation ratios are considered. Using experiments and theory, we examine the effects of fluctuations, produced by turbulence. Our measurements, from a multiphase, turbulent cloud chamber, show a clear transition from a regime in which the mean saturation ratio dominates to one in which the fluctuations determine cloud properties. The laboratory measurements demonstrate cloud formation in mean-subsaturated conditions (i.e., relative humidity <100%) in the fluctuation-dominant activation regime. The theoretical framework developed to interpret these measurements predicts a transition from a mean- to a fluctuation-dominated regime, based on the relative values of the mean and standard deviation of the environmental saturation ratio and the critical saturation ratio at which aerosol particles activate or become droplets. The theory is similar to the concept of stochastic condensation and can be used in the context of the atmosphere to explore the conditions under which droplet activation is driven by fluctuations as opposed to mean supersaturation. Overall, it provides a basis for future development of cloud droplet activation parameterizations that go beyond the internally homogeneous parcel calculations that have been used in the past.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Experimental and Computational Evaluation of Lipidomic In-Source Fragmentation as a Result of Postionization with Matrix-Assisted Laser Desorption/Ionization

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can provide spatially resolved molecular information about a sample. Recently, a postionization approach (MALDI-2) has been commercially integrated with MALDI-MSI, allowing for bettered sensitivity and consequent improved spatial resolution. While advantages of MALDI-2 have previously been established, we demonstrate here statistically increased in-source fragmentation (ISF) results from postionization with a commercial instrument. Via lipid standard analyses, known MALDI ISF pathways (e.g., loss of trimethylamine) were statistically increased in MALDI-2 compared to MALDI-1 (65–172% increase in fragmentation). Gas phase molecular modeling with density functional theory estimated that the most-weighted virtual orbitals to excite within lipids involve ester and phosphate bonds. Protonated lipid excitation energies are furthermore red-shifted compared to those of other adduct types [e.g., 254 nm for protonated PC(16:0/18:1)] and approach the MALDI-2 laser energy (266 nm). Analysis of rat brain homogenate detected statistically more positive-ion mode peaks with MALDI-2 (1090) than that with MALDI-1 (719), where Kernel density estimations showed that the majority of this enhancement occurs with low m/z ions (i.e., m/z 75–500). Taken together with the lipid standard data, these observations may indicate ISF due to postionization. Finally, while artifact contributions from matrix blanks were also noted, both experimental and computational data sets suggest that the overall extent of ISF is statistically increased in MALDI-2 compared to MALDI-1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Gravothermal phase transition, black holes and space dimensionality

In the framework of gravothermal evolution of an ideal monatomic fluid, I examine the dynamical instability of the fluid sphere in (N+1) dimensions by exploiting Chandrasekhar’s criterion to each quasi-static equilibrium along the sequence of the evolution. Once the instability is triggered, it would probably collapse into a black hole if no other interaction halts the process. From this viewpoint, the privilege of (3+1)-dimensional spacetime is manifest, as it is the marginal dimensionality in which the ideal monatomic fluid is stable but not too stable. Moreover, it is the unique dimensionality that allows stable hydrostatic equilibrium with positive cosmological constant. While all higher dimensional (N > 3) spheres are genuinely unstable. In contrast, in (2+1)-dimensional spacetime it is too stable either in the context of Newton’s theory of Gravity or Einstein’s General Relativity. It is well-known that the role of negative cosmological constant is crucial to have the Banados-Teitelboim-Zanelli (BTZ) black hole solution and the equilibrium configurations of a fluid disk. Owing to the negativeness of the cosmological constant, there is no unstable configuration for a homogeneous fluid disk to collapse into a naked singularity, which supports the Cosmic Censorship Conjecture. However, BTZ holes of mass M BTZ > 0 could emerge from collapsing fluid disks. The implications of spacetime dimensionality are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Tuning the Dimensionality of Excitons in Colloidal Quantum Dot Molecules

Electrically coupled quantum dots (QDs) can support unique optoelectronic properties arising from the superposition of single-particle excited states. Experimental methods for integrating colloidal QDs within the same nano-object, however, have remained elusive to the rational design. Here, we demonstrate a chemical strategy that allows for the assembling of colloidal QDs into coupled composites, where proximal interactions give rise to unique optoelectronic behavior. The assembly method employing “adhesive” surfactants was used to fabricate both homogeneous (e.g., CdS–CdS, PbS–PbS, CdSe–CdSe) and heterogeneous (e.g., PbS–CdS, CdS–CdSe) nanoparticle assemblies, exhibiting quasi-one-dimensional exciton fine structure. In addition, tunable mixing of single-particle exciton states was achieved for dimer-like assemblies of CdSe/CdS core–shell nanocrystals. The nanoparticle assembly mechanism was explained within the viscoelastic interaction theory adapted for molten-surface colloids. Furthermore, we expect that the present work will provide the synthetic and theoretical foundation needed for building assemblies of many inorganic nanocrystals.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Real-Frequency Response Functions at Finite Temperature

Building on previous developments [A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlanc, Phys. Rev. B 99, 035120 (2019); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B101, 125109 (2020); A. Taheridehkordi, S. H. Curnoe, and J. P. F. LeBlancPhys. Rev. B102, 045115 (2020), B. Holm and U. von Barth, Phys. Rev. B 57, 2108 (1998), J. Vičičevic and M. Ferrero, Phys. Rev. B 101, 075113 (2020)], we show that the diagrammatic Monte Carlo technique allows us to compute finite-temperature response functions directly on the real-frequency axis within any field-theoretical formulation of the interacting fermion problem. There are no limitations on the type and nature of the system’s action or whether partial summation and self-consistent treatment of certain diagram classes are used. In particular, by eliminating the need for numerical analytic continuation from a Matsubara representation, our scheme allows us to study spectral densities of arbitrary complexity with controlled accuracy in models with frequency-dependent effective interactions. Finally, for illustrative purposes we consider the problem of the plasmon linewidth in a homogeneous electron gas (jellium).

3-dimensional systems↗

DOE OSTI2022

Implementation of hybrid finite element method based transport solver in GRIFFIN

A new transport solver option based on the hybrid FEM (HFEM) was implemented in GRIFFIN, the MOOSE-based reactor analysis code, as an effort to support routine core design calculations for advanced reactor applications. The HFEM formulation with P{sub N} (spherical harmonics expansion), akin to the variational nodal method, is effective for solving a spatially homogenized problem with strong transport effect. The residual and Jacobian evaluations of the HFEM weak form were derived and successfully implemented in GRIFFIN, having the diffusion and the PN options available in the new HFEM based transport solver. The performance was tested with the simplified ABTR benchmark problems. The results indicate that the HFEM-based transport solver is a feasible option for solving problems with spatially homogenized and strong streaming by providing superior accuracy with a proper p-refinement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Direct Identification of Mixed-Metal Centers in Metal–Organic Frameworks: Cu 3 (BTC) 2 Transmetalated with Rh 2+ Ions

Raman spectroscopy was used to establish direct evidence of heterometallic metal centers in a metal–organic framework (MOF). The Cu 3 (BTC) 2 MOF HKUST-1 (BTC 3– = benzenetricarboxylate) was transmetalated by heating it in a solution of RhCl 3 to substitute Rh 2+ ions for Cu 2+ ions in the dinuclear paddlewheel nodes of the framework. In addition to the Cu–Cu and Rh–Rh stretching modes, Raman spectra of (Cu x Rh 1– x ) 3 (BTC) 2 show the Cu–Rh stretching mode, indicating that mixed-metal Cu–Rh nodes are formed after transmetalation. Density functional theory studies confirmed the assignment of a Raman peak at 285 cm –1 to the Cu–Rh stretching vibration. Electron paramagnetic resonance spectroscopy experiments further supported the conclusion that Rh 2+ ions are substituted into the paddlewheel nodes of Cu 3 (BTC) 2 to form an isostructural heterometallic MOF, and electron microscopy studies showed that Rh and Cu are homogeneously distributed in (Cu x Rh 1– x ) 3 (BTC) 2 on the nanoscale.

02 PETROLEUM↗

DOE OSTI2020

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

High-Rate, Selective Electrosynthesis of Cyclohexanone Oxime via In Situ Generation and Release of Hydroxylamine on Bismuth

Oxime compounds are key industrial intermediates for nylon precursors and commodity chemicals. However, conventional routes rely on multistep reactions and hydroxylamine (NH 2 OH) salts, raising significant safety and sustainability concerns. Although electrosynthesis offers an alternative, oxime formation on d-block transition metals suffers from poor selectivity, as nitrogen oxyanion intermediates bind strongly to the surface and are readily over-reduced to ammonia. Here, we report morphology-controlled p-block bismuth rhombic dodecahedra (Bi RDs) that promote in situ NH 2 OH generation and its desorption into the electrolyte, enabling an electrochemical-chemical decoupled route for cyclohexanone oxime (CHO) synthesis. Bi RDs deliver nearly 100% Faradaic efficiency (FE) at −0.5 V vs. RHE and a yield of 1.4 mmol h –1 cm –2 at −0.9 V vs. RHE in an H-cell, while maintaining a CHO selectivity of nearly 100% at 100 mA cm –2 in a flow cell. Under identical conditions, d-block electrodes (Cu, Pd, Ag) show FE below 30%. Density functional theory calculations reveal that Bi 6p orbital-derived surface states weaken intermediate binding and facilitate NH 2 OH desorption, suppressing over-reduction. Kinetic analysis, post-addition trapping experiments, and in situ ATR-FTIR and Raman spectroscopy suggest the following reaction mechanism: NH 2 OH is selectively generated at the electrode surface, released as a freely diffusing intermediate, and undergoes homogeneous condensation with cyclohexanone in the bulk electrolyte, bypassing the surface-confined Langmuir–Hinshelwood pathway. These findings demonstrate that regulating intermediate desorption through p-block orbital chemistry provides a general strategy for achieving high selectivity in electro-organic nitrogen synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Interplay between copper redox and transfer and support acidity and topology in low temperature NH 3 -SCR

Low-temperature standard NH 3 -SCR over copper-exchanged zeolite catalysts occurs on NH 3 -solvated Cu-ion active sites in a quasi-homogeneous manner. As key kinetically relevant reaction steps, the reaction intermediate Cu II (NH 3 ) 4 ion hydrolyzes to Cu II (OH)(NH 3 ) 3 ion to gain redox activity. The Cu II (OH)(NH 3 ) 3 ion also transfers between neighboring zeolite cages to form highly reactive reaction intermediates. Via operando electron paramagnetic resonance spectroscopy and SCR kinetic measurements and density functional theory calculations, we demonstrate here that such kinetically relevant steps become energetically more difficult with lower support Brønsted acid strength and density. Consequently, Cu/LTA displays lower Cu atomic efficiency than Cu/CHA and Cu/AEI, which can also be rationalized by considering differences in their support topology. By carrying out hydrothermal aging to eliminate support Brønsted acid sites, both Cu II (NH 3 ) 4 ion hydrolysis and Cu II (OH)(NH 3 ) 3 ion migration are hindered, leading to a marked decrease in Cu atomic efficiency for all catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗