Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Markovnikov alcohols via epoxide hydroboration by molecular alkali metal catalysts

Synthesis of branched “Markovnikov” alcohols is crucial to various chemical industries. The catalytic reduction of substituted epoxides under mild conditions is a highly attractive method for preparing such alcohols. Classical methods based on heterogeneous or homogeneous transition metal-catalyzed hydrogenation, hydroboration, or hydrosilylation usually suffer from poor selectivity, reverse regioselectivity, limited functional group compatibility, high cost, and/or low availability of the catalysts. Here we report the discovery of highly regioselective hydroboration of nonsymmetrical epoxides catalyzed by ligated archetypal reductants in organic chemistry – alkali metal triethylborohydrides. The chemoselectivity and turnover efficiencies of the present catalytic approach are excellent. Thus, terminal and internal epoxides with ene, yne, aryl, and halo groups were selectively and quantitatively reduced under a substrate-to-catalyst ratio (S/C) of up to 1000. Mechanistic investigations point to a mechanism reminiscent of frustrated Lewis pair action on substrates in which a nucleophile and Lewis acid act cooperatively on the substrate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

3D printed transparent ceramic YAG laser rods: Matching the core-clad refractive index

Yttrium Aluminum Garnet (YAG) solid state laser gain media rods with an active Neodymium-doped core and an optically-clear cladding region were additively manufactured via direct-ink-writing (DIW), followed by sintering and hot isostatic pressing to form fully dense optical ceramics. Lutetium and Gadolinium were chosen as optically-inert co-doping ions in the clad to match the increase in refractive index caused by the Neodymium in the core. Additionally, either 11.6% Lutetium or 3.8% Gadolinium can be used to match the index change from 2% neodymium; however, differences in diffusion distances across the core-clad interface lead to large fluctuations in index in that region. These index fluctuations can be minimized either by matching dopants with similar diffusion distances, or by implementing a gradual gradient in the doping profile, possible through DIW, rather than a sharp compositional interface attainable via more standard fabrication methods. This improvement in index homogeneity resulted in a 40% improvement in lasing performance compared with that of a core-clad rod fabricated with a sharp interface between the doped and undoped regions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interfacial Electrochemistry of Catalyst-Coordinated Graphene Nanoribbons

The immobilization of molecular electrocatalysts on conductive electrodes is an appealing strategy for enhancing their overall activity relative to those of analogous molecular compounds. Here, in this study, we report on the interfacial electrochemistry of self-assembled two-dimensional nanosheets of graphene nanoribbons (GNR-2DNS) and analogs containing a Rh-based hydrogen evolution reaction (HER) catalyst (RhGNR-2DNS) immobilized on conductive electrodes. Proton-coupled electron transfer (PCET) taking place at N-centers of the nanoribbons was utilized as an indirect reporter of the interfacial electric fields experienced by the monolayer nanosheet located within the electric double layer. The experimental Pourbaix diagrams were compared with a theoretical model, which derives the experimental Pourbaix slopes as a function of parameter f, a fraction of the interfacial potential drop experienced by the redox-active group. Interestingly, our study revealed that GNR-2DNS was strongly coupled to glassy carbon electrodes (f = 1), while RhGNR-2DNS was not (f = 0.15). We further investigated the HER mechanism by RhGNR-2DNS using electrochemical and X-ray absorption spectroelectrochemical methods and compared it to homogeneous molecular model compounds. RhGNR-2DNS was found to be an active HER electrocatalyst over a broader set of aqueous pH conditions than its molecular analogs. We find that the improved HER performance in the immobilized catalyst arises due to two factors. First, redox-active bipyrimidine-based ligands were shown to dramatically alter the activity of Rh sites by increasing the electron density at the active Rh center and providing RhGNR-2DNS with improved catalysis. Second, catalyst immobilization was found to prevent catalyst aggregation that was found to occur for the molecular analog in the basic pH. Overall, this study provides valuable insights into the mechanism by which catalyst immobilization can affect the overall electrocatalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

On the three-dimensional spatial correlations of curved dislocation systems

