Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigating the crust of neutron stars with neural-network quantum states

An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm −3 and 0.1 fm −3 , matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.

Nuclear astrophysics↗
DOE OSTI2025

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗
DOE OSTI2025

Boundary Effects in the Diffusion of New Products on Cartesian Networks

The Role of Boundaries in the Spreading of Solar Peer effects by neighbors play a key role in the spreading of residential solar. Thus, people are more likely to install a solar system on their roof if some of their neighbors have already done so. Because people who live near the municipality boundary have fewer neighbors, does this imply that they are less likely to adopt solar? In “Boundary Effects in the Diffusion of New Products on Cartesian Networks,” Fibich, Levin, and Gillingham analyze this problem analytically using the Bass model on two-dimensional networks and empirically using data on installations of solar systems. They show that boundaries have a significant impact on the adoption of residential units near the municipality boundary. Their effect on the aggregate adoption in the municipality, however, is negligible.

Business & Economics↗
DOE OSTI2020

Inpainting Galactic Foreground Intensity and Polarization Maps Using Convolutional Neural Networks

The Deep Convolutional Neural Networks (DCNNs) have been a popular tool for image generation and restoration. In this work, we applied DCNNs to the problem of inpainting non-Gaussian astrophysical signal, in the context of Galactic diffuse emissions at the millimetric and submillimetric regimes, specifically Synchrotron and Thermal Dust emissions. Both signals are affected by contamination at small angular scales due to extragalactic radio sources (the former) and dusty star-forming galaxies (the latter). Furthermore, we compare the performance of the standard diffusive inpainting with that of two novel methodologies relying on DCNNs, namely Generative Adversarial Networks and Deep-Prior. We show that the methods based on the DCNNs are able to reproduce the statistical properties of the ground-truth signal more consistently with a higher confidence level.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Single crystal tubes of beta alumina

Edge-defined, film-fed growth process allows both tubular shapes and single crystallinity to be achieved. Beta alumina in single crystal form makes possible membranes with improved conductivities. Single crystal membranes also eliminate problems associated with electrical short circuiting of membrane due to possible sodium metal diffusion.

Stormont, R. W.↗
NASA NTRS1973

NASA/MSFC multilayer diffusion models and computer program for operational prediction of toxic fuel hazards

The NASA/MSFC multilayer diffusion models are discribed which are used in applying meteorological information to the estimation of toxic fuel hazards resulting from the launch of rocket vehicle and from accidental cold spills and leaks of toxic fuels. Background information, definitions of terms, description of the multilayer concept are presented along with formulas for determining the buoyant rise of hot exhaust clouds or plumes from conflagrations, and descriptions of the multilayer diffusion models. A brief description of the computer program is given, and sample problems and their solutions are included. Derivations of the cloud rise formulas, users instructions, and computer program output lists are also included.

Dumbauld, R. K.↗
NASA NTRS1979

A new approach to the correlation of boundary layer mass transfer rates with thermal diffusion and/or variable properties

A rational approach to the correlation of boundary layer mass transport rates, applicable to many commonly encountered laminar flow conditions with thermal diffusion and/or variable properties, is outlined. The correlation scheme builds upon already available constant property blowing/suction solutions by introducing appropriate correction factors to account for the additional ('pseudo' blowing and source) effects identified with variable properties and thermal diffusion. Applications of the scheme to the particular laminar boundary layer mass transfer problems considered herein (alkali and transition metal compound vapor transport) indicates satisfactory accuracy up to effective blowing factors equivalent to about one third of the 'blow off' value. As a useful by-product of the variable property correlation, we extend the heat-mass transfer analogy, for a wide range of Lewis numbers, to include variable property effects.

Srivastava, R.↗
NASA NTRS1986

Modeling of infrared flux spectra from disk-shaped interstellar dust clouds

