Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flows”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hydrodynamic irreversibility of non-Brownian suspensions in highly confined duct flow

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

42 ENGINEERING↗
DOE OSTI2019

Numerical Simulation of the Three-Phase Flow of a Bubble Interacting with the Steel–Slag Interface During the Secondary Refining Process

Gas stirring is an effective method for the secondary refining process of steel. The microphenomena among the phases of an argon bubble floating in steel and motion across the steel–slag interface were simulated and discussed with a sharp interface tracking method. A water-model experiment was set to validate the numerical model, and the experimental results showed agreement with the numerical model. A critical size for bubbles to pass the interface was found using the model.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Computational study of the bubbling-to-slugging transition in a laboratory-scale fluidized bed

Here, we report results from a computational study of the transition from bubbling to slugging in a laboratory-scale fluidized-bed reactor with Geldart Group B glass particles. For simulating the three-dimensional fluidized-bed hydrodynamics, we employ MFiX, a widely studied multi-phase flow simulation tool, that uses a two-fluid Eulerian-Eulerian approximation of the particle and gas dynamics over a range of gas flows. We also utilize a previously published algorithm to generate bubble statistics that can be correlated with pressure fluctuations to reveal previously unreported details about the stages through which the hydrodynamics progress during the bubbling-to-slugging transition. We expect this new information will lead to improved approaches for on-line reactor diagnostics, as well as new approaches for validating the results of computational fluidized-bed simulations with experimental measurements.

42 ENGINEERING↗
DOE OSTI2017

Clustering instabilities in sedimenting fluid–solid systems: critical assessment of kinetic-theory-based predictions using direct numerical simulation data

In this paper the quantitative and qualitative ability of a kinetic-theory-based two-fluid model (KT-TFM) is assessed in a state of fully periodic sedimentation (fluidization), with a focus on statistically steady, unstable (clustered) states. The accuracy of KT-TFM predictions is evaluated via direct comparison to direct numerical simulation (DNS) data. The KT-TFM and DNS results span a rather wide parameter space: mean-flow Reynolds numbers on the order of 1 and 10, mean solid volume fractions from 0.1 to 0.4, solid-to-fluid density ratios from 10 to 1000 and elastic and moderately inelastic (restitution coefficient of 0.9) conditions. Data from both KT-TFM and DNS display a rich variety of statistically steady yet unstable structures (clusters). Instantaneous snapshots of KT-TFM and DNS demonstrate remarkable qualitative agreement. This qualitative agreement is quantified by calculating the critical density ratio at which the structure transitions from a chaotic, dynamic state to a regular, plug-flow state, with good overall comparisons. Further quantitative assessments of mean and fluctuating velocities show good agreement at high density ratios but weaker agreement at intermediate to low density ratios depending on the mean-flow Reynolds numbers and solid fractions. Deviations of the KT-TFM results from the DNS data were traced to a breakdown in one of the underlying assumptions of the kinetic theory derivation: high thermal Stokes number. Surprisingly, however, even though the low Knudsen number assumption, also associated with the kinetic theory derivation, is violated throughout most of the parameter space, it does not seem to affect the good quantitative accuracy of KT-TFM simulations.

42 ENGINEERING↗
DOE OSTI2017

Catalyst Residence Time Distributions in Riser Reactors for Catalytic Fast Pyrolysis. Part 2: Pilot-Scale Simulations and Operational Parameter Study

Here, wsing the validated simulation model developed in part one of this study for biomass catalytic fast pyrolysis (CFP), we assess the functional utility of using this validated model to assist in the development of CFP processes in fluidized catalytic cracking (FCC) reactors to a commercially viable state. Specifically, we examine the effects of mass flow rates, boundary conditions (BCs), pyrolysis vapor molecular weight variation, and the impact of the chemical cracking kinetics on the catalyst residence times. The factors that had the largest impact on the catalyst residence time included the feed stock molecular weight and the degree of chemical cracking as controlled by the catalyst activity. Lastly, because FCC reactors have primarily been developed and utilized for petroleum cracking, we perform a comparison analysis of CFP with petroleum and show the operating regimes are fundamentally different.

09 BIOMASS FUELS↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Potential seal bypass and caprock storage produced by deformation-band-to-opening-mode-fracture transition at the reservoir/caprock interface

Here, we examined the potential impact on CO 2 transport of zones of deformation bands in reservoir rock that transition to opening-mode fractures within overlying caprock. Sedimentological and petrophysical measurements were collected along an approximately 5 m × 5 m outcrop of the Slick Rock and Earthy Members of the Entrada Sandstone on the eastern flank of the San Rafael Swell, Utah, USA. Measured deformation band permeability (2 mD) within the reservoir facies is about three orders of magnitude lower than the host sandstone. Average permeability of the caprock facies (0.0005 mD) is about seven orders of magnitude lower than the host sandstone. Aperture-based permeability estimates of the opening-mode caprock fractures are high (3.3 × 10 7 mD). High-resolution CO 2 –H 2 O transport models incorporate these permeability data at the millimeter scale. We then varied fault properties at the reservoir/caprock interface between open fractures and deformation bands as part of a sensitivity study. Numerical modeling results suggest that zones of deformation bands within the reservoir strongly compartmentalize reservoir pressures largely blocking lateral, cross-fault flow of supercritical CO 2 . Significant vertical CO 2 transport into the caprock occurred in some scenarios along opening-mode fractures. The magnitude of this vertical CO 2 transport depends on the small-scale geometry of the contact between the opening-mode fracture and the zone of deformation bands, as well as the degree to which fractures penetrate caprock. Finally, the presence of relatively permeable units within the caprock allows storage of significant volumes of CO 2 , particularly when the fracture network does not extend all the way through the caprock.

