SEARCH · Search NASA

Results for “effective conductivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Code Benchmark of the HTTF Pressurized Conduction Cooldown Test Using SAM

The High Temperature Test Facility (HTTF) at Oregon State University is an integral system test facility to simulate postulated reactor transients of prismatic high-temperature gas-cooled reactors(HTGRs). A series of test campaigns was launched, providing abundant test data that could be used to benchmark reactor system analysis codes like the System Analysis Module (SAM). In this study, a SAM model of the facility is developed based on the two-dimensional (2D) ring model approach. All components including the ceramic matrix, graphite heaters, helium coolant channels, core barrel, upcomer, pressure vessel, and reactor cavity cooling system are modeled as concentric cylindrical rings. The model is used to simulate one of the benchmark problems-Pressurized Conduction Cooldown (PCC)-within the scope of the Organisation for Economic Co-operation and Development Nuclear Energy Agency International HTTF Benchmark. The simulations consist of two parts. In the first part, operating and boundary conditions as well as thermophysical properties of materials are specified for the benchmark problem. In this work, results from the first part will be used in code-to-code comparison. In the second part, the SAM model is used to simulate Test PG-27, which is the first PCC test carried out in the HTTF, with only two of the ten heater banks activated. The results in the second part are used for code-to-data comparison. Because the helium coolant flow rate is not measured in this facility, it is estimated using the input power and inlet/outlet coolant temperatures. Additionally, radial heat flow in the ceramic blocks is complicated by hundreds of cylindrical coolant channels and heater rods embedded in them. As such, it is necessary to deduce an effective thermal conductivity for the ceramic to analyze the core thermal behavior. SAM predictions of the helium coolant and ceramic temperatures are compared with test data measured in three equivalent sectors. Overall, the SAM results agree reasonably well with test data within the variation of data among the three sectors, which demonstrates SAM's capability in capturing transient effects in HTGR using the simplified 2D ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗

NASA NTRS2009

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

DOE OSTI2022

Giant resonances in topological spin Hall effect due to electron-skyrmion scattering in two-dimensional Rashba spin-orbit ferromagnets

Here we study a topological spin Hall effect where conduction electrons are scattered by Néel-type skyrmions in two-dimensional Rashba spin-orbit ferromagnets. We find the resonance structures in direct and Hall electric conductivities. The resonances strongly depend on the skyrmion size, the Rashba spin-orbit coupling vector value, and the relative direction between magnetization and Rashba vectors. The antiparallel arrangement of the Rashba spin-orbit coupling vector with respect to the magnetization determines the resonance structure and increases the direct conductivity by about two orders of magnitude. For the parallel arrangement the electric conductivities decrease by one to three orders of magnitude and the Hall conductivity changes the sign. The Rashba spin-orbit coupling also modifies electron bands depending on the ratio ε R /J, where ε R is Rashba energy and J is an exchange integral between conduction and localized electrons. If ε R /J ≤ 1, each energy band has the single minimum while for ε R /J > 1 the lower band has the “Mexican-hat” shape with the maximum at k = 0. We focus on the dependencies of direct and Hall electric conductivities on Fermi energy, ε F , skyrmion sizes, and Rashba spin-orbit coupling constant values and its relative direction with the magnetization. Analyzing both types of the energy bands, in general, we find two types of resonances: at the minimum of the upper band and behaviors, we employ the scattering pattern analysis. The resonance dependencies on skyrmion sizes and Rashba, a spin-orbit coupling contestant, can be qualitatively explained in terms of the Ramsauer-Townsend scattering of the upper band electrons by the skyrmion quantum well. The resonance properties can be used in spin transistors. To discover the resonances it is necessary to know the specific range of the parameters. At some values of the parameters the electric conductivity changes by about two orders of magnitude in the narrow range of ε F (< 0.01 eV). To detect the resonances it is also important to identify the relative direction of the Rashba spin-orbit coupling vector with respect to the magnetization.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗

