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,495 records · Page 83

RECONNECTION OF QUASI-SINGULAR CURRENT SHEETS: THE ''IDEAL'' TEARING MODE

A strong indication that fast reconnection regimes exist within resistive magnetohydrodynamics was given by the proof that the Sweet-Parker current sheet, maintained by a flow field with an appropriate inflow-outflow structure, could be unstable to a reconnecting instability which grows without bound as the Lundquist number, S, tends to infinity. The requirement of a minimum value for S in order for the plasmoid instability to kick in does little to resolve the paradoxical nature of the result. Here we argue against the realizability of Sweet-Parker current sheets in astrophysical plasmas with very large S by showing that an ''ideal'' tearing mode takes over before current sheets reach such a thickness. While the Sweet-Parker current sheet thickness scales as ∼S {sup –1/2}, the tearing mode becomes effectively ideal when a current sheet collapses to a thickness of the order of ∼S {sup –1/3}, up to 100 times thicker than S {sup –1/2}, when (as happens in many astrophysical environments) S is as large as 10{sup 12}. Such a sheet, while still diffusing over a very long time, is unstable to a tearing mode with multiple x-points: here we detail the characteristics of the instability and discuss how it may help solve the flare trigger problem and effectively initiate the turbulent disruption of the sheet.

79 ASTRONOMY AND ASTROPHYSICS↗

First principles modeling of microscopic scintillation mechanisms using a Many-Body Green’s function approach (Final Report)

In 2010, we started the first lifecycle of a project on first principles modeling of microscopic unit processes that occur during scintillation in order to bring about a capability of formulating physics-based reaction-diffusion models that can predict scintillator light yield and explain their non-proportional response to the energy of incident ionizing radiation, with the hope that this will accelerate materials discovery that can mitigate this problem. During the ensuing years we focused on investigating from first principles the elementary quantum mechanical processes leading to polaron formation and exciton self-trapping in scintillators and identify the microscopic pathways leading to their transport. The effort to derive the master equation for carrier and exciton transport in realistic scintillators materials from first principles constitutes a challenge that requires an extensive background in various areas of computational physics. Previous theoretical investigations involved significant approximations and were limited in their scope. Most notably, the standard workhorse of the computational material scientist, density-functional theory, fails in describing the localized nature of self-trapped carriers as well as the electronic structure of rare-earth dopants. Advances in both algorithms and computational resources however, have rendered a comprehensive and accurate approach feasible. This has been amply demonstrated by the results produced by us during the first lifecycle of this project, which impressed upon the reviewers of the Independent Project Review: “The Review Team agreed, and commented that ’amazing and excellent progress had been made to date.’ The researchers tried and did improve the technology to advance the state of the art and show that it was adequate for the project.” In the following, we detail the progress achieved during each year of the project. We conclude by a final section where all our main achievements are summarized.

36 MATERIALS SCIENCE↗

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗

A gallium phosphide high-temperature bipolar junction transistor

Preliminary results are reported on the development of a high temperature (350 C) gallium phosphide bipolar junction transistor (BJT) for geothermal and other energy applications. This four-layer p(+)n(-)pp(+) structure was formed by liquid phase epitaxy using a supercooling technique to insure uniform nucleation of the thin layers. Magnesium was used as the p-type dopant to avoid excessive out-diffusion into the lightly doped base. By appropriate choice of electrodes, the device may also be driven as an n-channel junction field-effect transistor. The initial design suffers from a series resistance problem which limits the transistor's usefulness at high temperatures.

Zipperian, T. E.↗

The anisotropic radiative transfer problem in optically thick, strongly magnetized plasma - A comparison of results

The accuracy of four different methods for solving the coherent radiative transfer problem in a strongly magnetized plasma is compared. Specific attention is given to the case of a semi-infinite homogeneous plasma with a magnetic field perpendicular to the surface. The four approaches evaluated were: diffusion approximation schemes; the generalized H-function method of Silant'ev (1982) the Feautrier method of Nagel (1981) and the integral equation approach of Meszares and Bonazzola (1981). It is found that the coupled diffusion approximation schemes of Nagel (1981) and Kaminker et al. (1982) provided satisfactory results and were more efficient in many cases than direct numerical methods.

