Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion equations”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Subgrid Approach for Modeling Microtopography Effects on Overland Flow

Microtopography, or heterogeneities in the elevation across scales much smaller than the domain of interest, plays a critical role in surface water retention, surface/subsurface interactions, and runoff. Resolving microtopographic influences on flow requires high–resolution simulations, which are computationally demanding even when considering the surface system in isolation and even more so when surface flow is one component in integrated simulations that couple surface flow with unsaturated subsurface flow. There is thus significant motivation for models that allow the effects of subgrid microtopography to be better represented. Subgrid models modify coarsened models to capture the microtopography–induced nonlinear effects on hydrologic processes. We present a subgrid model that alters the water storage and flow terms in the diffusion wave equation for surface flow. Stochastically generated microtopography with strongly contrasting spatial structure, high–resolution digital elevation maps from a polygonal tundra site on the North Slope of Alaska and a hummocky microtopography from a field site in Northern Minnesota are used to assess the accuracy and applicability of the subgrid model to disparate landscapes. Approaches for determining subgrid model parameters are tested and simulation results using the subgrid model are compared to benchmark fine–scale simulations and to coarse simulations that ignore microtopography. Our findings confirm that a properly parameterized subgrid model greatly improves the coarse–scale representation of hydrographs and total water content in the system. Furthermore, using the polygonal tundra example, we propose and test a strategy for moving to application–relevant spatial scales by combining microtopography classification and a few fine–scale simulations on small subdomains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observations of Offshore Internal Boundary Layers

The growth of the marine internal boundary layer (MIBL, height h i ) with the shore-normal distance x, is a topic of continuing interest because of its applications in coastal dispersion, offshore wind farm siting, coastal air-sea fluxes and in evaporative ducting. Available data on MIBL are only scarce, given its difficult to measure the variability of coastal winds. During Coupled Air-Sea Processes and Electromagnetic Research (CASPER) campaigns an array of instrumentation was deployed to measure offshore spatial variability and its effect on electromagnetic wave propagation. Meteorological sensors (flux towers and remote sensing) were deployed along the coast of Point Mugu, California, on a research vessel and FLoating Instrument Platform (FLIP) provided surface layer and boundary layer observations. In this article, measurements from multiple remote sensing instruments, such as synchronized triple Doppler lidars, small boat operations with tethered lifting system, and radiosondes provide a holistic view of the MIBL growth and its spatial variability from the coast. A convective and stable MIBL observed during two intensive operating period days showed distinct growth characteristics off the coast of Point-Mugu. During stable stratified atmospheric conditions, an MIBL was observed at least up to 47 km from the coast. The growth of the MIBL within the nearshore adjustment zone was influenced by surrounding atmospheric, oceanographic, and topographic conditions. A new parameterization scheme is developed based on advection-diffusion balance equations accounting for upstream turbulence and compared with h i observations from Doppler lidar and profiles from small boat. An evaluation of existing IBL theories was also conducted.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2020

Extended families of critical and stationary droplets for nonequilibrium phase transitions in spatially discrete bistable systems

Bistable nonequilibrium systems are realized in catalytic reaction-diffusion processes, biological transport and regulation, spatial epidemics, etc. Behavior in spatially continuous formulations, described at the mean-field level by reaction-diffusion type equations (RDEs), often mimics that of classic equilibrium van der Waals type systems. When accounting for noise, similarities include a discontinuous phase transition at some value, p eq , of a control parameter, p , with metastability and hysteresis around p eq . For each p , there is a unique critical droplet of the more stable phase embedded in the less stable or metastable phase which is stationary (neither shrinking nor growing), and with size diverging as p → p eq . Spatially discrete analogs of these mean-field formulations, described by lattice differential equations (LDEs), are more appropriate for some applications, but have received less attention. It is recognized that LDEs can exhibit richer behavior than RDEs, specifically propagation failure for planar interphases separating distinct phases. Herein, we show that this feature, together with an orientation dependence of planar interface propagation also deriving from spatial discreteness, results in the occurrence of entire families of stationary droplets. The extent of these families increases approaching the transition and can be infinite if propagation failure is realized. In addition, there can exist a regime of generic two-phase coexistence where arbitrarily large droplets of either phase always shrink. Such rich behavior is qualitatively distinct from that for classic nucleation in equilibrium and spatially continuous nonequilibrium systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An Illumination- and Temperature-Dependent Analytical Model for Copper Indium Gallium Diselenide (CIGS) Solar Cells

