scholarly journals An efficient multigrid-DEIM semi-reduced-order model for simulation of single-phase compressible flow in porous media

2020 ◽  
Author(s):  
Jing-Fa Li ◽  
Bo Yu ◽  
Dao-Bing Wang ◽  
Shu-Yu Sun ◽  
Dong-Liang Sun

Abstract In this paper, an efficient multigrid-DEIM semi-reduced-order model is developed to accelerate the simulation of unsteady single-phase compressible flow in porous media. The cornerstone of the proposed model is that the full approximate storage multigrid method is used to accelerate the solution of flow equation in original full-order space, and the discrete empirical interpolation method (DEIM) is applied to speed up the solution of Peng–Robinson equation of state in reduced-order subspace. The multigrid-DEIM semi-reduced-order model combines the computation both in full-order space and in reduced-order subspace, which not only preserves good prediction accuracy of full-order model, but also gains dramatic computational acceleration by multigrid and DEIM. Numerical performances including accuracy and acceleration of the proposed model are carefully evaluated by comparing with that of the standard semi-implicit method. In addition, the selection of interpolation points for constructing the low-dimensional subspace for solving the Peng–Robinson equation of state is demonstrated and carried out in detail. Comparison results indicate that the multigrid-DEIM semi-reduced-order model can speed up the simulation substantially at the same time preserve good computational accuracy with negligible errors. The general acceleration is up to 50–60 times faster than that of standard semi-implicit method in two-dimensional simulations, but the average relative errors of numerical results between these two methods only have the order of magnitude 10−4–10−6%.

2009 ◽  
Vol 131 (10) ◽  
Author(s):  
M. M. Awad ◽  
S. D. Butt

A simple semitheoretical method for calculating the two-phase frictional pressure gradient in porous media using asymptotic analysis is presented. The two-phase frictional pressure gradient is expressed in terms of the asymptotic single-phase frictional pressure gradients for liquid and gas flowing alone. In the present model, the two-phase frictional pressure gradient for x≅0 is nearly identical to the single-phase liquid frictional pressure gradient. Also, the two-phase frictional pressure gradient for x≅1 is nearly identical to the single-phase gas frictional pressure gradient. The proposed model can be transformed into either a two-phase frictional multiplier for liquid flowing alone (ϕl2) or a two-phase frictional multiplier for gas flowing alone (ϕg2) as a function of the Lockhart–Martinelli parameter X. The advantage of the new model is that it has only one fitting parameter (p), while the other existing correlations, such as the correlation of Larkins et al., Sato et al., and Goto and Gaspillo, have three constants. Therefore, calibration of the new model to the experimental data is greatly simplified. The new model is able to model the existing multiparameter correlations by fitting the single parameter p. Specifically, p=1/3.25 for the correlation of Midoux et al., p=1/3.25 for the correlation of Rao et al., p=1/3.5 for the Tosun correlation, p=1/3.25 for the correlation of Larkins et al., p=1/3.75 for the correlation of Sato et al., and p=1/3.5 for the Goto and Gaspillo correlation.


Open Physics ◽  
2016 ◽  
Vol 14 (1) ◽  
pp. 588-601 ◽  
Author(s):  
Yi Wang ◽  
Bo Yu ◽  
Shuyu Sun

AbstractFast prediction modeling via proper orthogonal decomposition method combined with Galerkin projection is applied to incompressible single-phase fluid flow in porous media. Cases for different configurations of porous media, boundary conditions and problem scales are designed to examine the fidelity and robustness of the model. High precision (relative deviation 1.0 × 10−4% ~ 2.3 × 10−1%) and large acceleration (speed-up 880 ~ 98454 times) of POD model are found in these cases. Moreover, the computational time of POD model is quite insensitive to the complexity of problems. These results indicate POD model is especially suitable for large-scale complex problems in engineering.


Author(s):  
M. M. Awad ◽  
S. D. Butt

A simple semi-theoretical method for calculating two-phase frictional pressure gradient in porous media using asymptotic analysis is presented. Two-phase frictional pressure gradient is expressed in terms of the asymptotic single-phase frictional pressure gradients for liquid and gas flowing alone. In the present model, the two-phase frictional pressure gradient for x ≅ 0 is nearly identical to single-phase liquid frictional pressure gradient. Also, the two-phase frictional pressure gradient for x ≅ 1 is nearly identical to single-phase gas frictional pressure gradient. The proposed model can be transformed into either a two-phase frictional multiplier for liquid flowing alone (φl2) or two-phase frictional multiplier for gas flowing alone (φg2) as a function of the Lockhart-Martinelli parameter, X. The advantage of the new model is that it has only one fitting parameter (p) while the other existing correlations such as Larkins et al. correlation, Sato et al. correlation, and Goto and Gaspillo correlation have three constants. Therefore, calibration of the new model to experimental data is greatly simplified. The new model is able to model the existing multi parameters correlations by fitting the single parameter p. Specifically, p = 1/3.25 for Midoux et al. correlation, p = 1/3.25 for Rao et al. correlation, p = 1/3.5 for Tosun correlation, p = 1/3.25 for Larkins et al. correlation, p = 1/3.75 for Sato et al. correlation, and p = 1/3.5 for Goto and Gaspillo correlation.


