Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal diffusion”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Silicon nanocluster anion-argon cation recombination via hybrid continuum-molecular dynamics calculations

Prediction of the charge distribution on particles in an aerosol is critical not only in electrical mobility based characterization methods, but also in understanding the role that charging plays in particle growth in gas phase synthesis reactors. The latter is particularly important in non-thermal plasma synthesis reactors, wherein nanoclusters form and grow from vapor phase precursors in a high electron and high ion density environment. In plasmas, because free electrons are much less massive and much more mobile than positive gas ions, nanocluster charge distributions are biased negative. Furthermore, neutral or even positively charged nanoclusters may exist, depending on the rate of nanocluster-gas ion recombination, and nanocluster-nanocluster collisions may greatly contribute to nanocluster growth if not Coulombically suppressed.

42 ENGINEERING↗
NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗
NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗
NASA NTRS2000

General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT

This update reports a state-of-the-art discrete ordinate algorithm for monochromatic unpolarized radiative transfer in non-isothermal, vertically inhomogeneous, but horizontally homogeneous media. The physical processes included are Planckian thermal emission, scattering with arbitrary phase function, absorption, and surface bidirectional reflection. The system may be driven by parallel or isotropic diffuse radiation incident at the top boundary, as well as by internal thermal sources and thermal emission from the boundaries. Radiances, fluxes, and mean intensities are returned at user-specified angles and levels. DISORT has enjoyed considerable popularity in the atmospheric science and other communities since its introduction in 1988. Several new DISORT features are described in this update: intensity correction algorithms designed to compensate for the 8-M forward-peak scaling and obtain accurate intensities even in low orders of approximation; a more general surface bidirectional reflection option; and an exponential-linear approximation of the Planck function allowing more accurate solutions in the presence of large temperature gradients. DISORT has been designed to be an exemplar of good scientific software as well as a program of intrinsic utility. An extraordinary effort has been made to make it numerically well-conditioned, error-resistant, and user-friendly, and to take advantage of robust existing software tools. A thorough test suite is provided to verify the program both against published results, and for consistency where there are no published results. This careful attention to software design has been just as important in DISORT's popularity as its powerful algorithmic content.

Tsay, Si-Chee↗
NASA NTRS2004

Domes on Europa: The Role of Thermally Induced Compositional Diapirism

The surface of Europa is peppered by topographic domes, interpreted as sites of intrusion and extrusion. Diapirism is consistent with dome morphology, but thermal buoyancy alone cannot produce sufficient driving pressures to create the observed dome elevations. Instead, we suggest that diapirs may initiate by thermal convection that induces compositional segregation within Europa's ice shell. This double-diffusive convection scenario allows sufficient buoyancy for icy plumes to create the observed surface topography, if the ice shell has a very small effective elastic thickness (approximately 0.1 to 0.5 km) and contains low-eutectic-point impurities at the percent level. Thermal buoyancy, compositional buoyancy and double-diffusive convection are discussed.

Pappalardo, R. T.↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Surface tension effects on the onset of double-diffusive convection

Experiments have been carried out to determine the critical thermal Rayleigh number for onset of convection in a horizontal layer of density-stratified fluid with a free surface when heated from below. Three different aqueous solutions were used: salt, glycerol, and acetic acid. The rates of change in surface tension with concentration for these three solutions are positive, nearly zero, and negative, respectively. Compared to the rigid-rigid boundaries, the critical thermal Rayleigh number was found to be larger by 11.2 percent for the salt solution and smaller by 10.0 percent for the glycerol solution. With the acetic acid solution, however, the effect of the free surface was found to be negligible.

Chen, C. F.↗
NASA NTRS1985

Current understanding of point defects and diffusion processes in silicon

The effects of oxidation of Si which established that vacancies (V) and Si self interstitials (I) coexist in Si at high temperatures under thermal equilibrium and oxidizing conditions are discussed. Some essential points associated with Au diffusion in Si are then discussed. Analysis of Au diffusion results allowed a determination of the I component and an estimate of the V component of the Si self diffusion coefficient. A discussion of theories on high concentration P diffusion into Si is then presented. Although presently there still is no theory that is completely satisfactory, significant progresses are recently made in treating some essential aspects of this subject.

Tan, T. Y.↗
NASA NTRS1993

Variable Refractive Index Effects on Radiation in Semitransparent Scattering Multilayered Regions

A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior is readily calculated from the analytical expressions that are obtained. By using many sublayers, expressions provide the temperature distribution and heat flow for a diffusing medium with a continually varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature and heat flux results are given to show the effect of variations in refractive index and optical thickness through the multilayer laminate.

