Search NASA⌕ Search

SEARCH · Search NASA

Results for “binary fluids”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Giant Planet Scatterings and Collisions: Hydrodynamics, Merger-Ejection Branching Ratio, and Properties of the Remnants

Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple “sticky-sphere” prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and N-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with MJ and 2MJ, we find that quick mergers occur when the impact parameter rp (the pericenter separation between the planets) is less than 2RJ, and the merger conserved at least 97% of the initial mass. Strong tidal effects can affect the “binary-planet” orbit when rp is between 2RJ and 4RJ. We quantify these effects using a set of fitting formulae that can be implemented in N-body codes. We run a suite of N-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular and nearly coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the branching ratio change. We also find that the merger products have broad distributions of spin magnitudes and obliquities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Efficient data-driven models for prediction and optimization of geothermal power plant operations

Increasing the capacity of geothermal energy as a renewable resource calls for development and deployment of efficient control and optimization technologies for geothermal power plants. A data-driven prediction and optimization model is presented as a cost-effective and efficient alternative to physics-based approach. The model predicts power output and operational cost by propagating the influence of control and disturbance variables within an artificial neural network (ANN). Numerical experiments with simulated and field data from a real geothermal power plant are first used to demonstrate the prediction performance of the ANN model. The model is then adopted to maximize the net predicted power production by automatically adjusting the working fluid circulation rate. The optimization performance of the model in evaluated using a thermodynamic flowsheet simulation model. The workflow is applied to model and control the effect of ambient temperature on an air-cooled binary cycle power plant, which is complex and costly to perform using a physics-based predictive model. As a result, the performance of the method is demonstrated by applying it to both simulated and field datasets from a binary cycle geothermal power plant.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2024

Accurate and Data‐Efficient Micro X‐ray Diffraction Phase Identification Using Multitask Learning: Application to Hydrothermal Fluids

Traditional analysis of highly distorted micro X‐ray diffraction (μ‐XRD) patterns from hydrothermal fluid environments is a time‐consuming process, often requiring substantial data preprocessing and labeled experimental data. Herein, the potential of deep learning with a multitask learning (MTL) architecture to overcome these limitations is demonstrated. MTL models are trained to identify phase information in μ‐XRD patterns, minimizing the need for labeled experimental data and masking preprocessing steps. Notably, MTL models show superior accuracy compared to binary classification convolutional neural networks. Additionally, introducing a tailored cross‐entropy loss function improves MTL model performance. Most significantly, MTL models tuned to analyze raw and unmasked XRD patterns achieve close performance to models analyzing preprocessed data, with minimal accuracy differences. This work indicates that advanced deep learning architectures like MTL can automate arduous data handling tasks, streamline the analysis of distorted XRD patterns, and reduce the reliance on labor‐intensive experimental datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Bimodal particle distributions with increased thermal conductivity for solid particles as heat transfer media and storage materials

Solid particles are being considered in several high temperature thermal energy storage systems and as heat transfer media in a variety of advanced power generation systems, particularly in concentrated solar power plants. The downside of such an approach is the low overall heat transfer coefficients caused by the inherently low thermal conductivity values of the low-cost solid media when coupled to heat exchanger for the power cycle working fluid. Choosing the right particle size distribution, emittance, and material of the solid media can all make a substantial difference in packed bed thermal conductivity. Current research though exclusively focuses on continuous unimodal distributions of particles. Here, in this work, we propose the use of a binary particle system with a bimodal size distribution to significantly increase packed bed thermal conductivity by reducing packed bed porosity. This is the first study related to ceramic solid particle heat transfer that has considered the thermal conductivity of non-unimodal size distributions at room and elevated temperatures. The following study found that for the binary particle system using Carbo particles CP 16/30 – CP 70/140 where the large particle volume fraction was 50% there was an 17–47% increase in packed bed thermal conductivity when compared with a nearly unimodal particle size distribution of CP 16/30 between 50 and 300 °C. Two different porosity and effective thermal conductivity models were studied, with one providing better prediction of porosity but both effective thermal conductivity models providing less predictive capacity. Importantly this approach can have a substantial impact of thermal performance, with little to no impact on the particle cost.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions

Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.

