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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Axisymmetric global gravitational equilibrium for magnetized, rotating hot plasma

In this work, we present analytic solutions for three-dimensional magnetized axisymmetric equilibria confining rotating hot plasma in a gravitational field. Our up–down symmetric solution to the full Grad–Shafranov equation can exhibit equatorial plane localization of the plasma density and current, resulting in disk equilibria for the plasma density. For very weak magnetic fields and high plasma pressure, we find strongly rotating thin plasma disk gravitational equilibria that satisfy strict Keplerian motion provided the gravitational energy is much larger than the plasma pressure, which must be large compared to the magnetic energy of the poloidal magnetic field. When the rotational energy exceeds the gravitational energy and it is larger than the plasma pressure, diffuse disk equilibrium solutions continue to exist provided the poloidal magnetic energy remains small. For stronger magnetic fields and lower plasma pressure and rotation, we can also find gravitational equilibria with strong localization to the equatorial plane. However, a toroidal magnetic field is almost always necessary to numerically verify these equilibria are valid solutions in the presence of gravity for the cases considered in Catto & Krasheninnikov (J. Plasma Phys., vol. 81, 2015, 105810301). In all cases both analytic and numerical results are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗
DOE OSTI2023

Double Diffusion Maps and their Latent Harmonics for scientific computations in latent space

In this work, we introduce a data-driven approach to building reduced dynamical models through manifold learning; the reduced latent space is discovered using Diffusion Maps (a manifold learning technique) on time series data. A second round of Diffusion Maps on those latent coordinates allows the approximation of the reduced dynamical models. This second round enables mapping the latent space coordinates back to the full ambient space (what is called lifting); it also enables the approximation of full state functions of interest in terms of the reduced coordinates. In our work, we develop and test three different reduced numerical simulation methodologies, either through pre-tabulation in the latent space and integration on the fly or by going back and forth between the ambient space and the latent space. The data-driven latent space simulation results, based on the three different approaches, are validated through (a) the latent space observation of the full simulation through the Nyström Extension formula, or through (b) lifting the reduced trajectory back to the full ambient space, via Latent Harmonics. Latent space modeling often involves additional regularization to favor certain properties of the space over others, and the mapping back to the ambient space is then constructed mostly independently from these properties; here, we use the same data-driven approach to construct the latent space and then map back to the ambient space.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Faster ablative Kelvin–Helmholtz instability growth in a magnetic field

Shear flows along a plasma interface will quickly grow unstable due to the Kelvin–Helmholtz instability. If there is a concurrent temperature gradient across the interface, higher modes are stabilized by the thermal diffusion. These ablative effects must be considered in, for example, jet features in inertial confinement fusion hot-spots, or plasma plumes in young supernovae. We show that magnetization of the plasma can greatly affect the instability, even if magnetic pressure is small. This is because electrons are localized by their gyromotion, reducing the heat flux and material ablation. We use a two-dimensional numerical extended-magnetohydrodynamics approach to assess this effect for dense fusion conditions. In comparison with the unmagnetized case, self-generated Biermann fields make only a minor difference to growth rates. However, simulations with a large 50 kT external field found that the growth rate of the least stable mode increased by 40%. In conclusion, this has implications for mix processes in Z-pinches and magnetized inertial confinement fusion concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Magnetically self-regulated formation of early protoplanetary disks

The formation of protoplanetary disks during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disk radius, r≃18 au (η{sub AD}/0.1s){sup 2/9}(B{sub z}/0.1G){sup −4/9}(M/0.1M{sub ⊙}){sup 1/3}, where M is the total disk plus protostar mass, η{sub AD} is the ambipolar diffusion coefficient, and B {sub z} is the magnetic field in the inner part of the core. This is significantly smaller than the disks that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disk radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disk radii upon the various relevant quantities, suggesting weak variations among class-0 disk sizes. In some cases, we note the onset of large spiral arms beyond this radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Novel Multi-scale Simulation Strategy for Turbulent Reacting Flows

In this project, a new methodology was proposed to bridge the gap between Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES). This novel methodology, titled Lattice-Based Multiscale Simulation (LBMS), creates a lattice structure of One-Dimensional Turbulence (ODT) models. This model has been shown to capture turbulent combustion with high fidelity by fully resolving interactions between turbulence and diffusion. By creating a lattice of ODT models, which are then coupled, LBMS overcomes the shortcomings of ODT, which are its inability to capture large scale three dimensional flow structures. However, by spacing these lattices significantly apart, LBMS can avoid the curse of dimensionality that creates untenable computational costs associated with DNS. This project has shown that LBMS is capable of reproducing statistics of isotropic turbulent flows while coarsening the spacing between lines significantly. It also investigates and resolves issues that arise when coupling ODT lines, such as flux reconstruction perpendicular to a given ODT line, preservation of conserved quantities when eddies cross a course cell volume and boundary condition application. Robust parallelization is also investigated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Performance Analysis of Transposition Models Simulating Solar Radiation on Inclined Surfaces

