Code Verification of Non-Newtonian Fluid Solvers for Single- and Two-Phase Laminar Flows

2021 ◽  
Vol 6 (2) ◽  
Author(s):  
Stefano Lovato ◽  
Serge L. Toxopeus ◽  
Just W. Settels ◽  
Geert H. Keetels ◽  
Guilherme Vaz
Keyword(s):  
Free Surface ◽  
Newtonian Fluid ◽  
Analytical Solutions ◽  
Two Phase Flow ◽  
Fluid Model ◽  
Laminar Flows ◽  
Phase Flow ◽  
Code Verification ◽  
Two Phase ◽  
Flow Solver

Abstract The presence of complex fluids in nature and industrial applications combined with the rapid growth of computer power over the past decades has led to an increasing number of numerical studies of non-Newtonian flows. In most cases, non-Newtonian models can be implemented in existing Newtonian solvers by relatively simple modifications of the viscosity. However, due to the scarcity of analytical solutions for non-Newtonian fluid flows and the widespread use of regularization methods, performing rigorous code verification is a challenging task. The method of manufactured solutions (MMS) is a powerful tool to generate analytical solutions for code verification. In this article, we present and discuss the results of three verification exercises based on MMS: (i) steady single-phase flow; (ii) unsteady two-phase flow with a smooth interface; (iii) unsteady two-phase flow with a free surface. The first and second exercises showed that rigorous verification of non-Newtonian fluid solvers is possible both on single- and two-phase flows. The third exercise revealed that “spurious velocities” typical of free-surface calculations with the Volume-of-Fluid model lead to “spurious viscosities” in the non-Newtonian fluid. The procedure is illustrated herein on a second-order finite volume flow solver, using the regularized Herschel-Bulkley fluid model as an example. The same methodology is however applicable to any flow solver and to all the rheological models falling under the class of generalized Newtonian fluid models.

Performance of a Two-Phase Flow Solver for the Simulation of Breaking Waves

2019 ◽  
Author(s):  
Qiu Jin ◽  
Dominic Hudson ◽  
Pandeli Temarel ◽  
W. Geraint Price
Keyword(s):  
Free Surface ◽  
Turbulent Flow ◽  
Cfd Simulation ◽  
Two Phase Flow ◽  
Breaking Waves ◽  
Phase Flow ◽  
Ocean Engineering ◽  
Two Phase ◽  
Flow Solver ◽  
Air Water Interface

Abstract Wave breaking is one of the most violent phenomena observed in air-water interface interactions. This phenomenon commonly occurs in real ship flows and is one of the main sources of underwater noise and white-water wakes. The investigation of this phenomenon is thus important in ship and ocean engineering. The performance of a two-phase flow solver is investigated for a simulation of spilling breaking waves generated by a shallowly submerged hydrofoil (NACA0024) in a uniform flow. An algebraic Volume of Fluid (AVOF) method is applied to capture the dynamic behaviour of the free surface and a standard k-ε turbulence model is selected to capture the turbulent flow around and downstream of the hydrofoil. The wave profiles, pressure and velocity contours are computed to investigate the overall flow conditions and a detailed analysis of the flow field downstream of the hydrofoil is conducted in terms of velocity components and turbulence intensities at six measurement sections. A comparison of the numerical and experimental results shows that an accurate representation of the free surface and the turbulent flow beneath it is obtained with the present numerical scheme. It is expected that the systematic documentation of the performance of the AVOF two-phase solver will enable its more accurate and optimal use for simulating ship-related flows, as well as increase awareness of its potential shortcomings for those interested in general CFD simulation of breaking waves.

scholarly journals Natural modes of the two-fluid model of two-phase flow

2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano
Keyword(s):  
Two Phase Flow ◽  
Fluid Model ◽  
Phase Flow ◽  
Two Phase ◽  
Natural Modes ◽  
Two Fluid Model ◽  
Two Fluid

Numerical Analysis of Two-Phase Pipe Flow of Liquid Helium Using Multi-Fluid Model

2001 ◽  
Vol 123 (4) ◽  
pp. 811-818 ◽  
Author(s):  
Jun Ishimoto ◽  
Mamoru Oike ◽  
Kenjiro Kamijo
Keyword(s):  
Phase Transition ◽  
Gas Phase ◽  
Liquid Helium ◽  
Two Phase Flow ◽  
Fluid Model ◽  
Numerical Results ◽  
Phase Flow ◽  
Two Dimensional ◽  
Two Phase ◽  
Normal Fluid

