Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

New discretization of Caputo–Fabrizio derivative

We derive a numerical approximation, namely L1–2 formula, to the Caputo–Fabrizio derivative by using a quadratic interpolation. Quadratic and cubic convergence rates are achieved for L1 and L1–2 formulas using Lagrange interpolation, respectively. We compute Caputo–Fabrizio derivatives of some known functions both theoretically and numerically. In addition, we solve non/linear sub-diffusion equations to test theoretical findings. Numerical results confirm the theoretically observed convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

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

Electrostatic Barriers to Nanoparticle Accessibility of a Porous Matrix

Translocation from one cavity to another through a narrow constriction (i.e., a “hole”) represents the fundamental elementary process underlying hindered mass transport of nanoparticles and macromolecules within many natural and synthetic porous materials, including intracellular environments. This process is complex and may be influenced by long-range (e.g., electrostatic) particle–wall interactions, transient adsorption/desorption, surface diffusion, and hydrodynamic effects. In this study, we used a three-dimensional (3D) tracking method to explicitly visualize the process of nanoparticle diffusion within periodic porous nanostructures, where electrostatic interactions were mediated via ionic strength. The effects of electrostatic interactions on nanoparticle transport were surprisingly large. For example, an increase in the Debye length of only a few nanometers (in a material with a hole diameter of ~100 nm) increased the mean cavity escape time 3-fold. A combination of computational and experimental analyses indicated that this hindered cavity escape was due to an electrostatic energy barrier in the region of the hole, which was quantitatively explained using DLVO theory. These findings explicitly demonstrate that the cavity escape process was barrier-limited and dominated by electrostatic effects.

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

Diffusive Escape of a Nanoparticle from a Porous Cavity

Narrow escape from confinement through a nano-channel is the critical step of complex transport processes including size exclusion-based separations, oil and gas extraction from the microporous subsurface environment, and RNA translocation through nuclear pore complex channels. While narrow escape has been studied using theoretical and computational methods, experimental quantification is rare because of the difficulty in confining a particle into a microscopic space through a nanoscale hole. In this paper, we studied narrow escape in the context of continuous nanoparticle diffusion within the liquid-filled void space of an ordered porous material. Specifically, we quantified the spatial dependence of nanoparticle motion and the sojourn times of individual particles in the interconnected confined cavities of a liquid-filled inverse opal film. We found that nanoparticle motion was inhibited near cavity walls and cavity escape was slower than predicted by existing theories and random-walk simulations. Here, a combined computational-experimental analysis indicated that translocation through a nano-channel is barrier-controlled rather than diffusion-controlled.

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

Features of the dynamics of particles and fields at cyclotron resonances

Some features of the dynamics of particles and fields at cyclotron resonances have been discussed with the focus on chaotic dynamical regimes. It has been shown that the known criterion of the transition of the regular dynamics of particles to chaotic dynamics at cyclotron resonances sometimes describes this transition incorrectly. The reason for such a feature of the criterion has been revealed. The anomalous sensitivity of the dynamics of particles to external fluctuations at autoresonance has been analyzed. A theory of excitation of electromagnetic waves by a beam of phased oscillators under the conditions of isolated nonlinear cyclotron resonance has been developed. It has been shown that the chaotic dynamical regime is due to the periodic change in the structure of the phase portrait of particles in the wave field. It has been shown that higher moments can play a more significant role than lower moments in almost all chaotic dynamical regimes at cyclotron resonances. In this case, the known kinetic diffusion equations should be generalized with the inclusion of these higher moments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Nuclear Effects in Neutrino Detection

Although the interactions of neutrinos with elementary particles are described with impressive precision by the Standard Model, neutrino-nucleus reactions remain less well understood. Improvements in our knowledge of neutrino-nucleus scattering will be an important requirement for the success of future experimental efforts involving neutrinos at both low ($\sim$\SI{10}{\MeV}) and medium ($\sim$\SI{1}{\GeV}) energies. In the low energy regime, the planned DUNE experiment will attempt to observe neutrinos from a nearby core-collapse supernova using liquid argon time projection chambers (LArTPCs). Unlike other large neutrino detectors, which are primarily sensitive to electron antineutrinos, DUNE will observe mostly charged current absorptions of electron neutrinos on \isotope[40]{Ar} in response to a supernova, providing a unique window into the physics of stellar collapse. Despite the importance of low-energy neutrino-nucleus reactions to DUNE's supernova physics goals, prior to the work presented in this thesis, no thorough consideration of the many possible final states generated by neutrino-argon scattering, including those involving the emission of nucleons or heavier nuclear fragments, had yet been attempted in the literature for the energy range of interest for supernova neutrinos. To aid DUNE's supernova physics program, this thesis presents a detailed theoretical model of low-energy neutrino-argon scattering. This model has been implemented within a new event generator called MARLEY (Model of Argon Reaction Low Energy Yields) which may be used to simulate realistic $\isotope[40]{Ar}(\nu_e,e^-)\isotope[40]{K}^*$ events for LArTPC supernova neutrino sensitivity studies. At medium energies, large theoretical uncertainties in predictions of neutrino-induced neutron production present a problem for precision neutrino oscillation experiments, attempts to discover the Diffuse Supernova Neutrino Background, and searches for proton decay. To const rain the widely-varying predictions of current nuclear model! s, the Accelerator Neutrino Neutron Interaction Experiment (ANNIE) collaboration plans to make a systematic measurement of neutron production by neutrino interactions in water. Because the ANNIE experimental hall is located only \SI{100}{\meter} away from the Fermilab Booster Neutrino Beam target, background neutrons correlated in time with the beam could potentially interfere with this proposed measurement. As a first step toward ANNIE's ultimate physics goals, this thesis presents an analysis of background neutron rates as a function of position within the ANNIE detector. These rates are found to be small enough for the neutron yield measurements to proceed as planned.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