PLoS ONE ◽  
2021 ◽  
Vol 16 (10) ◽  
pp. e0258047
Author(s):  
Konstantinos G. Lyras ◽  
Jack Lee

Quantification of pressure drop across stenotic arteries is a major element in the functional assessment of occlusive arterial disease. Accurate estimation of the pressure drop with a numerical model allows the calculation of Fractional Flow Reserve (FFR), which is a haemodynamic index employed for guiding coronary revascularisation. Its non-invasive evaluation would contribute to safer and cost-effective diseases management. In this work, we propose a new formulation of a reduced-order model of trans-stenotic pressure drop, based on a consistent theoretical analysis of the Navier-Stokes equation. The new formulation features a novel term that characterises the contribution of turbulence effect to pressure loss. Results from three-dimensional computational fluid dynamics (CFD) showed that the proposed model produces predictions that are significantly more accurate than the existing reduced-order models, for large and small symmetric and eccentric stenoses, covering mild to severe area reductions. FFR calculations based on the proposed model produced zero classification error for three classes comprising positive (≤ 0.75), negative (≥ 0.8) and intermediate (0.75 − 0.8) classes.


2019 ◽  
Vol 24 (1) ◽  
pp. 131-142 ◽  
Author(s):  
E. Tengs ◽  
F. Charrassier ◽  
M. Holst ◽  
Pål-Tore Storli

Abstract As part of an ongoing study into hydropower runner failure, a submerged, vibrating blade is investigated both experimentally and numerically. The numerical simulations performed are fully coupled acoustic-structural simulations in ANSYS Mechanical. In order to speed up the simulations, a model order reduction technique based on Krylov subspaces is implemented. This paper presents a comparison between the full ANSYS harmonic response and the reduced order model, and shows excellent agreement. The speedup factor obtained by using the reduced order model is shown to be between one and two orders of magnitude. The number of dimensions in the reduced subspace needed for accurate results is investigated, and confirms what is found in other studies on similar model order reduction applications. In addition, experimental results are available for validation, and show good match when not too far from the resonance peak.


2019 ◽  
Vol 29 (11) ◽  
pp. 4167-4204 ◽  
Author(s):  
Jingfa Li ◽  
Tao Zhang ◽  
Shuyu Sun ◽  
Bo Yu

Purpose This paper aims to present an efficient IMPES algorithm based on a global model order reduction method, proper orthogonal decomposition (POD), to achieve the fast solution and prediction of two-phase flows in porous media. Design/methodology/approach The key point of the proposed algorithm is to establish an accurate POD reduced-order model (ROM) for two-phase porous flows. To this end, two projection methods including projecting the original governing equations (Method I) and projecting the discrete form of original governing equations (Method II) are respectively applied to construct the POD-ROM, and their distinctions are compared and analyzed in detail. It is found the POD-ROM established by Method I is inapplicable to multiphase porous flows due to its failed introduction of fluid saturation and permeability that locate on the edge of grid cell, which would lead to unphysical results. Findings By using Method II, an efficient IMPES algorithm that can substantially speed up the simulation of two-phase porous flows is developed based on the POD-ROM. The computational efficiency and numerical accuracy of the proposed algorithm are validated through three numerical examples, and simulation results illustrate that the proposed algorithm displays satisfactory computational speed-up (one to two orders of magnitude) without sacrificing numerical accuracy obviously when comparing to the standard IMPES algorithm that without any acceleration technique. In addition, the determination of POD modes number, the relative errors of wetting phase pressure and saturation, and the influence of POD modes number on the overall performances of the proposed algorithm, are investigated. Originality/value 1. Two projection methods are applied to establish the POD-ROM for two-phase porous flows and their distinctions are analyzed. The reason why POD-ROM is difficult to be applied to multiphase porous flows is clarified firstly in this study. 2. A highly efficient IMPES algorithm based on the POD-ROM is proposed to accelerate the simulation of two-phase porous flows. 3. Satisfactory computational speed-up (one to two orders of magnitude) and prediction accuracy of the proposed algorithm are observed under different conditions.


Author(s):  
Joseph Bakarji ◽  
Khoudor Keniar ◽  
Mohammad Cheikh ◽  
Issam Lakkis

A reduced-order model of a Microfluidic Transistor is presented. The transistor is essentially a long micro channel between substrate and a membrane that is pressure actuated. The proposed model captures steady (DC) and small signal (AC) behavior of the device in a manner analogous to standard semiconductor transistor models. The model is based on steady and perturbed unsteady solutions of the conservation of mass and momentum, coupled with an elastic model for the membrane. To improve the accuracy and to enhance the range of validity, the model is enhanced by numerical simulations of the coupled fluid-structure problem. The model predicts dependence of the transconductance on the pressure differentials across the membrane and along the channel. The proposed model also investigates the impact of flow inertia, among other effects, on the dynamic behavior of the transistor.


Sign in / Sign up

Export Citation Format

Share Document