Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion model”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Modeling the low-altitude trapped radiation environment

The predictions of models of trapped energetic particle radiation at low altitudes, including the NASA trapped radiation models, are compared with low altitude measurements of energetic protons and electrons. Knowledge of the trapped energetic particle radiation at low altitudes in the vicinity of the south Atlantic magnetic anomaly is considered to be important for the analysis of the space environmental effects on space systems in this region. The long term trends in the observations are generally consistent with the NASA trapped radiation models. The discrepancies between the models and observations are currently being investigated using a 3D diffusive model and realistic atmospheric and magnetic field models. The requirements for future work and development of low altitude trapped radiation environment models are discussed.

Fung, Shing F.↗
NASA NTRS2004

Rotating Reverse Osmosis for Wastewater Reuse

Our previous work established the concept of a low-pressure rotating reverse osmosis membrane system. The rotation of the cylindrical RO filter produces shear and Taylor vortices in the annulus of the device that decrease the concentration polarization and fouling commonly seen with conventional RO filtration techniques. A mathematical model based on the film theory and the solution-diffusion model agrees well with the experimental results obtained using this first generation prototype. However, based on the model, the filtrate flux and contaminant rejection depend strongly on the transmembrane pressure. Therefore, the goal of our current work is to improve the flux of the device by increasing the transmembrane pressure by a factor of 3 to 4. In addition, the rejections for a wider variety of inorganic and organic compounds typically found in space mission wastewater are measured.

Lueptow, Richard M.↗
NASA NTRSDate not supplied

Nodal Numerical Modeling of Submerged Helium Injection in a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper presents a mathematical model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. This mathematical model was incorporated in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. The numerical model developed with GFSSP has been validated against two sets of experimental data and has been shown to predict both propellant subcooling and helium consumption to within 30% in most cases. The test data used for the model validation were taken in 1-g, but the mass diffusion model was developed to be applied in both 1-g and micro-g environments.

Cryogenics↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
DOE OSTI2021

A Flow Rate Dependent 1D Model for Thermally Stratified Hot-Water Energy Storage

Stratified tank models are used to simulate thermal storage in applications such as residential or commercial hot-water storage tanks, chilled-water storage tanks, and solar thermal systems. The energy efficiency of these applications relates to the system components and the level of stratification maintained during various flow events in the tank. One-dimensional (1D) models are used in building energy simulations because of the short computation time but often do not include flow-rate dependent mixing. The accuracy of 1D models for plug flow, plug flow with axial conduction, and two convection eddy-diffusivity models were compared with experimental data sets for discharging a 50-gal residential tank and recharging the tank with hot water from an external hot-water source. A minimum and maximum relationship for the eddy diffusivity factor were found at Re <2100 and >10,000 for recirculation of hot water to the top of the tank and vertical tubes inletting cold water at the bottom. The root mean square error decreased from >4 °C to near 2 °C when considering flow-based mixing models during heating, while the exponential decay of the eddy diffusion results in a root mean square error reduction of 1 °C for cone-shaped diffusers that begin to relaminarize flow at the inlet.

1D flow model↗
NASA NTRS2000

A Mathematical Model of Diffusion-Limited Gas Bubble Dynamics in Tissue with Varying Diffusion Region Thickness

A three-region mathematical model of gas bubble dynamics has been shown suitable for describing diffusion-limited dynamics of more than one bubble in a given volume of extravascular tissue. The model is based on the dynamics of gas exchange between a bubble and a well-stirred tissue region through an intervening unperfused diffusion region previously assumed to have constant thickness and uniform gas diffusivity. As a result, the gas content of the diffusion region remains constant as the volume of the region increases with bubble growth, causing dissolved gas in the region to violate Henry's law. Earlier work also neglected the relationship between the varying diffusion region volume and the fixed total tissue volume, because only cases in which the diffusion region volume is a small fraction of the overall tissue volume were considered. We herein extend the three-region model to correct these theoretical inconsistencies by allowing both the thickness and gas content of the diffusion region to vary during bubble evolution. A postulated difference in gas diffusivity between an infinitesimally thin layer at the bubble surface and the remainder of the diffusion region leads to variation in diffusion region gas content and thickness during bubble growth and resolution. This variable thickness, differential diffusivity (VTDD) model can yield bubble lifetimes considerably longer than those yielded by earlier three-region models for given model and decompression parameters, and meets a need for theoretically consistent but relatively simple bubble dynamics models for use in studies of decompression sickness (DCS) in human subjects, Keywords: decompression sickness, gas diffusion in tissue, diffusivity

Srinivasan, R. Srini↗
NASA NTRS1972

Model for radon diffusion through the lunar regolith.