Although there is growing observational and theoretical evidence for many disk-shaped objects of astrophysical interest, spherical geometry is assumed in most radiative transfer models. Recently we generalized the quasi-diffusion method developed by Leung (1975, 1976) for spherical geometry to solve the problem of scattering, absorption, and reemission by dust grains in a medium of disk geometry. The method is applicable to a variety of astronomical sources whose dynamics are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, bipolar-flow molecular clouds, accretion disks and disk galaxies. Using this technique and realistic grain opacities we construct theoretical models to determine self-consistently the dust temperature distribution and infrared emission from disk-shaped, quiescent dark globules heated externally by the intersteller radiation field. The effects of the following parameters on the temperature structure and the emergent spectrum are studied: grain type (graphite and silicate), optical depth, density inhomogeneity, and degree of disk flattening. The disk models are characterized by an aspect ratio and an optical depth at 0.55 microns. For inhomogeneous models, a gaussian density distribution is assumed such that the ratio of central to surface density is 100. To study the effects of source geometry, we also compare results for models with spherical and disk geometry. In this case both models have the same radius and central optical depth, the disk models having a 1:1 aspect ratio. While the dust temperature distributions in the two cases are very similar, the emergent flux for the disk model depends sensitively on the viewing angle. For clouds which are unresolved, one would expect, since the thermal emission is isotropic in the neighborhood of an emitting grain, and since the emission (in the far infrared) is optically thin, that the emergent flux spectrum should be characteristic only of the dust temperature, and independent of viewing angle. However, because of the lack of complete symmetry and the resulting radiation anisotropy, this is found not to be the case for the disk models. Angle is that it implies a large uncertainty in estimating the radiation dust mass.

Spagna, G. F., Jr.↗
NASA NTRS1986

CONDIF - A modified central-difference scheme with unconditional stability and very low numerical diffusion

This paper describes a method, called CONDIF, which retains the essential nature of the central difference scheme (CDS) but eliminates the over- and under-shoots. It modifies the CDS by introducing a controlled amount of numerical diffusion based on local gradients. The results obtained for a number of test problems show that the CONDIF is significantly more accurate than the hybrid scheme when the Peclet number is very high and the flow is at large angles to the grid. In the worst case, it yields results similar to those of the hybrid scheme. Though the CONDIF is unconditionally stable, under some conditions the rate of convergence deteriorates as the grid Peclet number is increased. Efforts are currently underway to develop a faster converging variation of CONDIF.

Runchal, Akshai K.↗
NASA NTRS1987

Approximation techniques for domain identification in two-dimensional parabolic systems under boundary observations

The identification of the geometrical structure of the system boundary for a two-dimensional diffusion system is considered. The domain identification treated is converted into an optimization problem based on a fit-to-data criterion and theoretical convergence results for approximate identification techniques are discussed. Results of numerical experiments to demonstrate the efficiency of the theoretical ideas are reported.

Banks, H. T.↗
NASA NTRS1989

Exponential Finite-Difference Technique

Report discusses use of explicit exponential finite-difference technique to solve various diffusion-type partial differential equations. Study extends technique to transient-heat-transfer problems in one dimensional cylindrical coordinates and two and three dimensional Cartesian coordinates and to some nonlinear problems in one or two Cartesian coordinates.

Handschuh, Robert F.↗
NASA NTRS1992

Onset of salt finger convection in a gravity gradient

The effects of gravity level elevation in a centrifuge on the onset of salt finger convection in a double-diffusive system are examined. It is shown that the steady onset of salt finger convection in a double-diffusive layer is exactly analogous to the small-gap Taylor-Couette problem. The results presented by Chandrasekhar (1961) for the Taylor-Couette problem are extended to cover a wider range of application in the present case.

Chen, C. F.↗
NASA NTRS1992

Electromagnetic scattering in clouds

Techniques used to explain the nature of the optical effects of clouds on the light produced by lightning include a Monte Carlo simulation, an equivalent medium approach, and methods based on Boltzmann transport theory. A cuboidal cloud has been considered using transform methods and a diffusion approximation. Many simplifying assumptions have been used by authors to make this problem tractable. In this report, the cloud will have a spherical shape and its interior will consist of a uniform distribution of identical spherical water droplets. The source will be modeled as a Hertz dipole, electric or magnetic, inside or outside the cloud. An impulsive source is used. Superposition may be employed to obtain a sinusoid within an envelope which describes a lightning event. The problem is investigated by transforming to the frequency domain, obtaining Green's functions, and then using the Cagniard-DeHoop method to symbolically recover the time domain solution.

