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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Relativistic electron precipitation events driven by electromagnetic ion-cyclotron waves

We adopt a canonical approach to describe the stochastic motion of relativistic belt electrons and their scattering into the loss cone by nonlinear EMIC waves. The estimated rate of scattering is sufficient to account for the rate and intensity of bursty electron precipitation. This interaction is shown to result in particle scattering into the loss cone, forming ∼10 s microbursts of precipitating electrons. These dynamics can account for the statistical correlations between processes of energization, pitch angle scattering, and relativistic electron precipitation events, that are manifested on large temporal scales of the order of the diffusion time ∼tens of minutes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Four decades of implicit Monte Carlo

In 1971, Fleck and Cummings derived a system of equations to enable robust Monte Carlo simulations of time-dependent, thermal radiative transfer problems. Denoted the “Implicit Monte Carlo” (IMC) equations, their solution remains the de facto standard of high-fidelity radiative transfer simulations. Over the course of 44 years, their numerical properties have become better understood, and accuracy enhancements, novel acceleration methods, and variance reduction techniques have been suggested. In this review, we rederive the IMC equations—explicitly highlighting assumptions as they are made—and outfit the equations with a Monte Carlo interpretation. We put the IMC equations in context with other approximate forms of the radiative transfer equations and present a new demonstration of their equivalence to another well-used linearization solved with deterministic transport methods for frequency-independent problems. We discuss physical and numerical limitations of the IMC equations for asymptotically small time steps, stability characteristics and the potential of maximum principle violations for large time steps, and solution behaviors in an asymptotically thick diffusive limit. We provide a new stability analysis for opacities with general monomial dependence on temperature. Here, we consider spatial accuracy limitations of the IMC equations and discussion acceleration and variance reduction techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Calculation of the RF Conductivity and Hall Constant of a Thin Metal Film

The kinetic problem of finding the rf conductivity and Hall constant of a thin metal film placed in a transverse steady magnetic field and a longitudinal variable electric field has been considered. It has been assumed that electrons diffusely reflect from the upper and lower surfaces of the film. No limitations have been imposed on the relationship between the thickness of the film and the electron free path. The dependences of the conductivity and Hall constant on dimensionless parameters, namely, electric field frequency, magnetic field induction, and thickness of the film, have been studied. Calculation results have been compared with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Modulation of solar flare particles and track density profiles in gas-rich meteorite grains

A solution is presented to the problem concerning the time-averaged solar flare particle flux as a function of kinetic energy and distance from the sun for a given particle injection spectrum at the sun within the framework of standard diffusion-convection-adiabatic deceleration theory with the diffusion coefficient independent of distance from the sun. Results of the calculations which give best agreement with observations at 1 AU are presented and discussed, with particular reference to their implications for gas-rich meteorites. Normalization at the orbit of earth is achieved via observed track density versus depth profiles in lunar vug crystals. It is shown that if gas-rich meteorite grains were irradiated in the asteroid belt and if source and modulation parameters have changed little since irradiation, the track density should be 'harder' than the lunar vug profile by about 0.2-0.3 in the index. Quantitative estimation of solar flare particle exposure ages is discussed.

Lee, M. A.↗
NASA NTRS1978

Experimental clean combustor program: Phase 3: Turbulence measurement addendum

Airflow turbulence parameters were measured in the high pressure, high temperature flow stream leaving the compressor of an operating gas turbine engine. A water cooled hot film turbulence probe was used to determine the turbulence intensity and length scale in the compressor exit flow stream of the CF6-50 engine. Data were obtained only at idle operating conditions. At engine power levels above idle conditions, durability problems and erratic data readings were encountered with the turbulence measurement probes. Turbulence measurements were made at three radial immersions at a point 15.5 cm downstream of the compressor outlet guide vanes in the compressor exit diffuser. The passage height at this point is 5.54 cm. Data reduction was accomplished with a time-data fast Fourier transform (FFT) system. This system finds the power spectral density distribution (PSD) of a large number of data samples using a direct Fourier transform algorithm and finds the autocorrelation parameter for these data by doing an FFT analysis of the PSD curves for a series of time delay intervals.

Taylor, J. R.↗
NASA NTRS1979

A throat-bypass stability-bleed system using relief valves to increase the transient stability of a mixed-compression inlet

A stability-bleed system was installed in a YF-12 flight inlet that was subjected to internal and external airflow disturbances in the NASA Lewis 10 by 10 foot supersonic wind tunnel. The purpose of the system is to allow higher inlet performance while maintaining a substantial tolerance (without unstart) to internal and external disturbances. At Mach numbers of 2.47 and 2.76, the inlet tolerance to decreases in diffuser-exit corrected airflow was increased by approximately 10 percent of the operating-point airflow. The stability-bleed system complemented the terminal-shock-control system of the inlet and did not show interaction problems. For disturbances which caused a combined decrease in Mach number and increase in angle of attack, the system with valves operative kept the inlet started 4 to 28 times longer than with the valves inoperative. Hence, the stability system provides additional time for the inlet control system to react and prevent unstart. This was observed for initial Mach numbers of 2.55 and 2.68. For slow increase in angle of attack at Mach 2.47 and 2.76, the system kept the inlet started beyond the steady-state unstart angle. However, the maximum transient angles of attack without unstart could not be determined because wind-tunnel mechanical-stop limits for angle of attack were reached.

