Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Analysis of operator splitting errors for near-limit flame simulations

High-fidelity simulations of ignition, extinction and oscillatory combustion processes are of practical interest in a broad range of combustion applications. Splitting schemes, widely employed in reactive flow simulations, could fail for stiff reaction–diffusion systems exhibiting near-limit flame phenomena. The present work first employs a model perfectly stirred reactor (PSR) problem with an Arrhenius reaction term and a linear mixing term to study the effects of splitting errors on the near-limit combustion phenomena. Analysis shows that the errors induced by decoupling of the fractional steps may result in unphysical extinction or ignition. The analysis is then extended to the prediction of ignition, extinction and oscillatory combustion in unsteady PSRs of various fuel/air mixtures with a 9-species detailed mechanism for hydrogen oxidation and an 88-species skeletal mechanism for n-heptane oxidation, together with a Jacobian-based analysis for the time scales. The tested schemes include the Strang splitting, the balanced splitting, and a newly developed semi-implicit midpoint method. Results show that the semi-implicit midpoint method can accurately reproduce the dynamics of the near-limit flame phenomena and it is second-order accurate over a wide range of time step size. For the extinction and ignition processes, both the balanced splitting and midpoint method can yield accurate predictions, whereas the Strang splitting can lead to significant shifts on the ignition/extinction processes or even unphysical results. With an enriched H radical source in the inflow stream, a delay of the ignition process and the deviation on the equilibrium temperature are observed for the Strang splitting. On the contrary, the midpoint method that solves reaction and diffusion together matches the fully implicit accurate solution. The balanced splitting predicts the temperature rise correctly but with an over-predicted peak. For the sustainable and decaying oscillatory combustion from cool flames, both the Strang splitting and the midpoint method can successfully capture the dynamic behavior, whereas the balanced splitting scheme results in significant errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Surprisingly High Conductivity and Efficient Exciton Blocking in Fullerene/Wide-Energy-Gap Small Molecule Mixtures

We find that mixtures of C 60 with the wide energy gap, small molecular weight semiconductor bathophenanthroline (BPhen) exhibit a combination of surprisingly high electron conductivity and efficient exciton blocking when employed as buffer layers in organic photovoltaic cells. Photoluminescence quenching measurements show that a 1:1 BPhen/C 60 mixed layer has an exciton blocking efficiency of 84 ± 5% compared to that of 100% for a neat BPhen layer. This high blocking efficiency is accompanied by a 100-fold increase in electron conductivity compared with neat BPhen. Transient photocurrent measurements show that charge transport through a neat BPhen buffer is dispersive, in contrast to nondispersive transport in the compound buffer. Interestingly, although the conductivity is high, there is no clearly defined insulating-to-conducting phase transition with increased insulating BPhen fraction. Furthermore, we infer that C 60 undergoes nanoscale (<10 nm domain size) phase segregation even at very high (>80%) BPhen fractions.

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

Measurements of atoms and metastable species in N 2 and H 2 –N 2 nanosecond pulse plasmas

