Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Topographic controls on a phreatomagmatic maar-diatreme eruption: field and numerical results from the Holocene Dotsero volcano (Colorado, USA)

We analyzed the eruptive products of Holocene Dotsero volcano (Colorado, USA; ~4.15 ka), which recorded evidence of a progression from effusive magmatic activity through explosive phreatomagmatic maar-forming activity to a final explosive magmatic phase. The nature of the deposits suggests that the irregular and mountainous pre-eruptive topography strongly influenced the formation of directed pyroclastic density currents (PDCs) during the phreatomagmatic phase, where the topographically high northern crater rim acted as a barrier and promoted transport to the south. Furthermore, the crater shape strongly controlled the grain size of the final deposits, providing further evidence that deposit grain size can be a misleading proxy for fragmentation or explosion intensity. We test these hypotheses, as gleaned from fieldwork, by implementing a multiphase numerical model to simulate discrete explosions and their resulting PDCs at maar-diatreme volcanoes. We confirm that a crater shape where one rim is higher than the other promotes multiple pulses of particle transport in the direction of the lower rim. The initial PDCs contain both coarse and fine particles, whereas PDCs generated by later pulses contain, predominately, finer particles. This mechanism provides an alternative explanation to the formation of lapilli tuff/tuff couplets often found at maar-diatreme tephra rings (including Dotsero volcano). Additional modeling shows how the crater depth can control grain size of the resulting PDC deposits. In the case of a deep crater and a relatively weak explosion, it is possible that only the fine particles escape following collapse of the eruptive column, resulting in a fine-grained deposit. This work allows us to improve models developed to describe emplacement of PDCs for different crater geometries and informs on the dynamics of maar-diatreme eruptions, but also highlights the limitations of our knowledge of key parameters, especially those related to the processes of magma-water interaction and syn-eruptive hydrogeology, during such events.

58 GEOSCIENCES↗
DOE OSTI2019

Modeling Three-Dimensional Fluid-Driven Propagation of Multiple Fractures using TOUGH-FEMM

Here, a novel numerical simulation tool named TOUGH-FEMM is presented and applied to model three-dimensional (3D) hydraulic fracturing in porous rock. The fluid flow in both fractures and porous rock is modeled using TOUGH2, which is a well-established code for analysis of multiphase and multi-component fluid flow. Rock deformations associated with fracture propagation are modeled using finite element-meshfree method (FEMM). FEMM is an approach to simulate fracture propagation without remeshing, in which the fracture path does not need to be predetermined. Fracture mechanics with mixed-mode stress intensity factors are employed to detect fracture instability and determine the direction of fracture propagation. TOUGH-FEMM is confirmed for modeling fluid-driven fracture propagation in 3D through a number of simulation examples, including modeling of hydraulic fracturing laboratory experiments and by comparison to independent numerical simulation results for multiple interacting hydraulic fractures at ten to hundred meter scale.

58 GEOSCIENCES↗
DOE OSTI2019

Eigenparticles: characterizing particles using eigenfaces

The shape characteristics of particles have a pinnacle role in microsopic and macroscopic features of a system. Several studies have highlighted the need for considering deviations from a spherical representation of particles for accurate modeling of granular and multiphase flow systems. Using a shape factor, sphericity or roundness parameter alone is proven to be inadequate to capture the physical phenomena. In the present study we propose a novel metric based on the pattern recognition method Eigenfaces, coining the technique ‘Eigenparticles’. Using this technique we create a single statistical distribution of basis shapes to describe the morphological composition. The proposed technique is successfully validated with test shapes and applied to real particles. As a result, when compared with a state-of-the-art Fourier based method, ‘Eigenparticles’ performs favorably, clearly distinguishing the different particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Computational Biorheology of Human Blood Flow in Health and Disease

Hematologic disorders arising from infectious diseases, hereditary factors and environmental influences can lead to, and can be influenced by, significant changes in the shape, mechanical and physical properties of red blood cells (RBCs), and the biorheology of blood flow. Hence, modeling of hematologic disorders should take into account the multiphase nature of blood flow, especially in arterioles and capillaries. We present here an overview of a general computational framework based on dissipative particle dynamics (DPD) which has broad applicability in cell biophysics with implications for diagnostics, therapeutics and drug efficacy assessments for a wide variety of human diseases. This computational approach, validated by independent experimental results, is capable of modeling the biorheology of whole blood and its individual components during blood flow so as to investigate cell mechanistic processes in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to arterioles and can also be used to model RBCs down to the spectrin level. We start from experimental measurements of a single RBC to extract the relevant biophysical parameters, using single-cell measurements involving such methods as optical tweezers, atomic force microscopy and micropipette aspiration, and cell-population experiments involving microfluidic devices. We then use these validated RBC models to predict the biorheological behavior of whole blood in healthy or pathological states, and compare the simulations with experimental results involving apparent viscosity and other relevant parameters. Finally, while the approach discussed here is sufficiently general to address a broad spectrum of hematologic disorders including certain types of cancer, this paper specifically deals with results obtained using this computational framework for blood flow in malaria and sickle cell anemia.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Reactive transport benchmarks for subsurface environmental simulation

