Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Onset of selective laser flash sintering of AlN

Abstract Flash sintering uses a combination of heating and electric fields to rapidly densify ceramics. Previously, it has been shown that a scanning laser can be used to initiate flash sintering in localized regions on an yttria‐stabilized zirconia (YSZ) sample in a process known as selective laser flash sintering (SLFS). In this work, we show using a combination of measurements of electric current flowing through the sample and observations of necks formed between powder particles that aluminum nitride (AlN) can also undergo SLFS. Scan conditions required to initiate SLFS are characterized over a range of laser powers and laser scan speeds in a dry nitrogen environment. It is shown that initiation of SLFS in AlN is governed by both the local input energy density per scan and heat dissipation and a numerical model is developed to predict temperatures during SLFS. Assuming the minimum temperature along the conductive path determines the onset of SLFS, the minimum temperature and time required is 450–670 K in 2–0.25 s for the pressed AlN pellets used in this study for laser scan speeds of 33–300 m/s, laser powers of 10–30 W, and an applied electric field of 3000 V/cm.

Gouws, Andre↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
NASA NTRS1981

Combustor liner durability analysis

An 18 month combustor liner durability analysis program was conducted to evaluate the use of advanced three dimensional transient heat transfer and nonlinear stress-strain analyses for modeling the cyclic thermomechanical response of a simulated combustor liner specimen. Cyclic life prediction technology for creep/fatigue interaction is evaluated for a variety of state-of-the-art tools for crack initiation and propagation. The sensitivity of the initiation models to a change in the operating conditions is also assessed.

Moreno, V.↗
NASA NTRS1987

Electrical conductivity of carbonaceous chondrites and electric heating of meteorite parent bodies

Electromagnetic heating of rock-forming materials most probably was an important process in the early history of the solar system. Electrical conductivity experiments of representative materials such as carbonaceous chondrites are necessary to obtain data for use in electromagnetic heating models. With the assumption that carbon was present at grain boundaries in the material that comprised the meteorite parent bodies, the electrical heating of such bodies was calculated as a function of body size and solar distance using the T-Tauri model of Sonett and Herbert (1977). The results are discussed.

Duba, AL↗
NASA NTRS1990

Nonisothermal Pluto atmosphere models

The present thermal profile calculation for a Pluto atmosphere model characterized by a high number fraction of CH4 molecules encompasses atmospheric heating by solar UV flux absorption and conductive transport cooling to the surface of Pluto. The stellar occultation curve predicted for an atmosphere of several-microbar surface pressures (which entail the existence of a substantial temperature gradient close to the surface) agrees with observations and implies that the normal and tangential optical depth of the atmosphere is almost negligible. The minimum period for atmospheric methane depletion is calculated to be 30 years.

Hubbard, W. B.↗
NASA NTRS1995

Cooling Panel Optimization for the Active Cooling System of a Hypersonic Aircraft

Optimization of cooling panels for an active cooling system of a hypersonic aircraft is explored. The flow passages are of rectangular cross section with one wall heated. An analytical fin-type model for incompressible flow in smooth-wall rectangular ducts with coupled wall conduction is proposed. Based on this model, the a flow rate of coolant to each design minimum mass flow rate or coolant for a single cooling panel is obtained by satisfying hydrodynamic, thermal, and Mach number constraints. Also, the sensitivity of the optimal mass flow rate of coolant to each design variable is investigated. In addition, numerical solutions for constant property flow in rectangular ducts, with one side rib-roughened and coupled wall conduction, are obtained using a k-epsilon and wall function turbulence model, these results are compared with predictions of the analytical model.

Youn, B.↗
NASA NTRS2003

Flamelet Formation In Hele-Shaw Flow

A Hele-Shaw flow apparatus constructed at Michigan State University (MSU) produces conditions that reduce influences of buoyancy-driven flows. In addition, in the MSU Hele-Shaw apparatus it is possible to adjust the heat losses from the fuel sample (0.001 in. thick cellulose) and the flow speed of the approaching oxidizer flow (air) so that the "flamelet regime of flame spread" is entered. In this regime various features of the flame-to-smolder (and vice versa) transition can be studied. For the relatively wide (approx. 17.5 cm) and long (approx. 20 cm) samples used, approximately ten flamelets existed at all times. The flamelet behavior was studied mechanistically and statistically. A heat transfer analysis of the dominant heat transfer mechanisms was conducted. Results indicate that radiation and conduction processes are important, and that a simple 1-D model using the Broido-Shafizadeh model for cellulose decomposition chemistry can describe aspects of the flamelet spread process. Introduction

Wichman, I. S.↗
NASA NTRS2011

Frost Growth and Densification in Laminar Flow Over Flat Surfaces

