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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

High-velocity interstellar clouds

Interstellar absorption line studies have revealed that rapidly moving clouds within one kiloparsec of the sun have considerably less depletion of heavy elements than in the low-velocity diffuse clouds, and their density is somewhat greater and their thickness smaller. Evaporation of interstellar grains, resulting from grain-grain collisions in interstellar shocks, may provide a possible explanation for the observations. Some problems involved in this mechanism are discussed, particularly the likely dominant importance of the interstellar magnetic field in confining the grains within the shocked layer and in accelerating them through the associated electric field, thus offsetting the frictional force of the gas.

Spitzer, L.↗
NASA NTRS1981

An extended classical solution of the droplet growth problem

Problems of applying the classical kinetic theory to the growth of small droplets from vapor are examined. A solution for the droplet growth equation is derived which is based on the assumption of a diffusive field extending to the drop surface. The method accounts for partial thermal and mass accommodation at the interface and the kinetic limit to the mass and heat fluxes, and it avoids introducing the artifact of a discontinuity in the thermal and vapor field near the droplet. Consideration of the environmental fields in spherical geometry utilizing directional fluxes yields boundary values in terms of known parameters and a new Laplace transform integral.

Anderson, B. J.↗
NASA NTRS1984

Combustor flame flashback

Flashback, a problem that occurs in premixed-prevaporized combustors, is the upstream propagation of the flame from the combustor into the premixing tubes. Not only does flashback change the combustion process from premixed burning to diffusion burning, thus creating more pollutants, but it also inflicts considerable damage to the fuel injector, premixing tube and other equipment upstream. The conditions at which flashback occurs in steady burning and the mechanism that causes flashback in both steady and transient flow are studied. The equivalence ratio at which flashback occurs is being measured for inlet temperatures of 600-950 K, premixer wall temperatures of 450-1050 K and premixer velocities of 40-80 ft/s. These data are presented.

Proctor, M. P.↗
NASA NTRS1987

A geostationary imaging spectrometer TOMS instrument

One design for a geostationary Total Ozone Mapping Spectrometer (TOMS) with many desirable features is an imaging spectrometer. A preliminary study makes use of a 0.25 m Czerny-Turner spectrometer with which the Earth is imaged on a charge-coupled device (CCD) in dispersed light. The wavelength is determined by a movable grating which can be set arbitrarily by ground control. The signal integration time depends on wavelength but this system allows arbitrary timing by command. Special circumstances such as a requirement to track a low-lying sulfur dioxide cloud or a need to discriminate high level ozone from total ozone at midlatitudes could be obtained by adding a particular wavelength to the normally pre-programmed time sequence. The incident solar irradiance is measured by deploying a diffuser plate in the field of view. Individual detector elements correspond to scene elements in which the several wavelengths are serially sampled and the Earth radiance is compared to the incident sunlight. Thus the problem of uncorrelated drift of multiple detectors is removed.

Krueger, Arlin J.↗
NASA NTRS1987

A parallelized elliptic solver for reacting flows

A modified Newton algorithm for the solution of nonlinear elliptic boundary value problems via finite discretization methods is presented. A serial implementation of this algorithm which has recently been applied successfully to the computation of an axisymmetric over-ventilated subsonic laminar methane-air jet diffusion flame is described. Parallel implementation issues and a complexity theory are presented. Included as well are actual performance data for model systems obtained on the Intel Hypercube and a discussion of its implications for modeling realistic systems.

Keyes, David E.↗
NASA NTRS1996

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗
NASA NTRS2000

The AstroBiology Explorer (ABE) MIDEX Mission

The Astrobiology Explorer (ABE) is a Medium-Class Explorer (MIDEX) mission concept currently under study at NASA's Ames Research Center. ABE will conduct infrared (IR) spectroscopic observations with much better sensitivity than Infrared Space Observatory (ISO) or the Stratospheric Observatory for Infrared Astronomy program (SOFIA) in order to address outstanding astrobiologically important problems in astrochemistry as well as important astrophysical investigations. The core observational astrobiology program would make fundamental scientific progress in understanding the cosmic history of molecular carbon, the distribution of organic matter in the diffuse interstellar medium, tracing the chemical history of complex organic molecules in the interstellar medium, and the evolution of organic ices in young planetary systems. The ABE instrument concept includes a 0.5 m aperture Cassegrain telescope and a suite of three moderate resolution (R = 1000 - 4000) spectrographs which cover the entire lambda = 2.5-20 micron spectral region. Use of large format (1024 x 1024 pixel or larger) IR detector arrays will allow each spectrograph to cover an entire octave of spectral range per exposure without any moving parts. The telescope is passively cooled by a sun shade to below 65 K, and the detectors are cooled with solid H2 cryogen to approximately 8 K. ABE will be placed in an Earth-trailing one AU solar orbit by a Delta II launch vehicle. This energetically favorable orbit provides a low thermal background, affords good access to the entire sky over the one year mission lifetime, and allows adequate communications bandwidth. The spacecraft will be stabilized in three axes and will be pointed to an accuracy of approximately one arcsecond at ABE's several thousand individual scientific targets.

