SEARCH · Search NASA

Results for “Distributed temperature”

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 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Low-frequency modes and nonbarotropic effects in pseudo-Newtonian accretion disks

A recently developed formalism is used to reexamine the question of the existence of hydrodynamical modes that pulsate with very low frequencies in the inner regions of accretion disks. The formalism is valid in an exact sense for the adiabatic pulsations of rotating Newtonian fluids that are generally nonbarotropic (such as those with 'nonadiabatic temperature gradients,' for example), and hence its application in the present context represents an improvement over previous analyses that are more approximate. The formalism is applied to thin non-self-gravitating disks, with the gravitational potential of the central source modified in the usual way in order to simulate relativistic effects. In the barotropic limit, the analyses indicate that in many cases nearly Keplerian disks exhibit nonaxisymmetric modes of pulsation that are trapped in the inner disk regions, with pulsation periods much longer than the dynamical timescale. These results are similar to those of earlier calculations that assume disks pulsate without changing the temperature distribution. A method is developed for including lowest order nonbarotropic effects. Previous analyses have been incapable of accurately treating the nonbarotropic regime. The application of the present method to the low-frequency modes reveals that, due to unexpected cancellations among terms, the nonbarotropic correction to the pusation frequency omega is only of order tilde-omega(sub BV exp 2) omega, where tilde-omega(sub BV) is the appropriate dimensionless Brunt-Vaisala frequency. This correction is much smaller than the expected correction of order tilde-omega(sub BV) Omega, where Omega is the rotation angular velocity. The important conclusion drawn from this is that nonbarotropic corrections are generally small and hence that low-frequency modes persist into the nonbarotropic regime. For disk temperatures appropriate to X-ray emission, the adiabatic frequencies of trapped modes are of the same order as the frequencies associated with quasi-periodic oscillations (QPOs) of certain low-mass X-ray binaries and other X-ray sources. It is suggested that further consideration be given to the possiblity that modulation of emission processes by these modes can explain at least some of the QPO phenomena.

Ipser, James R.↗

NASA NTRS2006

The X-ray Emitting Components towards l = 111 deg: The Local Hot Bubble and Beyond

We have obtained an XMM-Newton spectrum of the diffuse X-ray emission towards (l, b) = (111.14,1.11), a line of sight with a relatively simple distribution of absorbing clouds; > 9 x 10(exp 19)/sq cm at R>170 pc, a 6 x 10(exp 21)/sq cm molecular cloud at 2.5-3.3 kpc, and a total column of 1.2 x 10(exp 22)/sq cm. We find that the analysis of the XMM-Newton spectrum in conjunction with the RASS spectral energy distribution for the same direction requires three thermal components to be well fit: a "standard" Local Hot Bubble component with kT = 0.089, a component beyond the molecular cloud with kT = 0.59, and a component before the molecular cloud with kT = 0.21. The strength of the O VII 0.56 keV line from the Local Hot Bubble, 2.1+/-0.7 photons/sq cm/s/sr, is consistent with other recent measures. The 0.21 keV component has an emission measure of 0.0022+/-0.0006 pc and is not localized save as diffuse emission within the Galactic plane; it is the best candidate for a pervasive hot medium. The spatial separation of the approx. 0.2 keV component from the approx. 0.6 keV component suggests that the spectral decompositions of the emission from late-type spiral disks found in the literature do represent real temperature components rather than reflecting more complex temperature distributions.

Kuntz, K. D.↗

DOE OSTI2021

Transition to Online Cable Insulation Condition Monitoring

