Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “directional heat flow”

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 325 records · Page 18

Stability of compressible reacting mixing layer

Linear instability of compressible reacting mixing layers is analyzed with emphasis on the effects of heat release and compressibility. Laminar solutions of the compressible boundary-layer equations are used as the base flows. The parameters of this study are the adiabatic flame temperature, the Mach number of the upper stream, frequency, wavenumber, and the direction of propagation of the disturbance wave. Stability characteristics of the flow are presented. Three groups of unstable modes are found when the Mach number and/or heat release are large. Finally, it is shown that the unstable modes are two-dimensional for large heat release even in highly compressible flow.

Shin, D. S.↗

Small scale demand type neon liquefaction plant

Low-temperature measurement of the thermal conductivity of insulating materials is generally made using a boil-off calorimetry technique involving liquid hydrogen (LH2). Liquid neon (LNe) has nearly the same normal boiling point as LH2, but has a much larger heat of vaporization, allowing extended run times. The main drawback of using LNe has been its excessive cost; $170.00 versus $1.50/l for LH2 (1989 prices). A neon liquefaction plant has been designed and constructed to capture, purify, and refrigerate the neon boil-off from calorimetry experiments. Recycling the neon reduces operating costs to approximately $20/l. The system consists of a purification section, a heat exchanger, LNe and LH2 storage dewars, and a fully automated control system. After purification, neon is liquified in the heat exchanger by LH2 flowing countercurrently through stainless steel cooling coils. Hydrogen flow is automatically adjusted to keep the neon at its normal saturation temperature, 27 K. The liquid neon is then stored in a dewar placed directly below the heat exchanger.

Dube, W. P.↗

Approaches for Water Removal in Direct-Fired sCO 2 Power Cycles

There is interest in investigation of water removal processes in direct fired sCO 2 flows, as this may potentially lead to greater system efficiency as the removal of this contaminant will result in sCO 2 behaving close to idealized behaviors. Water removal should be split into a two-step process, condensation of the water, followed by separation of the liquid phase water from the sCO 2 . The two main avenues of condensation are manipulation of pressure and temperature for phase change. For this paper, temperature-based phase change is the primary focus through the implementation of heat exchangers. Of the heat exchangers investigated it was found that printed circuit heat exchangers (PCHEs) could be an alternative for this use case, though the specific design of flow channel geometry and flow direction depends on the specific system case and cannot be determined at this point. For water separation there were four processes identified, all of which already assume water is in liquid phase at that point in the system. Of these separation avenues the best candidate is the hydrocyclone as it has a proven history of separating liquid-liquid phase mixtures with small density differences in oilfield use, in addition they have been investigated and modeled specifically for water separation for sCO 2 flows and the footprint is relatively small.

20 FOSSIL-FUELED POWER PLANTS↗

Developments on two lithium vapor-box linear test-stand experiments

The lithium vapor-box divertor is a possible fusion power exhaust solution.It uses condensation pumping to create a gradient of vapor density in a divertor slot; this should allow a stable detachment front without active feedback.As initial explorations of the concept, two test stands which take the form of three connected cylindrical stainless steel boxes are being developed: one without plasma at PPPL, to test models of lithium evaporation and flow; and one for the linear plasma device Magnum-PSI (at DIFFER in Eindhoven, The Netherlands) to test the ability of a lithium vapor cloud to induce volumetric detachment and redistribute the plasma power.The first experiment uses boxes with diameters of 6 cm, joined by apertures with diameters of 2.2 cm. Up to 1 g of Li is placed in one box, which is heated to up to 600 degrees C. The Li evaporates, then flows to and condenses in the two other, cooler boxes over several minutes. The quantity of Li transported is assessed by weighing the boxes before and after the heating cycle, and is compared to the quantity predicted to flow for the box at its measured temperature using a Direct Simulation Monte Carlo code, SPARTA. With good experimental conditions, the two values agree to within 15%.The experiment on Magnum-PSI is in the conceptual design stage.The design is assessed by simulations using the code B2.5-Eunomia.They show that when the hydrogen-ion plasma beam, with n_e = 4e20 per cubic meter, T_e = 1.5 eV, and r = 1 cm, is passed through a 16 cm long, 12 Pa, 625 degree C Li vapor cloud, the plasma heat flux and pressure on the target are significantly reduced compared to the case without Li.With the Li present, the plasma is cooled by excitation of Li neutrals followed by radiation until it volumetrically recombines, lowering the heat flux from 3.7 MW/m^2 to 0.13 MW/m^2, and the pressure is reduced by 93%, largely by collisions of hydrogen ions with neutral Li.

