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 235 records · Page 13

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Association of qualitative and quantitative imaging features on multiphasic multidetector CT with tumor grade in clear cell renal cell carcinoma

Purpose: The purpose of the study was to determine if enhancement features and qualitative imaging features on multiphasic multidetector computed tomography (MDCT) were associated with tumor grade in patients with clear cell renal cell carcinoma (ccRCC). Methods: In this retrospective, IRB approved, HIPAA-compliant, institutional review board-approved study with waiver of informed consent, 127 consecutive patients with 89 low grade (LG) and 43 high grade (HG) ccRCCs underwent preoperative four-phase MDCT in unenhanced (UN), corticomedullary (CM), nephrographic (NP), and excretory (EX) phases. Previously published quantitative (absolute peak lesion enhancement, absolute peak lesion enhancement relative to normal enhancing renal cortex, 3D whole lesion enhancement and the wash-in/wash-out of enhancement within the 3D whole lesion ROI) and qualitative (enhancement pattern; presence of necrosis; pattern of; tumor margin; tumor–parenchymal interface, tumor–parenchymal interaction; intratumoral vascularity; collecting system infiltration; renal vein invasion; and calcification) assessments were obtained for each lesion independently by two fellowship-trained genitourinary radiologists. Comparisons between variables included χ{sup 2}, ANOVA, and student t test. p values less than 0.05 were considered to be significant. Inter-reader agreement was obtained with the Gwet agreement coefficient (AC1) and standard error (SE) was reported. Results: No significant differences were observed between the LG and HG ccRCC cohorts with respect to absolute peak lesion enhancement and relative lesion enhancement ratio. There was a significant inverse correlation between low and high grade ccRCC and tumor enhancement the NP (71 HU vs. 54 HU, p < 0.001) and EX (52 HU vs. 39 HU, p < 0.001) phases using the 3D whole lesion ROI method. The percent wash-in of 3D enhancement from the UN to the CM phase was also significantly different between LG and HG ccRCCs (352% vs. 255%, p = 0.003). HG lesions showed significantly more calcification, necrosis, collecting system infiltration and ill-defined tumor margins (p < 0.05). Overall agreement between the two readers had a mean AC1 of 0.8172 (SE 0.0235). Conclusions: Quantitatively, high grade ccRCC had significantly lower whole lesion enhancement in the NP and EX phases on MDCT. Qualitatively, high grade ccRCC were significantly more likely to be associated with calcifications, necrosis, collecting system infiltration, and an ill-defined tumor margin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Smoothed Particle Hydrodynamics and its applications for multiphase flow and reactive transport in porous media

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and Advection-Diffusion-Reaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-scale multiphase flows and reactive transport in porous and fractured media.

Tartakovsky, Alexandre M.↗

Modification of Inclusions in Molten Steel by Mg-Ca Transfer from Top Slag: Experimental Confirmation of the ‘Refractory-Slag-Metal-Inclusion (ReSMI)’ Multiphase Reaction Model

High-temperature experiments and Refractory-Slag-Metal-Inclusion (ReSMI) multiphase reaction simulations were carried out to determine the effect of the ladle slag composition on the formation behavior of non-metallic inclusions in molten steel. Immediately after the slag-metal reaction, magnesium migrated to the molten steel and a MgAl{sub 2}O{sub 4} spinel inclusion was formed due to a reaction between magnesium and alumina inclusions. However, the spinel inclusion changed entirely into a liquid oxide inclusion via the transfer of calcium from slag to metal in the final stage of the reaction. Calcium transfer from slag to metal was more enhanced for lower SiO{sub 2} content in the slag. Consequently, the spinel inclusion was modified to form a liquid CaO-Al{sub 2}O{sub 3}-MgO-SiO{sub 2} inclusion, which is harmless under steelmaking conditions. The modification reaction was more efficient as the SiO{sub 2} content in the slag decreases.

36 MATERIALS SCIENCE↗

Numerical Simulation of Multiphase Magnetohydrodynamic Flow and Deformation of Electrolyte–Metal Interface in Aluminum Electrolysis Cells

