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 361 records · Page 20

On the preferential acceleration and heating of solar wind heavy ions

It is pointed out that the properties of solar wind ions heavier than H(+) have presented a puzzle to theorists for some time. The most recent theoretical efforts to explain the observations are considered, taking into account studies conducted by Dusenbery and Hollweg (1981) and Marsch et al. (1982). The first model of a wave-driven, three-fluid, supersonic solar wind is considered. Details of the resonant wave-particle interaction are discussed and the relevant equations for the model are derived. An investigation is conducted of the effects of resonant cyclotron acceleration and heating by a spectrum of nondispersive ion cyclotron waves on the model. It is found that there are serious difficulties with the hypothesis that the resonant cyclotron interaction is responsible for the preferential heating and acceleration of heavy ions in the solar wind. The model presented here is unable to produce the required energization by this mechanism. However, this model should be viewed as an approximate one, subject to possible correction by more advanced calculations.

Isenberg, P. A.↗

Dynamical System Scaling Application to Zircaloy Cladding Thermal Response During Reactivity-Initiated Accident Experiment

New fuel design and development currently requires 20 to 25 years to be qualified for use by the nuclear power industry. The thermal-hydraulics community has taken advantage of scaling theory to design reduced scale experiments that correctly preserve dominant key phenomena while quantifying distorted phenomena. These techniques can be leveraged in the design and analysis of fuel performance experiments to help reduce the timeline associated with fuel design and development. This study uses the Dynamical System Scaling (DSS) method to analyze cladding temperature data from the recent SETH-C experiment in the TREAT facility and accompanying BISON simulations to assess dynamic distortions occurring throughout the fast power excursion transient. The DSS analysis revealed that on the cool down from peak cladding temperature that the fuel radial power profile is the most sensitive modeling parameter with a heterogenous radial peaking factor corresponding to the lowest distortion compared to a uniform energy deposition. For the heat up to peak cladding temperature the heterogeneous radial power profile corresponded to the shortest process action. Finally, for the heat up to peak cladding temperature, the gap conductance model sensitivity was quantified using process action and shows that the default Light Water Reactor gap conductance model corresponded to the longest process action.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Low-order aquifer thermal energy storage model for geothermal system simulation

This paper presents a low-order aquifer thermal energy storage (ATES) model for simulation of combined subsurface and above-surface energy systems. The model is included in the Modelica IBPSA Library, which is a free open-source library with basic models for building and district energy and control systems. The model uses a lumped-component method, in which the transient conductive-convective heat and mass transfer equation is radially discretized. To verify the accuracy of the model, we present an inter-model comparison from a simulation test suite. Results show that the Modelica ATES model is in good agreement, with a normalized mean bias error for yearly variation of aquifer temperatures of 1.6×10−2 and 9×10−5 at 1 m and 10 m distance from the well.

Maccarini, Alessandro↗

NO PLIF Visualizations of the Orion Capsule in LENS-I

Planar laser-induced fluorescence (PLIF) of nitric oxide (NO) was used to visualize the interaction of reaction-control-system (RCS) jet flows in the wake of a hypersonic capsule reentry vehicle. The tests were performed at the Calspan University at Buffalo Research Center's (CUBRC) LENS-I reflected shock tunnel facility. This was the first application of PLIF to study RCS jets in a large-scale pulsed hypersonic facility. The LENS-I facility allowed RCS jet flows to be studied while varying the flow enthalpy, Reynolds number, angle of attack and jet configuration. The interaction of pitch and roll jets with the flowfield was investigated. Additionally, thin film sensors were used to monitor heat transfer on the surface of the model to detect any localized heating resulting from the firing of the RCS jets. Tests were conducted with the model held at angles of attack of 18deg and 22deg. The nominal Mach number in all tests was 8, while Reynolds number based on model diameter ranged from 2.2x10(exp 6) - 1.5x10(exp 7). Images were processed using the Virtual Diagnostics Interface (ViDI) system developed at NASA Langley Research Center to provide a three-dimensional display of the experimental data.