The two-dimensional characteristics of the vapor-liquid two-phase flow of liquid helium in a pipe are numerically investigated to realize the further development and high performance of new cryogenic engineering applications. First, the governing equations of the two-phase flow of liquid helium based on the unsteady thermal nonequilibrium multi-fluid model are presented and several flow characteristics are numerically calculated, taking into account the effect of superfluidity. Based on the numerical results, the two-dimensional structure of the two-phase flow of liquid helium is shown in detail, and it is also found that the phase transition of the normal fluid to the superfluid and the generation of superfluid counterflow against normal fluid flow are conspicuous in the large gas phase volume fraction region where the liquid to gas phase change actively occurs. Furthermore, it is clarified that the mechanism of the He I to He II phase transition caused by the temperature decrease is due to the deprivation of latent heat for vaporization from the liquid phase. According to these theoretical results, the fundamental characteristics of the cryogenic two-phase flow are predicted. The numerical results obtained should contribute to the realization of advanced cryogenic industrial applications.

A Physically Based, One-Dimensional Two-Fluid Model for Direct Contact Condensation of Steam Jets Submerged in Subcooled Water

2015 ◽  
Vol 1 (2) ◽  
Author(s):  
David Heinze ◽  
Thomas Schulenberg ◽  
Lars Behnke
Keyword(s):  
Direct Contact ◽  
Two Phase Flow ◽  
Fluid Model ◽  
Interfacial Area ◽  
Phase Flow ◽  
Two Phase ◽  
One Dimensional ◽  
Two Fluid Model ◽  
Contact Condensation ◽  
Two Fluid

A simulation model for the direct contact condensation of steam in subcooled water is presented that allows determination of major parameters of the process, such as the jet penetration length. Entrainment of water by the steam jet is modeled based on the Kelvin–Helmholtz and Rayleigh–Taylor instability theories. Primary atomization due to acceleration of interfacial waves and secondary atomization due to aerodynamic forces account for the initial size of entrained droplets. The resulting steam-water two-phase flow is simulated based on a one-dimensional two-fluid model. An interfacial area transport equation is used to track changes of the interfacial area density due to droplet entrainment and steam condensation. Interfacial heat and mass transfer rates during condensation are calculated using the two-resistance model. The resulting two-phase flow equations constitute a system of ordinary differential equations, which is solved by means of the explicit Runge–Kutta–Fehlberg algorithm. The simulation results are in good qualitative agreement with published experimental data over a wide range of pool temperatures and mass flow rates.

Eulerian Two-Phase Flow Modeling of Steam Direct Contact Condensation for the Fukushima Accident Investigation

2014 ◽  
Author(s):  
Marco Pellegrini ◽  
Giulia Agostinelli ◽  
Hidetoshi Okada ◽  
Masanori Naitoh
Keyword(s):  
Two Phase Flow ◽  
Volume Fraction ◽  
Fluid Model ◽  
Fluid Dynamic ◽  
Interfacial Region ◽  
Phase Flow ◽  
Two Phase ◽  
Fukushima Accident ◽  
Average Cell Size ◽  
Average Cell

Steam condensation is characterized by a relatively large interfacial region between gas and liquid which, in computational fluid dynamic (CFD) analyses, allows the creation of a discretized domain whose average cell size is larger than the interface itself. For this reason generally one fluid model with interface tracking (e.g. volume of fluid method, VOF) is employed for its solution in CFD, since the solution of the interface requires a reasonable amount of cells, reducing the modeling efforts. However, for some particular condensation applications, requiring the computation of long transients or the steam ejected through a large number of holes, one-fluid model becomes computationally too expensive for providing engineering information, and a two-fluid model (i.e. Eulerian two-phase flow) is preferable. Eulerian two-phase flow requires the introduction of closure terms representing the interactions between the two fluids in particular, in the condensation case, drag and heat transfer. Both terms involve the description of the interaction area whose definition is different from the typical one adopted in the boiling analyses. In the present work a simple but effective formulation for the interaction area is given based on the volume fraction gradient and then applied to a validation test case of steam bubbling in various subcooling conditions. It has been shown that this method gives realistic values of bubble detachment time, bubble penetration for the cases of interest in the nuclear application and in the particular application to the Fukushima Daiichi accident.

scholarly journals TWO-PHASE FLOW SIMULATIONS OF SCOUR AROUND VERTICAL AND HORIZONTAL CYLINDERS

