SEARCH · Engineering Papers

Results for “Heat Flux”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1994

An implementation of a chemical and thermal nonequilibrium flow solver on unstructured meshes and application to blunt bodies

This paper presents a nonequilibrium flow solver, implementation of the algorithm on unstructured meshes, and application to hypersonic flow past blunt bodies. Air is modeled as a mixture of five chemical species, namely O2, N2, O, NO, and N, having two temperatures namely translational and vibrational. The solution algorithm is a cell centered, point implicit upwind scheme that employs Roe's flux difference splitting technique. Implementation of this algorithm on unstructured meshes is described. The computer code is applied to solve Mach 15 flow with and without a Type IV shock interference on a cylindrical body of 2.5mm radius representing a cowl lip. Adaptively generated meshes are employed, and the meshes are refined several times until the solution exhibits detailed flow features and surface pressure and heat flux distributions. Effects of a catalytic wall on surface heat flux distribution are studied. For the Mach 15 Type IV shock interference flow, present results showed a peak heat flux of 544 MW/m2 for a fully catalytic wall and 431 MW/m(exp 2) for a noncatalytic wall. Some of the results are compared with available computational data.

Prabhu, Ramadas K.↗

NASA NTRS1988

An attempt to remotely sense from space the surface heat budget over the Indian Ocean during the 1979 monsoon

Satellite data are used to estimate the net surface heat flux (Q) over the Indian Ocean during June 1979. Ten-day and monthly average fields of Q and its components are produced and analyzed in relation to monsoon dynamics. When compared to ship-based estimates, the satellite-derived fluxes exhibit correct orders of magnitude and their temporal evolution is consistent with the present monsoon knowledge. This study, which represents a first attempt to remotely sense from space the earth surface net heat flux, demonstrates that current satellite sensor data can be combined to accurately describe net heat flux changes in areas such as the Arabian Sea, where they are large, rapid, and spatially extended, and therefore not fully observed by ships.

Gautier, Catherine↗

NASA NTRS1992

Postflight aerothermodynamic analysis of Pegasus(tm) using computational fluid dynamic techniques

The objective was to validate the computational capability of the NASA Ames Navier-Stokes code, F3D, for flows at high Mach numbers using comparison flight test data from the Pegasus (tm) air launched, winged space booster. Comparisons were made with temperature and heat fluxes estimated from measurements on the wing surfaces and wing-fuselage fairings. Tests were conducted for solution convergence, sensitivity to grid density, and effects of distributing grid points to provide high density near temperature and heat flux sensors. The measured temperatures were from sensors embedded in the ablating thermal protection system. Surface heat fluxes were from plugs fabricated of highly insulative, nonablating material, and mounted level with the surface of the surrounding ablative material. As a preflight design tool, the F3D code produces accurate predictions of heat transfer and other aerodynamic properties, and it can provide detailed data for assessment of boundary layer separation, shock waves, and vortex formation. As a postflight analysis tool, the code provides a way to clarify and interpret the measured results.

Kuhn, Gary D.↗

NASA NTRSDate not supplied

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗

DOE OSTI2021

Divertor detachment in the pre-fusion power operation phase in ITER during application of resonant magnetic perturbations

Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Ammonium perchlorate gasification and combustion at high heating rates and low pressures.

Mass-spectrometric and linear regression rate characterizations are reported, derived from CO2 laser pyrolyses of pressed NH2ClO4 (AP) at incident heat fluxes ranging from 25 to 4000 cal/sq cm sec. Product evolution-rate histories were obtained in vacuo by time-resolved (5 msec) mass spectrometry during (1) transient heat-up, and (2) subsequent quasi-steady vaporization (QSV). Vaporization induction times were obtained for (1); these, coupled with heat-transfer approximations neglecting thermochemical heat release, indicated that optical absorption at 10.6 microns dominated over conduction for heat fluxes much greater than 300. Conclusions applying in vacuo were: preferential desorption of NH3, with net accumulation of adsorbed HClO4, occurred during transient heat-up and onset of condensed phase decomposition (CPD), but preferential decomposition of adsorbed HClO4 (compared to NH3) occurred during QSV when CPD was significant. CPD was the dominant mode of QSV at moderate heat fluxes.

Pellett, G. L.↗

DOE OSTI2021

ELM and inter-ELM heat and particle flux to a secondary divertor in the DIII-D tokamak

DIII-D ELMing discharges with a secondary X-point (XPT) inside the vessel were used to evaluate the power distribution between the primary and secondary divertors. The magnetic balance, defined by the distance between the primary and secondary separatrices at the outer midplane (dRsep), was varied from -5 mm (lower-single-null, LSN) to +16 mm (upper-single-null, USN). In this work, we find that the secondary divertor receives up to 1/3 of the total heat flux (primary + secondary divertors), and it is dominated by the ELM-driven flux. Both the radially-integrated and the peak heat flux to the secondary divertor decay below ~50% of the maximum as dRsep is varied from -5 to +6mm, but the integrated heat flux decay flattens out and it would take dRsep above ~25mm to reduce the ELM heat flux to the secondary divertor below 10% of the one deposited to a well-defined SN. Both the secondary inner and outer strike points receive heat flux during ELMs. The peak heat flux is comparable in both strike points for dRsep < 10mm, but the inner strike point receives less than 15% of the total flux. Furthermore, this is experimental evidence that the secondary inner divertor region receives significant ELM flux which should be taken into account for the design of future generation tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗

NASA NTRS1995

Measuring Fluxes Of Heat To A Plasma-Arc Anode

Three probes constructed to provide measurements indicative of conductive, convective, and radiative transfer of heat from free-burning plasma arc to water-cooled copper anode used in generating arc. Each probe consists mainly of copper body with two thermocouples embedded at locations 4 mm apart along length. Thermocouples provide measure of rate of conduction of heat along probe and transfers of heat from plasma to sensing surface at tip of probe. Probes identical except sensing surface of one uncoated and other two coated with different materials to make them sensitive to different components of overall flux of heat.

Sankovic, John M.↗

NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

NASA NTRS2011

Geomagnetic Secular Variation Prediction with Thermal Heterogeneous Boundary Conditions

It has long been conjectured that thermal heterogeneity at the core-mantle boundary (CMB) affects the geodynamo substantially. The observed two pairs of steady and strong magnetic flux lobes near the Polar Regions and the low secular variation in the Pacific over the past 400 years (and perhaps longer) are likely the consequences of this CMB thermal heterogeneity. There are several studies on the impact of the thermal heterogeneity with numerical geodynamo simulations. However, direct correlation between the numerical results and the observations is found very difficult, except qualitative comparisons of certain features in the radial component of the magnetic field at the CMB. This makes it difficult to assess accurately the impact of thermal heterogeneity on the geodynamo and the geomagnetic secular variation. We revisit this problem with our MoSST_DAS system in which geomagnetic data are assimilated with our geodynamo model to predict geomagnetic secular variations. In this study, we implement a heterogeneous heat flux across the CMB that is chosen based on the seismic tomography of the lowermost mantle. The amplitude of the heat flux (relative to the mean heat flux across the CMB) varies in the simulation. With these assimilation studies, we will examine the influences of the heterogeneity on the forecast accuracies, e.g. the accuracies as functions of the heterogeneity amplitude. With these, we could be able to assess the model errors to the true core state, and thus the thermal heterogeneity in geodynamo modeling.

Kuang, Weijia↗

NASA NTRS2013

Mars Science Laboratory Heatshield Aerothermodynamics: Design and Reconstruction

The Mars Science Laboratory heatshield was designed to withstand a fully turbulent heat pulse based on test results and computational analysis on a pre-flight design trajectory. Instrumentation on the flight heatshield measured in-depth temperatures in the thermal protection system. The data indicate that boundary layer transition occurred at 5 of 7 thermocouple locations prior to peak heating. Data oscillations at 3 pressure measurement locations may also indicate transition. This paper presents the heatshield temperature and pressure data, possible explanations for the timing of boundary layer transition, and a qualitative comparison of reconstructed and computational heating on the as-flown trajectory. Boundary layer Reynolds numbers that are typically used to predict transition are compared to observed transition at various heatshield locations. A uniform smooth-wall transition Reynolds number does not explain the timing of boundary layer transition observed during flight. A roughness-based Reynolds number supports the possibility of transition due to discrete or distributed roughness elements on the heatshield. However, the distributed roughness height would have needed to be larger than the pre-flight assumption. The instrumentation confirmed the predicted location of maximum turbulent heat flux near the leeside shoulder. The reconstructed heat flux at that location is bounded by smooth-wall turbulent calculations on the reconstructed trajectory, indicating that augmentation due to surface roughness probably did not occur. Turbulent heating on the downstream side of the heatshield nose exceeded smooth-wall computations, indicating that roughness may have augmented heating. The stagnation region also experienced heating that exceeded computational levels, but shock layer radiation does not fully explain the differences.

Edquist, Karl T.↗

DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2012

Space-based Observational Constraints for 1-D Plume Rise Models