One-dimensional frost growth and densification in laminar flow over flat surfaces has been theoretically investigated. Improved representations of frost density and effective thermal conductivity applicable to a wide range of frost circumstances have been incorporated. The validity of the proposed model considering heat and mass diffusion in the frost layer is tested by a comparison of the predictions with data from various investigators for frost parameters including frost thickness, frost surface temperature, frost density and heat flux. The test conditions cover a range of wall temperature, air humidity ratio, air velocity, and air temperature, and the effect of these variables on the frost parameters has been exemplified. Satisfactory agreement is achieved between the model predictions and the various test data considered. The prevailing uncertainties concerning the role air velocity and air temperature on frost development have been elucidated. It is concluded that that for flat surfaces increases in air velocity have no appreciable effect on frost thickness but contribute to significant frost densification, while increase in air temperatures results in a slight increase the frost thickness and appreciable frost densification.

Kandula, Max↗
NASA NTRS2011

A Rocket-Base Study of Auroral Electrodynamics Within the Current Closure Ionosphere

The Auroral Current and Electrodynamics Structure (ACES) mission consisted of two sounding rockets launched nearly simultaneously from Poker Flat Research Range, AK on January 29, 2009 into a dynamic multiple-arc aurora. The ACES rocket mission, in conjunction with the PFISR Radar, was designed to observe the three-dimensional current system of a stable auroral arc system. ACES utilized two well instrumented payloads flown along very similar magnetic field footprints, at various altitudes with small temporal separation between both payloads. ACES High, the higher altitude payload (apogee 360 km), took in-situ measurements of the plasma parameters above the current closure region to provide the input signature into the lower ionosphere. ACES Low, the low-altitude payload (apogee 130 km), took similar observations within the current closure region, where cross-field currents can flow. We present results comparing observations of the electric fields, magnetic fields, electron flux, and the electron temperature at similar magnetic footpoints between both payloads. We further present data from all-sky imagers and PFISR detailing the evolution of the auroral event as the payloads traversed regions connected by similar magnetic footpoints. Current measurements derived from the magnetometers on both payloads are further compared. We examine data from both PFISR and observations on the high-altitude payload which we interpreted as a signature of electron acceleration by means of Alfv n waves. We further examine all measurements to understand ionospheric conductivity and how energy is being deposited into the ionosphere through Joule heating. Data from ACES is compared against models of Joule heating to make inferences regarding the effect of collisions at various altitudes.

Kaeppler, Stephen R.↗
NASA NTRSDate not supplied

Preliminary Measurements on the BOLT Geometry in the Supersonic Low Disturbance Tunnel

Experiments were performed in the Mach 3.5 Supersonic Low Disturbance Tunnel on the BOundary-Layer Transition (BOLT) geometry. The goal of this campaign was to assess changes to the transition front by varying the freestream noise and the model surface quality. The model was printed of polycarbonate and is 30% scale of the flight geometry. It was tested under noisy and quiet conditions at different streamwise positions in the tunnel and was also tested before and after improvements were made to the model surface through sanding and gap reduction. Pitch and yaw angles were nominally zero. The model surface temperature was measured using infrared thermography and thermocouples as the unit Reynolds number was swept from 4.75 – 15.8 × 10 6 m −1 . Boundary-layer transition was observed near the midspan of the model for unit Reynolds numbers above 10×10 6 m −1 in quiet flow, while transition was already observed on the shoulders below 5×10 6 m −1 in noisy flow. The model surface enhancements improved the symmetry of the transition front for quiet flow but made only slight differences for noisy flow. In quiet flow, the impact of streamwise positioning on transition was more significant on the model shoulders compared to the central region. The difference in the shoulder heating was likely due to variation in noise radiated from the nozzle sidewalls. The total temperature and initial wall temperatures were near 300 K and 293 K, respectively, which generally resulted in negative convective heat-flux values into the model. Hot-wire anemometry was also utilized to characterize the freestream and to conduct a planar survey near the base of the model for the 7.92 × 10 6 m −1 quiet condition. The survey captured the thickening of the boundary layer at the centerline and the vortical nature of the flow outboard of it. Spectral analysis of the mass-flux fluctuations near the midspan demonstrates that the boundary layer is laminar and suggests that the thick boundary layer caused the relatively warm centerline observed in the infrared images.

boundary layer transition↗
NASA NTRSDate not supplied

Preliminary Measurements on the BOLT Geometry in the Supersonic Low Disturbance Tunnel

