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 487 records · Page 27

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

A multiscale linear-solver framework for the pressure equation associated with flow in highly heterogeneous porous formations was developed. The multiscale based approach is cast in a general algebraic form, which facilitates integration of the new scalable linear solver in existing flow simulators. The Algebraic Multiscale Solver (AMS) is employed as a preconditioner within a multi-stage strategy. The formulations investigated include the standard MultiScale Finite-Element (MSFE) andMultiScale Finite-Volume (MSFV) methods. The local-stage solvers include incomplete factorization and the so-called Correction Functions (CF) associated with the MSFV approach. Extensive testing of AMS, as an iterative linear solver, indicate excellent convergence rates and computational scalability. AMS compares favorably with advanced Algebraic MultiGrid (AMG) solvers for highly detailed three-dimensional heterogeneous models. Moreover, AMS is expected to be especially beneficial in solving time-dependent problems of coupled multiphase flow and transport in large-scale subsurface formations.

02 PETROLEUM↗

Evaluation of liquid-phase thermometry in impinging jet sprays using synchrotron x-ray scattering

Liquid thermometry during primary and secondary breakup of liquid sprays is challenging due to the presence of highly dynamic, optically complex flow features. This work evaluates the use of x-ray scattering from a focused, monochromatic beam of the Advanced Photon Source at Argonne National Laboratory for the measurement of liquid temperatures within the mixing zone of an impinging jet spray. The measured scattering profiles are converted to temperature through a previously developed two-component partial least squares (PLS) regression model. Transmitive mixing during jet merging is inferred through spatial mapping of temperatures within the impingement region. In this work, the technique exhibits uncertainties of ±2K in temperature and 2% in capturing the correct scattering profile, showing its potential utility for probing liquid temperature distributions in multiphase flows.

42 ENGINEERING↗

Cryogenic spray vaporization in high-velocity helium, argon and nitrogen gasflows

Effects of gas properties on cryogenic liquid-jet atomization in high-velocity helium, nitrogen, and argon gas flows were investigated. Volume median diameter, D(sub v.5e), data were obtained with a scattered-light scanning instrument. By calculating the change in spray drop size, -Delta D(sub v.5)(exp 2), due to droplet vaporization, it was possible to calculate D(sub v.5C). D(sub v.5C) is the unvaporized characteristic drop size formed at the fuel-nozzle orifice. This drop size was normalized with respect to liquid-jet diameter, D(sub O). It was then correlated with several dimensionless groups to give an expression for the volume median diameter of cryogenic LN2 sprays. This expression correlates drop size D(sub v.5c) with aerodynamic and liquid-surface forces so that it can be readily determined in the design of multiphase-flow propellant injectors for rocket combustors.

Ingebo, Robert D.↗

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗

High-resolution time-resolved experiments on mixing and entrainment of buoyant jets in stratified environments

