Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase change”

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 289 records · Page 16

The effect of gravity induced free convection upon the melting phenomena of a finite paraffin slab for thermal control

Two theoretical models were developed to predict the thermal response of the phase change material to a given hot plate temperature. A two-dimensional pure conduction model was developed to predict the melting of the phase change material when heat transfer was a function of conduction. A combined conduction-convection model, also two-dimensional, was developed to predict the phase change phenomena when heat transfer was a function of conduction and gravity-induced free convection. Both models were solved using explicit finite difference approximations on a digital computer. The experimental equipment consisted of a rectangular cell utilizing a heat chamber, an expansion chamber, and a test chamber; a sixteen channel multipoint recorder, and a fluid flow system. The recorder monitered hot and cold plate temperatures and interior node temperatures at two second intervals. A comparison of theoretical temperature profiles and experimental temperature profiles is presented for six runs at various angles of inclination of the test cell with respect to the horizontal direction. A detailed discussion of results is presented.

Bain, R. L.↗

Two-dimensional model of a Space Station Freedom thermal energy storage canister

The Solar Dynamic Power Module being developed for Space Station Freedom uses a eutectic mixture of LiF-CaF2 phase change salt contained in toroidal canisters for thermal energy storage. Results are presented from heat transfer analyses of the phase-change salt containment canister. A 2-D, axisymmetric finite-difference computer program which models the canister walls, salt, void, and heat engine working fluid coolant was developed. Analyses included effects of conduction in canister walls and solid salt, conduction and free convection in liquid salt, conduction and radiation across salt vapor filled void regions, and forced convection in the heat engine working fluid. Void shape, location, and growth or shrinkage (due to density difference between the solid and liquid salt phases) were prescribed based on engineering judgement. The salt phase change process was modeled using the enthalpy method. Discussion of results focuses on the role of free-convection in the liquid salt on canister heat transfer performance. This role is shown to be important for interpreting the relationship between groundbased canister performance (in 1-g) and expected on-orbit performance (in micro-g). Attention is also focused on the influence of void heat transfer on canister wall temperature distributions. The large thermal resistance of void regions is shown to accentuate canister hot spots and temperature gradients.

Kerslake, Thomas W.↗

Phase transition studies in barium and strontium titanates at microwave frequencies

The objectives were the following: to understand the phase transformations in barium and strontium titanates as the crystals go from one temperature to the other; and to study the dielectric behavior of barium and strontium titanate crystals at a microwave frequency of 9.12 GHz and as a function of temperature. Phase transition studies in barium and strontium titanate are conducted using a cylindrical microwave resonant cavity as a probe. The cavity technique is quite successful in establishing the phase changes in these crystals. It appears that dipole relaxation plays an important role in the behavior of the dielectric response of the medium loading the cavity as phase change takes place within the sample. The method of a loaded resonant microwave cavity as applied in this work has proven to be sensitive enough to monitor small phase changes of the cavity medium.

Dahiya, Jai N.↗

Thermally Actuated Hydraulic Pumps

