Validation of CFD and Simplified Models With Experimental Data for Multiphase Flow in Bends

2013 ◽  
Author(s):  
E. Nennie ◽  
S. P. C. Belfroid ◽  
T. S. D. O’Mahoney

In this paper details of the measurement results of the forces on the bends in a 4″ setup are compared to two models. The first model is a simple analytical model and is used to estimate the forces. In the second model, CFD is used. In the experiments only resulting forces, including upstream and downstream bends and mechanical resonance and interaction is measured. The goal of the CFD was to discriminate between flow and mechanics and to evaluate the influence of a flow disturbance as a result of a bend on the force on a downstream second bend. For the simplified analytical model the amplitudes are underestimated, but the frequency spectra look very reasonable in case of the slug flow regime. The main advantage of the simplified analytical model is that the computational time is in the order of seconds, but the accuracy is still reasonable for the use in an engineering approach of determining the structural integrity of the complete pipe system. For the CFD the shape of the force function is similar to the experiments. The CFD results indicate that the forces on the second downstream bend are also measured on the first upstream bend. The accuracy of the CFD simulation is the advantage of this model, but the computational time is very long, especially if the multiphase flow simulation is coupled to the structural model.

2021 ◽  
Vol 104 (1) ◽  
pp. 003685042110080
Author(s):  
Zheqin Yu ◽  
Jianping Tan ◽  
Shuai Wang

Shear stress is often present in the blood flow within blood-contacting devices, which is the leading cause of hemolysis. However, the simulation method for blood flow with shear stress is still not perfect, especially the multiphase flow model and experimental verification. In this regard, this study proposes an enhanced discrete phase model for multiphase flow simulation of blood flow with shear stress. This simulation is based on the discrete phase model (DPM). According to the multiphase flow characteristics of blood, a virtual mass force model and a pressure gradient influence model are added to the calculation of cell particle motion. In the experimental verification, nozzle models were designed to simulate the flow with shear stress, varying the degree of shear stress through different nozzle sizes. The microscopic flow was measured by the Particle Image Velocimetry (PIV) experimental method. The comparison of the turbulence models and the verification of the simulation accuracy were carried out based on the experimental results. The result demonstrates that the simulation effect of the SST k- ω model is better than other standard turbulence models. Accuracy analysis proves that the simulation results are accurate and can capture the movement of cell-level particles in the flow with shear stress. The results of the research are conducive to obtaining accurate and comprehensive analysis results in the equipment development phase.


Author(s):  
Seyed Kourosh Mahjour ◽  
Antonio Alberto Souza Santos ◽  
Manuel Gomes Correia ◽  
Denis José Schiozer

AbstractThe simulation process under uncertainty needs numerous reservoir models that can be very time-consuming. Hence, selecting representative models (RMs) that show the uncertainty space of the full ensemble is required. In this work, we compare two scenario reduction techniques: (1) Distance-based Clustering with Simple Matching Coefficient (DCSMC) applied before the simulation process using reservoir static data, and (2) metaheuristic algorithm (RMFinder technique) applied after the simulation process using reservoir dynamic data. We use these two methods as samples to investigate the effect of static and dynamic data usage on the accuracy and rate of the scenario reduction process focusing field development purposes. In this work, a synthetic benchmark case named UNISIM-II-D considering the flow unit modelling is used. The results showed both scenario reduction methods are reliable in selecting the RMs from a specific production strategy. However, the obtained RMs from a defined strategy using the DCSMC method can be applied to other strategies preserving the representativeness of the models, while the role of the strategy types to select the RMs using the metaheuristic method is substantial so that each strategy has its own set of RMs. Due to the field development workflow in which the metaheuristic algorithm is used, the number of required flow simulation models and the computational time are greater than the workflow in which the DCSMC method is applied. Hence, it can be concluded that static reservoir data usage on the scenario reduction process can be more reliable during the field development phase.


Author(s):  
Jian Liu ◽  
Yong Yu ◽  
Chenqi Zhu ◽  
Yu Zhang

