Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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

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

At least 1,549 records · Page 86

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Comparison of Heat Transfer from Airfoil in Natural and Simulated Icing Conditions

An investigation of the heat transfer from an airfoil in clear air and in simulated icing conditions was conducted in the NACA Lewis 6- by 9-foot icing-research tunnel in order to determine the validity of heat-transfer data as obtained in the tunnel. This investiation was made on the same model NACA 65,2-016 airfoil section used in a previous flight study, under similar heating, icing, and operating conditions. The effect of tunnel turbulence, in clear air and in icingwas indicated by the forward movement of transition from laminar to turbulent heat transfer. An analysis of the flight results showed the convective heat transfer in icing to be considerably different from that measured in clear air and. only slightly different from that obtained in the icing-research tunnel during simulated icing.

Gelder, Thomas F.↗

Experimental feasibility study of a thermoelectric heat flux gage

An experiment was conducted to determine the feasibility of using a commercially available thermoelectric device as a heat flux gage at near ambient conditions. In certain research applications, the thermoelectric heat flux gage can provide a relatively simple means to model a warm fluid cold wall convection environment. The experiment showed that heat flux through the gage could be correlated within 2.5 percent with a simple algebraic equation which considered the thermoelectric current through the device and the hot and cold side temperatures.

Vanfossen, G. J.↗

Effect of an oscillating flow direction on leading edge heat transfer

An experimental investigation was conducted to examine the effect of a periodic variation in the angle of attack on heat transfer at the leading edge of a gas turbine blade. A circular cylinder was used as a large-scale model of the leading edge region. The cylinder was placed in a wind tunnel and was oscillated rotationally about its axis. The incident flow Reynolds number and the Strouhal number of oscillation were chosen to model an actual turbine condition. Incident turbulence levels up to 4.9 percent were produced by grids placed upstream of the cylinder. The transfer rate was measured using a mass transfer technique and heat transfer rates inferred from the results. A direct comparison of the unsteady and steady results indicate that the effect is dependent on the Strouhal number, turbulence level, and the turbulence length scale, but that the largest observed effect was only a 10 percent augmentation at the nominal stagnation position.

Marziale, M. L.↗

A spacecraft cooling system for a charged coupled device

This paper describes the thermal analysis, design, and testing of a dedicated cooling system for a Spartan spacecraft payload. A simple reliable design that requires minimum power consumption and minimum weight was developed. The payload has a CCD detector that must be maintained at a temperature of approximately -40 C or colder. The cooling system consists of a fin radiator, dual redundant heat pipes, and a thermal electric device (TED). The system was analytically modeled through the use of the Simplified Shuttle Payload Thermal Analyzer (SSPTA) computer program. A thermal test of the system simulating flight conditions was conducted to correlate the computer model and verify performance specifications.

Walker, Mary S.↗

Natural melting within a spherical shell

Fundamental heat transfer experiments were performed on the melting of a phase change medium in a spherical shell. Free expansion of the medium into a void space within the sphere was permitted. A step function temperature jump on the outer shell wall was imposed and the timewise evolution of the melting process and the position of the solid-liquid interface was photographically recorded. Numerical integration of the interface position data yielded information about the melted mass and the energy of melting. It was found that the rate of melting and the heat transfer were significantly affected by the movement of the solid medium to the base of the sphere due to gravity. The energy transfer associated with melting was substantially higher than that predicted by the conduction model. Furthermore, the radio of the measured values of sensible energy in the liquid melt to the energy of melting were nearly proportional to the Stefan number. The experimental results are in agreement with a theory set forth in an earlier paper.

Bahrami, Parviz A.↗

Thermographic inspection of fatigue cracks in riveted plates

A theoretical and experimental analysis of the inspection of reiveted plates using the thermoelastic effect is presented. Theoretically, a model is employed to estimate the thermoelastic response expected from a cracked and uncracked rivet hole. The stress invariant, which causes the thermoelastic response, is calculated for the cracked case using an FEM model and for the uncracked case using a theoretical approximation. The differences in the two thermal responses are reduced when the effects of optical spatial averaging, thermal conduction, and a concentric blind area representing the rivet heat are considered. It is suggested that discrimination between the cracked and uncracked cases for short cracks will improve with increasing frequency of stress oscillation. Measurements are shown for aerospace aluminum samples with and without fatigue cracks propagating from holes in a plate. The effect of mean stress on the thermoelastic response is illustrated, and the results are qualitatively compared with the finite element analysis.

Welch, Christopher S.↗

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗

Method of Forming a Hot Film Sensor System on a Model

A method of forming a hot film sensor directly on a model is provided. A polyimide solution is sprayed onto the model. The model so sprayed is then heated in air. The steps of spraying and heating are repeated until a polyimide film of desired thickness is achieved on the model. The model with the polyimide film thereon is then thoroughly dried in air. One or more hot film sensors and corresponding electrical conducting leads are then applied directly onto the polyimide film.

Tran, Sang Q.↗

