Engineering Papers⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Literary review and modeling outlook for radiation damage in fused-silica optical fibers

The purpose of this document is to distill, from literature, the effects of radiation damage in fused-silica optical fibers. This information will be used to design a theoretical framework capable of predictive modeling that will be published elsewhere. By inspection of the relevant literature, the devised method is calculation of defect amorphous structures by molecular dynamic techniques, and computation of the macroscopic dielectric function using density functional theory. The dielectric function is then used to directly compute the transmittance of the core material before and after irradiation via solution of Maxwell’s equations in a homogeneous medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Stable Site‐Isolated Mono(phosphine)‐Rhodium Catalyst on a Metal‐Organic Layer for Highly Efficient Hydrogenation Reactions

Abstract Phosphine‐ligated transition metal complexes play a pivotal role in modern catalysis, but our understanding of the impact of ligand counts on the catalysis performance of the metal center is limited. Here we report the synthesis of a low‐coordinate mono(phosphine)‐Rh catalyst on a metal‐organic layer (MOL), P‐MOL • Rh, and its applications in the hydrogenation of mono‐, di‐, and tri‐substituted alkenes as well as aryl nitriles with turnover numbers (TONs) of up to 390000. Mechanistic investigations and density functional theory calculations revealed the lowering of reaction energy barriers by the low steric hindrance of site‐isolated mono(phosphine)‐Rh sites on the MOL to provide superior catalytic activity over homogeneous Rh catalysts. The MOL also prevents catalyst deactivation to enable recycle and reuse of P‐MOL • Rh in catalytic hydrogenation reactions.

Chemistry↗
DOE OSTI2022

Effect of loading path on grain misorientation and geometrically necessary dislocation density in polycrystalline aluminum under reciprocating shear

Solid phase processing (SPP) is a promising alloy fabrication technique to produce fine and homogeneous grain structures for high-performance alloys. However, there is very limited modeling capability to understand and predict the grain refinement during SPP. In this work, the crystal plasticity theory was used to study elastic-plastic deformation in polycrystalline aluminums under large shear deformation. Two approaches, kernel averaged misorientation (KAM) and grain reference orientation deviation (GROD), were used to assess the grain misorientations. The geometrically necessary dislocation (GND) density was computed with the plastic strain rate. The deformation simulations were carried out under two loading conditions to investigate the effect of loading paths on the evolutions of grain misorientation and GND density. The results show that the regions with high misorientation and GND density first appear near grain boundaries. These regions then extend toward interior grains. The loading path affects dislocation system activation and dislocation recovery, hence dislocation evolution and misorientation. In conclusion, both two- and three-dimensional simulations showed that the spatial and temporal evolutions of GROD, KAM, and GND density in are closely correlated, which indicates they all can be used as criteria of grain refinement or recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗
DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
DOE OSTI2021

Toward the Development of New Parameterizations for Surface Fluxes

For more than 50 years Monin-Obukhov Similarity Theory (MOST) has been used to represent near-surface exchanges of heat, moisture, and momentum in NWP models. However, limitations of MOST are well-known. In brief, MOST assumes a homogenous near-surface flux layer, which is oftentimes not valid, and suffers from autocorrelation. Limitations such as these motivate the need to modify the functional forms of the similarity equations used to represent fluxes in NWP models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Highly Asymmetric Graphene Layer Doping and Band Structure Manipulation in Rare Earth–Graphene Heterostructure by Targeted Bonding of the Intercalated Gadolinium

Heterostructures consisting of vertically stacked two-dimensional (2D) materials have recently gained large attention due to their highly controllable electronic properties and resulting quantum phases. In contrast to the mechanically stacked multilayered systems, which offer exceptional control over a stacking sequence or interlayer twist angles, the epitaxially grown 2D materials express unprecedented quality and stability over wafer-scale lengths. However, controlling the growth conditions remains a major obstacle toward the formation of complex, epitaxial heterostructures with well-defined electronic properties. Here, we synthesized a trilayer graphene heterostructure on the SiC(0001) substrate with two specific interlayer locations occupied by gadolinium. We applied multitechnique methodology based on low-temperature scanning tunneling microscopy/spectroscopy (STM/S) and angle-resolved photoelectron spectroscopy (ARPES) to determine the intercalant’s locations in the complex, epitaxial graphene heterostructure. Our approach relies on very high quality and large, micrometer-scale homogeneity of the synthesized system. The experimentally determined electronic structure is dominated by the two topmost graphene layers. Overall, our spectroscopic results show quantitative agreement between global ARPES, local STM/S, and density functional theory predictions. The characterized electronic properties primarily reflect highly anisotropic doping levels between the two corresponding graphene layers, which significantly affect the band structure topology. Two pairs of hybridized massive Dirac bands from our initial synthesis–the bilayer graphene on the SiC(0001) substrate–are transformed upon Gd intercalation into two pairs of massless Dirac bands with a new hybridization region in between. Our results open perspectives in the realization of exotic 2D quantum materials via atomically precise synthesis of epitaxial, multilayered graphene–rare earth heterostructures.

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

