Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal modeling techniques”

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 109 records · Page 6

Monte Carlo analysis of lobular gas-surface scattering in tubes applied to thermal transpiration

A model of rarefied gas flow in tubes was developed which combines a lobular distribution with diffuse reflection at the wall. The model with Monte Carlo techniques was used to explain previously observed deviations in the free molecular thermal transpiration ratio which suggest molecules can have a greater tube transmission probability in a hot-to-cold direction than in a cold-to-hot direction. The model yields correct magnitudes of transmission probability ratios for helium in Pyrex tubing (1.09 to 1.14), and some effects of wall-temperature distribution, tube surface roughness, tube dimensions, gas temperature, and gas molecular mass.

Smith, J. D.↗

Thermal modeling of nickel-hydrogen battery cells operating under transient orbital conditions

An analytical study of the thermal operating characteristics of nickel-hydrogen battery cells is presented. Combined finite-element and finite-difference techniques are employed to arrive at a computationally efficient composite thermal model representing a series-cell arrangement operating in conjunction with a radiately coupled baseplate and coldplate thermal bus. An aggressive, low-mass design approach indicates that thermal considerations can and should direct the design of the thermal bus arrangement. Special consideration is given to the potential for mixed conductive and convective processes across the hydrogen gap. Results of a compressible flow model are presented and indicate the transfer process is suitably represented by molecular conduction. A high-fidelity thermal model of the cell stack (and related components) indicates the presence of axial and radial temperature gradients. A detailed model of the thermal bus reveals the thermal interaction of individual cells and is imperative for assessing the intercell temperature gradients.

Schrage, Dean S.↗

The atmosphere of Saturn - An infrared perspective

A model of the radiative portion of the equatorial atmosphere of Saturn, constrained by the infrared data of various observers, has been constructed using a technique which includes the variation of thermal flux with depth. The model has a high-altitude temperature inversion due to the absorption of ultraviolet sunlight. The inversion causes the observed infrared emission peaks at 8 microns (methane) and 12 microns (ethane). Mixing ratios of these gases to hydrogen are computed from these emission features. The bottom of the modeled region occurs at the radiative convective boundary. At this level, an opaque cloud consisting of solid ammonia condensation particles is postulated. Above the cloud is a thin haze, also composed of ammonia particles. The haze is required to match infrared observations near 9.5 microns and hydrogen quadrupole equivalent widths near 0.64 micron. Predictions of the model are given for further observational tests.

Caldwell, J.↗

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗

Ion thruster thermal characteristics and performance.

Experimental and analytical studies of mercury-fed hollow cathode ion thrusters are presented. Pertinent thermal properties and plasma heating characteristics of the system have been determined through correlation and integration of temperature measurement data with a single stage Wiener-Kalman estimation technique. The resultant analytical network model is applied to predict the thermal performance of specified thruster systems under various operating conditions. The effects of solar intensity levels on the thruster thermal behavior, the dependency on the thruster dimension and multithruster design arrangements have been examined.

Wen, L.↗

Impact of uranium oxide (UO 2 ) fuel with molybdenum (Mo) inserts on pressurized water reactor performance and safety

This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal transport in defective and disordered materials

With significant recent advancements in thermal sciences—such as the development of new theoretical and experimental techniques, and the discovery of new transport mechanisms—it is helpful to revisit the fundamentals of vibrational heat conduction to formulate an updated and informed physical understanding. The increasing maturity of simulation and modeling methods sparks the desire to leverage these techniques to rapidly improve and develop technology through digital engineering and multi-scale, electro-thermal models. With that vision in mind, this review attempts to build a holistic understanding of thermal transport by focusing on the often unaddressed relationships between subfields, which can be critical for multi-scale modeling approaches. For example, we outline the relationship between mode-specific (computational) and spectral (analytical) models. We relate thermal boundary resistance models based on perturbation approaches and classic transmissivity based models. We discuss the relationship between lattice dynamics and molecular dynamics approaches along with two-channel transport frameworks that have emerged recently and that connect crystal-like and amorphous-like heat conduction. Throughout, we discuss best practices for modeling experimental data and outline how these models can guide material-level and system-level design.

42 ENGINEERING↗

Thermal Conductivity Measurement of Extrusion-Printed Silver Using Modulated Photothermal Radiometry

Flexible printed electronics is a rapidly growing field with applications in conformal and flexible devices. However, the physical properties of the films created by many state-of-the-art printing methods become highly dependent on printing parameters, resulting in varying thermal properties often differing significantly from their bulk ink components. To understand the influence of the printing process, we build upon our previous work, where a noncontact optical technique, known as modulated photothermal radiometry (MPTR), was used to measure the thermal conductivity of aerosol-jet-printed thin films. In this work, we use the method to study the thermal properties of extrusion-printed silver on glass and alumina substrates. A noise-resistant data analysis fitting technique is applied using a 2-D heat transfer model. Here, the thermal conductivity measurement is validated using the Weidemann-Franz (WF) relationship from measured electrical conductivity values.

42 ENGINEERING↗

High-Temperature Strain Sensing for Aerospace Applications

Thermal protection systems (TPS) and hot structures are utilizing advanced materials that operate at temperatures that exceed abilities to measure structural performance. Robust strain sensors that operate accurately and reliably beyond 1800 F are needed but do not exist. These shortcomings hinder the ability to validate analysis and modeling techniques and hinders the ability to optimize structural designs. This presentation examines high-temperature strain sensing for aerospace applications and, more specifically, seeks to provide strain data for validating finite element models and thermal-structural analyses. Efforts have been made to develop sensor attachment techniques for relevant structural materials at the small test specimen level and to perform laboratory tests to characterize sensor and generate corrections to apply to indicated strains. Areas highlighted in this presentation include sensors, sensor attachment techniques, laboratory evaluation/characterization of strain measurement, and sensor use in large-scale structures.