Fluid jets interacting with a stratified layer play an important role in the safety of several reactor designs. In the containment of nuclear power plants, fluid jets dominate the transport and mixing of gaseous species and consequent hydrogen distribution in case of a severe accident. The mixing phenomena in the containment are driven by buoyant high-momentum injections (jets) and low momentum injection plumes. Mixing near the postulated break is initially dominated by high flow velocities. Plumes with moderate flow velocities are instead relevant in the break compartment during the long-term pressurization phase, or in any of the apertures between two connected compartments if the mass flows are sufficiently high and the density differences between efflux and ambient are sufficiently low. Phenomena of interest include free plumes (as produced by the efflux from the break compartment in a larger room or directly from a break flow), wall plumes (such those produced by low mass flows through inter-compartment apertures), and propagating stratification fronts in the ambient (for any stably stratified conditions). These phenomena have been highly ranked about nuclear reactor design, especially regarding of safety protocols. During a Pressurized Thermal Shock (PTS) scenario, the interaction between the cold ECCS injection plume and the stratified fluid present in the cold (or hot) leg is important in order to determine the temperature at the time-dependent temperature at the inlet of the reactor pressure vessel (RPV) and the potential to cause a thermal shock on the RPV wall. In sodium-cooled fast reactors (SFRs), core channels are typically hydro-dynamically isolated so that there exists a considerable temperature variation at the exit of adjacent fuel assemblies. All the above phenomena are characterized by the interaction of buoyant jets with the stratified flow. In stratified layers baroclinic forces create significant redistribution of turbulent kinetic energy and scales, which leads to anisotropy. This important physical phenomenon is highly three dimensional and is challenging to capture even with high-fidelity CFD simulations, due in part to lack of sufficiently resolved validation data. Furthermore, the experimental data available in the open literature do not feature the level of fidelity needed for an extensive validation of turbulence models in lower order CFD. To shed new lights into the crucial phenomena object of the present research project, it was proposed to conduct coordinated experiments and simulations at the University of Michigan and the George Washington University. The project has resulted in an experimental database of high-resolution time-resolved measurements of jets in uniform and stratified environments. The novel experimental data will be used to validate computational fluid dynamic (CFD) codes, including both Large Eddy Simulations (LES) and unsteady Reynolds-averaged Navier-Stokes equations (URANS) methodologies. In the Experimental and Multiphase flow (ECMF) laboratory at Univerisity of Michigan, we built two experimental facilities to investigate also the effect of scaling. The first facility, DESTROJER (DEnsity Stratified Turbulent ROund free Jet ExpeRiment), featuring a contoured jet nozzle with a diameter of D=12.7mm and a 1m×1m×1m cubic tank, which is made of acrylic glass for optical access. The ratio between the tank width and the nozzle diameter is equal to 78, which ensures that there is no direct interaction between the jet and the side walls. A second, modular experimental facility, features three different tank sizes of size 10×10×30, 20×20×30, 30×30×30 cm3 respectively (all tanks have the same height of 30 cm), and a jet diameter of 2mm. For the modular facility, tank-to-nozzle ratio of 50, 100, and 150 are obtained respectively. Experiments with different-density jet impactions and sharp interface with a density difference of 3.16% In the Laboratory at George Washington University, complementary experiments have been carried out at a facility featuring the same dimensions as DESTROYER. Experiments with different jet diameters have been carried out as well. High-resolution time-resolved measurements have been performed at all facilities using high-speed particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF). In this way, simultaneous measurements of velocity and density fields have been carried out. At the University of Michigan, wire-mesh sensors have been designed and built to gather high-resolution, time-dependent data of the stratified front in the tank. These measurements are important to provide accurate initial conditions for the CFD simulations. An important outcome of the project has been the development of novel techniques to achieve refractive index matching (RIM) of acqueous solutions with high density differences. RIM is crucial if optical techniques such as PIV are used to measure velocity fields. Previous to this project, experiments up to 3% density difference had been reported in the literature. With the methodology developed in this project, we have been able to push the envelop to density differences up to 9%. The experiments have been used to validate Reynolds Averaged Navier-Stokes (RANS) turbulence models and Large Eddy Simulations (LES) models. Both NEK5000 and the commercial CFD code STAR-CCM+ have been used. The experimental data will be made available to the community through the NE-KAMS databasehave been successfully conducted in the scaled facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN

To better inform the subsurface scientist on the expected performance of parallel simulators, this work investigates performance of the reactive multiphase flow and multicomponent biogeochemical transport code PFLOTRAN as it is applied to several realistic modeling scenarios run on the Jaguar supercomputer. After a brief introduction to the code's parallel layout and code design, PFLOTRAN's parallel performance (measured through strong and weak scalability analyses) is evaluated in the context of conceptual model layout, software and algorithmic design, and known hardware limitations. PFLOTRAN scales well (with regard to strong scaling) for three realistic problem scenarios: (1) in situ leaching of copper from a mineral ore deposit within a 5-spot flow regime, (2) transient flow and solute transport within a regional doublet, and (3) a real-world problem involving uranium surface complexation within a heterogeneous and extremely dynamic variably saturated flow field. Weak scalability is discussed in detail for the regional doublet problem, and several difficulties with its interpretation are noted.

54 ENVIRONMENTAL SCIENCES↗

Transient seafloor venting on continental slopes from warming‐induced methane hydrate dissociation