Pavlov, G. G.↗

Inviscid spatial stability of a compressible mixing layer. Part 2: The flame sheet model

The results of an inviscid spatial calculation for a compressible reacting mixing layer are reported. The limit of infinitive activation energy is taken and the diffusion flame is approximated by a flame sheet. Results are reported for the phase speeds of the neutral waves and maximum growth rates of the unstable waves as a function of the parameters of the problem: the ratio of the temperature of the stationary stream to that of the moving stream, the Mach number of the moving streams, the heat release per unit mass fraction of the reactant, the equivalence ratio of the reaction, and the frequency of the disturbance. These results are compared to the phase speeds and growth rates of the corresponding nonreacting mixing layer. We show that the addition of combustion has important, and complex effects on the flow stability.

Jackson, T. L.↗

Inviscid spatial stability of a compressible mixing layer. II - The flame sheet model

The results of an inviscid spatial calculation for a compressible reacting mixing layer are reported. The limit of infinitive activation energy is taken and the diffusion flame is approximated by a flame sheet. Results are reported for the phase speeds of the neutral waves and maximum growth rates of the unstable waves as a function of the parameters of the problem: the ratio of the temperature of the stationary stream to that of the moving stream, the Mach number of the moving streams, the heat release per unit mass fraction of the reactant, the equivalence ratio of the reaction, and the frequency of the disturbance. These results are compared to the phase speeds and growth rates of the corresponding nonreacting mixing layer. We show that the addition of combustion has important and complex effects on the flow stability.

Jackson, T. L.↗

Nature of fluid flows in differentially heated cylindrical container filled with a stratified solution