Over the last 20 years, we have seen firsthand the evolution of multicomponent reactive transport modeling and the expanding range and increasing complexity of subsurface applications it is being used to address. There is a growing reliance on reactive transport modeling (RTM) to address some of the most compelling issues facing our planet: climate change, nuclear waste management, contaminant remediation, and pollution prevention. While these issues are motivating the development of new and improved capabilities for subsurface environmental modeling using RTM (e.g., biogeochemistry from cell-scale physiology to continental-scale terrestrial ecosystems, nonisothermal multiphase conditions, coupled geomechanics), there remain longstanding challenges in characterizing the natural variability of hydrological, biological, and geochemical properties in subsurface environments and limited success in transferring models between sites and across scales. An equally important trend over the last 20 years is the evolution of modeling from a service sought out after data has been collected to a multifaceted research approach that provides (1) an organizing principle for characterization and monitoring activities; (2) a systematic framework for identifying knowledge gaps, developing and integrating new knowledge; and (3) a mechanistic understanding that represents the collective wisdom of the participating scientists and engineers. There are now large multidisciplinary projects where the research approach is model-driven, and the principal product is a holistic predictive simulation capability that can be used as a test bed for alternative conceptualizations of processes, properties, and conditions. Much of the future growth and expanded role for RTM will depend on its continued ability to exploit technological advancements in the earth and environmental sciences. Advances in measurement technology, particularly in molecular biology (genomics), isotope fractionation, and high-resolution X-ray spectroscopy, have created new lines of research that can be used to inform the conceptualization of reactions and rate laws and validate mechanistic models. For example, spectroscopy has identified the oxidation states of key components and elemental distributions at increasingly smaller scales and lower concentrations; molecular biology has progressed from identifying the presence of microbes to characterization of which microbial communities are active and what they are doing (i.e., microbial function), which has led in turn to the identification of active processes under conditions beyond what analytical chemistry can discern; isotope ratios in pore water and solid phases that can be used to distinguish between biotic from abiotic processes, sorption from precipitation, and origin and age of groundwater. The other noteworthy development that is expanding the role of RTM in subsurface environmental modeling is he advance in computational technology that is enabling the simulation of more coupled processes with increasing mechanistic detail. In some cases, this involves the inclusion of more reactive species and/or microbial populations in the simulations; in other cases, the impact is through the ability to achieve high resolution of property distributions over longer simulated times. To achieve these ambitious objectives for subsurface reactive transport simulation, the subsurface science and engineering community is being driven to provide accurate assessments of engineering performance and risk for important issues with far-reaching consequences. As a result, the complexity and detail of subsurface processes, properties, and conditions that can be simulated have significantly expanded. This expansion was enabled, in part, by advances in measurement technology, computing technology, and numerical techniques.

benchmarking↗
DOE OSTI2026

Molten Halide Salt Surface Tension: Methods and Correlations

Here, this paper reviews various methods for studying surface tension and their applicability to fluoride and chloride molten salt systems, including a comparison of benefits and drawbacks. Such a comparison aids in experiment design based on desired factors such as scale, accuracy, and repeatability. A detailed review is presented for existing literature data regarding the surface tension of molten fluoride and chloride salts. These reference data were compiled and analyzed to determine cross-validated correlation equations for several alkali and alkaline earth fluoride and chloride salts as functions of temperature. These correlations are necessary for reliable multiphysics modeling approaches as well as accurate design and analysis of multiphase molten salt phenomena such as gas sparging and bubble formation/transport. This analysis supports the development of the thermophysical arm of the Molten Salt Thermal Properties Database (MSTDB-TP) managed by Oak Ridge National Laboratory.

Chloride↗
DOE OSTI2015

Conduction below 100 °C in nominal Li 6 ZnNb 4 O 14

