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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary-/dual- pore network: Part II - numerical model

To increase the market share of electric vehicles, it is desirable to reduce the battery charge times, which are significantly limited by poor electrolyte transport. A high rate charging is achievable by using expensive and low energy density cells with thin electrodes. For higher energy density cells, new electrolytes with improved conductivity and diffusivity and/or electrodes with advanced architecture are required to boost the electrolyte transport, leading to a more uniform utilization of active materials. In our previous work, an analytical model was developed to investigate the effect of secondary pore network (SPN) on electrolyte transport and the configuration of SPN was optimized by enforcing equal characteristic diffusion times in through-plane and in-plane directions. Here, to evaluate the effect of SPN on the fast-charging capability of lithium-ion batteries, a 2D physics-based electrochemical model is developed with SPN in either one or both electrodes. Additionally, the effect of SPN on cell energy density and lithium plating is investigated for cells with different loadings and electrode porosities. Combining SPN with elevated charging temperatures, the model predicts that the volumetric discharge energy density of a 3 mA h/cm2 cell can reach 270 Wh/L after a 6C constant-current charging.

25 ENERGY STORAGE↗
DOE OSTI2017

The Effects of Magnetic Fields on the Dynamics of Radiation Pressure–dominated Massive Star Envelopes

We use three-dimensional radiation magnetohydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure–dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100 to 1000 G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low-density regions. Magnetic buoyancy brings small-scale magnetic fields to the photosphere and increases the vertical energy transport, with the energy advection velocity proportional to the Alfvén velocity, although in all cases we study, photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared with the purely hydrodynamic case, reaching $\approx 100\,{{\rm{km}}{\rm{s}}}^{-1}$ at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Overlapping Schwarz Methods Are Not Anisotropy‐Robust Multigrid Smoothers

We analyze overlapping multiplicative Schwarz methods as smoothers in the geometric multigrid solution of two-dimensional anisotropic diffusion problems. For diffusion equations, it is well known that the smoothing properties of point-wise smoothers, such as Gauss Seidel, rapidly deteriorate as the strength of anisotropy increases. On the other hand, global smoothers based on line smoothing are known to generally provide good smoothing for diffusion problems, independent of the anisotropy strength. Here, a natural question is whether global methods are really necessary to achieve good smoothing in such problems, or whether it can be obtained with locally overlapping block smoothers using sufficiently large blocks and overlap. Through local Fourier analysis and careful numerical experimentation, we show that global methods are indeed necessary to achieve anisotropy-robust smoothing. Specifically, for any fixed block size bounded sufficiently far away from the global domain size, we find that the smoothing properties of overlapping multiplicative Schwarz rapidly deteriorate with increasing anisotropy, irrespective of the amount of overlap between blocks. Moreover, our results indicate that anisotropy-robust smoothing requires blocks of diameter 𝒪⁡(𝜖 −1/2 ) for anisotropy ratio 𝜖 ∈(0,1] .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Convergence analysis of two-node CMFD method for two-group neutron diffusion eigenvalue problem

In this paper, the nonlinear coarse-mesh finite difference method with two-node local problem (CMFD2N) is proven to be unconditionally stable for neutron diffusion eigenvalue problems. The explicit current correction factor (CCF) is derived based on the two-node analytic nodal method (ANM2N), and a Fourier stability analysis is applied to the linearized algorithm. It is shown that the analytic convergence rate obtained by the Fourier analysis compares very well with the numerically measured convergence rate. It is also shown that the theoretical convergence rate is only governed by the converged second harmonic buckling and the mesh size. It is also noted that the convergence rate of the CCF of the CMFD2N algorithm is dependent on the mesh size, but not on the total problem size. This is contrary to expectation for eigenvalue problem. The novel points of this paper are the analytical derivation of the convergence rate of the CMFD2N algorithm for eigenvalue problem, and the convergence analysis based on the analytic derivations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Rayleigh-Darcy convection with hydrodynamic dispersion

We investigate the effect of hydrodynamic dispersion on convection in porous media by performing direct numerical simulations (DNS) in a two-dimensional Rayleigh-Darcy domain. Scaling analysis of the governing equations shows that the dynamics of this system are not only controlled by the classical Rayleigh-Darcy number based on molecular diffusion, Ra m , and the domain aspect ratio, but also controlled by two other dimensionless parameters: the dispersive Rayleigh number Ra d =H/α t and the dispersivity ratio r=α l /α t , where H is the domain height and αt and αl are the transverse and longitudinal dispersivities, respectively. For Δ=Ra d /Ra m >O(1), the influence from the mechanical dispersion is minor; for Δ≲0.02, however, the flow pattern is determined by Rad while the convective flux is F~c(Ra d )Ra m for large Ram. Our DNS results also show that the increase of mechanical dispersion, i.e., decreasing Rad, will coarsen the convective pattern by increasing the plume spacing. Moreover, the inherent anisotropy of mechanical dispersion breaks the columnar structure of the megaplumes at large Ra m , if Ra d <5000. This results in a fan-flow geometry that reduces the convective flux.

