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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Convergence Study of LR-NDA Using Fourier Analysis

Numerical solution of neutron transport equations is challenging. It often requires acceleration techniques to improve computational efficiency. In the last two decades, various acceleration methods have been developed. Among those methods, the coarse-mesh nonlinear diffusion acceleration (CM-NDA) technique, or commonly known as coarse-mesh finite difference (CMFD), becomes one of the most popular methods due to its effectiveness and simplicity in recent years. However, CM-NDA will lose its effectiveness and could become unstable when the mesh optical thickness becomes large. There exists a number of stabilization techniques to improve the stability of CM-NDA, such as pCMFD and most recent odCMFD. Wang et al., proposed a new CM-NDA scheme, LR-NDA, which employs local refinement on the coarse mesh. The new scheme can significantly improve the effectiveness and stability of CM-NDA. In this paper, we perform a convergence study of LR-NDA using Fourier analysis. The spectral radius calculated from Fourier analysis is consistent with the numerical results for a 1-D fixed source model problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

CMB anisotropies from patchy reionisation and diffuse Sunyaev-Zel'dovich effects

Anisotropies in the Cosmic Microwave Background (CMB) can be induced during the later stages of cosmic evolution, and in particular during and after the Epoch of Reionisation. Inhomogeneities in the ionised fraction, but also in the baryon density, in the velocity fields and in the gravitational potentials are expected to generate correlated CMB perturbations. We present a complete relativistic treatment of all these effects, up to second order in perturbation theory, that we solve using the numerical Boltzmann code (\SONG). The physical origin and relevance of all second order terms are carefully discussed. In addition to collisional and gravitational contributions, we identify the diffuse analogue of the blurring and kinetic Sunyaev-Zel'dovich (SZ) effects. Our approach naturally includes the correlations between the imprint from patchy reionisation and the diffuse SZ effects thereby allowing us to derive reliable estimates of the induced temperature and polarisation CMB angular power spectra. In particular, we show that the B -modes generated at intermediate length-scales (ℓ ≅ 100) have the same amplitude as the B -modes coming from primordial gravitational waves with a tensor-to-scalar ratio r =10{sup −4}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A mixing timescale model for TPDF simulations of turbulent premixed flames

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priori assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive scalar mixing across different flame regimes are appropriately accounted for.

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

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Design, analysis and preliminary numerical results for the nonconforming VEM for parabolic problems

The virtual element method (VEM) is a numerical methodology for solving partial differential equations on unstructured polytopal meshes. In this work, we extend the nonconforming high-order formulation suitable to time-dependent diffusion problems with possibly nonlinear forcing terms. Coupling high-order virtual element formulations in space with time-marching schemes based on Backward Differention Formulas (BDFs) makes a fully discrete method with optimal order of convergence in the Method of Lines framework. We investigate and assess the performance of the method through theoretical analysis and a set of numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adapting confocal Raman microscopy for in situ studies of redox transformations at electrode-electrolyte interfaces

Confocal Raman microscopy was applied to quantify redox species present within the diffusion layer adjacent to an electrode surface under potentiostatic control. A glass microscope coverslip with a thin indium tin oxide (ITO) coating served as both the working electrode and optical window for a microscope-stage mountable spectroelectrochemical cell. A high numerical aperture objective mounted in an inverted microscope frame just below the stage brought excitation radiation through the coverslip window and to a tight focus a few micrometers above the ITO film surface. Species diffusing into the confocal probe volume defined by the excitation beam focus and the collected light region were detected, identified and quantified based on their Raman scattering frequencies and intensities. Here, in measurements that interrogated the interconversion of ferrocyanide and ferricyanide ions as a function of applied voltage, least-squares regression analysis of spectral datasets predicted the formal potential and relative surface concentrations of the ions in good agreement with the expected Nernstian response. Preliminary studies of methyl viologen reduction at an ITO film/Nafion membrane interface were conducted and show the possibility for estimation of mass transport coefficients of redox species within ionic polymer materials.

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

Reaction rates for a generalized reaction-diffusion master equation

It has been established that there is an inherent limit to the accuracy of the reaction-diffusion master equation. Specifically, there exists a fundamental lower bound on the mesh size, below which the accuracy deteriorates as the mesh is refined further. In this paper we extend the standard reaction-diffusion master equation to allow molecules occupying neighboring voxels to react, in contrast to the traditional approach in which molecules react only when occupying the same voxel. We derive reaction rates, in two dimensions as well as three dimensions, to obtain an optimal match to the more fine-grained Smoluchowski model, and show in two numerical examples that the extended algorithm is accurate for a wide range of mesh sizes, allowing us to simulate systems that are intractable with the standard reaction-diffusion master equation. In addition, we show that for mesh sizes above the fundamental lower limit of the standard algorithm, the generalized algorithm reduces to the standard algorithm. We derive a lower limit for the generalized algorithm which, in both two dimensions and three dimensions, is on the order of the reaction radius of a reacting pair of molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Shape effect on drying behavior of cement-based materials: Mechanisms and numerical analysis