The increasing demand for a safe rechargeable battery with a high energy density per cell is driving a search for a novel solid electrolyte with a high Li + or Na + conductivity that is chemically stable in a working Li-ion or Na-ion battery. Li 6 ZnNb 4 O 14 has been reported to exhibit a σ Li > 10 -2 S cm -1 at 250 °C, but to disproportionate into multiple phases on cooling from 850 °C to room temperature. An investigation of the room-temperature Li-ion conductivity in a porous pellet of a multiphase product of a nominal Li 6 ZnNb 4 O 14 composition is shown to have bulk σ Li 3.3 x 10 -5 S cm -1 at room temperature that increases to 1.4 x 10 -4 S cm -1 by 50 °C. 7Li MAS NMR spectra were fitted to two Lorentzian lines, one of which showed a dramatic increase with increasing temperature. As a result, a test for water stability indicates that Li + may move to the particle and grain surfaces to react with adsorbed water as occurs in the garnet Li + conductors.

25 ENERGY STORAGE↗
DOE OSTI2018

Unordinary Ferromagnetic Phase Transitions in Mn- and Fe-doped Half-Heusler CoTiSb

Unordinary ferromagnetic character and high Curie temperature around 900 K have been detected for CoTi {sub 1−x}M{sub x}Sb (M = Mn and Fe, x = 0.1) multiphases systems. The synthesis process was carried out at high temperature in evacuated quartz tubes. The x-ray powder diffraction analysis through Rietveld method was applied to refine the crystal structure. By manipulating the search match program, three nanocrystalline phases of CoTiSb (64.5%), CoSb (15.5%), and TiSb{sub 2} (20.5%) for Mn-doped system and CoTiSb (33%), Ti{sub 3}Sb (35.3%), and CoTi{sub 2} (30.7%) for Fe-doped system were recognized. The term “phase” was used to address the coexistence of different stable chemical compositions for the same half-Heusler alloy. The scanning electron microscope (SEM) and the high-resolution transmission electron microscope (HR-TEM) were applied to characterize the morphology, size, shape, crystallinity and lattice spacing. Good matching was detected between the obtained d values using electron diffraction data and those calculated from x-ray analysis. Well defined nanocrystalline structure with average interatomic distances equaling 0.229 nm, 0.286 nm together with many successive rings with sharp diffraction spots characterized the Fe-doped system. In contrast, partially distorted structured planes and smeared halo surrounded by weak rings were detected in the Mn-doped system. Thermo-magnetic measurements (M-T) were performed for both systems. In addition to Curie transition, many magnetic phase transitions were detected and confirmed through the ΔM/ ΔT vs. T relationship. The field-dependence magnetization (M-H) was measured at different temperatures. The M-H hysteresis loops showed the ferromagnetic character for both systems. Tentative magnetic phase diagrams are suggested as a guide for future research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Low-Temperature Short-Time SPS Processes to Produce Fine-Grained High-Coercivity Barium Hexaferrite Ceramics from Polyol Nanoparticles

Barium hexaferrite (BaFe{sub 12}O{sub 19}) ceramics were obtained by combination of polyol-synthesized precursors subsequently consolidated by spark plasma sintering (SPS) at a low temperature (800 {sup ∘}C). X-ray powder diffraction studies showed the influence of the experimental parameters to produce a virtually single- or multiphase material. The room temperature magnetic properties exhibited interesting variety of properties, including exchange coupling (spring magnet) between hexaferrite and other magnetic phases. The saturation magnetization and coercive field, at room temperature, are in the 34–70 emu/g and 4.1–5.0 kOe ranges, respectively. The combination of strong coercive field and magnetization resulted in a maximum energy product [(BH)max] = 10.9 kJ/m{sup 3} for the sample with the highest BaFe{sub 12}O{sub 19} content.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Phase Diagrams and Magnetic Properties of a Spin-3/2 Blume-Capel Nanoparticle

Within the framework of Monte Carlo simulation, the critical behavior and the magnetic properties of a spherical nanoparticle with the presence of a crystal field are studied. The particle is formed by a ferromagnetic core of spin-3/2 surrounded by a ferromagnetic shell of spin-3/2 with antiferromagnetic interface coupling located on a simple cubic lattice. The ground-state phase diagrams are obtained in the plane of the crystal field versus the shell coupling and the interface coupling versus the crystal field. The effects of the shell coupling and the crystal field on both the critical and compensation temperatures are investigated. The phase diagrams exhibit a rich variety of behaviors, which include the first- and second-order phase transitions, the multiphase critical point, the isolated critical point, and the compensation point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Spectral Analysis of Linear Models of Viscoelasticity

