Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Image Processing Methods for Coronal Hole Segmentation, Matching, and Map Classification

The paper presents the results from a multi-year effort to develop and validate image processing methods forselecting the best physical models based on solar image observations. The approach consists of selecting the physical models based on their agreement with coronal holes extracted from the images. Ultimately, the goal is to use physical models to predict geomagnetic storms. We decompose the problem into three subproblems: (i) coronal hole segmentation based on physical constraints, (ii) matching clusters of coronal holes between different maps, and (iii) physical map classification. For segmenting coronal holes, we develop a multi-modal method that uses segmentation maps from three different methods to initialize a level-set method that evolves the initial coronal hole segmentation to the magnetic boundary. Then, we introducea new method based on Linear Programming for matchingclusters of coronal holes. The final matching is then performedusing Random Forests. The methods were carefully validatedusing consensus maps derived from multiple readers, manualclustering, manual map classification, and method validation for50 maps. The proposed multi-modal segmentation method significantly outperformed SegNet, U-net, Henney-Harvey, and FCNby providing accurate boundary detection. Overall, the methodgave a 95.5% map classification accuracy.

Jatla, Venkatesh↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
NASA NTRS2003

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.↗
DOE OSTI2017

A finite element/level set model of polyurethane foam expansion and polymerization

Polyurethane foams are used widely for encapsulation and structural purposes because they are inexpensive, straightforward to process, amenable to a wide range of density variations (1 lb/ft3 - 50 lb/ft3), and able to fill complex molds quickly and effectively. Computational model of the filling and curing process are needed to reduce defects such as voids, out-of-specification density, density gradients, foam decomposition from high temperatures due to exotherms, and incomplete filling. This paper details the development of a computational fluid dynamics model of a moderate density PMDI structural foam, PMDI-10. PMDI is an isocyanate-based polyurethane foam, which is chemically blown with water. The polyol reacts with isocyanate to produces the polymer. PMDI- 10 is catalyzed giving it a short pot life: it foams and polymerizes to a solid within 5 minutes during normal processing. To achieve a higher density, the foam is over-packed to twice or more of its free rise density of 10 lb/ft3. The goal for modeling is to represent the expansion, filling of molds, and the polymerization of the foam. This will be used to reduce defects, optimize the mold design, troubleshoot the processed, and predict the final foam properties. A homogenized continuum model foaming and curing was developed based on reaction kinetics, documented in a recent paper; it uses a simplified mathematical formalism that decouples these two reactions. The chemo-rheology of PMDI is measured experimentally and fit to a generalized- Newtonian viscosity model that is dependent on the extent of cure, gas fraction, and temperature. The conservation equations, including the equations of motion, an energy balance, and three rate equations are solved via a stabilized finite element method. The equations are combined with a level set method to determine the location of the foam-gas interface as it evolves to fill the mold. Understanding the thermal history and loads on the foam due to exothermicity and oven curing is very important to the results, since the kinetics, viscosity, and other material properties are all sensitive to temperature. Results from the model are compared to experimental flow visualization data and post-test X-ray computed tomography (CT) data for the density. Several geometries are investigated including two configurations of a mock structural part and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Further model improvements are also discussed for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗
DOE OSTI2017

Simulating Bubble Deformations Using the Intersection Marker (ISM) Front-Tracking method

Surface and volume tracking methods have been under continuous development over the last 50+ years since the first MAC (Marker and Cell) and SLIC (Simple Line Interface Calculation) methods appeared in the 1960's. The sustained interest in simulating multiphase flow stem from their many practical applications in the physical world; from biological systems, to chemical systems, to boiling dynamics and atmospheric heat and mass transfer. The ability to capture the details of interface evolution as well as the physics of interfacial heat and mass transfer has thus been a tempting goal for many numerical researchers. Numerical analysts, apart from the well-established and popular VOF PLIC (Volume of Fluid, Piecewise Linear Interface Construction) method, can now call upon even more computationally intensive options such as the LS (Level Set) method, the FT (Front Tracking) method, the MoF (Moment of Fluid) method and hybrid methods.

42 ENGINEERING↗
DOE OSTI2016

Deformation of contact surfaces in a vacuum interrupter after high-current interruptions

In a high-current interruption, the contact surface in a vacuum interrupter might be severely damaged by constricted vacuum arcs causing a molten area on it. As a result, a protrusion will be initiated by a transient recovery voltage after current zero, enhancing the local electric field and making breakdowns occur easier. The objective of this paper is to simulate the deformation process on the molten area under a high electric field by adopting the finite element method. A time-dependent Electrohydrodynamic model was established, and the liquid-gas interface was tracked by the level-set method. From the results, the liquid metal can be deformed to a Taylor cone if the applied electric field is above a critical value. This value is correlated to the initial geometry of the liquid metal, which increases as the size of the liquid metal decreases. Moreover, the buildup time of a Taylor cone obeys the power law t = k × E{sup −3}, where E is the initial electric field and k is a coefficient related to the material property, indicating a temporal self-similar characteristic. In addition, the influence of temperature has little impact on the deformation but has great impact on electron emission. Finally, the possible reason to initiate a delayed breakdown is associated with the deformation. The breakdown does not occur immediately when the voltage is just applied upon the gap but is postponed to several milliseconds later when the tip is formed on the liquid metal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling curved interfaces without element‐partitioning in the extended finite element method