The shape of the specimen or structure has a direct impact on the behavior of drying concrete. This effect is not taken into account in most of shrinkage models and codes. Only the current B4 model and its previous version B3 model incorporate a shape factor. The determination of this coefficient is based on the nonlinear moisture diffusion analysis, which can explain the shape effect on the shrinkage rate. However, the model integrates this coefficient in the ultimate shrinkage values as well. An experimental study has been performed to confirm and understand this behavior. The mass loss and drying shrinkage of cylinders showed higher rate and ultimate values than the prisms, for the same effective cross section thickness. To explain this result, the drying, the differential shrinkage, the internal stress and the damaged depth were determined for both shapes through numerical modeling. The main input data of the model were the measured mass loss and shrinkage. The calculated desorption isotherm curve of the cylinder was below that of the prism, which is consistent with its higher ultimate mass-loss. The determination of stress profiles indicates slower stress mitigation for the cylinder. The cylinder develops a larger relative damaged area, defined as the part of the specimen where cracking is likely to occur due to tensile stresses. The results actually show that differential shrinkage between the edge and the core is more pronounced for the cylinder and decreases with time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical Modeling of Laser Heating and Evaporation of a Single Droplet

Laser technology is being widely studied for controlled energy deposition for a range of applications, including flow control, ignition, combustion, and diagnostics. The absorption and scattering of laser radiation by liquid droplets in aerosols affects propagation of the laser beam in the atmosphere, while the ignition and combustion characteristics in combustion chambers are influenced by the evaporation rate of the sprayed fuel. In this work, we present a mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion-dominated regime taking into account absorption of the laser radiation, evaporation process, and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for heating and evaporation of ethanol and water droplets. The two main regimes—continuous and pulsed laser heating—are explored. For continuous laser heating, the peak temperature is higher for larger droplets. For pulsed laser heating, when the peak irradiance is close to transition to the boiling regime, the temporal dynamics of the droplet temperature does not depend on the droplet size. With the empirical normalization of time, the dynamics of the droplet shrinkage and cooling are found to be independent of droplet sizes and peak laser intensities.

42 ENGINEERING↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Two-Relaxation-Time Lattice Boltzmann Method and its Application to Advective-Diffusive-Reactive Transport

The lattice Boltzmann method (LBM) based on single-relaxation-time (SRT) or multiple-relaxation-time (MRT) collision operators is widely used in simulating flow and transport phenomena. The LBM based on two-relaxation-time (TRT) collision operators possesses strengths from the SRT and MRT LBMs, such as its simple implementation and good numerical stability, although tedious mathematical derivations and presentations of the TRT LBM hinder its application to a broad range of flow and transport phenomena. This paper describes the TRT LBM clearly and provides a pseudocode for easy implementation. Various transport phenomena were simulated using the TRT LBM to illustrate its applications in subsurface environments. These phenomena include advection-diffusion in uniform flow, Taylor dispersion in a pipe, solute transport in a packed column, reactive transport in uniform flow, and bacterial chemotaxis in porous media. Finally, the TRT LBM demonstrated good numerical performance in terms of accuracy and stability in predicting these transport phenomena. Therefore, the TRT LBM is a powerful tool to simulate various geophysical and biogeochemical processes in subsurface environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

MODELING AND ANALYSIS OF FISSION PRODUCT TRANSPORT IN THE AGR-3/4 EXPERIMENT

In this work we describe the ongoing modeling and analysis efforts in support of the AGR-3/4 experiment. AGR-3/4 is intended to provide data to assess fission product retention and transport (e.g., diffusion coefficients) in fuel matrix and graphite materials. We describe a set of pre-test predictions that incorporate the results of detailed thermal and fission product release models into a coupled 1D radial diffusion model of the experiment, using diffusion coefficients reported in the literature for Ag, Cs, and Sr. We make some comparisons of the predicted Cs profiles to preliminary measured data for Cs and find these to be reasonable, in most cases within an order of magnitude. Our ultimate objective is to refine the diffusion coefficients using AGR-3/4 data, so we identify an analytical method for doing so and demonstrate its efficacy via a series of numerical experiments using the model predictions. Finally, we discuss development of a post-irradiation examination plan informed by the modeling effort and simulate some of the heating tests that are tentatively planned.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Theoretical Analysis of the EPR Spectrum Saturation Effect with Account for Spectral Diffusion in a System with a Gaussian Distribution of Spin Resonance Frequencies