Abstract Coarse-grained descriptions of dislocation motion in crystalline metals inherently represent a loss of information regarding dislocation-dislocation interactions. In the present work, we consider a coarse-graining framework capable of re-capturing these interactions by means of the dislocation-dislocation correlation functions. The framework depends on a convolution length to define slip-system-specific dislocation densities. Following a statistical definition of this coarse-graining process, we define a spatial correlation function which will allow the arrangement of the discrete line system at two points—and thus the strength of their interactions at short range—to be recaptured into a mean field description of dislocation dynamics. Through a statistical homogeneity argument, we present a method of evaluating this correlation function from discrete dislocation dynamics simulations. Finally, results of this evaluation are shown in the form of the correlation of dislocation densities on the same slip-system. These correlation functions are seen to depend weakly on plastic strain, and in turn, the dislocation density, but are seen to depend strongly on the convolution length. Implications of these correlation functions in regard to continuum dislocation dynamics as well as future directions of investigation are also discussed.

Anderson, Joseph Pierre (ORCID:0000000329760903)↗
DOE OSTI2024

Influence of Joule Heating on the Stability of High Temperature Irradiation-resistant Thermocouples

Development of in-core instrumentation is driven by the pursuit of safer, more economic energy production from the perspective of both existing nuclear reactors and Generation IV reactor designs. Idaho National Laboratory has developed high-temperature irradiation-resistant thermocouples (HTIR-TCs) for temperature sensing inside Generation IV nuclear reactors. These thermocouples are composed of phosphorus-doped niobium (Nb-P) and lanthana-doped molybdenum (Mo-LaO) thermoelements, an alumina (Al2O3) insulation, and a niobium sheath. HTIR-TCs require an initial heat treatment exceeding the maximum service temperature to stabilize the generated electromotive force (EMF) signal. The mechanism behind the stabilization of the HTIR-TCs through traditional heat treatment methods is understood; however, the traditional heat treatment method is expensive, time consuming, and results in a heterogeneous microstructure. Therefore, we investigated a rapid method for stabilization and microstructure homogeneity, through Joule heating. This work evaluates the impact of Joule heating on the thermoelements’ microstructures, chemical stability, and mechanical properties so as to determine the mechanisms by which stabilization of the EMF signal occurs. Accordingly, during the Joule heat treatment, a secondary Nb3P phase coarsened along the length of the Nb-P thermoelement, along with the formation of an interaction region at the Al2O3/niobium interface. The interaction between the alumina insulation and the Nb-P thermoelement was also observed within the Mo-LaO thermoelement. Joule heating induced stability within the generated HTIR-TC EMF signal via the formation of secondary phases within the Nb-P and the interaction between the alumina insulation and the thermoelements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2026

Subplane Decusping for BWRs in MPACT

Control blade cusping can introduce significant error in boiling water reactor (BWR) calculations with MPACT when blade tips fall partway within an axial method of characteristics (MOC) plane, requiring homogenization of controlled and uncontrolled regions. This work implements subplane decusping for BWRs in MPACT by enabling BWR-compatible subplane coarse mesh finite difference (CMFD) and extending the existing decusping framework to represent between-assembly control blades that insert from the bottom of the core. The method resolves axial heterogeneity on a refined subplane mesh in the low-order solve and uses the resulting subplane flux shape to form flux-volume homogenized transport cross sections for the partially rodded MOC plane. The capability is evaluated using a single physics General Electric (GE)-14 assembly and a multiphysics Peach Bottom Unit 2 Type 1 assembly (PB2T1A) with thermal hydraulic feedback. In both cases, coarse axial meshes with and without subplane decusping are compared against fine mesh reference solutions over the full range of blade withdrawal positions. Subplane decusping reduces maximum/average eigenvalue errors from 3,369/220 pcm to 303/36 pcm for GE-14 and from 6,689/585 pcm to 1,582/109 pcm for PB2T1A. Additionally, it reduces pin power root mean square errors from 5.2%/0.8% to 1.8%/0.2% for GE-14 and from 6.9%/0.6% to 2.3%/0.1% for PB2T1A. These results demonstrate an effective, practical correction for BWR blade cusping in VERA-MPACT.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Unveiling the Universe with emerging cosmological probes

