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,099 records · Page 61

Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)

Degradation of near-surface permafrost due to changes in the climate is expected to impact the hydrological, ecological and biogeochemical responses of the Arctic tundra. From a hydrological perspective, it is important to understand the movement of the various phases of water (gas, liquid and ice) during the freezing and thawing of near-surface soils. We present a new non-isothermal, single-component (water), three-phase formulation that treats air as an inactive component. This single component model works well and produces similar results to a more complete and computationally demanding two-component (air, water) formulation, and is able to reproduce results of previously published laboratory experiments. A proof-of-concept implementation in the massively parallel subsurface flow and reactive transport code PFLOTRAN is summarized, and parallel performance of that implementation is demonstrated. When water vapor diffusion is considered, a large effect on soil moisture dynamics is seen, which is due to dependence of thermal conductivity on ice content. A large three-dimensional simulation (with around 6 million degrees of freedom) of seasonal freezing and thawing is also presented.

54 ENVIRONMENTAL SCIENCES↗

An Investigation into Solution Verification for CFD-DEM

This report presents the study of the convergence behavior of the computational fluid dynamicsdiscrete element method (CFD-DEM) method, specifically National Energy Technology Laboratory’s (NETL) open source MFiX code (MFiX-DEM) with a diffusion based particle-tocontinuum filtering scheme. In particular, this study focused on determining if the numerical method had a solution in the high-resolution limit where the grid size is smaller than the particle size. To address this uncertainty, fixed particle beds of two primary configurations were studied: i) fictitious beds where the particles are seeded with a random particle generator, and ii) instantaneous snapshots from a transient simulation of an experimentally relevant problem. Both problems considered a uniform inlet boundary and a pressure outflow. The CFD grid was refined from a few particle diameters down to 1/6 th of a particle diameter. The pressure drop between two vertical elevations, averaged across the bed cross-section was considered as the system response quantity of interest. A least-squares regression method was used to extrapolate the grid-dependent results to an approximate “grid-free” solution in the limit of infinite resolution. The results show that the diffusion based scheme does yield a converging solution. However, the convergence is more complicated than encountered in simpler, single-phase flow problems showing strong oscillations and, at times, oscillations superimposed on top of globally non-monotonic behavior. The challenging convergence behavior highlights the importance of using at least four grid resolutions in solution verification problems so that (over-determined) regression-based extrapolation methods may be applied to approximate the grid-free solution. The grid-free solution is very important in solution verification and VVUQ exercise in general as the difference between it and the reference solution largely determines the numerical uncertainty. By testing different randomized particle configurations of the same general problem (for the fictitious case) or different instances of freezing a transient simulation, the numerical uncertainties appeared to be on the same order of magnitude as ensemble or time averaging uncertainties. By testing different drag laws, almost all cases studied show that model form uncertainty in this one, very important closure relation was larger than the numerical uncertainty, at least with a reasonable CFD grid, roughly five particle diameters. In this study, the diffusion width (filtering length scale) was mostly set at a constant of six particle diameters. A few exploratory tests were performed to show that similar convergence behavior was observed for diffusion widths greater than approximately two particle diameters. However, this subject was not investigated in great detail because determining an appropriate filter size is really a validation question which must be determined by comparison to experimental or highly accurate numerical data. Future studies are being considered targeting solution verification of transient simulations as well as validation of the filter size with direct numerical simulation data.

97 MATHEMATICS AND COMPUTING↗

Numerical and experimental studies on effects of moisture content on combustion characteristics of simulated municipal solid wastes in a fixed bed

Highlights: • The effects of moisture content on the burning process of MSW are investigated. • A two-dimensional mathematical model was built to simulate the combustion process. • Temperature distributions, process rates, gas species were measured and simulated. • The The conversion ratio of C/CO and N/NO in MSW are inverse to moisture content. - Abstract: In order to reveal the features of the combustion process in the porous bed of a waste incinerator, a two-dimensional unsteady state model and experimental study were employed to investigate the combustion process in a fixed bed of municipal solid waste (MSW) on the combustion process in a fixed bed reactor. Conservation equations of the waste bed were implemented to describe the incineration process. The gas phase turbulence was modeled using the k–ε turbulent model and the particle phase was modeled using the kinetic theory of granular flow. The rate of moisture evaporation, devolatilization rate, and char burnout was calculated according to the waste property characters. The simulation results were then compared with experimental data for different moisture content of MSW, which shows that the incineration process of waste in the fixed bed is reasonably simulated. The simulation results of solid temperature, gas species and process rate in the bed are accordant with experimental data. Due to the high moisture content of fuel, moisture evaporation consumes a vast amount of heat, and the evaporation takes up most of the combustion time (about 2/3 of the whole combustion process). The whole bed combustion process reduces greatly as MSW moisture content increases. The experimental and simulation results provide direction for design and optimization of the fixed bed of MSW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and validation of an enhanced filtered drag model for simulating gas-solid fluidization of Geldart A particles in all flow regimes

