Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 379 records · Page 21

Diffusion approximation for large absorption in radiative transfer

A new diffusion model is developed for radiative transfer in particulate media. It includes the effects of higher Legendre moments while avoiding the mathematical complexities of solving multiple coupled differential moment equations and satisfying higher-order boundary conditions. The method accurately extrapolates the conventional Eddington approximation to problems involving large absorption. Although the simplifying assumptions limit the model to nonbeam incidence and spatial homogeneity of the scattering material, they are nonrestrictive with regard to the size and shape of the medium, the character of background reflections, and the type of phase function.

Meador, W. E.↗

Interfacial geometry and D-variation effects in two-phase systems

Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.

Tenney, D. R.↗

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗

Fluorescence Studies of Protein Crystal Nucleation

One of the most powerful and versatile methods for studying molecules in solution is fluorescence. Crystallization typically takes place in a concentrated solution environment, whereas fluorescence typically has an upper concentration limit of approximately 1 x 10(exp -5)M, thus intrinsic fluorescence cannot be employed, but a fluorescent probe must be added to a sub population of the molecules. However the fluorescent species cannot interfere with the self-assembly process. This can be achieved with macromolecules, where fluorescent probes can be covalently attached to a sub population of molecules that are subsequently used to track the system as a whole. We are using fluorescence resonance energy transfer (FRET) to study the initial solution phase self-assembly process of tetragonal lysozyme crystal nucleation, using covalent fluorescent derivatives which crystallize in the characteristic P432121 space group. FRET studies are being carried out between cascade blue (CB-lys, donor, Ex 376 nm, Em 420 nm) and lucifer yellow (LY-lys, acceptor, Ex 425 nm, Em 520 nm) asp101 derivatives. The estimated R0 for this probe pair, the distance where 50% of the donor energy is transferred to the acceptor, is approximately 1.2 nm, compared to 2.2 nm between the side chain carboxyls of adjacent asp101's in the crystalline 43 helix. The short CB-lys lifetime (approximately 5 ns), coupled with the large average distances between the molecules ((sup 3) 50 nm) in solution, ensure that any energy transfer observed is not due to random diffusive interactions. Addition of LY-lys to CB-lys results in the appearance of a second, shorter lifetime (approximately 0.2 ns). Results from these and other ongoing studies will be discussed in conjunction with a model for how tetragonal lysozyme crystals nucleate and grow, and the relevance of that model to microgravity protein crystal growth

Pusey, Marc L.↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

The Atmospheric Effects of HSCT Emissions Simulated by a 3-Wave Interactive Model

An interactive model which couples a semi-spectral dynamical model, a radiative transfer code and a two-dimensional chemistry transport model (2-D CTM), is used to assess the atmospheric effects of the High-Speed Civil Transport (HSCT) engine emissions. The residual mean meridional circulation, the zonal-mean temperature and the eddy diffusion coefficients are calculated using zonal means and three longest zonal waves of dynamical variables integrated in the semi-spectral dynamical model. They are used in the 2-D CTM to simulate the distribution of trace gases in the atmosphere. The simulated ozone is sent to the radiative transfer code to calculate the heating rates, which drive the dynamics. This radiative coupling connects the dynamical and photochemical processes and creates feedback when the atmosphere is perturbed. It is found that in most areas the ozone depletion caused by HSCT emissions calculated using the 3-wave model has the features similar to, but with significantly larger magnitude than that calculated by the AER 2-D CTM with prescribed transport parameters and temperature. The difference is mostly due to the differences in the circulation in the two models. The radiative feedback effects are investigated by comparing the ozone depletion calculated with the baseline dynamics and with the dynamics perturbed by the HSCT emissions. The feedback through changes in the residual mean meridional circulation and the eddy diffusion coefficients has moderate effects on the simulated ozone depletion. It reduced the ozone depletion by 20-30% in northern mid and high-latitudes. However, the feedback through changes in the zonal- mean temperature is negligible.

Shia, R.-L.↗

Analysis of opposed jet hydrogen-air counter flow diffusion flame

A computational simulation of the opposed-jet diffusion flame is performed to study its structure and extinction limits. The present analysis concentrates on the nitrogen-diluted hydrogen-air diffusion flame, which provides the basic information for many vehicle designs such as the aerospace plane for which hydrogen is a candidate as the fuel. The computer program uses the time-marching technique to solve the energy and species equations coupled with the momentum equation solved by the collocation method. The procedure is implemented in two stages. In the first stage, a one-step forward overal chemical reaction is chosen with the gas phase chemical reaction rate determined by comparison with experimental data. In the second stage, a complete chemical reaction mechanism is introduced with detailed thermodynamic and transport property calculations. Comparison between experimental extinction data and theoretical predictions is discussed. The effects of thermal diffusion as well as Lewis number and Prandtl number variations on the diffusion flame are also presented.

