Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two Phase Flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Effect of Aluminum Alloying on the Hot Deformation Behavior of Nano-bainite Bearing Steel

Interest in using aluminum in nano-bainite steel, especially for high-carbon bearing steel, is gradually growing. In this study, GCr15SiMo and GCr15SiMoAl steels are introduced to investigate the effect of Al alloying on the hot deformation behavior of bearing steel. Results show that the addition of Al not only notably increases the flow stress of steel due to the strong strengthening effect of Al on austenite phase, but also accelerates the strain-softening rates for its increasing effect on stacking fault energy. Al alloying also increases the activation energy of deformation. Two constitutive equations with an accuracy of higher than 0.99 are proposed. The constructed processing maps show the expanded instability regions for GCr15SiMoAl steel as compared with GCr15SiMo steel. This finding is consistent with the occurrence of cracking on the GCr15SiMoAl specimens, revealing that Al alloying reduces the high-temperature plasticity of the bearing steel. On the contrary, GCr15SiMoAl steel possesses smaller grain size than GCr15SiMo steel, manifesting the positive effect of Al on bearing steel. Attention should be focused on the hot working process of bearing steel with Al.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A moist Boussinesq shallow water equations set for testing atmospheric models

The shallow water equations have long been used as an initial test for numerical methods applied to atmospheric models with the test suite of Williamson et al. being used extensively for validating new schemes and assessing their accuracy. However the lack of physics forcing within this simplified framework often requires numerical techniques to be reworked when applied to fully three dimensional models. In this paper a novel two-dimensional shallow water equations system that retains moist processes is derived. This system is derived from three-dimensional Boussinesq approximation of the hydrostatic Euler equations where, unlike the classical shallow water set, we allow the density to vary slightly with temperature. This results in extra (or buoyancy) terms for the momentum equations, through which a two-way moist-physics dynamics feedback is achieved. The temperature and moisture variables are advected as separate tracers with sources that interact with the mean-flow through a simplified yet realistic bulk moist-thermodynamic phase-change model. This moist shallow water system provides a unique tool to assess the usually complex and highly non-linear dynamics–physics interactions in atmospheric models in a simple yet realistic way. The full non-linear shallow water equations are solved numerically on several case studies and the results suggest quite realistic interaction between the dynamics and physics and in particular the generation of cloud and rain. - Highlights: • Novel shallow water equations which retains moist processes are derived from the three-dimensional hydrostatic Boussinesq equations. • The new shallow water set can be seen as a more general one, where the classical equations are a special case of these equations. • This moist shallow water system naturally allows a feedback mechanism from the moist physics increments to the momentum via buoyancy. • Like full models, temperature and moistures are advected as tracers that interact through a simplified yet realistic phase-change model. • This model is a unique tool to test numerical methods for atmospheric models, and physics–dynamics coupling, in a very realistic and simple way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Dominant two-fluid magnetohydrodynamic instabilities in CFETR upgrade phase-I scenario in presence of perfect conducting wall

Here, the dominant two-fluid magnetohydrodynamic (MHD) stabilities of China Fusion Engineering Test Reactor upgrade phase-I baseline equilibria in presence of perfect conducting wall have been evaluated using the initial value code NIMROD. Our MHD simulations show that for equilibrium with q axis < 1, high-n edge localized modes are slightly stabilized by toroidal shear flow, and n = 1 mode is found to be unstable in the core region with the characteristic of sawtooth instability. For q axis > 2 equilibrium, all modes are found to be unstable but edge localized both in single-fluid and two-fluid MHD models. Specifically for n ≤ 13 modes, the two-fluid MHD model gives a slightly higher growth rate than the single-fluid MHD model, while for n > 13 modes, this trend becomes the opposite, which means the two-fluid MHD model is needed for high-n mode analysis. In addition, n = 1–10 modes are found to be more unstable with increasing wall position and eventually their growth rates approach values in the no-wall limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multiscale characterization of phase change materials for building thermal energy storage applications