Mechanics↗
DOE OSTI2024

Geometric and thermodynamic characterization of binary neutron star accretion discs

Accretion disks formed in binary neutron star mergers play a central role in many astrophysical processes of interest, including the launching of relativistic jets or the ejection of neutron-rich matter hosting heavy element nucleosynthesis. In this work we analyze in detail the properties of accretion disks from 44 ab initio binary neutron star merger simulations for a large set of nuclear equations of state, binary mass ratios, and remnant fates, with the aim of furnishing reliable initial conditions for disk simulations and a comprehensive characterization of their properties. We find that the disks have a significant thermal support, with an aspect ratio decreasing with the mass ratio of the binary from ~0.7 to 0.3. Even if the disk sample spans a broad range in mass and angular momentum, their ratio is independent from the equation of state and from the mass ratio. This can be traced back to the rotational profile of the disk, characterized by a constant specific angular momentum (as opposed to a Keplerian one) of 3 – 5×10 16 cm 2 s –1 . The profiles of the entropy per baryon and of the electron fraction depend on the mass ratio of the binary. For more symmetric binaries, they follow a sigmoidal distribution as a function of the rest mass density, for which we provide a detailed description and a fit. Furthermore, the disk properties discussed in this work can be used as a robust set of initial conditions for future long-term simulations of accretion disks from binary neutron star mergers, posing the basis for progress in the quantitative study of the outflow properties.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
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↗
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↗
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 OSTI2024

EGS Collab Experiment 2: Microseismic Monitoring

This dataset contains continuous seismic waveform data recorded during stimulation and thermal circulation tests for the Enhanced Geothermal Systems (EGS) Collab Experiment #2, conducted from February to September 2022 at the Sanford Underground Research Facility in Lead, South Dakota. This experiment aimed to study and validate models of geothermal systems by injecting high-pressure fluids into rock formations 1200-1500 meters below the surface, inducing microseismic events. The seismic monitoring system included 16 three-component accelerometers and a 24-channel hydrophone array, installed in boreholes surrounding the test area. Data were recorded at high sampling rates using a continuous waveform recording system to monitor seismic activity in real time. The dataset contains the raw data stored in binary format, with files named based on timestamps, and includes calibration certificates for some sensors to facilitate corrections to real units. Users are strongly advised to consult the accompanying detailed report, which outlines the experimental setup, sensor specifications, installation procedures, and data processing methods. The report also describes important nuances, such as the hardware filters on hydrophones, sensor calibration details, and the naming conventions for the recorded data. Proper use of this dataset may require familiarity with seismic data analysis tools, such as the Obspy Python package, and an understanding of the SEED naming conventions used for channel identification.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2022

Vapor–liquid equilibrium estimation of n-alkane/nitrogen mixtures using neural networks

Understanding fluid phase behavior, like VLE, in high P&T conditions is crucial for developing high-fidelity simulations of chemically reacting flows in liquid-fueled combustion systems and also forms an integral part of the design-modeling of the control processes in chemical industries. Two data-driven models have been proposed here in this study, each of which was competent in estimating VLE for the Type III binary systems of C 10 /N 2 and C 12 /N 2 , at pressures ranging up to 50–60 MPa. Both models showed better performance in predicting equilibrium pressure as compared to VLE modeled using PR-EOS. A modified model has also been proposed, capable of estimating the full phase envelope for the binary systems of C 10 /N 2 and C 12 /N 2 across a wide range of temperatures, and thus exhibit the mixture critical pressure at the concerned temperature. The diverse applicability of the proposed network architecture was further exhibited while estimating the VLE of a ternary system of C 1 /C 10 /N 2 .

97 MATHEMATICS AND COMPUTING↗