Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Kinetic theory of phase space plateaux in a non-thermal energetic particle distribution

The transformation of kinetically unstable plasma eigenmodes into hole-clump pairs with temporally evolving carrier frequencies was recently attributed to the emergence of an intermediate stage in the mode evolution cycle, that of an unmodulated plateau in the phase space distribution of fast particles. The role of the plateau as the hole-clump breeding ground is further substantiated in this article via consideration of its linear and nonlinear stability in the presence of fast particle collisions and sources, which are known to affect the production rates and subsequent frequency sweeping of holes and clumps. In particular, collisional relaxation, as mediated by e.g. velocity space diffusion or even simple Krook-type collisions, is found to inhibit hole-clump generation and detachment from the plateau, as it should. On the other hand, slowing down of the fast particles turns out to have an asymmetrically destabilizing/stabilizing effect, which explains the well-known result that collisional drag enhances holes and their sweeping rates but suppresses clumps. It is further demonstrated that relaxation of the plateau edge gradients has only a minor quantitative effect and does not change the plateau stability qualitatively, unless the edge region extends far into the plateau shelf and the corresponding Landau pole needs to be taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Entropic and gradient flow formulations for nonlinear diffusion

Nonlinear diffusion ∂{sub t}ρ = Δ(Φ(ρ)) is considered for a class of nonlinearities Φ. It is shown that for suitable choices of Φ, an associated Lyapunov functional can be interpreted as thermodynamic entropy. This information is used to derive an associated metric, here called thermodynamic metric. The analysis is confined to nonlinear diffusion obtainable as hydrodynamic limit of a zero range process. The thermodynamic setting is linked to a large deviation principle for the underlying zero range process and the corresponding equation of fluctuating hydrodynamics. For the latter connections, the thermodynamic metric plays a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

DOUBLE-DIFFUSIVE MIXING IN STELLAR INTERIORS IN THE PRESENCE OF HORIZONTAL GRADIENTS

We have identified a potentially important source of mixing in stellar radiation zones which would arise whenever two conditions are satisfied: (1) the presence of an inverse vertical compositional gradient, and (2) the presence of density-compensating horizontal gradients of temperature and composition. The former can be caused by off-center nuclear burning, atomic diffusion, or surface accretion. The latter could be caused by rotation, tides, meridional flows, etc. The linear instability and its nonlinear development have been studied in an oceanographic context. It is known to drive the formation of stacks of fingering layers separated by diffusive interfaces, called intrusions. Using three-dimensional numerical simulations of the process in the astrophysically relevant region of parameter space, we find similar results and demonstrate that the material transport in the intrusive regime can be highly enhanced compared with pure diffusion, even in systems which would otherwise be stable to fingering (thermohaline) convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

STRAHL FORMATION IN THE SOLAR WIND ELECTRONS VIA WHISTLER INSTABILITY

This Letter puts forth a possible explanation for the formation of a solar wind strahl electron distribution function by means of local wave–particle interactions. A solar wind electron consists of a core and the hot “halo” electrons, which possess a net drift speed with respect to the core. According to the present model, pitch angle scattering of the initially isotropic drifting halo occurs when the enhanced whistler waves are excited by mildly anisotropic core electrons. The pitch angle scattering primarily affects the halo moving in the anti-sunward direction, resulting in pitch angle diffusion across the 90° gap through a nonlinear scattering process, and consequently leading to a reduction in the net drift speed of halo electrons. The remaining portion of the anti-sunward moving halo, which is not affected by pitch angle scattering, simply appears to form a field-aligned strahl in the electron velocity distribution. The present scenario of local generation of the strahl is demonstrated by the particle-in-cell simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2016

Review of mixing length estimates and effects of toroidicity in a fluid model for turbulent transport in tokamaks

Basic aspects of turbulent transport in toroidal magnetized plasmas are discussed. In particular the fluid closure has strong effects on zonal flows which are needed to create an absorbing boundary for long wave lengths and also to obtain the Dimits nonlinear upshift. The fluid resonance in the energy equation is found to be instrumental for generating the L–H transition, the spin-up of poloidal rotation in internal transport barriers, as well as the nonlinear Dimits upshift. The difference between the linearly fastest growing mode number and the corresponding longer nonlinear correlation length is also addressed. It is found that the Kadomtsev mixing length result is consistent with the non-Markovian diagonal limit of the transport at the nonlinearly obtained correlation length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear optical response of planar and spherical CdSe nanocrystals