Siegel, R.↗
NASA NTRS1981

Boundary-layer analysis for the convection/diffusion transition in dendritic growth

The supercooling dependence of dendritic growth kinetics under the influence of convective heat transport is investigated theoretically and experimentally with emphasis on theoretical prediction of the supercooling level at which the transition from diffusion-controlled to convection-controlled dendritic growth occurs. It is shown that the crossover between diffusive and convective transport depends on the relative thickness of the Stefan length compared with the thermal boundary layer. These lengths become equal at a supercooling which may be calculated from diffusion theory and fluid mechanics. It is also shown that the crossover supercooling varies weakly with the gravitational acceleration, melt viscosity, and the volumetric expansion coefficient.

Glicksman, M. E.↗
NASA NTRS2006

Particle Demagnetization in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS2006

How Does Collisionless Magnetic Reconnection Work in the Presence of a Guide Magnetic Field?

The dissipation mechanism of guide field magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide-field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS1994

Diffused ionized gas in the spiral galaxy M31

We present a study of diffuse ionized gas (DIG, or WIM for warm ionized medium) in the nearby spiral M31, based on charge coupled device (CCD) imagery in H alpha and (S II) emission lines. We find that the DIG in M31 contributes 40% of the observed total H alpha luminosity, and at least 20% of the total after correction for extinction. This fraction is constant across the disk of M31. The total H alpha luminosity of M31, corrected for extinction, is about 4.1 x 10(exp 40) ergs/s. This implies a low star formation rate of about 0.35 M solar mass/yr, in agreement with earlier estimates. DIG is concentrated near regions of star formation. The average emission measure (EM) perpendicular to the disk ranges from about 6 pc cm(exp-6) in weaker spiral arms to 15 pc cm(exp-6) in the main spiral arms, with (rare) peak values reaching 50 pc cm(exp-6). Overall, the EM in the main arms is substantially higher than that of DIG in the solar neighborhood, implying higher power requirements per unit area in these regions of M31. However, the strong concentration of DIG near regions of star formation keeps the total area covered by DIG in M31 relatively small, which makes it possible to sustain this gas in spite of the low star formation rate. We argue that the DIG extent in vertical direction in M31 may well be lower than in our Galaxy, because of its low star formation rate. A simple parameterization of the electron distribution then shows that the DIG in the spiral arms may contribute between 15% and 40% of the total atomic gas column. The average electron densities implied by the high EMs reached in some regions seem to indicate a higher thermal pressure in the DIG than in the diffuse H I gas in M31. This pressure in-equilibrium can be avoided if part of the diffuse emissions due to a denser component, as proposed most recently by Miller & Cox, because this would reduce the EM and implied pressures for the actual widespread DIG. Column densities of the ionized gas might then be lower as well. We also show, however, that it is not possible to attribute all of the diffuse emission to a relatively dense component, such as earlier proposed extended low-density H II regions. The DIG in M31 shows an overall (S II)/H alpha intensity ratio of 0.5, significantly higher than in discrete H II regions, and similar to that in other galaxies. This number does not vary with radius in M31. There is a suggestion that (S II)/H alpha intensity may be higher at the lowest detectable emission measures, up to 0.7 or even 1. The discrete source population in M31 also shows a trend toward increasing (S II)/H alpha flux ratios as the sources become more diffuse. Few sources (apart from supernova remnants), however, reach the value of 0.5 found in the DIG as a whole. The distinct spectral signature of the DIG in M31, and its energy requirements favor photoionization of the gas, although it is likely that shock ionization occurs in selected regions.

Walterbos, Rene A. M.↗
NASA NTRS2006

The Dissipation Mechanism in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗
NASA NTRS2009

What Turns Reconnection On and Off

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density. which matches key features of electron distribution functions.

Hesse, Michael↗
DOE OSTI2023

Diffuse neutrino background from past core collapse supernovae

Core collapse supernovae are among the most powerful explosions in the Universe, which emit thermal neutrinos that carry away most of the gravitational binding energy released. These neutrinos produce a diffuse supernova neutrino background (DSNB), which is one of the largest energy budgets among all radiation backgrounds. Detecting the DSNB is an important goal of modern high-energy astrophysics and particle physics, which provides valuable insights into core collapse modeling, neutrino physics, and cosmic supernova rate history. In this review, the key ingredients of DSNB calculation and what can be learned from future detections, including black hole formation and non-standard neutrino interactions are discussed. Moreover, an overview of the latest updates in neutrino experiments, which could lead to the detection of the DSNB in the next decade, is provided. With the promise of this breakthrough discovery on the horizon, the study of DSNB has great potential to further our understanding of the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