42 ENGINEERING↗
DOE OSTI2019

The Bloch equation for spin dynamics in electron storage rings: Computational and theoretical aspects

In this paper, we describe our work on spin polarization in high-energy electron storage rings which we base on the Full Bloch equation (FBE) for the polarization density and which aims towards the [Formula: see text] option of the proposed Future Circular Collider (FCC-ee) and the proposed Circular Electron Positron Collider (CEPC). The FBE takes into account non spin-flip and spin-flip effects due to synchrotron radiation including the spin-diffusion effects and the Sokolov–Ternov effect with its Baier–Katkov generalization as well as the kinetic-polarization effect. This mathematical model is an alternative to the standard mathematical model based on the Derbenev–Kondratenko formulas. For our numerical and analytical studies of the FBE, we develop an approximation to the latter to obtain an effective FBE. This is accomplished by finding a third mathematical model based on a system of stochastic differential equations (SDEs) underlying the FBE and by approximating that system via the method of averaging from perturbative ODE theory. We also give an overview of our algorithm for numerically integrating the effective FBE. This discretizes the phase space using spectral methods and discretizes time via the additive Runge–Kutta (ARK) method which is a high-order semi-implicit method. We also discuss the relevance of the third mathematical model for spin tracking.

Physics↗
DOE OSTI2016

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). III. THERMAL AND COMPOSITIONAL TRANSPORT THROUGH NON-LAYERED ODDC

Oscillatory double-diffusive convection (ODDC; also known as semi-convection) refers to a type of double-diffusive instability that occurs in regions of planetary and stellar interiors that have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. first showed that ODDC can either spontaneously form layers that significantly enhance the transport of heat and chemical species compared to microscopic transport or remain in a state dominated by large-scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series focused on identifying under what conditions layers form and quantifying transport through layered systems. Here we proceed to characterize transport through systems that are unstable to ODDC, but do not undergo spontaneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both two-dimensional (2D) and 3D numerical simulations, and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first-level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport systems unstable to ODDC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

COMBINING DIFFUSIVE SHOCK ACCELERATION WITH ACCELERATION BY CONTRACTING AND RECONNECTING SMALL-SCALE FLUX ROPES AT HELIOSPHERIC SHOCKS

Computational and observational evidence is accruing that heliospheric shocks, as emitters of vorticity, can produce downstream magnetic flux ropes and filaments. This led Zank et al. to investigate a new paradigm whereby energetic particle acceleration near shocks is a combination of diffusive shock acceleration (DSA) with downstream acceleration by many small-scale contracting and reconnecting (merging) flux ropes. Using a model where flux-rope acceleration involves a first-order Fermi mechanism due to the mean compression of numerous contracting flux ropes, Zank et al. provide theoretical support for observations that power-law spectra of energetic particles downstream of heliospheric shocks can be harder than predicted by DSA theory and that energetic particle intensities should peak behind shocks instead of at shocks as predicted by DSA theory. In this paper, a more extended formalism of kinetic transport theory developed by le Roux et al. is used to further explore this paradigm. We describe how second-order Fermi acceleration, related to the variance in the electromagnetic fields produced by downstream small-scale flux-rope dynamics, modifies the standard DSA model. The results show that (i) this approach can qualitatively reproduce observations of particle intensities peaking behind the shock, thus providing further support for the new paradigm, and (ii) stochastic acceleration by compressible flux ropes tends to be more efficient than incompressible flux ropes behind shocks in modifying the DSA spectrum of energetic particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Preliminary Evidence for a Virial Shock around the Coma Galaxy Cluster