Greene, Thomas↗
NASA NTRS2004

The Distribution of Cosmic-Ray Sources in the Galaxy, Gamma-Rays and the Gradient in the CO-to-H2 Relation

We present a solution to the apparent discrepancy between the radial gradient in the diffuse Galactic gamma-ray emissivity and the distribution of supernova remnants, believed to be the sources of cosmic rays. Recent determinations of the pulsar distribution have made the discrepancy even more apparent. The problem is shown to be plausibly solved by a variation in the Wco-to-N(H2) scaling factor. If this factor increases by a factor of 5-10 from the inner to the outer Galaxy, as expected from the Galactic metallicity gradient and supported by other evidence, we show that the source distribution required to match the radial gradient of gamma-rays can be reconciled with the distribution of supernova remnants as traced by current studies of pulsars. The resulting model fits the EGRET gamma-ray profiles extremely well in longitude, and reproduces the mid-latitude inner Galaxy intensities better than previous models.

Strong, A. W.↗
NASA NTRS2004

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗
DOE OSTI2018

Results of Robotic Evaluation of Uranium-235 in Gaseous Diffusion Piping Holdup Deposits - 18303

Measurement of U-235 in pipes within old uranium enrichment facilities is time-consuming and very costly. Radiometric assay from outside the pipes suffers from manual deployment challenges, attenuation of detection through pipe walls, long counting times, approximate modeling, and shortfalls associated with transcription and human interpretation. The D and D process of cutting into these pipes creates a unique possibility of robotic in-pipe assay to address these problems. This paper presents the results of testing a first-of-kind radiometric non-destructive assay (NDA) robot for high-cadence, accurate, certain, paperless assay of U-235 in holdup deposits within process piping. The robot's enabling capability is an innovative radiometric collimation method that views only a one-foot moving annulus of pipe wall. This exploits the internal axisymmetry of pipes to provide a direct measurement of industry-standard grams per foot U-235. The initial proof-of-application robot was successfully hot tested by operators at the Gaseous Diffusion Enrichment Facility in Piketon, Ohio, in late 2017. This testing achieved all goals of traversing pipe, measuring U-235, deployment by facility personnel, and displaying analyzed radiation, visual, and geometric information. The success was a Kitty Hawk moment. The system demonstrated robotic in-pipe measurement of holdup deposit to a resolution, accuracy, and speed unachievable by alternate methods. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1986

Observational discrimination between modes of shock propagation in interstellar clouds: Predictions of CH+ and SH+ column densities in diffuse clouds

Considerable effort in recent years has been devoted to the study of shocks in the diffuse interstellar medium. This work has been motivated partly by the observations of rotationally excited states of H2, and partly by the realization that species such as CH(+), OH and H2O might be formed preferentially in hot, post-shock gas. The problem of CH(+) and the difficulties encountered when trying to explain the high column densities, observed along lines of sight to certain hot stars, have been reviewed earlier. The importance of a transverse magnetic field on the structure of an interstellar shock was also demonstrated earlier. Transverse magnetic fields above a critical strength give rise to an acceleration zone or precursor, in which the parameters on the flow vary continuously. Chemical reactions, which change the degree of ionization of the gas, also modify the structure of the shock considerably. Recent work has shown that large column densities of CH(+) can be produced in magnetohydrodynamic shock models. Shock speeds U sub s approx. = 10 km/s and initial magnetic field strengths of a few micro G are sufficient to produce ion-neutral drift velocities which can drive the endothermic C(+)(H2,H)CH(+) reaction. It was also shown that single-fluid hydrodynamic models do not generate sufficiently large column densities of CH(+) unless unacceptably high shock velocities (u sub s approx. 20 km/s) are assumed in the models. Thus, the observed column densities of CH(+) provide a constraint on the mode of shock propagation in diffuse clouds. More precisely, they determine a lower limit to the ion-neutral drift velocity.

Flower, D. R.↗
DOE OSTI2021

Atomistic modeling of meso -timescale processes with $\mathrm{SEAKMC}$: A perspective and recent developments

On-the-fly kinetic Monte Carlo (kMC) methods have recently garnered significant attentions after successful applications to various atomic-scale problems using a timescale outside the reach of classical molecular dynamics. These methods play a critical role in modeling atomistic meso-timescale processes, and it is therefore essential to further improve their capabilities. In this report we review one of the on-the-fly kMC methods, Self-Evolving Atomistic kinetic Monte Carlo (SEAKMC) and propose two schemes that considerably enhance the efficiency of saddle point searches (SPSs) during the simulations. The performance of these schemes is tested using the diffusion of point defects in bcc Fe. In addition, we discuss approaches to significantly mitigate limitations of these schemes, which further improves their efficiencies. Importantly, these schemes improve the SPS efficiency not only for SEAKMC but also for other on-the-fly kMC methods, broadening the applications of on-the-fly kMC simulations to complex meso-timescale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Analytical Solutions of Transient Drift-Diffusion in P-N Junction Pixel Sensors

