Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “critical point”

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 271 records · Page 15

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

Science Requirements for a Space Flight Experiment Entitled Critical Viscosity of Xenon

We propose to measure in low gravity the viscosity of xenon close to its critical point. The accuracy will be sufficient to eliminate uncertainties currently associated with the analysis of l-g experiments. The measurements will provide the first direct observation of the predicted power-law divergence of viscosity in a pure fluid. The measurements will also strengthen Zeno's test of mode coupling theory by greatly increasing the reliability of the extrapolation of viscosity to low reduced temperatures. Our scientific objectives are described in more detail in one of the attached reports. The low-gravity experiment will be the final stage of a program whose completed ground-based stages are: (1) theoretical studies by one of the principal investigators (MRM) and coworkers, (2) critical viscosity measurements of binary liquid mixtures, (3) critical viscosity measurements of pure fluids in l-g, and development of a suitable vibration-insensitive viscometer. Our technical approach is described in the draft Science Requirements Document. One of us (MRM) has reviewed opportunities for critical phenomena research in low gravity. Both of us were co-principal investigators in the Thermal Equilibration Experiment in the Critical Point Facility, flown on IML-1 in 1992. From this experience, and from the technical maturity of our ground-based work, we believe our critical point viscometer is ready for development as a flight experiment.

Berg, Robert F.↗

Gallery of soft modes: Theory and experiment at a ferromagnetic quantum phase transition

We examine the low-energy excitations in the vicinity of the quantum critical point in LiHoF4, a physical realization of the Transverse Field Ising Model, focusing on the long-range fluctuations which soften to zero energy at the ferromagnetic quantum phase transition. Microwave spectroscopy in tunable loop-gap resonator structures identifies and characterizes the electronuclear soft mode and higher-energy electronuclear states as a function of frequency and magnetic fields applied transverse and parallel to the Ising axis. These are understood in the context of a theoretical model of a soft electronuclear mode that interacts with soft photons as well as soft phonons. We identify competing infrared divergences at the quantum critical point, coming from the photons and the electronuclear soft mode. It is an incomplete cancellation of these divergences that leads to the muted but distinct signatures observed in the experiments. The application of a longitudinal magnetic field gaps the soft mode. As a result, measurements well away from the quantum critical point reveal a set of "Walker'' modes associated with ferromagnetic domain dynamics.

Dipolar interaction↗

Low-gravity experiments in critical phenomena

Studies of anomalous thermodynamic, transport, and structural phenomena in multibody systems near critical points are reviewed. The nomenclature used to describe critical points is explained; theoretical predictions of the thermodynamic properties of bulk systems are presented; and experimental tests of these predictions systems are discussed, considering equilibration and gravity effects in fluid systems and emphasizing the value of experiments conducted in a reduced-gravity environment. Several such experiments are described, and the available academic-research opportunities are briefly surveyed.

Moldover, Michael R.↗

Uniformly-valid asymptotic solutions to the Orr-Sommerfeld equation using multiple scales

The Orr-Sommerfeld equation is studied using the method of multiple scales. The leading-order uniformly-valid asymptotic solutions are obtained without recourse to heuristic methods for two different types of velocity profiles. Solutions were obtained for arbitrary velocity profiles without inflection or critical points and for profiles with simultaneous inflection and critical points. These are shown to reproduce previously obtained solutions in the inner and outer regions.

Long, Lyle N.↗

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗

Critical phenomena experiments in space

The paper analyzes several types of critical phenomena in fluids, shows how they are affected by the presence of gravity, and describes how experiments conducted in an orbiting laboratory under low gravity conditions could extend the range of measurements needed to study critical phenomena. Future experiments are proposed. One would be a careful measurement of the dielectric constant in a low gravity environment. Two basic problems that can benefit especially from space experiments are the specific heat near the critical point and the shear viscosity at the gas-liquid critical point.

Sengers, J. V.↗

Synergetic Ferroelectricity and Superconductivity in Zero-Density Dirac Semimetals near Quantum Criticality