Phase change materials (PCMs) store and release large amounts of thermal energy because of their high latent energy storage capacity. However, long-term cyclic stability, supercooling and performance-scalability are some of the major challenges for their use in building thermal energy storage (TES) applications. Here, in this study, we present a comprehensive multiscale characterization of two commercially available organic PCMs, Puretemp 18 and Puretemp 23. At the microscale, differential scanning calorimetry (DSC) was used to characterize phase change temperature, specific heat, and latent heat. At the mesoscale, a heat flow meter apparatus (HFMA), following the ASTM C1784 standard, was employed to measure the phase change temperature, specific heat, and latent heat properties. A comparative analysis of latent heat as a function of temperature was conducted by integrating the DSC and HFMA results. At the macroscale, the thermal performance and cyclic stability of the TES system was evaluated using Puretemp 23. The TES system consisted of a finned tube heat exchanger with a storage volume of 0.0189 m 3 (5 gal), which represents a compact, real-world TES solution suitable for building energy storage. The results showed consistent thermal stability of the PCM over 200 cycles, and the supercooling temperature remained within 0.2 °C, which was not detected in smaller-scale characterization methods. Additionally, the macroscale testing methodology of the PCM revealed that the TES is able to charge and discharge stored latent energy within 2 h under a temperature differential of 16.67 °C measured between the inlet water temperature and the phase transition temperature of the PCM. The proposed multiscale PCM characterization method provides a systematic basis for comparing important thermal storage properties while also investigating the scalability, reliability and integration challenges in large scale TES applications.

Latent heat↗
DOE OSTI2020

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗
DOE OSTI2015

Richtmyer-Meshkov evolution under steady shock conditions in the high-energy-density regime

This work presents direct experimental evidence of long-predicted nonlinear aspects of the Richtmyer-Meshkov (RM) process, in which new modes first arise from the coupling of initially-present modes, and in which shorter-wavelength modes are eventually overtaken by longer-wavelength modes. This is accomplished using a technique we developed employing a long driving laser pulse to create a strong (Mach ~ 8) shock across a well-characterized material interface seeded by a two-mode sinusoidal perturbation. Furthermore, this technique further permits the shock to be sustained, without decay of the high-energy-density flow conditions, long enough for the system to evolve into the nonlinear phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Cross-Plane Acquisitions in Electrical Capacitance Volume Tomography

Electrical capacitance tomography (ECT) is a widely used imaging modality to image two-dimensional cross sections of multiphase flows. Recent developments in electrical capacitance volume tomography (ECVT) have made it possible to directly obtain volumetric images from measured data. An ECVT is instrumental for obtaining accurate phase hold up information and velocity information that are needed for the optimization of certain flow processes. However, compared to ECT, the high correlation between the measurements in ECVT exacerbates the ill-conditioning of the associated image reconstruction problem. Previous studies have suggested that neglecting mutual capacitance data between the ECVT electrodes located at cross-planes that are well separated along the sensor axis can be done without significantly affecting the reconstructed image. In addition, this may help constrain the ill-conditioning of the reconstruction problem. Here, we examine in detail and quantify the effect of reduced cross-plane acquisition strategies for optimizing the image reconstruction and constraining the ill-conditioning of typical ECVT settings.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2019

Asymmetries in adaptive optics point spread functions

An explanation for the origin of asymmetry along the preferential axis of the point spread function (PSF) of an AO system is developed. When phase errors from high-altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low toward the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
DOE OSTI2017

Generation Mechanisms of Quasi-parallel and Quasi-circular Flare Ribbons in a Confined Flare

