Parallel Programming Of A Reservoir Simulator

1992 ◽  
pp. 1-13
Author(s):  
Mariyamni Awang
Keyword(s):  
Parallel Programming ◽  
Gaussian Elimination ◽  
Parallel Computer ◽  
Two Phase ◽  
Fully Implicit ◽  
The Matrix ◽  
Speed Up ◽  
Black Oil Model ◽  
Reservoir Simulator ◽  
Raphson Method

This study concerns applying parallel programming to reservoir simulation using a 32-Mbyte, 12-processor parallel computer. The effects of number of processes, granularity, load balancing and program structure were studied. The model simulated was a two-dimensionals, two-phase, black oil model with a fully-implicit formulation. The differenced equations were solved by the Newton-Raphson method and, Gaussian elimination was used to solve the Jacobian matrix. Matrix generation was parallelized using monitors as macros to synchronize calculation. The performance of the simulator was measured by the speed up. The speed ups of the matrix generation time increased almost linearly with increasing number of processes. For all of the models tested, the speed ups ranged from 3.5 to 4.0 for four processes and 7.0 to 7.9 for eight proceses.

Two-Phase Flow (Oil-Water) in Petroleum Reservoir with Irregular Geometry Including Water Injection: Effect of Porosity on the Oil Recovery Factor

2012 ◽  
Vol 326-328 ◽  
pp. 181-186 ◽  
Author(s):  
F. Alves Batista ◽  
B. Gonçalves Coutinho ◽  
Severino Rodrigues de Farias Neto ◽  
Antônio Gilson Barbosa de Lima
Keyword(s):  
Oil Recovery ◽  
Two Phase Flow ◽  
Water Injection ◽  
Recovery Factor ◽  
Phase Flow ◽  
Petroleum Reservoir ◽  
Two Phase ◽  
Fully Implicit ◽  
Irregular Geometry ◽  
Black Oil Model

The aim of this work is to study theoretically the effect of porosity of an oil reservoir with arbitrary geometry on the oil recovery factor. A two-dimensional mathematical modeling (Black-oil model) and numerical solution applied to two-phase flow (water-oil) into the reservoir with irregular geometry including water injection is presented. The conservation equations written in generalized coordinates are solved using the finite volume method, with a fully implicit technique. Results of the pressure and saturation distributions and oil recovery factor over time are presented and evaluated for different values of porosity of the reservoir.

Transition alpha structure in beta-isomorphous Nb-Hf alloys

1987 ◽  
Vol 45 ◽  
pp. 232-233
Author(s):  
R.W. Carpenter ◽  
Changhai Li ◽  
David J. Smith
Keyword(s):  
Solid Solution ◽  
Solution Phase ◽  
Miscibility Gap ◽  
Large Strains ◽  
Solid Solution Phase ◽  
Two Phase ◽  
Diffuse Intensity ◽  
Metastable Form ◽  
The Matrix ◽  
Equilibrium Morphology

Binary Nb-Hf alloys exhibit a wide bcc solid solution phase field at temperatures above the Hfα→ß transition (2023K) and a two phase bcc+hcp field at lower temperatures. The β solvus exhibits a small slope above about 1500K, suggesting the possible existence of a miscibility gap. An earlier investigation showed that two morphological forms of precipitate occur during the bcc→hcp transformation. The equilibrium morphology is rod-type with axes along <113> bcc. The crystallographic habit of the rod precipitate follows the Burgers relations: {110}||{0001}, <112> || <1010>. The earlier metastable form, transition α, occurs as thin discs with {100} habit. The {100} discs induce large strains in the matrix. Selected area diffraction examination of regions ∼2 microns in diameter containing many disc precipitates showed that, a diffuse intensity distribution whose symmetry resembled the distribution of equilibrium α Bragg spots was associated with the disc precipitate.

Using crystallographic symmetry in diffraction and contrast analysis

1989 ◽  
Vol 47 ◽  
pp. 504-505
Author(s):  
U. Dahmen ◽  
K.H. Westmacott
Keyword(s):  
Electron Microscopy ◽  
High Symmetry ◽  
Two Phase ◽  
Tem Analysis ◽  
Crystallographic Symmetry ◽  
The Matrix ◽  
Contrast Analysis ◽  
Practical Tool ◽  
Convergent Beam ◽  
Beam Diffraction