Here, we study superconductivity in a three-dimensional zero-density Dirac semimetal in proximity to a ferroelectric quantum critical point. We find that the interplay of criticality, inversion-symmetry breaking, and Dirac dispersion gives rise to a robust superconducting state at the charge-neutrality point, where no Fermi surface is present. Using Eliashberg theory, we show that the ferroelectric quantum critical point is unstable against the formation of a ferroelectric density wave (FDW), whose fluctuations, in turn, lead to a first-order superconducting transition. Surprisingly, long-range superconducting and FDW orders are found to cooperate with each other, in contrast to the more usual scenario of phase competition. Therefore, we suggest that driving charge neutral Dirac materials, e.g., Pb x Sn 1–x Te, through a ferroelectric quantum critical point may lead to superconductivity intertwined with FDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dispersion-enhanced sequential batch sampling for adaptive contour estimation

In computer simulation and optimal design, sequential batch sampling offers an appealing way to iteratively stipulate optimal sampling points based upon existing selections and efficiently construct surrogate modeling. Nonetheless, the issue of near duplicates poses tremendous quandary for sequential learning. It refers to the situation that selected critical points cluster together in each sampling batch, which are individually but not collectively informative towards the optimal design. Near duplicates severely diminish the computational efficiency as they barely contribute extra information towards update of the surrogate. To address this issue, we impose a dispersion criterion on concurrent selection of sampling points, which essentially forces a sparse distribution of critical points in each batch, and demonstrate the effectiveness of this approach in adaptive contour estimation. Specifically, we adopt Gaussian process surrogate to emulate the simulator, acquire variance reduction of the critical region from new sampling points as a dispersion criterion, and combine it with the modified expected improvement (EI) function for critical batch selection. The critical region here is the proximity of the contour of interest. This proposed approach is vindicated in numerical examples of a two-dimensional four-branch function, a four-dimensional function with a disjoint contour of interest and a time-delay dynamic system.

97 MATHEMATICS AND COMPUTING↗

MISTE flight experiment status

The MISTE flight experiment has made significant scientific and technical progress in preparing for a future microgravity flight. We are collaborating with several theoretical modeling groups that have developed crossover (equation-of-state) models for predicting thermodynamic behavior near the liquid-gas critical point. Several of these models have already been used to test experimental measurements of the heat capacity at constant volume, isothermal susceptibility, and coexistence curve in the crossover region near the 3He liquid-gas critical point. A brief description of these models and a representative fit to experimental data will be presented. In collaboration with Mission Research Corporation, MISTE has been testing a new small pneumatic valve for use at low temperatures. The results of recent successful low temperature actuation tests will also be discussed.

critical phenomena liquid-gas critical point↗

Quantum criticality in the anisotropic spin-chain antiferromagnet Ca 2 ⁢CoH 2 ⁢(SeO 3 ) 4

Here, the quantum criticality of an XXZ-type spin chain antiferromagnet Ca 2 ⁢CoH 2 ⁢(SeO 3 ) 4 , subjected to an external magnetic field, was investigated through ultralow-temperature specific heat, magnetocaloric effect, and magnetic susceptibility measurements. In the absence of an external magnetic field, the system undergoes an antiferromagnetic phase transition at T N = 1.42K. The application of an external field gradually suppresses the antiferromagnetic order. Furthermore, an additional phase emerges in the intermediate field range when the zero-field antiferromagnetic order is suppressed, resulting in a quantum critical point at B s ≈ 2.55T. Near this field-induced quantum critical point, universal scaling behaviors are observed, characterized by the universality class parameters d=1, v=1/2, and z=2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lattice model theory of the equation of state covering the gas, liquid, and solid phases

The three stable states of matter and the corresponding phase transitions were obtained with a single model. Patterned after Lennard-Jones and Devonshires's theory, a simple cubic lattice model containing two fcc sublattices (alpha and beta) is adopted. The interatomic potential is taken to be the Lennard-Jones (6-12) potential. Employing the cluster variation method, the Weiss and the pair approximations on the lattice gas failed to give the correct phase diagrams. Hybrid approximations were devised to describe the lattice term in the free energy. A lattice vibration term corresponding to a free volume correction is included semi-phenomenologically. The combinations of the lattice part and the free volume part yield the three states and the proper phase diagrams. To determine the coexistence regions, the equalities of the pressure and Gibbs free energy per molecule of the coexisting phases were utilized. The ordered branch of the free energy gives rise to the solid phase while the disordered branch yields the gas and liquid phases. It is observed that the triple point and the critical point quantities, the phase diagrams and the coexistence regions plotted are in good agreement with the experimental values and graphs for argon.

