Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal & statistical models”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗
NASA NTRS2022

New Constraints on the Volatile Deposit in Mercury’s North Polar Crater, Prokofiev

We present new high-resolution topographic, illumination, and thermal models of Mercury’s 112 km-diameter north polar crater, Prokofiev. The new models confirm previous results that water ice is stable at the surface within the permanently shadowed region (PSR) of Prokofiev for geologic timescales. The largest radar-bright region in Prokofiev is confirmed to extend up to several kilometers past the boundary of its PSR making it unique on Mercury for hosting a significant radar-bright area outside a PSR. The near-infrared normal albedo distribution of Prokofiev’s PSR suggests the presence of a darkening agent rather than pure surface ice. Linear mixture models predict at least roughly half of the surface area to be covered with this dark material. Using improved altimetry in this crater, we place an upper limit of 26 m on its ice deposit thickness. The 1 km-baseline topographic slope and roughness of the radar-bright deposit are lower than the non-radar-bright floor although the difference is not statistically significant when compared to the non-radar-bright floor’s natural topographic variations. These results place new constraints on the nature of Prokofiev’s volatile deposit that will inform future missions, such as BepiColombo.

Michael K Barker↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
NASA NTRS2007

An Investigation of the Characterization of Cloud Contamination in Hyperspectral Radiances

In regions lacking direct observations, the assimilation of radiances from infrared and microwave sounders is the primary method for characterizing the atmosphere in the analysis process. In recent years, technological advances have led to the launching of more advanced sounders, particularly in the thermal infrared spectrum. With the advent of these hyperspectral sounders, the amount of data available for the analysis process has and will continue to be dramatically increased. However, the utilization of infrared radiances in variational assimilation can be problematic in the presence of clouds; specifically the assessment of the presence of clouds in an instantaneous field of view (IFOV) and the contamination in the individual channels within the IFOV. Various techniques have been developed to determine if a channel is contaminated by clouds. The work presented in this paper and subsequent presentation will investigate traditional techniques and compare them to a new technique, the C02 sorting technique, which utilizes the high spectral resolution of the Atmospheric Infrared Sounder (AIRS) within the framework of the Gridpoint Statistical Interpolation (GSI) 3DVAR system. Ultimately, this work is done in preparation for the assessment of short-term forecast impacts with the regional assimilation of AIRS radiances within the analysis fields of the Weather Research and Forecast Nonhydrostatic Mesoscale Model (WRF-NMM) at the NASA Short-term Prediction Research and Transition (SPORT) Center.

McCarty, William↗
NASA NTRS1986

Thermospheric nitric oxide and its role in thermospheric dynamics and composition

The global distribution of thermospheric nitric oxide has been measured under different geophysical conditions by the Atmosphere Explorer (AE) and Solar Mesosphere Explorer (SME) satellites. These observations show a great variability in the NO concentration but the existence of a latitudinal gradient is clearly evidenced by statistical maps of the NO meridional distribution. A two dimensional zonally averaged chemical-dynamical model was used to investigate the importance of nitric oxide 5.3 micrometer cooling and its role on thermospheric temperature, dynamics and major gas composition. For this purpose, a 2-D background atmosphere code and an odd nitrogen code were coupled and run to steady state. The NO distribution obtained is in good agreement with the AE-D global picture calculated at solstice for solar minimum activity. The importance of each term in the thermodynamic equation was studied spatially. It is found that the NO I.R. cooling term competes with conduction in the upper thermosphere and reaches its maximum value near 200 km at high summer latitudes. The primary effect of including the NO cooling term is to increase the temperature in the upper winter thermosphere and decrease it in the other hemisphere, if the global average temperature is fixed. In general, the strength of the circulation is decreased and the thermal gradient is smoother. As a consequence of these changes in the wind and temperature fields, the O2, N2 and O densities also react to the effect of the NO cooling. The importance of this process depends on the level of solar activity which controls the NO distribution.

Gerard, J. C.↗
NASA NTRS1992

One Year of Rayleigh Lidar Measurements at Toronto