Radiation detection in applications ranging from high energy physics to medical imaging rely on solid state detectors, often hybrid pixel detectors with (1) reverse biased p-n junction pixel sensors and (2) readout ASICs, attached by flip chip-bonding. Transient signals characteristics are important in, e.g., matching ASIC and sensor design, modeling and optimizing detector parameters and describing timing and charge sharing properties. Currently analytical forms of transient signals are available for only a few limited cases (e.g., drift or diffusion) or for the steady state (which is not relevant for high energy radiation detection). Tools are available for (relatively slow) numerical evaluation of the transient charge transport. We present here the first analytical solutions of partial differential equations de scribing drift-diffusion-recombination charge transport in planar p-n junction sensors in a variety of conditions: (1) undepleted, (2) fully depleted, (3) taking into account the gradual velocity saturation, and (4) overdepleted. We deduce the Green’s functions which can be applied to any detection problem through simple convolution with the initial conditions. We compare the analytical solutions with Monte Carlo simulations and industry standard simulations (Synopsys Sentaurus), demonstrating good agreement. Using the analytical equations enables fast modeling of the influence of various detector parameters on tracking, imaging and timing performance, describing performance and enabling optimizations for different applications. Finally, we illustrate this model with applications in 3D+T (x,y,z,time) photon tracking and 4D+T (x,y,θ,φ,time) relativistic charged particle tracking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2022

Grassmannian Diffusion Maps--Based Dimension Reduction and Classification for High-Dimensional Data

This work introduces the Grassmannian diffusion maps (GDMaps), a novel nonlinear dimensionality reduction technique that defines the affinity between points through their representation as low-dimensional subspaces corresponding to points on the Grassmann manifold. Here, the method is designed for applications, such as image recognition and data-based classification of constrained high-dimensional data where each data point itself is a high-dimensional object (i.e., a large matrix) that can be compactly represented in a lower-dimensional subspace. The GDMaps is composed of two stages. The first is a pointwise linear dimensionality reduction wherein each high-dimensional object is mapped onto the Grassmann manifold representing the low-dimensional subspace on which it resides. The second stage is a multipoint nonlinear kernel-based dimension reduction using diffusion maps to identify the subspace structure of the points on the Grassmann manifold. To this end, an appropriate Grassmannian kernel is used to construct the transition matrix of a random walk on a graph connecting points on the Grassmann manifold. Spectral analysis of the transition matrix yields low-dimensional Grassmannian diffusion coordinates embedding the data into a low-dimensional reproducing kernel Hilbert space. Further, a novel data classification/recognition technique is developed based on the construction of an overcomplete dictionary of reduced dimension whose atoms are given by the Grassmannian diffusion coordinates. Three examples are considered. First, a "toy" example shows that the GDMaps can identify an appropriate parametrization of structured points on the unit sphere. The second example demonstrates the ability of the GDMaps to revealing the intrinsic subspace structure of high-dimensional random field data. In the last ex- ample, a face recognition problem is solved considering face images subject to varying illumination conditions, changes in face expressions, and occurrence of occlusions. The technique presented high recognition rates (i.e., 95% in the best case) using a fraction of the data required by conventional methods.

42 ENGINEERING↗
DOE OSTI2013

A multiphase electrokinetic flow model for electrolytes with liquid/liquid interfaces

A numerical model for electrokinetic flow of multiphase systems with deformable interfaces is presented, based on a combined level set-volume of fluid technique. A new feature is a multiphase formulation of the Nernst–Planck transport equation for advection, diffusion and conduction of individual charge carrier species that ensures their conservation in each fluid phase. The numerical model is validated against the analytical results of Zholkovskij et al. (2002) [1], and results for the problem of two drops coalescing in the presence of mobile charge carriers are presented. The time taken for two drops containing ions to coalesce decreases with increasing ion concentration.

42 ENGINEERING↗
DOE OSTI2015

Path optimization with limited sensing ability

We propose a computational strategy to find the optimal path for a mobile sensor with limited coverage to traverse a cluttered region. The goal is to find one of the shortest feasible paths to achieve the complete scan of the environment. We pose the problem in the level set framework, and first consider a related question of placing multiple stationary sensors to obtain the full surveillance of the environment. By connecting the stationary locations using the nearest neighbor strategy, we form the initial guess for the path planning problem of the mobile sensor. Then the path is optimized by reducing its length, via solving a system of ordinary differential equations (ODEs), while maintaining the complete scan of the environment. Furthermore, we use intermittent diffusion, which converts the ODEs into stochastic differential equations (SDEs), to find an optimal path whose length is globally minimal. To improve the computation efficiency, we introduce two techniques, one to remove redundant connecting points to reduce the dimension of the system, and the other to deal with the entangled path so the solution can escape the local traps. Numerical examples are shown to illustrate the effectiveness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