The detection of the accelerated expansion of the Universe has been one of the major breakthroughs in modern cosmology. Several cosmological probes (Cosmic Microwave Background, Supernovae Type Ia, Baryon Acoustic Oscillations) have been studied in depth to better understand the nature of the mechanism driving this acceleration, and they are being currently pushed to their limits, obtaining remarkable constraints that allowed us to shape the standard cosmological model. In parallel to that, however, the percent precision achieved has recently revealed apparent tensions between measurements obtained from different methods. These are either indicating some unaccounted systematic effects, or are pointing toward new physics. Following the development of CMB, SNe, and BAO cosmology, it is critical to extend our selection of cosmological probes. Novel probes can be exploited to validate results, control or mitigate systematic effects, and, most importantly, to increase the accuracy and robustness of our results. This review is meant to provide a state-of-art benchmark of the latest advances in emerging “beyond-standard” cosmological probes. We present how several different methods can become a key resource for observational cosmology. In particular, we review cosmic chronometers, quasars, gamma-ray bursts, standard sirens, lensing time-delay with galaxies and clusters, cosmic voids, neutral hydrogen intensity mapping, surface brightness fluctuations, stellar ages of the oldest objects, secular redshift drift, and clustering of standard candles. The review describes the method, systematics, and results of each probe in a homogeneous way, giving the reader a clear picture of the available innovative methods that have been introduced in recent years and how to apply them. The review also discusses the potential synergies and complementarities between the various probes, exploring how they will contribute to the future of modern cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Lignin-based polymers with enhanced melt extrusion ability

A solid polymer blend material comprising: (i) a lignin-acrylonitrile component containing a homogeneous blend of a lignin component and an acrylonitrile-containing rubber component; and (ii) a styrene-containing thermoplastic component that is non-elastomeric; wherein components (i) and (ii) are homogeneously dispersed in the polymer blend material. Methods for producing the blend material are also described. Methods for producing objects made of the blend material by melt extrusion are also described, comprising: (a) melt blending components (i) and (ii) to form a polymer blend in which components (i) and (ii) are homogeneously blended, wherein the polymer blend exhibits a melt viscosity of no more than 2000 Pa·s at a shear rate of 100-1000 s −1 and when heated to a temperature of no more than 240° C.; and (b) forming an object made of said polymer blend material.

Naskar, Amit K.↗
DOE OSTI2023

Lignin-based polymers with enhanced melt extrusion ability

A solid polymer blend material comprising: (i) a lignin-acrylonitrile component containing a homogeneous blend of a lignin component and an acrylonitrile-containing rubber component; and (ii) a styrene-containing thermoplastic component that is non-elastomeric; wherein components (i) and (ii) are homogeneously dispersed in the polymer blend material. Methods for producing the blend material are also described. Methods for producing objects made of the blend material by melt extrusion are also described, comprising: (a) melt blending components (i) and (ii) to form a polymer blend in which components (i) and (ii) are homogeneously blended, wherein the polymer blend exhibits a melt viscosity of no more than 2000 Pa·s at a shear rate of 100-1000 s −1 and when heated to a temperature of no more than 240° C.; and (b) forming an object made of said polymer blend material.

Naskar, Amit K.↗
DOE OSTI2023

Computational design of metamaterials with self contact

Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.

42 ENGINEERING↗
DOE OSTI2024

Mechanical and Electrochemical Implications of Drying Temperature on Lithium‐Ion Battery Electrodes

Abstract Lithium‐ion battery (LIB) electrodes are typically produced with n‐methyl‐2‐pyrrolidone, a toxic solvent that is a known carcinogen and reproductive hazard. Accordingly, aqueous processing has been an expanding area of research interest in the field of LIB manufacturing. Although aqueous processing has been widely successful in anode processing, serious challenges remain in processing the cathode. In this work, the drying mechanics of cathode processed with both solvents is investigated though implementation of a chemical‐engineering‐based model to better understand the utilization of heat provided by experimentally determining the heat and mass transfer coefficients. Electrochemical performance is also evaluated to determine the impact of drying temperature on cycling performance. Binder distribution is determined via various methods to confirm differences in binder homogeneity as a function of both solvent and drying temperature. Identified is the large difference in the efficiency in which the heat is used as well as an ideal drying temperature for both aqueous and non‐aqueous processed cathodes. Also identified is the increased sensitivity to processing temperature for aqueous processed electrodes compared to non‐aqueous processed counterparts, pointing to the possibility of tuned drying regimes which would capitalize on the potential cost savings of aqueous processing for cathodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