We report time-resolved, absolute number densities of metastable N 2 (A 3 Σ u + , v = 0, 1) molecules, ground state N 2 and H atoms, and rotational–translational temperature have been measured by tunable diode laser absorption spectroscopy and two-photon absorption laser-induced fluorescence in diffuse N 2 and N 2 –H 2 plasmas during and after a nanosecond pulse discharge burst. Comparison of the measurement results with the kinetic modeling predictions, specifically the significant reduction of the N 2 (A 3 Σ u + ) populations and the rate of N atom generation during the burst, suggests that these two trends are related. The slow N atom decay in the afterglow, on a time scale longer than the discharge burst, demonstrates that the latter trend is not affected by N atom recombination, diffusion to the walls, or convection with the flow. This leads to the conclusion that the energy pooling in collisions of N 2 (A 3 Σ u + ) molecules is a major channel of N 2 dissociation in electric discharges where a significant fraction of the input energy goes to electronic excitation of N 2 . Additional measurements in a 1% H 2 –N 2 mixture demonstrate a further significant reduction of N 2 (A 3 Σ u + , v = 0, 1) populations, due to the rapid quenching by H atoms accumulating in the plasma. Comparison with the modeling predictions suggests that the N 2 (A 3 Σ u + ) molecules may be initially formed in the highly vibrationally excited states. The reduction of the N 2 (A 3 Σ u + ) number density also diminishes the contribution of the energy pooling process into N2 dissociation, thus reducing the N atom number density. The rate of N atom generation during the burst also decreases, due to its strong coupling to N 2 (A 3 Σ u + , v) populations. On the other hand, the rate of H atom generation, produced predominantly by the dissociative quenching of the excited electronic states of N 2 by H 2 , remains about the same during the burst, resulting in a nearly linear rise in the H atom number density. Comparison of the kinetic model predictions with the experimental results suggests that the yield of H atoms during the quenching of the excited electronic state of N 2 by molecular H 2 is significantly less than 100%. The present results quantify the yield of N and H atoms in high-pressure H 2 –N 2 plasmas, which have significant potential for ammonia generation using plasma-assisted catalysis

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

SU-E-T-90: Concrete Forward-Scatter Fractions for Radiotherapy Shielding Applications

Purpose: There is little instruction within the primary shielding guidance document NCRP 151 for vault designs where the primary beam intercepts the maze. We have conducted a Monte-Carlo study to characterize forward-scattered radiation from concrete barriers with the intent of quantifying what amount of additional shielding outside the primary beam is needed in this situation. Methods: We reproduced our vault in MCNP 5 and simulated spectra obtained from the literature and from our treatment planning system for 10 and 18 MV beams. Neutron and gamma-capture contributions were not simulated. Energy deposited was scored at isocenter in a water phantom, within various cells that comprised the maze, and within cells that comprised the vault door. Tracks were flagged that scattered from within the maze to the door and their contributions were tallied separately. Three different concrete mixtures found in the literature were simulated. An empirically derived analytic equation was used for comparison, utilizing patient scatter fractions to approximate the scatter from concrete. Results: Our simulated data confirms that maze-scattered radiation is a significant contribution to total photon dose at the door. It contributes between 20-35% of the photon shielding workload. Forward-scatter fractions for concrete were somewhat dependent on concrete composition and the relative abundance of higher-Z elements. Scatter fractions were relatively insensitive to changes in the primary photon spectrum. Analytic results were of the same magnitude as simulated results. Conclusions: Forward-scattered radiation from the maze barrier needs to be included in the photon workload for shielding calculations in non-standard vault designs. Scatter fractions will vary with concrete composition, but should be insensitive to spectral changes between machine manufacturers. Further plans for investigation include refined scatter fractions for various concrete compositions, scatter fraction measurement, and simulation of photo-nuclear interaction within the maze barrier for high-energy beams.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗
NASA NTRS1980

Ar-40/Ar-39 dating, Ar diffusion properties, and cooling rate determinations of severely shocked chondrites

The ages of Ar-40/Ar-39 chondrites were computed to be 4.29 to 1.0 Gyr, with degassing times of 0.5 to 1.0 Gyr. The values of the diffusion parameter for Ar in Arrhenius plots show linear relationships which correspond to the degassing of different mineral phases with distinct K/Ca ratios and different average temperatures for Ar release. The experimental values of the diffusion parameter for the high-temperature phase of severely shocked chondrites are 10 to the -7th to 10 to the -5th/s for the shock-heating temperatures in the 950-1200 C range; the inferred reheating temperatures and the fraction of the Ar-40 loss during the reheating event suggest post-shock cooling rates and burial depths of 0.01-0.0001 C/s and 0.5-2m, respectively.

Bogard, D. D.↗
NASA NTRS1980