DOE OSTI2022

Uncertainty in Thermal Modeling of Spent Nuclear Fuel Casks

Uncertainty is a key metric in computational modeling that must be evaluated for results to have wide ranging applicability. A well characterized uncertainty range is ideal with clear error bars on results that can be presented to stakeholders. In the field of spent fuel cask modeling, this ideal has been historically difficult to achieve in practice because of the computationally intensive nature of the models used and the difficulty assigning reasonable uncertainties to quantities in as-built systems. The work in this report has been conducted to evaluate the overall state of uncertainty and sensitivity in spent fuel cask models and develop methodologies for evaluating these uncertainties. These methodologies must be practical for engineering applications. They should not require excessive computational resources or calendar time to achieve results. In engineering, the model must be on a scale such that it can be changed and adapted throughout a project as new information is discovered and project goals evolve. This report covers three major modeling task areas that provide an overview of the types of sensitivity and uncertainty present in a spent fuel storage and transportation system. Section 3 discusses sensitivity and uncertainty analysis in the effective thermal conductivity model for the fuel region and applies these results to a single assembly model. Section 4 shows sensitivity analysis of a full cask model in the TN-32B and Section 5 demonstrates the overall uncertainty workflow using Coolant Boiling in Rod Arrays – Spent Fuel Storage and STAR-CCM+ developed from the sensitivity work in the preceding sections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1988

Effects of eddy viscosity and thermal conduction and Coriolis force in the dynamics of gravity wave driven fluctuations in the OH nightglow

The chemical-dynamical model of Walterscheid et al. (1987), which describes wave-driven fluctuations in OH nightglow, was modified to include the effects of both eddy thermal conduction and viscosity, as well as the Coriolis force (with the shallow atmosphere approximation). Using the new model, calculations were performed for the same nominal case as used by Walterscheid et al. but with only wave periods considered. For this case, the Coriolis force was found to be unimportant at any wave period. For wave periods greater than 2 or 3 hours, the inclusion of thermal conduction alone greatly modified the results (in terms of a complex ratio 'eta' which expresses the relationship between the intensity oscillation about the time-averaged intensity and the temperature oscillation about the time-averaged temperature); this effect was reduced with the further inclusion of the eddy viscosity.

Hickey, M. P.↗

DOE OSTI2023

Heat transfer model for moving packed-bed particle-to-sCO 2 heat exchangers integrated with metal foams

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchangers (HXs) play a vital role in coupling heat transfer fluid (HTF) from high-temperature thermal receivers to power cycle working fluids (WF). Heat transfer enhancement is essential for adopting particle-based moving packed-bed heat exchangers (MPBHXs) in next-generation thermal energy storage (TES) systems, as MPBHXs usually exhibit low particle bed-to-wall heat transfer coefficients. High-porosity metal foams have shown their effectiveness in heat transfer enhancement. Here, this work presents a continuum heat transfer model to demonstrate the heat transfer enhancement in MPBHXs when the particle bed channel is filled with high porosity metal foams compared to open channel MPBHXs. The presence of metal foams increases the effective thermal conductivity in the particle channel and enhances the interstitial heat transfer coefficient between the moving particle bed and the stationary metal foams. The present model considers coupled two-dimensional (2D) heat transfer in the particle channel with metal foams and 1D heat transfer in the sCO 2 channel and the dividing wall. The temperature profiles for the particle bed, metal foam, dividing wall, and sCO 2 stream, as well as the local heat flux profiles between those, are studied in detail for various foam porosities, from which the particle bed-to-wall heat transfer coefficient, overall heat transfer coefficient, and total heat transfer rate for the MPBHX are determined. The effects of major MPBHX design and operating parameters on the particle bed-to-wall heat transfer coefficient and overall heat exchange capacity are thoroughly examined; thus, the MPBHX performance improvement with metal foams is quantified for each case. The present heat transfer model can provide valuable insights into the metal-foam MPBHX design and optimization, scale-up, and operating parameters selection.