Coarse-grid two-fluid simulation of gas-solid fluidized bed reactors based on the kinetic theory of granular flow exhibits a significant dependence on drag models, especially for Geldart A particles. Many drag models are available in the literature, which have been reported to work for different systems. This work focused on the evaluation of an enhanced filtered drag model along with other different drag models derived from different methods for three-dimensional two-fluid model simulations of gas-solid fluidized beds of Geldart A particles covering a broad range of fluidization regimes, including bubbling fluidization, turbulent fluidization, fast fluidization, and dilute phase transport regimes. Eight drag models were selected, which included five heterogeneous drag models and three homogeneous drag models. Comparison with the available experimental data demonstrates the need for modification of homogeneous drag models to account for the effect of mesoscale structures (i.e., bubbles and clusters). The enhanced filtered drag model and energy-minimization multi-scale (EMMS) drag models were found to achieve superior predictions in all fluidization regimes, while the other drag models were only capable of predicting certain fluidization regimes. Finally, the results of this work provide a guideline for choosing appropriate drag models for simulating Geldart A particles and suggestions on developing more reliable and general drag models applicable in all flow regimes.

42 ENGINEERING↗

IMAGING AND SPECTROSCOPIC DIAGNOSTICS ON THE FORMATION OF TWO MAGNETIC FLUX ROPES REVEALED BY SDO/AIA AND IRIS

Helical magnetic flux rope (MFR) is a fundamental structure of coronal mass ejections (CMEs) and has been discovered recently to exist as a sigmoidal channel structure prior to its eruption in the EUV high-temperature passbands of the Atmospheric Imaging Assembly (AIA). However, when and where the MFR is built up are still elusive. In this paper, we investigate two MFRs (MFR1 and MFR2) in detail, whose eruptions produced two energetic solar flares and CMEs on 2014 April 18 and 2014 September 10, respectively. The AIA EUV images reveal that for a long time prior to their eruption, both MFR1 and MFR2 are under formation, which is probably through magnetic reconnection between two groups of sheared arcades driven by the shearing and converging flows in the photosphere near the polarity inversion line. At the footpoints of the MFR1, the Interface Region Imaging Spectrograph Si iv, C ii, and Mg ii lines exhibit weak to moderate redshifts and a non-thermal broadening in the pre-flare phase. However, a relatively large blueshift and an extremely strong non-thermal broadening are found at the formation site of the MFR2. These spectral features consolidate the proposition that the reconnection plays an important role in the formation of MFRs. For the MFR1, the reconnection outflow may propagate along its legs, penetrating into the transition region and the chromosphere at the footpoints. For the MFR2, the reconnection probably takes place in the lower atmosphere and results in the strong blueshift and non-thermal broadening for the Mg ii, C ii, and Si iv lines.

79 ASTRONOMY AND ASTROPHYSICS↗

Survivability of Autonomous Microgrid during Overload Events

Grid-forming sources are voltage sources that draw necessary currents to meet any load changes. A load step can cause part or all of these sources to become overloaded in a microgrid. This paper presents an overload mitigation controller that addresses the two overload issues in a microgrid by actively controlling the sources’ frequency. When part of the sources in a microgrid is overloaded, the controller autonomously transfers the extra load to other sources by rapidly reducing its frequency. The frequency difference between sources during transient results in a change of phase angle, which redistributes the power flow. When all sources in a microgrid are overloaded, each source keeps dropping the frequency. Therefore, under frequency load shedding can be used to trip the non-critical loads resulting in the survival of microgrid. The advantages of these concepts are that communications between sources are not needed during transient, and the robust voltage control is maintained. Lastly, simulation and field tests from CERTS/AEP microgrid test site verify that the control strategy is effective in both purely inverter-based microgrids and inverter & generator mixed microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Isothermal Microcalorimetric Analysis of Li/CF x Battery Discharge