Nuclear power plant cables were originally qualified for 40 year life and generally have not required specific test verification to assure service availability through the initial plant qualification period. However, license renewals to 60 and 80 years of operation require a cable aging management program that depends on some form of test and verification to assure fitness for service. Environmental stress (temperature, radiation, chemicals, water, and mechanical) varies dramatically within a nuclear power plant and, in some cases, cables have degraded and required repair or replacement before their qualified end-of-life period. In other cases, cable conditions have been mild and dependable cable performance confirmed to extend well beyond the initial qualified life. Most offline performance-based testing requires cables to be de-coupled and de-energized for specially trained technicians to perform testing. These offline tests constitute an expensive operational burden that limits the economic viability of nuclear power plants. Although initial investment may be higher, new online test practices are emerging as options or complements to offline testing that avoid or minimize the regularly scheduled offline test burden. These online methods include electrical and fiber-optic partial discharge measurement, spread spectrum time or frequency domain reflectometry, distributed temperature profile measurements, and local interdigital capacitance measurement of insulation characteristics. Introduction of these methods must be supported by research to confirm efficacy plus either publicly financed or market driven investment to support the start-up expense of cost-effective instrumentation to monitor cable condition and assure reliable operation. This work summarizes various online cable assessment technologies plus introduces a new cable motor test bed to assess some of these technologies in a controlled test environment.

Glass, Samuel W.↗

DOE OSTI2025

Heat load measurements for the PIP-II pHB650 cryomodule

This study presents a brief overview of the 1st and 2nd phases and an in-depth analysis of the 3rd phase heat load testing performed on the pHB650 (prototype High Beta 650 MHz) cryomodule at PIP2IT (PIP-II Injector Test Facility), with a focus on both the results and the methodological advancements that have improved testing efficiency and accuracy. A key challenge identified in the testing campaign is the higher-than-expected heat loads observed in the first PIP-II (Proton Improvement Plan II) prototype cryomodules (pSSR1 and pHB650) tested at PIP2IT. Elevated heat loads are concerning given the fixed capacity of the PIP-II cryoplant that is currently being installed at Fermilab. However, understanding the sources of these elevated heat loads offers a critical opportunity to implement effective heat load mitigations on upcoming PIP-II cryomodules to stay within the available capacity of the PIP-II cryoplant. The study includes a summary of test results, descriptions of measurement procedures, and key observations on parameters directly and indirectly related to heat load measurements. Direct observations include measured heat loads and the effectiveness of JT heat exchanger under varying conditions, while indirect observation analyze factors such as the temperature distribution on the two-phase pipe and relief piping under varying conditions. Thermal acoustic oscillations (TAO) were identified during testing, which was mitigated by replacing the original G10 stem with a stainless steel stem equipped with wipers for the cryomodule cooldown valve. A major innovation during pHB650 Phase 3 testing was the development of an automated Python script to streamline data acquisition, analysis, and reporting of heat load results. This script automatically retrieved data from ACNET (Accelerator Control Network), performed heat load calculations, and generated detailed reports featuring plots and tables. This advancement significantly reduced manual labor and enhanced the thoroughness of data analysis compared to earlier campaigns. The heat load test reports were promptly uploaded to the electronic logbook shortly after each test, enabling rapid feedback and collaboration between the SRF and cryogenic teams. The heat load measurements included various components: HTTS (high-temperature thermal shield), LTTS (low-temperature thermal shield), 2K isothermal and non-isothermal heat loads. Results were recorded both within the cryomodule and between the bayonet can supply and return. Measurements were conducted under different operating conditions such as "standard", "linac", and "simulated dynamic". Additionally, HTTS and LTTS heat loads were calculated in real time, allowing for the tracking of thermal stability and identification of changes during testing, both in steady-state and transient conditions. The results of this testing campaign not only provide valuable insights into the performance of the pHB650 cryomodule but also highlight best practices and lessons learned that will inform future cryomodule testing at PIP2IT. These include adopting automated tools for data analysis, refining real-time measurement capabilities, and emphasizing detailed pre-test planning. The framework established in this campaign aims to set an improved standard for cryomodule testing and heat load reporting in future cryomodule test campaigns.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.]↗

NASA NTRS1972

Survey of wing and flap lower-surface temperatures and pressures during full-scale ground tests of an externally blown flap system