Thermally actuated hydraulic pumps have been proposed for diverse applications in which direct electrical or mechanical actuation is undesirable and the relative slowness of thermal actuation can be tolerated. The proposed pumps would not contain any sliding (wearing) parts in their compressors and, hence, could have long operational lifetimes. The basic principle of a pump according to the proposal is to utilize the thermal expansion and contraction of a wax or other phase-change material in contact with a hydraulic fluid in a rigid chamber. Heating the chamber and its contents from below to above the melting temperature of the phase-change material would cause the material to expand significantly, thus causing a substantial increase in hydraulic pressure and/or a substantial displacement of hydraulic fluid out of the chamber. Similarly, cooling the chamber and its contents from above to below the melting temperature of the phase-change material would cause the material to contract significantly, thus causing a substantial decrease in hydraulic pressure and/or a substantial displacement of hydraulic fluid into the chamber. The displacement of the hydraulic fluid could be used to drive a piston. The figure illustrates a simple example of a hydraulic jack driven by a thermally actuated hydraulic pump. The pump chamber would be a cylinder containing encapsulated wax pellets and containing radial fins to facilitate transfer of heat to and from the wax. The plastic encapsulation would serve as an oil/wax barrier and the remaining interior space could be filled with hydraulic oil. A filter would retain the encapsulated wax particles in the pump chamber while allowing the hydraulic oil to flow into and out of the chamber. In one important class of potential applications, thermally actuated hydraulic pumps, exploiting vertical ocean temperature gradients for heating and cooling as needed, would be used to vary hydraulic pressures to control buoyancy in undersea research vessels. Heretofore, electrically actuated hydraulic pumps have been used for this purpose. By eliminating the demand for electrical energy for pumping, the use of the thermally actuated hydraulic pumps could prolong the intervals between battery charges, thus making it possible to greatly increase the durations of undersea exploratory missions.

Jones, Jack↗

Subduction and volatile recycling in Earth's mantle

The subduction of water and other volatiles into the mantle from oceanic sediments and altered oceanic crust is the major source of volatile recycling in the mantle. Until now, the geotherms that have been used to estimate the amount of volatiles that are recycled at subduction zones have been produced using the hypothesis that the slab is rigid and undergoes no internal deformation. On the other hand, most fluid dynamical mantle flow calculations assume that the slab has no greater strength than the surrounding mantle. Both of these views are inconsistent with laboratory work on the deformation of mantle minerals at high pressures. We consider the effects of the strength of the slab using two-dimensional calculations of a slab-like thermal downwelling with an endothermic phase change. Because the rheology and composition of subducting slabs are uncertain, we consider a range of Clapeyron slopes which bound current laboratory estimates of the spinel to perovskite plus magnesiowustite phase transition and simple temperature-dependent rheologies based on an Arrhenius law diffusion mechanism. In uniform viscosity convection models, subducted material piles up above the phase change until the pile becomes gravitationally unstable and sinks into the lower mantle (the avalanche). Strong slabs moderate the 'catastrophic' effects of the instabilities seen in many constant-viscosity convection calculations; however, even in the strongest slabs we consider, there is some retardation of the slab descent due to the presence of the phase change.

King, S. D.↗

Aerodynamic heat transfer to a hypersonic research aircraft model /X-24C/ at Mach 6

The paper reports on results of heat-transfer tests conducted on a 1/29-scale model of the X-24C-12I hypersonic research aircraft configuration in a Mach 6 tunnel at a Reynolds number of thirteen million using the phase-change heat transfer technique. Sequences of phase-change heat transfer pattern photographs are presented showing windward side and leeward side heating processes. Theoretical predictions of dimensionless heat transfer coefficients along a data line on lower fuselage and on fuselage side bracket the experimental values. A turbulent heating theory gives good agreement with data when shifted to a new virtual origin.

Lawing, P. L.↗

Characteristics of Whipple Shield Performance in the Shatter Regime