Transposition models have been widely used in the solar energy industry to simulate solar radiation on inclined photovoltaic panels. Following numerous studies comparing the performance of transposition models, this work aims to understand the quantitative uncertainty in state-of-the-art transposition models and the sources leading to the uncertainty. Our results show significant differences between two highly used isotropic transposition models, with one substantially underestimating the diffuse plane-of-array irradiances when diffuse radiation is perfectly isotropic. In the empirical transposition models, the selection of the empirical coefficients and land surface albedo can both result in uncertainty in the output. This study can be used as a guide for the future development of physics-based transposition models and evaluations of system performance.

14 SOLAR ENERGY↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2025

Enabling fast-charging of lithium-ion batteries through printed electrodes

It has been well recognized that introducing secondary porous networks (SPNs) into the electrodes can effectively improve the electrochemical performance of lithium-ion batteries (LIBs), especially under fast-charging operations. However, the process complexity and high cost limit the commercial success of advanced electrodes with SPNs. To address this issue, we developed a facile screen-printing process to produce structured graphite electrodes with SPNs. The experimental results demonstrated that, by tuning the diameter and center-to-center (C2C) distance of emulsion dots on the stencil screen, the pore diameters and C2C pore distances of SPNs in screenprinted electrodes can be precisely controlled in the range of 100 mu m to 1 mm and 100 mu m to 3 mm respectively. In addition, the SPNs with hexagonal and square-shape pore alignments have also been imprinted onto the electrode coatings through adjusting the patterns of screen stencils. Used as anodes, the printed graphite electrodes demonstrated significantly reduced overpotential and voltage fluctuation under fast-charging operations from 2C to 6C. Coupled with LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) cathodes, the full cells with printed graphite anodes exhibited an unprecedently stable performance with almost no capacity decay up to 170 cycles when charged to 80 % SOC at 2C. Observations from electron microscopy showed plated lithium undetectable at the surface of printed graphite electrodes after numerous cycles. The electrochemical analysis on the voltage evolution during the cell rest period indicated the significantly delayed onset of lithium plating in the presence of printed graphite electrodes. In conclusion, all these results suggest that the significantly improved cell performance is associated with the shortened Li-ion diffusion distance, reduced polarization and suppressed Li plating in the printed electrodes with patterned SPNs.

25 ENERGY STORAGE↗
DOE OSTI2018

Systematic search for tidal features around nearby galaxies: I. Enhanced SDSS imaging of the Local Volume

In hierarchical models of galaxy formation, stellar tidal streams are expected around most, if not all, galaxies. Although these features may provide useful diagnostics of the Λ CDM model, their observational properties remain poorly constrained because they are challenging to detect and interpret and have been studied in detail for only a sparse sampling of galaxy population. More quantitative, systematic approaches are required. We advocate statistical analysis of the counts and properties of such features in archival wide-field imaging surveys for a direct comparison against results from numerical simulations. We aim to study systematically the frequency of occurrence and other observational properties of tidal features around nearby galaxies. The sample we construct will act as a foundational dataset for statistical comparison with cosmological models of galaxy formation. Our approach is based on a visual classification of diffuse features around a volume-limited sample of nearby galaxies, using a post-processing of Sloan Digital Syk Survey (SDSS) imaging optimized for the detection of stellar structure with low surface brightness. At a limiting surface brightness of 28 mag arcsec –2 , 14% of the galaxies in our sample exhibit evidence of diffuse features likely to have arisen from minor merging events. Our technique recovers all previously known streams in our sample and yields a number of new candidates. Consistent with previous studies, coherent arc-like features and shells are the most common type of tidal structures found in this study. Here, we conclude that although some detections are ambiguous and could be corroborated or refuted with deeper imaging, our technique provides a reliable foundation for the statistical analysis of diffuse circumgalactic features in wide-area imaging surveys, and for the identification of targets for follow-up studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Collisional dependence of Alfvén mode saturation in tokamaks

Saturation of Alfvén modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Finally, saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Collisional dependence of Alfven mode saturation in tokamaks

Saturation of \alfven modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A numerical study of ignition and flame development characteristics in GCI combustion using large eddy simulations and chemical explosive mode analysis