We analyze a confined multiple-ribbon M2.1 flare (SOL2015-01-29T11:42) that originated from a fan-spine coronal magnetic field configuration, within active region NOAA 12268. The observed ribbons form in two steps. First, two primary ribbons form at the main flare site, followed by the formation of secondary ribbons at remote locations. We observe a number of plasma flows at extreme-ultraviolet temperatures during the early phase of the flare (as early as 15 minutes before the onset) propagating toward the formation site of the secondary ribbons. The secondary ribbon formation is co-temporal with the arrival of the pre-flare generated plasma flows. The primary ribbons are co-spatial with Ramaty High Energy Spectroscopic Imager ( RHESSI ) hard X-ray sources, whereas no enhanced X-ray emission is detected at the secondary ribbon sites. The (E)UV emission, associated with the secondary ribbons, peaks ∼1 minute after the last RHESSI hard X-ray enhancement. A nonlinear force-free model of the coronal magnetic field reveals that the secondary flare ribbons are not directly connected to the primary ribbons, but to regions nearby. Detailed analysis suggests that the secondary brightenings are produced due to dissipation of kinetic energy of the plasma flows (heating due to compression), and not due to non-thermal particles accelerated by magnetic reconnection, as is the case for the primary ribbons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Electrical impedance tomography spectroscopy method for characterising particles in solid-liquid phase

Electrical impedance tomography (EIT) is one of the process tomography techniques to provide an on-line non-invasive imaging for multiphase flow measurement. With EIT measurements, the images of impedance real part, impedance imaginary part, phase angle, and magnitude can be obtained. However, most of the applications of EIT in the process industries rely on the conductivity difference between two phases in fluids to obtain the concentration profiles. It is not common to use the imaginary part or phase angle due to the dominant change in conductivity or complication in the use of other impedance information. In a solid-liquid two phases system involving nano- or submicro-particles, characterisation of particles (e.g. particle size and concentration) have to rely on the measurement of impedance phase angle or imaginary part. Particles in a solution usually have an electrical double layer associated with their surfaces and can form an induced electrical dipole moment due to the polarization of the electrical double layer under the influence of an alternating electric field. Similar to EIT, electrical impedance spectroscopy (EIS) measurement can record the electrical impedance data, including impedance real part, imaginary part and phase angle (θ), which are caused by the polarization of the electrical double layer. These impedance data are related to the particle characteristics e.g. particle size, particle and ionic concentrations in the aqueous medium, therefore EIS method provides a capability for characterising the particles in suspensions. Electrical impedance tomography based on EIS measurement or namely, electrical impedance tomography spectroscopy (EITS) could image the spatial distribution of particle characteristics. In this paper, a new method, including test set-up and data analysis, for characterisation of particles in suspensions are developed through the experimental approach. The experimental results on tomographic imaging of colloidal particles based on EIS measurement using a sensor of 8 electrodes are reported. Results have demonstrated the potential as well as revealed the challenge in the use of EIS and EITS for characterisation of particle in suspension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

Performance of a Cross-Flow Humidifier with a High Flux Water Vapor Transport Membrane

Water vapor transport (WVT) flux across a composite membrane that consists of a very thin perfluorosulfonic acid (PFSA) ionomer layer sandwiched between two expanded polytetrafluoroethylene (PTFE) microporous layers is investigated. Static and dynamic tests are conducted to measure WVT flux for different composite structures; a transport model shows that the underlying individual resistances for water diffusion in the gas phase and microporous and ionomer layers and for interfacial kinetics of water uptake at the ionomer surface are equally important under different conditions. A finite-difference model is formulated to determine water transport in a full-scale (2-m2 active membrane area) planar cross-flow humidifier module assembled using pleats of the optimized composite membrane. In agreement with the experimental data, the modeled WVT flux in the module increases at higher inlet relative humidity (RH) of the wet stream and at lower pressures, but the mass transfer effectiveness is higher at higher pressures. The model indicates that the WVT flux is highest under conditions that maintain the wet stream at close to 100% RH while preventing the dry stream from becoming saturated. The overall water transport is determined by the gradient in RH of the wet and dry streams but is also affected by vapor diffusion in the gas layer and the microporous layer.

combined thermal and mass transfer↗
DOE OSTI2020

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Self-diffusion in the non-Newtonian regime of shearing liquid crystal model systems based on the Gay-Berne potential