Superfluid dark matter around black holes

Superfluid dark matter, consisting of self-interacting light particles that thermalize and condense to form a superfluid in galaxies, provides a novel theory that matches the success of the standard ΛCDM model on cosmological scales while simultaneously offering a rich phenomenology on galactic scales. Within galaxies, the dark matter density profile consists of a nearly homogeneous superfluid core surrounded by an isothermal envelope. In this work we compute the density profile of superfluid dark matter around supermassive black holes at the center of galaxies. Here, in this work, we show that, depending on the fluid equation of state, the dark matter profile presents distinct power-law behaviors, which can be used to distinguish it from the standard results for collisionless dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Mixed stochastic-deterministic time-dependent density functional theory: application to stopping power of warm dense carbon

Warm dense matter (WDM) describes an intermediate phase, between condensed matter and classical plasmas, found in natural and man-made systems. In a laboratory setting, WDM is often created dynamically. It is typically laser or pulse-power generated and can be difficult to characterize experimentally. Measuring the energy loss of high energy ions, caused by a WDM target, is both a promising diagnostic and of fundamental importance to inertial confinement fusion research. However, electron coupling, degeneracy, and quantum effects limit the accuracy of easily calculable kinetic models for stopping power, while high temperatures make the traditional tools of condensed matter, e.g. time-dependent density functional theory (TD-DFT), often intractable. In this study, we have developed a mixed stochastic-deterministic approach to TD-DFT which provides more efficient computation while maintaining the required precision for model discrimination. Recently, this approach showed significant improvement compared to models when compared to experimental energy loss measurements in WDM carbon. Here, we describe this approach and demonstrate its application to warm dense carbon stopping across a range of projectile velocities. We compare direct stopping-power calculation to approaches based on combining homogeneous electron gas response with bound electrons, with parameters extracted from our TD-DFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mechanistic Insights into Co and Fe Quaterpyridine-Based CO 2 Reduction Catalysts: Metal–Ligand Orbital Interaction as the Key Driving Force for Distinct Pathways

