SEARCH · Engineering Papers
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.
Advanced manned spacecraft heat shield thermal strain analysis - equilibrium and stress equations, singular point, and computer program
Navier-Stokes equation in predicting solar wind flow pattern
A method of determining molecular eigen-frequencies based on the function of Einstein expressing the variation of the atomic heat of various elements is proposed. It is shown that the same equation can be utilized to calculate both atomic heat and optically identifiably eigen-frequencies - at least to an order of magnitude - suggesting that in both cases the same oscillating structure is responsible.
Provided in this study are predicted in-depth temperature and pyrolysis gas pressure distributions for carbon phenolic materials that are externally heated with a laser source. Governing equations, numerical techniques and comparisons to measured temperature data are also presented. Surface thermochemical conditions were determined using the Aerotherm Chemical Equilibrium (ACE) program. Surface heating simulation used facility calibrated radiative and convective flux levels. Temperatures and pyrolysis gas pressures are predicted using an upgraded form of the SINDA/CMA program that was developed by NASA during the Solid Propulsion Integrity Program (SPIP). Multispecie mass balance, tracking of condensable vapors, high heat rate kinetics, real gas compressibility and reduced mixture viscosity's have been added to the algorithm. In general, surface and in-depth temperature comparisons are very good. Specie partial pressures calculations show that a saturated water-vapor mixture is the main contributor to peak in-depth total pressure. Further, for most of the cases studied, the water-vapor mixture is driven near the critical point and is believed to significantly increase the local heat capacity of the composite material. This phenomenon if not accounted for in analysis models may lead to an over prediction in temperature response in charring regions of the material.
Transport property difference effect on flow prediction of blunt body stagnation region heat transfer, noting laminar flow and compressible boundary layer equations
Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass heat transfer model with both radially inward and outward flow. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature, rotation number, Reynolds number, and radius-to-passage hydraulic diameter ratio. Results were correlated and compared to previous results from similar stationary and rotating models with smooth walls and with trip strips normal to the flow direction. It was concluded that (1) both Coriolis and buoyancy must be considered in turbine blade cooling designs with trip strips, (2) the effects of rotation are markedly different depending upon the flow direction, and (3) the heat transfer with skewed trip strips is less sensitive to buoyancy than the heat transfer models with either smooth or normal trips. Therefore, skewed trip strips rather than normal trip strips are recommended and geometry-specific tests are required for accurate design.
The National Aeronautics and Space Administration Armstrong Flight Research Center (Edwards, California) Flight Loads Laboratory is a unique national laboratory that supports thermal, mechanical, thermal/mechanical, and structural dynamics research and testing. A Proportional Integral thermal control system was designed and implemented to support thermal tests. A thermal control algorithm supporting a quartz lamp heater was developed based on the Proportional Integral control concept and a linearized heating process. The thermal control equations were derived and expressed in terms of power levels, integral gain, proportional gain, and differences between thermal setpoints and skin temperatures. Besides the derived equations, user's predefined thermal test information generated in the form of thermal maps was used to implement the thermal control system capabilities. Graphite heater closed-loop thermal control and graphite heater open-loop power level were added later to fulfill the demand for higher temperature tests. Verification and validation tests were performed to ensure that the thermal control system requirements were achieved. This thermal control system has successfully supported many milestone thermal and thermal/mechanical tests for almost a decade with temperatures ranging from 50 F to 3000 F and temperature rise rates from -10 F/s to 70 F/s for a variety of test articles having unique thermal profiles and test setups.
Solutions are presented for the equations describing natural convection in a planar, horizontal layer of liquid with a constant, linear temperature gradient along the unbounded top and bottom surfaces. Both buoyancy and capillary forces are considered to drive thermal convection. The results are compared to earlier analyses for bounded liquid layers and then applied to containerless metal melts undergoing solidification in the microgravity environment of near-earth orbit.
An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.
A near-wall four-equation turbulence model is developed for the calculation of high-speed compressible turbulent boundary layers. The four equations used are the k-epsilon equations and the theta(exp 2)-epsilon(sub theta) equations. These equations are used to define the turbulent diffusivities for momentum and heat fluxes, thus allowing the assumption of dynamic similarity between momentum and heat transport to be relaxed. The Favre-averaged equations of motion are solved in conjunction with the four transport equations. Calculations are compared with measurements and with another model's predictions where the assumption of the constant turbulent Prandtl number is invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant temperature wall boundary conditions are considered. Results for the range of low Mach numbers and temperature ratios investigated are essentially the same as those obtained using an identical near-wall k-epsilon model. In general, the numerical predictions are in very good agreement with measurements and there are significant improvements in the predictions of mean flow properties at high Mach numbers.
A near-wall four-equation turbulence model is developed for the calculation of high-speed compressible turbulent boundary layers. The four equations used are the k-epsilon equations and the theta(exp 2)-epsilon (sub theta) equations. These equations are used to define the turbulent diffusivities for momentum and heat fluxes, thus allowing the assumption of dynamic similarity between momentum and heat transport to be relaxed. The Favre-averaged equations of motion are solved in conjunction with the four transport equations. Calculations are compared with measurements and with another model's predictions where the assumption of the constant turbulent Prandtl number is invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant temperature wall boundary conditions are considered. Results for the range of low Mach numbers and temperature ratios investigated are essentially the same as those obtained using an identical near-wall k-epsilon model. In general, there are significant improvements in the predictions of mean flow properties at high Mach numbers.
A fully coupled magneto-thermo-elastokinetic model of laminated composite, finitely electroconductive plates incorporating geometrical nonlinearities and subjected to a combination of magnetic and thermal fields, as well as carrying an electrical current is developed, In this context. the first-order transversely shearable plate theory in conjunction with von-Karman geometrically nonlinear strain concept is adopted. Related to the distribution of electric and magnetic field disturbances within the plate, the assumptions proposed by Ambartsumyan and his collaborators are adopted. Based on the electromagnetic equations (i.e. the ones by Faraday, Ampere, Ohm, Maxwell and Lorentz), the modified Fourier's law of heat conduction and on the elastokinetic field equations, the 3-D coupled problem is reduced to an equivalent 2- D one. The theory developed herein provides a foundation for the investigation, both analytical and numerical, of the interacting effects among the magnetic, thermal and elastic fields in multi-layered thin plates made of anisotropic materials.
In the present article, we derive a new dispersion relation for slow magnetoacoustic waves invoking the effect of thermal conductivity, compressive viscosity, radiation, and an unknown heating term along with the consideration of heating–cooling imbalance from linearized MHD equations. We solve the general dispersion relation to understand the role of compressive viscosity and thermal conductivity in the damping of slow waves in coronal loops with and without heating–cooling imbalance. We have analyzed the wave damping for the range of loop length L=50--500 Mm, temperature T=5--30 MK, and density ρ=10^(−11)--10^(−9) kg/cu.m. It was found that the inclusion of compressive viscosity along with thermal conductivity significantly enhances the damping of the fundamental mode oscillations in shorter (e.g. L=50 Mm) and super-hot (T>10 MK) loops. However, the role of viscosity in the damping is insignificant in longer (e.g. L=500 Mm) and hot loops (T≤10 MK) where, instead, thermal conductivity along with the presence of heating–cooling imbalance plays a dominant role. For shorter loops at a super-hot regime of temperature, the increment in the loop density substantially enhances the damping of the fundamental modes due to thermal conductivity when viscosity is absent, however, when the compressive viscosity is added the increase in density substantially weakens the damping. Thermal conductivity alone is found to play a dominant role in longer loops at lower temperatures (T≤10 MK), while compressive viscosity dominates the damping at super-hot temperatures (T>10 MK) in shorter loops. The predicted scaling law between damping time (τ) and wave period (P) is found to better match the observed SUMER (Solar Ultraviolet Measurements of Emitted Radiation) oscillations when the heating–cooling imbalance is taken into account in addition to thermal conductivity and compressive viscosity for the damping of the fundamental slow mode oscillations.
In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.
The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.
A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.
Membrane-based dehumidification is currently being considered as a promising solution for the building application due to its low cost and very limited energy consumption. Developing a simple and efficient open-source code simulation tool is important for boosting the optimization and evaluation of such device in HVAC community. This paper reports a first-order physics based model which accounts for the fundamental heat and mass transfer of humid-air vapor at feed side to flow stream at permeate side. The current model comprises two membrane mass transfer submodels (i.e. microstructure model and performance map model); and it adopts a segment-by-segment methodology for discretizing heat and mass transfer governing equations. The model is capable of simulating both dehumidifiers and energy recovery ventilators with parallel-flow cross-flow, and counter-flow configurations. The model was validated with the measurements at appropriate device. The practices in dehumidification and energy recovery exchangers are also discussed. The model and open-source codes are expected to become a solid fundament for developing a more comprehensive and accurate membrane-based dehumidification in the future.
In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.