divertor↗

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗

Internal-Film Cooling of Rocket Nozzles

Experiments were conducted with 1000-pound-thrust rocket engine to determine feasibility of cooling convergent-divergent nozzle by internal film of water introduced at nozzle entrance. Water flow of 3 percent of propellant flow reduced heat flow into nozzle to 55 percent of uncooled heat flow. Introduction of water by porous ring before nozzle resulted in more uniform coverage of nozzle than water introduced by single arrangement of 36 jets directed along nozzle wall. Water flow through porous ring of 3.5 percent of propellant flow stabilized wall temperature in convergent section but did not adequately cool throat or divergent sections.

HEAT TRANSFER↗

Analytical studies of hypersonic viscous dissociated flows

This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Inger, George R.↗

Numerical Studies of Turbulent Supersonic Plane-Channel Flows

This two-part talk begins with an overview of the `turbulence problem’, the computational strategies used to address it, and the role of direct numerical simulation (DNS). The second part presents DNS results to determine the turbulent Prandtl number Pr t above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number Reτw from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of Pr t away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall Pr t profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗wall scaling proposed by Huang et al. (1995), with only a weak dependence on R eτw . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.

Direct numerical simulation↗

Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility

In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Verification and Validation of Spectral Element Code for Supercritical CO2 Flow in Vertical Heated Tubes

The investigation of heat transfer in supercritical CO2 (sCO2) has garnered considerable attention in recent decades, given sCO2's potential as a promising working fluid for advanced power conversion cycles. Despite previous research efforts, there are still gaps in our understanding of sCO2 heat transfer, particularly in conditions associated with heat transfer deterioration. To delve into sCO2 heat transfer more comprehensively, we propose employing the high-fidelity computational fluid dynamics code NekRS to simulate sCO2 flow using the large eddy simulation technique. Through graphics processing unit acceleration, NekRS achieves a higher computational speed than traditional CPU-based systems. However, before using NekRS in practical applications involving sCO2, it is imperative to perform verification and validation. Here, this paper presents our efforts to verify and validate the NekRS code's capability for simulating sCO2 using heated vertical tubes, where heat transfer deterioration usually happens. To accommodate the unique properties of sCO2, we have modified the NekRS code by integrating third-party property modules, such as REFPROP and PROPATH. Our simulations are compared with experimental and numerical data from the literature, instilling confidence in leveraging NekRS for future engineering applications. Our simulations also reveal that the accuracy of the property module significantly impacts the results, with REFPROP outperforming PROPATH for sCO2 properties. Additionally, we observed that, depending on the flow direction, buoyancy can either enhance or suppress turbulence in sCO2 flow. In upward flow, under certain conditions, the suppressed turbulence leads to heat transfer deterioration, resulting in elevated wall temperatures.

NekRS↗

A simple method for obtaining a three-dimensional proton distribution function from Voyager plasma data

The main sensor of the Vogager plasma experiment consists of a cluster of three, modulated-grid Faraday cups whose normals are arranged symmetrically about the symmetry axis of the cluster at an angle of 20 degrees to that axis. In interplanetary space, each cup explores the positive ion distribution by accepting particles from contiguous slices in velocity space. The slices are narrow in the direction of the normal to the modulating grid but are broad in planes parallel to that grid. The resulting three sets of measurements can be combined to yield the three-dimensional distribution function in the following way: the distribution function is assumed to be gyrotropic. For each value of speed in a frame of reference moving with the bulk velocity of the solar wind, the variation of the distribution function with angle from the field direction is represented by a series of Legendre polynomials. Effects such as double-streaming and heat flow can be well represented by using only the first three terms of the series which are fully specified by the measurements. Examples of the use of this method in the analysis of Voyager data are shown.

Olbert, S.↗

A hydromagnetic vortex seen by ISEE-1 and 2

Magnetometer and plasma data from the dual ISEE spacecraft are combined in a study of the initial plasma vortex event reported by Hones et al. (1978) in the dawn plasma sheet. The event is a transient hydromagnetic wave of two cycles duration with a six minute period. Large amplitude compressional and transverse magnetic components were present. Particle and magnetic pressure oscillations were in strict antiphase, but did not balance. When combined with the plasma velocity data these properties show that substantial Earthward field-aligned flows of electromagnetic energy and heat flux occurred during the vortex. The net energy flow perpendicular to B was in the antisolar direction. This event possesses hydromagnetic features unique to a hot plasma environment.

Saunders, M. A.↗

Small particle transport across turbulent nonisothermal boundary layers

The interaction between turbulent diffusion, Brownian diffusion, and particle thermophoresis in the limit of vanishing particle inertial effects is quantitatively modeled for applications in gas turbines. The model is initiated with consideration of the particle phase mass conservation equation for a two-dimensional boundary layer, including the thermophoretic flux term directed toward the cold wall. A formalism of a turbulent flow near a flat plate in a heat transfer problem is adopted, and variable property effects are neglected. Attention is given to the limit of very large Schmidt numbers and the particle concentration depletion outside of the Brownian sublayer. It is concluded that, in the parameter range of interest, thermophoresis augments the high Schmidt number mass-transfer coefficient by a factor equal to the product of the outer sink and the thermophoretic suction.

Rosner, D. E.↗

A three-dimensional spectral algorithm for simulations of transition and turbulence

A spectral algorithm for simulating three dimensional, incompressible, parallel shear flows is described. It applies to the channel, to the parallel boundary layer, and to other shear flows with one wall bounded and two periodic directions. Representative applications to the channel and to the heated boundary layer are presented.

Zang, T. A.↗

Fuel-Cell Structure Prevents Membrane Drying

Embossed plates direct flows of reactants and coolant. Membrane-type fuel-cell battery has improved reactant flow and heat removal. Compact, lightweight battery produces high current and power without drying of membranes.

Mcelroy, J.↗

Two-phase pressure drop of ammonia in small diameter horizontal tubes

Data for pressure drop in adiabatic two-phase ammonia flows in small diameter horizontal tubes are presented. The data has direct application to the sizing of the flow-through radiator tubes in the Space Station Freedom heat rejection system. The data are compared to existing correlations for pressure drop and are found to be significantly lower than the most commonly used correlations. However, several of the less commonly used correlations predict the data accurately. Alternate pressure drop prediction methods are explored and a recommendation is made for a method to accurately predict the pressure drop in two-phase ammonia flows in small horizontal tubes.

Ungar, Eugene K.↗

Pool and flow boiling in variable and microgravity

As is well known, boiling is an effective mode of heat transfer in that high heat flux levels are possible with relatively small temperature differences. Its optimal application requires that the process be adequately understood. A measure of the understanding of any physical event lies in the ability to predict its behavior in terms of the relevant parameters. Despite many years of research the predictability of boiling is currently possible only for quite specialized circumstances, e.g., the critical heat flux and film boiling for the pool boiling case, and then only with special geometries. Variable gravity down to microgravity provides the opportunity to test this understanding, but possibly more important, by changing the dimensional and time scales involved permits more detailed observations of elements involved in the boiling process, and perhaps discloses phenomena heretofore unknown. The focus here is on nucleate boiling although, as will be demonstrated below, under but certain circumstances in microgravity it can take place concurrently with the dryout process. In the presence of earth gravity or forced convection effects, the latter process is usually referred to as film boiling. However, no vapor film as such forms with pool boiling in microgravity, only dryout. Initial results are presented here for pool boiling in microgravity, and were made possible at such an early date by the availability of the Get-Away-Specials (GAS). Also presented here are some results of ground testing of a flow loop for the study of low velocity boiling, eventually to take place also in microgravity. In the interim, variable buoyancy normal to the heater surface is achieved by rotation of the entire loop relative to earth gravity. Of course, this is at the expense of varying the buoyancy parallel to the heater surface. Two questions which must be resolved early in the study of flow boiling in microgravity are (1) the lower limits of liquid flow velocity where buoyancy effects become significant to the boiling process (2) the effect of lower liquid flow velocities on the Critical Heat Flux when buoyancy is removed. Results of initial efforts in these directions are presented, albeit restricted currently to the ever present earth gravity.

Merte, Herman, Jr.↗

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 K) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

thermal protection↗