42 ENGINEERING↗

NASA NTRS2003

Electric Field Effects in Self-Propagating High-Temperature Synthesis under Microgravity Conditions

Self-propagating high-temperature synthesis (SHS) has been used to form many materials. SHS generally involves mixing reactants together (e.g., metal powders) and igniting the mixture such that a combustion (deflagration) wave passes though the mixture. The imposition of an electric field (AC or DC) across SHS reactants has been shown to have a marked effect on the dynamics of wave propagation and on the nature, composition, and homogeneity of the product . The use of an electric field with SHS has been termed "field-assisted SHS". Combustion wave velocities and temperatures are directly affected by the field, which is typically perpendicular to the average wave velocity. The degree of activation by the field (e.g., combustion rate) is related to the current density distribution within the sample, and is therefore related to the temperature-dependent spatial distribution of the effective electrical conductivity of reactants and products. Furthermore, the field can influence other important SHS-related phenomena including capillary flow, mass-transport in porous media, and Marangoni flows. These phenomena are influenced by gravity in conventional SHS processes (i.e., without electric fields). As a result the influence of the field on SHS under reduced gravity is expected to be different than under normal gravity. It is also known that heat loss rates from samples, which can depend significantly on gravity, can influence final products in SHS. This research program is focused on studying field-assisted SHS under reduced gravity conditions. The broad objective of this research program is to understand the role of an electric field in SHS reactions under conditions where gravity-related effects are suppressed. The research will allow increased understanding of fundamental aspects of field-assisted SHS processes as well as synthesis of materials that cannot be formed in normal gravity.

Unuvar, C.↗

NASA NTRS2015

Single-Event Effect Performance of a Conductive-Bridge Memory EEPROM

We investigated the heavy ion SEE characteristics of an EEPROM based on CBRAM technology. SEFI is the dominant type of SEE for each operating mode (standby, read-only, write/read). We also observed single bit upsets in the CBRAM cell, during write/read tests. the SEULET threshold is between 10 and 20 MeV * sq cm/mg, with an upper fluence limit of 3 × 10(exp 6) cm(exp -2) at 10 MeV * sq cm/mg. In the stand by mode, the CBRAM array appears immune to bit upsets.

single event effects↗

DOE OSTI2022

Linking the Salmonella enterica 1,2-Propanediol Utilization Bacterial Microcompartment Shell to the Enzymatic Core via the Shell Protein PduB

Bacterial microcompartments (MCPs) are protein-based organelles that house the enzymatic machinery for metabolism of niche carbon sources, allowing enteric pathogens to outcompete native microbiota during host colonization. While much progress has been made toward understanding MCP biogenesis, questions still remain regarding the mechanism by which core MCP enzymes are enveloped within the MCP protein shell. Here, we explore the hypothesis that the shell protein PduB is responsible for linking the shell of the 1,2-propanediol utilization (Pdu) MCP from Salmonella enterica serovar Typhimurium LT2 to its enzymatic core. Using fluorescent reporters, we demonstrate that all members of the Pdu enzymatic core are encapsulated in Pdu MCPs. We also demonstrate that PduB is critical for linking the entire Pdu enzyme core to the MCP shell. Using MCP purifications, transmission electron microscopy, and fluorescence microscopy, we find that shell assembly can be decoupled from the enzymatic core, as apparently empty MCPs are formed in Salmonella strains lacking PduB. Mutagenesis studies reveal that PduB is incorporated into the Pdu MCP shell via a conserved, lysine-mediated hydrogen bonding mechanism. Finally, growth assays and system-level pathway modeling reveal that unencapsulated pathway performance is strongly impacted by enzyme concentration, highlighting the importance of minimizing polar effects when conducting these functional assays. Together, these results provide insight into the mechanism of enzyme encapsulation within Pdu MCPs and demonstrate that the process of enzyme encapsulation and shell assembly are separate processes in this system, a finding that will aid future efforts to understand MCP biogenesis.

