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 1,657 records · Page 92

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↗

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↗

Arcjet load characteristics

Experiments were conducted to define the interface characteristics and constraints of 1 kW class arcjets run on simulated decomposition products of hydrazine and power processors. The impacts of power supply output current ripple on arcjet performance were assessed by variation of the ripple frequency from 100 Hz to 100 kHz with 10 percent peak-to-peak ripple amplitude at 1.2 kW. Ripple had no significant effects on thrust, specific impulse or efficiency. The impact of output ripple on thruster lifetime was not assessed. The static and dynamic impedances of the arcjet were quantified with two thrusters of nearly identical configuration. Superposition of an AC component on the DC arc current was used to characterize the dynamic impedance as a function of flow rate and DC current level. A mathematical model was formulated from these data. Both the static and dynamic impedance magnitude were found to be dependent on mass flow rate. The amplitude of the AC component was found to have little effect on the dynamic impedance. Reducing the DC level from 10 to 8 amps led to a large change in the magnitude of the dynamic impedance with no observable phase change. The impedance data compared favorably between the two thrusters.

Hamley, John A.↗

Arcjet load characteristics

Experiments were conducted to define the interface characteristics and constraints of 1 kW class arcjets run on simulated decomposition products of hydrazine and power processors. The impacts of power supply output current ripple on arcjet performance were assessed by variation of the ripple frequency from 100 Hz to 100 kHz with 10 percent peak-to-peak ripple amplitude at 1.2 kW. Ripple had no significant effects on thrust, specific impulse or efficiency. The impact of output ripple on thruster lifetime was not assessed. The static and dynamic impedances of the arcjet were quantitied with two thrusters of nearly identical configuration. Superposition of an AC component on the DC arc current was used to characterize the dynamic impedance as a function of flow rate and DC current level. A mathematical model was formulated from these data. Both the static and dynamic impedance magnitude were found to be dependent on mass flow rate. The amplitude of the AC component was found to have little effect on the dynamic impedance. Reducing the DC level from 10 to 8 amps led to a large change in the magnitude of the dynamic impedance with no observable phase change. The impedance data compared favorably between the two thrusters.

Hamley, John A.↗

A pictorial study of an invariant torus in phase space of four dimensions.

An investigation was conducted with the aid of a computer graphics device at Goddard Space Flight Center to study the behavior of the invariant manifolds of a particular fourth-order equation, as a parameter in the equation is varied over the interval from 0 to 1. The equation consists of two coupled Van der Pol equations. For a small parameter value, the manifold is an asymptotically stable torus, where the flow on the torus is simply a rotation. As the value of the parameter is increased, the only thing that changes is the nature of the flow on the torus, which itself persists throughout the parameter variation. It is shown that ultimately the four periodic cycles which appear play a more significant part in the phase profile of the system than does the torus itself.

Baxter, R.↗

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↗

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↗

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↗

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↗

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.↗

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↗

Automatic thermal switch

An automatic thermal switch to control heat flow includes two thermally conductive plates and a thermally conductive switch saddle pivotally mounted to the first plate. A flexible heat carrier is connected between the switch saddle and the second plate. A phase-change power unit, including a piston coupled to the switch saddle, is in thermal contact with the first thermally conductive plate. A biasing element biases the switch saddle in a predetermined position with respect to the first plate. When the phase-change power unit is actuated by an increase in heat transmitted through the first place, the piston extends and causes the switch saddle to pivot, thereby varying the thermal conduction between the two plates through the switch saddle and flexible heat carrier. The biasing element, switch saddle, and piston can be arranged to provide either a normally closed or normally opened thermally conductive path between the two plates.

Cunningham, J. W.↗

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.↗

Preliminary design studies of an advanced general aviation aircraft

The preliminary design results are presented of the advanced aircraft design project. The goal was to take a revolutionary look into the design of a general aviation aircraft. Phase 1 of the project included the preliminary design of two configurations, a pusher, and a tractor. Phase 2 included the selection of only one configuration for further study. The pusher configuration was selected on the basis of performance characteristics, cabin noise, natural laminar flow, and system layouts. The design was then iterated to achieve higher levels of performance.

Barrett, Ron↗

X-Ray Emission from Protostellar Jets

The goals of this program were to identify the termination shocks in several parsec-scale protostellar jets through the thermal soft X-rays emitted by the high-velocity terminal shock waves, and to assess the impact these outflows have on the surrounding interstellar and intercloud medium. The terminal shock waves that plow into the undisturbed interstellar medium may have shock velocities commensurate with the observed space motions of several hundred km s(exp -1). Such shocks will heat and ionize the interstellar medium, perhaps creating large ovoid bubbles of hot gas. Identifying the location of the terminal shocks in these outflows would also allow us to place constraints on the ages of these outflows, and hence the duration of the mass-loss phase that accompanies the formation of a star. We targeted four outflows where the outer optical shock waves are projected against low-extinction backgrounds. From the two dozen or so parsec-scale flows known, we chose those that have the highest velocities, brightest optical counterparts, and/or the richest concentration of shock waves in a small area on the sky. Four giant HH flows have been observed with the ROSAT HRI at Priority A for 30 ksec each. Results: Unfortunately, we failed to detect X-ray emission from the terminal bow shocks of the giant HH flows using the ROSAT HRI. The reasons are likely to be: (1) The shock velocities may be too low to emit in the soft X-rays. (2) The sensitivity of the HRI is too low to detect the diffuse emission. The fields that we observed were chosen to be low extinction sight-lines, with the consequence that the tenuous media into which the outer bow shocks are propagating produce low fluxes. Despite the failure to achieve the primary science goals, we have identified in each image a half a dozen or so point-source young stellar objects. The HRI images have been very useful for discerning potential outflow sources and are being combined with optical, near-IR, and radio imaging data of each outflow field. The first paper on L1551 NE is currently in preparation.

Morse, Jon A.↗

OASIS Observation and Analysis of Smectic Islands in Space

The Observation and Analysis of Smectic Islands in Space (OASIS) project comprises a series of experiments that will probe the interfacial and hydrodynamic behavior of freely suspended liquid crystal films in space. These are the thinnest known stable condensed phase structures, making them ideal for studies of fluctuation and interface phenomena. The experiments seek to verify theories of coarsening dynamics, hydrodynamic flow, relaxation of hydrodynamic perturbations, and hydrodynamic interactions of a near two-dimensional structure. The effects of introducing islands or droplets on a very thin bubble will be studied, both as controllable inclusions that modify the flow and as markers of flow.

liquid crystals↗

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↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for representing the phase boundary and the associated interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed in the K-site flightweight hydrogen storage tank at NASA Glenn Research Center. Laminar and turbulent simulations together with conjugated heat transfer analysis are performed. Effects of turbulence, as well as tank wall conduction are presented and discussed. Predicted tank pressures and fluid temperatures are compared with the experimental data at two different heat loads for validating the CFD model.

phase change↗