Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Potential–Dependent BDAC Adsorption on Zinc Enabling Selective Suppression of Zinc Corrosion for Energy Storage Applications

Utility-scale zinc (Zn) batteries are a promising solution to address the problem of intermittency of renewable energy sources; however, Zn-metal anodes in these batteries suffer from capacity loss due to spontaneous corrosion of the Zn especially when high-surface area anode configurations are employed. Additionally, Zn dendrites are known to form during battery charging limiting the cycle-life of these batteries. Electrolyte additives have been explored that prevent aforementioned issues, but these too come at a cost, i.e., surface-blocking additives polarize the electrode surface leading to loss in the voltaic and energy efficiencies of the battery. In this contribution, a novel electrolyte additive, benzyldimethylhexadecylammonium chloride (BDAC), is investigated for its ability to suppresses corrosion of Zn in an acidic (pH = 3) electrolyte. An attribute of BDAC distinct from previously studied additives is that it selectively suppresses electrochemical activity when the Zn electrode is at its corrosion potential; however, during high-rate Zn deposition (charging) or stripping (discharging), BDAC is essentially deactivated and thus it does not appreciably polarize the electrode surface, thus minimizing voltaic efficiency losses. This selective corrosion suppression behavior is explored using slow-scan voltammetry, which reveals hysteresis implying a potential- or current-dependent BDAC adsorption mechanism in which BDAC reaches higher surface coverages when the partial currents at the Zn surface are low (e.g., at or near the corrosion potential), but BDAC coverage is reduced considerably when the Zn deposition or stripping rates are increased. Numerical simulations of the BDAC diffusion-adsorption process corroborate this mechanism. Ramifications of our approach to the selective suppression of Zn dendrites are discussed.

25 ENERGY STORAGE↗
DOE OSTI2024

M3SF-24LL010302052-Comprehensive Analysis of Radionuclide Interaction with Hydrothermally Altered Repository Materials