Semiconductor crystals such as Hg(1-x)Cd(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. Because of the higher thermal conductivity in the melt than that in the crystal there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart (1971), Thorpe, Hutt and Soulsby (1969) have shown that under such condition a fluid will become convectively unstable as a result of different diffusivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositional segregation in the unidirectionally solidified crystals. To reach this goal, we start with a simplified problem. We study the nature of fluid flows of a stratified solution in a cylindrical container with a radial temperature gradient. The cylindrical container wall is considered to be maintained at a higher temperature than that at the center of the solution and the solution in the lower gravitational direction has higher solute concentration which decrease linearly to a lower concentration and then remain constant to the top of the solution. The sample solution is taken to be salt water.

Wang, Jai-Ching↗

Holographic interferometric tomography for reconstructing flow fields

Holographic interferometric tomography is a technique for instantaneously capturing and quantitatively reconstructing three-dimensional flow fields. It has a very useful application potential for high-speed aerodynamics. However, three major challenging tasks need to be accomplished before its practical applications. First, fluid flows are mostly unsteady or at least non repeatable. Consequently, a means for Instantaneously recording three-dimensional flow fields, that is, a simple holographic technique for simultaneously recording multi-directional projections, needs to be developed. Second, while holographic interferometry provides enormous data storage capabilities, expeditious data extraction from complicated interferograms is very important for timely near real-time applications. Third, unlike medical applications, flow tomography does not provide complete data sets but instead involves ill-posed reconstruction problems of incomplete projection and limited angular scanning. During this summer research period, new experimental techniques and corresponding hardware were developed and tested to address the above mentioned tasks. The first task was achieved by diffuser illumination. This concept allows instantaneous capture of many projections with a conventional setup for single-projection recording. For the second task, a phase-shifting technique was incorporated. This technique allows one to acquire multiple phase-stepped interferograms for a single projection and thus to extract phase information from intensity data almost at real-time. For the third task, the research that has been extensively conducted previously was utilized. In this research period, a complete experimental setup that provides the above three major capabilities was designed, built, and tested by integrating all the techniques. A simple laboratory experiment for simulating wind-tunnel testing was then conducted. A test flow was produced by employing a relatively simple device that generated a gravity-driven flow. The flow was then experimentally investigated to check the viability of the holographic interferometric tomographic technique before wind-tunnel application.

Cha, Soyoung S.↗

The morphology of faint galaxies in Medium Deep Survey images using WFPC2

First results from Hubble Space Telescope (HST) Medium Deep Survey images taken with Wide Field/Planetary Camera-2 (WFPC2) demonstrate that galaxy classifications can be reliably performed to magnitudes I814 approximately less than 22.0 in the F815W band. Published spectroscopic surveys to this depth indicate a mean redshift of bar-z approximately 0.5. We have classified over 200 galaxies in nine WFPC2 fields according to a basic morphological scheme. The majority of these faint galaxies appear to be similar to regular Hubble-sequence examples observed at low redshift. To the precision of our classification scheme, the relative proportion of spheroidal and disk systems of normal appearance is as expected from nearby samples, indicating that the bulk of the local galaxy population was in place at half the Hubble time. However, the most intriguing result is the relatively high proportion (approximately 40%) of objects which are in some way anomalous, and which may be of relevance in understanding the origin of the familiar excess population of faint galaxies established by others. These diverse objects include apparently interacting pairs whose multiple structure is only revealed with HST's angular resolution, galaxies with superluminous star-forming regions, diffuse low surface brightness galaxies of various forms, and compact galaxies. These anomalous galaxies contribute a substantial fraction of the excess counts at our limiting magnitude, and may provide insights into the 'faint blue galaxy' problem.

Griffiths, R. E.↗

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 studies (telescopic and laboratory) have demonstrated that much of the carbon in the interstellar medium (ISM) is in complex organic species of a variety of types, but the distribution, abundance, and evolutionary relationships of these materials are not well understood. The Astrobiology Explorer (ABE) is a MIDEAST mission concept designed to conduct IR spectroscopic observations to detect and identify these materials to address outstanding important problems in astrobiology, astrochemistry, and astrophysics. Systematic studies include the observation of planetary nebulae and stellar outflows, protostellar objects, Solar System Objects, and galaxies, and multiple lines of sight through dense molecular clouds and the diffuse ISM. ABE will also search for evidence of D enrichment in complex molecules in all these environments. The mission is currently under study at NASA's Ames Research Center in collaboration with Ball Aerospace and Technologies Corp. ABE is a cryogenically-cooled 60 cm diameter space telescope equipped with 3 cryogenic cross-dispersed spectrographs that share a single common slit. The 3 spectrometers each measure single spectral octaves (2.5-5, 5-10, 10-20 microns) and together cover the entire 2.5 - 20 micron region simultaneously. The spectrometers use state-of-the-art 1024x1024 pixel detectors, with a single InSb array for the 2.5-5 micron region and two Si:As arrays for the 5-10 and 10-20 micron regions. The spectral resolution is wavelength dependent but is greater than 2000 across the entire spectral range. ABE would operate in a heliocentric, Earth drift-away orbit and is designed to take maximum advantage of this environment for cooling, thermal stability, and mission lifetime. ABE would have a core science mission lasting approximately 1.5 years.

Sandford, Scott A.↗

Towards a 3D Space Radiation Transport Code

High-speed computational procedures for space radiation shielding have relied on asymptotic expansions in terms of the off-axis scatter and replacement of the general geometry problem by a collection of flat plates. This type of solution was derived for application to human rated systems in which the radius of the shielded volume is large compared to the off-axis diffusion limiting leakage at lateral boundaries. Over the decades these computational codes are relatively complete and lateral diffusion effects are now being added. The analysis for developing a practical full 3D space shielding code is presented.

Wilson, J. W.↗

Development of Multilayered Coatings for Nerva-Type Graphite Fuel Elements- An Overview

“Mid-passage corrosion” (MPC) was observed in NbC and ZrC-coated graphite (Gr)-based nuclear fuels during the NERVA/Rover programs. A new multilayered coating concept was proposed to solve this problem. The concept envisions designing a compliant multilayered coatings architecture to accommodate the thermal expansion mismatch strains between the ZrC outer coating and the Gr/(U,Zr)C fuel matrix. Additionally, the coating architecture uses diffusion barriers to minimize carbon diffusion from the substrate, and hydrogen diffusion from the flowing gas stream, into the multilayer coating. The present paper provides a brief overview of the development of multilayered coatings deposited on Gr disks, as well as within the inner channels of a19-hole hexagonal rod 152.4 mm long, by chemical vapor deposition (CVD) in proof-of concept studies. The cross-sectional microstructural observations are discussed. These trials have established that it is possible to deposit multilayered coatings by CVD in 19-hole hexagonal rods similar to the NERVA fuel rod design although the quality of the coatings and the uniformity of their thicknesses require further improvement.

Multilayer coatings, graphite, nuclear thermal pro↗

Interface design of VSOP'94 computer code for safety analysis

Today, most software applications, also in the nuclear field, come with a graphical user interface. VSOP'94 (Very Superior Old Program), was designed to simplify the process of performing reactor simulation. VSOP is a integrated code system to simulate the life history of a nuclear reactor that is devoted in education and research. One advantage of VSOP program is its ability to calculate the neutron spectrum estimation, fuel cycle, 2-D diffusion, resonance integral, estimation of reactors fuel costs, and integrated thermal hydraulics. VSOP also can be used to comparative studies and simulation of reactor safety. However, existing VSOP is a conventional program, which was developed using Fortran 65 and have several problems in using it, for example, it is only operated on Dec Alpha mainframe platforms and provide text-based output, difficult to use, especially in data preparation and interpretation of results. We develop a GUI-VSOP, which is an interface program to facilitate the preparation of data, run the VSOP code and read the results in a more user friendly way and useable on the Personal 'Computer (PC). Modifications include the development of interfaces on preprocessing, processing and postprocessing. GUI-based interface for preprocessing aims to provide a convenience way in preparing data. Processing interface is intended to provide convenience in configuring input files and libraries and do compiling VSOP code. Postprocessing interface designed to visualized the VSOP output in table and graphic forms. GUI-VSOP expected to be useful to simplify and speed up the process and analysis of safety aspects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Derivation of the Point-Kinetics Parameters for the Coupled Reactors

The design of coupled reactors, which provide two distinct spectral zones (fast and thermal), has recently drawn attention because it can provide configurations that could not be met in current thermal neutron test reactors and thus allows versatile irradiation tasks. A coupled reactor core is characterized by the neutron multiplicity of each zone as well as by the coupling between the two zones. It is desired to evaluate the integral point kinetics parameters that characterize the neutron physics of the coupled system. This summary provides a derivation of the equations for evaluating these point kinetics parameters. Following the strong and the weak form of the multigroup diffusion equation, the weak forms for the partial forward and adjoint solutions, each for a particular active region of the reactor core, are derived. These partial solutions are then used to collapse the weak form into the coupled-reactor point kinetics equations with the approximation of equilibrium. It is noted that these parameters are not meant for transient analysis because of the equilibrium approximation. They are rather used for characterizing the coupled system. The derived reactivity equations are then compared with those in Ref. 1. The equations are implemented in Rattlesnake, the INL MOOSE-based radiation transport application for multi-physics simulations. Numerical results are presented for a simple two-group problem. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-Driven First-Principles Methods for the Study and Design of Alkali Superionic Conductors

We present a detailed exposition of how first-principles methods can be used to guide alkali superionic conductor (ASIC) study and design. Using the argyrodite Li 6 PS 5 Cl as a case study, we demonstrate how modern information technology (IT) infrastructure and software tools can facilitate the assessment of alkali superionic conductors in terms of various critical properties of interest such as phase and electrochemical stability and ionic conductivity. The emphasis is on well-documented, reproducible analysis code that can be readily generalized to other material systems and design problems. For our chosen Li 6 PS 5 Cl case study material, we show that Li excess is crucial to enhancing its conductivity by increasing the occupancy of interstitial sites that promote long-range Li + diffusion between cage-like frameworks. The predicted room-temperature conductivities and activation barriers are in reasonably good agreement with experimental values.

42 ENGINEERING↗