Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BINARY FLUID”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Solutal separation in a binary nanofluid due to thermodiffusion

Transport phenomena in porous media have received considerable attention due to an increasing interest in geothermal processes, chemical catalytic reactors, waste storage (especially geological or ocean storage of carbon dioxide), etc. Among others, oil industry has shown an increasing interest in studying diffusion phenomenon. Nanofluid is a term used to describe the suspension of low concentration of metallic and non-metallic nanoparticles in a base fluid. The size of a nanoparticle ranges from 10 to 100nm, and the conventional fluids used are water, ethylene glycol (C{sub 2}H{sub 6}O{sub 2}) or engine oil. Various studies have proven that nanoparticles improve the heat transfer of a base fluid. However, using various nanofluids it has been shown that the results could vary depending on different initial concentrations. The main objective of this paper is to study the diffusion and the thermodiffusion effect in a nanofluid for different fluid/porous media configurations. In this configuration, a liquid layer surrounds a porous layer. The full Brinkman equation coupled with the heat and mass transfer equations have been solved numerically for the porous layer using the finite element technique. The full Navier stokes equation coupled with heat and mass transfer equations have been solved for the liquid layer using the finite element method. A constraint between the liquid and porous layer has been applied to ensure heat flow and mass transfer continuity is maintained. A square cavity filled with hydrocarbon nanofluid of a mixture of fullerene-toluene with varying concentration of fullerene has been subject to different heating conditions. The entire cavity has been considered to be fully wetted with nanofluid. Results have confirmed that in the presence of a nanofluid a heat transfer enhancement is present up to certain initial concentration of the fullerene. The heat convection coefficient has been found to be 16% higher when a nanofluid is used as the working fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Main-ion intrinsic toroidal rotation across the ITG/TEM boundary in DIII-D discharges during ohmic and electron cyclotron heating

Direct measurements of deuterium main-ion toroidal rotation spanning the linear ohmic to saturated ohmic confinement (LOC-SOC) regime and with additional electron cyclotron heating (ECH) are proposed and compared with the more commonly measured impurity (carbon) ion rotation in DIII-D. Main ions carry the bulk of the plasma toroidal momentum, and hence, the shape of the main-ion rotation is more relevant to the study of angular momentum transport in tokamaks. Both in the LOC regime and with ECH, the main-ion toroidal rotation frequency is flat across the profile from the sawtooth region to the plasma separatrix. However, the impurity rotation profile possesses a rotation gradient, with the rotation frequency being lower near the plasma edge, implying a momentum pinch or negative residual stress inferred from the impurity rotation that differs from the main-ion rotation. In the SOC regime, both the main-ion and impurity rotation profiles develop a deeply hollow feature near the midradius while maintaining the offset in the edge rotation, both implying a positive core residual stress. In the radial region where the rotation gradient changes most dramatically, turbulence measurements show that density fluctuations near the trapped electron mode (TEM) scale are higher when the rotation profile is flat and drop significantly when the plasma density is raised and the rotation profile hollows, consistent with instabilities damped by collisions. Linear initial value gyrokinetic simulations with GYRO indicate that the transition from LOC-SOC in DIII-D occurs as TEMs are replaced by ion temperature gradient (ITG) driven modes from the outer radii inwards as the plasma collisionality increases, Z eff decreases, and the power flow through the ion channel progressively increases due to the electron-ion energy exchange. Gyrofluid modeling with trap gyro-Landau fluid (TGLF) successfully reproduces the plasma profiles at key times in the discharge and in time dependent simulations with predictive TRANSP. TGLF indicates that in the LOC and SOC regimes as well as with ECH, subdominant modes are present and that the plasma is not in a pure TEM or ITG binary state, but rather a more subtle mixed state. Predictions of the main-ion rotation profiles are observed with global nonlinear gyrokinetic simulations using GTS and reveal that the flat rotation is due to oscillatory variation of the turbulent residual stress across the profile, whereas the deeply hollow rotation profile is due to a larger-scale, dipole-like stress profile. In these cases, the predicted and observed main-ion rotation profile is consistent with the balance of turbulent residual stress and momentum diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Optimized collision-specific parameters for binary mixtures of nitrogen, oxygen, argon, and helium

Recently proposed collision-specific parameters for direct simulation Monte Carlo simulations are tested for binary mixtures of nitrogen, oxygen, and argon. Near ambient conditions, the traditional collision-averaged parameters are highly accurate, whereas the collision-specific parameters are not. The simulated transport using the collision-averaged parameters for mixtures with helium, however, is found to be inaccurate. Therefore, we propose a novel method to determine molecular parameters by combining the Chapman–Enskog theory with empirical mixing rules and experimental data. The optimized parameters are highly accurate for the binary mixtures of nitrogen, oxygen, and argon and greatly improve the simulated transport for the helium mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Modeling the use of a binary mixture as a control scheme for two-phase thermal systems