In this study, we present a physics-based analytical model for copper indium gallium diselenide (CIGS) solar cells that describes the illumination- and temperature-dependent current-voltage (I-V) characteristics and accounts for the statistical shunt variation of each cell. The model is derived by solving the drift-diffusion transport equation so that its parameters are physical and, therefore, can be obtained from independent characterization experiments. The model is validated against CIGS I-V characteristics as a function of temperature and illumination intensity. This physics-based model can be integrated into a large-scale simulation framework to optimize the performance of solar modules, as well as predict the long-term output yields of photovoltaic farms under different environmental conditions.

14 SOLAR ENERGY↗
DOE OSTI2016

Effective conditions for the neutron flux density at axial boundaries of the core

Analytical expressions for elements of the triangular matrix of effective conditions at the boundary of the core with a multiregion reflector are derived in the few-group diffusion approximation. The developed technique is verified using the example of fuel assemblies of a light-water reactor with an intermediate neutron spectrum.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Charon v. 2.0

Charon is used to solve the drift-diffusion PDE Equations associated with semi-conductor physics. it is used to stimulate the operation of semi-conductor devices such as transistors and diodes. Additionally, the ode includes the functionality, but no the restricted data, to model these devices under the effects of neutron radiation.

Hennigan, GaryL.↗
DOE OSTI2018

Theoretical Analysis of the Entrainment–Mixing Process at Cloud Boundaries. Part I: Droplet Size Distributions and Humidity within the Interface Zone

Abstract The problem of a complex entrainment–mixing process is analyzed by solving a diffusion–evaporation equation for an open region in the vicinity of the cloud–dry air interface. Upon normalization the problem is reduced to a one-parametric one, the governing parameter being the potential evaporation parameter R proportional to the ratio of saturation deficit in the dry air to the available liquid water content in the cloud air. As distinct from previous multiple studies analyzing mixing within closed adiabatic volumes, we consider a principally nonstationary problem that never leads to a homogeneous equilibrium state. It is shown that at R < −1 the cloud edge shifts toward the cloud; that is, the cloud dissipates due to mixing with dry air, and the cloud volume decreases. If R > −1, the cloud edge shifts outside; that is, the mixing leads to an increase in the cloud volume. The time evolution of droplet size distribution and its moments, as well as the relative humidity within the expanding cloud–dry air interface, are calculated and analyzed. It is shown that the values of the mean volume radii rapidly decrease within the interface zone in the direction away from the cloud, indicating significant changes in the cloud edge microstructure. Scattering diagrams plotted for the cloud edge agree well with high-frequency in situ measurements, corroborating the reliability of the proposed approach. It is shown that the humidity front moves toward dry air faster than the front of liquid water content. As a result, the mixing leads to formation of a humid air shell around the cloud. The widths of the interface zone and humid shell are evaluated.

Pinsky, Mark↗
DOE OSTI2019

Theoretical Analysis of the Entrainment–Mixing Process at Cloud Boundaries. Part II: Motion of Cloud Interface

The entrainment–mixing process leads to formation of a cloud–environment interface zone, consisting of the cloud dilution zone, where the liquid water content decreases as the distance from the cloud interior increases, and a humid shell with higher humidity in comparison with the dry environment. The time evolution of the location and width of these zones is analyzed by solving the diffusion–evaporation equation for the open region in the vicinity of the cloud–dry air interface. Upon normalization, the problem is reduced to a one-parametric one, the governing parameter being the potential evaporation parameter R < 0, which is proportional to the ratio of saturation deficit in the dry air to the available liquid water content in the cloud air. It is shown that the widths of the dilution zone and the humid shell increase with time. At R < −1, the interface between the dilution zone and the humid shell (i.e., the cloud edge) moves toward the cloud interior (i.e., the cloud dissipates). If R > −1, the cloud edge moves outward; that is, the mixing leads to an increase in the cloud volume. It is shown that the relative humidity remains high within most of the cloud dilution zone, and the effective radius changes only slightly, especially in growing clouds, where the effective radius changes by less than 10% within more than 80% of the cloud dilution zone. The strong decrease in the LWC is accompanied by a corresponding decrease in the droplet concentration. The study is illustrated by means of two examples of clouds: a small cumulus (Cu) typical of the trade wind zone and a deep convective cloud typical of the monsoon period over India. It is shown that while the small Cu tends to dissipate, the deep convective cloud expands. The characteristic sizes of the dilution zone and the humid shell are several tens of meters for the small Cu and several hundred meters for the deep convective cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE

Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Adjoint sensitivity analysis and data assimilation for verification of dry storage cask contents