The self-diffusion coefficients of nematic phases of various model systems consisting of regular convex calamitic and discotic ellipsoids and non-convex bodies such as bent-core molecules and soft ellipsoid strings have been obtained as functions of the shear rate in a shear flow. Then the self-diffusion coefficient is a second rank tensor with three different diagonal components and two off-diagonal components. These coefficients were found to be determined by a combination of two mechanisms, which previously have been found to govern the self-diffusion of shearing isotropic liquids, namely, (i) shear alignment enhancing the diffusion in the direction parallel to the streamlines and hindering the diffusion in the perpendicular directions and (ii) the distortion of the shell structure in the liquid whereby a molecule more readily can escape from a surrounding shell of nearest neighbors, so that the mobility increases in every direction. Thus, the diffusion parallel to the streamlines always increases with the shear rate since these mechanisms cooperate in this direction. In the perpendicular directions, these mechanisms counteract each other so that the behaviour becomes less regular. In the case of the nematic phases of the calamitic and discotic ellipsoids and of the bent core molecules, mechanism (ii) prevails so that the diffusion coefficients increase. However, the diffusion coefficients of the soft ellipsoid strings decrease in the direction of the velocity gradient because the broadsides of these molecules are oriented perpendicularly to this direction due the shear alignment (i). The cross coupling coefficient relating a gradient of tracer particles in the direction of the velocity gradient and their flow in the direction of the streamlines is negative and rather large, whereas the other coupling coefficient relating a gradient in the direction of the streamlines and a flow in the direction of the velocity gradient is very small.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Parametric study of flow patterns behind the standing accretion shock wave for core-collapse supernovae

In this study, we conduct three-dimensional hydrodynamic simulations systematically to investigate the flow patterns behind the accretion shock waves that are commonly formed in the post-bounce phase of core-collapse supernovae. Adding small perturbations to spherically symmetric, steady, shocked accretion flows, we compute the subsequent evolutions to find what flow pattern emerges as a consequence of hydrodynamical instabilities such as convection and standing accretion shock instability for different neutrino luminosities and mass accretion rates. Depending on these two controlling parameters, various flow patterns are indeed realized. We classify them into three basic patterns and two intermediate ones; the former includes sloshing motion (SL), spiral motion (SP), and multiple buoyant bubble formation (BB); the latter consists of spiral motion with buoyant-bubble formation (SPB) and spiral motion with pulsationally changing rotational velocities (SPP). Although the post-shock flow is highly chaotic, there is a clear trend in the pattern realization. The sloshing and spiral motions tend to be dominant for high accretion rates and low neutrino luminosities, and multiple buoyant bubbles prevail for low accretion rates and high neutrino luminosities. It is interesting that the dominant pattern is not always identical between the semi-nonlinear and nonlinear phases near the critical luminosity; the intermediate cases are realized in the latter case. Running several simulations with different random perturbations, we confirm that the realization of flow pattern is robust in most cases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Two-current spin-dependent conduction in polycrystalline LaMnO{sub 3} produced under oxygen gas flow

We have studied two-current spin-dependent conduction in polycrystalline LaMnO{sub 3} (LMO) produced under the oxygen gas flow (OGF). The polycrystalline La{sub 1−x}Sr{sub x}MnO{sub 3} (LSMO) samples were prepared with x = 0, 0.1, 0.125, 0.15, 0.175, 0.2 by use of a solid-state-reaction method. The LMO sample produced under the OGF showed large electrical conduction due to the self-hole-doping caused by the excess oxygen ions in it. The electrical resistivity and magnetoresistance (MR) of the LSMO samples were measured as a function of temperature (4K-300K). With increasing the temperature, we have observed that the MR ratios of the LSMO samples with x ≠ 0 have one maximum, while that of LMO sample has two maxima. The temperature at the maximum of the MR ratio corresponds to the magnetic phase transition temperature. The existence of two MR maxima for the LMO sample is considered to imply that the LMO sample has two regions with different doping levels; the one is the crystalline grain region with low doping and the other is the amorphous-like grain-boundary region with high doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