Two-phase thermal loops using mechanical pumps, capillary pumps, or a combination of the two have been chosen as the main heat transfer systems for the space station. For these systems to operate optimally, the flow rate in the loop should be controlled in response to the vapor/liquid ratio leaving the evaporator. By substituting a mixture of two non-azeotropic fluids in place of the single fluid normally used in these systems, it may be possible to monitor the temperature of the exiting vapor and determine the vapor/liquid ratio. The flow rate would then be adjusted to maximize the load capability with minimum energy input. A FLUINT model was developed to study the system dynamics of a hybrid capillary pumped loop using this type of control and was found to be stable under all the test conditions.

Benner, S. M.↗
DOE OSTI2019

Implications of real-gas behavior on refractive index calculations for optical diagnostics of fuel–air mixing at high pressures

Three models to compute the refractive index of gaseous mixtures at real gas conditions are presented with the purpose to improve the accuracy of state relationships between refractive index and thermodynamic properties. Models are compared with experimental data to determine one that is applicable to high pressure mixtures with non-ideal thermodynamic behavior near or above a fluids’ critical point. The optimal model is applied to analyze adiabatic thermal mixing of fuel and air at typical diesel engine conditions. Our results show that the ideal gas mixture law is appropriate in the vapor region, assuming it is known, for example, from experimental measurements. Finally, the model is applied for binary fuel-air mixing at supercritical conditions to demonstrate its full potential.

02 PETROLEUM↗
NASA NTRS1977

Cloud physics laboratory project science and applications working group

The conditions of the expansion chamber under zero gravity environment were simulated. The following three branches of fluid mechanics simulation under low gravity environment were accomplished: (1) oscillation of the water droplet which characterizes the nuclear oscillation in nuclear physics, bubble oscillation of two phase flow in chemical engineering, and water drop oscillation in meteorology; (2) rotation of the droplet which characterizes nuclear fission in nuclear physics, formation of binary stars and rotating stars in astrophysics, and breakup of the water droplet in meteorology; and (3) collision and coalescence of the water droplets which characterizes nuclear fusion in nuclear physics and processes of rain formation in meteorology.

Hung, R. J.↗
DOE OSTI2024

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Revealing 3D Morphological and Chemical Evolution Mechanisms of Metals in Molten Salt by Multimodal Microscopy

Growing interest in molten salts as effective high-temperature heat-transfer fluids for sustainable energy systems drives a critical need to fundamentally understand the interactions between metals and molten salts. This work utilizes the multimodal microscopy methods of synchrotron X-ray nanotomography and electron microscopy to investigate the 3D morphological and chemical evolution of two-model systems, pure nickel metal and Ni-20Cr binary alloy, in a representative molten salt (KCl-MgCl2 50–50 mol %, 800 °C). In both systems, unexpected shell-like structures formed because of the presence of more noble tungsten, suggesting a potential route of using Ni–W alloys for enhanced molten-salt corrosion resistance. The binary alloy Ni-20Cr developed a bicontinuous porous structure, reassembling functional porous metals manufactured by dealloying. This work elucidates better mechanistic understanding of corrosion in molten salts, which can contribute to the design of more reliable alloys for molten salt applications including next-generation nuclear and solar power plants and opens the possibility of using molten salts to fabricate functional porous materials.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Thermosolutal convection during directional solidification

During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly - a phenomenon also indicated by linear stability analysis. The lateral inhomogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time.

Mcfadden, G. B.↗
NASA NTRS1985

Double-diffusive convection with sidewalls

Stommel et al. (1956) have first described an instability, known as thermosolutal convection, thermohaline convection, or double-diffusive convection. This instability may occur in the case of a fluid in a gravitational field with two diffusing components present. The present study is concerned with the effect of sidewalls on flow in the fingering regime in the absence of applied horizontal gradients. The work was motivated by numerical results obtained on the basis of a simulation of thermosolutal convection occurring during the unidirectional solidification of a binary alloy. In this case, the unperturbed solute field in the liquid ahead of the solidifying planar interface has an exponential vertical profile because of the rejection or preferential incorporation of solute by the solid phase.

Mcfadden, G. B.↗
NASA NTRS1993

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗
NASA NTRS1996

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗
DOE OSTI2017

Species Separation and Other Kinetic Effects at Shock Wave Fronts in Plasma Mixtures Assessed using Molecular Dynamics

