SEARCH · Search NASA

Results for “Heat flux calculation”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1994

The Martian atmospheric planetary boundary layer stability, fluxes, spectra, and similarity

This is the first analysis of the high frequency data from the Viking lander and spectra of wind, in the Martian atmospheric surface layer, along with the diurnal variation of the height of the mixed surface layer, are calculated for the first time for Mars. Heat and momentum fluxes, stability, and z(sub O) are estimated for early spring, from a surface temperature model and from Viking Lander 2 temperatures and winds at 44 deg N, using Monin-Obukhov similarity theory. The afternoon maximum height of the mixed layer for these seasons and conditions is estimated to lie between 3.6 and 9.2 km. Estimations of this height is of primary importance to all models of the boundary layer and Martian General Circulation Models (GCM's). Model spectra for two measuring heights and three surface roughnesses are calculated using the depth of the mixed layer, and the surface layer parameters and flow distortion by the lander is also taken into account. These experiments indicate that z(sub O), probably lies between 1.0 and 3.0 cm, and most likely is closer to 1.0 cm. The spectra are adjusted to simulate aliasing and high frequency rolloff, the latter caused both by the sensor response and the large Kolmogorov length on Mars. Since the spectral models depend on the surface parameters, including the estimated surface temperature, their agreement with the calculated spectra indicates that the surface layer estimates are self consistent. This agreement is especially noteworthy in that the inertial subrange is virtually absent in the Martian atmosphere at this height, due to the large Kolmogorov length scale. These analyses extend the range of applicability of terrestrial results and demonstrate that it is possible to estimate the effects of severe aliasing of wind measurements, to produce a models which agree well with the measured spectra. The results show that similarity theory developed for Earth applies to Mars, and that the spectral models are universal.

Tillman, James E.↗

NASA NTRS1972

A computer program for the calculation of laminar and turbulent boundary layer flows

The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.

Dwyer, H. A.↗

NASA NTRS1977

The effect of sea ice extent on the climatology of the Giss general circulation model

Two model simulations were calculated with the GISS general circulation model corresponding to hypothetical maximum and minimum January sea ice conditions in the north and south polar regions. Results indicate that there were large differences in the Northern Hemisphere circulation between maximum and minimum ice conditions in zonally averaged temperature, vertically averaged eddy sensible heat flux, and mean 500 mb geopotential height. The calculated differences are found to be greater than the inherent variability of the model.

Herman, G. F.↗

NASA NTRS1980

The evolution of active region loop plasma

The adjustment of coronal active-region loops to changes in their heating rate is investigated numerically. The one-dimensional hydrodynamic equations are solved subject to boundary conditions in which heat flux-induced mass exchange between coronal and chromospheric components is allowed. The calculated evolution of physical parameters suggests that (1) mass supplied during chromospheric evaporation is much more effective in moderating coronal temperature excursions than when downward heat flux is dissipated by a static chromosphere, and (2) the method by which the chromosphere responds to changing coronal conditions can significantly influence coronal readjustment time scales. Observations are cited which illustrate the range of possible fluctuations in the heating rates.

Krall, K. R.↗

NASA NTRS1981

Liquid water on Mars - An energy balance climate model for CO2/H2O atmospheres

A simple climatic model is developed for a Mars atmosphere containing CO2 and sufficient liquid water to account for the observed hydrologic surface features by the existence of a CO2/H2O greenhouse effect. A latitude-resolved climate model originally devised for terrestrial climate studies is applied to Martian conditions, with the difference between absorbed solar flux and emitted long-wave flux to space per unit area attributed to the divergence of the meridional heat flux and the poleward heat flux assumed to equal the atmospheric eddy heat flux. The global mean energy balance is calculated as a function of atmospheric pressure to assess the CO2/H2O greenhouse liquid water hypothesis, and some latitude-resolved cases are examined in detail in order to clarify the role of atmospheric transport and temperature-albedo feedback. It is shown that the combined CO2/H2O greenhouse at plausible early surface pressures may account for climates hot enough to support a hydrological cycle and running water at present-day insolation and visible albedo levels.

Hoffert, M. I.↗

NASA NTRS1984

A seasonal model of the Saturnian upper troposphere Comparison with Voyager infrared measurements

In the present extension of Cess and Caldwell's (1979) seasonal climate model to the upper troposphere of Saturn, the ring-modulated latitudinal dependence of the insolation, the ring thermal emission, the oblateness and orbital eccentricity of the planet, and the latitudinal variation of the internal heat flux, are taken into account. While these calculations agree with the temperature-latitude profiles retrieved from Voyager IRIS measurements above 0.2-bar pressure level atmospheric strata, the model fails to predict the retrieved temperature-latitude profiles below the 0.3-bar level. This discrepancy may be due to the existence of clouds at these levels.

Bezard, B.↗

NASA NTRS1999

Cryogenic Fuel Tank Draining Analysis Model

One of the technological challenges in designing advanced hypersonic aircraft and the next generation of spacecraft is developing reusable flight-weight cryogenic fuel tanks. As an aid in the design and analysis of these cryogenic tanks, a computational fluid dynamics (CFD) model has been developed specifically for the analysis of flow in a cryogenic fuel tank. This model employs the full set of Navier-Stokes equations, except that viscous dissipation is neglected in the energy equation. An explicit finite difference technique in two-dimensional generalized coordinates, approximated to second-order accuracy in both space and time is used. The stiffness resulting from the low Mach number is resolved by using artificial compressibility. The model simulates the transient, two-dimensional draining of a fuel tank cross section. To calculate the slosh wave dynamics the interface between the ullage gas and liquid fuel is modeled as a free surface. Then, experimental data for free convection inside a horizontal cylinder are compared with model results. Finally, cryogenic tank draining calculations are performed with three different wall heat fluxes to demonstrate the effect of wall heat flux on the internal tank flow field.

Greer, Donald↗

NASA NTRS1977

Vector sweep

Program is numerical procedure for calculating radiation/geometrical configuration factors and/or space-environmental heat-flux histories. System can accomodate complex systems of diverse geometric shapes residing in unlimited variety of mutual shadow conditions. Basic analytical device is finite-difference employment of unit-sphere method of calculating configuration factors.

Moutrie, C. L.↗

NASA NTRS1992

The theory of magnetohydrodynamic wave generation by localized sources. III - Efficiency of plasma heating by dissipation of far-field waves

The fraction of radiation emitted by Alfven waves is calculated by using two separate methods to determine whether the Alfven flux generated in the photosphere is sufficient to heat the corona. One method employs a set of scaling laws for the fluxes as functions of plasma and source parameters; the second method consist of a procedure for calculating the flux in each waveband from the interaction of vector-harmonic components of an arbitrary applied forcing. Both methods indicate that the Alfven flux accounts roughly for half of the total emission. The need to reexamine estimates of the amount of Alfven flux reaching the corona based on observations of plasma disturbances in the photosphere is emphasized.

Collins, William↗

NASA NTRS1995

Detailed modeling analysis for soot formation and radiation in microgravity gas jet diffusion flames

Radiation heat transfer in combustion systems has been receiving increasing interest. In the case of hydrocarbon fuels, a significant portion of the radiation comes from soot particles, justifying the need for detailed soot formation model and radiation transfer calculations. For laminar gas jet diffusion flames, results from this project (4/1/91 8/22/95) and another NASA study show that flame shape, soot concentration, and radiation heat fluxes are substantially different under microgravity conditions. Our emphasis is on including detailed soot transport models and a detailed solution for radiation heat transfer, and on coupling them with the flame structure calculations. In this paper, we will discuss the following three specific areas: (1) Comparing two existing soot formation models, and identifying possible improvements; (2) A simple yet reasonably accurate approach to calculating total radiative properties and/or fluxes over the spectral range; and (3) Investigating the convergence of iterations between the flame structure solver and the radiation heat transfer solver.

Ku, Jerry C.↗

NASA NTRSDate not supplied

Geochemical Insights Into Volcanic and Lithospheric Evolution of Mars

Introduction: The Martian lithosphere plays a pivotal role in volcanic processes, altering surface features, climate shifts, and the planet's early habitability. Crucially, the rigid lithosphere significantly influences Mars' long-term thermal conditions. Studying various geological eons’ volcanic compositions unveils how the lithosphere has evolved, even though understanding the early Noachian volcanic chemistry remains challenging due to weathering and resurfacing. In this study, newly discovered Noachian volcanic terranes[1-2] along with Hesperian and Amazonian volcanic terranes[3] are considered(Figure 1), to infer the evolution of the Martian lithosphere. Figure 1:Mars topographic map derived from MOLA-HRSC data [4] highlights the study regions including the Noachian volcanic terranes, Arabia Terra (AT), Thaumasia Minor (TM), and the Southwest Margin of Icaria Planum (SWIP)[5].Additionally, the Hesperian volcanic terranes such as Syrtis Major (SM), Hesperia Planum (HP), and, the Amazonian volcanic terranes encompass Elysium Mons (EM), Alba Patera (AP), and Olympus Mons (OM) are marked. Methods and Datasets: We use the most recent elemental mass fraction maps derived from gamma spectroscopy data from the Mars Odyssey 2001mission [6]to ascertain the bulk composition of the Martian landscape. Utilizing the updated geochemical provinces dataset of Mars[7], we consider bulk chemistry variations within our study regions. The maps, at a resolution of 50×50, encompass elements like Al, Ca, Fe, Si, K, Th, H2O, Cl, and S, detected through Gamma-Ray Spectroscopy (GRS) with decimeter scale depth sensitivity[6].Additional major and minor elements like Mg, Na, Ti, P, and Mn, using mass balance methods [8]serve as inputs for melting and crystallization simulations. We employ the pMELTS model with FMQ-3to FMQ oxygen fugacity to estimate the Pressure(P)and melting degree(F)in various eons of volcanic terranes. We have also used conventional thermobarometric calculations to determine their formation P-T conditions, using silica activity and olivine-melt Mg-exchange thermometry[9-10]. Results: Chemical Composition of Studied Volcanic Terranes: To assess formation conditions, we have investigated if GRS-measured compositions signify primary igneous processes by evaluating weathering through the Chemical Index of Alteration (CIA), K/Th ratios, and ternary plots for aqueous alteration. Consistent K/Th ratios and low CIA (<50) across volcanic regions suggest minimal large-scale weathering, affirming well-preserved igneous material extending at decimeter depths. In addition, the calculated bulk compositions align with Martian mantle equilibrium, confirming studied volcanic terranes-Arabia Terra, Thaumasia Minor, and SWIP -as primary or less altered compositions, reinforcing their representativeness at a regional GRS scale on Mars. We have also estimated the bulk compositions of Hesperian and Amazonian volcanic terranes shown in Figure 1. All volcanic terranes represent the primary or minimally altered composition. Temporal Evolution of Martian Volcanic Terranes: Our results show P-T conditions ranging from 1.3-1.6 GPa and 1350-1390°C for Noachian terranes, while Hesperian volcanic terranes vary from 1.6-1.7 GPa and 1370-1395°C, and Amazonian volcanic terranes range between 1.9-2.8 GPa and 1380-1415°C. Our analyses correspond to lithospheric thicknesses/ depth of melting for Noachian which ranges from110-135 km, 120-235 km for Hesperian terranes, and 160-235 km for Amazonian terranes. Partial melting percentages (F) range from 8-12% for Noachian, 10-11% for Hesperian, and 10-12% for Amazonian. We have also calculated mantle potential temperatures (Tp) to show variations across epochs(Figure 2). Heat flux estimation shows variations from 51-65 mW/m2 for Noachian, 38-45 mW/m2for Hesperian, and 27-39 mW/m2for Amazonian terranes, implying temporal changes in lithospheric thickness and heat flow. Comparison of Noachian volcanic terranes with Hesperian and Amazonian reveals variations in lithospheric thickness and heat flux, indicating potential implications for Martian surface conditions, including volcanic activity, climate evolution, and early habitability. Figure 2: P-T diagram displays formation conditions of studied Martian volcanic regions using GRS data, alongside Noachian-age surface basalt samples from Gusev, Gale, and Jezero. Dashed lines represent Martian mantle potential temperatures, calculated by adding the latent heat of fusion to the equilibrium temperature at which melt and solid mantle coexist and subtracting the Martian adiabatic gradient corresponding to the depth of melt. The solid line shows the Martian mantle solidus [11].Discussion& Implication: The study explores Martian volcanic activity across eons. Uniform lithospheric thickness from the Noachian to Hesperian eons suggests persistent weak plumes until Hesperian, influencing heat flow variations. Noachian's high heat flux, likely due to concentrated radioactive elements, shows spatial heterogeneity, indicating early Mars' differentiated mantle. Volcanic outgassing plays avital role in the formation of the Martian climate, likely sustained water on the surface, influencing warm-wet or cold-dry conditions. Hesperian's favorable lithospheric heat flow (45 mW/m2) supports extended water existence, resembling Earth's conditions. Whereas, Amazonian's diminished heat flux (25 mW/m2) corresponds to water loss and arid climate. Hesperian's similar lithospheric thickness to mid-Noachian implies prolonged water retention, influencing climate and potential for astrobiological research, expanding perspectives for future Martian missions. Future work: As volcanism was active throughout its history. The emergence of Mars' extensive volcanic regions, known as large igneous provinces, which are linked to intermittent volcanic activity driven by mantle plumes, poses challenges to understanding the scale and persistence of these phenomena in Martian convection models. Therefore, understanding the regional-scale changes in eruptive processes within Martian volcanic provinces over time remains unclear and crucial for deciphering the evolution of the Martian interior without Earth-like plate tectonics. We will continue this work with petrological and thermoelastic analyses to study the compositional and spatiotemporal evolution of the Elysium Volcanic Province.

A Rani↗

NASA NTRS2021

LOFTID Heat Flux Gauge Calibration: What is Truth?

The Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) is a demonstration of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology, which may enable the delivery of heavy payloads to Mars, Venus, and Titan, as well as return to Earth. Unlike rigid aeroshells that are constrained by the size of the rocket’s shroud, inflatable aeroshells can be deployed to a much larger scale, thus allowing a spacecraft to begin its deceleration earlier and experience less heating. On LOFTID, there will be 4 total heat flux gauges (HFG) with a range of 70 W/cm2 and 1 radiometer with a range of 3 W/cm2, arranged as shown in Fig. 1. Both the radiometer and total HFGs are Schmidt-Boelter gauges purchased from an external vendor. Radiative calibrations were performed in-house at NASA Ames’ Sensors and TPS Advanced Research Laboratories (STAR Labs) before and after environmental testing to investigate how the testing affected the sensors' response. Additional rounds of radiative calibration at STAR Labs were also performed in order to investigate the large uncertainties associated with these tests. For example, a survey of multiple calibration facilities concluded that the uncertainty within a given facility was +/-3% [1]. An additional NIST study that calibrated heat flux gauges at 7 different facilities also found the variation in calibration coefficients to be up to ~3% within a given facility, but up to 15% between facilities, suggesting systematic differences between test setups [2]. Finally, the response of heat flux gauges to radiative versus convective heat flux has shown to differ by up to 20% [3],[4]. Because the heat flux gauges on LOFTID will predominantly experience convective heat flux during flight, a convective calibration study was performed at Boeing's Large-Core Arc Tunnel (LCAT) facility. Radiative Calibration Procedure The calibrations performed at STAR Labs utilize a quartz lamp bank (QLB) that provides a maximum heat flux of 50 W/cm2, which bounds the expected LOFTID flight environment. The calibration involves exposing a water-cooled Gardon gauge (reference) and then the unit-under-test (UUT) to 5 different heat fluxes multiple times for 10 seconds each, and then calculating a linear fit. The test setup is shown in Fig. 2. The total HFGs were calibrated at STAR Labs 3 times, denoted as STAR 1 (before environmental testing), STAR 2 (after protoflight vibration and thermal-vacuum testing), and STAR 3 (no change from previous test). All 8 flight-lot total HFGs showed a decrease in full-scale output from STAR 1 to STAR 2 by between 0.5% and 10. The first portion of this investigation was to determine whether the change could be due to differences in temperature between the two calibration runs. A typical linear fit to the calibration data was performed using Eq. 1 where q’ is the heat flux in W/cm2, c is the calibration coefficient, and mV is the sensor output. To account for temperature, the data were fit to a nonlinear function that included both the sensor output (mV) and the temperature from the thermocouple embedded inside the HFG near the surface (T): q'=mV/(c1* T + c0}. The residuals between the fits and the actual data points were calculated for every point, and proven to be much smaller for the temperature-compensated fits than for the linear fits for all sensors. An example is shown in Fig. 3. When the temperature-compensated fits from STAR 1 were applied to the STAR 2 data, the residuals did not improve, suggesting that the change in sensitivity between these two calibration runs was not due to temperature. A third round of calibration (STAR 3) was conducted to further address the temperature dependence of the total HFGs, and the resulting sensitivities matched closely to STAR 2 (within 2%). Temperature-compensated calibration curves were once again fit to the data. In this case, when the temperature-compensated fits from STAR 3 were applied to STAR 2 data, the residuals between the fits and STAR 2 data were much lower than the residuals due to the linear fits. This suggests that the changes seen between STAR 1 and STAR 2 were likely due to actual changes in the sensors caused by the environmental testing between the two calibrations. A modification of the original calibration process, in which the UUT was exposed to each heat flux for just 3 seconds (instead of 10) to reduce the temperature increase during the test, was additionally performed on several of the HFGs. In general, the sensitivities were 1-1.5% lower than from the 10-second tests, but the temperatures were also significantly lower. When the temperature-compensated fits from the 10-second tests were applied to the 3-second test data, the residuals were greatly improved than when just using the linear fits, further suggesting that the temperature-compensated fits may lead to better accuracy than the linear fits in flight. Convective Calibration The second portion of this study was to create a mapping between the radiative and convective calibration coefficients. The majority of the heating during flight will be convective, so it is important to understand how the HFG response differs under these conditions. However, there are no standardized methods for convective calibration [5]. Because the TPS aerothermal response models were validated at LCAT, the same facility was chosen for convective calibration of two of the total HFGs (Fig. 4). Preliminary results showed that the full-scale output was 3% and 8% higher in convective heat flux as compared to radiative heat flux. However, tunnel variation may have contributed to noise and uncertainty in the measurements, and more testing and analysis remains to be done. Scope of Presentation The presentation will include an overview of the changes seen in HFG calibration before and after environmental testing, differences between radiative and convective calibrations, the modeling work done to aid in understanding the sensor response to varying environments, and recommended future work.

