SEARCH · Engineering Papers

Results for “Multiphase 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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Learning Distribution Grid Topologies: A Tutorial

Unveiling feeder topologies from data is of paramount importance to advance situational awareness and proper utilization of smart resources in power distribution grids. This tutorial summarizes, contrasts, and establishes useful links between recent works on topology identification and detection schemes that have been proposed for power distribution grids. The primary focus is to highlight methods that overcome the limited availability of measurement devices in distribution grids, while enhancing topology estimates using conservation laws of power-flow physics and structural properties of feeders. Grid data from phasor measurement units or smart meters can be collected either passively in the traditional way, or actively, upon actuating grid resources and measuring the feeder's voltage response. Analytical claims on feeder identifiability and detectability are reviewed under disparate meter placement scenarios. Such topology learning claims can be attained exactly or approximately so via algorithmic solutions with various levels of computational complexity, ranging from least-squares fits to convex optimization problems, and from polynomial-time searches over graphs to mixed-integer programs. Although the emphasis is on radial single-phase feeders, extensions to meshed and/or multiphase circuits are sometimes possible and discussed. Here this tutorial aspires to provide researchers and engineers with knowledge of the current state-of-the-art in tractable distribution grid learning and insights into future directions of work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗

NASA NTRS1997

Evidence for a Multiphase ISM in Early Type Galaxies and Elliptical Galaxies with Strong Radio Continuum

We have observed NGC 1316 (Fornax A) with the ROSAT HRI. In this paper, we present the results of these observations and we complement them with the spectral analysis of the archival PSPC data. The spectral properties suggest the presence of a significant component of thermal X-ray emission (greater than 60%), amounting to approx. 10(exp 9) solar mass of hot ISM. Within 3 feet from the nucleus of NGC 1316, the HRI X-ray surface brightness falls as r(exp -2) following the stellar light. In the inner approx. 30 inch., however, the X-ray surface brightness is significantly elongated, contrary to the distribution of stellar light, which is significantly rounder within 10 inch. This again argues for a non-stellar origin of the X-ray emission. This flattened X-ray feature is suggestive of either the disk-like geometry of a rotating cooling flow and/or the presence of extended, elongated dark matter. By comparing the morphology of the X-ray emission with the distribution of optical dust patches, we find that the X-ray emission is significantly reduced at the locations where the dust patches are more pronounced, indicating that at least some of the X-ray photons are absorbed by the cold ISM. We also compare the distribution of the hot and cold ISM with that of the ionized gas, using recently obtained H(sub alpha) CCD data. We find that the ionized gas is distributed roughly along the dust patches and follows the large scale X-ray distribution at r greater than 1 foot from the nucleus. However, there is no one-to-one correspondence between ionized gas and hot gas. Both morphological relations and kinematics suggest different origins for hot and cold ISM. The radio jets in projection appear to pass perpendicularly through the central X-ray ellipsoid. Comparison of thermal and radio pressures suggests that the radio jets are confined by the surrounding hot gaseous medium.

Kim, Dong Woo↗

NASA NTRSDate not supplied

Parabolic Flights in Support of the Ring-Sheared Drop (RSD) Module Aboard the ISS

The Ring-Sheared Drop module was conceived as a centimeter-scale, containerless biochemical reactor to study shearing flow at gas-liquid interfaces without complications associated with solid walls. The drop is constrained by a thin, stationary contact ring in one hemisphere and sheared by the steady rotation of another ring in the other hemisphere. The original experiment considers a 2.5 cm diameter drop of human insulin solution in order to examine the formation of amyloid fibrils without wall nucleation. Shear-induced amyloidogenesis is of wide interest due to the pathogenic role amyloid proteins play in diseases such as Alzheimer's, Parkinson's, prion diseases, and type-2 diabetes. Pinning full-scale drops of native protein solutions - prior to any fibrillization - was challenging and required a collaborative effort. In this talk we report on laboratory experiments, multiphase CFD simulations, and several parabolic flights to find strategies that are effective for constraining large drops of protein in microgravity.

Amir Hirsa↗

DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Turbulence kinetic energy equation for dilute suspensions

A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Abou-Arab, T. W.↗

NASA NTRS1995

Microgravity combustion of dust clouds: Quenching distance measurements