Lithium/carbon monofluoride (Li/CF x ) batteries generate substantial amounts of heat during discharge in part due to the large deviation of the loaded voltage from the measured and thermodynamically predicted open circuit voltage. Here, in this study, we further analyze this system by estimating the equilibrium voltage (V eq ) and its temperature dependence (dV eq /dT) over the entire discharge range. Based on these results, the ohmic and entropic heat contributions to the overall heat flow are calculated from experimental data. The ohmic heat flow is consistent with an electrochemical impedance spectroscopy model and electrode density measurements, indicating decreasing porosity of the electrode during discharge. Entropic heating increases during discharge as the cell reaction becomes increasingly entropically unfavorable. The total energy dissipated by the cell (electrical and thermal) remains similar over the entire discharge, evidence for a two-phase reaction with a constant, and rate-independent, ΔG rxn of about -465 kJ mol -1 . However, near the end of discharge (>90% DOD), ΔG rxn changes significantly indicating possible secondary reactions.

25 ENERGY STORAGE↗

Flow Boiling Pressure Drop Characteristics of Next-generation Refrigerants in a Micro‑fin Copper Tube

This paper presents experimental frictional pressure-drop data for flow boiling of R-410A, R-134a, and next-generation alternatives R-454C, R-455A, R-1234yf, and R-1234ze(E) in a horizontal micro-fin copper tube. Tests were conducted over a range of mass fluxes and evaporation temperatures to characterize refrigerant-dependent two-phase pressure-drop behavior. The results show that frictional pressure gradient increased with mass flux and vapor quality and generally increased as evaporation temperature decreased, with liquid viscosity strongly affecting the observed trends. Among the evaluated correlations, the Goto et al. (2001) model gave the best overall agreement with the measurements before optimization. Further optimization of the Kuo and Wang (1996) and Goto et al. (2001) models reduced the overall mean absolute deviation to below 15%, with the optimized Goto (2001) model providing the most consistent predictions across all six refrigerants. The results support improved pressure-drop prediction and evaporator design for next-generation refrigerants in micro-fin tubes.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational study of the shock driven instability of a multiphase particle-gas system

This paper considers the interaction of a shock wave with a multiphase particle-gas system which creates an instability somewhat similar to the Richtmyer-Meshkov instability but with a larger parameter space. Because this parameter space is large, we only present an introductory survey of the effects of many of these parameters. We highlight the effects of particle-gas coupling, incident shock strength, particle size, effective system density differences, and multiple particle relaxation time effects. We focus on dilute flows with mass loading up to 40% and do not attempt to cover all parametric combinations. Instead, we vary one parameter at a time leaving additional parametric combinations for future work. In this work, the simulations are run with the Ares code, developed at Lawrence Livermore National Laboratory, which uses a multiphase particulate transport method to model two-way momentum and energy coupling. A brief validation of these models is presented and coupling effects are explored. It is shown that even for small particles, on the order of 1μm, multi-phase coupling effects are important and diminish the circulation deposition on the interface by up to 25%. These coupling effects are shown to create large temperature deviations from the dusty gas approximation, up to 20% greater, especially at higher shock strengths. It is also found that for a multiphase instability, the vortex sheet deposited at the interface separates into two sheets. In conclusion, depending on the particle and particle-gas Atwood numbers, the instability may be suppressed or enhanced by the interactions of these two vortex sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

TEA of a Unique Two-Pathways Process for Post-Combustion CO 2 Capture

A unique two-Pathways process using aqueous sodium glycinate for CO 2 capture from a split flue gas stream emitted from 600 MWe post-combustion coal power plant was developed in Aspen Plus v.10. The split gas flow rate used was 44.75 ton/h and contained 0.0023 mol% SO 2 and 13.33 mol% CO 2 . The process includes a washing unit, a CO 2 absorption unit, a reverse osmosis unit, and a solvent regeneration unit or an ultrafiltration unit. The washing unit uses deionized water to completely remove SO 2 and the CO 2 absorption unit uses SGS to capture at least 90 mol% of the CO 2 in the split flue gas stream. Upon CO 2 and SGS reactions, the resulting liquid products exhibit phase-separation into CO 2 -lean phase and CO 2 -rich phase, allow two distinct pathways. Pathway (i) is to regenerate mostly the CO 2 -rich phase, collect the released CO 2 , and compress it for sequestration purposes. Pathway (ii) is to send the liquid stream from the CO 2 absorption unit to the ultrafiltration unit to separate the solid nanomaterials. The hydraulics and mass transfer characteristics in the washing and CO 2 absorption units were obtained; and techno-economic analysis (TEA) for Pathways (i) and (ii), including Capital Expenditure (CAPEX), Operating Expenditure (OPEX), and Levelized Cost of CO 2 Captured (LCOC), were calculated and compared. The simulation results revealed that the CAPEX, OPEX, and LCOC for Pathway (i) were ($\$12,039,251$), (261 dollar/h), and (54.01 dollar/ton-CO 2 captured), respectively, and those for Pathway (ii) were ($\$5,908,000$), (237.2 dollar/h), and (39.90 dollar/ton-CO 2 captured), respectively. Moreover, in Pathway (ii), 8.19 ton/h of CO 2 were captured to produce 15.62 ton/h NaHCO 3 nanomaterials, which were sold to offset the overall process cost. The LCOC values indicate that Pathway (ii) is more cost-effective than Pathway (i) because LCOC values for Pathway (ii) are much lower than those for Pathway (i).COC values for Pathway (ii) are much lower than those for Pathway (i).