Dry cask storage is a method for interim storage of spent fuel assemblies which contain fissile isotopes of uranium and plutonium. These can present a proliferation concern and consequently there is a need for non-destructive testing methods to verify a dry cask's contents for proliferation protection. We present an application of adjoint sensitivity analysis and data assimilation to a multigroup diffusion model of dry cask storage. Adjoint sensitivity analysis allows the efficient calculation of sensitivities for use in data assimilation to calibrate imprecisely known parameter values and data consistency tests to detect diversion scenarios. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Analyzing the 2019 Chilean social outbreak: Modelling Latin American economies

In this work, we propose a quantitative model for the 2019 Chilean protests. We utilize public data for the consumer price index, the gross domestic product, and the employee and per capita income distributions as inputs for a nonlinear diffusion-reaction equation, the solutions to which provide an in-depth analysis of the population dynamics. Specifically, the per capita income distribution stands out as a solution to the extended Fisher-Kolmogorov equation. According to our results, the concavity of employee income distribution is a decisive input parameter and, in contrast to the distributions typically observed for Chile and other countries in Latin America, should ideally be non-negative. Based on the results of our model, we advocate for the implementation of social policies designed to stimulate social mobility by broadening the distribution of higher salaries.

Curilef, Sergio (ORCID:0000000283183206)↗
DOE OSTI2016

Report on Evaluation, Verification, and Assessment of Porosity Migration Model in Fast Reactor MOX Fuel

Abstract This report documents the progress of simulating pore migration in ceramic (UO 2 and mixed oxide or MOX) fuel using BISON. The porosity field is treated as a function of space and time whose evolution is governed by a custom convection-diffusion-reaction equation (described here) which is coupled to the heat transfer equation via the temperature field. The porosity is initialized to a constant value at every point in the domain, and as the temperature (and its gradient) are increased by application of a heat source, the pores move up the thermal gradient and accumulate at the center of the fuel in a time-frame that is consistent with observations from experiments. There is an inverse dependence of the fuel’s thermal conductivity on porosity (increasing porosity decreases thermal conductivity, and vice-versa) which is also accounted for, allowing the porosity equation to couple back into the heat transfer equation. Results from an example simulation are shown to demonstrate the new capability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Simultaneous Momentum and Position Measurement and the Instrumental Weyl-Heisenberg Group

The canonical commutation relation, [Q,P]=iℏ, stands at the foundation of quantum theory and the original Hilbert space. The interpretation of P and Q as observables has always relied on the analogies that exist between the unitary transformations of Hilbert space and the canonical (also known as contact) transformations of classical phase space. Now that the theory of quantum measurement is essentially complete (this took a while), it is possible to revisit the canonical commutation relation in a way that sets the foundation of quantum theory not on unitary transformations but on positive transformations. This paper shows how the concept of simultaneous measurement leads to a fundamental differential geometric problem whose solution shows us the following. The simultaneous P and Q measurement (SPQM) defines a universal measuring instrument, which takes the shape of a seven-dimensional manifold, a universal covering group we call the instrumental Weyl-Heisenberg (IWH) group. The group IWH connects the identity to classical phase space in unexpected ways that are significant enough that the positive-operator-valued measure (POVM) offers a complete alternative to energy quantization. Five of the dimensions define processes that can be easily recognized and understood. The other two dimensions, the normalization and phase in the center of the IWH group, are less familiar. The normalization, in particular, requires special handling in order to describe and understand the SPQM instrument.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Solar Energetic Particle Acceleration by a Shock Wave Accompanying a Coronal Mass Ejection in the Solar Atmosphere

Solar energetic particle acceleration by a shock wave accompanying a coronal mass ejection (CME) is studied. The description of the accelerated particle spectrum evolution is based on the numerical calculation of the diffusive transport equation with a set of realistic parameters. The relation between the CME and shock speeds, which depend on the initial CME radius, is determined. Depending on the initial CME radius, its speed, and the magnetic energy of the scattering Alfvén waves, the accelerated particle spectrum is established 10–60 minutes from the beginning of CME motion. The maximum energies of particles reach 0.1–10 GeV. The CME radii of 3–5 R{sub ⊙} and the shock radii of 5–10 R{sub ⊙} agree with observations. The calculated particle spectra agree with the observed ones in events registered by ground-based detectors if the turbulence spectrum in the solar corona significantly differs from the Kolmogorov one.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Numerical investigation of microbial quorum sensing under various flow conditions

Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.

59 BASIC BIOLOGICAL SCIENCES↗