The current level of physical understanding of dust combustion phenomena is still in a rudimentary state compared with the understanding of gas combustion processes. The reason for such a lack of fundamental understanding is partially based on the complexity of multiphase combustion and the enormous diversity of chemical-physical properties of heterogeneous combustible mixtures but is largely due to difficulties in the experimental investigation of dust combustion. The influence of gravity on a dust suspension is the main reason. First of all, when particulates (either solid particles or liquid droplets) with a characteristic size of the order of tens of microns are suspended, they rapidly settle in the gravitational field. To maintain a particulate suspension for a time duration adequate to carry out combustion experiments invariably requires continuous convection of particulates at or in excess of the gravitational settling velocity. Of necessity, this makes the experiments turbulent in character and makes it impossible to study laminar dust flames. For particle sizes of the order of microns a stable laminar dust flow can be maintained only for relatively small dust concentrations (e.g., for low fuel equivalence ratios) at normal gravity conditions. High dust loading leads to gravitational instability of the dust cloud and to the formation of recirculation cells in a dust suspension in a confined volume, or to the rapid sedimentation of the dense dust cloud as a whole in an unconfined volume. In addition, many important solid fuels such as low volatile coal, carbon, and boron have low laminar flame speeds (of the order of several centimeters per second). Gravitational convection that occurs in combustion products due to the buoyancy forces disrupts low speed dust flames and, therefore, makes observation of such flames at normal gravity impossible. The only way to carry out 'clean' fundamental experiments in dust combustion over a wide range of dust cloud parameters is in a gravity-free environment. Access to the microgravity environment provided by the use of large-scale drop towers, parabolic flights of aircraft and rockets, and shuttle and space station orbits has permitted now to proceed with a systematic program of dust combustion microgravity research. For example, the NASA-Lewis drop tower and a Lear jet parabolic flight aircraft were used by Ross et al. and by Berlad and Tangirala for experiments with Iycopodium/air mixtures. The Japan Microgravity Center drop shaft (JAMIC) where a microgravity condition of 10(exp -4) g for 10 s is available, was recently used by Kobayashi, Niioka et al. for measuring flame propagation velocities in polymethyl methacrylate dust/air suspensions. Microgravity dust combustion experiments were started at McGill University in the early 90's under the sponsorship of the Canadian Space Agency. Several generations of dust combustion platforms permitting dust combustion microgravity experiments to be carried out on board a parabolic flight aircraft (KC-135, NASA) have been designed and tested. The experimental data and experience gained from this research allowed us to design and build in a current phase of this program the microgravity apparatus for the visual observation of freely propagating constant pressure laminar dust flames. Quenching distances in aluminum dust suspensions have been measured in a wide range of dust cloud parameters in ground-based experiments and in recent microgravity experiments (KC-135 parabolic flights, Houston, February 1995).

Goroshin, S.↗

DOE OSTI2025

Scaling microstructural processes in the sintering of ionic ceramics

A multi-scale framework, combining a multiphase field formulation and large deformation mechanics, was developed as a stepping stone to perform the data analytics of the microstructural level kinetics of a sintering solid. Relevant microstructural information from this framework, such as grain, stress, and porosity statistics, was scaled up to describe the macroscopic level sintering kinetics. Here, the developed formulation was applied to describe the electric field assisted sintering of Y 2 O 3 . Microstructural inhomogeneities in a multi-granular solid result in the formation of a field of compressive stress networks, which interleave with low compression and weakly tensile regions, defining a scaffolding for sintering concentration regions to develop. A Poisson effect-induced lateral stress network is also naturally self-induced as a result of the mechanical constraints imposed by the sintering apparatus. For long sintering times, localized shear stresses enhancing mass flow along grain boundaries and internal surfaces develop. Three-sided pores are removed by either vacancy transport to the surrounding pores, or move towards the external surfaces through grain boundary diffusion. Four- and higher order-sided pores stabilize because an equal amount of vacancies are gained and lost through the connecting grain boundaries. Grain dewetting contributes to pore coalescence, suggesting that pore kinetics and grain growth are coupled and should be analyzed in concert. The combined sintering and grain growth kinetics define six regimes of sintering behavior: (1) T, the transient regime; (2) E$_Υ$, the surface energy dominated, early sintering regime, where the grain growth exponent, p = 1, and the stress concentration factor, $f$ ~ $1/\hat{ρ}^{4.6}$; (3) E S , the stress dominated, early sintering regime, where p = 1 and $f$ ~ $1/\hat{ρ}^{4}$; (4) I$_Υ$, the surface energy dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4.6}$; (5) I S , the stress dominated, intermediate sintering regime, where p = 2 and $f$ ~ $1/\hat{ρ}^{4}$; and (6) L, the late sintering regime, where p = 3 and $f$ ~ 1. At the macroscopic level, the rapid densification and suppression of grain growth observed in the electric field assisted sintering process is a consequence of the compounding effects of the underlying stress-, transport-, and interfacial-energy-induced energy minimization kinetics, as predicted by the multi-scale framework.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2018

Columbus IFHX Ammonia Leak Analysis

After the Columbus Moderate Temperature Loop (MTL) InterFace Heat eXchanger (IFHX) low temperature event of GMT 345-2013, NASA investigated relevant transient scenarios involving IFHX rupture after water freezing and subsequent thawing. NASA recommended development of a Fault Detection Isolation and Recovery (FDIR) plan that would, in the event of a heat exchanger freeze event, close the Water On/Off Valves (WOOVs) to isolate the heat exchanger and prevent ammonia from the external flow loops from spreading into the cabin. NASA performed a preliminary simplified analysis for the reference case of IFHX rupture, but for a deeper understanding TAS developed detailed SINDA-FLUINT models of the Columbus ITCS that were built and run through the SINAPS GUI. This allowed simulation of the ammonia leakage physics including the variation of environmental parameters, thus providing more accurate and specific input to the FDIR under development. The result was finalization of the IFHX WOOVs closure sequence and wait times to contain the ammonia propagation to Columbus and allow identification of the leaking IFHX. In addition, the analysis results provided reference pressure profiles to be used on console and by the Engineering as support for the telemetry data assessment in case of failure.This paper gives an overview on the issue and focuses on the analytical aspects of the multiphase fluid dynamics involved.