The nonlinear optical response of a colloidal solution of planar CdSe semiconductor nanocrystals (nanoplatelets) is studied for the first time. The nonlinear optical response of these nanoparticles is compared to that of spherical CdSe nanocrystals (quantum dots). The photoinduced nonlinearity is attributed to the optical generation of long-lived charge carriers in the nanoobjects under study. It is shown that, upon the exposure of a cell with the solution of nanoparticles to focused continuous-wave (cw) laser radiation with a wavelength of 473 nm, the nonlinear optical responses of CdSe nanoplatelets and quantum dots are somewhat different at identical optical densities at the above-indicated wavelength. The differences are supposedly associated with a higher diffusion rate of spherical nanoparticles in the solution because of their smaller size compared to that of nanoplatelets.

36 MATERIALS SCIENCE↗
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

Save the Planet, Feed the Star: How Super-Earths Survive Migration and Drive Disk Accretion

Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse—and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity (“viscosity”). In inviscid, laminar disks, a planet readily repels gas away from its orbit. On short timescales, zero viscosity gas accumulates inside a planet’s orbit to slow Type I migration by orders of magnitude. On longer timescales, multiple super-Earths (distributed between, say, ∼0.1–10 au) can torque inviscid gas out of interplanetary space, either inward to feed their stars, or outward to be blown away in a wind. We explore this picture with 2D hydrodynamics simulations of Earths and super-Earths embedded in inviscid disks, confirming their slow/stalled migration even under gas-rich conditions, and showing that disk transport rates range up to ∼10{sup −7} M{sub ⊙} yr{sup −1} and scale as M-dot ∝ΣM{sub p}{sup 3/2}, where Σ is the disk surface density and M {sub p} is the planet mass. Gas initially sandwiched between two planets is torqued past both into the inner and outer disks. In sum, sufficiently compact systems of super-Earths can clear their natal disk gas in a dispersal history that may be complicated and non-steady but which conceivably leads over Myr timescales to large gas depletions similar to those characterizing transition disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model

We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatment of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2015

A new approach to enforcing discrete maximum principles in continuous Galerkin methods for convection-dominated transport equations

This work introduces a set of design principles and new algorithmic tools for enforcing maximum principles and/or positivity preservation in continuous finite element approximations to convection-dominated transport problems. Enabling a linear first-order advection equation as a model problem, we address the design of first-order artificial diffusion operators and their higherorder counterparts at the element matrix level. The proposed methodology leads to a nonlinear high-resolution scheme capable of resolving moving fronts and internal/boundary layers as sharp localized nonoscillatory features. The amount of numerical dissipation depends on the difference between the solution value at a given node and a local maximum or minimum. The shockcapturing numerical diffusion coefficient is designed to vanish as the nodal values approach a mass-weighted or linearity-preserving average. The universal applicability and simplicity of the element-based limiting procedure makes it an attractive alternative to edge-based algebraic flux correction.

Kuzmin, D.↗
DOE OSTI2016

Solving Drift-flux Model Using High-Resolution Spatial Discretization Scheme

Accurate modeling and simulation of the two-phase flow phenomena are critical to the system and safety analysis of light water reactors. Compared to the more complex two-fluid model, the drift-flux model treats the two phases as a mixture, and thus the models are formulated to consider the conservation laws of the mixture. The relative motion between the two phases is realized by constitutive correlations. Although the drift flux models have limitations in certain applications, they are widely used in many applications due to their simplicity and applicability to a wide range of two-phase flow problems. First-order numerical schemes are widely used in existing codes, which introduces excessive numerical diffusions. The objective of this work is to investigate a high-resolution spatial scheme in the application to solve drift-flux two-phase flow problems. In this study, a high-resolution spatial discretization scheme has been investigated in the application of solving drift-flux two-phase flow problem. The scheme is validated using existing experimental results. Compared to the first order scheme, it demonstrated excellent capability in reducing numerical error. Extended work will focus on the comparison of numerical errors from discretization schemes and model uncertainties. (authors)

42 ENGINEERING↗
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 OSTI2017

Variationally Derived Discontinuity Factors for the Asymptotic Homogenized Diffusion Equation

Here in this work, we derive and test variational discontinuity factors (DFs) for the asymptotic homogenized diffusion equation. We begin with a functional for optimally estimating the reactor multiplication factor, then introduce asymptotic expressions for the forward and adjoint angular fluxes, and finally require that all first-order error terms vanish. In this way, the reactor multiplication factor can be calculated with second-order error. The analysis leads to (1) an alternate derivation of the asymptotic homogenized diffusion equation, (2) variational boundary conditions for large periodic systems, and (3) variational DFs to be applied between adjacent periodic regions (e.g., fuel assemblies). Numerical tests show that applying the variational DFs to the asymptotic homogenized diffusion equation yields the most accurate estimates of the reactor multiplication factor compared to other DFs for a wide range of problems. However, the resulting assembly powers are less accurate than those obtained using other DFs for many realistic problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