Combs, C.↗

Numerical modeling of laser thermal propulsion flows

An review of the problems associated with modeling laser thermal propulsion flows, a synopsis of the status of such models, and the attributes of a successful model are presented. The continuous gaseous hydrogen laser-supported combustion wave (LSCW) thruster, for which a high-energy laser system (preferably space-based) should exist by the time the propulsion technology is developed, is considered in particular. The model proposed by Raizer (1970) is based on the assumptions of one-dimensional flow at constant pressure with heat conduction as the principal heat transfer mechanism. Consideration is given to subsequent models which account for radiative transfer into the ambient gas; provide a two-dimensional generalization of Raizer's analysis for the subsonic propagation of laser sparks in air; include the effect of forward plasma radiation in a one-dimensional model; and attempt a time-dependent (elliptic) solution of the full Navier-Stokes equations for the flow in a simple axisymmetric thruster. Attention is also given to thruster and nozzle flow models and thermodynamic and transport properties.

Mccay, T. D.↗

Verification of thermal analysis codes for modeling solid rocket nozzles

One of the objectives of the Solid Propulsion Integrity Program (SPIP) at Marshall Space Flight Center (MSFC) is development of thermal analysis codes capable of accurately predicting the temperature field, pore pressure field and the surface recession experienced by decomposing polymers which are used as thermal barriers in solid rocket nozzles. The objective of this study is to provide means for verifications of thermal analysis codes developed for modeling of flow and heat transfer in solid rocket nozzles. In order to meet the stated objective, a test facility was designed and constructed for measurement of the transient temperature field in a sample composite subjected to a constant heat flux boundary condition. The heating was provided via a steel thin-foil with a thickness of 0.025 mm. The designed electrical circuit can provide a heating rate of 1800 W. The heater was sandwiched between two identical samples, and thus ensure equal power distribution between them. The samples were fitted with Type K thermocouples, and the exact location of the thermocouples were determined via X-rays. The experiments were modeled via a one-dimensional code (UT1D) as a conduction and phase change heat transfer process. Since the pyrolysis gas flow was in the direction normal to the heat flow, the numerical model could not account for the convection cooling effect of the pyrolysis gas flow. Therefore, the predicted values in the decomposition zone are considered to be an upper estimate of the temperature. From the analysis of the experimental and the numerical results the following are concluded: (1) The virgin and char specific heat data for FM 5055 as reported by SoRI can not be used to obtain any reasonable agreement between the measured temperatures and the predictions. However, use of virgin and char specific heat data given in Acurex report produced good agreement for most of the measured temperatures. (2) Constant heat flux heating process can produce a much higher heating rate than the radiative heating process. The results show that heating rates of 125 C/s (225 F/s) can be achieved. (3) The electrical resistance of the FM 5055 samples were about 150 Omega in the virgin state, and decreased to about 35 Omega when charred. A reliable scheme must be developed to electrically insulate the composite from the heater foil in order to prevent any current leakage into the sample which can result in volumetric heating of the char zone of the composite.

Keyhani, M.↗

Magnetization around mix jets entering inertial confinement fusion fuel

Engineering features are known to cause jets of ablator materials to enter the fuel hot-spot in inertial confinement fusion implosions. The Biermann battery mechanism wraps them in a self-generated magnetic field. Here we show that higher-Z jets have an additional thermoelectric magnetic source term that is not present for hydrogen jets, verified here through a kinetic simulation. It has similar magnitude to the Biermann term. We then include this in an extended magneto-hydrodynamic approach to post process an xRAGE radiation-hydrodynamic implosion simulation. The simulation includes an accurate model for the capsule fill tube, producing a dense carbon jet that becomes wrapped in a 4000 T magnetic field. A simple spherical carbon mix model shows that this insulates the electron heat conduction enough to cause contraction of the jet to an optically thick equilibrium. The denser magnetized jet hydrodynamics could change its core penetration and therefore the final mix mass, which is known to be well correlated with fusion yield degradation. Fully exploring this will require self-consistent magneto-hydrodynamic simulations. Experimental signatures of this self-magnetization may emerge in the high energy neutron spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗

Cryogenic System Preliminary Design for a 0.5m-Long, Conduction-Cooled Nb3Sn Undulator Magnet Prototype

There are several NbTi superconducting undulator (SCU) magnets currently in operation at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). The development on Nb3Sn-based superconducting undulator magnets at APS is underway due to the potential to further enhance the performance of the SCUs. Superconducting undulator magnets need to keep temperature gradients minimized in order to retain thermal and operating current margin. We have designed the thermal links for efficient heat conduction using 3D finite element analysis (FEA) simulation in COMSOL Multiphysics software, which was later used for the calculation of the temperature distribution across a 0.5 m long, conduction-cooled, Nb3Sn undulator magnet prototype that includes both conductive heat transfer and radiative heating components. We have modelled the evolution of the thermal properties of the magnet winding as well as other cold parts during cool-down from ambient temperature with the operation of a SHI RDK-415D cryocooler to examine the estimated time that is needed for the cooldown as well as the baseline temperature we could achieve. A key result was that a maximum coil delta T = 0.11 K along the designed thermal links and a temperature range of 3.42 K-3.53 K for the winding were predicted at steady state. A pair of G10 support rods for carrying the magnet as well as a pair of current leads made of metals were also designed and optimized for the testing system, targeting for an operation temperature of 4 K and a coil current of 950 A DC.

Zhang, Danlu↗

Models of the Venus ionosphere

Venus daytime ionospheric models using electron, ion and neutral gas heat conduction with momentum and chemical equations for charged particle densities

Bauer, S. J.↗

Heat transfer in rotating serpentine passages with trips normal to the flow

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multipass, heat transfer model with both radially inward and outward flow. Trip strips on the leading and trailing surfaces of the radial coolant passages were used to produce the rough walls. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature ratio, Rossby number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from stationary and rotating similar models with trip strips. The heat transfer coefficients on surfaces, where the heat increased with rotation and buoyancy, varied by as much as a factor of four. Maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels obtained with the smooth wall model. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation, varied by as much as a factor of three due to rotation and buoyancy. It was concluded that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips and that the effects of rotation were markedly different depending upon the flow direction.

Wagner, J. H.↗

Heat transfer in rotating serpentine passages with trips normal to the flow

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multipass, heat transfer model with both radially inward and outward flow. Trip strips on the leading and trailing surfaces of the radial coolant passages were used to produce the rough walls. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature ratio, Rossby number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from stationary and rotating similar models with trip strips. The heat transfer coefficients on surfaces, where the heat increased with rotation and buoyancy, varied by as much as a factor of four. Maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels obtained with the smooth wall model. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation, varied by as much as a factor of three due to rotation and buoyancy. It was concluded that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips and that the effects of rotation were markedly different depending upon the flow direction.

Wagner, J. H.↗

Preliminary SP-100/Stirling heat exchanger designs

Analytic modeling of several heat exchanger concepts to couple the SP-100 nuclear reactor primary lithium loop and the Space Stirling Power Convertor (SSPC) was performed. Four 25 kWe SSPC's are used to produce the required 100 kW of electrical power. This design work focused on the interface between a single SSPC and the primary lithium loop. Manifolding to separate and collect the four channel flow was not modeled. This work modeled two separate types of heat exchanger interfaces (conductive coupling and radiative coupling) to explore their relative advantages and disadvantages. The minimum mass design of the conductively coupled concepts was 18 kg or 0.73 kg/kWe for a single 25 kWe convertor. The minimum mass radiatively coupled concept was 41 kg or 1.64 kg/kWe. The direct conduction heat exchanger provides a lighter weight system because of its ability to operate the Stirling convertor evaporator at higher heat fluxes than those attainable by the radiatively coupled systems. Additionally the conductively coupled concepts had relatively small volumes and provide potentially simpler assembly. Their disadvantages were the tight tolerances and material joining problems associated with this refractory to superalloy interface. The advantages of the radiatively coupled designs were the minimal material interface problems.