Solakiewicz, Richard↗
NASA NTRS1958

Application of the Method of Coordinate Perturbation to Unsteady Duct Flow

The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.

Himmel, Seymour C.↗
NASA NTRS1961

Near-Field and Far-Field Noise Measurements for a Blowdown-Wind-Tunnel Supersonic Exhaust Jet having about 475,000 Pounds of Thrust

Near-field and far-field noise surveys were made of the supersonic The exhaust jet of the Langley 9- by 6-foot thermal structures tunnel. The jet had a thrust rating of approximately 475,000 pounds. The sound power radiated was found to be about 3.6 x 10(exp 6) watts, and on an acoustical-mechanical efficiency basis this value is in reasonable agreement with data for smaller supersonic jets and for rocket engines of other investigations. Octave-band analyses of the near-field noise show that the maximum sound pressure levels in the low-frequency bands are greatest downstream, whereas maximum sound pressure levels in the high-frequency bands were greatest near the jet exit. A comparison of near-field noise measurements is made with data previously obtained for rocket engines. Noise survey measurements of the original jet are compared with similar data obtained after the addition of a 97-foot-long exit diffuser section, and an example of the application of this facility to the problem of acoustic environmental testing of a large space capsule is cited.

Mayes, William H.↗
DOE OSTI2020

A fast implicit solver for semiconductor models in one space dimension

Several different approaches are proposed for solving fully implicit discretizations of a simplified Boltzmann-Poisson system with a linear relaxation-type collision kernel. This system models the evolution of free electrons in semiconductor devices under a low-density assumption. At each implicit time step, the discretized system is formulated as a fixed-point problem, which can then be solved with a variety of methods. A key algorithmic component in all the approaches considered here is a recently developed sweeping algorithm for Vlasov-Poisson systems. A synthetic acceleration scheme has been implemented to accelerate the convergence of iterative solvers by using the solution to a drift-diffusion equation as a preconditioner. The performance of four iterative solvers and their accelerated variants has been compared on problems modeling semiconductor devices with various electron mean-free-path.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

PolyPole-1: An accurate numerical algorithm for intra-granular fission gas release

This paper describes the development of a new numerical algorithm (called PolyPole-1) to efficiently solve the equation for intra-granular fission gas release in nuclear fuel. The work was carried out in collaboration with Politecnico di Milano and Institute for Transuranium Elements. The PolyPole-1 algorithms is being implemented in INL's fuels code BISON code as part of BISON's fission gas release model. The transport of fission gas from within the fuel grains to the grain boundaries (intra-granular fission gas release) is a fundamental controlling mechanism of fission gas release and gaseous swelling in nuclear fuel. Hence, accurate numerical solution of the corresponding mathematical problem needs to be included in fission gas behaviour models used in fuel performance codes. Under the assumption of equilibrium between trapping and resolution, the process can be described mathematically by a single diffusion equation for the gas atom concentration in a grain. In this work, we propose a new numerical algorithm (PolyPole-1) to efficiently solve the fission gas diffusion equation in time-varying conditions. The PolyPole-1 algorithm is based on the analytic modal solution of the diffusion equation for constant conditions, with the addition of polynomial corrective terms that embody the information on the deviation from constant conditions. The new algorithm is verified by comparing the results to a finite difference solution over a large number of randomly generated operation histories. Furthermore, comparison to state-of-the-art algorithms used in fuel performance codes demonstrates that the accuracy of the PolyPole-1 solution is superior to other algorithms, with similar computational effort. Finally, the concept of PolyPole-1 may be extended to the solution of the general problem of intra-granular fission gas diffusion during non-equilibrium trapping and resolution, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