CEED-MS8:Initial Integration of CEED Software in ECP Applications

Key components of CEED software development involve fast finite element operator storage and evaluation, architecture optimizations, performant algorithms for all orders, global kernels for finite element operators, efficient use of the memory sub-system, optimal data locality and motion, enhanced scalability and parallelism, and fast tensor contractions. As ultimate goals of the project, CEED is exploring and identifying the best algorithms for the full range of discretizations and applying algorithmic and software development to support ECP applications’ needs. In this milestone we report on the initial integration of CEED software in the selected first-wave ECP/CEED apps and the continued exploration of other ECP applications as second-wave (milestone CEED-MS23) and third-wave (milestone CEED-MS35) ECP/CEED application candidates. An extensive evaluation of various approaches for GPU acceleration of the CEED computational motives is also a major focus of the report, both in the settings of CEED’s Nekbone and Laghos miniapps, as well as versions of CEED’s bake-off problems (introduced in milestone CEED-MS6). The first-wave ECP applications targeted for coupling with CEED software include the ExaSMR and MARBL applications. We also report on near-term collaboration with the seed-funded Urban Systems project, as well as engagements with ACME, ExaAM, GEOS and the Application Assessment and Proxy App projects. The ExaSMR project involves coupled thermal-hydraulics/neutronics, with the latter based on Monte Carlo methods and the former based on high-order discretizations of the Navier-Stokes and energy equations for accurate simulation of turbulent transport. The project requires fast and efficient turbulence simulation capabilities that will scale to exascale platforms. Work within CEED will ensure that all aspects of the Nek5000 workflow scale to the target problem space, which involves > 10 10 degrees of freedom in typical production runs. CEED’s goal is to provide implementations that perform with maximum attainable efficiency and deliver low turn-around times. A deliverable for this milestone is beta-release of the GPU-based variant of Nek5000. MARBL is targeting high-energy-density physics problems, including single fluid multi-material hydrodynamics and radiation/magnetic diffusion simulation, with applications in inertial confined fusion, pulsed power experiments, and equation of state/material strength analysis. The code uses high-order methods based arbitrary Lagrangian-Eulerian, direct Eulerian, and unstructured adaptive mesh refinement. A deliverable for this milestone is collaboration with MARBL for improved performance based on partially-assembled operators. The Urban exascale project is planning to couple real-time data with high performance simulations to study particulate transport in cities. Nek5000 will be used for large eddy and RANS (Reynolds-averaged Navier-Stokes) simulations of turbulent transport in these domains. Areas where the CEED team is aiding the Urban Canyon group include creating high-order meshes for urban geometries and developing boundary conditions appropriate for simulations close to real-life urban conditions. The software artifacts delivered as part of this milestone include GPU-enabled versions of the Nekbone and Laghos miniapps, which are provided through the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q4 of FY17, including: the CEED first annual meeting and participation of CEED researchers in a variety of outreach activities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Perpendicular Diffusion Coefficient of Comic Rays: The Presence of Weak Adiabatic Focusing

The influence of adiabatic focusing on particle diffusion is an important topic in astrophysics and plasma physics. In the past, several authors have explored the influence of along-field adiabatic focusing on the parallel diffusion of charged energetic particles. In this paper, using the unified nonlinear transport theory developed by Shalchi and the method of He and Schlickeiser, we derive a new nonlinear perpendicular diffusion coefficient for a non-uniform background magnetic field. This formula demonstrates that the particle perpendicular diffusion coefficient is modified by along-field adiabatic focusing. For isotropic pitch-angle scattering and the weak adiabatic focusing limit, the derived perpendicular diffusion coefficient is independent of the sign of adiabatic focusing characteristic length. For the two-component model, we simplify the perpendicular diffusion coefficient up to the second order of the power series of the adiabatic focusing characteristic quantity. We find that the first-order modifying factor is equal to zero and that the sign of the second order is determined by the energy of the particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Simple Analytical Forms of the Perpendicular Diffusion Coefficient for Two-component Turbulence. III. Damping Model of Dynamical Turbulence