Experiments were performed in the Mach 3.5 Supersonic Low Disturbance Tunnel on the BOundary-Layer Transition (BOLT) geometry. The goal of this campaign was to assess changes to the transition front by varying the freestream noise and the model surface quality. The model was printed of polycarbonate and is 30% scale of the flight geometry. It was tested under noisy and quiet conditions at different streamwise positions in the tunnel and was also tested before and after improvements were made to the model surface through sanding and gap reduction. Pitch and yaw angles were nominally zero. The model surface temperature was measured using infrared thermography and thermocouples as the unit Reynolds number was swept from 4.75 -- 15.8 x 10^6 m^{-1}. Boundary-layer transition was observed near the midspan of the model for unit Reynolds numbers above 10 x 10^6 m^{-1} in quiet flow, while transition was already observed on the shoulders below 5 x 10^6 m^{-1} in noisy flow. The model surface enhancements improved the symmetry of the transition front for quiet flow but made only slight differences for noisy flow. In quiet flow, the impact of streamwise positioning on transition was more significant on the model shoulders compared to the central region. The difference in the shoulder heating was likely due to variation in noise radiated from the nozzle sidewalls. The total temperature and initial wall temperatures were near 300 K and 293 K, respectively, which generally resulted in negative convective heat-flux values into the model. Hot-wire anemometry was also utilized to characterize the freestream and to conduct a planar survey near the base of the model for the 7.92 x 10^6 m^{-1} quiet condition. The survey captured the thickening of the boundary layer at the centerline and the vortical nature of the flow outboard of it. Spectral analysis of the mass-flux fluctuations near the midspan demonstrates that the boundary layer is laminar and suggests that the thick boundary layer caused the relatively warm centerline observed in the infrared images.

boundary layer transition↗
NASA NTRS1993

The effect of SST and soil moisture anomalies on GLA model simulations of the 1988 U.S. summer drought

A series of simulations of the late spring and early summer of 1988 were conducted in order to study the relative importance of different boundary forcings to the Goddard Laboratory for Atmospheres (GLA) model's simulation of the heat wave and drought over the Great Plains of the United States during this time period. Separate 60-day simulations were generated from 10, 20, and 30 May 1988 with a variety of boundary condition datasets. For the control experiment, climatological boundary conditions were used. This was followed by experiments in which either the observed 1988 sea surface temperatures (SST) or derived 1988 soil moisture values, or both, were used in place of the climatological fields. Additional experiments were conducted in which only tropical or midlatitude SST anomalies were used. The impact of the different boundary forcings was evaluated relative to the control simulations of the precipitation and surface air temperature over the Great Plains. It was found that the tropical SST anomalies had a significant effect in reducing precipitation in this area, while the midlatitude anomalies did not. Due to the prescribed climatological soil moistures for the SST experiments, a significant increase in surface temperature did not occur in these simulations. In contrast, the simulations with the anomalous 1988 soil moistures produced both a larger reduction of precipitation and a significant increase in surface temperature over the Great Plains. The simulations with both anomalous SST and soil moisture showed only a slight augmentation of the heat wave and drought relative to the experiments with anomalous soil moisture alone.

Atlas, R.↗
DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Parametric Analysis of Cyclic Phase Change and Energy Storage in Solar Heat Receivers

A parametric study on cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, has been performed. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results in the two-phase regime. Results indicate that parametric changes in receiver gas inlet temperature and receiver heat input effects higher sensitivity to changes in receiver gas exit temperatures.

Hall, Carsie A., III↗
DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗
NASA NTRS2017

Experimental Study of Convective Heating on the Back Face and Payload of a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Aeroshell

A wind tunnel test program has been conducted to define convective heating environments on the back-face of a Hypersonic Inflatable Aerodynamic Decelerator aeroshell. Wind tunnel testing was conducted at Mach 6 and Mach 10 at unit Reynolds numbers from 0.5×10(exp 6)/ft to 3.9×10(exp 6)/ft on a 6.3088 in diameter aeroshell model. Global heating data were obtained through phosphor thermography on the aeroshell back face, as well as on the payload and the aeroshell front face. For all test conditions, laminar flow was produced on the aeroshell front face, while the separated wake shear layer and aeroshell back-face boundary layer were transitional or turbulent. Along the leeward centerline of the aeroshell back face and payload centerbody, heating levels increased with both free stream Reynolds number and angle of attack. The Reynolds number dependency was due to increasing strength of wake turbulence with Reynolds number. The angle-of-attack dependency was due to movement of the wake-vortex reattachment point on the aeroshell back face. The maximum heating levels on the aeroshell back face and payload were approximately 5% to 6%, respectively, of the aeroshell front-face stagnation point. To allow for extrapolation of the ground test data to flight conditions, the back face and payload heating levels were correlated as a function of aeroshell front-face peak momentum thickness Reynolds numbers.

Hollis, Brian R.↗
DOE OSTI2024

ECAR-7932 Rev 0 Large Eddy Simulation of MARVEL Reactor Core Subchannel to Evaluate Model Uncertainty of Reynolds-Averaged Navier-Stokes Equation Based Computational Fluid Dynamics Analysis

In the previous work (ECAR-7210), the peak cladding temperature of the MARVEL microreactor has been evaluated by steady-state Reynolds-Averaged Navier-Stokes (RANS) based computational fluid dynamics (CFD) simulations. Although numerical uncertainties of RANS-based CFD simulations has been assessed in ECAR-7210, the model uncertainty of RANS turbulence models must be investigated to resolve the issues related to inaccurate prediction of turbulent heat flux and flow pulsation in a tight lattice rod bundle using the steady-state RANS simulations. Consequently, this ECAR conducted a high-fidelity CFD analysis utilizing Large Eddy Simulation (LES) to generate reference data and investigated the model uncertainty of RANS-based CFD simulations.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