Neiner, G. H.↗
NASA NTRS1981

Symposium /International/ on Combustion, 18th, University of Waterloo, Waterloo, Ontario, Canada, August 17-22, 1980, Proceedings

Problems related to combustion generated pollution are explored, taking into account the mechanism of NO formation from nitrogen compounds in hydrogen flames studied by laser fluorescence, the structure and similarity of nitric oxide production in turbulent diffusion flames, the effect of steam addition on NO formation, and the formation of NO2 by laminar flames. Other topics considered are concerned with propellant combustion, fluidized bed combustion, the combustion of droplets and sprays, premixed flame studies, fire studies, and flame stabilization. Attention is also given to coal flammability, chemical kinetics, turbulent combustion, soot, coal combustion, the modeling of combustion processes, combustion diagnostics, detonations and explosions, ignition, internal combustion engines, combustion studies, and furnaces.

Source record↗
NASA NTRS1987

Spacecraft particulate contaminant redistribution

This paper describes analyses that were performed in support of the Hubble Space Telescope (HST)) particulate contamination control effort. The specific problems addressed include extension of available particle removal data to launch acoustic and random vibration conditions, development of an engineering model for transport of suspended particulates by airflow and in the presence of vehicle acceleration, turbulent diffusion, and migration of particulates over vibrating surfaces, and integration of the various models into a code that could be used to generate contamination level estimates for the HST primary mirror and other critical surfaces for the HST mission phases. The overall redistribution calculations were made assuming a specified initial contaminant distribution in terms of MIL STD 1246A levels, and using predicted vibration data for the various HST surfaces and mission phases. As expected, the effects of airflow were found to be significant, particularly for the larger particles. Particles smaller than about 20 microns did not participate appreciably in the redistribution.

Klavins, A.↗
NASA NTRS2002

The AstroBiology Explorer (ABE) MIDEX Mission Concept: Using Infrared Spectroscopy to Identify Organic Molecules in Space

One of the principal means by which organic compounds are detected and identified in space is by infrared spectroscopy. Past IR telescopic and laboratory studies have shown that much of the carbon in the interstellar medium (ISM) is in complex organic species but the distribution, abundance and evolutionary relationships of these materials are not well understood. The Astrobiology Explorer (ABE) is a MIDEX mission concept designed to conduct IR spectroscopic observations to detect and identify these materials and address outstanding problems in astrobiology, astrochemistry, and astrophysics. ABE's core science program includes observations of planetary nebulae and stellar outflows, protostellar objects, Solar System objects, and galaxies, and lines of sight through dense molecular clouds and the diffuse ISM. ABE is a cryogenically-cooled 60 cm diameter space telescope equipped with 3 cross-dispersed R-2000 spectrometers that share a single common slit. Each spectrometer measures one spectral octave and together cover the entire 2.5-20 micron region simultaneously. The spectrometers use state-of-the-art InSb and Si:As 1024x1024 pixel detectors. ABE would operate in a heliocentric, Earth drift-away orbit and have a core science mission lasting approximately 1.5 years. ABE is currently under study at NASA's Ames Research Center in collaboration with Ball Aerospace and Technologies Corp.

Sandford, Scott A.↗
NASA NTRS2014

Capillary Flows Along Open Channel Conduits: The Open-Star Section

Capillary rise in tubes, channels, and grooves has received significant attention in the literature for over 100 years. In yet another incremental extension of such work, a transient capillary rise problem is solved for spontaneous flow along an interconnected array of open channels forming what is referred to as an 'open-star' section. This geometry possesses several attractive characteristics including passive phase separations and high diffusive gas transport. Despite the complex geometry, novel and convenient approximations for capillary pressure and viscous resistance enable closed form predictions of the flow. As part of the solution, a combined scaling approach is applied that identifies unsteady-inertial-capillary, convective-inertial-capillary, and visco-capillary transient regimes in a single parameter. Drop tower experiments are performed employing 3-D printed conduits to corroborate all findings.

Weislogel, Mark↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗
NASA NTRS1997

Measurement of Atomic Oxygen in Diffuse Aurora and Ion Density in the E-Region

An ion mass spectrometer (IMS) was refurbished, calibrated and supplied to the University of Colorado payload (Dr. Charles Barth, P.I.) which was launched from White Sands in September of 1993 as NASA 33.062. The nose cone failed to deploy and their were problems with the ACS so the mission was declared a failure. However, the door covering the IMS deployed and the instrument obtained data. The launch occurred shortly after a payload carrying solar x-ray detectors was launched. Thus a small portion of the Colorado payload science was salvaged; namely, the NO(+)/O2(+) ratio to compare with the measured x-ray flux. Figure I shows the NO(+) to O2(+) ratio vs. altitude. The behavior is typical of the E-region.