In this paper, we examine Volterra integrodifferential equations with unbounded operator coefficients in Hilbert spaces. Equations considered are abstract hyperbolic equations perturbed by terms containing Volterra integral operators. These equations can be realized as partial integrodifferential equations that appear in the theory of viscoelasticity (see [2, 5]), as Gurtin–Pipkin integrodifferential equations (see [1, 7]) that describe finite-speed heat transfer in materials with memory. They also appear in averaging problems for multiphase media (Darcy’s law.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Revisiting Wagner’s Criteria to Predict Establishment and Retention of Alumina Scales on Ternary NiCrAl Alloys

Despite the advent of numerical modeling approaches and high-performance computing infrastructure, the design and development of corrosion-resistant high temperature alloys (> 500 °C) continue to be largely empirical and typically involve extensive experimentation. This is mainly due to the lack of a single unified physics-based model that can address the impact of multiple competing factors such as time, environment, alloy composition, microstructure, and geometry. The classical Wagner’s criteria have been foundational to estimate the minimum concentrations required of an oxide-forming element to establish and sustain a protective oxide scale. However, the formulation is primarily limited to lower-order alloy systems (binary alloys) and ignores the time dependence of subsurface compositional changes in the alloy. The lack of key data on the temperature and composition dependence of the solubility and transport of oxidants in multicomponent-multiphase alloys further exacerbates the problem. In the present work, a few of these limitations were addressed using a flux-based approach (FLAP) which tracks the spatiotemporal evolution of the fundamental flux balance between oxygen and the oxide-forming elements to enable the prediction of the formation of an external alumina scale in ternary NiCrAl alloys. In conclusion, the modeling results were validated with the literature findings and additional experimental work conducted in the present work.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Revealing pMDI Spray Initial Conditions: Flashing, Atomisation and the Effect of Ethanol

Sprays from pressurised metered-dose inhalers are produced by a transient discharge of a multiphase mixture. Small length and short time scales have made the investigation of the governing processes difficult. Consequently, a deep understanding of the physical processes that govern atomisation and drug particle formation has been elusive. X-ray phase contrast imaging and quantitative radiography were used to reveal the internal flow structure and measure the time-variant nozzle exit mass density of 50 µL metered sprays of HFA134a, with and without ethanol cosolvent. Internal flow patterns were imaged at a magnification of 194 pixels/mm and 7759 frames per second with 150 ps temporal resolution. Spray projected mass was measured with temporal resolution of 1 ms and spatial resolution 6 µm × 5 µm. The flow upstream of the nozzle comprised large volumes of vapour at all times throughout the injection. The inclusion of ethanol prevented bubble coalescence, altering the internal flow structure and discharge. Radiography measurements confirmed that the nozzle exit area is dominantly occupied by vapour, with a peak liquid volume fraction of 13%. Vapour generation in pMDIs occurs upstream of the sump, and the dominant volume component in the nozzle exit orifice is vapour at all times in the injection. Furthermore, the flow in ethanol-containing pMDIs has a bubbly structure resulting in a comparatively stable discharge, whereas the binary structure of propellant-only flows results in unsteady discharge and the production of unrespirable liquid masses.

42 ENGINEERING↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Theory and Applications of Macroscale Models in Porous Media

Systems dominated by heterogeneity over a multiplicity of scales, like porous media, still challenge our modeling efforts. The presence of disparate length- and time-scales that control dynamical processes in porous media hinders not only models predictive capabilities, but also their computational efficiency. Macrosopic models, i.e., averaged representations of pore-scale processes, are computationally efficient alternatives to microscale models in the study of transport phenomena in porous media at the system, field or device scale (i.e., at a scale much larger than a characteristic pore size). We present an overview of common upscaling methods used to formally derive macroscale equations from pore-scale (mass, momentum and energy) conservation laws. This review includes the volume averaging method, mixture theory, thermodynamically constrained averaging, homogenization, and renormalization group techniques. We apply these methods to a number of specific problems ranging from food processing to human bronchial system, and from diffusion to multiphase flow, to demonstrate the methods generality and flexibility in handling different applications. The primary intent of such an overview is not to provide a thorough review of all currently available upscaling techniques, nor a complete mathematical treatment of the ones presented, but rather a primer on some of the tools available for upscaling, the basic principles they are based upon, and their specific advantages and drawbacks, so to guide the reader in the choice of the most appropriate method for particular applications and of the most relevant technical literature.

42 ENGINEERING↗
DOE OSTI2019

Aspects of Solvers for Large-Scale Coupled Problems in Porous Media

This report summarizes solution strategies for discrete systems occurring in the simulation of processes in the subsurface. The focus is on scalable solvers for large and coupled systems. The goal of this work is to enable researchers to select suitable algorithms and parameter settings to efficiently solve their problems. The work provides an overview of existing methods, highlighting their features, potential, and also frequent pitfalls. Numerical examples are provided for single phase flow, density driven flow and poroelasticity.Aspects of multiphase flow are discussed briefly; a detailed discussion of reactive transport is beyond the scope of the article. Future trends are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