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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗
DOE OSTI2017

Rosseland and Flux Mean Opacities for Compton Scattering

Rosseland mean opacity plays an important role in theories of stellar evolution and X-ray burst models. In the high-temperature regime, when most of the gas is completely ionized, the opacity is dominated by Compton scattering. Our aim here is to critically evaluate previous works on this subject and to compute the exact Rosseland mean opacity for Compton scattering over a broad range of temperature and electron degeneracy parameter. We use relativistic kinetic equations for Compton scattering and compute the photon mean free path as a function of photon energy by solving the corresponding integral equation in the diffusion limit. As a byproduct we also demonstrate the way to compute photon redistribution functions in the case of degenerate electrons. We then compute the Rosseland mean opacity as a function of temperature and electron degeneracy and present useful approximate expressions. We compare our results to previous calculations and find a significant difference in the low-temperature regime and strong degeneracy. We then proceed to compute the flux mean opacity in both free-streaming and diffusion approximations, and show that the latter is nearly identical to the Rosseland mean opacity. We also provide a simple way to account for the true absorption in evaluating the Rosseland and flux mean opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Propagation of Solar Energetic Particles in Three-dimensional Interplanetary Magnetic Fields: Radial Dependence of Peak Intensities

A functional form I{sub max}(R)=kR{sup −α}, where R is the radial distance of a spacecraft, was usually used to model the radial dependence of peak intensities I{sub max}(R) of solar energetic particles (SEPs). In this work, the five-dimensional Fokker–Planck transport equation incorporating perpendicular diffusion is numerically solved to investigate the radial dependence of SEP peak intensities. We consider two different scenarios for the distribution of a spacecraft fleet: (1) along the radial direction line and (2) along the Parker magnetic field line. We find that the index α in the above expression varies in a wide range, primarily depending on the properties (e.g., location and coverage) of SEP sources and on the longitudinal and latitudinal separations between the sources and the magnetic foot points of the observers. Particularly, whether the magnetic foot point of the observer is located inside or outside the SEP source is a crucial factor determining the values of index α. A two-phase phenomenon is found in the radial dependence of peak intensities. The “position” of the break point (transition point/critical point) is determined by the magnetic connection status of the observers. This finding suggests that a very careful examination of the magnetic connection between the SEP source and each spacecraft should be taken in the observational studies. We obtain a lower limit of R{sup −1.7±0.1} for empirically modeling the radial dependence of SEP peak intensities. Our findings in this work can be used to explain the majority of the previous multispacecraft survey results, and especially to reconcile the different or conflicting empirical values of the index α in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The Acceleration of Energetic Particles at Coronal Shocks and Emergence of a Double Power-law Feature in Particle Energy Spectra

We present numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that the location on the shock where the high-energy particle intensity is the largest, depends on the energy of the particles and on time. The acceleration of particles to more than 100 MeV mainly occurs in the shock-streamer interaction region, due to perpendicular shock geometry and the trapping effect of closed magnetic fields. A comparison of the particle spectra to that in a radial magnetic field shows that the intensity at 100 MeV (200 MeV) is enhanced by more than one order (two orders) of magnitude. This indicates that the streamer-like magnetic field can be an important factor in producing large solar energetic particle events. We also show that the energy spectrum integrated over the simulation domain consists of two different power laws. Further analysis suggests that it may be a mixture of two distinct populations accelerated in the streamer and open field regions, where the acceleration rate differs substantially. Our calculations also show that the particle spectra are affected considerably by a number of parameters, such as the streamer tilt angle, particle spatial diffusion coefficient, and shock compression ratio. While the low-energy spectra agree well with standard diffusive shock acceleration theory, the break energy ranges from ~1 MeV to ~90 MeV and the high-energy spectra can extend to ~1 GeV with a slope of ~2–3.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Semi-Algebraic Two Level Solver

