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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

The Impacts of Dry Dynamic Cores on Asymmetric Hurricane Intensification

The fundamental pathways for tropical cyclone (TC) intensification are explored by considering axisymmetric and asymmetric impulsive thermal perturbations to balanced, TC-like vortices using the dynamic cores of three different nonlinear numerical models. Attempts at reproducing the results of previous work, which used the community WRF Model, revealed a discrepancy with the impacts of purely asymmetric thermal forcing. The current study finds that thermal asymmetries can have an important, largely positive role on the vortex intensification, whereas other studies find that asymmetric impacts are negligible. Analysis of the spectral energetics of each numerical model indicates that the vortex response to asymmetric thermal perturbations is significantly damped in WRF relative to the other models. Spectral kinetic energy budgets show that this anomalous damping is primarily due to the increased removal of kinetic energy from the vertical divergence of the vertical pressure flux, which is related to the flux of inertia-gravity wave energy. The increased kinetic energy in the other two models is shown to originate around the scales of the heating and propagate upscale with time from nonlinear effects. For very large thermal amplitudes (50 K), the anomalous removal of kinetic energy due to inertia-gravity wave activity is much smaller, resulting in good agreement between models. The results of this paper indicate that the numerical treatment of small-scale processes that project strongly onto inertia-gravity wave energy can lead to significant differences in asymmetric TC intensification. Sensitivity tests with different time integration schemes suggest that diffusion entering into the implicit solution procedure is partly responsible for the anomalous damping of energy.

Guimond, Stephen R.↗
NASA NTRSDate not supplied

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗
NASA NTRSDate not supplied

Large Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulation↗
DOE OSTI2021

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗
NASA NTRS1984

The dayside Venus ionosphere. II - Combined numerical model of ion and neutral composition above 120 km

The median vertical profile of ion densities for the dayside Venusian ionosphere obtained by the orbiter retarding potential analyzer (ORPA) is simulated by one-dimensional model calculations. The model includes both neutral and ionic chemistry, eddy and molecular diffusion for neutral constituents, and ion-plasma diffusion for ionic constituents. The electron and ion temperatures measured by the ORPA are used to calculate the plasma diffusion coefficients and scale heights for ions. The predicted O2(+) densities below about 200 km agree particularly well with observations by the ORPA, but the model values are significantly less than those measured by the orbiter ion mass spectrometer. The observed ion composition is interpreted in terms of densities of the neutral atmosphere and its composition.

Shimazaki, T.↗
NASA NTRS1999

Experimental Observations on a Low Strain Counter-Flow Diffusion Flame: Flow and Bouyancy Effects

Diffusion flames are of great interest in fire safety and many industrial processes. The counter-flow configuration provides a constant strain flow, and therefore is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in microgravity conditions have begun to explore the low strain regimes. Numerical work has shown the coupling between gas phase reaction rates, soot reaction rates, and radiation. For these programs, size, geometry and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. When the physical boundaries can not be considered infinitely far from the reaction zone discrepancies arise. A computational study that includes boundary effects and accounts for the deviations occurring when the major potential flow assumptions are relaxed was presented by Borlik et al. This development properly incorporates all heat loss terms and shows the possibility of extinction in the low strain regime. A major constraint of studying the low strain regime is buoyancy. Buoyant instabilities have been shown to have a significant effect on the nature of reactants and heat transport, and can introduce instabilities on the flow that result in phenomena such as flickering or fingering. The counter-flow configuration has been shown to provide a flame with no symmetry disrupting instabilities for inlet velocities greater than 50 mm/s. As the velocity approaches this limit, the characteristic length of the experiment has to be reduced to a few millimetres so as to keep the Rayleigh number (Ra L = (β g 0 L 3 ∆T)/(𝑎v)) below 2000. In this work, a rectangular counter-flow burner was used to study a two-dimensional counter-flow diffusion flame. Flow visualisation and Particle Image Velocimetry served to describe the nature of the stagnation plane for strain rates smaller than 100 (1/s). These experiments were conducted with a non-reacting flow. Video images of a propane air diffusion flame were used to describe the behaviour of a diffusion flame in this regime. Flame geometry and pulsation frequency are described.

J. A. Sutula↗
NASA NTRS2016

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗
NASA NTRS2016

Product of the SNPP VIIRS SD Screen Transmittance and the SD BRDF (RSB) From Both Yaw Maneuver and Regular On-Orbit Data