Full-scale ground tests of an externally blown flap system were made using the wing of an F-111B airplane and a CF700 engine. Pressure and temperature distributions were determined on the undersurface of the wing, vane, and flap for two engine exhaust nozzles (conical and daisy) at several engine power and engine/wing positions. The tests were made with no airflow over the wing. The leading-edge wing sweep angle was fixed at 26 deg, the angle of incidence between the engine and the wing was fixed at 3 deg, and the tests were conducted with the flap retracted, extended and deflected 35 deg, and extended and deflected 60 deg. The integrated local pressures on the undersurface of the flap produced loads approximately three times as great at the 60 deg flap position as at the 35 deg flap position. With both nozzle configurations, more than 90 percent of the integrated pressure loads were contained within plus or minus 20 percent of the flap span centered around the engine exhaust centerline. The maximum temperature recorded on the flaps was 218 C (424 F) for the conical nozzle and 180 C (356 F) for the daisy nozzle.

Hughes, D. L.↗

NASA NTRS2005

Navier-Stokes Simulation of Airconditioning Facility of a Large Modem Computer Room

NASA recently assembled one of the world's fastest operational supercomputers to meet the agency's new high performance computing needs. This large-scale system, named Columbia, consists of 20 interconnected SGI Altix 512-processor systems, for a total of 10,240 Intel Itanium-2 processors. High-fidelity CFD simulations were performed for the NASA Advanced Supercomputing (NAS) computer room at Ames Research Center. The purpose of the simulations was to assess the adequacy of the existing air handling and conditioning system and make recommendations for changes in the design of the system if needed. The simulations were performed with NASA's OVERFLOW-2 CFD code which utilizes overset structured grids. A new set of boundary conditions were developed and added to the flow solver for modeling the roomls air-conditioning and proper cooling of the equipment. Boundary condition parameters for the flow solver are based on cooler CFM (flow rate) ratings and some reasonable assumptions of flow and heat transfer data for the floor and central processing units (CPU) . The geometry modeling from blue prints and grid generation were handled by the NASA Ames software package Chimera Grid Tools (CGT). This geometric model was developed as a CGT-scripted template, which can be easily modified to accommodate any changes in shape and size of the room, locations and dimensions of the CPU racks, disk racks, coolers, power distribution units, and mass-storage system. The compute nodes are grouped in pairs of racks with an aisle in the middle. High-speed connection cables connect the racks with overhead cable trays. The cool air from the cooling units is pumped into the computer room from a sub-floor through perforated floor tiles. The CPU cooling fans draw cool air from the floor tiles, which run along the outside length of each rack, and eject warm air into the center isle between the racks. This warm air is eventually drawn into the cooling units located near the walls of the room. One major concern is that the hot air ejected to the middle isle might recirculate back into the cool rack side and cause thermal short-cycling. The simulations analyzed and addressed the following important elements of the computer room: 1) High-temperature build-up in certain regions of the room; 2) Areas of low air circulation in the room; 3) Potential short-cycling of the computer rack cooling system; 4) Effectiveness of the perforated cooling floor tiles; 5) Effect of changes in various aspects of the cooling units. Detailed flow visualization is performed to show temperature distribution, air-flow streamlines and velocities in the computer room.

Source record↗

DOE OSTI2023

High-Resolution Thermal Profiling of a High-Power Diode Laser Facet During Aging

Here, we study the facet temperature distribution of a high-power diode laser over its lifetime using a noncontact, high spatial resolution CCD-based thermoreflectance technique. Based on the known correlation between non-radiative defects and heating, thermal maps can provide valuable information regarding the formation and evolution of small point defects that are at or near the facet during aging. In the laser under study in this work we measure the appearance of local hot spots on the facet, including concentrated hot spots that appear just before or just after COD and are correlated with loss of local light emission. The locations of these hot spots do not exhibit morphology changes in high-resolution SEM imaging of the facet, indicating that the related defects are too small to be observable in SEM or are located at some depth under the facet. Prior to COD, we measure a gradual facet temperature increase accompanied by a gradual optical power decrease and gradual facet optical absorption increase, indicating gradual degradation of the laser.

47 OTHER INSTRUMENTATION↗

NASA NTRS1989

Infrared Emission from Comets