Methane held in frozen hydrate cages within marine sediment comprises one of the largest carbon reservoirs on the planet. Recent submarine observations of widespread methane seepage may record hydrate dissociation due to oceanic warming, which consequently may further amplify climate change. Here, in this work, we simulate the effect of seafloor warming on marine hydrate deposits using a multiphase flow model. We show that hydrate dissociation, gas migration, and subsequent hydrate formation can generate temporary methane venting into the ocean through the hydrate stability zone. Methane seeps venting through the hydrate stability zone on the eastern Atlantic margin may record this process due to warming begun thousands of years ago. Our results contrast with the traditional view that venting occurs only updip of the hydrate stability zone.

hydrate↗

Experimental analysis of spatial correlation effects on capillary trapping of supercritical CO 2 at the intermediate laboratory scale in heterogeneous porous media

Several numerical studies have demonstrated that the heterogeneous nature of typical sedimentary formations can favorably dampen the accumulation of mobile CO 2 phase underneath the caprock. Core flooding experiments have also shown that contrasts in capillary entry pressure can lead to buildup of nonwetting fluid phase (NWP) at interfaces between facies. Explicit representation of geological heterogeneity at the intermediate (cm-to-m) scale is a powerful approach to identify the key mechanisms that control multiphase flow dynamics in porous media. The ability to carefully control flow regime and permeability contrast at a scale that is relevant to CO 2 plume dynamics in saline formations offers valuable information to understand immiscible displacement processes and provides a benchmark for mathematical models. To provide insight into the impact of capillary heterogeneity on flow dynamics and trapping efficiency of supercritical CO 2 under successive drainage and imbibition conditions, we present an experimental investigation conducted in a synthetic sand reservoir. By mimicking the interplay of governing forces at reservoir conditions via application of surrogate fluids, we performed three immiscible displacement experiments to observe the entrapment of NWP in heterogeneous porous media. Capillary trapping performance is evaluated for each scenario through spatial and temporal variations of NWP saturation; for this reason we adopted X-ray attenuation to precisely measure phase saturation throughout the flow domain and apply spatial moment analysis. Furthermore, the sweeping performance of two different permeability fields with comparable variance but distinct spatial correlation was compared against a homogeneous base case with equivalent mean permeability by means of spatial moment analysis.

58 GEOSCIENCES↗

An extended JBN method of determining unsteady-state two-phase relative permeability

Relative permeability is the reduction of permeability of porous media when subjected to multiphase flow and a key parameter in subsurface hydrology. The JBN method is a well-known method of obtaining relative permeability, which measures the overall pressure drop and the effluent phase ratio versus time during two-phase displacements. By assuming no capillary pressure or gravity, the JBN method obtains the relative permeabilities to both phases at the core outlet. Since data across a range of saturations are acquired in a relatively short time, this method is widely used. This work extends the JBN method by having (1) section-wise pressure drop measurements between the core inlet, four pressure taps on the core and the outlet, (2) local saturation measurements, and (3) local phase fluxes. With these data, the extended JBN method can determine relative permeabilities to both phases at each pressure tap of the core (not just at the core outlet). The JBN extension is shown using a data set where CO 2 invades a brine-filled core. From this it is found that the advantages of the extended JBN method over the regular JBN method are: (1) four times more data are obtained, and (2) data are more accurate because the capillary end effect is experimentally avoided. Avoiding the end effect results in tripling the saturation range, and obtaining relative permeabilities that are consistent with steady-state measurements and roughly 40% higher than those from the regular JBN method.

54 ENVIRONMENTAL SCIENCES↗

Length scale selects directionality of droplets on vibrating pillar ratchet

Directional control of droplet motion at room temperature is of interest for applications such as microfluidic devices, self-cleaning coatings, and directional adhesives. Here, arrays of tilted pillars ranging in height from the nanoscale to the microscale are used as structural ratchets to directionally transport water at room temperature. Water droplets deposited on vibrating chips with a nanostructured ratchet move preferentially in the direction of the feature tilt while the opposite directionality is observed in the case of microstructured ratchets. This remarkable switch in directionality is consistent with changes in the contact angle hysteresis. To glean further insights into the length scale dependent asymmetric contact angle hysteresis, the contact lines formed by a nonvolatile room temperature ionic liquid placed onto the tilted pillar arrays were visualized and analyzed in situ in a scanning electron microscope. As a result, the ability to tune droplet directionality by merely changing the length scale of surface features all etched at the same tilt angle would be a versatile tool for manipulating multiphase flows and for selecting droplet directionality in other lap-on-chip applications.