This progress report (Level 3 Milestone Number M3SF-24LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Work Package Number SF-24LL01030205. The research is focused on actinide and radionuclide sequestration in hydrothermally altered repository materials. In FY24, we completed a rigorous analysis of Se sorption to iron oxide phases using our L-SCIE sorption database. This effort explicitly accounts for surface titration behavior of oxide surfaces that was the subject of a recently published manuscript (Han et al., 2023). With this in mind, our Se sorption analysis now yields a more robust workflow for developing self-consistent surface complexation modeling approaches that can be adapted to specific SCM conceptual and numerical approaches (i.e. non-electrostatic, diffuse layer, triple layer models). In late FY24, we plan to publish the results of our comprehensive surface complexation modeling of Se(IV) and Se(VI) sorption to iron oxide mineral phases and provide a path forward to developing robust radionuclide sorption models for use in performance assessment. In FY24, we also submitted a manuscript summarizing our approach to integrating radionuclide sorption and coprecipitation phenomena and evaluation of radionuclide partitioning values across a range of radionuclides relevant to performance assessment. We demonstrated our approach in detail using Se sorption and coprecipitation with iron oxide minerals as a test case. This manuscript was recently accepted for publication in Applied Geochemistry (Balboni et al., Accepted). We also continued experiments to identify radionuclide interaction with hydrothermally altered crystalline repository and backfill materials. Recent research performed at Los Alamos National Laboratory (LANL) and Sandia National Laboratory (SNL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are now examining how mineral alteration affects retardation behavior of a suite of radionuclides of interest to repository performance assessment. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database, developed using our L-SCIE approach, can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. A short update of results to date is presented below.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗
NASA NTRS1987

Theoretical approach to oxygen atom degradation of silver

Based on available Rutherford backscattering spectrometry (RBS), proton induced X-ray emission (PIXE) and ellipsometry data obtained on silver specimens subjected to atomic oxygen attack in low Earth orbit STS flight 41-G, a theory was developed to model the oxygen atom degradation of silver. The diffusion of atomic oxygen in a microscopically nonuniform medium is an essential constituent of the theory. The driving force for diffusion is the macroscopic electrochemical potential gradient developed between the specimen surface exposed to the ambient and the bulk of the silver specimen. The longitudinal electric effect developed parallel to the gradient is modified by space charge of the diffusing charged species. Lateral electric fields and concentration differences also exist due to the nonuniform nature of the medium. The lateral concentration differences are found to be more important than the lateral electric fields in modifying the diffusion rate. The model was evaluated numerically. Qualitative agreement exists between the kinetics predicted by the theory and kinetic data taken in ground-based experiments utilizing a plasma asher.

Fromhold, Albert T., Jr.↗
NASA NTRS1987

Computation of diffuse sky irradiance from multidirectional radiance measurements

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Ahmad, Suraiya P.↗
NASA NTRS1988

Turbulent diffusivity

A method for determining the turbulent diffusivity tensor is proposed which has application to problems such as the mixing in the interior of the sun for the solar neutrino problem, the chemical evolution of the galaxy, and grain sedimentation on the primitive solar nebula. A model for determining the turbulent energy spectral function, the eddy decorrelation rate, and the growth rate of the underlying instability is discussed. Good agreement is obtained between the present values for the coefficient of turbulent diffusivity and those obtained from both numerical simulation studies and laboratory data on shear turbulence.

Canuto, V. M.↗
NASA NTRS1990

Process modelling for Space Station experiments

Examined here is the sensitivity of a variety of space experiments to residual accelerations. In all the cases discussed the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion, or fluid motion must be suppressed or eliminated. In the latter case, the experiments are sensitive to steady and low frequency accelerations. For experiments such as the directional solidification of melts with two or more components, determination of the velocity response alone is insufficient to assess the sensitivity. The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt-crystal interface. As far as the response to transient disturbances is concerned, the sensitivity is determined by both the magnitude and frequency of the acceleration and the characteristic momentum and solute diffusion times. The microgravity environment, a numerical analysis of low gravity tolerance of the Bridgman-Stockbarger technique, and modeling crystal growth by physical vapor transport in closed ampoules are discussed.

Alexander, J. Iwan D.↗
NASA NTRS1996

The Effects of Flame Structure on Extinction of CH4-O2-N2 Diffusion Flames

The effects of flame structure on the extinction limits of CH4-O2-N2 counterflow diffusion flames were investigated experimentally and numerically by varying the stoichiometric mixture fraction Z(sub st), Z(sub st) was varied by varying free-stream concentrations, while the adiabatic flame temperature T(sub ad) was held fixed by maintaining a fixed amount of nitrogen at the flame. Z(sub st) was varied between 0.055 (methane-air flame) and 0.78 (diluted- methane-oxygen flame). The experimental results yielded an extinction strain rate K(sub ext) of 375/s for the methane-air flame, increasing monotonically to 1042/s for the diluted-methane-oxygen flame. Numerical results with a 58-step Cl mechanism yielded 494/s and 1488/s, respectively. The increase in K(sub ext) with Z(sub st) for a fixed T(sub ad) is explained by the shift in the O2 profile toward the region of maximum temperature and the subsequent increase in rates for chain-branching reactions. The flame temperature at extinction reached a minimum at Z(sub st) = 0.65, where it was 200 C lower than that of the methane-air flame. This significant increase in resistance to extinction is seen to correspond to the condition in which the OH and O production zones are centered on the location of maximum temperature.

Du, J.↗
NASA NTRS2011

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗
NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗
DOE OSTI2022

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1980

Transient and diffusion analysis of Hg Cd Te

The solution of the one dimensional planar interface solidification problem is presented. Numerical solutions are applied to experimental solute profiles to determine the effective diffusion constant for Hg Co Te. Progress on the design and construction of a high gradient furnace is reported.

Clayton, J. C.↗
NASA NTRS1987

Numerical simulation of cometary nuclei. III - Internal temperatures of cometary nuclei

The thermal diffusion equation for the internal temperature of cometary nuclei is exactly solved by means of a one-dimensional numerical model in order to shed light on the complex behavior of these temperatures with varying orbital and thermal parameters and in order to consider possible cometary nucleus thermal evolution targets for comet rendezvous and/or sample-return missions. The concept of 'new' and 'old' comets, classified in terms of how many passages around the sun have been made, may take on new meaning in view of the present demonstration that differences in eccentricity between two comets having the same period and number of apparitions may yield significantly different internal temperature profiles.

Herman, Gary↗
DOE OSTI2022

The role of geographical spreaders in infectious pattern formation and front propagation speeds

The pattern formation and spatial spread of infectious populations are investigated using a kernel-based Susceptible–Infectious–Recovered (SIR) model applicable across a wide range of basic reproduction numbers R O . The goal is to examine the role of geographical spreaders on transient spatial pattern formation of infectious populations and the associated maximum invasive front speeds c max . In the simulations conducted here, geographical spreaders are defined as a portion of the infected population Φ experiencing high mobility between identical communities. The spatial organization of the infected population and c max are determined when the infections are randomly initiated in space within multiple communities. For small but finite , scaling analysis and numerical simulations in 1-dimension suggest that when the spreading kernel is Gaussian-shaped, where is the inverse of the infectious duration. This finding for agrees with a diffusion-based representation of mobility in 1-D. Numerical simulations in 2-D across wide-ranging suggest that , the variance of the spatial kernel describing mobility of long-distance geographical spreaders across communities, determines the spatial organization of infections across communities. When (long-distance mobility, where is the minimum spatial extent defining adjacent communities), the infectious population will experience a transient but spatially coherent pattern with a wavelength that can be derived from the spreading kernel properties. Moreover, the 2-D simulations for the bounded kernel suggest that attainment of is also dictated by but the magnitude is not sensitive to unlike diffusion-based models.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Leveraging inter-firm influence in the diffusion of energy efficiency technologies: An agent-based model

Energy efficiency technologies (EETs) are crucial for saving energy and reducing carbon dioxide emissions. However, the diffusion of EETs in small and medium-sized enterprises is rather slow. Literature shows the interactions between innovation adopters and potential adopters have significant impacts on innovation diffusion. Enterprises lack the motivation to share information, and EETs usually lack observability, which suppress the inter-firm influence. Thus, an information platform, together with proper policies encouraging or forcing enterprises to disclose EET-related information, should help harness inter-firm influence to accelerate EETs' diffusion. To explore whether and how such an information platform affects EETs' diffusion in small and medium-sized enterprises, this report builds an agent-based model to mimic EET diffusion processes. Based on a series of controlled numerical experiments, some counter-intuitive phenomena are discovered and explained. The results show that the information platform is a double-edged sword that notably accelerates EETs' diffusion by approximately 47% but may also boost negative information to diffuse even faster and delay massive adoption of EETs. Increasing network density and the intensity of inter-firm influence are effective to speed EET diffusion, but their impacts diminish drastically after reaching some critical values (0.05 and 0.15 respectively) and eventually harm the stability of the system. Ultimately, the findings implicate that EET suppliers should carefully launch their promising but immature products; policies that can reduce the perceived risk by enterprises and the effort to maintain an informative rather than judgmental information platform can prominently mitigate the negative side effects brought by high fluidity of information.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1994

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.↗