Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat equation”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The effect of micro-encapsulation on thermal characteristics of metallic phase change materials

Micro-encapsulated phase change material with metallic core can work under high temperature for energy storage purpose, making it an appealing candidate for renewable energy technologies, such as thermal energy storage for concentrating solar power plants. Here, the objective of this work is to study the impact of encapsulation on the thermal performance of the micro-encapsulated phase change materials under various operating conditions. To evaluate the thermal performance of the micro-encapsulated phase change material particles, the governing equation of transient heat conduction with phase change in spheres with composite walls is solved numerically. The impact of encapsulation shell on the convective heat transfer rate, total energy absorption, latent heat ratio, energy density, and duration of phase change is presented and analyzed under different operating conditions. Three different core materials including tin, aluminum and copper with distinct Stephan numbers are investigated. In addition to the encapsulation shell effect, this model can be used to study micro-encapsulated phase change materials with different material compositions, core to shell ratios, and it can be used to analyze the impact of air gaps in between materials. Understanding the thermal performance at particle level is essential, and the results of this study could potentially be used as input for thermal energy storage system level analysis, such as the charging and discharging processes.

42 ENGINEERING↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Validation of implosion modeling through direct-drive shock timing experiments at the National Ignition Facility

Precise modeling of shocks in inertial confinement fusion implosions is critical for obtaining the desired compression in experiments. Shock velocities and post-shock conditions are determined by laser-energy deposition, heat conduction, and equations of state. This paper describes experiments at the National Ignition Facility (NIF) where multiple shocks are launched into a cone-in-shell target made of polystyrene, using laser pulse shapes with two or three pickets and varying on-target intensities. Shocks are diagnosed using the Velocity Interferometric System for Any Reflector diagnostic. Simulated and inferred shock velocities agree well for the range of intensities studied in this work. These results are the first shock-velocity measurements using direct drive on the NIF and provide a good measure of early-time, directly driven laser-energy coupling. Furthermore, the validated models add to the credibility of direct-drive-ignition designs at the megajoule scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Photochemically Induced Acousto-optics Fluid Simulations

PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.