Both [Co II (qpy)(H 2 O) 2 ] 2+ and [Fe II (qpy)(H 2 O) 2 ] 2+ (with qpy = 2,2':6',2'':6'',2'''-quaterpyridine) are efficient homogeneous electrocatalysts and photoelectrocatalysts for the reduction of CO 2 to CO. The Co catalyst is more efficient in the electrochemical reduction, while the Fe catalyst is an excellent photoelectrocatalyst. This work uses density functional theory to shed light on the contrasting catalytic pathways. While both catalysts experience primarily ligand-based reductions, the second reduction in the Co catalyst is delocalized onto the metal via a metal-ligand bonding interaction, causing a spin transition and a distorted ligand framework. This orbital interaction explains the experimentally observed mild reduction potential and slow kinetics of the second reduction. The decreased hardness and doubly occupied d z 2 -orbital facilitate a σ-bond with the CO 2 -π* in an η 1 - κC binding mode. CO 2 binding is only possible after two reductions resulting in an EEC mechanism (E = electron transfer, C = chemical reaction), and the second protonation is rate-limiting. In contrast, the Fe catalyst maintains a Lewis acidic metal center throughout the reduction process because the metal orbitals do not strongly mix with the qpy-π* orbitals. Further, this allows binding of the activated CO 2 in an η 2 -binding mode. This interaction stabilizes the activated CO 2 via a π-type interaction of a Fe-t 2g orbital and the CO 2 -π* and a dative bond of the oxygen lone pair. This facilitates CO 2 binding to a singly reduced catalyst resulting in an ECE mechanism. The barrier for CO 2 addition and the second protonation are higher than those for the Co catalyst and rate-limiting.

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

Alluvial fan morphology: A self-similar free boundary problem description

In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.

58 GEOSCIENCES↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Space-Dependent Calculation of the Multiplicity Moments for Shells With the Inclusion of Scattering

In recent work, we extended the methodology of multiplicity counting in nuclear safeguards by elaborating the one-speed stochastic transport theory of the calculation of the so-called multiplicity moments, i.e., the factorial moments of the number of neutrons emitted from a fissile item, following a source event from an internal neutron source [spontaneous fission and (α; n) reactions]. Calculations were made for solid spheres and cylinders, with the source being homogeneously distributed within the item. Recent measurements of the Rocky Flats Shells during the Measurement of Uranium Subcritical and Critical (MUSIC) campaign conducted by Los Alamos National Laboratory and assisted by the University of Michigan inspired us to extend the model to spherical shell geometry with a point source in the middle of the central cavity. Comparison of the calculated results with the experimental ones indicated that accounting for fission as the only neutron reaction (the standard procedure in the point model, adapted also in our work so far) was not sufficient for reaching good agreement with measurements. The model was therefore extended to include elastic scattering into the one speed formalism, whereas the effect of inelastic scattering was accounted for in an empirical way. After these extensions, good agreement was found between the calculated and the measured values. The paper describes the extension of the theory and provides concrete quantitative results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Carbene Functionalization of Monolayer Tungsten Disulfide for Enhanced Quantum Emission

Semiconducting two-dimensional (2D) transition metal dichalcogenides (TMDs) are promising materials for an array of applications, ranging from conventional field-effect transistors, photodetectors, and light-emitting diodes to their more recent use in quantum photonic technologies. Chemical functionalization of 2D TMDs with organic ligands and adlayers provides an additional means for customizing their electronic and optical properties. While many pathways have been reported for the chemical functionalization of 2D TMDs, their frequent reliance on solution-based methods results in limited control over adlayer thickness and coverage, thus hindering utility in high-performance applications. Here, in this study, we describe the vapor-phase functionalization of a 2D TMD with carbene ligands, specifically tungsten disulfide (WS 2 ) with N-heterocyclic carbenes (NHCs), resulting in molecularly smooth, thin, and uniform adlayers. Reacting NHCs with monolayer WS 2 reduces the broad photoluminescence background observed at cryogenic temperatures by 58%, which facilitates the detection of single-photon emitters from strained monolayer WS 2 , as indicated by second order correlation values ( g (2) ) as low as 0.17 ± 0.07. Chemical characterization coupled with density functional theory calculations suggests that the NHC adlayer has a dual defect-passivation and doping effect on monolayer WS 2 that results in enhanced single-photon emission. Overall, this study establishes vapor-phase carbene functionalization as a homogeneous surface modification scheme for tailoring the quantum emission properties of semiconducting 2D TMDs.

N-heterocyclic carbenes↗
DOE OSTI2025

Electron Transfer Theory Elucidates the Hidden Role Played by Triethylamine and Triethanolamine during Photocatalysis

Triethylamine (TEA) and triethanolamine (TEOA) are renowned, in part, for their ability to reductively quench excited states by outer-sphere electron transfer with vast and still growing applications as sacrificial electron donors for photocatalytic systems. Upon amine oxidation, the resulting TEA •+ and TEOA •+ radical cations undergo proton transfer (or hydrogen atom transfer), resulting in the formation of a chemical reductant that has an α-carbon centered radical adjacent to the nitrogen center (TEA • and TEOA • ). In this contribution, we have electrochemically and spectroscopically characterized a set of electron acceptors which, upon accepting an electron, are a series of photocatalysts, [ReCl(R 1 R 2 -bpy)(CO) 3 ] •– , where R 1 and R 2 are electron-donating and electron-withdrawing groups in the 4,4′- and 5,5′-positions on the bipyridyl ligand. Here, we substantiated the formation of the electron donors, TE(O)A • , by spin trapping electron paramagnetic resonance spectroscopy, where TE(O)A • reacts with 2,4,6-tri-tert-butylnitrosobenzene to generate N-centered and O-centered radical adducts. Having established the chemical behaviors of the electron acceptors and donors individually, the electron transfer rate constants were determined across a 1.43 V range in driving force. The redox potential of TEA • was benchmarked to within ±80 mV on an absolute scale in V vs Fc + /Fc in CH 3 CN by using an empirical rate vs free-energy correlation, electron transfer theory, and density functional theory calculations. The equilibrium potentials for TEA • and TEOA • were determined to be −1.98 V and −1.76 V, respectively. Based on the kinetic and thermochemical analysis presented for TEA • and TEOA • , these transient radicals can be broadly considered strong homogeneous chemical reductants within the wider context of photoredox potentials. Thus, this work clarifies a frequently unnoticed secondary function for these sacrificial electron donors during photocatalysis and rationalizes the possibility of a one-photon/two-electron conversion process that is dependent on the free-energy exchange between TE(O)A • and photocatalysts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Copper carbon composite wire with a uniform carbon dispersion made by friction extrusion

Copper-carbon composites are a group of materials with excellent mechanical, electrical, thermal, and tribological properties. However, bulk size copper-carbon composites made by the traditional manufacturing processes, like rolling or extrusion, fall short of reaching some of these properties predicted by theory or demonstrated only by samples at centimeter scale or smaller. The two main challenges to the successful scaling-up are: 1) to uniformly disperse carbon in the metal matrix; 2) to avoid degradation due to oxidation or reaction from overheating. In this work, we first demonstrate friction extrusion as a new method to make bulk-size void-free copper-carbon composite wires with homogenized carbon dispersion. Three different carbon varieties, graphite powder, graphene nanopowder, and carbon nanotubes, were added to the copper matrix with the concentration ranging from 0.5 wt% to 15 wt%. Special tooling, processing parameters, and procedures were developed, especially for high carbon content samples. Ten-fold reductions of both copper grain size and carbon particle size were achieved and attributed to the high shear deformation. Energy dispersive X-ray spectrometry indicates the carbon powder was refined to a sub-micron level and uniformly dispersed in the copper matrix. Compared with that of pure copper, the thermal capacity of the composite wire increases by 30 % while density reduces by 29 %.

Friction extrusion, Metal matrix composite, Graphe↗