A computational fluid dynamics based multiphase magnetohydrodynamic (MHD) flow model for simulating the melt flow and bath–metal interface deformation in realistic aluminum reduction cells is presented. The model accounts for the complex physics of the MHD problem in aluminum reduction cells by coupling two immiscible fluids, electromagnetic field, Lorentz force, flow turbulence, and complex cell geometry with large length scale. Especially, the deformation of bath–metal interface is tracked directly in the simulation, and the condition of constant anode–cathode distance (ACD) is maintained by moving anode bottom dynamically with the deforming bath–metal interface. The metal pad deformation and melt flow predicted by the current model are compared to the predictions using a simplified model where the bath–metal interface is assumed flat. The effects of the induced electric current due to fluid flow and the magnetic field due to the interior cell current on the metal pad deformation and melt flow are investigated. The presented model extends the conventional simplified box model by including detailed cell geometry such as the ledge profile and all channels (side, central, and cross-channels). The simulations show the model sensitivity to different side ledge profiles and the cross-channel width by comparing the predicted melt flow and metal pad heaving. In addition, the model dependencies upon the reduction cell operation conditions such as ACD, current distribution on cathode surface and open/closed channel top, are discussed.

42 ENGINEERING↗

Multiphase Model of Semisolid Slurry Generation and Isothermal Holding During Cooling Slope Rheoprocessing of A356 Al Alloy

In the present paper, we present an experimentally validated 3D multiphase and multiscale solidification model to understand the transport processes involved during slurry generation with a cooling slope. In this process, superheated liquid alloy is poured at the top of the cooling slope and allowed to flow along the slope under the influence of gravity. As the melt flows down the slope, it progressively loses its superheat, starts solidifying at the melt/slope interface with formation of solid crystals, and eventually exits the slope as semisolid slurry. In the present simulation, the three phases considered are the parent melt as the primary phase, and the solid grains and air as secondary phases. The air phase forms a definable air/liquid melt interface as the free surface. After exiting the slope, the slurry fills an isothermal holding bath maintained at the slope exit temperature, which promotes further globularization of microstructure. The outcomes of the present model include prediction of volume fractions of the three different phases considered, grain evolution, grain growth, size, sphericity and distribution of solid grains, temperature field, velocity field, macrosegregation and microsegregation. In addition, the model is found to be capable of making predictions of morphological evolution of primary grains at the onset of isothermal coarsening. The results obtained from the present simulations are validated by performing quantitative image analysis of micrographs of the rapidly oil-quenched semisolid slurry samples, collected from strategic locations along the slope and from the isothermal slurry holding bath.

36 MATERIALS SCIENCE↗

Numerical Prediction of Oscillation Behaviors of a Multiphase Core–Shell Droplet During Interfacial Tension Measurement

Interfacial tension between molten high-temperature materials is to be measured using the electrostatic levitation furnace (ELF) and the electromagnetic levitator (ISS-EML) aboard the International Space Station. A levitated compound droplet of a concentric core–shell structure is excited either by an impulsive electromagnetic field or by a superimposed electrostatic field. The oscillation behavior of the compound droplet is analyzed to measure interfacial tension at the interface of the two phases. In support of the space experiments, a computational fluid dynamic model was developed to characterize the oscillation behavior of multiphase core–shell droplet. The developed model predicted the interfacial tension between molten copper-rich and cobalt-rich phases with a difference of 3.3 pct compared to the values reported in literature. The developed model is being utilized to investigate the influence of various test parameters on measured surface tension and also being extended for molten steel–molten slag systems.

36 MATERIALS SCIENCE↗

Nanoscale Solute Partitioning and Carbide Precipitation in a Multiphase TRIP Steel Analyzed by Atom Probe Tomography

Nanoscale solute partitioning across multiple constituent phases in a 980-grade quenched and partitioned (Q&P) steel was analyzed using atom probe tomography (APT). The Q&P process was used to increase the C content in the retained austenite phase thereby improving its stability under plastic straining. Significant carbon enrichment of austenite was measured with decreased levels of C in martensite and almost depleted C content in ferrite, supporting the C partitioning mechanism in the literature. The APT analysis of retained austenite surrounded by martensite demonstrated a higher amount of C content compared with retained austenite surrounded by the ferrite phase. Lath and discrete carbide particle precipitation was also observed inside martensite colonies, tying up C and reducing the total amount of C available for austenite stabilization. In addition, the partitioning of Mn and other minor elements was quantitatively investigated by correlating APT and SEM-EBSD. Furthermore, these techniques provide a robust methodology for analyzing nanoscale compositional partitioning in multiphase steels, TRIP steels in particular, which can be used to better explain their microstructure-mechanical property relationships.

