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 757 records · Page 42

Photometric and spectroscopic observations of 5145 Pholus

Light curve observations carried out in January and February 1992 are presented. We find an amplitude of 0.2 magnitude and a period of 0.4157 day. Fine structures of the light curve and their variability with time are discussed in terms of the changing phase angle and deviations of the object's shape from a regular ellipsoid. The present orbital and physical properties alone for 5145 Pholus do not provide a clue for either a cometary or an asteroidal origin. It seems to have properties in common with both and thus could form a bridge between these two groups.

Hoffmann, Martin↗

ASDA - Advanced Suit Design Analyzer computer program

An ASDA model developed to evaluate the heat and mass transfer characteristics of advanced pressurized suit design concepts for low pressure or vacuum planetary applications is presented. The model is based on a generalized 3-layer suit that uses the Systems Integrated Numerical Differencing Analyzer '85 in conjunction with a 41-node FORTRAN routine. The latter simulates the transient heat transfer and respiratory processes of a human body in a suited environment. The user options for the suit encompass a liquid cooled garment, a removable jacket, a CO2/H2O permeable layer, and a phase change layer.

Bue, Grant C.↗

A finite element analysis of the freeze/thaw behavior of external artery heat pipes

A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.

Lu, X. J.↗

Incineration for resource recovery in a closed ecological life support system

A functional schematic, including mass and energy balance, of a solid waste processing system for a controlled ecological life support system (CELSS) was developed using Aspen Plus, a commercial computer simulation program. The primary processor in this system is an incinerator for oxidizing organic wastes. The major products derived from the incinerator are carbon dioxide and water, which can be recycled to a crop growth chamber (CGC) for food production. The majority of soluble inorganics are extracted or leached from the inedible biomass before they reach the incinerator, so that they can be returned directly to the CGC and reused as nutrients. The heat derived from combustion of organic compounds in the incinerator was used for phase-change water purification. The waste streams treated by the incinerator system conceptualized in this work are inedible biomass from a CGC, human urine (including urinal flush water) and feces, humidity condensate, shower water, and trash. It is estimated that the theoretical minimum surface area required for the radiator to reject the unusable heat output from this system would be 0.72 sq m/person at 298 K.

Upadhye, R. S.↗

Rates of fluvio-thermal erosion on Mars

In order to get some ideas about the order of magnitude of fluvio-thermal recession rate, a mathematical thermal model for Martian outflow channels is proposed. This model corresponds to a system undergoing a permanent thermal regime where the surface temperature is constant and equal to the phase change temperature (due to the immediate removal of melted materials). This is an ablation model. For its application, estimations of the heat transfer coefficient h and thermal flux q are necessary. Determination of these coefficients needs the calculation of dimensionless numbers (Reynolds, Prandtl, and Nusselt), and the consideration of turbulent regime of the flow.

Aguirre-Puente, J.↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

A new release rate equation to model the phase change of water to steam in composite materials was derived from the theory of molecular diffusion and equilibrium moisture concentration. The new model is dependent on internal pressure, the microstructure of the voids and channels in the composite materials, and the diffusion properties of the matrix material. Hence, it is more fundamental and accurate than the empirical Arrhenius rate equation currently in use. The model was mathematically formalized and integrated into the thermostructural analysis code CHAR. Parametric studies on variation of several parameters have been done. Comparisons to Arrhenius and straight-line models show that the new model produces physically realistic results under all conditions.

Mcmanus, Hugh L.↗

Dynamics of face seals for high speed turbomachinery

Face seals in rocket engine fuel and oxidizer turbopumps have been the subject of intense investigation for over 25 years. While advances have been made in the understanding of thin film lubrication between seal faces, valuable data has been produced on the friction and wear of material pairs in cryogenic environments; pioneering work has been done on the effect of lubricant phase change in seals, and many improvements have been made in mechanical seal design. Relatively superficial attention has been given to the vibrational dynamics of face seals in high-speed turbomachinery. BHR Group Ltd. (formerly BHRA) has recently completed the first stage of a study, commissioned by the European Space Agency, to investigate this area. This has involved the development of a two-dimensional adiabatic, turbulent lubrication model for thick gas film applications, the production of an integrated mathematical model of gas seal vibrational dynamics for thin film applications, implementation in software, the undertaking of an experimental program to validate software against variations in operating conditions and design variables, and suggestions for improved seal design.

Leefe, Simon↗

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

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.↗

Thermocapillary flow with evaporation and condensation and its effect on liquid retention in low-G fluid acquisition devices

The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.

Schmidt, George R.↗

Personal cooling systems: Possibilities and limitations