58 GEOSCIENCES↗
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 OSTI2017

Dynamic reduced-order models of integrated physics-specific systems for carbon sequestration

In this study, CO 2 leakage from a deep storage reservoir into a shallow groundwater aquifer through a wellbore is simulated using physics-specific models conceptualized as nonisothermal two-phase flow in a deep storage reservoir and through the wellbore, and reactive transport in the aquifer. Supercritical CO 2 is injected into the reservoir subsystem with uncertain permeabilities of reservoir and injection strategies (rate and duration as decision variables). The simulated pressure and CO 2 saturation at an uncertain location are connected to the wellbore model as a boundary condition to output CO 2 flux to the reactive transport model in the shallow aquifer. Uncertainties are propagated from the deep reservoir model, to the wellbore leakage model, and eventually to the reactive transport model, thus contributing to system performance. Finally, to minimize the risk associated with the CO 2 leakage, we develop dynamic reduced order models for quantifying and managing the uncertainty propagation.

58 GEOSCIENCES↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Experimental Investigation of the Effect of Atomizing Electrospray Nozzles on the Cooling and Reduced Humidification of Air

Separating sensible and latent cooling have shown the potential to reduce energy consumption in A/C systems for building HVAC applications. Several technologies exist based on vapor-compression refrigeration, enthalpy wheels, chemical adsorption and absorption materials, and mechanical cooling. However, their thermodynamics limits and high energy consumption hinder their deployment in hot and humid climatic regions. This paper focuses on separating sensible and latent cooling by using a new thermodynamics process and in-kind (i.e., non-vapor compression-based) technology. Highly electrically charged water droplets were sprayed in the airflow. These droplets attracted water vapor molecules to their surfaces and promoted condensation. The phenomenon was the result of simultaneous dielectrophoresis and electro-diffusion interactions. Studies in the literature have shown that using multiple capillary electrodes reduced air moisture by up to 5% when using nanometer-size droplets in the spray. Unfortunately, these studies were limited to low airflow rates, and the objective of this paper was to investigate how to scale up this in-kind approach to airflows typical of buildings. In the present paper, droplets of micrometer size were utilized to control the humidity for a 5-cfm flow rate. While this airflow was still low for building applications, it was 100 times fold the airflows in the literature studies. The air was tested at 20°C and at 50 and 80% relative humidity. A two-fluid atomizing nozzle produced the droplets in the spray, and high DC electric potential, up to 25 kV, was used to charge the fine droplets electrically. The air atomizing nozzle with high voltage potential resembled an evaporative cooler process. However, a measurable reduction of the absolute humidity of up to 2% was observed compared to the case of the nozzle with no high voltage potential. The entire device had one small nozzle selected from off-the-shelf components and had less than a 9 cm2 footprint area.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Development of Multiple Solid Components Capability in GOTHIC Droplet Fields

GOTHIC is a versatile and generally applicable software package that solves complex thermal hydraulics problems. GOTHIC solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and/or multi-dimensional geometries. GOTHIC version 8.2 includes a new feature to model one or more solid components in each of the modeled droplet fields. A solid component is a material present in the water which has volume, density, and heat capacity. GOTHIC has supported modeling multiple droplet fields since version 8.0 - each field has its own mass, energy, and momentum equation. Recent work by Lane et. al. allowed multiple solid components to be tracked in the continuous liquid field. With the release of version 8.2 the droplet fields can now include one or more solid components within each droplet field. The addition of solid components to the droplet field allows faithful simulation of one or more solid components throughout the liquid and droplet fields. Preliminary expected applications are for scenarios such as GSI-191 and longer term desired applications are aerosol simulation. (authors)

42 ENGINEERING↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2024

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗
DOE OSTI2016

Enhanced Performance Assessment Models for Generic Deep Geologic Repositories for High-Level Waste and Spent Nuclear Fuel - 16223

Development of an enhanced performance assessment (PA) capability for geologic disposal of spent nuclear fuel and high-level waste has been ongoing for several years in the U.S. repository program. The new Generic Disposal System Analysis (GDSA) modeling and software framework is intended to be flexible enough to evolve through the various phases of repository activities, beginning with generic PA activities in the current Concept Evaluation phase to site-specific PA modeling in the Repository Development phase. The GDSA Framework utilizes modern software and hardware capabilities by being based on open-source software architecture and being configured to run in a massively parallel, high-performance computing (HPC) environment. It consists of two main components, the open-source Dakota uncertainty sampling and analysis software and the PFLOTRAN reactive multi-phase flow and transport simulator. Reference cases or 'generic repositories' have been, and are being developed, based on typical properties for potential salt, clay, and granite host-rock formations and corresponding engineered design concepts for each medium. Past simulations have focused on a generic repository in bedded-salt host rock, while the most recent research has focused on a reference case for a typical clay/shale host rock. A variety of single-realization (i.e., deterministic) and multi-realization (probabilistic) results for the new clay reference case are presented, including an analysis of the effects of heat generation on repository performance, assuming a 100-year out-of-reactor commercial SNF waste form. Order-of-magnitude differences between predicted radionuclide concentrations in thermal versus isothermal simulations imply that mechanistic, coupled-process modeling in three-dimensional (3-D) domains can be important for building confidence in post-closure performance assessments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