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
NASA NTRS1999

Radiative Smoothing in Clouds at Transparent and Absorbing Wavelengths

For absorbing and transparent wavelengths, we discuss the effect of horizontal solar radiative fluxes in clouds on the accuracy of a conventional plane-parallel radiative transfer calculations for a single pixel, known as the Independent Pixel Approximation (IPA). We address the question of correlations between horizontal fluxes, IPA accuracies and radiative smoothing. By smoothing we understand a radiative transfer process whereby radiation does not follow the small-scale fluctuations of cloud structure, producing much smoother radiation fields. The scale eta that characterizes this process is called "radiative smoothing scale." We relate radiative smoothing to the photon's horizontal displacement that characterizes a "spot" of reflected light associated with a point source. We generalize the "spot-size" estimate derived for conservative scattering using the diffusion theory to the case of non-conservative scattering. For reflected light, theoretical results are confirmed with numerical simulations. The radiative smoothing scale eta is a critical value where IPA effectively breaks down; for scales smaller than TI, real radiation field are much smoother than their IPA counterparts for the same cloud structure. In addition to the estimate of il for absorbing wavelengths, we show that: (1) with more absorption, the scale break determined by eta in a log-log plot of wavenumber spectra moves towards smaller scales and (2) the smaller eta the flatter the small-scale slope which means less radiative smoothing, thus more accuracy in the IPA reflection.

Marshak, A.↗
NASA NTRS2015

Badhwar-O'Neill 2014 Galactic Cosmic Ray Flux Model Description

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) model is based on GCR measurements from particle detectors. The model has mainly been used by NASA to certify microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BON14 model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model also incorporates an empirical time delay function to account for the lag of the solar activity to reach the boundary of the heliosphere. This technical paper describes the most recent improvements in parameter fits to the BON model (BON14). Using a comprehensive measurement database, it is shown that BON14 is significantly improved over the previous version, BON11.

O'Neill, P. M.↗
NASA NTRS2015

NASA Galactic Cosmic Radiation Environment Model: Badhwar - O'Neill (2014)

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) flux model has been used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a comprehensive database of GCR measurements from various particle detectors to determine boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

Golge, S.↗
NASA NTRS2016

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

shock wave/boundary-layer interaction↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2017

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A channel flow cell with double disk electrodes for oxygen electroreduction study at elevated temperatures and pressures: Theory

We report a new double-disk-electrode channel flow cell has been designed, constructed and tested. A theoretical treatment and implicit finite difference numerical simulations are presented, which permit the computation of the hydrodynamic behavior of this cell under steady state conditions, i.e., the diffusion limited current at the working electrode and the collection efficiency of the detector. The simulated results are consistent with those measured experimentally, indicating that our simulations appropriately describe the hydrodynamic behavior of this cell. The oxygen reduction reaction (ORR) on a carbon-supported nanoparticle Pt catalyst at both room temperature and 110 °C were preliminarily studied with this cell.

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

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

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

Large-eddy simulation of a multi-injection flame in a diesel engine environment using an unsteady flamelet/progress variable approach

In this work, large-eddy simulations (LESs) are conducted for a multiple-injection flame in a diesel engine environment using an unsteady flamelet/progress variable (UFPV) approach in which differential diffusion is considered. The suitability of the UFPV tabulation approach is first evaluated through a priori analyses using the state-of-the-art direct numerical simulation (DNS) dataset. Both the instantaneous data and the conditional values for the major and minor species' mass fractions are compared between the UFPV and the DNS. The comparisons show that the proposed UFPV tabulation approach can give good predictions for the multiple-injection flame at different injection phases. While the gas temperature and major species mass fractions can be accurately predicted with or without differential diffusion being considered in the UFPV flamelet library, the prediction accuracy for the highly diffusive species (e.g., hydrogen) in the main injection phase can be noticeably improved when differential diffusion is taken into account. The fully coupled LES/UFPV simulations show that the overall structure of the multiple-injection flame can be predicted, and the conditional thermo-chemical values are close to the filtered DNS dataset. The reasons for the remaining discrepancies found in the a priori analyses and the a posteriori simulations using the UFPV approach are analyzed.

Wen, Xu (ORCID:0000000331334292)↗
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 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 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↗
NASA NTRS2012

The Extended Pulsar Magnetosphere

We present the structure of the 3D ideal MHD pulsar magnetosphere to a radius ten times that of the light cylinder, a distance about an order of magnitude larger than any previous such numerical treatment. Its overall structure exhibits a stable, smooth, well-defined undulating current sheet which approaches the kinematic split monopole solution of Bogovalov 1999 only after a careful introduction of diffusivity even in the highest resolution simulations. It also exhibits an intriguing spiral region at the crossing of two zero charge surfaces on the current sheet, which shows a destabilizing behavior more prominent in higher resolution simulations. We discuss the possibility that this region is physically (and not numerically) unstable. Finally, we present the spiral pulsar antenna radiation pattern.

Constantinos, Kalapotharakos↗