Ho, Y. H.↗

Euler solutions for the flow around a hovering helicopter rotor

A method for computing the flow around a hovering helicopter rotor is presented. The solution procedure is divided into two parts. The first part is a finite volume solution of the Euler equations for the near field flow around a rotor blade. The second part is a free wake approach for determining the wake geometry and induced velocities. These two parts are solved in a coupled fashion. The rolled-up vortex wake of the rotor is included in the Euler computational domain through a nonlinear perturbation technique. This eliminates the numerical diffusion of the wake vorticity due to truncation error and artificial viscosity without the need for excessive grid resolution near the vortex core. The method is used to compute both hovering rotor and wing/vortex interaction flows, and the results are compared to experiment.

Roberts, T. W.↗

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗

Phosphate Reactions as Mechanisms of High-Temperature Lubrication

One of the major problems preventing the operation of advanced gas turbine engines at higher temperatures is the inability of currently used liquid lubricants to survive at these higher temperatures under friction and wear conditions. Current state-of-the-art organic liquid lubricants rapidly degrade at temperatures above 300 C; hence some other form of lubrication is necessary. Vapor-phase lubrication is a promising new technology for high-temperature lubrication. This lubrication method employs a liquid phosphate ester that is vaporized and delivered to bearings or gears; the vapor reacts with the metal surfaces, generating a solid lubricious film that has proven very stable at high temperatures. In this study, solid lubricious films were grown on cast-iron foils in order to obtain reaction and diffusion rate data to help characterize the growth mechanism. A phenomenological mathematical model of the film deposition process was derived incorporating transport and kinetic parameters that were coupled to the experimental data. This phenomenological model can now be reliably used as a predictive and scale-up tool for future vapor-phase lubrication studies.

Nagarajan, Anitha↗

Energetic electrons in Jupiter's magnetosphere

A theoretical model for the energetic electron fluxes in the Jovian magnetosphere is developed. Electrons are transported inward from the solar wind or Jovian magnetospheric tail by radial diffusion. The radial diffusion is driven by fluctuating ionospheric dynamo electric fields associated with a neutral-wind tidal eigenmode at ionospheric altitudes. The tidal mode is excited by the electromagnetic coupling of the solar wind to the polar ionosphere. Two injection models are considered: (1) electron penetration through the dayside magnetopause - low-energy model; and (2) injection of electrons from an assumed magnetospheric tail - high-energy model. Both thermal solar-wind electrons and energetic solar-flare electrons are considered.

Coroniti, F. V.↗

The effect of increased carbon dioxide concentrations on stratospheric ozone

Consideration is given to the influence on ozone of an increased carbon dioxide concentration, for which a measurable growth has been observed in the recent past. Increased carbon dioxide can indirectly affect ozone by perturbing atmospheric temperatures, which will alter ozone production, whose rate displays a fairly strong temperature dependence. This paper presents one-dimensional model results for the steady state ozone behavior when the CO2 concentration is twice its ambient level; the results account for coupling between chemistry and temperature. When the CO2 level doubled, the total ozone burden increased in relation to the ambient burden by 1.2-2.5%, depending on the vertical diffusion coefficient used. Above 30 km, ozone concentrations were larger than the ambient values, a maximum increase of 16% being reached at 43 km. In this region the relative variations were insensitive to the choice of diffusion coefficient. Below 30 km, ozone concentrations were smaller than the unperturbed values and were sensitive to the vertical diffusion profile in this region (10-30 km).

Boughner, R. E.↗

A kinetic approach to compressible turbulence

A kinetic theory of strong turbulence is developed by decomposing fluctuation of distribution function into scaled components, or ranks, to represent three transport processes: macroscopic evolution, transport property and relaxation. Closure is obtained by formulating a relaxation process for the approach of eddy diffusivity to equilibrium. The kinetic equation is solved and three transport functions of turbulence in adiabatic gas are derived: two transfer functions describe turbulent mode transfer, and a coupling function describes the build-up of kinetic energy at the expense of internal energy. The energy balance in the inertial subrange of homogeneous and isotropic turbulence in an adiabatic gas yields a k to the -3rd law, in agreement with numerical computations from one dimensional model of compressible turbulence.