A brief discussion of the infrared observations from 4 to 20 micrometers of seven comets is presented. The observed infrared emission from comets depends primarily on their heliocentric distance. A model based on grain populations composed of a mixture of silicate and amorphous carbon particles in the mass ratio of about 40 to 1, with a power-law size distribution similar to that inferred for comet Halley, is applied to the observations. The model provides a good match to the observed heliocentric variation of both the 10 micrometers feature and the overall thermal emission from comets West and Halley. Matches to the observations of comet IRAS-Araki-Alcock and the antitail of comet Kohoutek require slightly larger grains. While the model does not match the exact profile and position of the 3.4 micrometers feature discovered in comet Halley, it does produce a qualitative fit to the observed variation of the feature's strength as a function of heliocentric distance. The calculations predict that the continuum under the 3.4 micrometers feature is due primarily to thermal emission from the comet dust when the comet is close to the Sun and to scattered solar radiation at large heliocentric distances, as is observed. A brief discussion of the determination of cometary grain temperatures from the observed infrared emission is presented. It is found that the observed shape of the emission curve from about 4 to 8 micrometers provides the best spectral region for estimating the cometary grain temperature distribution.

Swamy, K. S. Krishna↗

DOE OSTI2023

Computational fluid dynamics investigations of flow, heat transfer, and oxidation in heat recovery steam generator

Modern heat recovery steam generators (HRSGs) operate at elevated temperatures, leading to the formation of oxides inside the tubes of heat exchangers (HXs). This oxide growth reduces the heat recovery efficiency. Moreover, after reaching a certain critical thickness, some oxide scales detach from the tube surface (exfoliation), causing erosion damage to the components downstream. Predicting the metal temperature distribution and associated oxide thickness in the HXs of an HRSG can aid in mitigating these problems. Here, a computational fluid dynamics (CFD) model was developed within the commercial code STAR-CCM+ for the prediction of fluid flow, conjugate heat transfer, and associated oxidation in HRSGs. Moreover, a new Porous Media Model (PMM) method was developed to model the fin effect on the heat transfer in HX, which can substantially reduce the prohibitive computational costs of fin meshes. The developed CFD model was used to conduct a high-fidelity simulation of a real-scale HRSG to investigate flow, heat transfer, and oxide growth. The calculated oxide thickness on different tubes can be used to identify HX regions that require oxide-resistant coatings to prevent exfoliation and ensuing damages. Furthermore, this CFD framework can serve as a reference for future studies that intend to model and investigate high-temperature oxidation in HXs used for any applications.

42 ENGINEERING↗

DOE OSTI2018

Experimental Study on Portable Air-Conditioning System with Enhanced PCM Condenser

A portable personal air-conditioning system with a vapor compression cycle (VCC), electrical battery, and thermal battery is an attractive option for localized cooling as compared with stationary air-conditioning systems. A phase change material (PCM) can be used for storing waste heat from the VCC condenser, thereby acting as a thermal battery. In the PCM condenser, copper tubes are submerged in the PCM container. Due to the low thermal conductivity of pure PCM, the overall heat transfer coefficient of the PCM condenser is limited resulting in high VCC condensing pressure and thus, low VCC coefficient of performance (COP). Therefore, it is necessary to improve VCC system performances by applying heat transfer augmentation in PCM. In this work, a portable air-conditioning system with heat transfer enhanced PCM that stores heat from condenser during VCC operation and is regenerated by a thermosyphon operation is proposed. The system also includes a receiver to address the charge difference between the thermosyphon and the VCC operation. The paper experimentally studies the effects of the following three types of enhanced PCM condensers on the system performance: (1) copper mesh enhanced PCM with helically coiled copper tube branches and large refrigerant distribution headers (CMPHX), (2) graphite matrix enhanced PCM with straight copper tube branches and small headers (GMPHX1), and (3) graphite matrix enhanced PCM with straight copper tube branches and large headers (GMPHX2). Results show that this system operates properly and effectively during both VCC and PCM regeneration operations. The system can provide four hours of continuous cooling, and the PCM can be fully regenerated after 7.4 hours by thermosyphon operation. The comparisons of three PCM condensers indicate that the system with GMPHX2 performs the best with the average VCC COP of 4.7. The VCC system with CMPHX has an average COP of 4.2, and GMPHX1 has an average COP of 3.1. GMPHX2 has the best heat transfer enhancement with an overall heat transfer coefficient of 632.1 W∙m-2∙K-1, which is more than twice of GMPHX1 and five times of CMPHX. Moreover, GMPHX2 can achieve more uniform refrigerant flows and PCM temperature distributions than GMPHX1. This leads to the higher utilization ratio of PCM latent heat. CMPHX is the second best but needs larger heat transfer surface area and inner volume as well as more refrigerant charge. In conclusion, the system design with GMPHX2 is recommended due to its small refrigerant charge, high VCC COP, and the use of good heat transfer enhancement material.