Personal thermal control by means of gas- or liquid-conditioned garments was developed during the 1960s and has been applied in a variety of aerospace and industrial settings. Both USAF fighter pilots and astronauts are required to wear heavy protective clothing which insulates them from the environment and thus creates stress through storage of metabolic heat. The problem is particularly severe in astronauts who perform heavy physical work during extra-vehicular activity (EVA); without artificial cooling they could reach incapacitating hyperthermia in a matter of minutes. This paper reviews the factors which influence the design of personal cooling systems. An important early step is determination of acceptable heat stress level, taking into account possible interactions of heat stress with other physiological problems such as motion sickness, diminished plasma volume, decompression sickness and acceleration tolerance. Other factors which require consideration include the work schedule, the area to be covered by the cooling garment, and the practicalities of a fixed or body-mounted heat sink and its power source. Nearly every imaginable heat sink has been proposed or tried over the past 30 years, including direct gas systems, phase-change systems with open or closed loops and thermoelectric heat sinks. The latter are now the system of choice for aircraft.

Nunneley, Sarah A.↗

Effect of Microgravity on Material Undergoing Melting and Freezing: the TES Experiment

This experiment is the first to melt and freeze a high temperature thermal energy storage (TES) material under an extended duration of microgravity. It is one of a series to validate an analytical computer program that predicts void behavior of substances undergoing phase change under microgravity. Two flight experiments were launched in STS-62. The first, TES-1, containing lithium fluoride in an annular volume, performed flawlessly in the 22 hours of its operation. Results are reported in this paper. A software failure in TES-2 caused its shutdown after 4 seconds. A computer program, TESSIM, for thermal energy storage simulation is being developed to analyze the phenomena occurring within the TES containment vessel. The first order effects, particularly the surface tension forces, have been incorporated into TESSIM. TESSIM validation is based on two types of results. First is the temperature history of various points of the containment structure, and second, upon return from flight, the distribution of the TES material within the containment vessel following the last freeze cycle. The temperature data over the four cycles showed a repetition of results over the third and fourth cycles. This result is a confirmation that any initial conditions prior to the first cycle had been damped out by the third cycle. The TESSIM simulation showed a close comparison with the flight data. The solidified TES material distribution within the containment vessel was obtained by a tomography imaging process. The frozen material was concentrated toward the colder end of the annular volume. The TESSIM prediction showed the same pattern. With the general agreement of TESSIM and the data, a computerized visual representation can be shown which accurately shows the movement and behavior of the void during the entire freezing and melting cycles.

Namkoong, David↗

Nucleation and growth of crystals under cirrus and polar stratospheric cloud conditions

Laboratory studies examine phase changes of hygroscopic substances which occur as aerosol in stratosphere and troposphere (sodium chloride, ammonium sulfate, ammonium bisulfate, nitric acid, sulfuric acid), under controlled conditions, in samples volume 1 to 10(exp -4) ml. Crystallization of salts from supersaturated solutions is examined by slowly evaporating a solution drop on a substrate, under controlled relative humidity, until self nucleation occurs; controlled nucleation of ice in a mm capillary U-tube gives a measured ice crystallization velocity at known supercooling. Two states of crystallization occur for regions where hydrates exist. It is inferred that all of the materials readily exist as supersaturated/supercooled solutions; the degree of metastability appears to be slightly enhanced by inclusion of aircraft produced soot. The crystallization velocity is taken as a measure of viscosity. Results suggest an approach to a glass transition at high molality, supersaturation and/or supercooling within the range of atmospheric interest. It is hypothesized that surface reactions occur more readily on solidified particles - either crystalline or glass, whereas volume reactions are more important on droplets with sufficiently low viscosity and volume diffusivity. Implications are examined for optical properties of such particles in the atmosphere. In a separate experiment, crystal growth was examined in a modified thermal vapor diffusion chamber over the range of cirrus temperature (-30 to -70 C) and under controlled supersaturation and air pressure. The crystals grew at a velocity of 1-2 microns/s, thickness 60-70 micron, in the form of thin column crystals. Design criteria are given for a system to investigate particle growth down to -100 C, (PSC temperatures) where nitric acid particles can be grown under similar control and in the form of hydrate crystals.

Hallett, John↗

Numerical modeling of Bridgman growth of PbSnTe in a magnetic field