42 ENGINEERING↗

Experimental investigation on multiscale hydrodynamics in a novel gas–Liquid–Solid three phase jet-Loop reactor

Multiphase flow hydrodynamics in a novel gas–liquid–solid jet-loop reactor (JLR) were experimentally investigated at the macroscales and mesoscales. The chord length distribution was measured by an optical fiber probe and transformed for bubble size distribution through the maximum entropy method. The impacts of key operating conditions (superficial gas and liquid velocity, solid loading) on hydrodynamics at different axial and radial locations were comprehensively investigated. JLR was found to have good solid suspension ability owing to the internal circulation of bubbles and liquid flow. The gas holdup, axial liquid velocity, and bubble velocity increase with gas velocity, while liquid velocity has little influence on them. Compared with the gas–liquid JLRs, solids decrease the gas holdup and liquid circulation, reduces the bubble velocity and delays the flow development due to the enhanced interaction between bubbles and particles (Stokes number >1). This work also provides a benchmark data for computational fluid dynamics (CFD) model validation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of drilling fluid rheology and wellbore pressure on rock cuttings removal performance: numerical investigation

Drilling fluids exhibit non-Newtonian behavior which is generally described in the literature by two rheological models: Power Law and Herschel-Bulkley models. We have used Eulerian- Eulerian multiphase flow model to evaluate the drilling fluid flow patterns and filter cake formation patterns for both models for ultra-deep drilling conditions. Simulated drilling fluid flow patterns clearly show that the Herschel - Bulkley model exhibits fewer vortices around the drilling pipe, indicating better rock cuttings removal performance. In addition to mud rheology, we have analyzed the pressure drop variation with the filter cake formation on the wellbore during ultradeep drilling process. We have shown that filter cake thickness at the wall linearly increases with differential pressure between the well and formation. This implies that the differential pressure needs to be optimized. Very thin filter cake may lead to excessive fluid seepage into the permeable subterranean vertical zone resulting in a loss of fluid circulation and hence in a loss of rock cuttings removal. Yet, very thick filter cake reduces the effective diameter of the hole which may lead to a poor rock cuttings removal performance and stuck pipe. The mathematical model presented here needs to be experimentally validated when filter cake thickness data in deep drilling conditions become available. The model could be useful to provide drilling engineers insight on the filter cake formation, drilling fluid flow pattern, and give guide-line to optimize the filter cake as field experimental data are scarce.

02 PETROLEUM↗

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↗

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↗

Characterization of immiscible fluid displacement processes with various capillary numbers and viscosity ratios in 3D natural sandstone

We characterize the influence of reservoir conditions upon multiphase flow, we calculated fluid displacements (drainage processes) in 3D pore spaces of Berea sandstone using two-phase lattice Boltzmann (LB) simulations. The results of simulations under various conditions were used to classify the resulting two-phase flow behavior into three typical fluid displacement patterns on the diagram of capillary number (Ca) and viscosity ratio of the two fluids (M). In addition, the saturation of the nonwetting phase was calculated and mapped on the Ca–M diagram. We then characterized dynamic pore-filling events (i.e., Haines jumps) from the pressure variation of the nonwetting phase, and linked this behavior to the occurrence of capillary fingering. The results revealed the onset of capillary fingering in 3D natural rock at a higher Ca than in 2D homogeneous granular models, with the crossover region between typical displacement patterns broader than in the homogeneous granular model. Furthermore, saturation of the nonwetting phase mapped on the Ca–M diagram significantly depends on the rock models. These important differences between two-phase flow in 3D natural rock and in 2D homogeneous models could be due to the heterogeneity of pore geometry in the natural rock and differences in pore connectivity. By quantifying two-phase fluid behavior in the target reservoir rock under various conditions (e.g., saturation mapping on the Ca–M diagram), our approach could provide useful information for investigating suitable reservoir conditions for geo-fluid management (e.g., high CO 2 saturation in CO 2 storage).

58 GEOSCIENCES↗