Qiao, Yiyuan↗

DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗

NASA NTRS2002

Thermostructural Analysis of Carbon Cloth Phenolic Material Tested at the Laser Hardened Material Evaluation Laboratory

RSRM nozzle liner components have been analyzed and tested to explore the occurrence of anomalous material performance known as pocketing erosion. Primary physical factors that contribute to pocketing seem to include the geometric permeability, which governs pore pressure magnitudes and hence load, and carbon fiber high temperature tensile strength, which defines a material limiting capability. The study reports on the results of a coupled thermostructural finite element analysis of Carbon Cloth Phenolic (CCP) material tested at the Laser Hardened Material Evaluation Laboratory (the LHMEL facility). Modeled test configurations will be limited to the special case of where temperature gradients are oriented perpendicular to the composite material ply angle. Analyses were conducted using a transient, one-dimensional flow/thermal finite element code that models pore pressure and temperature distributions and in an explicitly coupled formulation, passes this information to a 2-dimensional finite element structural model for determination of the stress/deformation behavior of the orthotropic fiber/matrix CCP. Pore pressures are generated by thermal decomposition of the phenolic resin which evolve as a multi-component gas phase which is partially trapped in the porous microstructure of the composite. The nature of resultant pressures are described by using the Darcy relationships which have been modified to permit a multi-specie mass and momentum balance including water vapor condensation. Solution to the conjugate flow/thermal equations were performed using the SINDA code. Of particular importance to this problem was the implementation of a char and deformation state dependent (geometric) permeability as describing a first order interaction between the flow/thermal and structural models. Material property models are used to characterize the solid phase mechanical stiffness and failure. Structural calculations were performed using the ABAQUS code. Iterations were made between the two codes involving the dependent variables temperature, pressure and across-ply strain level. Model results comparisons are made for three different surface heat rates and dependent variable sensitivities discussed for the various cases.

Clayton, J. Louie↗

NASA NTRS2001

Measurement of distributed strain and temperature based on higher order and higher mode Bragg conditions

A Bragg grating sensor for measuring distributed strain and temperature at the same time comprises an optical fiber having a single mode operating wavelength region and below a cutoff wavelength of the fiber having a multimode operating wavelength region. A saturated, higher order Bragg grating having first and second order Bragg conditions is fabricated in the optical fiber. The first order of Bragg resonance wavelength of the Bragg grating is within the single mode operating wavelength region of the optical fiber and the second order of Bragg resonance wavelength is below the cutoff wavelength of the fiber within the multimode operating wavelength region. The reflectivities of the saturated Bragg grating at the first and second order Bragg conditions are less than two orders of magnitude of one another. In use, the first and second order Bragg conditions are simultaneously created in the sensor at the respective wavelengths and a signal from the sensor is demodulated with respect to each of the wavelengths corresponding to the first and second order Bragg conditions. Two Bragg conditions have different responsivities to strain and temperature, thus allowing two equations for axial strain and temperature to be found in terms of the measure shifts in the primary and second order Bragg wavelengths. This system of equations can be solved for strain and temperature.

Sirkis, James S.↗

NASA NTRS1973