Between the onset of projectile fragmentation and the assumption of rear wall failure due to an impulsive load, multi-wall ballistic limit equations are linearly interpolated to provide reasonable yet conservative predictions of perforation thresholds with conveniently simple mathematics. Although low velocity and hypervelocity regime predictions are based on analytical expressions, there is no such scientific foundation for predictions in the intermediate (or shatter) regime. As the debris flux in low earth orbit (LEO) becomes increasingly dominated by manmade pollution, the profile of micrometeoroid and orbital debris (MMOD) risk shifts continually towards lower velocities. For the International Space Station (ISS), encounter velocities below 7 km/s now constitute approximately 50% of the penetration risk. Considering that the transition velocity from shatter to hypervelocity impact regimes described by common ballistic limit equations (e.g. new non-optimum Whipple shield equation [1]) occurs at 7 km/s, 50% of station risk is now calculated based on failure limit equations with little analytical foundation. To investigate projectile and shield behavior for impact conditions leading to projectile fragmentation and melt, a series of hypervelocity impact tests have been performed on aluminum Whipple shields. In the experiments projectile diameter, bumper thickness, and shield spacing were kept constant, while rear wall thickness was adjusted to determine spallation and perforation limits at various impact velocities and angles. The results, shown in Figure 1 for normal and 45 impacts, demonstrated behavior that was not sufficiently described by the simplified linear interpolation of the NNO equation (also shown in Figure 1). Hopkins et al. [2] investigated the performance of a nominally-identical aluminum Whipple shield, identifying the effects of phase change in the shatter regime. The results (conceptually represented in Figure 2) were found to agree well with those obtained in this study at normal incidence, suggesting that shielding performance in the shatter regime could be well described by considering more complex phase conditions than currently implemented in most BLEs. Furthermore, evidence of these phase effects were found in the oblique test results, providing the basis for an empirical description of these effects that can be applied in MMOD risk assessment software. In this paper, results of the impact experiments are presented, and characteristics of target damage are evaluated. A comparison of intermediate velocity impact failure mechanisms in current BLEs are discussed and compared to the findings of the experimental study. Risk assessment calculations have been made on a simplified structure using currently implemented penetration equations and predicted limits from the experimental program, and the variation in perceived mission risk is discussed. It was found that ballistic limit curves that explicitly incorporated phase change effects within the intermediate regime lead to a decrease in predicted MMOD risk for ISS-representative orbits. When considered for all Whipple-based shielding configurations onboard the ISS, intermediate phase change effects could lead to significant variations in predicted mission risk.

Ryan, Shannon↗

Selection of high temperature thermal energy storage materials for advanced solar dynamic space power systems

Under the direction of NASA's Office of Aeronautics and Technology (OAST), the NASA Lewis Research Center has initiated an in-house thermal energy storage program to identify combinations of phase change thermal energy storage media for use with a Brayton and Stirling Advanced Solar Dynamic (ASD) space power system operating between 1070 and 1400 K. A study has been initiated to determine suitable combinations of thermal energy storage (TES) phase change materials (PCM) that result in the smallest and lightest weight ASD power system possible. To date the heats of fusion of several fluoride salt mixtures with melting points greater than 1025 K have been verified experimentally. The study has indicated that these salt systems produce large ASD systems because of their inherent low thermal conductivity and low density. It is desirable to have PCMs with high densities and high thermal conductivities. Therefore, alternate phase change materials based on metallic alloy systems are also being considered as possible TES candidates for future ASD space power systems.

Lacy, Dovie E.↗

Selection of high temperature thermal energy storage materials for advanced solar dynamic space power systems

Under the direction of NASA's Office of Aeronautics and Technology (OAST), the NASA Lewis Research Center has initiated an in-house thermal energy storage program to identify combinations of phase change thermal energy storage media for use with a Brayton and Stirling Advanced Solar Dynamic (ASD) space power system operating between 1070 and 1400 K. A study has been initiated to determine suitable combinations of thermal energy storage (TES) phase change materials (PCM) that result in the smallest and lightest weight ASD power system possible. To date the heats of fusion of several fluoride salt mixtures with melting points greater than 1025 K have been verified experimentally. The study has indicated that these salt systems produce large ASD systems because of their inherent low thermal conductivity and low density. It is desirable to have PCMs with high densities and high thermal conductivities. Therefore, alternate phase change materials based on metallic alloy systems are also being considered as possible TES candidates for future ASD space power systems.

Lacy, Dovie E.↗

Phase Time and Envelope Time in Time-Distance Analysis and Acoustic Imaging