In this work we study heat and mass transport, fluid motion, and solid/liquid phase change in the process of steady Bridgman growth of Pb(.8)Sn(.2)Te (LTT) in an axially-imposed uniform magnetic field under terrestrial and microgravity conditions. In particular, this research is concerned with the interrelationships among segregation, buoyancy-driven convection, and magnetic damping in the LTT melt. The main objectives are to provide a quantitative understanding of the complex transport phenomena during solidification of the nondilute binary of LTT, to provide estimates of the strength of magnetic field required to achieve the desired diffusion-dominated growth, and to assess the role of magnetic damping for space and earth based control of the buoyancy-induced convection. The problem was solved by using FIDAP and numerical results for both vertical and horizontal growth configurations with respect to the acceleration of gravity vector are presented.

Yao, Minwu↗

Numerical simulation of heat and mass transport during space crystal growth with MEPHISTO

The MEPHISTO space experiments are collaborative United States and French investigations aimed at understanding the fundamentals of crystal growth. Microgravity experiments were conducted aboard the USMP-1 and -2 missions on STS-52 and 62 in October 1992 and March 1994 respectively. MEPHISTO is a French designed and built Bridgman type furnace which uses the Seebeck technique to monitor the solid/liquid interface temperature and Peltier pulsing to mark the location and shape of the solid/liquid interface. In this paper the Bridgman growth of Sn-Bi and Bi-Sn under terrestrial and microgravity conditions is modeled using the finite element code, FIDAP*. The numerical model considers fully coupled heat and mass transport, fluid motion and solid/liquid phase changes in the crystal growth process. The primary goals of this work are: to provide a quantitative study of the thermal buoyancy-induced convection in the melt for the two flight experiments; to compare the vertical and horizontal growth configurations and systematically evaluate the effects of various gravity levels on the solute segregation. Numerical results of the vertical and horizontal Bridgman growth configurations are presented.

Yao, Minwu↗

Microphysical modeling of cirrus. 2: Sensitivity studies

The one-dimensional cirrus model described in part 1 of this issue has been used to study the sensitivity of simulated cirrus microphysical and radiative properties to poorly known model parameters, poorly understood physical processes, and environmental conditions. Model parameters and physical processes investigated include nucleation rate, mode of nucleation (e.g., homogeneous freezing of aerosols and liquid droplets or heterogeneous deposition), ice crystal shape, and coagulation. These studies suggest that the leading sources of uncertainty in the model are the phase change (liquid-solid) energy barrier and the ice-water surface energy which dominate the homogeneous freezing nucleation rate and the coagulation sticking efficiency at low temperatures which controls the production of large ice crystals (radii greater than 100 mcirons). Environmental conditions considered in sensitivity tests were CN size distribution, vertical wind speed, and cloud height. We found that (unlike stratus clouds) variations in the total number of condensation nuclei (NC) have little effect on cirrus microphysical and radiative properties, since nucleation occurs only on the largest CN at the tail of the size distribution. The total number of ice crystals which nucleate has little or no relationship to the number of CN present and depends primarily on the temperature and the cooling rate. Stronger updrafts (more rapid cooling) generate higher ice number densities, ice water content, cloud optical depth, and net radiative forcing. Increasing the height of the clouds in the model leads to an increase in ice number density, a decrease in effective radius, and a decrease in ice water content. The most prominent effect of increasing cloud height was a rapid increase in the net cloud radiative forcing which can be attributed to the change in cloud temperature as well as change in cloud ice size distributions. It has long been recognized that changes in cloud height or cloud area have the greatest potential for causing feedbacks on climate change. Our results suggest that variations in vertical velocity or cloud microphysical changes associatd with cloud height changes may also be important.

Jensen, Eric J.↗

A two-dimensional model study of the QBO signal in SAGE II NO2 and O3

Calculations of the quasi biennial oscillation (QBO) signal in Stratospheric Aerosol and Gas Experiment (SAGE) II O3 and NO2 data between 1984 and 1991 are presented and have been investigated by using a two-dimensional model. The isentropic 2D model is a fully interactive radiative-dynamical-chemical model in which the eddy fluxes of chemical species are calculated in a consistent manner. The QBO in the model has been forced by relaxing the equatorial zonal wind toward the observations at Singapore allowing the comparison of the model with observations from specific years. The model reproduces the observed vertical structure of the equatorial ozone anomaly with the well-known transition from dynamical to photochemical control at around 28km. The model also reproduces the observed vertical structure of the SAGE II observed NO2 anomaly. The model studies have shown that it is the QBO modulation of NO2 which the main cause of QBO signal in O3 above 30 km. The model also reproduces the observed latitudinal structure of the QBO signals in O3 and NO2. Due to the differing horizontal distribution of O3 and NO(y) the ozone signal shows a distinct phase change in the subtropics whereas the NO2 anomaly gives a broader signal.

Chipperfield, M. P.↗