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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1984

Intensified silicon photodiode array detector linearity Application to coherent anti-Stokes Raman spectroscopy

The linearity of an intensified silicon photodiode array multichannel detector is studied with coherent anti-Stokes Raman and other similar signals. Studies with diffuse and focused (spherically and cylindrically) signals resolved apparent saturation problems which limit the dynamic range of the detector. In addition, it has been shown that there is no short-range wavelength (473-532-nm) dependence on this saturation. Theoretical explanations for these phenomena are also included.

Antcliff, R. R.↗
NASA NTRS1987

CONDIF - A modified central-difference scheme for convective flows

The paper presents a method, called CONDIF, which modifies the CDS (central-difference scheme) by introducing a controlled amount of numerical diffusion based on the local gradients. The numerical diffusion can be adjusted to be negligibly low for most problems. CONDIF results are significantly more accurate than those obtained from the hybrid scheme when the Peclet number is very high and the flow is at large angles to the grid.

Runchal, Akshai K.↗
NASA NTRS1993

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

The hypothesis is explored that the crust-mantle boundary of Venus is not in phase equilibrium but rather is rate-limited by the temperature-dependent volume diffusion of the slowest ionic species. The 1D thermal evolution problem is solved assuming that the mountains formed by uniform horizontal shortening of the crust and the lithospheric mantle at a constant rate. The time-dependent density structure and surface elevation are calculated by assuming a temperature-dependent reaction rate and local Airy isostatic compensation. For a horizontal strain rate of 10 exp -15/s or greater, the temperature increase at the base of the crust during mountain formation is modest to negligible, the deepening lower crust is metastable, and the surface elevation increases as the crust thickens. For strain rates less than 10 exp -16/s, crustal temperature increases with time because of internal heat production and the lower crust is more readily transformed to the dense eclogite assemblage. For such models, a maximum elevation is reached during crustal shortening.

Namiki, Noriyuki↗
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.↗
NASA NTRS1993

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗
NASA NTRS1993

Why does preferential diffusion strongly affect premixed turbulent combustion?

Combustion of premixed reactants in a turbulent flow is a classical but unresolved problem. The key problem is to explain the following data: the maximal turbulent and laminar burning velocities u(sub t) and u(sub L) occur at different equivalence ratios Phi. It is known that the equivalence ratio varies along a curved flame if molecular diffusivity D(sub fuel) does not equal D(sub oxygen). However, the mean flame radius of curvature is much larger than the laminar flame thickness delta-L. Therefore, significant influence of preferential diffusion should occur only if the flame propagation speed varies with flame curvature. This conclusion agrees with Zel'dovich's long-standing idea about the important role of leading points of a flame. The main objective of this paper is to prove Zel'dovich's hypothesis. An equation for the mean flame surface area density (MFSAD) is employed for this purpose. The second objective of this paper is to suggest a different approach to the derivation of the equation for MFSAD. It is based on the pdf equation for the reaction progress variable C and the relation between the pdf and MFSAD. This treatment suggests an entirely different closure assumption.

Kuznetsov, Vadim R.↗