Time-distance analysis and acoustic imaging are two related techniques to probe the local properties of solar interior. In this study, we discuss the relation of phase time and envelope time between the two techniques. The location of the envelope peak of the cross correlation function in time-distance analysis is identified as the travel time of the wave packet formed by modes with the same w/l. The phase time of the cross correlation function provides information of the phase change accumulated along the wave path, including the phase change at the boundaries of the mode cavity. The acoustic signals constructed with the technique of acoustic imaging contain both phase and intensity information. The phase of constructed signals can be studied by computing the cross correlation function between time series constructed with ingoing and outgoing waves. In this study, we use the data taken with the Taiwan Oscillation Network (TON) instrument and the Michelson Doppler Imager (MDI) instrument. The analysis is carried out for the quiet Sun. We use the relation of envelope time versus distance measured in time-distance analyses to construct the acoustic signals in acoustic imaging analyses. The phase time of the cross correlation function of constructed ingoing and outgoing time series is twice the difference between the phase time and envelope time in time-distance analyses as predicted. The envelope peak of the cross correlation function between constructed ingoing and outgoing time series is located at zero time as predicted for results of one-bounce at 3 mHz for all four data sets and two-bounce at 3 mHz for two TON data sets. But it is different from zero for other cases. The cause of the deviation of the envelope peak from zero is not known.

Chou, Dean-Yi↗

Differential scanning calorimetric survey of brominated PAN, pitch-based and vapor-grown fibers

A differential scanning calorimetric survey of brominated PAN, pitch-based and vapor-grown fibers revealed that only vapor-grown fibers heat-treated to temperatures of 2250 C or above exhibited the phase change at 100 C associated with the two-dimensional melting of bromine, similar to that observed in HOPG. All the pitch-based fibers tested, as well as vapor-grown fibers annealed to 2250 C, showed a substantial decrease in resistivity upon bromination but gave no indication of a phase change, suggesting that the bromine in these fibers does not have the long-range order needed to undergo two-dimensional melting. The brominated PAN fibers and unannealed vapor-grown fibers showed no phase change and no change in resistivity, indicating that these fibers do not intercalate bromine.

Jaworske, D. A.↗

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗

Ice Particle Analysis of the Honeywell AL502 Engine Booster

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell ALF502 engine. The analysis focused on two closely related conditions one of which produced an icing event and another which did not during testing of the ALF502 engine in the Propulsion Systems Lab (PSL) at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.63 ice accretion software. The inflow conditions for the two conditions were similar with the main differences being that the condition that produced the icing event was 6.8 K colder than the non-icing event case and the inflow ice water content (IWC) for the non-icing event case was 50% less than for the icing event case. The particle analysis, which considered sublimation, evaporation and phase change, was generated for a 5 micron ice particle with a sticky impact model and for a 24 micron median volume diameter (MVD), 7 bin ice particle distribution with a supercooled large droplet (SLD) splash model used to simulate ice particle breakup. The particle analysis did not consider the effect of the runback and re-impingement of water resulting from the heated spinner and anti-icing system. The results from the analysis showed that the amount of impingement for the components were similar for the same particle size and impact model for the icing and non-icing event conditions. This was attributed to the similar aerodynamic conditions in the booster for the two cases. The particle temperature and melt fraction were higher at the same location and particle size for the non-icing event than for the icing event case due to the higher incoming inflow temperature for the non-event case. The 5 micron ice particle case produced higher impact temperatures and higher melt fractions on the components downstream of the fan than the 24 micron MVD case because the average particle size generated by the particle breakup was larger than 5 microns which yielded less warming and melting. The analysis also showed that the melt fraction and wet bulb temperature icing criterion developed during tests in the Research Altitude Test Facility (RATFac) at the National Research Council (NRC) of Canada were useful in predicting icing events in the ALF502 engine. The development of an ice particle impact model which includes the effects of particle breakup, phase change, and surface state is necessary to further improve the prediction of ice particle transport with phase change through turbomachinery.

Ice Accretion Codes↗

An implicit formulation for the one-dimensional two-phase multi-interface Stefan Problem