Galaxy clusters, the largest gravitationally bound objects in the universe, are thought to grow by accreting mass from their surroundings through large-scale virial shocks. Due to electron acceleration in such a shock, it should appear as a γ-ray, hard X-ray, and radio ring, elongated toward the large-scale filaments feeding the cluster, coincident with a cutoff in the thermal Sunyaev–Zel’dovich (SZ) signal. However, no such signature was found until now, and the very existence of cluster virial shocks has remained a theory. We find preliminary evidence for a large γ-ray ring of ∼5 Mpc minor axis around the Coma cluster, elongated toward the large-scale filament connecting Coma and Abell 1367, detected at the nominal 2.7σ confidence level (5.1σ using control signal simulations). The γ-ray ring correlates both with a synchrotron signal and with the SZ cutoff, but not with Galactic tracers. The γ-ray and radio signatures agree with analytic and numerical predictions if the shock deposits ∼1% of the thermal energy in relativistic electrons over a Hubble time and ∼1% in magnetic fields. The implied inverse Compton and synchrotron cumulative emission from similar shocks can contribute significantly to the diffuse extragalactic γ-ray and low-frequency radio backgrounds. Our results, if confirmed, reveal the prolate structure of the hot gas in Coma, the feeding pattern of the cluster, and properties of the surrounding large-scale voids and filaments. The anticipated detection of such shocks around other clusters would provide a powerful new cosmological probe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

THE EFFECT OF ROTATION ON OSCILLATORY DOUBLE-DIFFUSIVE CONVECTION (SEMICONVECTION)

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermocompositional layers that significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and nonlayered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work, we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations, we then analyze the effect of rotation on properties of layered and nonlayered ODDC, and we study how the angle of the rotation axis with respect to the direction of gravity affects layering. We find that rotating systems can be broadly grouped into two categories based on the strength of rotation. The qualitative behavior in the more weakly rotating group is similar to nonrotating ODDC, but strongly rotating systems become dominated by vortices that are invariant in the direction of the rotation vector and strongly influence transport. We find that whenever layers form, rotation always acts to reduce thermal and compositional transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Prediction of Cooling Curves for Controlled Unidirectional Solidification Under the Influence of Shrinkage: A Semi-analytical Approach

Quality of casting products can be significantly improved by controlling the crystal growth rate. Both dendritic and equiaxed grain structures depend strongly on crystal growth rate along the crystallographic direction. Uncontrolled solidification processes may lead to the formation of porous and columnar material with an extremely non-homogeneous composition distribution. These undesirable micro-structures can be avoided by adopting an optimal cooling rate (K/s) during the crystal growth. The present work focuses on determining suitable cooling curves to obtain desired unidirectional crystal growth rates, by implementing a semi-analytical model. Here, we defined cooling curve as the time history of temperature evolution at the cold boundary. The proposed semi-analytical heat transfer model is diffusion driven and accounts for shrinkage effect during solidification. The proposed model is found capable of predicting cooling curves for prescribed interface velocities. The model is self-sufficient in validating itself by means of solving inverse and forward problems. The cooling curve prediction by the model is validated by using it as the transient boundary condition for an existing enthalpy updating scheme-based numerical model. The proposed model is also applied to empirical data reported in literature, to validate the model against experimental results. All the validations showcased excellent fit. The proposed model will facilitate controlled unidirectional crystal growth rate by providing cooling rate (K/s) as an input for the experiments involving directional solidification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatio-temporal control of phenylpropanoid biosynthesis by inducible complementation of a cinnamate 4-hydroxylase mutant

Abstract Cinnamate 4-hydroxylase (C4H) is a cytochrome P450-dependent monooxygenase that catalyzes the second step of the general phenylpropanoid pathway. Arabidopsis reduced epidermal fluorescence 3 (ref3) mutants, which carry hypomorphic mutations in C4H, exhibit global alterations in phenylpropanoid biosynthesis and have developmental abnormalities including dwarfing. Here we report the characterization of a conditional Arabidopsis C4H line (ref3-2pOpC4H), in which wild-type C4H is expressed in the ref3-2 background. Expression of C4H in plants with well-developed primary inflorescence stems resulted in restoration of fertility and the production of substantial amounts of lignin, revealing that the developmental window for lignification is remarkably plastic. Following induction of C4H expression in ref3-2pOpC4H, we observed rapid and significant reductions in the levels of numerous metabolites, including several benzoyl and cinnamoyl esters and amino acid conjugates. These atypical conjugates were quickly replaced with their sinapoylated equivalents, suggesting that phenolic esters are subjected to substantial amounts of turnover in wild-type plants. Furthermore, using localized application of dexamethasone to ref3-2pOpC4H, we show that phenylpropanoids are not transported appreciably from their site of synthesis. Finally, we identified a defective Casparian strip diffusion barrier in the ref3-2 mutant root endodermis, which is restored by induction of C4H expression.

Plant Sciences↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2015

FCRD Transmutation Fuels Handbook 2015