Solar wind sputtering effects in the atmospheres of Mars and Venus

It is found through an investigation, combining Monte Carlo simulations and analytical techniques, of the direct collisional interaction of an energetic particle flux with the neutral components of a planetary atmosphere, that solar wind sputtering could provide an important exospheric mass sink on both Mars and Venus under appropriate conditions. The computed rates of helium loss in the Venusian atmosphere and of carbon, nitrogen and oxygen in the Martian atmosphere imply that sputtering would have a significant effect on the noble gas budget of Venus and dominate the chemical and photochemical loss processes of Mars. Because of diffusive separation of lighter elements and isotopes, and because the gravitational binding energy is proportional to the mass, the erosion process preferentially removes the lighter components of the atmosphere. Solar wind sputtering could therefore compete with other erosion mechanisms in generating substantial fractionation effects.

Watson, C. C.↗
NASA NTRS1994

Papers presented to the Conference on Chondrules and the Protoplanetary Disk

The following topics are covered in the presented papers: (1) producing chondrules; (2) carbons, CAI's, and chondrules; (3) large scale processes in the solar nebula; (4) chondrule-matrix relationships in chondritic meteorites; (5) overview of nebula models; (6) constraints placed on the nature of chondrule precursors; (7) turbulent diffusion and concentration of chondrules in the protoplanetary nebula; (8) heating and cooling in the solar nebula; (9) crystallization trends of precursor pyroxene in ordinary chondrites; (10) precipitation induced vertical lightning in the protoplanetary nebula; (11) the role of chondrules in nebular fractionations of volatiles and other elements; (12) astronomical observations of phenomena in disks; (13) experimental constraints on models for origins of chondrules, and various other topics.

Source record↗
NASA NTRS1994

Weakly ionized cosmic gas: Ionization and characterization

Since collective plasma behavior may determine important transport processes (e.g., plasma diffusion across a magnetic field) in certain cosmic environments, it is important to delineate the parameter space in which weakly ionized cosmic gases may be characterized as plasmas. In this short note, we do so. First, we use values for the ionization fraction given in the literature, wherein the ionization is generally assumed to be due primarily to ionization by cosmic rays. We also discuss an additional mechanism for ionization in such environments, namely, the photoelectric emission of electrons from cosmic dust grains in an interstellar Far Ultra Violet (FUV) radiation field. Simple estimates suggest that under certain conditions this mechanism may dominate cosmic ray ionization, and possibly also the photoionization of metal atoms by the interstellar FUV field, and thereby lead to an enhanced ionization level.

Rosenberg, M.↗
NASA NTRS1994

Ion injection and acceleration at quasi-perpendicular shocks

We present and discuss results of a new model for ion injections and acceleration at quasi-perpendicular collisionless shocks. We use the one-dimensional hybrid simulation (kinetic ions/fluid electrons) and impose an assumption on the ion motion so that diffusion across the magnetic field (normal to the shock front) is possible. These motions are otherwise suppressed by both one- and two-dimensional simulations. We find that, even in strictly perpendicular shocks, when scattering normal to the field is included, a fraction of the incident ions are accelerated to suprathermal energies. When reasonable scattering times are considered, only pickup ions are injected, whereas thermal solar wind ions are not. The acceleration of these ions is very rapid. We have found that a few of the initially low-energy pickup ions can reach many tens to a few hundred times the plasma ramming energy in less than 100 gyroperiods. Furthermore, highly field aligned energetic ions are found to exist upstream of the slightly oblique shock. The most direct applications of this study is toward the interpretation of observations of both solar wind and interstellar pickup ion distributions in the vicinity of interplanetary shocks which are most often quasi-perpendicular. This work also directly addresses a fundamental issue with regard to our current understanding of the anomalous component of cosmic rays.

Giacalone, J.↗
DOE OSTI2020

Bio-inspired redox mediated electrolyte for high performance flexible supercapacitor applications over broad temperature domain

