Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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 505 records · Page 28

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗

Nonlinear theory of electron neutralization waves in ions beams with dissipation

An analytical theory of nonlinear neutralization waves generated by injection of electrons from a grid in the direction of a homogeneous ion beam of uniform velocity and infinite extension is presented. The electrons are assumed to interact with the ions through the self-consistent space charge field and by strong collective interactions, while diffusion in the pressure gradient is disregarded (zero-temperature approximation). The associated nonlinear boundary-value problem is solved in closed form by means of a von Mises transformation. It is shown that the electron gas moves into the ion space in the form of a discontinuous neutralization wave, which exhibits a periodic field structure (incomplete neutralization). This periodic wave structure is damped out by intercomponent momentum transfer - i.e., after a few relaxation lengths a quasi-neutral plasma results.

Wilhelm, H. E.↗

Cosmic-ray transport theory and out-of-the-ecliptic exploration

The reasons for studying cosmic-ray transport theory are summarized and the fundamentally three-dimensional nature of the process is pointed out. It is shown that observations in the solar ecliptic plane cannot unambiguously test transport theories since the solutions to the transport equations depend critically on boundary conditions and variation of parameters such as diffusion tensor out of the ecliptic. Sample calculations (Fokker-Planck coefficient) are shown which illustrate the problem. It is concluded that out-of-the-ecliptic observations are essential to further test transport theory.

Jokipii, J. R.↗

Evolution of weak disturbances in inert binary mixtures

The evolution of weak disturbances in inert binary mixtures is determined for the one-dimensional piston problem. The interaction of the dissipative and nonlinear mechanisms is described by Burgers' equation. The binary mixture diffusion mechanisms enter as an additive term in an effective diffusivity. Results for the impulsive motion of a piston moving into an ambient medium and the sinusoidally oscillating piston are used to illustrate the results and elucidate the incorrect behavior pertaining to the associated linear theory.

Rasmussen, M. L.↗

Development of low cost contacts to silicon solar cells

The results of the second phase of the program of developing low cost contacts to silicon solar cells using copper are presented. Phase 1 yielded the development of a plated Pd-Cr-Cu contact system. This process produced cells with shunting problems when they were heated to 400 C for 5 minutes. Means of stopping the identified copper diffusion which caused the shunting were investigated. A contact heat treatment study was conducted with Pd-Ag, Ci-Ag, Pd-Cu, Cu-Cr, and Ci-Ni-Cu. Nickel is shown to be an effective diffusion barrier to copper.

Tanner, D. P.↗

A theoretical study of the effects of interacting grain boundaries on electron-beam-induced currents

The influence of several grain boundaries on the electron-beam-induced current in polycrystalline silicon is investigated. The simple geometry contemplated here, consisting of two plane parallel grain boundaries possessing different effective, albeit constant, surface recombination velocities, intersecting a planar junction at right angles allows for an exact solution of the corresponding boundary value problem. It is shown that if the distance between the grain boundaries w becomes less than the diffusion length for minority carriers L within the grain, the grain boundary with the smaller effective recombination velocity becomes less and less visible the smaller the ratio w/L.

Luke, K. L.↗

Solvability condition for needle crystals at large undercooling in a nonlocal model of solidification

It is explicitly shown that, in a realistic model of diffusion-controlled dendritic solidification, Ivantsov's continuous family of steady-state needle crystals is destroyed by the addition of surface tension. The starting point is in the exact integro-differential equation for the one-sided model, in two dimensions, in a moving frame of reference. In the limit of large undercooling, where the range of the diffusion field is much smaller than the radius of curvature of the tip of the needle, this problem is reduced to a linear, inhomogeneous differential equation of infinite order. A solvability condition for this equation is derived and it is shown that solutions cease to exist for arbitrarily small but finite isotropic surface tension.

Caroli, B.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Obliquity variation in a Mars climate evolution model