Using In-Situ Deposition of Metallic Thin Films on Mars to Monitor Atmospheric H20

We describe the use of a novel sensor, designed to characterize the reactive nature of the Martian atmosphere, to also characterize the instantaneous abundance of atmospheric H,O. The sensor deposits, in situ, a thin silver film onto a sapphire substrate and monitors oxidation by measuring resistance of the Ag film during both deposition and subsequent oxidation [1,2]. . The evaporation source is placed in the center of a hollow, open-ended tube. On flash heating the source, evaporated metallic silver rapidly reacts with oxidizing gases in the tube, depositing on the tube walls as a resisitive Ag oxide film. Unoxidized silver deposits over the Ag oxides, once the oxidizing gases in the vicinity of the source have been consumed. Since the ends of the tube are open, atmospheric gases in the tube cause a time delay as the evaporated silver initially reacts with the available O2, and H2O; once oxidizing gases in the tube are significantly depleted, a silver film closes the chemiresistor circuit. In principle this operation mode provides a repeatable measure of the variable H2O abundance, since O2, levels are constant in the atmosphere. To explore the utility of the instrument as an H2O sensor, predict its behavior, and identify possible failure modes, we developed a numerical model of the instrument, and exercised it for Mars conditions. The model predicts the behavior of the electrical circuit, the temperature dependence of resistivity for each component, and the resultant Joule heating. The model balances Joule heating against radiation, sensible heat loss to ambient CO2, and, latent heat loss from sublimating Ag. The flux of Ag atoms from the source is tracked continuously. In the gas phase, the number density of Ag, O2, and H2O is calculated, along with their reaction rates on collision. The model suggests that Ag substantially depletes H2O only near the center of the tube. The model supports the hypothesis that the onset of high conductivity between the A u electrodes is a function of the H2O abundance in the ambient gas, but only to a limit; when H2O abundances are buffered at or above 240K (an H2O abundance that would not have occurred on Mars), metallic Ag never deposits in excess of Ago, and the circuit does not close. A significant difference in the Mars simulation is that the Ag source sublimates more quickly, primarily because of differences in sensible heat fluxes from the source. This allows subsequent heating of the W filament used to support the Ag source. Had the originally-designed experimental sequence been carried out on Mars, the filament might have melted, making subsequent investigations impossible.

Zent, A. P.↗

Process Heating Assessments Using DOE’s Manufacturing Energy Assessment Software for Utility Reduction (MEASUR) Tool Suite

Process heating is the most energy-intensive manufacturing process for most sectors of industry. To quantify energy savings from various energy conservation measures, the Department of Energy (DOE) sponsored the development of the Process Heating Assessment and Survey Tool (PHAST) and similar tools for other industrial systems in the early 2000s. It has been used extensively in the Save Energy Now Program’s Energy Savings Assessments and the Better Plants Program’s In-Plant Trainings. Since the initial development of the legacy tools, both computer operating systems and software development have evolved significantly. Thus, DOE has invested in the modernization of PHAST and other legacy software tools to create the Manufacturing Energy Assessment Software for Utility Reduction (MEASUR) tool suite. MEASUR offers a collection of software tools that can aid manufacturing facilities in improving the efficiency of energy systems and equipment (specifically pumps, fans, steam, and process heating) and in conducting “Energy Treasure Hunts”. Eventually, the tool will also add compressed air and process cooling systems. The Process Heating Assessment (PHA) module of MEASUR is an upgrade of the PHAST tool. PHA provides the means to model fuel-fired, steam-based, and electric process heating systems, covering process heating for most industrial plants in manufacturing sector. It also includes several key upgrades, including the ability to consider multi-component charge loads and account for several different areas of energy losses. The new tool includes a comprehensive flue gas calculator to quantify available heat and heat loss for various gaseous, liquid, and solid fuels and new heat loss calculators. It generates a report and a dynamic Sankey diagram to show the energy consumption in various areas of energy use. MEASUR has significantly improved the user experience by adopting a modern software design. This paper details the structure and workflow of PHA and presents a real-world case study to demonstrate energy savings quantification and MEASUR’s outstanding reporting capabilities.

Nimbalkar, Sachin U.↗

Morphological-informed Thermal Property Prediction in the Engineering Domain

Using first-principle, atomic scale methods to predict thermal carrier (electron and phonon) behavior at nanoscale yields highly resolved thermal properties. Density functional theory simulations can model myriad effects on the transport of phonons and electrons (e.g., phonon and electron coupling, electron correlation, defect presence, carrier scattering), and yield thermal properties for small geometric domains such as thermal and electrical conductivity, or heat capacity. While these thermal properties are accurate in atomistic and nanoscopic systems, they must be scaled up to the microstructural domain to be a useful predictor for experimental basis, as the microscale is where myriad changes and physical phenomenon occur within a material (e.g., grain boundaries, precipitate aggregations, interfaces, defect clusters). My talk will discuss the necessity of thermal property predictions for nuclear fuels applications and more. Additionally I will discuss the necessity in developing multi-physics, multi-scale methods for future applications in predictions for not only nuclear fuels, but other materials applications.