1, 2-propanediol utilization↗

DOE OSTI2023

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗

DOE OSTI2022

Synthesis and electrical behavior of VO 2 thin films grown on SrRuO 3 electrode layers

VO 2 thin films were grown on conducting oxide underlayer SrRuO 3 buffered SrTiO 3 (111) and Si/SiO 2 substrates, respectively, using sputtering. X-ray diffraction phi-scans revealed the epitaxial nature of the VO 2 films grown on SrRuO 3 buffered SrTiO 3 and polycrystalline structure for films grown on SrRuO 3 buffered Si/SiO 2 . X-ray photoelectron spectroscopy confirms a dominant presence of V 4+ in both films and establishes a high-quality growth of single-phase VO 2 films. Temperature and electric-field driven metal-insulator-transition in both the in-plane and out-of-plane configurations were investigated. Depending on the configuration, the resistance change across the metal-insulator-transition varies from a factor of 1.57–3. The measured resistance in each state as well as the magnitude of resistance change were similar during temperature and electric-field driven metal-insulator-transition. To shed light on the suppressed metal-insulator-transition characteristics due to the current shunting effect from conducting SrRuO 3 bottom electrode, a distributed resistance network model is proposed and benchmarked against reports from the literature. Here, the results demonstrate the growth of high-quality VO 2 on conducting SrRuO 3 layers and their electrical behavior, which is of particular interest for all-oxide electronic devices utilizing phase transitions such as resistive memory and neuromorphic oscillators.

36 MATERIALS SCIENCE↗

NASA NTRS2013

Silicon Heat Pipe Array

Improved methods of heat dissipation are required for modern, high-power density electronic systems. As increased functionality is progressively compacted into decreasing volumes, this need will be exacerbated. High-performance chip power is predicted to increase monotonically and rapidly with time. Systems utilizing these chips are currently reliant upon decades of old cooling technology. Heat pipes offer a solution to this problem. Heat pipes are passive, self-contained, two-phase heat dissipation devices. Heat conducted into the device through a wick structure converts the working fluid into a vapor, which then releases the heat via condensation after being transported away from the heat source. Heat pipes have high thermal conductivities, are inexpensive, and have been utilized in previous space missions. However, the cylindrical geometry of commercial heat pipes is a poor fit to the planar geometries of microelectronic assemblies, the copper that commercial heat pipes are typically constructed of is a poor CTE (coefficient of thermal expansion) match to the semiconductor die utilized in these assemblies, and the functionality and reliability of heat pipes in general is strongly dependent on the orientation of the assembly with respect to the gravity vector. What is needed is a planar, semiconductor-based heat pipe array that can be used for cooling of generic MCM (multichip module) assemblies that can also function in all orientations. Such a structure would not only have applications in the cooling of space electronics, but would have commercial applications as well (e.g. cooling of microprocessors and high-power laser diodes). This technology is an improvement over existing heat pipe designs due to the finer porosity of the wick, which enhances capillary pumping pressure, resulting in greater effective thermal conductivity and performance in any orientation with respect to the gravity vector. In addition, it is constructed of silicon, and thus is better suited for the cooling of semiconductor devices.

Yee, Karl Y.↗

DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1981

The effect of interstitial gaseous pressure on the thermal conductivity of a simulated Apollo 12 lunar soil sample