Tchen, C. M.↗

Two-Dimensional Array Beam Scanning Via Externally and Mutually Injection Locked Coupled Oscillators

Some years ago, Stephan proposed an approach to one dimensional (linear) phased array beam steering which requires only a single phase shifter. This involves the use of a linear array of voltage-controlled electronic oscillators coupled to nearest neighbors. The oscillators are mutually injection locked by controlling their coupling and tuning appropriately. Stephan's approach consists of deriving two signals from a master oscillator, one signal phase shifted with respect to the other by means of a single phase shifter. These two signals are injected into the end oscillators of the array. The result is a linear phase progression across the oscillator array. Thus, if radiating elements are connected to each oscillator and spaced uniformly along a line, they will radiate a beam at an angle to that line determined by the phase gradient which is, in turn, determined by the phase difference between the injection signals.The beam direction is therefore controlled by adjusting this phase difference. Recently, Pogorzelski and York presented a formulation which facilitates theoretical analysis of the above beam steering technique. This was subsequently applied by Pogorzelski in analysis of two dimensional beam steering using perimeter detuning of a coupled oscillator array. The formulation is based on a continuum model in which the oscillator phases are represented by a continuous function satisfying a partial differential equation of diffusion type. This equation can be solved via the Laplace transform and the resulting solution exhibits the dynamic behavior of the array as the beam is steered. Stephan's beam steering technique can be similarly generalized to two-dimensional arrays in which the beam control signals are applied to the oscillators on the perimeter of the array. In this paper the continuum model for this two-dimensional case is developed and the dynamic solution for the corresponding aperture phase function is obtained. The corresponding behavior of the resulting far-zone radiation pattern is displayed as well.

Pogorzelski, Ronald J.↗

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗

Alluvial Fans on Dunes in Kaiser Crater Suggest Niveo-Aeolian and Denivation Processes on Mars

On Earth, cold region sand dunes often contain inter-bedded sand, snow, and ice. These mixed deposits of wind-driven snow, sand, silt, vegetal debris, or other detritus have been termed Niveo-aeolian deposits. These deposits are often coupled with features that are due to melting or sublimation of snow, called denivation features. Snow and ice may be incorporated into dunes on Mars in three ways. Diffusion of water vapour into pore spaces is the widely accepted mechanism for the accretion of premafrost ice. Additional mechanisms may include the burial by sand of snow that has fallen on the dune surface or the synchronous transportation and deposition of snow, sand and ice. Both of these mechanisms have been reported for polar dunes on Earth. Niveo-aeolian deposits in polar deserts on Earth have unique morphologies and sedimentary structures that are generally not found in warm desert dunes. Recent analysis of MOC-scale data have found evidence for potential niveo-aeolian and denivation deposits in sand dunes on Mars.

Bourke, M. C.↗

Reliable Spectroscopic Constants for CCH-, NH2- and Their Isotopomers from an Accurate Potential Energy Function

Accurate quartic force fields have been determined for the CCH- and NH2- molecular anions using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). Very large one-particle basis sets have been used including diffuse functions and up through g-type functions. Correlation of the nitrogen and carbon core electrons has been included, as well as other "small" effects, such as the diagonal Born-Oppenheimer correction, and basis set extrapolation, and corrections for higher-order correlation effects and scalar relativistic effects. Fundamental vibrational frequencies have been computed using standard second-order perturbation theory as well as variational methods. Comparison with the available experimental data is presented and discussed. The implications of our research for the astronomical observation of molecular anions will be discussed.

Lee, Timothy J.↗

The Electron Heat Fluxes Associated with Electron Precipitation in the Region of Diffuse Aurora

Electron heat flux that comes from the magnetosphere to the upper ionospheric altitudes controls the value of electron temperature in the core plasma, and, as a result, the total electron density content that is required for different kinds of space weather applications. Knowing the thermal electron heat flux at the upper ionospheric boundaries is the Achilles' heel of all ionospheric models. Such a thermal heat flux setting is especially difficult to justify in the region of the diffuse aurora that is connected to a large energy reservoir of electrons with energies of a few kiloelectron volts, the Earth's plasma sheet, where magnetosphere-ionosphere coupling processes are strongly interconnected. We use the simulated heat flux provided by SuperThermal Electron Transport (STET) code to estimate electron temperatures at the upper ionospheric altitudes and compare these results with corresponding observations from the Defense Meteorological Satellite Program satellite during Saint Patrick's Days 2013 and 2015 Geomagnetic Storms.

George V. Khazanov↗