Schmitz, Paul↗

CISCC Countermeasures for Spent Fuel Canisters in Concrete Casks

Aging management of spent fuel is an important issue worldwide. Chloride-induced stress corrosion cracking (CISCC) of welded canisters is particularly important for concrete casks, where it occurs in response to a combination of environmental, material, and residual stress factors. In this paper, we introduce a countermeasure to prevent the occurrence of CISCC by adjusting the flow rate of the cooling air into the concrete cask so that the canister surface temperature does not decrease below 100°C, which is a conservative deliquesce temperature threshold for CISCC occurrence. The flow adjustment device was installed in a cask model, and two tests were conducted under the heat rate conditions corresponding to those of an actual canister. In the first test, the device was operated under a condition where the canister surface temperature was below 100°C to verify whether it could increase the canister surface temperature to above 100°C. In the second test, the heat rate was decreased from a state where the canister surface temperature was above 100°C, to verify whether the device could maintain the canister surface temperature above 100°C. Good results were obtained in both cases. However, in the case of a low decay heat rate of spent fuel after extended dry storage, the temperature of the canister surface may decrease below 100°C even if no air supply is provided. Therefore, in this case, it may be necessary to take measures such as blocking air outlets.

Takeda, Hirofumi↗

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

Community Geothermal: Soil Conductivity, Borehole Design, Energy Models, and Load Data for a Residential System Development - Hinesburg, VT

This dataset contains materials from the Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES) project, which evaluated the techno-economic feasibility of a community geothermal system for a residential development in Hinesburg, VT. The dataset includes detailed soil conductivity test reports, energy models, borehole design reports, hourly energy loads for heating, cooling, and hot water, and design layouts. EnergyPlus was used to model building energy loads, and Modelica software was applied for geothermal loop sizing based on these loads and soil conductivity results. Python scripts for network design further refined the models. Key files include PDF reports on borehole design (with projections for 1-year, 15-year, and 30-year systems), soil conductivity test results, EnergyPlus modeling outputs, and 2D/3D design drawings in PDF, DWG, and DXF formats. Python notebooks for network design and OnePipe model files are also provided, with Modelica required for viewing certain files. Outputs and modeling data are in various formats including CSV, JPG, HTML, and IDF, with units and data clearly labeled to support understanding of system design and performance for the proposed geothermal solution.

15 GEOTHERMAL ENERGY↗

Thermal transport in defective and disordered materials

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

42 ENGINEERING↗

Heat transfer tests of an 0.006-scale thin-skin space shuttle thermocouple model (41-OTS) in the Langley Research Center unitary plan wind tunnel at M equals 3.7 (IH16)

The results are presented of supersonic heat transfer tests performed on the .006 scale space shuttle vehicle model (41-OTS) in the Langley Research Center Unitary Plan Wind Tunnel. These tests were conducted to parametrically investigate ascent heating of the integrated vehicle and its components. The tests were conducted at a nominal Mach number of 3.7 and Reynolds numbers per foot of 2 and 5 million. The model configurations investigated were the integrated vehicle and each component alone (i.e. orbiter, tank and SRB). All the configurations were run with and without transition strips and through an angle of attack range of 0 deg to minus 5 deg with the exception of the SRB which was tested through an angle of attack range of minus 5 deg to 90 deg. The heat transfer data were obtained from 223 iron constantan thermocouples attached to stainless steel thin-skin areas of the model.

Walstad, D. G.↗