Summary In this paper, we model holes and material interfaces (weak discontinuities) in two‐dimensional linear elastic continua using the extended finite element method on higher‐order (spectral) finite element meshes. Arbitrary parametric curves such as rational Bézier curves and cubic Hermite curves are adopted in conjunction with the level set method to represent curved interfaces. Efficient computation of weak form integrals with polynomial integrands is realized via the homogeneous numerical integration scheme—a method that uses Euler's homogeneous function theorem and Stokes' theorem to reduce integration to the boundary of the domain. Numerical integration on cut elements requires the evaluation of a one‐dimensional integral over a parametric curve, and hence, the need to partition curved elements is eliminated. To improve stiffness matrix conditioning, ghost penalty stabilization and the Jacobi preconditioner are used. For material interface problems, we develop an enrichment function that captures weak discontinuities on spectral meshes. Taken together, we show through numerical experiments that these advances deliver optimal algebraic rates of convergence with h ‐refinement ( p =1,2,…,5) and exponential rates of convergence with p ‐refinement ( p =1,2,…,7) for elastostatic problems with holes and material inclusions on Cartesian p th‐order spectral finite element meshes.

Chin, Eric B.↗
DOE OSTI2017

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A sharp interface method for compressible liquid–vapor flow with phase transition and surface tension

The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevant physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Level set topology optimization of structural problems with interface cohesion

Summary This paper presents a finite element topology optimization framework for the design of two‐phase structural systems considering contact and cohesion phenomena along the interface. The geometry of the material interface is described by an explicit level set method, and the structural response is predicted by the extended finite element method. In this work, the interface condition is described by a bilinear cohesive zone model on the basis of the traction‐separation constitutive relation. The non‐penetration condition in the presence of compressive interface forces is enforced by a stabilized Lagrange multiplier method. The mechanical model assumes a linear elastic isotropic material, infinitesimal strain theory, and a quasi‐static response. The optimization problem is solved by a nonlinear programming method, and the design sensitivities are computed by the adjoint method. The performance of the presented method is evaluated by 2D and 3D numerical examples. The results obtained from topology optimization reveal distinct design characteristics for the various interface phenomena considered. In addition, 3D examples demonstrate optimal geometries that cannot be fully captured by reduced dimensionality. The optimization framework presented is limited to two‐phase structural systems where the material interface is coincident in the undeformed configuration, and to structural responses that remain valid considering small strain kinematics. Copyright © 2017 John Wiley & Sons, Ltd.

Behrou, Reza↗
DOE OSTI2017

Numerical comparison of bubbling in a waste glass melter

Radioactive tank waste is scheduled for vitrification at the Hanford Tank Waste Treatment and Immobilization Plant, also known as the Vit Plant, being constructed at the Hanford Site. Testing of the pilot-scale, joule-heated DuraMelter 1200 at the Vitreous State Laboratory has shown that bubbling increases the melt rate, and as a result, melter throughput. Computational fluid dynamics (CFD) models of this pilot-scale waste glass melter were developed to improve our understanding of the processes that occur within the melter to aid in process optimization and troubleshooting of the Vit Plant melters. Unfortunately, model validation is complicated by the difficulty of obtaining suitable experimental data for operational melters attributed to the large-scale, radioactivity and high temperatures at which these melters are operated. This study focuses on confirming the fidelity of the CFD models to accurately reproduce the bubbling behavior. Because of the paucity of experimental data at the resolution required for CFD validation, a code-to-code comparison was used to evaluate two common approaches for simulating flows of immiscible fluids on numerical grids and resolving multiphase interfaces. Here, the volume of fluid and level set methods are used to track the dynamically evolving interfaces between the molten glass and the air bubbles. To aid in the verification of the results of these codes, a comparison of the bubble behavior, growth, and frequency of bubble generation are presented and a grid convergence study is performed for the two approaches. In addition to comparing the numerical results to available experimental data and observations from pilot-scale testing, the results are evaluated in relation to laboratory experiments with bubbles injected into viscous fluids, empirical correlations obtained from the published literature, and the basic laws of fluid dynamics. Furthermore, the results of this work can assist in validating waste glass melter models and provide a better understanding of the flow patterns within the WTP melters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