2018 ◽  
pp. 22
Author(s):  
Tim Nagel ◽  
Julien Chauchat ◽  
Cyrille Bonamy ◽  
Antoine Mathieu ◽  
Xiaofeng Liu ◽  
...  
Keyword(s):  
Sediment Transport ◽  
Two Phase Flow ◽  
Numerical Data ◽  
Horizontal Cylinder ◽  
Transport Processes ◽  
Phase Flow ◽  
Two Phase ◽  
Flow Solver ◽  
Sediment Bed ◽  
Horizontal Cylinders

Scour around structures is a major engineering issue that requires a detailed description of the flow field as well as sediment transport processes. Due to enhanced suspended load associated with vortices generated around structures, sediment transport cannot be solely related to bed shear stress, such as Shields parameter based formula. In order to address this issue, we used a multi-dimensional two-phase flow solver, sedFoam-2.0 (Chauchat et al., GMD 2017) implemented under the open-source CFD toolbox OpenFOAM. Three configurations are studied and compared with experimental and numerical data from the literature. First, the 2D configurations of an horizontal cylinder lying on a sediment bed (Mao, 1986; Sumer et al., 2001) are investigated. Then, the 3D configuration of the scour around a vertical cylindrical pile reported by Roulund et al. (2005) for rigid-bed and live bed cases is investigated.

Analysis of Flow-Induced Vibration Due to Stratified Wavy Two-Phase Flow

2016 ◽  
Vol 138 (9) ◽  
Author(s):  
Shuichiro Miwa ◽  
Takashi Hibiki ◽  
Michitsugu Mori
Keyword(s):  
Two Phase Flow ◽  
Dynamic Properties ◽  
Fluid Model ◽  
Dominant Frequency ◽  
Phase Flow ◽  
Flow Induced Vibration ◽  
Two Phase ◽  
Pipe Bend ◽  
Experimental Database ◽  
Force Fluctuation

Fluctuating force induced by horizontal gas–liquid two-phase flow on 90 deg pipe bend at atmospheric pressure condition is considered. Analysis was conducted to develop a model which is capable of predicting the peak force fluctuation frequency and magnitudes, particularly at the stratified wavy two-phase flow regime. The proposed model was developed from the local instantaneous two-fluid model, and adopting guided acoustic theory and dynamic properties of one-dimensional (1D) waves to consider the collisional force due to the interaction between dynamic waves and structure. Comparing the developed model with experimental database, it was found that the main contribution of the force fluctuation due to stratified wavy flow is from the momentum and pressure fluctuations, and collisional effects. The collisional effect is due to the fluid–solid interaction of dynamic wave, which is named as the wave collision force. Newly developed model is capable of predicting the force fluctuations and dominant frequency range with satisfactory accuracy for the flow induced vibration (FIV) caused by stratified wavy two-phase flow in 52.5 mm inner diameter (ID) pipe bend.

Gradient-Augmented Level Set Two-Phase Flow Method With Pretreated Reinitialization for Three-Dimensional Violent Sloshing

2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Jianjian Xin ◽  
Fulong Shi ◽  
Qiu Jin ◽  
Lin Ma
Keyword(s):  
Free Surface ◽  
Level Set ◽  
Two Phase Flow ◽  
Three Dimensional ◽  
Surface Shape ◽  
Phase Flow ◽  
Two Phase ◽  
Correction Technique ◽  
Highly Nonlinear ◽  
Violent Sloshing

Abstract A three-dimensional (3D) gradient-augmented level set (GALS) two-phase flow model with a pretreated reinitialization procedure is developed to simulate violent sloshing in a cuboid tank. Based on a two-dimensional (2D) GALS method, 3D Hermite, and 3D Lagrange polynomial schemes are derived to interpolate the level set function and the velocity field at arbitrary positions over a cell, respectively. A reinitialization procedure is performed on a 3D narrow band to treat the strongly distorted interface and improve computational efficiency. In addition, an identification-correction technique is proposed and incorporated into the reinitialization procedure to treat the tiny droplet which can distort the free surface shape, even lead to computation failure. To validate the accuracy of the present GALS method and the effectiveness of the proposed identification-correction technique, a 3D velocity advection case is first simulated. The present method is validated to have better mass conservation property than the classical level set and original GALS methods. Also, distorted and thin interfaces are well captured on all grid resolutions by the present GALS method. Then, sloshing under coupled surge and sway excitation, sloshing under rotational excitation are simulated. Good agreements are obtained when the present wave and pressure results are compared with the experimental and numerical results. In addition, the highly nonlinear free surface is observed, and the relationship between the excitation frequency and the impulsive pressure is investigated.

Sign in / Sign up

Export Citation Format

Share Document