Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high pressure diffusion coefficients”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
NASA NTRS1992

Constraints on the thermal evolution of Venus inferred from Magellan data

The impact craters with diameters from 1.5 to 280 km compiled from Magellan observations indicate that the crater population on Venus has a completely spatially random distribution and the size/density distribution of craters with diameters greater than or equal to 35 km is consistent with a 'production' population with an age of 500 plus or minus 250 m.y. The similarity in size distribution from area to area indicates that the crater distribution is independent of crater size. Also, the forms of the modified craters are virtually identical to those of the pristine craters. These observations imply that Venus reset its cratering record by global resurfacing 500 m.y. ago, and resurfacing declined relatively fast. The fact that less than 40 percent of all craters have been modified and that the few volcanically embayed craters are located on localized tectonic regions indicate that only minor and localized volcanism and tectonism have occurred since the latest vigorous resurfacing event approximately 500 m.y. ago and the interior of Venus has been solid and possibly colder than Earth's. This is because the high-temperature lithosphere of Venus would facilitate upward ascending of mantle plumes and result in extensive volcanism if the venusian upper mantle were as hot as or hotter than Earth's. Therefore, the present surface morphology of Venus may provide useful constraints on the pattern of that vigorous convection, and possibly on the thermal state of the venusian mantle. We examine this possibility through numerical calculations of three-dimensional thermal convection models in a spherical shell with temperature- and pressure-dependent Newtonian viscosity, temperature-dependent thermal diffusivity, pressure-dependent thermal expansion coefficient, and time-dependent internal heat production rate solar magnitude.

Arkani-Hamed, Jafar↗
NASA NTRS1983

High altitude Venus haze from Pioneer Venus limb scans

High-resolution limb scans of Venus made with the Orbiter Cloud Photopolarimeter aboard Pioneer Venus have been used to observe the high altitude Venus haze and its variation over almost one half a Venusian year. A model for a spherical atmosphere with tangential viewing geometry was used to analyze the observations and derive the haze properties. Values are found for the particle size and refractive index, particle number density and scale height, and eddy diffusion coefficient. Single scattering computations show good agreement with observed intensities for particles smaller than 0.3 micron radius and refractive index less than 1.7, consistent with, but not limited to, concentrated sulfuric acid. Particle scale height in the 0.5 to 2 mbar pressure regions varies between 1 and 3 km over the season, latitude, and local time ranges of the observations. Detached layers of haze are sometimes present. An average particle scale height of 2.2 km at 84 km altitude yields an eddy diffusion coefficient of 130,000 sq cm/sec.

Lane, W. A.↗
NASA NTRS1986

Experimental evaluation of two turning vane designs for high-speed corner of 0.1-scale model of NASA Lewis Research Center's proposed altitude wind tunnel

Two turning vane designs were experimentally evaluated for corner 1 (downstream of the test section) of a 0.1-scale model of the NASA Lewis Research Center's proposed Altitude Wind Tunnel (AWT). Vane A was a controlled-diffusion airfoil shape; vane B was a circular-arc airfoil shape. The vane designs were tested over corner inlet Mach numbers from 0.16 to 0.465. Several modifications in vane setting angle and vane spacing were also evaluated for vane A. The overall performance obtained from total pressure rakes indicated that vane B had a slightly lower loss coefficient than vane A. At Mach 0.35 (the design Mach number without the engine exhaust removal scoop), the loss coefficients were 0.150 and 0.178 for vanes B and A, respectively. Resetting the vane A angle by -5 deg. (vane A10) to turn the flow toward the outside corner reduced the loss coefficient to 0.119. The best configuration (vane A10) was also tested with a simulated engine exhaust removal scoop. The loss coefficient for that configuration was 0.164 at Mach 0.41 (the approximate design Mach number with the scoop).

Moore, R. D.↗
NASA NTRS2001

The Cool Flames Experiment

A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.

Pearlman, Howard↗
NASA NTRS2019

Two-Stage Autoignition Dynamics of N-Dodecane Droplets Under Normal Gravity at High Pressures