For practical supercapacitors, recent advances in electrolytes have diminished long-standing problems such chemical stability, low ionic conductivity and limited electrochemical performance. Herewith, we communicate a facile method to produce anhydrous redox-active Bio-inspired electrolytes for flexible supercapacitors with application over a broad temperature range. The glycerol (Gly) gel electrolyte consists of various fractions of potassium hydroxide (KOH) and redox active unit, i.e., ammonium molybdate, (NH 4 ) 2 MoO 4 (abbreviated as Mo), in the anhydrous state. The doped glycerol matrix acts as a framework that provides a stable gel electrolyte, forming hierarchical pathway for ion diffusion. The developed gel network with Gly-KOH-MoX (X = 3, 5, 10, and 15%, w/w) provides an effective route to carry the ions onto the pores of the electrodes. The supercapacitor with Gly-KOH5-Mo10 gel provides a specific capacitance of 328 F g –1 , and the same device achieved an energy density of 45.6 Wh kg –1 at a power density of 497 W kg –1 . In addition, the supercapacitor maintained extraordinary cyclic robustness and wide temperature tolerance. Thus, the low cost and environmentally benign, cyclic durability of Gly-KOH-Mo makes these gel electrolytes a potential candidate for use in various energy storage devices.

25 ENERGY STORAGE↗
DOE OSTI2015

A spectral analysis of the domain decomposed Monte Carlo method for linear systems

The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear oper- ator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approxi- mation and the mean chord approximation are applied to estimate the leakage frac- tion of random walks from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem in numerical experiments to test the models for symmetric operators with spectral qualities similar to light water reactor problems. We find, in general, the derived approximations show good agreement with random walk lengths and leakage fractions computed by the numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Where is the cosmic-ray modulation boundary of the heliosphere?

The intensity of Galactic cosmic rays in the heliosphere is modulated by solar activities. The outer boundary where the solar modulation begins has always been a subject matter of debate in the cosmic-ray and heliophysics community. Various experimental methods and theoretical model calculations have been used to determine the boundary. Although the heliopause was always suspected to be the boundary, it is only until very recently after Voyager 1 had crossed the heliopause did we confirm that the boundary is indeed the heliopause. In this paper, we use a model simulation and detailed Voyager observation of cosmic rays at the heliopause crossing to show that the modulation boundary, in fact, is a fraction of an AU beyond the heliopause. Such a conclusion requires a very low turbulence level of the interstellar magnetic field in the outer heliosheath. According to the quasi-linear theory, a low level of turbulence should result in a very large diffusion coefficient parallel to the magnetic field and a very small perpendicular diffusion coefficient. For the first time, we are confident that Voyager 1 has obtained the truly pristine local interstellar cosmic-ray spectra down to the energies below 1 MeV. The cosmic-ray intensity is rapidly filtered by a thin layer of the interstellar magnetic field immediately outside of the heliopause. Its filtration amount depends on the conditions of magnetic field turbulence on the both sides of the heliopause, thus making it solar-cycle dependent as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantifying the hydrodynamic contribution to electrical transport in non-Brownian suspensions

Electrical transport in semiconducting and metallic particle suspensions is an enabling feature of emerging grid-scale battery technologies. Although the physics of the transport process plays a key role in these technologies, no universal framework has yet emerged. Here, we examine the important contribution of shear flow to the electrical transport of non-Brownian suspensions. We find that these suspensions exhibit a strong dependence of the transport rate on the particle volume fraction and applied shear rate, which enables the conductivity to be dynamically changed by over 10 7 decades based on the applied shear rate. We combine experiments and simulations to conclude that the transport process relies on a combination of charge and particle diffusion with a rate that can be predicted using a quantitative physical model that incorporates the self-diffusion of the particles.

42 ENGINEERING↗
NASA NTRS2002

Radio Sources Associated with Intermediate X-ray Luminosity Objects in Merging Galaxy Systems