36 MATERIALS SCIENCE↗

CFD modeling of multiphase reacting flow in blast furnace shaft with layered burden

The ironmaking blast furnace is a counter-current chemical reactor which includes the ascending gas flow and the counter-current descending porous bed (burden). A Computational Fluid Dynamics (CFD) model has been developed to simulate the multiphase reacting flow in blast furnace shaft. The gas flow dynamics, burden movement, chemical reactions, heat and mass transfer between the gas phase and burden phase are included in the CFD model. The blast furnace burden consists of alternative layers of iron ore and coke. We propose a novel methodology in order to efficiently model the effects of alternative burden layer structure on gas flow, heat transfer, mass transfer and chemical reactions. Different reactions and heat transfer characteristics are applied for difference types of layer. In addition, the layered CFD model accurately predicts the Cohesive Zone (CZ) shape where the melting of solid burden taking place. Furthermore, we determine the shape and location of the CZ by an iterative method based on the ore temperature distribution. The theoretical formation and the methodology of the CFD model are presented and the model is applied to simulate industry blast furnaces. Finally, the proposed method can be applied to investigate the blast furnace shaft process and other moving bed system with periodic burden structure configuration.

42 ENGINEERING↗

TOUGH3: A new efficient version of the TOUGH suite of multiphase flow and transport simulators

The TOUGH suite of nonisothermal multiphase flow and transport simulators has been updated by various developers over many years to address a vast range of challenging subsurface problems. The increasing complexity of the simulated processes as well as the growing size of model domains that need to be handled call for an improvement in the simulator's computational robustness and efficiency. Moreover, modifications have been frequently introduced independently, resulting in multiple versions of TOUGH that (1) led to inconsistencies in feature implementation and usage, (2) made code maintenance and development inefficient, and (3) caused confusion to users and developers. TOUGH3—a new base version of TOUGH—addresses these issues. It consolidates both the serial (TOUGH2 V2.1) and parallel (TOUGH2-MP V2.0) implementations, enabling simulations to be performed on desktop computers and supercomputers using a single code. New PETSc parallel linear solvers are added to the existing serial solvers of TOUGH2 and the Aztec solver used in TOUGH2-MP. The PETSc solvers generally perform better than the Aztec solvers in parallel and the internal TOUGH3 linear solver in serial. TOUGH3 also incorporates many new features, addresses bugs, and improves the flexibility of data handling. So due to the improved capabilities and usability, TOUGH3 is more robust and efficient for solving tough and computationally demanding problems in diverse scientific and practical applications related to subsurface flow modeling.

97 MATHEMATICS AND COMPUTING↗

A coupled, multiphase heat flux boundary condition for the discrete element method

Similar to single-phase flows, multiphase systems with a constant heat flux at the wall are common in practice. Unlike single-phase flows, however, the numerical implementation of a constant boundary flux is non-trivial due to the coupling between phases – i.e., the partition of the total flux between the phases can vary in space and time. A numerical technique for modeling such a boundary condition is proposed here and verified via simulations of gas-solid flows. While discrete-particle simulations of monodisperse particles are considered in this report, the technique can be extended to include radiation, polydisperse systems, and/or a continuum representation of the solids phase.

42 ENGINEERING↗

Multiphase flow studies for microscale hydrodynamics in the structured packed column