Pressures and temperatures on the lower surfaces of an externally blown flap system during full-scale ground tests

Full-scale ground tests of an externally blown flap system were made using the wing of an F-111B airplane and a CF700 engine. Pressure and temperature distributions were determined on the undersurface of the wing, vane, and flap for two engine exhaust nozzles (conical and daisy) at several engine power levels and engine/wing positions. The test were made with no airflow over the wing. The wing sweep angle was fixed at e6 deg; and the angle of incidence between the engine and the wing was fixed at 3 deg; and the flap was in the retracted, deflected 35 deg, and deflected 60 deg positions. The pressure load obtained by integrating the local pressures on the undersurface of the flap, F sub p was approximately three times greater at the 60 deg flap position than at the 35 deg flap position. At the 60 deg flap position, F sub p was between 40 percent and 55 percent of the engine thrust over the measured range of thrust. More than 90 percent of F sub p was contained within plus or minus 20 percent of the flap span centered around the engine exhaust centerline with both nozzle configurations. Maximum temperatures recorded on the flaps were 218 C (424 F) and 180 C (356 F) for the conical and daisy nozzles, repectively,

Hughes, D. L.↗

NASA NTRS2014

Transient Convection Due to Imposed Heat Flux: Application to Liquid-Acquisition Devices

A model problem is considered that addresses the effect of heat load from an ambient laboratory environment on the temperature rise of liquid nitrogen inside an enclosure. This model has applications to liquid acquisition devices inside the cryogenic storage tanks used to transport vapor-free propellant to the main engine. We show that heat loads from Q = 0.001 to 10 W, with corresponding Rayleigh numbers from Ra = 109 to 1013, yield a range of unsteady convective states and temperature rise in the liquid. The results show that Q = 1 to 10 W (Ra = 1012 to 1013) yield temperature distributions along the enclosure height that are similar in trend to experimental measurements. Unsteady convection, which shows selfsimilarity in its planforms, is predicted for the range of heat-load conditions. The onset of convection occurs from a free-convection-dominated base flow that becomes unstable against convective instability generated at the bottom of the enclosure while the top of the enclosure is convectively stable. A number of modes are generated with small-scale thermals at the bottom of the enclosure in which the flow selforganizes into two symmetric modes prior to the onset of the propagation of the instability. These symmetric vertical modes transition to asymmetric modes that propagate as a traveling-wave-type motion of convective modes and are representative of the asymptotic convective state of the flow field. Intense vorticity production is created in the core of the flow field due to the fact that there is shear instability between the vertical and horizontal modes. For the higher Rayleigh numbers, 1012 to 1013, there is a transition from a stationary to a nonstationary response time signal of the flow and temperature fields with a mean value that increases with time over various time bands and regions of the enclosure.

Duval, Walter M. B.↗

DOE OSTI2019

A Forward Time Stepping Heat Conduction Model for Spot Melt Additive Manufacturing

Solidification dynamics are crucial for determining microstructure development in additively manufactured parts. Multiphysics models based on finite element or finite volume methods may help gain insight for complicated phenomena such as fluid flow, keyholing, and porosity but are too computationally expensive to use for simulating actual builds. Recent analytic and semi-analytic solutions for moving heat sources in a semi-infinite three-dimensional space provide a way to accurately estimate the solidification conditions for entire builds. The downside to these methods is that, unlike finite element or finite volume methods, they cannot use the temperature distribution of the previous timesteps to march the solution forward in time. This paper provides the mathematical formulation and implementation of a forward time stepping (FTS) approach to an existing semi-analytic solution. The speed and accuracy of the two methods are then compared for various scan patterns. The result is that, for spot melts, the forward time-stepping model provides improvements in both speed and accuracy. This is especially true for longer simulations, where the simulation can be orders of magnitude faster. Finally, the longest simulation analyzed in this paper was roughly 30× faster when using the forward time-stepping model versus the straightforward implementation of the semi-analytic solution.

36 MATERIALS SCIENCE↗

NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗

DOE OSTI2022

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