H S Alpert↗

NASA NTRS1984

A theoretical model for the cross spectra between pressure and temperature downstream of a combustor

A theoretical model developed to calculate pressure-temperature cross spectra, pressure spectra, temperature spectra and pressure cross spectra in a ducted combustion system is presented. The model assumes the presence of a fluctuating-volumetric-heat-release-rate disk source and takes into account the spatial distribution of the steady-state volumetric-heat flux. Using the model, pressure, velocity, and temperature perturbation relationships can be obtained. The theoretical results show that, at a given air mass flow rate, the calculated pressure-temperature cross spectra phase angle at the combustor exit depends on the model selected for the steady-state volumetric-heat flux in the combustor. Using measurements of the phase angle, an appropriate source region model was selected. The model calculations are compared with the data. The comparison shows good agreement and indicates that with the use of this model the pressure-temperature cross spectra measurements provide useful information on the physical mechanisms active at the combustion noise source.

Miles, J. H.↗

NASA NTRS1996

An Investigation of the Compatibility of Radiation and Convection Heat Flux Measurements

A method for determining time-resolved absorbed surface heat flux and surface temperature in radiation and convection environments is described. The method is useful for verification of aerodynamic, heat transfer and durability models. A practical heat flux gage fabrication procedure and a simple one-dimensional inverse heat conduction model and calculation procedure are incorporated in this method. The model provides an estimate of the temperature and heat flux gradient in the direction of heat transfer through the gage. This paper discusses several successful time-resolved tests of this method in hostile convective heating and cooling environments.