The thermal conductivity of a simulated Apollo 12 soil sample is measured as a function of interstitial gas density, and implications for the thermal properties of lunar and Martian regolith are discussed. Measurements were performed for samples consisting of a mixture of Knippa and Berkely basalt powders with a grain size distribution identical to that of Apollo 12 lunar soil samples by the needle probe technique at interstitial pressures of He, N2, Ar and CO2 from 133,000 to 0.0133 Pa. It is shown that sample thermal conductivity decreases with decreasing interstitial gas pressure down to 1.0 Pa, due to the decreasing effective thermal conductivity of interstitial gas with decreasing gas pressure. Constant thermal conductivity values of 8.8 mW/m per K and 10.9 mW/m per K are obtained for sample densities of 1.70 and 1.85 g/cu cm, respectively, in agreement with in situ lunar regolith measurements. The results, which are greater than those obtained in previous soil studies, are explained by the dense packing of soil particles and enhanced intergranular thermal contact in the present experimental configuration, rather than the influence of interstitial gas pressure. The differences in conductivity between loose soils and packed regolith may also be used to account for the two peaks observed in Martian surface thermal inertia data.

Horai, K.-I.↗

DOE OSTI2024

Comparison of Results between the Legacy and Refined RELAP5-3D Models of the High Temperature Test Facility in Exercises 1 and 2 of the HTTF Benchmark

Work conducted in FY23 identified that RELAP5-3D was capable of reproducing trends in HTTF data during experiment PG-27 but was incapable of reproducing measured values. The primary cause of this discrepancy between RELAP5-3D results and experimental data was hypothesized to be a distortion in power density that was introduced by the radial nodalization of the model. We further hypothesized that a new model would provide better results when compared to the experiments PG-27 and PG-29. Work this FY developed a new model that is better capable of capturing local heat generation rates and contains a representation of each 1/6 azimuthal sector of the core. We used this model to develop a new set of solutions to Exercises 1 and 2 of Problems 2 and 3 in the benchmark. In this report, we present the first comprehensive comparison of the results between the two models. We see that in Exercise 1A, which is common between problems 2 and 3, the results are similar, though the results from the new model show greater detail than those from the legacy model. In Problem 2 Exercise 1B and Problem 3 Exercise 1B, we see that heat removal is slower in the new model than the legacy model. Problem 3 Exercise 1C shows temperatures that are lower in most places in the new model than the legacy model, but the area with active heat generation has higher block temperatures in the new model than the legacy model. Problem 3 Exercise 1D further shows that long-term heat removal is lower in the new model. Problem 2 Exercise 1C demonstrated that the new model observes higher temperatures in the core regions than the legacy model, justifying the need to preserve the power density in HTTF. The validation of PG-27 and PG-29 also demonstrated the improved temperature agreement in the core regions, particularly with a calibrated model that implements an effective thermal conductivity for the core material. Overall, PG-27 models show reasonable to excellent agreement for steady-state temperatures and minimal to reasonable agreement for transients. PG-29 models showed minimal to insufficient agreement with the data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Squeezing a Prism into a Surface: Emulating Bulk Optics with Achromatic Metasurfaces

Metasurfaces promise to replace bulky prisms and lenses with 2D surfaces, revolutionizing wavefront control with technologically significant advantages in size, weight, and planar fabrication. However, conventional implementations suffer from large chromatic aberrations and cannot sustain performance over practical bandwidths of real-world signals because of the limited phase modulation margin available in a surface. How can an infinitely thin surface generate the arbitrarily large, broadband phase delay that bulk phase accumulation can provide? Here, equivalence between bulk optics and certain multiresonant metasurfaces is demonstrated, where phase delay arises from trains of multiple resonances in the effective sheet conductivities of the surface itself instead of accumulation of propagation phase. The fundamentally required electromagnetic surface conductivities of a purely achromatic metasurface are derived and general design rules for arbitrarily broadband beam steering and lensing are obtained. Both operation in transmission and reflection can be achieved in a unified way by proper alignment of the resonances, enabling 360° directionality. The arbitrary spectral bandwidth, dual reflection/transmission operation, and complete dispensability of propagation phase constitute a major advance in the state of the art of achromatic metasurfaces for wavefront manipulation.

36 MATERIALS SCIENCE↗