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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗
NASA NTRS1996

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗
NASA NTRS2004

Renormalization Group Tutorial

Complex physical systems sometimes have statistical behavior characterized by power- law dependence on the parameters of the system and spatial variability with no particular characteristic scale as the parameters approach critical values. The renormalization group (RG) approach was developed in the fields of statistical mechanics and quantum field theory to derive quantitative predictions of such behavior in cases where conventional methods of analysis fail. Techniques based on these ideas have since been extended to treat problems in many different fields, and in particular, the behavior of turbulent fluids. This lecture will describe a relatively simple but nontrivial example of the RG approach applied to the diffusion of photons out of a stellar medium when the photons have wavelengths near that of an emission line of atoms in the medium.

Bell, Thomas L.↗
NASA NTRS1991

The chaotic dynamics of comets and the problems of the Oort cloud

The dynamic properties of comets entering the planetary zone from the Oort cloud are discussed. Even a very slight influence of the large planets can trigger stochastic cometary dynamics. Multiple interactions of comets with the large planets produce diffusion of the parameters of cometary orbits and a mean increase in the semi-major axis of comets. Comets are lifted towards the Oort cloud, where collisions with stars begin to play a substantial role. The transport of comets differs greatly from the customary law of diffusion and noticeably alter cometary distribution.

Sagdeev, Roald Z.↗
NASA NTRS1992

The 1991 AVIRIS/POLDER experiment in Camargue, France

Airborne campaigns during the eighties provided high spectral resolution data, collected with imaging instruments such as AIS, AVIRIS, FLI, CAESI, and ISM, in order to investigate the relationship with canopy biophysical characteristics. The statistical approaches used to analyze these data do not allow investigation of the causality and the applicability of the observed correlations. Further, statistical studies demonstrated the high degree of redundancy of the spectral information amongst many others. And for retrieving vegetation biophysical characteristics, few results demonstrate the real information gain attributable to the high spectral resolution capability as compared to the use of a few wide wavelength bands. With several new imaging spectrophotometers scheduled for launch during the next 10 years (MERIS, MODIS, HIRIS), progress in the description and understanding of the mechanisms that drive the spectral variation of canopy reflectance is required. Most of these new sensor systems will also have the capability to observe the target under differing view directions. The problem of the combination and the use of the synergy between both the spectral and the directional sources of canopy reflectance variations has to be addressed. Apart from the atmospheric effects, the spectral variation of the light reflected by canopies originates from the leaves, the soil or the other vegetation elements such as branches and fruits. At leaf level, both diffuse reflectance and transmittance may be simulated by simple models, although no accurate information exists on the absorption features of the biochemicals (except water) in the 900-2500nm wavelength range. Many models mimic the directional variability of canopy reflectance at a given wavelength. Combining a leaf spectral model with a canopy directional model provides a powerful tool to analyze this problem. Some of us have initiated such a study, but our approach and theory remain to be tested using canopy data with their complexity, associated experimental error, and atmospheric effects. The main objective of the 1990 POLDER/AVIRIS experiment in Camargue was to provide a consistent data set over various canopies in order to test the applicability of the theory. The experiment, part of the 1991 MAC Europe experiment, involved simultaneous data collection using two sensors: AVIRIS and POLDER which measures the bidirectional and polarization properties of the targets at 670 and 880nm wavebands.

Baret, F.↗
NASA NTRS1992

Detection of x ray sources in PROS

The problem of detecting discrete sources in x-ray images has much in common with the problem of automatic source detection at other wavelengths. In all cases, one searches for positive brightness enhancements exceeding a certain threshold, which appear consistent with what one expects for a point source, in the presence of a (possibly) spatially variable background. Multidimensional point spread functions (e.g., dependent on detector position and photon energy) are also common. At the same time, the problem in x-ray astronomy has some unique aspects. For example, for typical x-ray exposures in current or recent observatories, the number of available pixels far exceeds the number of actual x-ray events, so Poisson, rather than Gaussian statistics apply. Further, extended cosmic x-ray sources are common, and one often desires to detect point sources in the vicinity or even within bright, diffuse x-ray emission. Finally, support structures in x-ray detectors often cast sharp shadows in x-ray images making it necessary to detect sources in a region of rapidly varying exposure. We have developed a source detection package within the IRAF/PROS environment which attempts to deal with some of the problems of x-ray source detection. We have patterned our package after the successful Einstein Observatory x-ray source detection programs. However, we have attempted to improve the flexibility and accessibility of the functions and to provide a graphical front-end for the user. Our philosophy has been to use standard IRAF tasks whenever possible for image manipulation and to separate general functions from mission-specific ones. We will report on the current status of the package and discuss future developments, including simulation tasks, to allow the user to assess detection efficiency and source significance, tasks to determine source intensity, and alternative detection algorithms.

