Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic diffusion”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Effect of magnetic helicity upon rectilinear propagation of charged particles in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by any random field components that are superposed on the guiding field. If the random field configuration embodies helicity, the scattering is asymmetrical with respect to a plane perpendicular to the guiding field, for particles moving into the forward hemisphere are scattered at different rates from those moving into the backward hemisphere. This asymmetry gives rise to new terms in the transport equations that describe propagation of charged particles. Helicity has virtually no impact on qualitative features of the diffusive mode of propagation. However, characteristic velocities of the coherent modes that appear after a highly anisotropic injection exhibit an asymmetry related to helicity. Explicit formulas, which embody the effects of helicity, are given for the anisotropies, the coefficient diffusion, and the coherent velocities. Predictions derived from these expressions are in good agreement with Monte Carlo simulations of particle transport, but the simulations reveal certain phenomena whose explanation calls for further analytical work.

Earl, James A.↗
DOE OSTI2019

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2022

Multi-scale boron penetration toward stabilizing nickel-rich cathode

Nickel-rich layered oxides LiNi x Co y Mn 1-x-y O 2 (x≥0.8) have been recognized as the preferred cathode materials to develop lithium-ion batteries with high energy density (>300Wh kg -1 ). However, the poor cycling stability and rate capability stemming from intergranular cracks and sluggish kinetics hinder their commercialization. To address such issues, a multi-scale boron penetration strategy is designed and applied on the polycrystalline LiNi 0.83 Co 0.11 Mn 0.06 O 2 particles that are pre-treated with pore construction. The lithium-ion conductive lithium borate in grain gaps functions as the grain binder that can bear the strain/stress from anisotropic contraction/expansion, and provides more pathways for lithium-ion diffusion. As a result, the intergranular cracks are ameliorated and the lithium-ion diffusion kinetics is improved. Moreover, the coating layer separates the sensitive cathode surface and electrolyte, helping to suppress the parasitic reactions and related gas evolution. In addition, the enhanced structural stability is acquired by strong B-O bonds with trace boron doping. As a result, the boron-modified sample with an optimized boron content of 0.5% (B5-NCM) exhibits a higher initial discharge capacity of 205.5mAh g -1 at 0.1C (1C=200mA g -1 ) and improved capacity retention of 81.7% after 100 cycles at 1C. Furthermore, the rate performance is distinctly enhanced by high lithium-ion conductive LBO (175.6mAh g -1 for B5-NCM and 154.6mAh g -1 for B0-NCM at 5C).

25 ENERGY STORAGE↗
DOE OSTI2018

Tunable Assembly of Gold Nanorods in Polymer Solutions To Generate Controlled Nanostructured Materials

In this study, gold nanorods grafted with short chain polymers are assembled into controlled open structures using polymer-induced depletion interactions and structurally characterized using small angle x-ray scattering. When the nanorod diameter is smaller than the radius of gyration of the depletant polymer, the depletion interaction depends solely on the correlation length of the polymer solution and not directly on the polymer molecular weight. As the polymer concentration increases, the stronger depletion interactions increasingly compress the grafted chains and push the gold nanorods closer together. By contrast, other structural characteristics such as the number of nearest neighbors and fractal dimension exhibit a non-monotonic dependence on polymer concentration. These parameters are maximal at intermediate concentrations, which are attributed to a crossover from reaction-limited to diffusion-limited aggregation. Finally, the control over structural properties of anisotropic nanoscale building blocks demonstrated here will be beneficial to designing and producing materials in situ with specific direction-dependent nanoscale properties and provides a crucial route for advances in additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Theoretical prediction and atomic kinetic Monte Carlo simulations of void superlattice self-organization under irradiation

Nano-structured superlattices may have novel physical properties and irradiation is a powerful mean to drive their self-organization. However, the formation mechanism of superlattice under irradiation is still open for debate. Here we use atomic kinetic Monte Carlo simulations in conjunction with a theoretical analysis to understand and predict the self-organization of nano-void superlattices under irradiation, which have been observed in various types of materials for more than 40 years but yet to be well understood. The superlattice is found to be a result of spontaneous precipitation of voids from the matrix, a process similar to phase separation in regular solid solution, with the symmetry dictated by anisotropic materials properties such as one-dimensional interstitial atom diffusion. This discovery challenges the widely accepted empirical rule of the coherency between the superlattice and host matrix crystal lattice. The atomic scale perspective has enabled a new theoretical analysis to successfully predict the superlattice parameters, which are in good agreement with independent experiments. The theory developed in this work can provide guidelines for designing target experiments to tailor desired microstructure under irradiation. Furthermore, it may also be generalized for situations beyond irradiation, such as spontaneous phase separation with reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Local charge-screening and polarization switching in a ferroelectric ionic-superconductor