Liebert, Curt H.↗

NASA NTRSDate not supplied

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

NASA NTRS2016

Simulation of Thermographic Responses of Delaminations in Composites with Quadrupole Method

The application of the quadrupole method for simulating thermal responses of delaminations in carbon fiber reinforced epoxy composites materials is presented. The method solves for the flux at the interface containing the delamination. From the interface flux, the temperature at the surface is calculated. While the results presented are for single sided measurements, with ash heating, expansion of the technique to arbitrary temporal flux heating or through transmission measurements is simple. The quadrupole method is shown to have two distinct advantages relative to finite element or finite difference techniques. First, it is straight forward to incorporate arbitrary shaped delaminations into the simulation. Second, the quadrupole method enables calculation of the thermal response at only the times of interest. This, combined with a significant reduction in the number of degrees of freedom for the same simulation quality, results in a reduction of the computation time by at least an order of magnitude. Therefore, it is a more viable technique for model based inversion of thermographic data. Results for simulations of delaminations in composites are presented and compared to measurements and finite element method results.

Winfree, William P.↗

NASA NTRS1994

Numerical modeling studies of wake vortex transport and evolution within the planetary boundary layer

The proposed research involves four tasks. The first of these is to simulate accurately the turbulent processes in the atmospheric boundary layer. TASS was originally developed to study meso-gamma scale phenomena, such as tornadic storms, microbursts and windshear effects in terminal areas. Simulation of wake vortex evolution, however, will rely on appropriate representation of the physical processes in the surface layer and mixed layer. This involves two parts. First, a specified heat flux boundary condition must be implemented at the surface. Using this boundary condition, simulation results will be compared to experimental data and to other model results for validation. At this point, any necessary changes to the model will be implemented. Next, a surface energy budget parameterization will be added to the model. This will enable calculation of the surface fluxes by accounting for the radiative heat transfer to and from the ground and heat loss to the soil rather than simple specification of the fluxes. The second task involves running TASS with prescribed wake vortices in the initial condition. The vortex models will be supplied by NASA Langley Research Center. Sensitivity tests will be performed on different meteorological environments in the atmospheric boundary layer, which include stable, neutral, and unstable stratifications, calm and severe wind conditions, and dry and wet conditions. Vortex strength may be varied as well. Relevant non-dimensional parameters will include the following: Richardson number or Froude number, Bowen ratio, and height to length scale ratios. The model output will be analyzed and visualized to better understand the transport, decay, and growth rates of the wake vortices. The third task involves running simulations using observed data. MIT Lincoln Labs is currently planning field experiments at the Memphis airport to measure both meteorological conditions and wake vortex characteristics. Once this data becomes available, it can be used to validate the model for vortex behavior under different atmospheric conditions. The fourth task will be to simulate the wake in a more realistic environment covering a wider area. This will involve grid nesting, since high resolution will be required in the wake region but a larger total domain will be used. During the first allocation year, most of the first task will be accomplished.

Lin, Yuh-Lang↗