Deponte, J.↗
DOE OSTI2018

Prediction of Cooling Curves for Controlled Unidirectional Solidification Under the Influence of Shrinkage: A Semi-analytical Approach

Quality of casting products can be significantly improved by controlling the crystal growth rate. Both dendritic and equiaxed grain structures depend strongly on crystal growth rate along the crystallographic direction. Uncontrolled solidification processes may lead to the formation of porous and columnar material with an extremely non-homogeneous composition distribution. These undesirable micro-structures can be avoided by adopting an optimal cooling rate (K/s) during the crystal growth. The present work focuses on determining suitable cooling curves to obtain desired unidirectional crystal growth rates, by implementing a semi-analytical model. Here, we defined cooling curve as the time history of temperature evolution at the cold boundary. The proposed semi-analytical heat transfer model is diffusion driven and accounts for shrinkage effect during solidification. The proposed model is found capable of predicting cooling curves for prescribed interface velocities. The model is self-sufficient in validating itself by means of solving inverse and forward problems. The cooling curve prediction by the model is validated by using it as the transient boundary condition for an existing enthalpy updating scheme-based numerical model. The proposed model is also applied to empirical data reported in literature, to validate the model against experimental results. All the validations showcased excellent fit. The proposed model will facilitate controlled unidirectional crystal growth rate by providing cooling rate (K/s) as an input for the experiments involving directional solidification.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2013

Convection of a stratified colloidal suspension

The convection of a colloidal suspension, which is a binary mixture of a carrier medium with an admixture of nanoparticles having a large positive thermal diffusion parameter, has been studied for the case of the heating of a horizontal cell from below and periodic conditions at the vertical boundaries corresponding to the experimental situation of ring channels. Bifurcation diagrams have been constructed for vibrational and monotonic regimes of the convection of the colloidal mixture. The time dependences of the maximum stream function and the stream function at a fixed point of the cell, as well as the spatial distributions of the concentration field of the colloid admixture, have been obtained. It has been shown that a stable regime of traveling waves exists in a certain region of the parameters of the problem (Boltzmann and Rayleigh numbers characterizing the gravitational stratification and intensity of the thermal effect, respectively)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exponential Runge-Kutta Parareal for non-diffusive equations

Parareal is a well-known parallel-in-time algorithm that combines a coarse and fine propagator within a parallel iteration. It allows for large-scale parallelism that leads to significantly reduced computational time compared to serial time-stepping methods. However, like many parallel-in-time methods it can fail to converge when applied to non-diffusive equations such as hyperbolic systems or dispersive nonlinear wave equations. Here, this paper explores the use of exponential integrators within the Parareal iteration. Exponential integrators are particularly interesting candidates for Parareal because of their ability to resolve fast-moving waves, even at the large stepsizes used by coarse propagators. This work begins with an introduction to exponential Parareal integrators followed by several motivating numerical experiments involving the nonlinear Schrödinger equation. These experiments are then analyzed using linear analysis that approximates the stability and convergence properties of the exponential Parareal iteration on nonlinear problems. The paper concludes with two additional numerical experiments involving the dispersive Kadomtsev-Petviashvili equation and the hyperbolic Vlasov-Poisson equation. These experiments demonstrate that exponential Parareal methods offer improved time-to-solution compared to serial exponential integrators when solving certain non-diffusive equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

High-resolution methods for preserving the sum of mass fractions: improved χ-scheme and an alternative

When high resolution convection schemes are used for discretizing chemical species mass balance equations, the mass fractions are not guaranteed to add to one. We show that a proposed remedy called χ-scheme (Darwish and Moukalled, Comput.Methods Appl.Mech. Engrg. 192 (2003): 1711) will degrade to a diffusive first-order scheme when a chemical species vanishes from the mixture, for example, because of chemical reactions. We propose an improvement to the χ-scheme to overcome this problem. Furthermore, a computationally efficient alternative scheme is proposed and evaluated with several examples, to quantify the improvements in the accuracy and the computational time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗
DOE OSTI2019

Phosphate amendments for chemical immobilization of uranium in contaminated soil