To assure data quality, the Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) regularly performs on-orbit radiometric calibrations of its 22 spectral bands. The primary calibration radiance source for the reflective solar bands (RSBs) is a sunlit solar diffuser (SD). During the calibration process, sunlight goes through a perforated plate (the SD screen) and then strikes the SD. The SD scattered sunlight is used for the calibration, with the spectral radiance proportional to the product of the SD screen transmittance and the SD bidirectional reflectance distribution function (BRDF). The BRDF is decomposed to the product of its value at launch and a numerical factor quantifying its change since launch. Therefore, the RSB calibration requires accurate knowledge of the product of the SD screen transmittance and the BRDF (RSB; launch time). Previously, we calculated the product with yaw maneuver data and found that the product had improved accuracy over the prelaunch one. With both yaw maneuver and regular on orbit data, we were able to improve the accuracy of the SDSM screen transmittance and the product for the solar diffuser stability monitor SD view. In this study, we use both yaw maneuver and a small portion of regular on-orbit data to determine the product for the RSB SD view.

Lei, Ning↗
DOE OSTI2016

Are oxygen limitations under recognized regulators of organic carbon turnover in upland soils?

Understanding the processes controlling organic matter (OM) stocks in upland soils, and the ability to management them, is crucial for maintaining soil fertility and carbon (C) storage as well as projecting change with time. OM inputs are balanced by the mineralization (oxidation) rate, with the difference determining whether the system is aggrading, degrading or at equilibrium with reference to its C storage. In upland soils, it is well recognized that the rate and extent of OM mineralization is affected by climatic factors (particularly temperature and rainfall) in combination with OM chemistry, mineral–organic associations, and physical protection. Here we examine evidence for the existence of persistent anaerobic microsites in upland soils and their effect on microbially mediated OM mineralization rates. Further, we corroborate long-standing assumptions that residence times of OM tend to be greater in soil domains with limited oxygen supply (aggregates or peds). Moreover, the particularly long residence times of reduced organic compounds (e.g., aliphatics) are consistent with thermodynamic constraints on their oxidation under anaerobic conditions. Incorporating (i) pore length and connectivity governing oxygen diffusion rates (and thus oxygen supply) with (ii) ‘hot spots’ of microbial OM decomposition (and thus oxygen consumption), and (iii) kinetic and thermodynamic constraints on OM metabolism under anaerobic conditions will thus improve conceptual and numerical models of C cycling in upland soils. We conclude that constraints on microbial metabolism induced by oxygen limitations act as a largely unrecognized and greatly underestimated control on overall rates of C oxidation in upland soils.

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

A comparison of Redlich-Kister polynomial and cubic spline representations of the chemical potential in phase field computations

Free energies play a central role in many descriptions of equilibrium and non-equilibrium properties of solids. Continuum partial differential equations (PDEs) of atomic transport, phase transformations and mechanics often rely on first and second derivatives of a free energy function. The stability, accuracy and robustness of numerical methods to solve these PDEs are sensitive to the particular functional representations of the free energy. In this communication we investigate the influence of different representations of thermodynamic data on phase field computations of diffusion and two-phase reactions in the solid state. First-principles statistical mechanics methods were used to generate realistic free energy data for HCP titanium with interstitially dissolved oxygen. While Redlich-Kister polynomials have formed the mainstay of thermodynamic descriptions of multi-component solids, they require high order terms to fit oscillations in chemical potentials around phase transitions. Here, we demonstrate that high fidelity fits to rapidly fluctuating free energy functions are obtained with spline functions. As a result, spline functions that are many degrees lower than Redlich-Kister polynomials provide equal or superior fits to chemical potential data and, when used in phase field computations, result in solution times approaching an order of magnitude speed up relative to the use of Redlich-Kister polynomials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Era of entropy: Synthesis, structure, properties, and applications of high-entropy materials

The field of high-entropy materials (HEMs) has emerged as a dynamic area of scientific exploration, driven by the exceptional properties arising from their compositional complexity. Encompassing both high-entropy alloys (HEAs) and high-entropy ceramics (HECs), these materials have garnered significant attention across diverse research domains. From investigations into phase evolution and mechanical characteristics to studies of ionic, electronic, and magnetic behaviors, HEMs demonstrate remarkable potential for a wide array of applications. These range from catalysis and tribology to energy storage and superconductivity. Fundamental research has shed light on crucial phenomena such as configurational entropy, lattice distortion, and sluggish diffusion. These discoveries are paving the way for materials design strategies that enable new functional tunability and resistance to application-specific harsh environments. This burgeoning field promises to revolutionize material design and performance across numerous technological sectors. Here, this special collection between Applied Physics Letters and the Journal of Applied Physics provides a timely overview of the latest research in this area. It highlights the growing interest in understanding the impact of high compositional complexity on conventional structure–process–property–performance relationships in HEMs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generative unfolding with distribution mapping