Experimental observations of two-stage, autoignition characteristics of fiber-supported, n-dodecane fuel droplets under normal gravity are reported for a wide range of ambient temperatures, and pressures varying from 1 atm to 25 atm, which encompasses the fuel critical pressure (~ 18 atm). A high-speed camera is used to obtain backlit images of the autoignition process. During two-stage ignition, a cool-flame front forms below the droplet in the natural convective plume of the fuel vapor originating from the vaporizing droplet, and propagates upward, followed by the formation of a hot-flame kernel that expands rapidly and surrounds the droplet, eventually leading to the classical diffusion flame. Experimentally measured cool, and hot flame formation locations and propagation speeds are presented for varying initial temperatures and pressures. Measured first, and second stage ignition delay times are also reported. The influences of negative temperature coefficient (NTC) regime chemistry on flame starting location, propagation speeds, and ignition delay times are described. The first-stage cool flame characteristics are found to be controlled by both the chemistry, and physical aspects such as droplet evaporation and fuel-oxidizer mixing, while the second-stage hot flame ignition is primarily controlled by chemical kinetics. Further work is currently underway to develop models to quantitatively predict the observed behavior.

Rose, Evan↗
NASA NTRS1990

Solubility and diffusion of hydrogen in IN-100 and Waspaloy

The solubility of hydrogen, at a fixed hydrogen pressure, was measured in two different heat treatment conditions in both IN-100 and Waspaloy by equilibrium and analysis. The diffusion coefficients were determined by measuring the rate of hydrogen evolution in a vacuum. Hydrogen evolution appeared to be a diffusion-controlled process in all cases. For hydrogen diffusion in Waspaloy, Delta H was 38.5 kJ and D0 was 2.6 x 10 to the -7th sq/s. For IN-100, Delta H was 68.2 kJ and D0 was 59 x 10 to the -7th sq/s. The values for Waspaloy are close to those for pure nickel. The values for IN-100 are more similar to high nickel alloys that show trapping or surface layer impedance.

Khan, A. S.↗
NASA NTRS1972

Analysis of inner and outer zone: OGO-1 and OGO-2 electron spectrometer and ion chamber data

The dynamic processes governing the acceleration and loss of electrons in the radiation zones are investigated. The radial diffusion coefficient was determined for a McIlwain parameter between 1.6 and 2.2 for electrons having a first adiabatic invariant of 12 MeV/gauss. The coefficient is larger than earlier values and suggests that there exists a lower limit to the fluxes in the inner zone. The agreement between observed and calculated magnetic fields and particle fluxes is improved by using solar wind pressure as input to the magnetic field models. Changes in the plasma pressure can cause apparent local time asymmetries in particle flux. A comparison of the magnetic field models with observed location of the trapping boundary also indicates the need for including distributed currents within the magnetosphere. The high latitude trapping boundary is only weakly dependent on A sub p, and the trapping boundary data are improved by including in the models a stand-off distance which varies with the plasma pressure.

Pfitzer, K. A.↗
NASA NTRS1993

Generic hypersonic vehicle performance model

An integrated computational model of a generic hypersonic vehicle was developed for the purpose of determining the vehicle's performance characteristics, which include the lift, drag, thrust, and moment acting on the vehicle at specified altitude, flight condition, and vehicular configuration. The lift, drag, thrust, and moment are developed for the body fixed coordinate system. These forces and moments arise from both aerodynamic and propulsive sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be determined. The vehicle is assumed to be a lifting body with a single aerodynamic control surface. The body shape and control surface location are arbitrary and must be defined. The aerodynamics are calculated using either 2-dimensional Newtonian or modified Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. Skin-friction drag was also accounted for. The skin-friction drag coefficient is a function of the freestream Mach number. The data for the skin-friction drag coefficient values were taken from NASA Technical Memorandum 102610. The modeling of the vehicle's SCRAMjet engine is based on quasi 1-dimensional gas dynamics for the engine diffuser, nozzle, and the combustor with heat addition. The engine has three variable inputs for control: the engine inlet diffuser area ratio, the total temperature rise through the combustor due to combustion of the fuel, and the engine internal expansion nozzle area ratio. The pressure distribution over the vehicle's lower aft body surface, which acts as an external nozzle, is calculated using a combination of quasi 1-dimensional gas dynamic theory and Newtonian or modified Newtonian theory. The exhaust plume shape is determined by matching the pressure inside the plume, calculated from the gas dynamic equations, with the freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this manner, the pressure distribution along the vehicle after body expansion surface is then determined. The aerodynamic modeling, the engine modeling, and the exhaust plume analysis are described in more detail. A description of the computer code used to perform the above calculations is given and an input/output example is then given. The computer code is available on a Macintosh floppy disk.