Uranium (U) contamination is a major environmental problem associated with the mining and processing of nuclear materials for both weapons and power production. When possible, in situ soil remediation techniques are preferable for reducing the risk associated with diffuse low-level U contamination. Uranium is known to form sparingly soluble phosphate compounds that persist in the environment. Therefore, batch experiments were performed to evaluate the efficacy of three phosphate amendments, hydroxyapatite (HA), sodium phytate (IP6) and sodium tripolyphosphate (TPP), to immobilize U in contaminated sediments. The amendments were added at equivalent phosphorus (P) concentrations and then equilibrated under a range of test conditions, with changes in soluble U and Ptotal monitored at pre-set time intervals. Only HA was effective at reducing the soluble U soil fraction when compared to the control, with IP6 and TPP increasing the soluble U soil fraction. After equilibration, changes in contaminant partitioning in the amended sediments were evaluated using operational extraction methods. Sequential extraction results for HA generally indicated a transfer of U from labile to more recalcitrant phases, while the results for IP6 and TPP were more ambiguous.

Uranyl↗
DOE OSTI2025

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mitigating spectral bias in neural operators via high-frequency scaling for physical systems

Neural operators have emerged as powerful surrogates for modeling complex physical problems. However, they suffer from spectral bias making them oblivious to high-frequency modes, which are present in multiscale physical systems. Therefore, they tend to produce over-smoothed solutions, which is particularly problematic in modeling turbulence and for systems with intricate patterns and sharp gradients such as multi-phase flow systems. In this work, we introduce a new approach named high-frequency scaling (HFS) to mitigate spectral bias in convolutional-based neural operators. By integrating HFS with proper variants of UNet, we demonstrate a higher prediction accuracy by mitigating spectral bias in single and two-phase flow problems. Unlike Fourierbased techniques, HFS is directly applied to the latent space, thus eliminating the computational cost associated with the Fourier transform. Additionally, we investigate alternative spectral bias mitigation through a diffusion model conditioned on neural operators. While the diffusion model integrated with the standard neural operator may still suffer from significant errors, these errors are substantially reduced when the diffusion model is integrated with a HFS-enhanced neural operator.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Multidimensional Fuel Performance Code: BISON

BISON is a finite element based nuclear fuel performance code applicable to a variety of fuel forms including light water reactor fuel rods, TRISO fuel particles, and metallic rod and plate fuel (Refs. [a, b, c]). It solves the fully-coupled equations of thermo-mechanics and species diffusion and includes important fuel physics such as fission gas release and material property degradation with burnup. BISON is based on the MOOSE framework (Ref. [d]) and can therefore efficiently solve problems on 1-, 2- or 3-D meshes using standard workstations or large high performance computers. BISON is also coupled to a MOOSE-based mesoscale phase field material property simulation capability (Refs. [e, f]). As described here, BISON includes the code library named FOX, which was developed concurrent with BISON. FOX contains material and behavioral models that are specific to oxide fuels.

Hales, Jason↗
NASA NTRS1991

A synchronous chopper mechanism for use at cryogenic temperature

A mechanically resonant synchronous chopper mechanism is described for use at cryogenic temperatures. The mechanism is a critical optical component of the Diffuse Infrared Background Experiment (DIRBE) and has been operating on orbit without incident since Nov. 1989. The requirements, electromechanical design, and testing of the mechanism are described. A description of the problems encountered and solutions implemented during the development of the mechanism is provided. Finally, a modified chopper design, which incorporated lessons learned and that has several advantages over the flight chopper design, is offered.

Hakun, Claef↗
NASA NTRS1996

Particle Acceleration in Active Galactic Nuclei

The investigation of stochastic particle acceleration through resonant interactions with plasma waves that populate the magnetosphere surrounding an accreting black hole is presented. Stochastic acceleration has been successfully applied to the problem of ion and electron energization in solar flares, and is capable of accounting for a wide range of both neutral and charged particle emissions. It is also a component in diffusive shock acceleration, since pitch-angle scattering (which is necessary for multiple shock crossings) is accompanied by diffusion in momentum space, which in turn yields a net systematic energy gain; however, stochastic energization will dominate the first-order shock process only in certain parameter regimes. Although stochastic acceleration has been applied to particle energization in the lobes of radio galaxies, its application to the central regions of AGNs (active galactic nuclei) has only recently been considered, but not in detail. We proposed to systematically investigate the plasma processes responsible for stochastic particle acceleration in black hole magnetospheres along with the energy-loss processes which impede particle energization. To this end, we calculated acceleration rates and escape time scales for protons and electrons resonating with Alfven waves, and for electrons resonating with whistlers. We also considered the "hot" topic of gamma-ray line emission from the Orion complex. We proposed that the observed gamma-ray lines are produced by energetic ions that are stochastically accelerated by cascading Alfven waves in the accretion plasma near a black hole. Related research papers that were published in journals are listed.

Miller, James A.↗