20 FOSSIL-FUELED POWER PLANTS↗

TH-CD-207A-05: Lung Surface Deformation Vector Fields Prediction by Monitoring Respiratory Surrogate Signals

Purpose: In this study, we developed and evaluated a method for predicting lung surface deformation vector fields (SDVFs) based on surrogate signals such as chest and abdomen motion at selected locations and spirometry measurements. Methods: A Patient-specific 3D triangular surface mesh of the lung region at end-expiration (EE) phase was obtained by threshold-based segmentation method. For each patient, a spirometer recorded the flow volume changes of the lungs; and 192 selected points at a regular spacing of 2cm X 2cm matrix points over a total area of 34cm X 24cm on the surface of chest and abdomen was used to detect chest wall motions. Preprocessing techniques such as QR factorization with column pivoting (QRCP) were employed to remove redundant observations of the chest and abdominal area. To create a statistical model between the lung surface and the corresponding surrogate signals, we developed a predictive model based on canonical ridge regression (CRR). Two unique weighting vectors were selected for each vertex on the surface of the lung, and they were optimized during the training process using the all other phases of 4D-CT except the end-inspiration (EI) phase. These parameters were employed to predict the vertices locations of a testing data set, which was the EI phase of 4D-CT. Results: For ten lung cancer patients, the deformation vector field of each vertex of lung surface mesh was estimated from the external motion at selected positions on the chest wall surface plus spirometry measurements. The average estimation of 98th percentile of error was less than 1 mm (AP= 0.85, RL= 0.61, and SI= 0.82). Conclusion: The developed predictive model provides a non-invasive approach to derive lung boundary condition. Together with personalized biomechanical respiration modelling, the proposed model can be used to derive the lung tumor motion during radiation therapy accurately from non-invasive measurements.

60 APPLIED LIFE SCIENCES↗

Pulsed laminar falling films in vertical tubes: Maintaining a continuous liquid film with reduced film thickness

Falling-liquid-film heat transfer is important to a wide variety of industrial applications, including two-phase thermosyphons, falling-film evaporators in the food industry, heat exchangers in fuel-ethanol processing, and the chemical service industry. Particularly for laminar falling films, the heat transfer increases as the film thickness decreases, hence it is desirable to keep the films as thin as possible. This article presents a periodically-pulsed flow concept to reduce the average mass flow inside a vertical tube. The technique provides thinner films than can be achieved with continuous flow. A drainage model is adopted to predict the film thickness as a function of flow rate and location along the tube. The concept is confirmed experimentally by measuring the film thickness in real time in a circular tube using a laser-based absorption technique. A new heat-transfer-based film continuity measurement is also integrated into the optical film measurement setup. Water at ambient temperature and pressure is used for this work. Film thickness reductions of up to 60%, and liquid-film volume-flow reductions of up to 80% compared to the minimum continuous-flow film have been obtained. The time-varying film thickness profile can also be controlled to a limited extent, which may be of interest for process industries.

42 ENGINEERING↗

Electromagnetic physics vectorization in the GeantV transport framework

The development of the GeantV Electromagnetic (EM) physics package has evolved following two necessary paths towards code modernization. A first phase required the revision of the main electromagnetic physics models and their implementation. The main objectives were to improve their accuracy, extend them to the new high-energy frontier posed by the Future Circular Collider (FCC) programme and allow a better adaptation to a multi-particle flow. Most of the EM physics models in GeantV have been reviewed from theoretical perspective and rewritten with vector-friendly implementations, being now available in scalar mode in the alpha release. The second phase consists of a thorough investigation on the possibility to vectorise the most CPU-intensive physics code parts, such as final state sampling. We have shown the feasibility of implementing electromagnetic physics models that take advantage of SIMD/SIMT architectures, thus obtaining gains in performance. After this phase, the time has come for the GeantV project to take a step forward towards the final proof of concept. This takes shape through the testing of the full simulation chain (transport + physics + geometry) running in vectorized mode. In this paper we will present the first benchmark results obtained after vectorizing a full set of electromagnetic physics models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Economic evaluation of distribution system smart grid investments