The propagation of a shock through a binary plasma mixture of hydrogen isotopes is simulated with molecular dynamics (MD) and a multispecies continuum hydrodynamics model. Separation of the two species due to their mass difference is observed at the shock front in both simulations. We analyze the shock structure in detail to determine which aspects of the wave propagation and shock profile obey single-fluid hydrodynamics and which require multi-component hydrodynamics with Fickian and barodiffusion. In order to reduce possible sources of uncertainty due to transport models used in the continuum simulations, we have performed separate, equilibrium MD simulations which verify the accuracy of the recently proposed Kinetic Molecular Dynamics (KMD) viscosity model, even at the relatively cold pre-shock temperature of 1 eV. We find the post-shock temperature and density are generally consistent between continuum and atomic-scale calculations, but the MD simulations predict a significantly larger thickness of the leading light-ion shock front, with some associated differences in the species separation. We analyze the shape and anisotropy of the MD velocity distributions in detail and compute the electric field at the shock interface from both the MD and continuum simulations. We find that the continuum hydrodynamics calculation captures the major features of the electric fi eld through the electron pressure term in the barodiffusion model and conclude that the differences observed in the shock structures from the two simulations arise largely from non-Maxwellian effects in the longitudinal velocities (the velocities in the direction that the shock propagates) that are not captured in the hydrodynamic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

EXPERIMENTAL AND COMPUTATIONAL STUDY OF DIRECT CONTACT MEMBRANE DISTILLATION

Experimental and computational studies of direct contact membrane distillation (DCMD) are conducted. The permeate flux is measured in a module containing a flat sheet membrane. Polytetrafluoroethylene (PTFE) hydrophobic membrane (Membrane Solutions, FPB045A16) with a nominal pore size of 0.45 μm, the thickness of 179~239 μm, the porosity of 80% is used in the experiments. The feed solution is a binary mixture of water and NaCl. The inlet feed temperature is fixed at 70°C, and the feed Reynolds number is varied between 500 and 2000. Two inlet feed concentrations, 0 and 35g/L, are considered. The permeate flow rate was kept constant at Re of 1300 and temperature of 20°C. Three-dimensional steady-state computational fluid dynamics (CFD) simulations are performed using the laminar model. The Navier’s Stokes, energy, and mass transport equations in each channel coupled with flux boundary conditions imposed at the membrane surface are solved. The predicted module-averaged flux is compared to measured data for all operating conditions. Measured and predicted flux agrees well; validating both experiments and model, Temperature polarization in the feed and permeate channel and concentration polarization in the feed channel are characterized by validated CFD simulations.

Oztekin, Alparslan↗
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermodynamic modeling of double azeotropy systems

Azeotropes are boiling liquid mixtures which have identical vapor and liquid compositions. Known as polyazeotropy, azeotropic mixtures may exhibit more than one azeotrope at a given temperature or pressure. While the thermodynamic condition for azeotrope formation in vapor-liquid equilibrium is well established, thermodynamic modeling of mixtures exhibiting polyazeotropy remains challenging and prior modeling studies often failed to satisfactorily correlate the polyazeotropic behavior of associating systems. Based on the association NRTL model, this study presents successful correlation and analysis of double azeotrope formation for three diverse systems: strongly associating methanol + diethylamine binary, weakly associating ammonia + pentafluoroethane binary, and non-associating benzene + hexafluorobenzene binary systems. Furthermore, these systems exhibit polyazeotropic behavior as the temperature increases and approaches their Bancroft points.

42 ENGINEERING↗
DOE OSTI2016

SPREADING LAYERS IN ACCRETING OBJECTS: ROLE OF ACOUSTIC WAVES FOR ANGULAR MOMENTUM TRANSPORT, MIXING, AND THERMODYNAMICS

Disk accretion at a high rate onto a white dwarf (WD) or a neutron star has been suggested to result in the formation of a spreading layer (SL)—a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics, we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large-scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g., α-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low-mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective rearrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin–Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal nonuniformity of mixing in accreting WDs may be linked to the observed bipolar morphology of classical nova ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Triaxial degenerate dwarfs and pulsing X-ray sources

It is shown that triaxial configurations for rotating degenerate dwarfs consisting entirely of fluid are not limited to forms analogous to the Jacobi ellipsoids. Arguments are presented which rule out only the specialized class of uniformly rotating analogs of the Jacobi ellipsoid for some degenerate dwarfs. Several ways of forming triaxial degenerate dwarfs are examined along with the characteristic lifetimes of such objects. Noting that some of the alternative forms may exist in mass-transfer binary systems, consideration is given to the secondaries in Her X-1 and Cen X-3. It is suggested that the secondary in Her X-1 must be close to an analog of a Dedekind ellipsoid, and plausible triaxial white-dwarf configurations are proposed for the secondary in Cen X-3.

Durisen, R. H.↗