The finite volume method (FVM)-based computational fluid dynamics (CFD) technology has been applied in the non-invasive diagnosis of coronary artery stenosis. Nonetheless, FVM is a time-consuming process. In addition to FVM, the lattice Boltzmann method (LBM) is used in fluid flow simulation. Unlike FVM solving the Navier–Stokes equations, LBM directly solves the simplified Boltzmann equation, thus saving computational time. In this study, 12 patients with left anterior descending (LAD) stenosis, diagnosed by CTA, are analysed using FVM and LBM. The velocities, pressures, and wall shear stress (WSS) predicted using FVM and LBM for each patient is compared. In particular, the ratio of the average and maximum speed at the stenotic part characterising the degree of stenosis is compared. Finally, the golden standard of LAD stenosis, invasive fractional flow reserve (FFR), is applied to justify the simulation results. Our results show that LBM and FVM are consistent in blood flow simulation. In the region with a high degree of stenosis, the local flow patterns in those two solvers are slightly different, resulting in minor differences in local WSS estimation and blood speed ratio estimation. Notably, these differences do not result in an inconsistent estimation. Comparison with invasive FFR shows that, in most cases, the non-invasive diagnosis is consistent with FFR measurements. However, in some cases, the non-invasive diagnosis either underestimates or overestimates the degree of stenosis. This deviation is caused by the difference between physiological and simulation conditions that remains the biggest challenge faced by all CFD-based non-invasive diagnostic methods.


Author(s):  
Nor Eddine Laghzale ◽  
Abdel-Hakim Bouzid

Steam generators are the subject of major concern in nuclear power plant safety. Within these generators, in addition to the structural integrity, the gross tightness barrier, which separates the primary and secondary circuits, is primarily ensured by the presence of a residual contact pressure at the tube-to-tubesheet joint interface. Any leakage is unacceptable, and its consequences are very heavy in terms of the human and environmental safety as well as maintenance cost. Some studies have been conducted to understand the main reasons for such a failure. However, no analytical model able to predict the attenuation of the residual contact pressure under the effect of material creep relaxation behavior. The development of a simple analytical model able to predict the change of the residual contact pressure as a function of time is laid out in this paper. The results from the analytical model are checked and compared with those of finite elements.


2007 ◽  
Vol 30 (7) ◽  
pp. 640-648 ◽  
Author(s):  
R. Kaminsky ◽  
K. Dumont ◽  
H. Weber ◽  
M. Schroll ◽  
P. Verdonck

The aim of this study was to validate the 2D computational fluid dynamics (CFD) results of a moving heart valve based on a fluid-structure interaction (FSI) algorithm with experimental measurements. Firstly, a pulsatile laminar flow through a monoleaflet valve model with a stiff leaflet was visualized by means of Particle Image Velocimetry (PIV). The inflow data sets were applied to a CFD simulation including blood-leaflet interaction. The measurement section with a fixed leaflet was enclosed into a standard mock loop in series with a Harvard Apparatus Pulsatile Blood Pump, a compliance chamber and a reservoir. Standard 2D PIV measurements were made at a frequency of 60 bpm. Average velocity magnitude results of 36 phase-locked measurements were evaluated at every 10° of the pump cycle. For the CFD flow simulation, a commercially available package from Fluent Inc. was used in combination with in-house developed FSI code based on the Arbitrary Lagrangian-Eulerian (ALE) method. Then the CFD code was applied to the leaflet to quantify the shear stress on it. Generally, the CFD results are in agreement with the PIV evaluated data in major flow regions, thereby validating the FSI simulation of a monoleaflet valve with a flexible leaflet. The applicability of the new CFD code for quantifying the shear stress on a flexible leaflet is thus demonstrated. (Int J Artif Organs 2007; 30: 640–8)


Author(s):  
Wei Li ◽  
Daoming Liu ◽  
Mathieu Desbrun ◽  
Jin Huang ◽  
Xiaopei Liu

Author(s):  
Guomin Ji ◽  
Nabila Berchiche ◽  
Sébastien Fouques ◽  
Thomas Sauder ◽  
Svein-Arne Reinholdtsen