This paper investigates economic benefits of smart grid automation investments. A system consisting of 7 substations and 14 feeders is used in the evaluation. Here benefits that can be quantified in terms of dollar savings are considered, termed “hard dollar” benefits. Smart Grid investment evaluations to be considered include investments in improved efficiency, more cost effective use of existing system capacity with automated switches, and coordinated control of capacitor banks and voltage regulators. These Smart Grid evaluations are sequentially ordered, resulting in a series of incremental hard dollar benefits. Hard dollar benefits come from improved efficiency, delaying large capital equipment investments, shortened storm restoration times, and reduced customer energy use. Analyses used in the evaluation involve hourly power flow analysis over multiple years and Monte Carlo simulations of switching operations during storms using a reconfiguration for restoration algorithm. The economic analysis uses the time varying value of the Locational Marginal Price. Algorithms used include reconfiguration for restoration involving either manual or automated switches and coordinated control involving two modes of control. Field validations of phase balancing and capacitor design results are presented. The evaluation shows that investments in automation can improve performance while at the same time lowering costs.

25 ENERGY STORAGE↗

Dynamic structural and microstructural responses of a metal–organic framework type material to carbon dioxide under dual gas flow and supercritical conditions

The structural and microstructural responses of a model metal–organic framework material, Ni(3-methyl-4,4'-bipyridine)[Ni(CN) 4 ] (Ni-BpyMe or PICNIC-21), to CO 2 adsorption and desorption are reported for in situ small-angle X-ray scattering and X-ray diffraction measurements under different gas pressure conditions for two technologically important cases. These conditions are single or dual gas flow (CO 2 with N 2 , CH 4 or H 2 at sub-critical CO 2 partial pressures and ambient temperatures) and supercritical CO 2 (with static pressures and temperatures adjusted to explore the gas, liquid and supercritical fluid regimes on the CO 2 phase diagram). The experimental results are compared with density functional theory calculations that seek to predict where CO 2 and other gas molecules are accommodated within the sorbent structure as a function of gas pressure conditions, and hence the degree of swelling and contraction in the associated structure spacings and void spaces. Furthermore, these predictions illustrate the insights that can be gained concerning how such sorbents can be designed or modified to optimize the desired gas sorption properties relevant to enhanced gas recovery or to addressing carbon dioxide reduction through carbon mitigation, or even direct air capture of CO 2 .

36 MATERIALS SCIENCE↗

Scaling, rotation, and channeling behavior of helical and skyrmion spin textures in thin films of Te-doped Cu 2 OSeO 3

Topologically nontrivial spin textures such as vortices, skyrmions, and monopoles are promising candidates as information carriers for future quantum information science. Their controlled manipulation including creation and annihilation remains an important challenge toward practical applications and further exploration of their emergent phenomena. Here, we report controlled evolution of the helical and skyrmion phases in thin films of multiferroic Te-doped Cu 2 OSeO 3 as a function of material thickness, dopant, temperature, and magnetic field using in situ Lorentz phase microscopy. We report two previously unknown phenomena in chiral spin textures in multiferroic Cu 2 OSeO 3 : anisotropic scaling and channeling with a fixed-Q state. The skyrmion channeling effectively suppresses the recently reported second skyrmion phase formation at low temperature. Our study provides a viable way toward controlled manipulation of skyrmion lattices, envisaging chirality-controlled skyrmion flow circuits and enabling precise measurement of emergent electromagnetic induction and topological Hall effects in skyrmion lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Concept Design of Active Shielding for Dynamic Wireless Charging of Light-duty EV

Dynamic wireless charging of electric vehicles is a flexible and state of the art charging technology with the potential capability of enabling fully automated in-motion charging. With charging power increasing to more than 100 kW for light duty vehicles, electromagnetic field (EMF) emission becomes a critical challenge. Due to the high costs of ferrite materials, this paper proposes an active shielding solution with multiple canceling coils installed on the ground side to supplement ferrite passive shielding to ensure electromagnetic safety. Two canceling coils are designed on two sides of the ground side coil. The canceling coils are small in size and 180 degree opposite in phase to the ground coil. Simulation and modeling shows that the canceling coils can reduce EMF emission from 37.2 µT to 18.2 µT at 0.8 m during 100 kW operation with only 2.5% of ground side current flowing in the canceling coils. These results have been preliminarily verified by inductance measurements and magnetic field measurement at 1.1 m distance without canceling coils. By further increasing the canceling coils’ currents to 4%, EMF emission can also be mitigated at 200 kW, although the canceling coils’ shape, position, and phase angle can be further optimized to improve the three-dimensional field distribution.

25 ENERGY STORAGE↗