SEARCH · Engineering Papers
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.
Linear mathematical model predicts increase in temperature in each segment of 15-segment resistive structure used to test electromigration. Assumption of linearity based on fact: equations that govern flow of heat are linear and coefficients in equations (heat conductivities and capacities) depend only weakly on temperature and considered constant over limited range of temperature.
Thermal analyses are performed on the liquid hydrogen (LH2) tank designed for an unmanned aerial vehicle (UAV) powered by solar arrays and a regenerative proton-exchange membrane (PEM) fuel cell. A 14-day cruise mission at a 65,000 ft altitude is considered. Thermal analysis provides the thermal loads on the tank system and the boiling-off rates of LH2. Different approaches are being considered to minimize the boiling-off rates of the LH2. It includes an evacuated multilayer insulation (MLI) versus aerogel insulation on the LH2 tank and aluminum versus stainless steel spacer rings between the inner and outer tank. The resulting boil-off rates of LH2 provided by the one-dimensional model and three-dimensional finite element analysis (FEA) on the tank system are presented and compared to validate the results of the three-dimensional FEA. It concludes that heat flux through penetrations by conduction is as significant as that through insulation around the tank. The tank system with MLI insulation and stainless steel spacer rings result in the lowest boiling-off rate of LH2.
Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.
An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.
Membrane distillation (MD) is a membrane-based thermal desalination process capable of treating hypersaline brines. Standard MD systems rely on preheating the feed to drive the desalination process. However, relying on the feed to carry thermal energy is limited by a decline of the thermal driving force as the water moves across the membrane, and temperature polarization. In contrast, supplying heat directly into the feed channel, either through the membrane or other channel surfaces, has the potential of minimizing temperature polarization, increasing single-pass water recoveries, and decreasing the number of heat exchangers in the system. When solar thermal energy can be utilized, particularly if the solar heat is optimally delivered to enhance water evaporation and process performance, MD processes can potentially be improved in terms of energy efficiency, environmental sustainability, or operating costs. Here we describe an MD process using layered composite membranes that include a high-thermal-conductivity layer for supplying heat directly to the membrane-water interface and the flow channel. Here, the MD system showed stable performance with water flux up to 9 L/m 2 /hr, and salt rejection >99.9% over hours of desalinating hypersaline feed (100 g/L NaCl). In addition to bench-scale system, we developed a computational fluid dynamics model that successfully described the transport phenomena in the system.
The primary purpose of this project was to evaluate and improve models for the bulk-skin temperature difference to the point where they could accurately and reliably apply under a wide variety of environmental conditions. To accomplish this goal, work was conducted in three primary areas. These included production of an archive of available data sets containing measurements of the skin and bulk temperatures and associated environmental conditions, evaluation of existing skin layer models using the compiled data archive, and additional theoretical work on the development of an improved model using the data collected under diverse environmental conditions. In this work we set the basis for a new physical model of renewal type, and propose a parameterization for the temperature difference across the cool skin of the ocean in which the effects of thermal buoyancy, wind stress, and microscale breaking are all integrated by means of the appropriate renewal time scales. Ideally, we seek to obtain a model that will accurately apply under a wide variety of environmental conditions. A summary of the work in each of these areas is included in this report. A large amount of work was accomplished under the support of this grant. The grant supported the graduate studies of Sandra Castro and the preparation of her thesis which will be completed later this year. This work led to poster presentations at the 1999 American Geophysical Union Fall Meeting and 2000 IGARSS meeting. Additional work will be presented in a talk at this year's American Meteorological Society Air-Sea Interaction Meeting this May. The grant also supported Sandra Castro during a two week experiment aboard the R/P Flip (led by Dr. Andrew Jessup of the Applied Physics Laboratory) to help obtain additional shared data sets and to provide Sandra with a fundamental understanding of the physical processes needed in the models. In a related area, the funding also partially supported Dr. William Emery and Daniel Baldwin in the preparation of their publication "Accuracy of in situ sea surface temperatures used to calibrate infrared satellite measurements". The remainder of this report is drawn from these publications and presentations.
An electrically conductive spherical sample located in an electromagnetic field excited by rf (radio frequency) current in a system of coaxial coils is treated theoretically. Maxwell's equations are solved exactly and all integrals in the formulas for the fields are evaluated analytically for the case where the sphere is on the axis and the coil system is modeled by a stack of filamentary circular loops. Formulas are also derived for electromagnetic force exerted on the sphere, excess impedance in the coil system due to the presence of the sphere, and power absorbed by the sphere. All integrals in those formulas have been evaluated analytically. Force measurements are presented and they are in excellent agreement with the new theory. A low-power electromagnetic levitator that is accurately described by the theory has been demonstrated and is discussed. Experimental measurements of excess impedance are presented and compared with theory, and those results are used to demonstrate an accurate noncontact method for determining electrical conductivity. Theoretical formulas for power absorption are evaluated numerically and their usefulness in both rf heating and in making noncontact measurements of a number of thermophysical properties of materials is discussed.
The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.
Heat transfer across nuclear fuels and structural interfaces is an important factor for evaluating the performance of nuclear power systems. Specifically, heat generated as nuclear fuel fissions must be transported through the cladding material and through the reactor to reach the steam turbine for power generation. As new microreactor designs emerge, maximizing the efficiency of this heat transfer process becomes crucial to make them commercially viable. This article examines thermal diffusivity and gap conductance in uranium nitride (UN) fuel and Zircaloy-4 (Zry4) cladding using light flash analysis (LFA). Thermal diffusivity measurements were made on monolithic UN pellets and Zry4 exposed to carbon at peak operating temperatures of microreactors and show that carbon ingress has a minimal effect on thermal diffusivity when compared with identical materials not exposed to carbon. Evaluation of gap conductance at the UN-Zry4 interface was done using one-dimensional two-layer thermal transport models as a function of applied pressure. Here the results show that increasing pressure on the UN-Zry4 interface leads to gains in gap conductance per unit area in fuel-cladding assemblies at microreactor operating temperatures. While many other variables are expected to influence UN-Zry4 interfacial gap conductance (e.g. contact surface roughness, porosity, localized heating, environmental gas pressure), the work offers a demonstration of using a conventional LFA apparatus to determine this parameter at elevated temperatures.
We have developed a comprehensive model to study the dynamics and energetics of the ionosphere of Titan. We solved the one-dimensional, time-dependent, coupled continuity and momentum equations for several ion species, together with single ion and electron energy equations, in order to calculate density, velocity, and temperature profiles. Calculations were carried out for several cases corresponding to different local times and configurations of the Titan-Saturn system. In our model the effects of horizontal magnetic fields were assumed to be negligible, except for their effect on reducing the electron and ion thermal conductivities and inhibiting vertical transport in the subram region. The ionospheric density peak was found to be at an altitude of about 1100 km, in accordance with earlier model calculations. The ionosphere is chemically controlled below an altitude of about 1500 km. Above this level, ion densities differ significantly from their chemical equilibrium values due to strong upward ion velocities. Heat is deposited in a narrow region around the ionospheric peak, resulting in temperature profiles increasing sharply and reaching nearly constant values of 800-1000 deg K for electrons and 300 deg K for ions in the topside, assuming conditions appropriate for the wake region. In the subram region magnetic correction factors make the electron heat conductivities negligible, resulting in electron temperatures increasing strongly with altitude and reaching values in the order of 5000 deg K at our upper boundary located at 2200 km. Ion chemical heating is found to play an important role in shaping the ion energy balance in Titan's ionosphere.
Protective coatings can revolutionize the aerospace industry by providing a way to build a fuel-efficient and low-emission gas turbine engine operating at elevated temperature. Currently, significant challenges persist in the multi-physics modeling of damage processes in protective coatings exposed to the hot corrosive environment of gas turbine engine. Examples of such detrimental processes include an interface damage between different material layers forming protective coating, a formation of the through-the-thickness vertical cracks, a diffusion of oxygen and moisture through those cracks, an oxidation, and a delamination. In this paper, a feasibility of a multi-physics peridynamic model to provide insight to these detrimental processes has been demonstrated. 3-Dimensional peridynamic equations which are valid everywhere, including dynamically evolving discontinuities, have been presented. In conclusion, validation and demonstration have been shown for conduction heat transfer, reactive oxidation, and delamination of protective coating with a composition borrowed from the open literature.
A parametric study of the effect of bulk and interface properties on device characteristics as a function of heat and light stress under open-circuit and short-circuit conditions is developed for Cu(In,Ga)(S,Se) 2 solar cells using SCAPS-1D simulator. The variables and interdependencies modeled include: (1) conduction band offset (CBO) between buffer and absorber; (2) buffer ionized donor density; (3) CIGS shallow acceptor density; (4) CIGS deep acceptor density; (5) CIGS shallow donor density; (6) ionized acceptor concentration in the ordered vacancy compound (OVC) CIGS/buffer interface; and (7) back contact work-function. An increase in absorber shallow acceptor or a decrease in shallow donor density results in an increase in net carrier density and open-circuit voltage (V OC ). Further, a decrease in V OC is observed with an increase in the deep acceptors or neutral midgap defects in the CIGS due to higher charge carrier recombination. A change of CBO from spike to cliff condition is a dominant mechanism of decreasing VOC with an interface defect density of 10 12 cm -2 . The minimum depletion width is found to be specifically sensitive to CBO and interface defect density. An inflection in the capacitance–voltage curves (in reverse bias) is observed with simultaneous increase in bulk deep acceptor density and shallow donor density near the back contact.
High-Harmonic Fast Wave (HHFW) heating experiments in NSTX have shown that up to 60% of the injected power can be lost in the Scrape-Off Layer (SOL) when the density is above the fast wave cutoff density in front of the antenna, for which the fast wave is able to propagate into the plasma. This work models HHFW propagation in the SOL plasmas of NSTX-U using a 2D divertor SOL profile derived from the pressure assumption and the finite element temperature solution, which accounts for the high anisotropy of heat conduction in a 2D axisymmetric geometry. In this work, the two-dimensional axisymmetric SOL temperature profile is first evaluated by solving the steady-state non-linear heat conduction equation, in which thermal conductivity depends on temperature, using a finite element approach in the Petra-M workbench. A 2D density profile is then obtained from a prescribed density profile at the outer midplane assuming pressure is uniform along a flux tube in the SOL. This approach results in density and temperature profiles in which the strong asymmetric nature of heat conduction between the parallel and perpendicular background magnetic field is successfully captured. Furthermore, this work focuses on investigating the effect of the SOL plasma density profile on parasitic HHFW propagation in the SOL. The simulation results show that the radial gradient of the density profile affects the wavefield propagation in the SOL. As the density profile broadens, the wavefield intensity and its poloidal extent are reduced in the SOL, and the core coupling increases. Similarly, as shown by collisional power deposition, a proxy for power absorption, the fraction of power deposited in the SOL decreases with a broadening profile.
The proposed electron lens for Fermilab's Integrable Optics Test Accelerator (IOTA) will broaden its capabilities by enabling new research in nonlinear integrable optics, space-charge compensation, proton beam cooling, and more. The electron lens is based on a 5–10 keV, 1–2 A electron beam, shaped using a 0.7 m long, 0.8 T solenoidal magnetic field. A cryogen-free superconducting solenoid has been designed to provide this solenoidal field, taking into consideration the constraints on space, utilities, and infrastructure in the IOTA experimental hall. The solenoid is made of copper stabilized niobium-titanium conductor, conduction-cooled using 4 K closed-cycle cryocoolers. This paper describes the overall design of the solenoid encompassing its mechanical construction, current leads optimization, cryogenic thermal modeling that provides estimates of cooldown time and static/dynamic heat loads, and quench analysis.
Summary Early thermal breakthrough in enhanced geothermal systems (EGS) due to the presence of preferential flow channels is a major challenge that endangers efficient and economic heat extraction in such systems. Previous studies mainly focused on adjusting circulation rates of the working fluid, which still leaves significant amounts of untapped heat behind. Currently, there is a lack of technologies for altering flow distribution within the fracture network to achieve uniform heat sweeping in the reservoir. This work presents a novel concept for making proppants to autonomously control fracture conductivity based on the surrounding temperature. Here, proppants with negative thermal expansion coefficients have demonstrated the capability for appropriate fracture conductivity adjustment as a function of temperature to achieve uniform flow across the fracture network. Particle-particle interactions governing such functions are explicitly modeled, and then the Lattice Boltzmann methods (LBM) is used to determine the potential impact of closure stress and temperature changes on the permeability of the proposed proppant packs. Microscale analyses are further used to determine the required material properties to achieve a certain improvement in the permeability of the proppant pack. Our analyses show an enhancement in permeability and the associated fracture conductivity by half of their initial values. Field-scale analysis further confirms the effectiveness of the proposed concept as 31.4% more heat can be extracted from EGS over 50 years of production when the proposed proppants are used. Such novel proppants may effectively delay thermal breakthrough, sweep heat from larger rock volumes, and elongate the life span of the EGS project.
Knowing thermophysical and thermochemical property values of salts over the expected range of operating temperatures is essential to modeling and simulation efforts that support the commercialization of molten salt reactor (MSR) technologies. Salt properties being measured at Argonne include thermal transitions, phase behavior, heat capacity, density, surface tension, volumetric thermal expansion, thermal diffusivity, thermal conductivity, and viscosity. The properties of several salts of interest including NaCl-UCl 3 and LiF-NaF-KF (FLiNaK) have been measured at Argonne and reported previously. These property measurements are suitable for use in evaluating reactor performance during startup and the early operating life of the reactor. The physical and chemical behavior of the fuel salt is altered by the buildup of fission products and the resulting changes in heat transfer and other properties must be understood by reactor developers. This report summarizes melting behavior, heat capacity, and thermal diffusivity values that were measured for mixtures of FLiNaK and NaCl-UCl 3 salts that have been doped with surrogate fission products. Additional heat capacity measurements were also performed on a separate mixture of FLiNaK without dopants. Measured properties are suitable for incorporation into the Molten Salt Thermal Properties Database–Thermochemical (MSTDB–TC) and the Molten Salt Thermal Properties Database–Thermophysical (MSTDB–TP) to support the development of MSRs. Salts were prepared by doping eutectic mixtures of NaCl-UCl 3 and FLiNaK with surrogate fission products. These salts were prepared for use in salt spill experiments being performed at Argonne and were subsequently used for thermal property measurements. Thermal transitions and heat capacities were measured by using differential scanning calorimetry (DSC) and thermal diffusivity was measured by using laser flash analysis (LFA) at temperatures spanning the typical MSR operating range. The measured property values of the doped salts were compared to measured property values of salts without fission products to quantify the effects of the doped fission products.
A wind tunnel test program to determine aerodynamic heat transfer distributions on the McDonnell-Douglas configurations is reported. The tests were conducted at the Arnold Engineering Development Center (AEDC) in Tunnel B of the von Karman Gas Dynamics Facility (VKF). Heat-transfer rates were determined by the phase-change paint technique on 0.011-scale Stycast models using Tempilaq as the surface temperature indicator. The nominal test conditions were; Mach 8, freestream unit Reynolds numbers of 0.8 x one million, 2.5 x one million, and 3.7 x one million, and angles of attack of -5 deg, 0 deg, +5deg. Model details, test conditions, phase-change paint photographs and reduced heat-transfer coefficients are presented.
Results from wind tunnel tests to determine hypersonic aerodynamic heating rates on a NASA/Rockwell Space Shuttle Orbiter are reported. The tests were to determine Mach number effects, if any, and to obtain overall heating rate data at high Mach numbers from 10.5 to 16. The model used was a 0.0175-scale model built to Rockwell Orbiter lines VL70-000139. The model identity number is 29-0. These tests, designated OH11, were conducted in the AEDC Tunnel F.