The paper addresses the structural integrity assessment of lifeboat launched from floating production, storage and offloading (FPSO) vessels. The study is based on long-term drop lifeboat simulations accounting for more than 50 years of hindcast data of metocean conditions and corresponding FPSO motions. Selection of the load cases and strength analyses with high computational time is a challenge. The load cases analyzed are those corresponding to the 99th percentile of long term distribution of indicators for large slamming loads (CARXZ) or large submergence (Imaxsub). For six selected cases, the time-varying pressure distribution on the lifeboat hull during and after water impact is calculated by CFD simulations using StarCCM+. The finite element model (FEM) of the composite structure of the lifeboat is modelled by ABAQUS. Quasi-static finite element (FE) analyses are performed for the selected load cases. The structural integrity is assessed by the maximum stress and Tsai-Wu failure measure. In the present study, the load and resistance factors are combined and applied to the response. A sensitivity study is performed to investigate the non-linear load/response effects when the load factor is applied to the load. In addition, dynamic analysis is performed with the time-varying pressure distribution for selected case and the dynamic effect is investigated.


2021 ◽  
Author(s):  
Rafael M. D. Rosa ◽  
Arthur B. Soprana ◽  
Vinicius Girardi ◽  
Fernando M. Villagra

Abstract This work presents a numerical assessment of chemical inhibitor injection to mitigate wax deposition in unconventional wells. The goal of this study is to simulate the deposition of wax under several operational conditions and later optimize the chemical inhibitor injection position, using two different types of numerical simulations. A transient one-dimensional multiphase flow simulator - ALFAsim, with a dedicated wax model, was used to predict flow conditions such pressure, temperature, holdup and flow pattern profiles, as well the position and rates that wax accumulates. The results from the 1D simulation were then used as boundary conditions in a 3D CFD simulator, which aimed to assess how long it would take to a satisfactory homogenization of the inhibitor with the flow and what would be the minimum depth for the injector should be installed. In this work, a 1D multiphase flow simulator with wax deposition model was used to identify on which operational conditions (flow rates and environmental temperatures) an unconventional well would start to present wax deposition on its tubing walls. After defining the susceptible region where the paraffin could deposit, it was important to verify if the inhibitor would be well homogenized with the stream when reaching this region. For that, a 3D CFD simulation was performed, using information obtained directly from the 1D simulator as boundary conditions. The CFD model was capable to show the mixing evolution of the inhibitor with the stream and it was possible to determine the minimum distance where the injector should be placed to guarantee such homogeneity. A real well was selected to provide comparisons between field observations and simulated data, in order to validate the model assumptions and accuracy.


Author(s):  
Shoichi Yoshida

Floating roofs are widely used to prevent evaporation of content in large cylindrical aboveground oil storage tanks. The 2003 Hokkaido Earthquake caused severe damages to the floating roofs due to sloshing. These accidents became a cause to establish structural integrity of the floating roof tanks in sloshing. However, many designers do not have a solution for the sloshing of floating roof tanks except for three-dimensional FEA computer codes. The three-dimensional FEA requires a long computational time and expenses. The sloshing of floating roof tanks is a coupling vibration problem with fluid and structure. The simplified and convenient method has been desired for this solution. This paper presents a simplified development method of a FEA code in the axisymmetric linear problem. It is performed to modify an existing structural analysis code. The fluid behavior is formulated in terms of displacement as the Lagrangian approach.


Author(s):  
Mrinalgouda Patil ◽  
Anubhav Datta

A time-parallel algorithm is developed for large-scale three-dimensional rotor dynamic analysis. A modified harmonic balance method with a scalable skyline solver forms the kernel of this algorithm. The algorithm is equipped with a solution procedure suitable for large-scale structures that have lightly damped modes near resonance. The algorithm is integrated in X3D, implemented on a hybrid shared and distributed memory architecture, and demonstrated on a three-dimensional structural model of a UH-60A-like fully articulated rotor. Flight-test data from UH-60A Airloads Program transition flight C8513 are used for validation. The key conclusion is that the new solver converges to the time marching solution more than 50 times faster and achieves a performance greater than 1 teraFLOPS. The significance of this conclusion is that the principal barrier of computational time for trim solution using high-fidelity three-dimensional structures can be overcome with the scalable harmonic balance method demonstrated in this paper.


Sign in / Sign up

Export Citation Format

Share Document