We develop a simple semi-algebraic 2-level solver built on traditional multigrid ideas. It is designed to be easily incorporated into existing simulation software. It exhibits good convergence for many classes of challenging problems including discontinuous diffusion, convection- diffusion, and Helmholtz equations. It has built-in structure that makes it simple to generalize in several interesting directions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Evaluating foam heterogeneity

New analytical tool is available to calculate the degree of foam heterogeneity based on the measurement of gas diffusivity values. Diffusion characteristics of plastic foam are described by a system of differential equations based on conventional diffusion theory. This approach saves research and computation time in studying mass or heat diffusion problems.

Liou, D. W.↗
NASA NTRS1972

Observed effects of earth-reflected radiation and hydrogen drag on the orbital accelerations of balloon satellites.

From March through May 1967, the high-altitude balloon satellites Pageos (1966 56A) and Dash 2 (1963 30D) were both in continuous sunlight orbits. Observed variations in the orbital accelerations of both satellites are found to be due primarily to a combination of the perturbing effects of earth albedo radiation and atmospheric drag. An equation based on diffuse Lambertian reflection of sunlight from the earth predicts quite well the observed trend of orbital accelerations for both satellites due to the force of earth albedo radiation. After correcting for this effect, the remaining orbital energy changes are found to be consistent with atmospheric drag perturbations, making it possible to deduce mean exospheric densities above 2300 km. The results indicate that hydrogen concentrations during the observation period were about three times greater than that given by the U.S. Standard Atmosphere Supplements (1966).

Prior, E. J.↗
NASA NTRS1975

Wave theory of turbulence in compressible media

An acoustical theory of turbulence was developed to aid in the study of the generation of sound in turbulent flows. The statistical framework adopted is a quantum-like wave dynamical formulation in terms of complex distribution functions. This formulation results in nonlinear diffusion-type transport equations for the probability densities of the five modes of wave propagation: two vorticity modes, one entropy mode, and two acoustic modes. This system of nonlinear equations is closed and complete. The technique of analysis was chosen such that direct applications to practical problems can be obtained with relative ease.

Kentzer, C. P.↗
NASA NTRS1975

Interpretation of recent positron-electron measurements between 20 and 800 MeV

Recently measured positron and negatron spectra are discussed with regard to the problem of solar modulation. At energies above 180 MeV, the spherically symmetric Fokker-Planck equation with a diffusion coefficient proportional to particle rigidity provides reasonable fits to both the positron and total electron data. At energies below 180 MeV, the data are consistent with a continuation of the same diffusion coefficient and a local source of negatrons or with a change in the diffusion coefficient to a constant value.

Pellerin, C. J.↗
NASA NTRS1976

Cosmic-ray positron and negatron spectra between 20 and 800 MeV measured in 1974

The flux density of cosmic-ray positrons and electrons was measured by a balloon-borne spark chamber magnetic spectrometer in two flights in an attempt to study the solar modulation of galactic cosmic rays. The total electron flux was found to be about 0.03/cu m/s/sr/MeV between 70 and 800 MeV and to increase toward lower energies, while the positron flux decreased sharply toward lower energies from about 0.008/cu m/s/sr/MeV at 650 MeV and only upper limits were obtained for positrons below 200 MeV. At energies above 180 MeV, the spherically symmetrical Fokker-Planck equation with a diffusion coefficient proportional to the particle rigidity provided a reasonable fit to both positron and total electron data, while at energies below 180 MeV, the data were consistent with a continuation of the same diffusion coefficient and local source of electrons, or a change in the diffusion coefficient to a constant value.

Hartman, R. C.↗
NASA NTRS1977

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Pavish, D. L.↗
NASA NTRS1982

Estimation of atmospheric species concentrations from remote sensing data

The atmospheric species concentration distribution-estimation problem that arises from the measurement configuration of remote sensing, which involves a field of view that moves across the region, is solved by means of the partial differential equation for atmospheric diffusion and the Wiener-Hopf theory. A numerical study is also undertaken for the estimation of concentration distribution downwind of a hypothetical, continuous ground-level source of pollutants.