We present new, high-resolution 6, 3.6, and 2 cm radio images of a time-ordered sequence of merging galaxy systems. The new data have a resolution of less than 100pc and a sensitivity comparable to a few x Cas A. We detect compact radio sources in all systems, generally embedded in more diffuse radio emission at the longer wavelengths. Several of the compact radio sources are coincident with compact Intermediate-luminosity X-ray Objects (IXOs) in these systems, and many more are within the 3$/sigma$ Chandra position errors for other IXOs. The fraction of radio identifications and the nature of the radio sources changes as a function of merger stage. These data suggest that the IXOs are associated with complexes of supernova remnants, and therefore with star formation that has occurred within the last $/sim$10$circumflex7$ yr, but are not located in HII regions where copious star formation is occurring currently.

Neff, S. G.↗
DOE OSTI2013

Effect of magnetic fluctuations on the confinement and dynamics of runaway electrons in the HT-7 tokamak

The nature of runaway electrons is such that the confinement and dynamics of the electrons can be strongly affected by magnetic fluctuations in plasma. Experimental results in the HT-7 tokamak indicated significant losses of runaway electrons due to magnetic fluctuations, but the loss processes did not only rely on the fluctuation amplitude. Efficient radial runaway transport required that there were no more than small regions of the plasma volume in which there was very low transport of runaways. A radial runaway diffusion coefficient of D{sub r} Almost-Equal-To 10 m{sup 2}s{sup -1} was derived for the loss processes, and diffusion coefficient near the resonant magnetic surfaces and shielding factor #Greek Upsilon With Hook Symbol#=0.8 were deduced. Test particle equations were used to analyze the effect of magnetic fluctuations on runaway dynamics. It was found that the maximum energy that runaways can gain is very sensitive to the value of {alpha}{sub s} (i.e., the fraction of plasma volume with reduced transport). {alpha}{sub s}=(0.28-0.33) was found for the loss processes in the experiment, and maximum runaway energy could be controlled in the range of E=(4 MeV-6 MeV) in this case. Additionally, to control the maximum runaway energy below 5 MeV, the normalized electric field needed to be under a critical value D{sub {alpha}}=6.8, and the amplitude normalized magnetic fluctuations b(tilde sign) needed to be at least of the order of b(tilde sign) Almost-Equal-To 3 Multiplication-Sign 10{sup -5}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Thermodynamic and kinetic modeling of Mn-Ni-Si precipitates in low-Cu reactor pressure vessel steels

Formation of large volume fractions of Mn-Ni-Si precipitates (MNSPs) causes excess irradiation embrittlement of reactor pressure vessel (RPV) steels at high, extended-life fluences. Thus, a new and unique, semi-empirical cluster dynamics model was developed to study the evolution of MNSPs in low-Cu RPV steels. The model is based on CALPHAD thermodynamics and radiation enhanced diffusion kinetics. The thermodynamics dictates the compositional and temperature dependence of the free energy reductions that drive precipitation. The model treats both homogeneous and heterogeneous nucleation, where the latter occurs on cascade damage, like dislocation loops. The model has only four adjustable parameters that were fit to an atom probe tomography (APT) database. The model predictions are in semi-quantitative agreement with systematic Mn, Ni and Si composition variations in alloys characterized by APT, including a sensitivity to local tip-to-tip variations even in the same steel. The model predicts that heterogeneous nucleation plays a critical role in MNSP formation in lower alloy Ni contents. Single variable assessments of compositional effects show that Ni plays a dominant role, while even small variations in irradiation temperature can have a large effect on the MNSP evolution. Within typical RPV steel ranges, Mn and Si have smaller effects. Furthermore, the delayed but then rapid growth of MNSPs to large volume fractions at high fluence is well predicted by the model. For purposes of illustration, the effect of MNSPs on transition temperature shifts are presented based on well-established microstructure-property and property-property models.

36 MATERIALS SCIENCE↗