For ferroelectric ionic-conductors, polarization switching is complicated by the interplay between ion mobility and charge screening effects. When the ionic charge carriers also play a key role in the domain reversal, such as in Rb-doped KTiOPO 4 (RKTP), a higher level of complexity is introduced. RKTP provides an ideal platform for investigating the relationship between ionic conductivity and polarization reversal because its highly anisotropic crystal properties allow selective modification of material characteristics through diffusive cation doping. Here, we use indiffused Ba/K doping to create a significant increase in the ionic conductivity. Time-of-flight secondary ion mass spectrometry is employed to map Ba/K doping within the RKTP crystal and correlate it to changes in ionic mobility and polarization switching characteristics under an external field applied to the nonpolar face. Using band-excitation piezoresponse force microscopy, we demonstrate a selective switching-inhibition mechanism driven by the enhanced charge screening.

36 MATERIALS SCIENCE↗
NASA NTRS2022

The Extended Field-Aligned Suprathermal Proton Beam and Long-Lasting Trapped Energetic Particle Population Observed Upstream of a Transient Interplanetary Shock

The properties of the suprathermal particle distributions observed upstream of interplanetary shocks depend not only on the properties of the shocks but also on the transport conditions encountered by the particles as they propagate away from the shocks. The confinement of particles in close proximity to the shocks, as well as particle scattering processes during propagation to the spacecraft, lead to the common observation of upstream diffuse particle distributions. We present observations of a rare extended anisotropic low-energy(≤30 keV)proton beam together with a trapped ≥500 keV proton population observed in association with the arrival of an oblique interplanetary shock at the Advanced Composition Explorer, the Interplanetary Monitoring Platform-8, and the Wind spacecraft on 2001 January 31. Continuous injection of particles by the traveling shock into a smooth radial magnetic field region formed in the tail of a modest high-speed solar wind stream produced an extended foreshock region of energetic particles. The absence of enhanced magnetic field fluctuations upstream of the shock results in the observation of a prolonged anisotropic field-aligned beam of ≤30 keV protons as well as a population of higher-energy(≥500 keV)protons with small pitch-angle cosine(μ∼0)extending far from the shock.

D. Lario↗
DOE OSTI2018

Net-baryon diffusion in fluid-dynamic simulations of relativistic heavy-ion collisions

A hybrid (hydrodynamics + hadronic transport) theoretical framework is assembled to model the bulk dynamics of relativistic heavy-ion collisions at energies accessible in the Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC) and the NA61/SHINE experiment at CERN. The system's energy-momentum tensor and net baryon current are evolved according to relativistic hydrodynamics with finite shear viscosity and non-zero net baryon diffusion. Our hydrodynamic description is matched to a hadronic transport model in the dilute region. With this fully integrated theoretical framework, we present a pilot study of the hadronic chemistry, particle spectra, and anisotropic flow. Phenomenological effects of a non-zero net-baryon current and its diffusion on hadronic observables are presented for the first time. Furthermore, the importance of the hadronic transport phase is also investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-coercivity ThMn 12 -type monocrystalline Sm–Zr–Fe–Co–Ti particles by high-temperature reduction diffusion

The ThMn 12 -type (Sm,Zr)1(Fe,Co,Ti) 12 compounds have the potential of powerful permanent magnets. Magnetically hard and anisotropic powders of such compound have been prepared by subjecting elemental oxides and Co mixed with Ca and CaO dispersant to a succession of high-energy ball-milling, reduction diffusion at 990–1220 °C and repeated washing. The size of the resulting ThMn 12 -type crystallites, their coercivity and fraction of monocrystalline particles were all found to increase with the reduction-diffusion temperature. Particles synthesized at 1220 °C were highly monocrystalline with a mean size of 0.54 μm and, after a magnetic-field alignment, exhibited a coercivity of 1.26 T and a full-density-projected maximum energy product of at least 209 kJ/m 3 (26.3 MGOe). Here, the strong positive effect of the reduction-diffusion temperature on the coercivity has been attributed to separation of the crystallites and to decrease in the incidence of structural defects.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Effective diffusion equation in a random velocity field

The effects are studied of assumed random velocity fields on diffusion in a binary fluid. Random velocity fields can result, for example, from the high-frequency components of residual accelerations onboard spacecraft (often called g-jitter). An effective diffusion equation is derived for an average concentration which includes spatial and temporal correlations induced by the fluctuating velocity fields assumed to be Gaussianly distributed. The resulting equation becomes nonlocal, and if correlations between different components of the velocity field exist, it is also anisotropic. The simple limiting case of short correlation times is discussed and an effective diffusivity is obtained which reflects the enhanced mixing caused by the velocity fields. The results obtained in the limit of short correlation times are valid even if the probability distribution of the velocity field is not Gaussian.

Vinals, Jorge↗
NASA NTRS1974