The Rayleigh lidar technique is being used at York University to study middle atmospheric thermal structure with particular emphasis on the properties and influence of internal gravity waves. These measurements are now carried out on a routine basis (every clear night) in order to determine seasonal variability and, eventually, long term trends. The transmitter consists of a doubled Nd:YAG laser with an output of 0.4 Joules per pulse at 20Hz (8 Watts) which is expanded to obtain a divergence of 0.4 milliradians. Backscattered photons are collected by a half meter cassegrain telescope, filtered to a bandwidth of 1nm, detected by a PMT, and counted with a multi-channel scaling averager. Strong backscatter from lower altitudes is blocked using a rotating chopper (25000rpm) which triggers the laser. The range gate and hardware averaging interval are variable but typically 300m and 5 minutes. Relative density profiles are obtained by simply normalizing to a model (eg. CIRA) at some altitude. Absolute temperature is derived using the hydrostatic equation, ideal gas law, and a model pressure value at the altitude where the photocount statistical uncertainty is about 15 percent - typically 90km for the entire night's average. The lower limit is currently about 30-35km due the presence aerosol from the volcanic eruption of Mt. Pinatubo.

Whiteway, James A.↗
NASA NTRS1994

The chromospheric structure of the cool giant star g Herculis

Non-Local Thermodynamic Equilibrium (LTE) calculations of semiempirical chromospheric models are presented for 30 g Her (M6 III). This star is one of the coolest (T(sub eff) = 3250 K) SRb (semiregular) variable stars and has a mass perhaps as great as 4 solar mass. Chromospheric features we have observed in its spectrum include Mg II h and k; C II) UV0.01, which is sensitive to electron density; Mg I lambda 2852; Ca II H, K, and IRT; Ca I lambda 4227 and lambda 6573; Al II) UV 1; and H alpha. We pay special attention to fitting the C II intersystem lines and the Mg II resonance lines but use all the other features as constraints to some extent. The equations of radiative transfer and statistical equilibrium are solved self-consistently for H I, H(-), H2, He I, C I, C II, Na I, Mg I, Mg II, Al I, Al II, Ca I, and Ca II with the equivalent two-level technique. To simplify these calculations, a one-dimensional hydrostatic, plane-parallel atmosphere is assumed. We investigate 10 separate 'classical' chromospheric models, differing most importantly in total mass column density above the temperature minimum. Synthetic spectra from these models fit some but not all of the observations. These comparisons are discussed in detail. However, we find that no single-component classical model in hydrostatic equilibrium is able to reproduce both the Mg II line profiles and the relative strengths of the CII) lines. In all these models, chromospheric emission features are formed relatively close to the star (approximately less than 0.05 R(sub *). The circumstellar environment has a thick, cool component overlying the Mg II emission region, which is relatively static and very turbulent. Finally, we find that thermalization in the Mg II h and k lines in the coolest giant stars is controlled by continuum absorption from Ca I 4p 4p3 P0 bound-free opacity and not collisional de-excitation as is the case for warmer K giants.

Luttermoser, Donald G.↗
NASA NTRSDate not supplied

Parallelized Carbon Oxidation and Surface Recession Model in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that often consist of a carbon-based reinforcement and a polymeric matrix. During Earth re-entry, they endure high-temperature oxidation and surface recession. Oxidation is an important mechanism for ablation, sometimes leading to the weakening, spallation, or failure of the oxidized fibers at the surface and in the char layer. However, more details are required including accurate material properties of the fiber microstructure, whether the fibers recede homogenously or localized at pits, and the role of pyrolysis outgassing in order to evaluate the role oxidation plays in the degradation and failure mechanisms of these materials. In this work, we demonstrate a parallelizable oxidation-driven ablation model developed for detailed, large-scale simulations in the DSMC code SPARTA. We also develop robust algorithms which handle the conservation of the surface state after an ablation step. Finally, we verify our model for both simple and more complex chemistry as well as microstructures with reaction statistics, oxidation depth calculations, and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗
DOE OSTI2022

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

The RHESSI Microflare Height Distribution

We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6-10 keV count-rate when RHESSI's full sensitivity was available in order to find the smallest events (Christe et al., 2008). Between March 2002 and March 2007, a total of 25,006 events were found. Source locations were determined in the 4-10 keV, 10-15 keV, and 15-30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g (h) oc exp(-h/lambda) where lambda = 6.1 plus or minus 0.3 Mm is the scale height. A power law height distribution with a negative power law index, gamma = 3.1 plus or minus 0.3 is also consistent with the data. Interpreted as thermal loop top sources, these heights are compared to loops generated by a potential field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential field loop height distribution which may be a signature of the flare energization process.

Christe, P.↗
DOE OSTI2020

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Using optical-microwave synergy for estimating surface energy fluxes over semi-arid rangeland

First results of the Walnut Gulch '92 (Arizona, U.S.) experiment concerning the combined use of radar backscattering (ERS-1) and thermal infrared (LANDSAT TM (Thematic Mapper)) data to estimate surface sensible heat flux are reported. The first step investigates the potential use of ERS-1 SAR (Synthetic Aperture Radar) images for surface soil moisture monitoring of the watershed using five calibrated images acquired during the year 1992 (dry to wet conditions). Results show that despite the typical low level of biomass of semi arid rangeland, an attenuation of the soil backscatter (up to 2 dB) can occur during the rainy season mainly due to the vegetation characteristics. A statistical relationship is then used to retrieve the volumetric surface soil moisture from ERS-1 backscattering (sensitivity of 0.23 dB/% moisture) with a resulting root mean square error of 1.3% of soil moisture. In a second step a semi empirical approach based on energy balance relates soil temperatures to this estimated surface soil moisture. Vegetation temperature is then deduced from soil temperatures and LANDSAT TM composite temperature in order to estimate sensible heat flux according to a two layer type model providing an RMSE of 29 W/sq m.

Troufleau, D.↗
NASA NTRS1996

Dynamics of the Impacts

The SL9 impacts are best known by their plumes, several of which were imaged towering 3000 km above Jupiter's limb. The heat released when they fen back on the atmosphere produced the famous infrared main events. The reentry shock was observed directly in CO emission and indirectly through thermally glowing dust. The relatively low color temperature of the latter implies that most of the energy released on reentry was radiated. If so, the R impact released some 3 x 10(exp 26) ergs, and the larger L, K, and G impacts as much as 3 x 10(exp 27) ergs. These estimates agree well with impact energies deduced from tidal disruption computations. The duration of the infrared events measures time aloft and hence provides a second measure of plume height. These indicate that plume height was only weakly dependent on impact energy. Evidently all the plumes were launched at roughly 10-13 km/s. Using a semi-analytic model for the deceleration, disintegration, and destruction of intruding bodies by an ever-vigilant atmosphere, we find that similar plume heights is a direct consequence of smaller impactors exploding at higher altitudes, in such a way that the different explosions were geometrically similar. We then compare the predictions of our model to the Venerian cratering record, which provides a good statistical test of impact deceleration by a thick atmosphere. Chemistry should have provided an independent measure of explosion altitude: abundant shock-generated CS, CS2 and HCN indicates a source in dry jovian air, above the putative water clouds. However, the Galileo Probe results seem to imply that we should expect no more. Observed water and S2 are consistent with a somewhat oxidized gas (presumably the comet itself), but the absence of SO2 and CO2 shows that conditions were neither too oxidizing nor the shocks too hot.

Zahnle, Kevin↗
DOE OSTI2021

Melt extruded versus extrusion compression molded glass-polypropylene long fiber thermoplastic composites

Long fiber thermoplastic (LFT) composites are processed either with extrusion compression or fiber injection processes. The properties of LFT materials are highly influenced by processing techniques, which lead to different porosity contents, fiber length distributions, and fiber orientation distributions. It is important to understand the various LFT processing techniques and their effects on mechanical, thermal, and microscopic properties. This work considered LFT sheets processed via extrusion, which offers highly aligned fibers (referred to as “Tecnogor composites”), and LFT composites produced via extrusion compression molding (ECM), which offers a random distribution of fibers. Tecnogor composites exhibited higher flexural strength (35–65%), flexural modulus (132–172%), tensile strength (39–52%), tensile modulus (67–75%), and Izod impact resistance (195–220%) than the random LFT composites. This response was attributed to the aligned fibers in Tecnogor composites. Mathematical models including Halpin-Tsai and Lavengood-Goettler were used to predict and compare the Young’s modulus of Tecnogor and ECM composites, respectively.

mechanical testing↗
DOE OSTI2025

Fundamental advantages of multijunction thermoradiative cells

Thermoradiative (TR) cells convert heat to work through emission of thermal radiation. Multijunction thermoradiative cells have received little research interest due to the apparent overlap with energy harvesting limits. Through detailed balance formalism, the present study models both single- and multi-diode TR devices, establishing their performance limits and identifying key factors that influence efficiency. For single-junction TR cells, we derive a relationship for the optimal bandgap and show that higher emitter temperatures increase efficiency. When the receiver is at absolute zero temperature, multijunction TR cells have little advantage over the single-junction cell in terms of the maximum output power density. However, for higher receiver temperatures, the multijunction configuration demonstrates significant improvements in both power and efficiency through the optimization of the chemical potential and bandgap of each diode in the complete ensemble. A 500 K emitter and 300 K receiver TR system can achieve a 21.25% increase in efficiency and a 10% improvement in the power density with multijunction architecture compared to a single junction through chemical potential optimization. These findings suggest that multijunction TR cells offer a promising approach to advancing low-grade heat recovery technologies for efficient heat-to-work conversion.

30 DIRECT ENERGY CONVERSION↗
NASA NTRSDate not supplied

Isolating the Surface Type Influence on Arctic Cloud Properties

In response to anthropogenic climate change, substantial declines in sea ice have been observed. An urgent question is whether and how clouds respond to the changing Arctic surface. Due to the important radiative influence of clouds, the response of clouds to sea ice cover change constitutes a potentially important climate feedback within the Arctic. The sign of this feedback is strongly influenced by the seasonality of the cloud response to sea ice loss. Given the importance of this phenomenon, previous research focused on using observations of inter-annual variability to quantify the cloud-sea ice relationship. However, the reliance on inter-annual variability to assess the covariance between clouds and sea ice make it challenging to control for the influence of large-scale meteorology on the clouds. We have devised an approach that applies an in situ, airborne observational strategy to quantify the influence of surface type on Arctic clouds to satellite data. Our event-based method composites cloud property differences for ocean, marginal ice zone, and sea ice surface types for >5000 MIZ crossing events. We find that cloud properties in non-summer months below ~1.5 km are influenced by surface type such that greater cloud fraction and water content are found over ocean relative to sea ice regions. Results are statistically significantly at the 95% confidence level. During summer, the surface type does not influence the cloud properties. The primary cause of the surface type dependent cloud property differences is the thermodynamic profiles differences that occur between the two surface types: namely, ocean footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their sea ice counterparts. We find that the differences in thermodynamic profiles by surface type that correspond to the cloud property differences also explains a significant portion of the variability across events. Our analysis provides evidence that the influence of surface type on cloud properties fundamentally occurs due to the inherent surface temperature differences between the surface types. Significant surface type influences on cloud properties are only found in the presence of surface temperature differences between the surface types and only non-summer months meet this criterion. A conceptual model of this process is as follows, fundamental properties of the ocean and sea ice surface types (e.g., surface albedo, heat capacity/thermal inertia, surface roughness, and thermal conductance) lead to systematic differences in the surface temperature by surface type. In response to this surface temperature difference, atmosphere boundary layer processes (including turbulent mixing and radiation) communicate the surface temperature differences to the lower atmosphere driving the greater stability over sea ice than ocean. These differences in the background thermodynamic state profiles between the surface types yield atmospheric conditions that are more conducive to cloud development over ocean than sea ice. The details of these boundary layer processes are sensitive to the surface type but may not be fundamental to the resulting thermodynamic structure and cloud differences. Thus, our results indicate that to accurately model that cloud response to sea ice loss in climate models the surface type properties that yield the systematic differences in ocean and sea ice temperature (e.g., surface albedo, surface turbulent fluxes parameters, and thermal conductance) must be accurately represented.

Arctic radiation budget↗
DOE OSTI2025

Coherency-Constrained Spectral Clustering for Power Network Reduction

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

42 ENGINEERING↗