Description of a model for radon diffusion through the lunar regolith in which the atom migrates by random walk. The regolith is represented by a system of randomly oriented baffles in which the mean distance which the atom travels between two collisions takes on the role of a mean free path. The effective mean time between two collisions depends on two entities: the actual mean time-of-flight and the mean sticking time on grain surfaces for one collision. The latter depends strongly on the temperature and the heat of adsorption of radon on regolith materials. Both the mean free path as well as the heat of adsorption are either poorly known or unknown for the lunar regolith; hence these quantities are treated as free parameters. Because of the greatly different mean lifetimes against radioactive decay of Rn219, Rn220, and Rn222, the regolith acts as a powerful 'filter' for these species. Rn222 escape is significant (32%) even for a mean free path of 1 micron, a heat of adsorption of 7.0 kcal/mole and a regolith depth of 4 m. Calculations of radon escape from a 4 m thick regolith, using mean free paths of 1, 10, and 80 microns and heats of adsorption of 4.0, 5.2, and 7.0 kcal/mole show that the Rn222/Rn220 escape ratio can be as small as 7.7 and as large as, or larger than 47. The small value of 7.7 is of particular interest, because it is nearly equal to the escape ratio inferred by Turkevich et al. (1970) from their Surveyor 5 results.

Friesen, L. J.↗
DOE OSTI2024

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗
DOE OSTI2022

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗
NASA NTRS1991

Analysis of redox additive-based overcharge protection for rechargeable lithium batteries

The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection, has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. Digital simulation of the overcharge experiment leads to numerical representation of the potential transients, and estimate of the influence of diffusion coefficient and interelectrode distance on the transient attainment of the steady state during overcharge. The model has been experimentally verified using 1,1-prime-dimethyl ferrocene as a redox additive. The analysis of the experimental results in terms of the theory allows the calculation of the diffusion coefficient and the formal potential of the redox couple. The model and the theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.

Narayanan, S. R.↗
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Electron transport and the magnetic field in the Crab Nebula

Observational data obtained at radio, optical, and X-ray frequencies on the size of the emission region in the Crab Nebula are used to examine various energy-transport models. It is shown that the data cannot be fitted by diffusion models, even when the energy-injection region is extended. A more favorable comparison is made with an extended source model in which bulk motion is the transport mechanism. Assuming a point source of electrons and a uniform nebular magnetic field, a perpendicular magnetic-field component of less than 0.0008 gauss is determined. Combining this with the lower limit obtained from gamma-ray observations, a perpendicular field component of between 0.0005 and 0.0008 gauss is obtained. A break in the size dependence is predicted at about 12 microns and 45 keV. Above 45 keV, the size should decrease as the inverse square root of proton energy.

Wolff, R. S.↗
NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain Corrected Satellite Sea Surface Salinity Using the GMAO S2S Forecast System

During the past years, we have seen that the La Nina to El Nino transition has had a significant global socio-economic impact and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help to improve the initialization of the thermocline, while satellite Sea Surface Temperature (SST) aids in constraining surface heat-fluxes, leading to improved subseasonal to seasonal forecasts of the coupled system. However, few studies have focused on improving the fresh-water flux and near-surface density and mixing through assimilation of satellite sea surface salinity (SSS). For expediency, the few ocean models that do assimilate SSS bias-correct the satellite SSS data to normalize towards the near-surface Argo data. However, in rainy regions, where buoyant water sits as a fresh lens at the surface, this assumption is likely inadequate. In previous work, we have shown that adjusting SSS data to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 has improved the near-surface density and mixed layer depth, leading to deeper thermocline and improved the NINO3.4 SST forecasts. Now we utilize the Soil Moisture/Ocean Salinity, Rain Corrected (SMOS_RC) SSS product provided by the Centre Aval de Traitement des données SMOS (CATDS CPDC) to represent the Sb more accurately at first model layer (in our case 5 m). Rather than using a diffusivity model as with RIM, SMOS_RC relies on an observed relationship between the spatial heterogeneity of SMOS SSS and instantaneous rain rate (RR) (Supply et al., 2020). In order to test the impact of SMOS_RC versus SMOS, we compare two reanalyses over the period 2014 to 2021. For both reanalysis experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is (i.e., with the fresh bias), and a separate reanalysis is performed assimilating the SMOS_RC data. We assess the impact for near-surface and subsurface dynamics within ocean reanalyses by validating against observations and explore how SSS assimilation (SMOS versus SMOS_RC) impacts dynamical ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (GEOS S2S-3, Molod et al., 2020, Hackert et al., 2023). We will show that improved SSS estimates and near-surface density and mixing led to more accurate coupled air/sea interaction and better ENSO forecasts.

Eric Hackert↗
DOE OSTI2022

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

A two component model of the diurnal variations in the thermospheric composition

A self-consistent, two-dimensional, two-component diffusion model is presented to describe the diurnal component of thermosphere dynamics. In the regions where oxygen and helium are the minor constituents, diffusion is shown to increase the amplitudes of these constituents and to advance the times of their diurnal maxima by several hours with respect to the gas temperature and the major species. Both effects are in basic agreement with satellite observations. The diffusion process is relatively insignificant for temperature, total mass density, and wind fields. The magnitude of the temperature-density phase delay is discussed with respect to its dependence on the electron density distribution.

Mayr, H. G.↗