MFiX-Multiphase Flow with Interphase Exchanges
Presented by LRST at the NETL 2020 FE R&D Virtual Project Review Meeting, 8/24
SEARCH · Engineering Papers
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.
Presented by LRST at the NETL 2020 FE R&D Virtual Project Review Meeting, 8/24
Apparatus and methods for the measurement of gas volume fraction of produced oil are described. A first method measures the response of a pipe containing the produced oil excited by a source of vibration in the form of an acoustic frequency chirp containing a linearly varying range of frequencies in the tens of kilohertz range encompassing at least one resonant mode of the pipe. As the gas volume fraction increases, the location of the peak maximum of the measured frequency spectrum responsive to the excitation increases in frequency, and the height of the peak maximum increases, thereby permitting a linear calibration curve to be obtained. A second method measures the response of a pipe containing the produced oil to excitation by a continuous source of vibration having a chosen frequency above those which excite flexural vibrations in the pipe and simultaneously excite acoustic waves in the fluid contained in the pipe, known as the coincidence frequency. Gas present in the fluid will interrupt sound propagation or reverberation, thereby generating fluctuations in the amplitude of the measured vibrations of the pipe. The amplitude fluctuation level provides a measure of the gas volume present inside the pipe. A third method measures the response of a pipe containing the produced oil to excitation by a high-bandwidth, short pulse having a chosen center frequency above the coincidence frequency. Gas present in the fluid will interrupt pulse propagation, thereby generating fluctuations in the amplitude of the measured vibrations of the pipe.
Abstract not provided.
Abstract not provided.
2021 Spring FE R&D Virtual Project Review Meeting, Virtual, May 24, 2021
FY22 FECM Spring R&D Project Review Meeting, Virtual, May 6, 2022
Abstract not provided.
For presentation at the Clearwater Clean Energy Conference, Clearwater, FL, July 23-28, 2023.
General overview of MFIX-Exa and weak scaling data from Frontier.
Explore the source record for details and available documents.
This paper provides an overview of the microgravity fluid physics and transport phenomena experiments planned for the International Spare Station. NASA's Office of Life and Microgravity Science and Applications has established a world-class research program in fluid physics and transport phenomena. This program combines the vast expertise of the world research community with NASA's unique microgravity facilities with the objectives of gaining new insight into fluid phenomena by removing the confounding effect of gravity. Due to its criticality to many terrestrial and space-based processes and phenomena, fluid physics and transport phenomena play a central role in the NASA's Microgravity Program. Through widely publicized research announcement and well established peer-reviews, the program has been able to attract a number of world-class researchers and acquired a critical mass of investigations that is now adding rapidly to this field. Currently there arc a total of 106 ground-based and 20 candidate flight principal investigators conducting research in four major thrust areas in the program: complex flows, multiphase flow and phase change, interfacial phenomena, and dynamics and instabilities. The International Space Station (ISS) to be launched in 1998, provides the microgravity research community with a unprecedented opportunity to conduct long-duration microgravity experiments which can be controlled and operated from the Principal Investigators' own laboratory. Frequent planned shuttle flights to the Station will provide opportunities to conduct many more experiments than were previously possible. NASA Lewis Research Center is in the process of designing a Fluids and Combustion Facility (FCF) to be located in the Laboratory Module of the ISS that will not only accommodate multiple users but, allow a broad range of fluid physics and transport phenomena experiments to be conducted in a cost effective manner.
Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.
A numerical model for electrokinetic flow of multiphase systems with deformable interfaces is presented, based on a combined level set-volume of fluid technique. A new feature is a multiphase formulation of the Nernst–Planck transport equation for advection, diffusion and conduction of individual charge carrier species that ensures their conservation in each fluid phase. The numerical model is validated against the analytical results of Zholkovskij et al. (2002) [1], and results for the problem of two drops coalescing in the presence of mobile charge carriers are presented. The time taken for two drops containing ions to coalesce decreases with increasing ion concentration.
In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.
Simulation of multiphase flow systems are of critical importance in managing hydrological systems. Flow simulations are affected by a number of factors including structure and flow properties including porosity and permeability as well as the anisotropy and heterogeneity of these properties. In many cases traditional hydrological and reservoir data are highly affected by these parameters, but are not directly sensitive to them. As such modellers often adjust these parameters in an ad hoc manner until solutions numerically converge. Simulation models are generally based on structural data from reflection seismics whose physical flow properties are then populated using geostatistical extrapolation techniques utilizing a sparse number of borehole logs and core analysis. In multiphase systems including enhanced oil recovery and carbon capture and sequestration uncertainties regarding phase-dependent physical properties confounds this challenge further. Geophysical methods provide a means by which to gain an improved understanding of phase distributions in the subsurface. In this paper we will look at applications from active carbon capture and sequestration and enhanced oil recovery applications, as well as synthetic examples. Geophysical data including electromagnetic and gravity are inverted using structural constraints from the reservoir model. Inversions are then mapped into flow properties using calibrated relations such as Archie’s Equation. The coupled models can then be used to both verify and improve on the reservoir flow model which improves it’s predictive power and utility as a management tool.
A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).
The simulation of multiphase flow is a difficult physical, computational, and numerical problem. Sophisticated tools are necessary to provide results encompassing a wide range of physical phenomena. The development of these tools must be supported by rigorous analysis and ongoing efforts of verification and validation to ensure the accuracy of simulations. The Nek-2P multiphase code is being developed at Argonne with the general goal of simulating a wide range of multiphase flows and the primary goal of simulating boiling flows in LWR fuel assemblies. As part of this, there is an ongoing effort to verify the code and validate the underlying physical models. The objective of this work is to present recent results in the verification of Nek-2P. As a starting point, this work is focused on adiabatic flows for the sake of simplicity with the immediate next planned step to expand the effort to include heat and mass transfer. Details of the verification of the implemented mass fraction transport, momentum conservation, and pressure equations in Nek-2P are presented. Results of simulations are compared to analytic solutions when possible, however due to the complexity of multiphase flows, these are seldom available. To overcome this difficulty, an independent model is used to provide solutions for comparison.