Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion equation”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An experimental comparison of a space-time multigrid method with PFASST for a reaction-diffusion problem

We consider two parallel-in-time approaches applied to a (reaction) diffusion problem, possibly non-linear. In particular, we consider PFASST (Parallel Full Approximation Scheme in Space and Time) and space-time multigrid strategies. For both approaches, we start from an integral formulation of the continuous time dependent problem. Then, a collocation form for PFASST and a discontinuous Galerkin discretization in time for the space-time multi-grid are employed, resulting in the same discrete solution at the time nodes. Strong and weak scaling of both multilevel strategies are compared for varying orders of the temporal discretization. Moreover, we investigate the respective convergence behavior for non-linear problems and highlight quantitative differences in execution times

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Extension of modified power method to two-dimensional problems

In this study, the generalized modified power method was extended to two-dimensional problems. A direct application of the method to two-dimensional problems was shown to be unstable when the number of requested eigenmodes is larger than a certain problem dependent number. The root cause of this instability has been identified as the degeneracy of the transfer matrix. In order to resolve this instability, the number of sub-regions for the transfer matrix was increased to be larger than the number of requested eigenmodes; and a new transfer matrix was introduced accordingly which can be calculated by the least square method. The stability of the new method has been successfully demonstrated with a neutron diffusion eigenvalue problem and the 2D C5G7 benchmark problem. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Approximate static condensation algorithm for solving multi-material diffusion problems on meshes non-aligned with material interfaces

In this study, we describe a new algorithm for solving multi-material diffusion problem when material interfaces are not aligned with the mesh. In this case interface reconstruction methods are used to construct approximate representation of interfaces between materials. They produce so-called multi-material cells, in which materials are represented by material polygons that contain only one material. The reconstructed interface is not continuous between cells. Finally, we suggest the new method for solving multi-material diffusion problems on such meshes and compare its performance with known homogenization methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Salt and ion transport in a series of crosslinked AMPS/PEGDA hydrogel membranes

Salt permeability and ionic conductivity are critical properties for membranes in water purification and energy applications. Both properties depend on individual ion sorption and diffusion coefficients, which are significantly influenced by polymer chemical and physical parameters such as fixed charge concentration and membrane water content. However, systematic studies connecting polymer structure to ion transport properties are still lacking. In this study, a series of uncharged and charged membranes were synthesized using poly(ethylene glycol) diacrylate (PEGDA) as a cross-linker and 2-acrylamido-2-methyl-1-propanesulfonic acid (AMPS) as a charged monomer. Membrane fixed charge concentration and water uptake were systematically varied by adjusting AMPS content in the pre-polymerization mixture. Salt sorption and permeability coefficients and ionic conductivity of these membranes were measured as a function of NaCl solution concentration (0.01–1 M). Combining the solution-diffusion model and Nernst-Planck equation, individual ion diffusion coefficients were calculated. Experimental Na + diffusion coefficients for all materials were well described by the Mackie and Meares tortuosity model, highlighting the strong influence of water content on ion diffusivity in both uncharged and charged polymers. Model predictions for Cl – diffusion coefficients agree reasonably well with experimental values, with some deviations occurring in more highly charged membranes. Finally, this discrepancy might result from interactions not captured by the Mackie and Meares model (e.g., fixed charge-ion interactions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