Despite the increased use of convergent beam diffraction, symmetry concepts in their more general form are not commonly applied as a practical tool in electron microscopy. Crystal symmetry provides an abundance of information that can be used to facilitate and improve the TEM analysis of crystalline solids. This paper draws attention to some aspects of symmetry that can be put to practical use in the analysis of structures and morphologies of two-phase materials.It has been shown that the symmetry of the matrix that relates different variants of a precipitate can be used to determine the axis of needle- or lath-shaped precipitates or the habit plane of plate-shaped precipitates. By tilting to a special high symmetry orientation of the matrix and by measuring angles between symmetry-related variants of the precipitate it is possible to find their habit from a single micrograph.

Microstructure of Polycrystalline MgO Penetrated by a Silicate Liquid

1996 ◽  
Vol 2 (3) ◽  
pp. 113-128 ◽  
Author(s):  
Sundar Ramamurthy ◽  
Michael P. Mallamaci ◽  
Catherine M. Zimmerman ◽  
C. Barry Carter ◽  
Peter R. Duncombe ◽  
...  
Keyword(s):  
Grain Boundaries ◽  
Grain Growth ◽  
Glassy Phase ◽  
Silicate Liquid ◽  
Two Phase ◽  
The Matrix ◽  
Devitrification Process ◽  
Phase Mixture

Dense, polycrystalline MgO was infiltrated with monticellite (CaMgSiO4) liquid to study the penetration of liquid along the grain boundaries of MgO. Grain growth was found to be restricted with increasing amounts of liquid. The inter-granular regions were generally found to be comprised of a two-phase mixture: crystalline monticellite and a glassy phase rich in the impurities present in the starting MgO material. MgO grains act as seeding agents for the crystallization of monticellite. The location and composition of the glassy phase with respect to the MgO grains emphasizes the role of intergranular liquid during the devitrification process in “snowplowing” impurities present in the matrix.

COMPUTATIONS OF PULSATILE AORTIC BLOOD FLOW PROBLEMS ON PARALLEL COMPUTERS

2003 ◽  
Vol 15 (03) ◽  
pp. 109-114
Author(s):  
YANG-YAO NIU ◽  
SHOU-CHENG TCHENG
Keyword(s):  
Blood Flow ◽  
Stokes Equations ◽  
Computing Time ◽  
Time Integration ◽  
Parallel Computer ◽  
Aortic Blood Flow ◽  
Pc Cluster ◽  
Flow Problems ◽  
Aortic Blood ◽  
Speed Up

In this study, a parallel computing technology is applied on the simulation of aortic blood flow problems. A third-order upwind flux extrapolation with a dual-time integration method based on artificial compressibility solver is used to solve the Navier-Stokes equations. The original FORTRAN code is converted to the MPI code and tested on a 64-CPU IBM SP2 parallel computer and a 32-node PC Cluster. The test results show that a significant reduction of computing time in running the model and a super-linear speed up rate is achieved up to 32 CPUs at PC cluster. The speed up rate is as high as 49 for using IBM SP2 64 processors. The test shows very promising potential of parallel processing to provide prompt simulation of the current aortic flow problems.

Failure Mechanisms of Particulate Two-Phase Composites

1998 ◽  
Vol 529 ◽  
Author(s):  
T. Antretter ◽  
E D. Fischer
Keyword(s):  
Residual Stresses ◽  
Crack Growth ◽  
Opening Angle ◽  
Geometrical Parameters ◽  
Two Phase ◽  
Absolute Size ◽  
The Matrix ◽  
Unit Cells ◽  
Angle Criterion ◽  
Probability Of Fracture

AbstractIn many composites consisting of hard and brittle inclusions embedded in a ductile matrix failure can be attributed to particle cleavage followed by ductile crack growth in the matrix. Both mechanisms are significantly sensitive towards the presence of residual stresses.On the one hand particle failure depends on the stress distribution inside the inclusion, which, in turn, is a function of various geometrical parameters such as the aspect ratio and the position relative to adjacent particles as well as the external load. On the other hand it has been observed that the absolute size of each particle plays a role as well and will, therefore, be taken into account in this work by means of the Weibull theory. Unit cells containing a number of quasi-randomly oriented elliptical inclusions serve as the basis for the finite element calculations. The numerical results are then correlated to the geometrical parameters defining the inclusions. The probability of fracture has been evaluated for a large number of inclusions and plotted versus the particle size. The parameters of the fitting curves to the resulting data points depend on the choice of the Weibull parameters.A crack tip opening angle criterion (CTOA) is used to describe crack growth in the matrix emanating from a broken particle. It turns out that the crack resistance of the matrix largely depends on the distance from an adjacent particle. Residual stresses due to quenching of the material tend to reduce the risk of particle cleavage but promote crack propagation in the matrix.