Post combustion carbon capture by solvent absorption in a structured packed column is a promising technology for mitigating greenhouse gas emissions. Computational fluid dynamic (CFD) modelling of such a column is a challenging multiscale problem due to the range in length and time scales. The microscale hydrodynamics play a key role in the overall column efficiency with the interfacial area significantly influencing the mass transfer between the gas and liquid phases. In this context, multiphase flow simulations using the volume of fluid (VOF) method in a representative elementary unit (REU) of the packed column can provide fundamental insights into the microscale hydrodynamics, such as, interfacial area and liquid holdup. The present paper systematically examines the impact of various factors (e.g., physical properties and contact angle) on the interfacial area. The results are compared with existing correlations and a scaling analysis is also performed. The solvent physical properties are characterized by the Kapitza number (Ka), a dimensionless number that depends only on fluid properties. So that the Ka number decreases with increasing viscosity. At a fixed liquid load, the interfacial area and liquid holdup are observed to increase with decreasing Ka number. The impact of contact angle (i.e., solid surface characteristics) is effectively investigated by modifying the wall boundary conditions. The interfacial area and liquid holdup are found to decrease with increasing contact angle. Subsequently, a phenomenological correlation for interfacial area is proposed that includes the impact of these parameters. Finally, this correlation may be used to predict the interfacial area for gas-liquid flow in a structured packing for rivulet to fully wetted flow regimes.

42 ENGINEERING↗

High-performance organic photodetectors enabled by a refined fibrillar multiphase morphology

Film morphology of donor:acceptor blend layers plays a critical role in photon-to-current efficiency and dark/noise current of organic photodetectors (OPDs). One effective approach to manipulate crystallization and mesoscale phase separation of such blend layers is sequential casting (SC). However, the guiding strategies to control the morphology and the impacts on OPD performance of the SC films remain elusive, as the film structural evolution during SC is different from conventional blend casting (BC). In this report a refined fibrillar multiphase morphology is demonstrated by SC in an NT812:IEICO-4F blend film, where the high-quality IEICO-4F crystallites infiltrate the robust NT812 fibril network from the surface to the bulk. Such a morphology leads to improved charge generation and collection, reduced trap states, and enhanced charge block capability of resultant device, enabling simultaneous achievement of high external quantum efficiency and low dark/noise current. A maximum special detectivity of 5 × 10 13 Jones is achieved at 860 nm under –0.1 V, which is among the highest detectivities for vis-to-NIR OPDs. The linear dynamic range and response speed are also improved. Such enhancements are parallelly observed from OPD devices based on other blend systems with similar fibrillar refinement, which provides guidelines for film structure manipulation towards high OPD performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid multiphase-field model development using a modular free energy based approach with automatic differentiation in MOOSE/MARMOT

In this paper, we present a novel rapid phase-field model development capability using a modular free energy based approach in the finite element based MOOSE framework. The phase-field equations are implemented in a general form as separate code form the model specific data such as the thermodynamic free energy density and mobility functions. This isolates researchers from complicated code and allows them to focus on implementing new material models. In addition, using a native symbolic differentiation algorithm, all required derivatives of the free energy density are automatically determined, removing sources for potential human error, and guaranteeing optimally preconditioned equation systems through fully populated Jacobian matrices. Free energy terms contributing to a phase-field model are abstracted into self contained objects that can be combined at will at simulation run time. Combining multiple chemical and mechanical free energy contributions allows the construction of coupled phase-field, mechanics, and multiphase models. Through just-in-time compilation, we greatly reduce the computational time required to evaluate the differentiated expressions. Lastly, we demonstrate the new capability on a variety of example simulations.

97 MATHEMATICS AND COMPUTING↗

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Data-driven reduced order models for effective yield strength and partitioning of strain in multiphase materials

There is a critical need for the development and verification of practically useful multiscale modeling strategies for simulating the mechanical response of multiphase metallic materials with heterogeneous microstructures. In this contribution, we present data-driven reduced order models for effective yield strength and strain partitioning in such microstructures. These models are built employing the recently developed framework of Materials Knowledge Systems that employ 2-point spatial correlations (or 2-point statistics) for the quantification of the heterostructures and principal component analyses for their low-dimensional representation. The models are calibrated to a large collection of finite element (FE) results obtained for a diverse range of microstructures with various sizes, shapes, and volume fractions of the phases. The performance of the models is evaluated by comparing the predictions of yield strength and strain partitioning in two-phase materials with the corresponding predictions from a classical self-consistent model as well as results of full-field FE simulations. The reduced-order models developed in this work show an excellent combination of accuracy and computational efficiency, and therefore present an important advance towards computationally efficient microstructure-sensitive multiscale modeling frameworks. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