Sharp, William E.↗
DOE OSTI2022

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Using Diffusion Bonding in Making Piezoelectric Actuators

A technique for the fabrication of piezoelectric actuators that generate acceptably large forces and deflections at relatively low applied voltages involves the stacking and diffusion bonding of multiple thin piezoelectric layers coated with film electrodes. The present technique stands in contrast to an older technique in which the layers are bonded chemically, by use of urethane or epoxy agents. The older chemical-bonding technique entails several disadvantages, including the following: It is difficult to apply the bonding agents to the piezoelectric layers. It is difficult to position the layers accurately and without making mistakes. There is a problem of disposal of hazardous urethane and epoxy wastes. The urethane and epoxy agents are nonpiezoelectric materials. As such, they contribute to the thickness of a piezoelectric laminate without contributing to its performance; conversely, for a given total thickness, the performance of the laminate is below that of a unitary piezoelectric plate of the same thickness. The figure depicts some aspects of the fabrication of a laminated piezoelectric actuator by the present diffusion- bonding technique. First, stock sheets of the piezoelectric material are inspected and tested. Next, the hole pattern shown in the figure is punched into the sheets. Alternatively, if the piezoelectric material is not a polymer, then the holes are punched in thermoplastic films. Then both faces of each punched piezoelectric sheet or thermoplastic film are coated with a silver-ink electrode material by use of a silkscreen printer. The electrode and hole patterns are designed for minimal complexity and minimal waste of material. After a final electrical test, all the coated piezoelectric layers (or piezoelectric layers and coated thermoplastic films) are stacked in an alignment jig, which, in turn, is placed in a curved press for the diffusion-bonding process. In this process, the stack is pressed and heated at a specified curing temperature and pressure for a specified curing time. The pressure, temperature, and time depend on the piezoelectric material selected. At the end of the diffusion-bonding process, the resulting laminated piezoelectric actuator is tested to verify the adequacy of the mechanical output as a function of an applied DC voltage.

Sager, Frank E.↗
DOE OSTI2016

Preliminary Compendium of the Initial Conditions and Evolution of the Fundamental Interactions of SiC with Ag and Pd

Tri-structural isotropic (TRISO) particles are the fuel of choice for very-high temperature reactors. At the core, these spherical particles consist of an inner fuel kernel of UO{sub 2}, UC{sub 2} or a combination of both, which is coated radially outward by successive layers of low density pyrolytic graphite, an inner high density pyrolytic graphite layer, silicon carbide (SiC) and an outer high density pyrolytic graphite layer. SiC is the main diffusion barrier against the release of fission products due to its hardness and high melting point. Irradiation testing of TRISO particles indicates that trace amounts of metallic fission products can damage or potentially escape through the barrier layer. Surrogate experimental studies have previously been performed by various methods seeking to shed light on and explain the diffusion behavior of Ag and Pd. Although some of these experiments have yielded realistic diffusion coefficients the discrepancy of results between different groups is such that the essential details of the diffusive behavior are still to be fully understood. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. Theoretical simulations of the diffusion of Ag through SiC have been exhaustive but typically yield bulk diffusion coefficients that due to the nature of bulk diffusion phenomena (i.e. slow diffusion) the diffusion coefficients are calculated to be orders of magnitude lower than what experimental observations suggest. This in part is due to the difficulty in simulating realistic models that include various polymorphs of SiC, and the dynamical structural transformations that SiC undergoes in the TRISO environment (which are also to be explained). The result is a particularly challenging problem to solve given that the initial conditions of the interactions of Ag and Pd with SiC are difficult to address, both in theoretical models and surrogate experiments. In order to improve these models the approach taken here is to consider fundamental interactions between Ag and SiC that take place in short time scales by varying the initial environmental conditions, i.e. annealing temperature and time, initial amounts of the metallic species (M) Ag, Pd in SiC/M/SiC multilayers. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. In order to do so we have synthesized SiC/Ag/SiC and SiC/Pd/SiC multilayers of various thicknesses by PLD and characterized them using SEM/EDS. Initial results show that upon annealing Ag agglomerates and ruptures the SiC layer to be released. In the case of Pd a strong chemical interaction is deduced from the slow interlayer diffusion and limited release. Further studies using photoelectron spectroscopy and RHEED will be carried out in order obtain detailed information on the chemical and physical transitions that undergo upon annealing of these systems. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2017

Surface tension models for a multi-material ALE code with AMR

A number of surface tension models have been implemented in a 3D multi-physics multi-material code, ALE–AMR, which combines Arbitrary Lagrangian Eulerian (ALE) hydrodynamics with Adaptive Mesh Refinement (AMR). ALE–AMR is unique in its ability to model hot radiating plasmas, cold fragmenting solids, and most recently, the deformation of molten material. The surface tension models implemented include a diffuse interface approach with special numerical techniques to remove parasitic flow and a height function approach in conjunction with a volume-fraction interface reconstruction package. These surface tension models are benchmarked with a variety of test problems. In conclusion, based on the results, the height function approach using volume fractions was chosen to simulate droplet dynamics associated with extreme ultraviolet (EUV) lithography.

42 ENGINEERING↗