We have analyzed the effect of saturation of the EPR spectrum for the model situation when a set of resonance frequencies of spin packets is described by a Gaussian distribution and spectral diffusion, which is a random process without correlation, exists. We have obtained the form of the spectrum for arbitrary values of magnetic–resonant parameters: relaxation times for the longitudinal and transverse components of the magnetization vector of the system, the spectral diffusion rate, the dispersion of the Gaussian distribution of resonance frequency, and the magnetic induction of the microwave field. We have derived an analytic expression for the saturation curve of the EPR spectrum. Detailed analysis of the dependence of the spectrum saturation curve on the spin system parameters is performed for the optimal case when the microwave field frequency coincides with the average frequency of the Gaussian distribution for resonance frequencies. The value of the microwave field induction at which the peak of the saturation curve is attained has been determined for arbitrary values of the paramagnetic relaxation time, spectral diffusion rate, and dispersion of the Gaussian spin frequency distribution. We have formulated the algorithm for determining the spin–lattice relaxation time using the experimentally determined microwave field induction for which the saturation curve passes through the peak. It should be noted that numerical calculations have been performed for parameters typical of EPR spectroscopy. The general analytic expression derived for the shape of the spectrum in the saturation conditions and qualitative conclusions concerning the effect of spectral diffusion on the saturation of stationary spectra obtained in this study can also be used in other fields of spectroscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

IMEX Runge-Kutta Parareal for Non-diffusive Equations

Parareal is a widely studied parallel-in-time method that can achieve meaningful speedup on certain problems. However, it is well known that the method typically performs poorly on non-diffusive equations. This paper analyzes linear stability and convergence for IMEX Runge-Kutta Parareal methods on non-diffusive equations. By combining standard linear stability analysis with a simple convergence analysis, we find that certain Parareal configurations can achieve parallel speedup on non-diffusive equations. These stable configurations all posses low iteration counts, large block sizes, and a large number of processors. Numerical examples using the nonlinear Schrodinger equation demonstrate the analytical conclusions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Fast Multilevel Solvers for a Class of Discrete Fourth Order Parabolic Problems

In this paper, we study fast iterative solvers for the solution of fourth order parabolic equations discretized by mixed finite element methods. We propose to use consistent mass matrix in the discretization and use lumped mass matrix to construct efficient preconditioners. We provide eigenvalue analysis for the preconditioned system and estimate the convergence rate of the preconditioned GMRes method. Furthermore, we show that these preconditioners only need to be solved inexactly by optimal multigrid algorithms. Our numerical examples indicate that the proposed preconditioners are very efficient and robust with respect to both discretization parameters and diffusion coefficients. We also investigate the performance of multigrid algorithms with either collective smoothers or distributive smoothers when solving the preconditioner systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗
DOE OSTI2015

Speeding up biomolecular interactions by molecular sledding

In numerous biological processes associations involve a protein with its binding partner, an event that is preceded by a diffusion-mediated search bringing the two partners together. Often hindered by crowding in biologically relevant environments, three-dimensional diffusion can be slow and result in long bimolecular association times. Moreover, the initial association step between two binding partners often represents a rate-limiting step in biotechnologically relevant reactions. We also demonstrate the practical use of an 11-a.a. DNA-interacting peptide derived from adenovirus to reduce the dimensionality of diffusional search processes and speed up associations between biological macromolecules. We functionalize binding partners with the peptide and demonstrate that the ability of the peptide to one-dimensionally diffuse along DNA results in a 20-fold reduction in reaction time. We also show that modifying PCR primers with the peptide sled enables significant acceleration of standard PCR reactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Charged particle dynamics in the presence of non-Gaussian Lévy electrostatic fluctuations

Full orbit dynamics of charged particles in a 3-dimensional helical magnetic field in the presence of -stable Levy electrostatic fluctuations and linear friction modeling collisional Coulomb drag is studied via Monte Carlo numerical simulations. The Levy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space resulting from intermittent electrostatic turbulence. The probability distribution functions of energy, particle displacements, and Larmor radii are computed and showed to exhibit a transition from exponential decay, in the case of Gaussian fluctuations, to power law decay in the case of Levy fluctuations. The absolute value of the power law decay exponents are linearly proportional to the Levy index. Furthermore, the observed anomalous non-Gaussian statistics of the particles' Larmor radii (resulting from outlier transport events) indicate that, when electrostatic turbulent fluctuations exhibit non-Gaussian Levy statistics, gyro-averaging and guiding centre approximations might face limitations and full particle orbit effects should be taken into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