The protons and electrons trapped in the Jovian dipole magnetic field region and their interaction with Io

Detailed analysis of electrons equal to or greater than 3 MeV and of protons 0.5 to 1.8 MeV and equal to or greater than 35 MeV for both the inbound and the outbound passes of the Pioneer 10 spacecraft. Conclusive evidence is obtained that the trapped radiation in Jupiter's inner magnetosphere is maintained and supplied by inward diffusion from the outer regions of the trapped radiation zone. It is shown that the time required for isotropization of an anisotropic flux by pitch angle scattering inside L approximately equal to 6 is long in comparison with the time required for particles to diffuse inward from L approximately equal to 6 to L approximately equal to 3, that the high-energy protons were not injected at high energies by the Crand (cosmic ray albedo neutron decay) process but were accelerated in the magnetosphere of Jupiter, and that the main conclusions of this analysis are unaffected by use of either the D sub 1 or the D sub 2 magnetic field models. Theoretical studies of the capture of trapped electrons and protons by Io have been carried out, and it is found that the probability of capture by Io depends strongly upon the particle species and kinetic energy.

Simpson, J. A.↗
DOE OSTI2015

The phase diagram of high-pressure superionic ice

Superionic ice is a special group of ice phases at high temperature and pressure, which may exist in ice-rich planets and exoplanets. In superionic ice liquid hydrogen coexists with a crystalline oxygen sublattice. At high pressures, the properties of superionic ice are largely unknown. Here we report evidence that from 280 GPa to 1.3 TPa, there are several competing phases within the close-packed oxygen sublattice. At even higher pressure, the close-packed structure of the oxygen sublattice becomes unstable to a new unusual superionic phase in which the oxygen sublattice takes the P2 1 /c symmetry. We also discover that higher pressure phases have lower transition temperatures. The diffusive hydrogen in the P2 1 /c superionic phase shows strong anisotropic behaviour and forms a quasi-two-dimensional liquid. The ionic conductivity changes abruptly in the solid to close-packed superionic phase transition, but continuously in the solid to P2 1 /c superionic phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Why is anatase a better photocatalyst than rutile? - Model studies on epitaxial TiO2 films

The prototypical photocatalyst TiO2 exists in different polymorphs, the most common forms are the anatase- and rutile-crystal structures. Generally, anatase is more active than rutile, but no consensus exists to explain this difference. Here we demonstrate that it is the bulk transport of excitons to the surface that contributes to the difference. Utilizing high –quality epitaxial TiO2 films of the two polymorphs we evaluate the photocatalytic activity as a function of TiO2-film thickness. For anatase the activity increases for films up to ,5 nm thick, while rutile films reach their maximum activity for ,2.5 nm films already. This shows that charge carriers excited deeper in the bulk contribute to surface reactions in anatase than in rutile. Furthermore, we measure surface orientation dependent activity on rutile single crystals. The pronounced orientation-dependent activity can also be correlated to anisotropic bulk charge carrier mobility, suggesting general importance of bulk charge diffusion for explaining photocatalytic anisotropies.

POLLUTION REMEDIATION↗
DOE OSTI2018

Influence of surface nano-patterning on the placement of InAs quantum dots

We have examined the influence of spontaneous nano-patterning on the placement of InAs quantum dots (QDs) on (Al)GaAs surfaces using an experimental-computational approach. Both atomically flat and mounded surfaces, generated via a surface instability induced by the Ehrlich- Schwoebel barrier, are employed as templates for the subsequent deposition of InAs QDs. Using height profiles from atomic-force micrographs, we simulate QD deposition with a 2D phase field model, which describes the time evolution of the InAs layer driven by a chemical potential gradient. For flat surfaces, phase-field simulations result in QD densities comparable to experimental observations. For mounded surfaces, the simulations reveal QDs preferentially positioned in regions of positive curvature (substrate valleys), e.g., at the edge of surface mounds, consistent with the anisotropic QD placement observed experimentally. We discuss the role of curvaturedriven diffusion in the spontaneous ordering of QDs, demonstrating the applicability of this mechanism to AlGaAs mounds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment and Improvement of Fission Product Transport Predictions of Particle Fuel in BISON

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This study was motivated by the need to incorporate more physics-based models in BISON in order to foster tri-structural isotropic (TRISO) applications. This document details the integration of new modeling capabilities in BISON, including (1) development of pyrolytic carbon (PyC) and silicon carbide (SiC) layer anisotropic thermal and mass transport capabilities, (2) verification of the mass diffusion solution in TRISO modeling, (3) calibration of fission product diffusivity using Advanced Gas Reactor (AGR) experiments, (4) improved fission product release modeling by developing compact diffusion modeling capabilities, and (5) documentation of accelerated failure analysis on the BISON website. Improvements made to the diffusion models and parameters were documented and validated against AGR-1 and -2 experiment data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