In recent years there has been an intensified effort in the analysis of the Stefan Problem. Particular interest is directed towards determining the location of the interface between the solid and liquid phases of the substance being analyzed. A method has been developed using an implicit finite difference scheme which provides an accurate calculation of the phase change interface locations for a problem with simultaneous multiple phases. This is achieved through an iterative technique which requires the phase change energy to be correctly distributed within the nodes which are undergoing a phase transition. The new method predicts the interface locations without requiring time step size control or node size modification at the interfaces.

Williams, S. D.↗

Thermophysical Property Measurements of Indium Iodide Crystals

InI single crystals are a promising room temperature detector material for X-rays and γ-rays. To improve crystal growth of the material by simulations, knowledge of thermophysical properties is essential, and since InI is orthorhombic, the anisotropy has to be taken into account. The temperature dependence of several thermophysical properties have been measured for InI, including the anisotropic thermal expansion, specific heat, and the thermal diffusivity in the b direction. The anisotropic thermal expansion coefficients, determined by X-ray diffraction, were α11=1.03·10-5 K-1, α22= 3.77·10-5 K-1, and α33 = 6.26·10-5 K-1. The specific heat, measured by DSC, was 0.226 J g-1 K-1 at 335K, with a temperature dependence of 9.582·10-5 J g-1 K-2. In the course of the X-ray diffraction and DSC measurements, it was shown that supposed phase changes, reported in older literature, are actually not phase changes but oxidation effects. The thermal diffusivity in the b direction, measured by the Xenon Flash method, was 0.288·10-6 m2s-1 at RT, decreasing to 0.253·10-6 m2 s-1 at 450K. In addition, the volume increase upon melting and the thermal expansion of the melt have been determined.

Martin Volz↗

A Hydraulic Motor-Alternator System for Ocean-Submersible Vehicles

An ocean-submersible vehicle has been developed at JPL that moves back and forth between sea level and a depth of a few hundred meters. A liquid volumetric change at a pressure of 70 bars is created by means of thermal phase change. During vehicle ascent, the phase-change material (PCM) is melted by the circulation of warm water and thus pressure is increased. During vehicle descent, the PCM is cooled resulting in reduced pressure. This pressure change is used to generate electric power by means of a hydraulic pump that drives a permanent magnet (PM) alternator. The output energy of the alternator is stored in a rechargeable battery that powers an on-board computer, instrumentation and other peripherals.The focus of this paper is the performance evaluation of a specific hydraulic motor-alternator system. Experimental and theoretical efficiency data of the hydraulic motor and the alternator are presented. The results are used to evaluate the optimization of the hydraulic motor-alternator system. The integrated submersible vehicle was successfully operated in the Pacific Ocean near Hawaii. A brief overview of the actual test results is presented.

alternator↗

The seasonal variation of the D region as inferred from propagation characteristics of LF radio waves

The propagation data of JG2AS 40 kHz (Japanese Standard Frequency), Loran C 100 kHz radio waves, and meteorological data were analyzed to study the association of propagation characteristics of LF radio waves with the atmospheric circulation in the mesosphere. The monthly averaged electric fields were depicted on the complex plane for typical summer and winter months, June and November. The locus traced out by the electric field vector during daytime is nearly circular. This is because during daytime the amplitude of the sky wave remains nearly constant while its phase changes in accord with the height change of the reflection layer, and thus the electric field vector traces out a circular locus with its center at the tip of the supposed ground wave vector. The locus has a loop during the sunrise or sunset period, which seems to arise from interference of two waves reflected by two different layers. In June the amplitude of the sky wave decreases rapidly before the dawn or increases after the dusk. In November such rapid change is not observed. During nighttime, the sky wave phase changes in such a way as to suggest that the reflection height moves upwards with time before midnight or lowers after midnight in November. In June it changes similarly before midnight, but after midnight it varies erratically. These characteristics are closely related to the structure of the D region, which is clearly shown by simulating the loci traced out by electric fields.

Ishimine, T.↗