The existence of layered terrain in both polar regions of Mars is strong evidence supporting a cyclic variation in climate. It has been suggested that periods of net deposition have alternated with periods of net erosion in creating the layered structure that is seen today. The cause for this cyclic climatic behavior is variation in the annually averaged latitudinal distribution of solar insolation in response to obliquity cycles. For Mars, obliquity variation leads to major climatological excursion due to the condensation and sublimation of the major atmospheric constituent, CO2. The atmosphere will collapse into the polar caps, or existing caps will rapidly sublimate into the atmosphere, dependent upon the polar surface heat balance and the direction of the change in obliquity. It has been argued that variations in the obliquity of Mars cause substantial departures from the current climatological values of the surface pressure and the amount of CO2 stored in both the planetary regolith and polar caps. In this new work we have modified the Haberle et al. model to incorporate variable obliquity by allowing the polar and equatorial insolation to become functions of obliquity, which we assume to vary sinusoidally in time. As obliquity varies in the model, there can be discontinuities in the time evolution of the model equilibrium values for surface pressure, regolith, and polar cap storage. The time constant, tau r, for the regolith to find equilibrium with the climate is estimated--depending on the depth, thermal conductivity, and porosity of the regolith--between 10(exp 4) and 10(exp 6) yr. Thus, using 2000-yr timesteps to move smoothly through the 0.1250 m.y. obliquity cycles, we have an atmosphere/regolith system that cannot be assumed in equilibrium. We have dealt with this problem by limiting the rate at which CO2, can move between the atmosphere and regolith, mimicking the diffusive nature and effects of the temperature and pressure waves, by setting the time rate of change of regolith storage proportional to the difference between equilibrium storage and current storage.

Tyler, D.↗

Filtration Combustion in Smoldering and SHS

Smolder waves and SHS (self-propagating high-temperature synthesis) waves are both examples of filtration combustion waves propagating in porous media. Smoldering combustion is important for the study of fire safety. Smoldering itself can cause damage, its products are toxic and it can also lead to the more dangerous gas phase combustion which corresponds to faster propagation at higher temperatures. In SHS , a porous solid sample, consisting of a finely ground powder mixture of reactants, is ignited at one end. A high temperature thermal wave, having a frontal structure, then propagates through the sample converting reactants to products. The SHS technology appears to enjoy a number of advantages over the conventional technology, in which the sample is placed in a furnace and "baked" until it is "well done". The advantages include shorter synthesis times, greater economy, in that the internal energy of the reactions is employed rather than the costly external energy of the furnace, purer products, simpler equipment and no intrinsic limitation on the size of the sample to be synthesized as exists in the conventional technology. When delivery of reactants through the pores to the reaction site is an important aspect of the combustion process, it is referred to as filtration combustion. The two types of filtration combustion have a similar mathematical formulation, describing the ignition, propagation and extinction of combustion waves in porous media. The goal in each case, however, is different. In smoldering the desired goal is to prevent propagation, whereas in SHS the goal is to ensure propagation of the combustion wave, leading to the synthesis of desired products. In addition, the scales in the two areas of application differ. Smoldering generally occurs at lower temperatures and propagation velocities than in SHS nevertheless, the two applications have much in common so that what is learned fit make application can be used to advantage in the other. In porous media, melting often occurs ahead of the propagating combustion wave. In certain cases there is so much melting that the porous solid structure is destroyed, e.g., by melting and a suspension arises, consisting of a liquid bath containing solid particles and/or gas bubbles. The resulting combustion wave is referred to as a liquid flame. We have considered a number of problems involving filtration combustion. Here, we describe four such studies: (A) rapid buoyant filtration combustion waves; (B) diffusion driven combustion waves; (C) rapidly propagating liquid flames in gravitational fields; and (D) gas-phase influence on liquid flames in gravitational fields.

Matkowsky, Bernard J.↗

L Band Brightness Temperature from Forest: Comparison of Approximate Techniques