Chavez, Frank R.↗
NASA NTRS1943

Heat Transfer Over the Circumference of a Heated Cylinder in Transverse Flow

A method for recording the local heat-transfer coefficients on bodies in flow was developed. The cylinder surface was kept at constant temperature by the condensation of vapor except for a narrow strip which is heated separately to the same temperature by electricity. The heat-transfer coefficient at each point was determined from the electric heat output and the temperature increase. The distribution of the heat transfer along the circumference of cylinders was recorded over a range of Reynolds numbers of from 5000 to 426,000. The pressure distribution was measured at the same time. At Reynolds numbers up to around 100,000 high maximums of the heat transfer occurred in the forward stagnation point at and on the rear side at 180C, while at around 80 the heat-transfer coefficient on both sides of the cylinder behind the forward stagnation point manifested distinct minimums. Two other maximums occurred at around 115 C behind the forward stagnation point between 170,000 and 426,000. At 426,000 the heat transfer at the location of those maximums was almost twice as great as in the forward stagnation point, and the rear half of the cylinder diffused about 60 percent of the entire heat, The tests are compared with the results of other experimental and theoretical investigations.

Schmidt, Ernst↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
NASA NTRS1980

A physical model of Titan's clouds

A physical model of the formation and growth of aerosols in the atmosphere of Titan has been constructed in light of the observed correlation between variations in Titan's albedo and the sunspot cycle. The model was developed to fit spectral observations of deep methane bands, pressures, temperature distributions, and cloud structure, and is based on a one-dimensional physical-chemical model developed to simulate the earth's stratospheric aerosol layer. Sensitivity tests reveal the model parameters to be relatively insensitive to particle shape but sensitive to particle density, with high particle densities requiring larger aerosol mass production rates to produce compatible clouds. Solution of the aerosol continuity equations for particles of sizes 13 A to about 3 microns indicates the importance of a warm upper atmosphere and a high-altitude mass injection layer, and the production of aerosols at very low aerosol optical depths. Limits are obtained for the chemical production of aerosol mass and the eddy diffusion coefficient, and it is found that an increase in mass input causes a decrease in mean particle size.

Toon, O. B.↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗
DOE OSTI2022

Evidence for superionic H 2 O and diffusive He–H 2 O at high temperature and high pressure

In this work, we present the evidence of superionic phase formed in H 2 O and, for the first time, diffusive H 2 O–He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H 2 O and diffusive He–H 2 O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H 2 O phase forms gradually in pure H 2 O over the temperature range of 1350–1400 K at 23 GPa, but the diffusive He–H 2 O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He–H 2 O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He–H 2 O.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1985

On boundary layer transition in high-subsonic and transonic flow under the influence of acoustic disturbances and free-stream turbulence

The existing data base which quantifies the roles played by the Mach number, turbulence level, and acoustic noise in the flow quality of wind tunnels is reviewed. Attention is given to the effects on end (ET) and beginning (BT) transitions in flows over a 10 deg cone in transonic channels. The Re of both ET and BT flows were significantly affected by the Mach number, acoustic noise level and the noise frequency content. A unit increment in the Mach number increased the Re beyond 3,000,000, while noise below the 1 percent level in the rms pressure coefficient did not alter the transition phenomenon. More data are required on the effects of turbulence, which was a significant factor, and in the settling chamber, test section, plenum and diffuser to trace the locations and magnitudes of the acoustic sources.

Murthy, S. V.↗
NASA NTRS1972

Transport properties of liquid metal hydrogen under high pressures

A theory is developed for the compressibility and transport properties of liquid metallic hydrogen, near to its melting point and under high pressure. The interionic force law is assumed to be of the screened Coulomb type, because hydrogen has no core electrons. The random phase approximation is used to obtain the structure factor S(k) of the system in terms of the Fourier transform of this force law. The long wavelenth limit of the structure factor S(o) is related to the compressibility, which is much lower than that of alkali metals at their melting points. The diffusion constant at the melting point is obtained in terms of the Debye frequency, using a frequency spectrum analogous with the phonon spectrum of a solid. A similar argument is used to obtain the combined shear and bulk viscosities, but these depend also on S(o). The transport coefficients are found to be about the same size as those of alkali metals at their melting points.

Brown, R. C.↗