In several astrophysical applications one needs analytical forms of cosmic-ray diffusion parameters. Some examples are studies of diffusive shock acceleration and solar modulation. In the current article we explore perpendicular diffusion based on the unified nonlinear transport theory. While we focused on magnetostatic turbulence in Paper I, we included the effect of dynamical turbulence in Paper II of the series. In the latter paper we assumed that the temporal correlation time does not depend on the wavenumber. More realistic models have been proposed in the past, such as the so-called damping model of dynamical turbulence. In the present paper we derive analytical forms for the perpendicular diffusion coefficient of energetic particles in two-component turbulence for this type of time-dependent turbulence. We present new formulas for the perpendicular diffusion coefficient and we derive a condition for which the magnetostatic result is recovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Discontinuous Finite Element Quasidiffusion Methods

Here in this paper, two-level methods for solving transport problems in one-dimensional slab geometry based on the quasi-diffusion (QD) method are developed. A linear discontinuous finite element method (LDFEM) is derived for the spatial discretization of the low-order QD (LOQD) equations. It involves special interface conditions at the cell edges based on the idea of QD boundary conditions (BCs). We consider different kinds of QD BCs to formulate the necessary cell-interface conditions. We develop two-level methods with independent discretization of the high-order transport equation and LOQD equations, where the transport equation is discretized using the method of characteristics and the LDFEM is applied to the LOQD equations. We also formulate closures that lead to the discretization consistent with a LDFEM discretization of the transport equation. The proposed methods are studied by means of test problems formulated with the method of manufactured solutions. Numerical experiments are presented demonstrating the performance of the proposed methods. Lastly, we also show that the method with independent discretization has the asymptotic diffusion limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Overcoming the Intrinsic Limitations of Fast Charging Lithium‐Ion Batteries Using Integrated Acoustic Streaming

A lithium‐ion battery's maximum charge rate and energy density are intrinsically limited by the ion diffusion rate in the electrolyte. Most research focuses on materials science solutions to this problem, with gradual improvement over the years. A mechanical solution is proposed to integrate an MHz‐order frequency surface acoustic wave (SAW) device into an existing 1.8 Ah multilayered Li‐ion pouch cell to enhance the ion diffusion rate and the overall battery performance. Both the charging rate and cycling lifetime are improved from SAW. At a 6C (10 min) charge and C/3 discharge rate, typical of electric vehicle applications, integrating SAW into the Li‐ion cell doubles the energy density and maintains at least 72% of the battery's initial capacity after 2000 cycles. Moreover, using SAW quantifiably reduces battery degradation in these conditions as determined by optical imaging, scanning electron microscopy, X‐ray diffraction, and neutron diffraction. The use of SAW appears to offer a method to avoid undesirable Li metal plating on the graphite anode during charging, and leads to a much longer battery lifetime and good charge capacity, all despite rapid charging.

Huang, An↗
DOE OSTI2018

Symmetry and separability of the neutron diffusion equation

Separation of variables is one of the oldest techniques for solving certain classes of partial differential equations (PDEs). As is the case with many other solution techniques for differential equations, separation of variables may be codified within the broader framework of symmetry analysis. Though the separation of variables technique is frequently used in the nuclear engineering context with various equations describing neutron transport, its connection to the symmetries of those equations has not yet been thoroughly established. It is thus the purpose of this work to establish that connection using neutron diffusion as both an initial step toward analysis of more generally applicable equations, and as a connection to previous results in related problems. Using Lie group analysis, it is found that the traditional space-time separable solution of the neutron diffusion equation (featuring a single α-eigenvalue) corresponds to time translation and flux scaling symmetries. Additional solutions of this equation are also constructed using its broader symmetry set.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Beyond the universal Dyson singularity for 1-D chains with hopping disorder

We study a simple non-interacting nearest neighbor tight-binding model in one dimension with disorder, where the hopping terms are chosen randomly. This model exhibits a well-known singularity at the band center both in the density of states and localization length. If the probability distribution of the hopping terms is well-behaved, then the singularities exhibit universal behavior, the functional form of which was first discovered by Freeman Dyson in the context of a chain of classical harmonic oscillators. We show here that this universal form can be violated in a tunable manner if the hopping elements are chosen from a divergent probability distribution. We also demonstrate a connection between a breakdown of universality in this quantum problem and an analogous scenario in the classical domain — that of random walks and diffusion with anomalous exponents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Emulators for Scarce and Noisy Data: Application to Auxiliary-Field Diffusion Monte Carlo for Neutron Matter

Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multimessenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as auxiliary-field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS, but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the effective field theory truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. In conclusion, this Letter will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

The arbitrary order mimetic finite difference method for a diffusion equation with a non-symmetric diffusion tensor

Here, we present the arbitrary order mimetic finite difference (MFD) discretization for the diffusion equation with non-symmetric tensorial diffusion coefficient in a mixed formulation on general polygonal meshes. The diffusion tensor is assumed to be positive definite. The asymmetry of the diffusion tensor requires changes to the standard MFD construction. We also present new approach for the construction that guarantees positive definiteness of the non-symmetric mass matrix in the space of discrete velocities. The numerically observed convergence rate for the scalar quantity matches the predicted one in the case of the lowest order mimetic scheme. For higher orders schemes, we observed super-convergence by one order for the scalar variable which is consistent with the previously published result for a symmetric diffusion tensor. The new scheme was also tested on a time-dependent problem modeling the Hall effect in the resistive magnetohydrodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