Omatu, S.↗
NASA NTRS1984

Calculation of vaporization rates assuming various rate determining steps: Application to the resistojet

The various steps that could control the vaporization rate of a material are discussed. These steps include the actual vaporization, flow rate of matrix gas, chemical reaction, gas diffusion, and solid state diffusion. The applicable equations have been collected from diverse appropriate sources, and their use is explained. Rate equations are derived for conditions where more than one step is rate controlling. Calculations are made for two model materials: rhenium which vaporizes congruently, and tantalum carbide which vaporizes incongruently. The case of vaporization under thermal gradient conditions is also treated. The existence of a thermal gradient in the resistojet means that the vaporization rate of a material may be only one thousandth of that predicted under isothermal conditions. Calculations show that rhenium might have a 100,000 hr lifetime at temperature in a 2500 C resistojet. Tantalum carbide would have a life of only 660 sec under similar conditions.

May, C. E.↗
NASA NTRS1983

How do cosmic rays change their energy in the solar wind?

The diffusion-convection (modulation) equation is derived directly from the Boltzmann equation on the basis of a minimum number of assumptions concerning the scattering process, among which are: (1) that the scattered particles undergo no energy change, and (2) that isotropy is an equilibrium state. It is noted that, in the event that the background plasma contains a magnetic field and the flow speeds of the plasma and scattering centers are different, additional terms arise that will modify the equations. If, moreover, the scatterers have individual motions relative to their average flow, the second-order Fermi acceleration term will appear.

Jones, F. C.↗
NASA NTRS1986

Applications of a new wall function to turbulent flow computations

A new wall function approach is developed based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients. This wall law was derived from a one-dimensional analysis of the turbulent kinetic energy equation with gradient diffusion concept employed in modeling the near-wall shear stress gradient. Numerical testing cases for the present wall functions include turbulent separating flows around an airfoil and turbulent recirculating flows in several confined regions. Improvements on the predictions using the present wall functions are illustrated. For cases of internal recirculating flows, one modification factor for improving the performance of the k-epsilon turbulence model in the flow recirculation regions is also included.

Chen, Y. S.↗
NASA NTRS1987

Turbulence models and their applications to the prediction of internal flows - A review

The paper presents a brief account of various turbulence models employed in the computation of turbulent flows, and evaluates the application of these models to selected internal flow configurations. The main conclusions of this analysis are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly (the algebraic stress model performs better in this case); and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well. Two important factors in the numerical solution of the model equations, namely false diffusion and inlet boundary conditions, are discussed. The existence of countergradient transport and its implications in turbulence modeling are examined. Finally, some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.

Nallasamy, M.↗
NASA NTRS1987

A multispacecraft study of the injection and transport of solar energetic particles

The November 22, 1977 and December 27, 1977 solar energetic particle (SEP) events observed simultaneously by the Helios and the Voyager space probes and the ISEE 1 satellite are described. The spacecraft instruments indicated flux increases for H, He, and (for November 22) for O and Fe, with the energy range of 0.5-20 MeV per nucleon. It was found that the data can be well fitted by a model based on the spherically symmetric Fokker-Planck equation including convection, diffusion, and adiabatic deceleration, assuming that the radial interplanetary diffusion coefficient increases in proportion to the heliocentric radius to the power of 0.6 and that there is an extended particle injection at the sun. This extended injection is of crucial importance for the event of November 22, where the injection functions are strongly dependent on times associated with a soft energy spectrum.

Beeck, J.↗
NASA NTRS1989

Galerkin/least-squares procedures in computational fluid dynamics

The application of Galerkin/least-squares methodology to computational fluid dynamics is reviewed. Applications to advective-diffusive systems of equations and incompressible and compressible Euler and Navier-Stokes equations are described. Space-time finite-element methodology for time-dependent problems is emphasized and the status and role of so-called 'discontinuity-capturing operators' is reviewed.

Hughes, Thomas J. R.↗