Machine learning enables unbinned, highly-differential cross section measurements. A recent idea uses generative models to morph a starting simulation into the unfolded data. We show how to extend two morphing techniques, Schrödinger Bridges and Direct Diffusion, in order to ensure that the models learn the correct conditional probabilities. This brings distribution mapping (DM) to a similar level of accuracy as the state-of-the-art conditional generative unfolding methods. Numerical results are presented with a standard benchmark dataset of single jet substructure as well as for a new dataset describing a 22-dimensional phase space of Z+2 -jets.

Butter, Anja↗
NASA NTRS1974

Analysis of chromatograph systems using orthogonal collocation

Research is generating fundamental engineering design techniques and concepts for the chromatographic separator of a chemical analysis system for an unmanned, Martian roving vehicle. A chromatograph model is developed which incorporates previously neglected transport mechanisms. The numerical technique of orthogonal collocation is studied. To establish the utility of the method, three models of increasing complexity are considered, the latter two being limiting cases of the derived model: (1) a simple, diffusion-convection model; (2) a rate of adsorption limited, inter-intraparticle model; and (3) an inter-intraparticle model with negligible mass transfer resistance.

Woodrow, P. T.↗
NASA NTRS1974

Radiation from an array of gray circular fins of trapezoidal profile

A numerical method is developed to calculate the temperature distribution and radiation heat transfer for an annular fin and tube radiator, with fins having trapezoidal profiles. All surfaces are assumed gray and to emit and reflect diffusely. Radiative interactions between adjacent fins and between the fins and tube are included. The thermal conductivity of the fin material may vary linearly with temperature. Results of a parametric study of the special case of circular fins of triangular profile having constant thermal conductivity are presented and used to optimize a fin array with respect to minimum weight.

Schnurr, N. M.↗
NASA NTRS1979

Anomalous growth of single ice crystals in solution

It is shown that major discrepancies exist between experiments and theory for ice crystal growth from solution. Accurate data, taken in a microgravity environment, approximate analytical models, and exact (probably numerical) models all are needed to advance our understanding of ice crystal growth phenomena. A new approximate semi-empirical theory is presented which predicts that a relatively sharp transition from natural convection control to diffusion control for ice growth in pure water occurs at a subcooling of about 10 C (a reduced temperature difference of about 0.125). No reliable data exist to test this prediction. The theory also predicts qualitatively the growth of ice in NaCl solution in which maxima in the growth rates are observed at various levels of subcooling.

Gill, W. N.↗
NASA NTRS1984

Turbulent compressible convection in a deep atmosphere. II - Two-dimensional results for main-sequence A5 and F0 type envelopes

In the present two-dimensional numerical study of turbulent compressible convection in the A5 and F0 main-sequence envelope types, ionization effects are included in the equation of state of the gas, and radiative transfer is modeled in a diffusive process with tabulated gas opacities. It is noted that the thermal effects of ionization significantly affect the dynamics of the flows and that an inversion of the mean density can be created and sustained in a dynamical situation. The substantial differences in the flows of A5 and F0 indicate development trends in the transition from the radiative to the convective mode of energy transport. As convection becomes more effective, the flow becomes more turbulent and the scaling effects of local scale heights become more significant.

Sofia, S.↗
NASA NTRS1986

Measurement and prediction of flow through a replica segment of a mildly atherosclerotic coronary artery of man

Pressure distributions were measured along a hollow vascular axisymmetric replica of a segment of the left circumflex coronary artery of man with mildly atherosclerotic diffuse disease. A large range of physiological Reynolds numbers from about 60 to 500, including hyperemic response, was spanned in the flows investigation using a fluid simulating blood kinematic viscosity. Predicted pressure distributions from the numerical solution of the Navier-Stokes equations were similar in trend and magnitude to the measurements. Large variations in the predicted velocity profiles occurred along the lumen. The influence of the smaller scale multiple flow obstacles along the wall (lesion variations) led to sharp spikes in the predicted wall shear stresses. Reynolds number similarity was discussed, and estimates of what time averaged in vivo pressure drop and shear stress might be were given for a vessel segment.

Back, L. H.↗