We use a space-based plume height climatology derived from observations made by the Multi-angle Imaging SpectroRadiometer (MISR) instrument aboard the NASA Terra satellite to evaluate the ability of a plume-rise model currently embedded in several atmospheric chemical transport models (CTMs) to produce accurate smoke injection heights. We initialize the plume-rise model with assimilated meteorological fields from the NASA Goddard Earth Observing System and estimated fuel moisture content at the location and time of the MISR measurements. Fire properties that drive the plume-rise model are difficult to estimate and we test the model with four estimates for active fire area and four for total heat flux, obtained using empirical data and Moderate Resolution Imaging Spectroradiometer (MODIS) re radiative power (FRP) thermal anomalies available for each MISR plume. We show that the model is not able to reproduce the plume heights observed by MISR over the range of conditions studied (maximum r2 obtained in all configurations is 0.3). The model also fails to determine which plumes are in the free troposphere (according to MISR), key information needed for atmospheric models to simulate properly smoke dispersion. We conclude that embedding a plume-rise model using currently available re constraints in large-scale atmospheric studies remains a difficult proposition. However, we demonstrate the degree to which the fire dynamical heat flux (related to active fire area and sensible heat flux), and atmospheric stability structure influence plume rise, although other factors less well constrained (e.g., entrainment) may also be significant. Using atmospheric stability conditions, MODIS FRP, and MISR plume heights, we offer some constraints on the main physical factors that drive smoke plume rise. We find that smoke plumes reaching high altitudes are characterized by higher FRP and weaker atmospheric stability conditions than those at low altitude, which tend to remain confined below the BL, consistent with earlier results. We propose two simplified parameterizations for computing injection heights for fires in CTMs and discuss current challenges to representing plume injection heights in large scale atmospheric models.

Martin, Maria Val↗

DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2001

Theoretical Prediction of Microgravity Ignition Delay of Polymeric Fuels in Low Velocity Flows

A new flammability apparatus and protocol, FIST (Forced Flow Ignition and Flame Spread Test), is under development. Based on the LIFT (Lateral Ignition and Flame Spread Test) protocol, FIST better reflects the environments expected in spacebased facilities. The final objective of the FIST research is to provide NASA with a test methodology that complements the existing protocol and provides a more comprehensive assessment of material flammability of practical materials for space applications. Theoretical modeling, an extensive normal gravity data bank and a few validation space experiments will support the testing methodology. The objective of the work presented here is to predict the ignition delay and critical heat flux for ignition of solid fuels in microgravity at airflow velocities below those induced in normal gravity. This is achieved through the application of a numerical model previously developed of piloted ignition of solid polymeric materials exposed to an external radiant heat flux. The model predictions will provide quantitative results about ignition of practical materials in the limiting conditions expected in space facilities. Experimental data of surface temperature histories and ignition delay obtained in the KC-135 aircraft are used to determine the critical pyrolysate mass flux for ignition and this value is subsequently used to predict the ignition delay and the critical heat flux for ignition of the material. Surface temperature and piloted ignition delay calculations for Polymethylmethacrylate (PMMA) and a Polypropylene/Fiberglass (PP/GL) composite were conducted under both reduced and normal gravity conditions. It was found that ignition delay times are significantly shorter at velocities below those induced by natural convection.

Fernandez-Pello, A. C.↗

NASA NTRS2017

The Role of Local Heating in the 2015 Indian Heat Wave

India faced a major heat wave during the summer of 2015. Temperature anomalies peaked in the dry period before the onset of the summer monsoon, suggesting that local land-atmosphere feedbacks involving desiccated soils and vegetation might have played a role in driving the heat extreme. Upon examination of in situ data, reanalysis, satellite observations, and land surface models, we find that the heat wave included two distinct peaks: one in late May, and a second in early June. During the first peak we find that clear skies led to a positive net radiation anomaly at the surface, but there is no significant sensible heat flux anomaly within the core of the heat wave affected region. By the time of the second peak, however, soil moisture had dropped to anomalously low levels in the core heat wave region, net surface radiation was anomalously high, and a significant positive sensible heat flux anomaly developed. This led to a substantial local forcing on air temperature that contributed to the intensity of the event. The analysis indicates that the highly agricultural landscape of North and Central India can reinforce heat extremes under dry conditions.

Debjani Ghatak↗

NASA NTRS1986

Turbulence in a convective marine atmospheric boundary layer

The structure and kinetic energy budget of turbulence in the convective marine atmospheric boundary layer as observed by aircraft during a cold air outbreak have been studied using mixed layer scaling. The results are significantly different from those of previous studies under conditions closer to free convection. The normalized turbulent kinetic energy and turbulent transport are about twice those found during the Air Mass Transformation Experiment (AMTEX). This implies that for a given surface heating the present case is dynamically more active. The difference is mainly due to the greater importance of wind shear in the present case. This case is closer to the roll vortex regime, whereas AMTEX observed mesoscale cellular convection which is closer to free convection. Shear generation is found to provide a significant energy source, in addition to buoyancy production, to maintain a larger normalized turbulent kinetic energy and to balance a larger normalized dissipation. The interaction between turbulent pressure and divergence (i.e., pressure scrambling) is also found to transfer energy from the vertical to the horizontal components, and is expected to be stronger in roll vortices than in m esoscale cells. The sensible heat flux is found to fit well with a linear vertical profile in a clear or subcloud planetary boundary layer (PBL), in good agreement with the results of Lenschow et al., (1980). The heat flux ratio between the PBL top and the surface, derived from the linear fitted curve, is approximately -0.14, in good agreement with that derived from the lidar data for the same case. Near the PBL top, the heat flux profiles are consistent with those of Deardoff (1979) and Deardorff et al. (1980).

Chou, S.-H.↗