Study on Expander-Auxiliary Compressor Unit Used in Refrigeration System

2011 ◽  
Vol 354-355 ◽  
pp. 811-818
Author(s):  
Tao He ◽  
Zhi Yuan Wang ◽  
Yu Qing Xue
Keyword(s):  
Energy Recovery ◽  
Rotation Speed ◽  
Cooling Capacity ◽  
Compressor Unit ◽  
Performance Parameters ◽  
Refrigeration System ◽  
Two Phase ◽  
Speed Up ◽  
Testing Condition ◽  
Auxiliary Unit

The throttling valve used in the refrigeration system always causes energy loss. In this paper, an energy recovery device in the refrigeration system which was composed of an expander-auxiliary compressor unit to replace the throttling valve was investigated. On the basis of thermodynamic analysis, two typical arrangements which the auxiliary compressor was connected to the main compressor of the refrigeration system were compared and the system performance parameters were discussed. A prototype of an expander-auxiliary unit was manufactured to observe the expander performance using R410A as refrigerant. The results showed the reliability of the unit working in the two-phase flow condition with the rotation speed up to 21020 rpm. And the maximum increases in the cooling capacity by 3.9% and COP by 3.2% could be obtained under the testing condition.

scholarly journals MC2: a two-phase algorithm for leveraged matrix completion

2018 ◽  
Vol 7 (3) ◽  
pp. 581-604 ◽  
Author(s):  
Armin Eftekhari ◽  
Michael B Wakin ◽  
Rachel A Ward
Keyword(s):  
Broad Class ◽  
Matrix Completion ◽  
A Priori ◽  
Total Sample ◽  
Second Phase ◽  
Sample Complexity ◽  
Two Phase ◽  
The Matrix ◽  
Leverage Scores ◽  
Uniform Random Sampling

Abstract Leverage scores, loosely speaking, reflect the importance of the rows and columns of a matrix. Ideally, given the leverage scores of a rank-r matrix $M\in \mathbb{R}^{n\times n}$, that matrix can be reliably completed from just $O (rn\log ^{2}n )$ samples if the samples are chosen randomly from a non-uniform distribution induced by the leverage scores. In practice, however, the leverage scores are often unknown a priori. As such, the sample complexity in uniform matrix completion—using uniform random sampling—increases to $O(\eta (M)\cdot rn\log ^{2}n)$, where η(M) is the largest leverage score of M. In this paper, we propose a two-phase algorithm called MC2 for matrix completion: in the first phase, the leverage scores are estimated based on uniform random samples, and then in the second phase the matrix is resampled non-uniformly based on the estimated leverage scores and then completed. For well-conditioned matrices, the total sample complexity of MC2 is no worse than uniform matrix completion, and for certain classes of well-conditioned matrices—namely, reasonably coherent matrices whose leverage scores exhibit mild decay—MC2 requires substantially fewer samples. Numerical simulations suggest that the algorithm outperforms uniform matrix completion in a broad class of matrices and, in particular, is much less sensitive to the condition number than our theory currently requires.

MRHS solver based on linear algebra and exhaustive search

2018 ◽  
Vol 12 (3) ◽  
pp. 143-157 ◽  
Author(s):  
Håvard Raddum ◽  
Pavol Zajac
Keyword(s):  
Linear Algebra ◽  
Matrix Representation ◽  
Equation System ◽  
Exhaustive Search ◽  
Binary Matrix ◽  
Algebraic Attacks ◽  
The Matrix ◽  
Speed Up ◽  
Symmetric Key

Abstract We show how to build a binary matrix from the MRHS representation of a symmetric-key cipher. The matrix contains the cipher represented as an equation system and can be used to assess a cipher’s resistance against algebraic attacks. We give an algorithm for solving the system and compute its complexity. The complexity is normally close to exhaustive search on the variables representing the user-selected key. Finally, we show that for some variants of LowMC, the joined MRHS matrix representation can be used to speed up regular encryption in addition to exhaustive key search.

Sign in / Sign up

Export Citation Format

Share Document