Transmutation of minor actinides such as Np, Am, and Cm in spent nuclear fuel is of international interest because of its potential for reducing the long-term health and safety hazards caused by the radioactivity of the spent fuel. One important approach to transmutation (currently being pursued by the DOE Fuel Cycle Research & Development Advanced Fuels Campaign) involves incorporating the minor actinides into U-Pu-Zr alloys, which can be used as fuel in fast reactors. It is, therefore, important to understand the properties of U-Pu-Zr alloys, both with and without minor actinide additions. In addition to requiring extensive safety precautions, alloys containing U and Pu are difficult to study for numerous reasons, including their complex phase transformations, characteristically sluggish phase-transformation kinetics, tendency to produce experimental results that vary depending on the histories of individual samples, and sensitivity to contaminants such as oxygen in concentrations below a hundred parts per million. Many of the experimental measurements were made before 1980, and the level of documentation for experimental methods and results varies widely. It is, therefore, not surprising that little is known with certainty about U-Pu-Zr alloys, and that general acceptance of results sometimes indicates that there is only a single measurement for a particular property. This handbook summarizes currently available information about U, Pu, Zr, and alloys of two or three of these elements. It contains information about phase diagrams and related information (including phases and phase transformations); heat capacity, entropy, and enthalpy; thermal expansion; and thermal conductivity and diffusivity. In addition to presenting information about materials properties, it attempts to provide information about how well the property is known and how much variation exists between measurements. Although the handbook includes some references to publications about modeling, its primary focus is experimental data. Most of the data has been published elsewhere (although scattered throughout numerous references, some quite obscure); however, some data is presented here for the first time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Numerical evaluation of AGR-2 fission product release

The AGR-2 experiment produced normal operation and accident condition fuel performance data for tri-structural isotropic (TRISO) particles with UCO and UO 2 kernels. Data for compacts with no failed particles and for compacts with one or more failed particles are available. To model the fission product diffusion and release from these compacts, it is important to consider the computational mesh as well as the heat conduction and mass diffusion properties of the materials in a TRISO particle. Code and solution verification studies with the Bison fuel performance code were performed to show that the code is computing correct solutions. Bison comparisons to post-irradiation examination data and PARFUME predictions were made for several sets of compacts. Furthermore, these comparisons considered silver, cesium, strontium, and krypton with and without failed particles. Bison predictions closely match those of PARFUME, while both codes generally overpredict post-irradiation examination data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Frequency spirals

We study the dynamics of coupled phase oscillators on a two-dimensional Kuramoto lattice with periodic boundary conditions. For coupling strengths just below the transition to global phase-locking, we find localized spatiotemporal patterns that we call “frequency spirals.” These patterns cannot be seen under time averaging; they become visible only when we examine the spatial variation of the oscillators' instantaneous frequencies, where they manifest themselves as two-armed rotating spirals. In the more familiar phase representation, they appear as wobbly periodic patterns surrounding a phase vortex. Unlike the stationary phase vortices seen in magnetic spin systems, or the rotating spiral waves seen in reaction-diffusion systems, frequency spirals librate: the phases of the oscillators surrounding the central vortex move forward and then backward, executing a periodic motion with zero winding number. We construct the simplest frequency spiral and characterize its properties using analytical and numerical methods. Simulations show that frequency spirals in large lattices behave much like this simple prototype.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Role of the edge stochastic layer in density pump-out by resonant magnetic perturbations

The role of the edge stochastic layer on particle transport is addressed in DIII-D plasmas with applied resonant magnetic perturbations (RMPs) causing density pump-out (a term used to denote the density reduction at the top of the pedestal due to the applied RMPs). Using an analytical model that adds the stochastic parallel transport of electrons (Harvey et al 1981 Phys. Rev. Lett. 47 102) to the fluid equations, the ambipolar radial electric field and particle flux are calculated simultaneously. In this model the nonambipolar electron particle flux, driven by the stochastic magnetic field, is predominantly balanced by the nonambipolar perpendicular ion flux, driven by toroidal viscosity, across a narrow stochastic layer of order two percent of the plasma radius (0.98 < ψn < 1). The model reproduces the level of density pump-out and its dependence on RMP amplitude near the separatrix for three plasma discharges for which the pedestal foot is at medium to high collisionality $(2\leqslant {\nu }_{\text{sep}}^{\ast }\leqslant 30)$. In this work, the experimentally observed and numerically tested inverse density dependence of density pump-out (Hu et al 2020 Nucl. Fusion 60 076001) is accurately captured by the model: an increase in collisionality with density in the pedestal foot results in a decrease in stochastic diffusivity, and hence a decrease in the level of pump-out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