Here, this work investigates the ignition and flame development processes of low reactivity fuel combustion under compression ignition conditions based on the large eddy simulation approach. The chemical explosive mode analysis (CEMA) is employed to characterize the local combustion features, including gas-liquid fuel zone, auto-ignition, diffusion-assisted, extinction, cool flame and post-ignition zone, among which auto-ignition and post-ignition are found to play a key role in the overall heat release process. The local flame propagation modes in gasoline compression ignition (GCI) are determined by quantifying the relative magnitude of diffusion/chemistry at a representative progress variable in the pre-ignition zone. The results show that autoignition fronts and deflagration waves exist simultaneously in the ignition and intense high temperature heat release (HTHR) stages, but autoignition fronts dominate. In addition, the chemical kinetic processes of four heat release periods are analyzed. The heat release during the ignition period is found to be dominated by the reactions CH 3 + H (+M) <=> CH 4 (+M) and CH 3 CHO + H <=> CH 2 CHO + H 2 . The reaction CH 2 OH + OH <=> CH 2 O + H 2 O always plays an important role in the heat releases during the other three combustion stages including intense HTHR, moderate HTHR and post-combustion.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

MODELLING FIBER ORIENTATION DURING ADDITIVE MANUFACTURING COMPRESSION MOLDING PROCESSES

The production of high-performance thermoplastic composites reinforced with short carbon fibers can be achieved by a novel “additive manufacturing-compression molding” technique. An advantage of such combination is two-fold: controlled fiber orientation in additive manufacturing and less void content by compression molding. In this study, a computational fluid dynamics model has been developed to predict the behavior of printed layers during fiber-reinforced thermoplastic extrusion and subsequent compression molding. The fiber orientation was modelled with a simple quadratic closure model. The interaction between the fibers is included using a rotary diffusion coefficient which becomes significant in concentrated regimes. Finally, the second order orientation tensor is coupled with the momentum equation as an anisotropic part of the stress term. The effect of processing parameters on the behavior of printed layers was investigated to determine the favorable printing scenarios. The developed numerical model enables design of high-performance composites with tunable mechanical properties.

Seta, Berin↗
DOE OSTI2016

Measurement of thermal gap conductance using the laser flash method

Thermal gap conductance refers to overall heat conductance through a gas-filled gap as interface thermal resistance between the gas and surface has to be considered. The analysis of new ATR (Advanced Test Reactor) fuel experiments requires experimental data and models for the gap conductance of the cladding/fuel capsule gap, which typically varies from 5 μm to 50 μm and is filled with He at 1-3 atm. Although past research has studied the gap conductance between UO{sub 2} fuel and Zircaloy cladding, new experiments are needed to provide the gap conductance data between HT9 and stainless steel 316, and between advanced fuels and advanced cladding in the Accident Tolerant Fuels program. In this paper, the laser flash method is modified to measure the gap conductance. The gap is formed by welding two stainless steel 316 disks with three spacers. The temperature at the rear surface, after exposure to the laser pulse at the front surface, is modeled with 1D heat transfer through three layers. The analytic solution is then fitted with the measured temperature to determine gap conductance. The sensitivity analysis based on the 1D assumption and 3D numerical simulations were performed to assess the sensitivity and uncertainty of the measurement. It is concluded the heat capacity of the gas in the gap has minimal impact on the measurement, which indicates that the thermal property determined by this method is the thermal conductance, instead of the thermal diffusivity. Convective heat loss through sidewalls under the current experimental design has insignificant influences. Heat transfer through spacers may cause up to 10% overestimation of thermal gap conductance, but the effects can be mostly eliminated by the design of grooves on samples. The result established a pathway for gap conductance determination. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Fast and accurate implementation of Fourier spectral approximations of nonlocal diffusion operators and its applications

This work is concerned with the Fourier spectral approximation of various integral differential equations associated with some linear nonlocal diffusion and peridynamic operators under periodic boundary conditions. For radially symmetric kernels, the nonlocal operators under consideration are diagonalizable in the Fourier space so that the main computational challenge is on the accurate and fast evaluation of their eigenvalues or Fourier symbols consisting of possibly singular and highly oscillatory integrals. For a large class of fractional power-like kernels, we propose a new approach based on reformulating the Fourier symbols both as coefficients of a series expansion and solutions of some simple ODE models. We then propose a hybrid algorithm that utilizes both truncated series expansions and high order Runge–Kutta ODE solvers to provide fast evaluation of Fourier symbols in both one and higher dimensional spaces. It is shown that this hybrid algorithm is robust, efficient and accurate. As applications, we combine this hybrid spectral discretization in the spatial variables and the fourth-order exponential time differencing Runge–Kutta for temporal discretization to offer high order approximations of some nonlocal gradient dynamics including nonlocal Allen–Cahn equations, nonlocal Cahn–Hilliard equations, and nonlocal phase-field crystal models. Numerical results show the accuracy and effectiveness of the fully discrete scheme and illustrate some interesting phenomena associated with the nonlocal models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation

In this report we propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite-volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