Oudin, AlbertineN [Lawrence Livermore National Lab↗
DOE OSTI2022

The Application of a Genetic Algorithm to the Optimization of a Mesoscale Model for Emergency Response

Besides solving the equations of momentum, heat, and moisture transport on the model grid, mesoscale weather models must account for subgrid-scale processes that affect the resolved model variables. These are simulated with model parameterizations, which often rely on values preset by the user. Such “free” model parameters, along with others set to initialize the model, are often poorly constrained, requiring that a user select each from a range of plausible values. Finding the values to optimize any forecasting tool can be accomplished with a search algorithm, and one such process—the genetic algorithm (GA)—has become especially popular. As applied to modeling, GAs represent a Darwinian process: an ensemble of simulations is run with a different set of parameter values for each member, and the members subsequently judged to be most accurate are selected as “parents” who pass their parameters onto a new generation. At the U.S. Department of Energy’s Savannah River Site in South Carolina, we are applying a GA to the Regional Atmospheric Modeling System (RAMS) mesoscale weather model, which supplies input to a model to simulate the dispersion of an airborne contaminant as part of the site’s emergency response preparations. An ensemble of forecasts is run each day, weather data are used to “score” the individual members of the ensemble, and the parameters from the best members are used for the next day’s forecasts. As meteorological conditions change, the parameters change as well, maintaining a model configuration that is best adapted to atmospheric conditions. Significance Statement We wanted to develop a forecasting system by which a weather model is run over the Savannah River Site each day and repeatedly adjusted according to how well it performed the previous day. To run the model, a series of values (parameters) must be set to control how the model will calculate winds, temperatures, and other desired variables. Each day the model was run several times using different combinations of these parameters and later compared with observed meteorological conditions. Parameters that produced the most accurate forecasts were preferentially reused to create the forecasts for the next day. The process was tested for the summer of 2020 and exhibited lower errors than forecasts produced by the model using default values of the parameters.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1962

Measurement of Local Heat Transfer Coefficients for Flow of Hydrogen and Helium in a Smooth Tube at High Surface to Fluid Bulk Temperature Ratios

Local values of heat transfer coefficients have been measured experimentally for helium and hydrogen gas flowing through an electrically heated inconel tube. The experiment was conducted pri­marily to determine the effect on the heat transfer coefficient of a large density change, radially, in the heat transfer fluid. This large density change was accomplished with relatively high surface temperatures as compared to fluid bulk temperatures or more commonly referred to as high surface to fluid bulk tempera­ture ratio. The large temperature ratio was achieved by precooling the gas with liquid nitrogen. Data were measured for local values of surface to fluid bulk temperature ratios up to 4.5, Reynolds numbers in the turbulent flow region, surface temperatures up to 2300° R, heat flux up to 1,600,000 Btu/(hr)(sq ft) and length to diameter ratio of 250. A comparison of this data with the conventional heat transfer· correlation equation (Dittus-Boelter eq.) is shown on a curve of Nusselt number divided by Prandtl number versus the Reynolds number. The gas properties were evaluated at the film temperature and the Reynolds number was modified by evaluating the velocity term at the fluid bulk temperature and the density at the film temperature.

Weiland, Walter F.↗
NASA NTRS1968

Photoelectron escape from the ionosphere.

Photoelectron escape from ionosphere to magnetosphere, noting flux dependence on energy and factors affecting transport equations and protonosphere heating

Nisbet, J. S.↗
NASA NTRS1974

Data correlation and analysis of arc tunnel and wind tunnel tests of RSI joints and gaps. Volume 1: Technical report

Heat transfer data measured in gaps typical of those under consideration for joints in space shuttle reusable surface insulation protection systems have been assimilated, analyzed and correlated. The data were obtained in four NASA facilities. Several types of gaps were investigated with emphasis on simple butt joints. Gap widths ranged from 0.07 to 0.7 cm and depths ranged from 1 to 6 cm. Laminar, transitional and turbulent boundary layer flows over the gap opening were investigated. Three-dimensional heating variations were observed within gaps in the absence of external flow pressure gradients. Heat transfer correlation equations were obtained for several of the tests. Thermal protection system performance with and without gaps was compared for a representative shuttle entry trajectory.

Christensen, H. E.↗
NASA NTRS1975

Auger spectroscopy of HCl interaction with Al at low pressures

The purpose of the present study was to determine some adsorption characteristics of HCl with partially clean polycrystalline Al by Auger electron spectroscopy. The adsorption was found to exhibit a Freundlich-type behavior pattern in which the number of molecules adsorbed varied with HCl pressure linearly on a log-log plot. Assuming that the adsorption behaves according to the Clausius-Clapeyron equation, the isosteric heat of adsorption increases with coverage at low pressures, appearing to approach a constant value of 4-5 kcal/mole.

Lassiter, W. S.↗
NASA NTRS1984

Comparison of predicted and experimental thermal performance of angular-contact ball bearings

Predicted bearing heat generation and bearing temperature were verified by experimental data for ball bearings over a range of sizes, shaft speeds, and lubricant flow rates. The computer program Shaberth requires, as input, a factor which describes the air-oil mixture in the bearing cavity for calculation of the ball drag contribution to bearing heat generation. An equation for this lubricant percent volume in the bearing cavity was derived and appears to be valid over the range of test conditions including bearing bore sizes from 35 to 167 mm and shaft speeds from 1.0 to 3.0 million DN.

Parker, R. J.↗
NASA NTRS1985

The Martian climate: Energy balance models with CO2/H2O atmospheres

Coupled equations are developed for mass and heat transport in a seasonal Mars model with condensation and sublimation of CO2 at the polar caps. Topics covered include physical considerations of planetary as mass and energy balance; effects of phase changes at the surface on mass and heat flux; atmospheric transport and governing equations; and numerical analysis.

Hoffert, M. I.↗
NASA NTRS1985

Geocoronal structure - The effects of solar radiation pressure and the plasmasphere interaction

The theory of planetary exospheres is extended to incorporate solar radiation pressure in a rigorous manner, and an evaporative geocoronal prototype (classical, motionless exobase) is constructed using Liouville's theorem. Model calculations for density and kinetic temperature at points along the earth-sun axis (solar and antisolar directions) reveal an extensive satellite component, comprising approximately 2/3 of the total hydrogen density near 10 earth radii, and a temperature profile suggestive of an isotropic quasi-Maxwellian velocity distribution for the bound component. A geotail is also evident as an enhancement of the density at local midnight compared to local noon that increases outward (from approximately 25 percent at 10 earth radii to over 60 percent at 20 earth radii). Additional mechanisms acting upon the geocorona alter the basic evaporative case in notable ways. Solar ionization has been included in a simple fashion; the effect is to partially deplete the density without otherwise altering the structure. Interaction with a simple plasmasphere via the Boltzmann equation results in 'heating' the geocorona and enhancing the escape flux at the expense of the density of the bound component, an effect not appreciated in earlier studies; the geotail survives this interaction.

Bishop, J.↗
NASA NTRS1985

Computation of flow regimes in parameter space for the AGCE

This report describes the results of a small study program in support of the design studies for NASA's proposed Atmospheric General Circulation Experiment (AGCE). The proposed experiment will model the atmosphere using a hemispherical layer of a dielectric fluid such as silicone oil, heated at the equator, and with a large radial AC electric field producing a temperature-dependent radial body force similar to radial gravity. The effect of terrestrial gravity on the experiment can be eliminated by doing the experiment in space flight. The author developed a series of three computer models to support these design studies. The first two calculate axisymmetric solutions and their stability to small non-axisymmetric perturbations. The third computes three-dimensional solutions. These codes allow the option of solving problems in a cylindrical geometry as well as a rather generally defined spherical layer.

Roberts, G. O.↗
NASA NTRS1986

Thermal stability analysis of the fine structure of solar prominences

The linear thermal stability of a 2D periodic structure (alternatively hot and cold) in a uniform magnetic field is analyzed. The energy equation includes wave heating (assumed proportional to density), radiative cooling and both conduction parallel and orthogonal to magnetic lines. The equilibrium is perturbed at constant gas pressure. With parallel conduction only, it is found to be unstable when the length scale 1// is greater than 45 Mn. In that case, orthogonal conduction becomes important and stabilizes the structure when the length scale is smaller than 5 km. On the other hand, when the length scale is greater than 5 km, the thermal equilibrium is unstable, and the corresponding time scale is about 10,000 s: this result may be compared to observations showing that the lifetime of the fine structure of solar prominences is about one hour; consequently, our computations suggest that the size of the unresolved threads could be of the order of 10 km only.

Demoulin, Pascal↗
NASA NTRS1987

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