In this paper, three approximate physical microwave radiometry models have been used to calculate brightness temperatures from a forest canopy at L-band. These models are (1) tau-omega model (zero order scattering approximation to radiative transfer equations), (2) successive order of scattering model up to first order (first order scattering approximation to the radiative transfer equations), and (3) Peak technique utilizing the active solution obtained from the Distorted Born Approximation (DBA). These models are physically-based and treat vegetation as a layer of discrete scatterers over a rough surface. Vegetation components within the canopy are represented by canonical shapes such as dielectric discs and cylinders. The tau-omega model is based on a zero-order solution to the radiative transfer (RT) equations. The model ignores scattering except for the effect of the scatterers in the attenuation of the emission through the vegetation. Application of the tau-omega model to data acquired during airborne and ground-based campaigns over the years has solidified scientific understanding of microwave interactions with different landscapes. In particular, shrubland, grasslands, agricultural crops, and light to moderate vegetation have been investigated. Its applicability to areas with a significant tree fraction is unknown. The first order scattering model is based on an iterative solution of the RT equation up to the first order. The first order solution is obtained by substituting the zeroth-order solution into the scattering source term and then solving the resulting radiative transfer equations. This formulation adds a new scattering term to the tau-omega model. It represents emission by particles in the layer and emission by the ground that is scattered once by particles in the layer. The resulting model represents an improvement over the standard zero-order solution (the tau-omega model) since it accounts for the scattered vegetation and ground radiation that can have a pronounced effect on the observed brightness temperature. The third model is based on the Peake formulation in conjunction with the DBA. The procedure for calculation of forest emission is accomplished by first calculating the bistatic scattering cross section for each type of scatterer, then by using the DBA to calculate specular albedo of the ground and the diffused albedo of the layer. Once the albedos are determined, Peake s principle relating active and passive problems can be used to determine the effective emissivity of the forest layer.

Kurum, Mehmet↗

VIIRS On-Orbit Optical Anomaly - Investigation, Analysis, Root Cause Determination and Lessons Learned

A gradual, but persistent, decrease in the optical throughput was detected during the early commissioning phase for the Suomi National Polar-Orbiting Partnership (SNPP) Visible Infrared Imager Radiometer Suite (VIIRS) Near Infrared (NIR) bands. Its initial rate and unknown cause were coincidently coupled with a decrease in sensitivity in the same spectral wavelength of the Solar Diffuser Stability Monitor (SDSM) raising concerns about contamination or the possibility of a system-level satellite problem. An anomaly team was formed to investigate and provide recommendations before commissioning could resume. With few hard facts in hand, there was much speculation about possible causes and consequences of the degradation. Two different causes were determined as will be explained in this paper. This paper will describe the build and test history of VIIRS, why there were no indicators, even with hindsight, of an on-orbit problem, the appearance of the on-orbit anomaly, the initial work attempting to understand and determine the cause, the discovery of the root cause and what Test-As-You-Fly (TAYF) activities, can be done in the future to greatly reduce the likelihood of similar optical anomalies. These TAYF activities are captured in the lessons learned section of this paper.

Iona, Glenn↗

Survey of the Relationship Between Theory and Experiment for Convective Heat Transfer From Rocket Combustion Gases

The problem of predicting the distribution of heat flux in a rocket engine is complicated by a flow field and thermodynamics so perturbed by as yet indescribable combustion effects that it continues to be necessary to utilize a simplified theoretical model against which to compare experimental results in order to gain insight. The most successful model applied to this problem is that of a compressible, turbulent, accelerating boundary layer flow with an arbitrary point of origin, and chemical reaction restricted to equilibrium recombination or diffusion-limited recombination on a catalytic wall. Analysis of this model was not possible without making several assumptions about the nature of the skin friction coefficient relationship and Reynold's analogy in an accelerating boundary layer flow. The paper summarizes the analysis of this model (together with detailed treatment of the special problems of variable properties, chemical reaction, and sources of transport properties), and describes several closed-form equations that have also been developed and applied to the problem. Described are results of several experiments with hot air flow made for the purpose of checking the validity of the assumptions in the boundary layer analysis. Most significant, however, are presentations of sample distributions of local heat flux measured in rocket engines operating with a wide range of conditions, propellants, injectors and sizes. When compared with predictions, these data show the significant range of deviation in heat flux that results from the combustion effects and which precludes a generally accurate prediction theory or correlation equation until the combustion process itself is correlated.

D R Bartz↗

Effect of shear stress on cell cultures and other reactor problems

Anchorage dependent cell cultures in fluidized beds are tested. Feasibility calculations indicate the allowed parameters and estimate the shear stresses therein. In addition, the diffusion equation with first order reaction is solved for the spherical shell (double bubble) reactor with various constraints.

Schleier, H.↗

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