Bonavito, N. L.↗

Partition function corrections

The vibration-rotation interaction coupling in diatomic molecules is analyzed and first-order expressions for the energy levels are deduced. Corrections to the partition function to account for this coupling are derived. The effects of long-range intermolecular forces due to induced dipole-induced dipole interactions, dipole-induced dipole interactions, charge-dipole interactions, and charge-charge interactions are considered. London's quantized oscillator model of the induced dipole-induced dipole interaction is developed. The partition function including effects of such dispersion forces is derived and is shown to lead to the virial equation of state. For a model combining the hard sphere and long-range attractive type potentials, this is found to lead to the van der Waals equation of state. The equation of state near the critical point is presented in reduced form (thermodynamic variables expressed in units of critical point values). Finally, approximate corrections to the electronic partition functions are derived to account for perturbations of excited electronic states induced by neighboring neutral particles, heavy ions, and fast electron collisions.

Source record↗

Critical exponent for the viscosity of carbon dioxide and xenon

The viscosities of carbon dioxide and xenon have been measured near their critical points and the critical exponent y characterizing the asymptotic divergence has been determined. Both fluids yielded exponents in the range y = 0.041 + or - 0.001 and thus also fell in the range y = 0.042 + or - 0.002 from an earlier study of four binary liquids. This agreement between experiments is the first evidence that pure fluids and binary liquids are in the same dynamic universality class. A recent theoretical value for y is 0.032. The 30 percent discrepancy is much greater than the combined errors from experiment and theory. The torsion oscillator viscometer operated at low frequency and low shear rate to avoid systematic errors caused by critical slowing down. Far from T(c) the analysis accounted for the crossover from critical to noncritical temperature dependence, where the latter was obtained from previously published correlations. Corrections for gravitational stratification were included close to T(c).

Berg, R. F.↗

Investigation of two and three parameter equations of state for cryogenic fluids

Two-phase flows are a common occurrence in cryogenic engines and an accurate evaluation of the heat-transfer coefficient in two-phase flow is of significant importance in their analysis and design. The thermodynamic equation of state plays a key role in calculating the heat transfer coefficient which is a function of thermodynamic and thermophysical properties. An investigation has been performed to study the performance of two- and three-parameter equations of state to calculate the compressibility factor of cryogenic fluids along the saturation loci. The two-parameter equations considered here are van der Waals and Redlich-Kwong equations of state. The three-parameter equation represented here is the generalized Benedict-Webb-Rubin (BWR) equation of Lee and Kesler. Results have been compared with the modified BWR equation of Bender and the extended BWR equations of Stewart. Seven cryogenic fluids have been tested; oxygen, hydrogen, helium, nitrogen, argon, neon, and air. The performance of the generalized BWR equation is poor for hydrogen and helium. The van der Waals equation is found to be inaccurate for air near the critical point. For helium, all three equations of state become inaccurate near the critical point.

Jenkins, Susan L.↗

Automated analysis for fluid flow topology

A new approach for visualizing vector data sets was developed by reducing the original vector field to a set of critical points and their connections, and was applied to fluid flow data sets. The critical point representation allows for considerable reduction in the data complexity. The representations are displayed as surfaces which are much simpler than the original data set, yet retain all the pertinent flow topology information. It is suggested that topological representations may be useful for database comparison.

Helman, James↗

Critical fluid thermal equilibration experiment (19-IML-1)

Gravity sometimes blocks all experimental techniques of making a desired measurement. Any pure fluid possesses a liquid-vapor critical point. It is defined by a temperature, pressure, and density state in thermodynamics. The critical issue that this experiment attempts to understand is the time it takes for a sample to reach temperature and density equilibrium as the critical point is approached; is it infinity due to mass and thermal diffusion, or do pressure waves speed up energy transport while mass is still under diffusion control. The objectives are to observe: (1) large phase domain homogenization without and with stirring; (2) time evolution of heat and mass after temperature step is applied to a one phase equilibrium sample; (3) phase evolution and configuration upon going two phase from a one phase equilibrium state; (4) effects of stirring on a low g two phase configuration; (5) two phase to one phase healing dynamics starting from a two phase low g configuration; and (6) effects of shuttle acceleration events on spatially and temporally varying compressible critical fluid dynamics.

Wilkinson, R. Allen↗