Piazza, Anthony↗

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

3D Characterization of the Columnar-to-Equiaxed Transition in Additively Manufactured Inconel 718

Additive manufacturing (AM) provides enormous processing flexibility, enabling novel part geometries and optimized designs. Access to a local heat source further permits the potential for local microstructure control on the scale of individual melt pools, which can enable local control of part properties. In order to design tailored processing strategies for target microstructures, models predicting the columnar-to-equiaxed transition must be extended to the high solidification velocities and complex thermal histories present in AM. Here, we combine 3D characterization with advanced modeling techniques to develop a more complete understanding of the solidification process and evolution of microstructure during electron beam melting (EBM) of Inconel 718. Full calibration of existing microstructure prediction models demonstrates the differences between AM processes and more conventional welding techniques, underlying the need for accurate determination of key parameters that can only be measured directly in 3D. The ability to combine multisensor data in a consistent 3D framework via data fusion algorithms is essential to fully leverage these advanced characterization approaches. Thermal modeling provides insight on microstructure development within isolated solidification events and demonstrates the role of Marangoni effects on controlling solidification behavior.

Polonsky, Andrew↗

Relationships between the energetics of impulsive and gradual emissions from solar flares

The gradual soft X-ray emissions associated with a homogeneous set of solar flares have been investigated in the context of a thermal model proposed to explain the impulsive components. The parametric techniques which successfully characterized the hard X-ray and microwave observations are employed in an event-by-event analysis to test for quantitative and correlative relationships between the impulsive and gradual emissions. The results of this investigation are consistent with the hypothesis that the hard X-ray and microwave emissions are produced by bulk heating of a common thermal source. The quantitative relationships require an additional source to explain the soft X-ray observations, consistent with previous results. Correlations between the energetics of the impulsive and gradual emissions, identified in the present work, provide the first clear evidence that their energizing mechanisms are related.

Crannell, C. J.↗

Thermal stress effects in intermetallic matrix composites

Intermetallic matrix composites develop residual stresses from the large thermal expansion mismatch (delta-alpha) between the fibers and matrix. This work was undertaken to: establish improved techniques to measure these thermal stresses in IMC's; determine residual stresses in a variety of IMC systems by experiments and modeling; and, determine the effect of residual stresses on selected mechanical properties of an IMC. X ray diffraction (XRD), neutron diffraction (ND), synchrotron XRD (SXRD), and ultrasonics (US) techniques for measuring thermal stresses in IMC were examined and ND was selected as the most promising technique. ND was demonstrated on a variety of IMC systems encompassing Ti- and Ni-base matrices, SiC, W, and Al2O3 fibers, and different fiber fractions (Vf). Experimental results on these systems agreed with predictions of a concentric cylinder model. In SiC/Ti-base systems, little yielding was found and stresses were controlled primarily by delta-alpha and Vf. In Ni-base matrix systems, yield strength of the matrix and Vf controlled stress levels. The longitudinal residual stresses in SCS-6/Ti-24Al-llNb composite were modified by thermomechanical processing. Increasing residual stress decreased ultimate tensile strength in agreement with model predictions. Fiber pushout strength showed an unexpected inverse correlation with residual stress. In-plane shear yield strength showed no dependence on residual stress. Higher levels of residual tension led to higher fatigue crack growth rates, as suggested by matrix mean stress effects.

Wright, P. K.↗

Urban and regional land use analysis: CARETS and Census Cities experiment package

The author has identified the following significant results. The arrival of the so-called energy crisis makes the portion of this experiment dealing with land use climatology of more immediate significance than before, since in addition to helping to understand the processes of climatic change associated with urbanization, the knowledge obtained may be useful in assigning an energy balance impact factor to proposed changes in land use in and around cities. Thermal maps derived from S-192 data are to be used as a measure of the energy being radiated into space from the mosaic of different surfaces in and around the city. While presenting excellent spatial sampling potential for a metropolitan area tests site, the Skylab data permit a very poor temporal sampling opportunity, owing to the large number of factors beyond the investigator's control that determine when data will be taken over a given test site. The strategy is to augment the thermal maps derived from S-192 with a modeling technique which enables the simulation of a number of components of the surface energy balance, calculated at regular time intervals throughout the day or year. Preliminary tests on the performance of the model are still underway, using airborne MSS data from NASA aircraft flights. Results look extremely promising.

Alexander, R. H.↗

Modeling the curing process of thick-section autoclave cured composites

Temperature gradients are significant during cure of large area, thick-section composites. Such temperature gradients result in nonuniformly cured parts with high void contents, poor ply compaction, and variations in the fiber/resin distribution. A model was developed to determine the temperature distribution in thick-section autoclave cured composites. Using the model, long with temperature measurements obtained from the thick-section composites, the effects of various processing parameters on the thermal response of the composites were examined. A one-dimensional heat transfer model was constructed for the composite-tool assembly. The governing differential equations and associated boundary conditions describing one-dimensional unsteady heat-conduction in the composite, tool plate, and pressure plate are given. Solution of the thermal model was obtained using an implicit finite difference technique.

Loos, A. C.↗

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