Tilloca, A.↗

DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Description of a 20 kilohertz power distribution system

A single phase, 440 VRMS, 20 kHz power distribution system with a regulated sinusoidal wave form is discussed. A single phase power system minimizes the wiring, sensing, and control complexities required in a multi-sourced redundantly distributed power system. The single phase addresses only the distribution links multiphase lower frequency inputs and outputs accommodation techniques are described. While the 440 V operating potential was initially selected for aircraft operating below 50,000 ft, this potential also appears suitable for space power systems. This voltage choice recognizes a reasonable upper limit for semiconductor ratings, yet will direct synthesis of 220 V, 3 power. A 20 kHz operating frequency was selected to be above the range of audibility, minimize the weight of reactive components, yet allow the construction of single power stages of 25 to 30 kW. The regulated sinusoidal distribution system has several advantages. With a regulated voltage, most ac/dc conversions involve rather simple transformer rectifier applications. A sinusoidal distribution system, when used in conjunction with zero crossing switching, represents a minimal source of EMI. The present state of 20 kHz power technology includes computer controls of voltage and/or frequency, low inductance cable, current limiting circuit protection, bi-directional power flow, and motor/generator operating using standard induction machines. A status update and description of each of these items and their significance is presented.

Hansen, I. G.↗

DOE OSTI2021

Model investigation of the longitudinal broadening of the transverse momentum two-particle correlator

Here, the multiphase transport model is used to investigate the longitudinal broadening of the transverse momentum two-particle correlator C 2 (Δη,Δφ), and its utility to extract the specific shear viscosity, η/s, of the quark-gluon plasma formed in ultrarelativistic heavy ion collisions. The results from these model studies indicate that the longitudinal broadening of C 2 (Δη,Δφ) is sensitive to the value of η/s. However, reliable extraction of the longitudinal broadening of the correlator requires the suppression of possible self-correlations associated with the definition of the collision centrality.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Characterizing initial- and final-state effects of relativistic nuclear collisions

Here, the multiphase transport model (AMPT) is used to study the final-state effects on the symmetric correlations (SC), asymmetric correlations (ASC), normalized symmetric correlations (NSC), and normalized asymmetric correlations (NASC) in Au+Au collisions at 200 GeV. The correlators' sensitivity to nonflow effects associated with long- and short-range nonflow correlations are also discussed using the HIJING model. The results indicate that SC, ASC, NSC, and NASC can give accompanying constraints for initial- and final-state effects. In addition, conducting further detailed experimental measurements spanning a broad range of collision systems and beam energies will serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Comprehensive Digital Twin of a Macro-fluidic Electrochemical Reactor to Optimize the Electrochemical-based Recovery of O2 from Metabolic CO2

Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery from metabolic CO 2 to achieve self-sufficiency for long space missions beyond Earth's low orbit. A Macro-fluidic Electrochemical Reactor (MFECR) development effort to electrochemically recover O 2 from CO 2 is underway at NASA Marshall Space Flight Center (MSFC) to increase current O 2 recovery efficiency and reduce air revitalization (AR) system complexity at the International Space Station (ISS) habitat and future long missions. The authors have developed and deployed a digital twin of an actual single cell of the MFECR via a comprehensive 3D multiphysics model that thoroughly replicates the exact configuration and fluid/material domains of the MFECR. This model's electrochemical physics consists of multicomponent-multiphase electrochemical-driven reactions leading to CO 2 conversion to C 2 H 4 and CO along with the formation of H 2 on the cathode in parallel with the generation of O 2 and H 2 O on the anode. This electrochemical model is coupled with all the physics phenomena involved in the process, including but not limited to fluid and non-ideal mass transfer of reactant and product species in free and porous media, convective/conduction/radiative heat transfer, and conduction of DC electrical current with Joule heating generation. The digital twin has proved to be an essential tool for performing different qualitative studies, including the effect of reducing the height (increasing the respect ratio) of the serpentine walls, leading to a further redesign of the EDU (the newly redesigned MFECR EDU has been fabricated and is expected to be used in future tests), the evaluation of setting different MFECRs connected in series, and the assessment of feeding air directly to the MFECR skipping the preceding metabolic CO 2 separation from air. The MFECR's test stand is fully automated and equipped with several inline measurements (flow, pressure, temperature, pH, component concentration) systems on all six MFECR's IO streams, allowing reliable experimental validation of the model, parametric determination of all electrochemical reactions, and process optimization.

Jesus A Dominguez↗

NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Computational Fluid Dynamics↗