36 MATERIALS SCIENCE↗

Measurement of Effective Canopy Temperature: The Missing Link to Modeling Transpiration in Controlled Environments

The objectives were to apply energy balance principles to plant canopies, and to determine which parameters are essential for predicting plant canopy transpiration (E) in controlled environments. Transpiration was accurately measured in a gas-exchange system. Absorbed radiation (R(sub abs)) by the canopy was measured with a net radiometer and calculated from short and long-wave radiation components. Average canopy foliar temperature T(sub L) can be measured with an infrared radiometer, but since T(sub L) is seldom uniform, a weighed average measurement of T(sub L) must be made. The effective canopy temperature T(sub C) is that temperature that balances the energy flux between absorbed radiation and latent heat L(sub E) and sensible heat (H) fluxes. TC should exactly equal air temperature T(sub A) when L(sub E) equals R(sub abs). When unnecessary thermal radiation from the lighting system is removed by a water filter, the magnitude of L(sub E) from transpiration approaches Rabs and T(sub C) is close to T(sub A). Unlike field models, we included the energy used in photosynthesis and found that up to 10% of Rabs was used in photosynthesis. We calculated aerodynamic conductance for H from measurements of wind speed and canopy height using the wind profile equation. Canopy aerodynamic conductance ranged from.03 to.04 m/s for wind speeds from.6 to 1 m/s; thus a 0.1 C canopy to air temperature difference results in a sensible heat flux of about 4 W/sq m, which is only 1% of R(sub abs). We examined the ability of wide angle infrared transducers to accurately integrate T(sub L) from the top to the bottom of the canopy. We measured evaporation from the hydroponic media to be approximately 1 micro mol/sq m s or 10% of R(sub abs). This result indicates that separating evaporation from transpiration is more important than exact measurement of canopy temperature.

Monje, O. A.↗

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↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗

DEM Modeling and Validation of Pebble Bed Packing Using Chrono::GPU

Accurate prediction of pebble packing structure is important for pebble bed reactors because the spatial distribution of void fraction directly affects coolant flow, pressure drop, heat transfer, and neutronic behavior. However, experimentally validated DEM studies that directly evaluate local void-fraction structure in reactor-relevant pebble beds remain limited. In this work, the pebble bed experiment conducted at Missouri University of Science and Technology is simulated using the graphics processing unit (GPU)-based discrete element method (DEM) code Chrono::GPU. The study focuses on evaluating the ability of Chrono::GPU to reproduce the packing arrangement and void-fraction distribution of a randomly packed spherical pebble bed. The DEM results are first verified against established radial void-fraction correlations, including the Mueller and Vortmeyer-Schuster models, to assess the predicted bulk porosity, near-wall behavior, and oscillatory packing structure. The simulation is then verified against reference DEM data and validated against gamma-ray computed tomography (CT) experimental data at three axial locations. The Chrono::GPU results reproduce the main features of the experimental packing, including the high void fraction near the wall, the first near-wall trough, and the damped oscillatory radial profile caused by wall-induced ordering. Quantitative comparison with DEM data and the CT-based radial profiles shows good agreement, with mean absolute errors on the order of 0.07 and root-mean-square errors below 0.09 for the averaged profiles. These results demonstrate that Chrono::GPU can accurately capture the void-fraction structure of spherical pebble beds and provides a reliable DEM framework for future pebble bed reactor packing, recycling, and thermal-hydraulic studies.

97 - MATHEMATICS AND COMPUTING↗

Temperature distribution in a laser-heated diamond anvil cell as described by finite element analysis

Finite element analysis (FEA) is a powerful tool for numerically solving partial differential equations over complex geometries and is thus useful for analyzing heat transport in laser-heated diamond anvil cell (LHDAC) experiments. Our models expand on previously published simulations by calculating the volume-averaged temperatures of both the sample and insulation/pressure media under steady-state heating to determine the thermal pressure of the hot sample. Our goal is to produce an accurate relationship between the measured surface temperature of the absorbing sample and the temperature of the transparent insulating media, which is used to determine thermal pressure but susceptible to steep temperature gradients. We find that in doing so, our FEA models of temperature within the pressure/insulation media can differ from simplified estimates of temperature gradients by more than a factor of 2. We also explore temperature-dependent and temperature-independent thermal conductivity models and find that the volume-averaged temperatures differ by up to a factor of 1.3, forcing the predicted thermal pressures determined to also differ by up to a factor of 1.5 at a temperature of 2000 K at 50 GPa for neon. Higher temperatures exacerbate this difference. We also find that unintentional asymmetric sample insertion and sample heating, which are common in LHDAC experiments, do not have a first-order effect on volume-averaged temperatures. The FEA models, available in both Python and FlexPDE, are versatile across different sample geometries, materials, and heat source laser shapes.

Farah, Frederick↗