Comments (6)
Enhancing TOPS to support non-constant top surface requires us to define just what we mean by it. For a simple example: given a grid of two cells with varying TOPS, should the grid become (1) a stair-step grid, (2) a conforming grid with sloping top/bottom faces, (3) or a stair-step grid where the connectivities are (magically) the same as for the conforming grid? Since we have not decided what it should mean, we do not support it. If you have suggestions as to the handling/meaning of TOPS, I am listening...
For now, I suggest using DEPTHZ (giving pillar top z values) or COORDS/ZCORN (full corner-point grid).
For the particular case of SPE9, we have made a version of that (very simple) grid using COORDS/ZCORN, this version makes the assumption that what you want is option (2) -- a conforming grid with parallelepiped cells. Please send me an email if this is what you want.
from opm-core.
@atgeirr Just as you said, enhancing TOPS is an undefined problems which I don't think much more about. What I want is trying to use TOPS create a conforming grid with sloping top/bottom faces especial for SPE9. DEPTHZ may be a good choice. By the way, could you send me your SPE9 corner grid file? And usually, which tools is used to create corner grid mesh file? Thanks.
from opm-core.
Just a small remark. All definitions I've ever seen of the SPE9 grid use the DX
, DY
, and DZ
keywords in addition to TOPS
. Strictly speaking, these keywords will (probably, I haven't tested it) allow specifying a geometry that cannot possibly exist in three space dimensions. The format does, in principle, allow specifying
DX(i) = i
and similarly for DY
and DZ
.
Therefore, whenever I have a choice, I personally prefer the more FrontSim-like keywords DXV
, DYV
, DZV
and DEPTHZ
. In the case of SPE9, these would be
DXV
24*300
/
DYV
25*300
/
DZV
20 15 26 15 16 14 8
8 18 12 19 18 20 50 100
/
and then the DEPTHZ
would be the output of the following programme
#include <math.h>
#include <stddef.h>
#include <stdio.h>
int
main(void)
{
size_t i, j;
double x, dip;
dip = tan(10.0 / 180.0); /* 10 degree dip angle */
printf("DEPTHZ\n");
for (j = 0; j < 26; j++) {
for (i = 0; i < 25; i++) {
x = 300.0 * i;
printf(" %.10e", 9000.0 + x*dip);
if ((i + 1) % 5 == 0) {
printf("\n");
}
}
}
printf("/\n");
return 0;
}
from opm-core.
dip = tan(10.0 / 180.0); /* 10 degree dip angle */
Well, almost. Of course, the real value would be something along the lines of
dip = tan(1.745329251994329e-01); /* approximately pi/18 = 10*pi/180 */
Sorry about that.
from opm-core.
@bska Thank you for your comments, that's a another way for giving the grid file for SPE9, when you don't have the COORD/ZCORN mesh file. Thanks. :)
from opm-core.
There's been no activity on this issue for eight days and it appears to be resolved. I'll close the issue, but please reopen if the problem resurfaces.
from opm-core.
Related Issues (20)
- Some tests break after OPM/opm-parser#768 HOT 1
- Initialization and restart HOT 1
- Failing test_relpermdiagnostics HOT 4
- Dereferencing of end iterator happens in initStateEquil.hpp
- test_pinchprocessor is broken? HOT 3
- opm-core build fails HOT 4
- parameter's ouput does not show up in .CASE_NAME.DEBUG HOT 2
- Allow passing the size of the initial time step to adaptive time stepper HOT 1
- Runtime failure in "linsolver" test - umfpack not configured HOT 2
- Add restart test with no wells present.
- WellState::init has memory problems HOT 7
- Don't we need to keep track of (artificially) shut completions in the simulator for output? HOT 3
- Windows Installer (guide) failed HOT 3
- Compilation error due to not finding opm/output/data/Wells.hpp HOT 2
- Broken test HOT 3
- Multiple definitions of OPM_THROW HOT 1
- opm-core broken related to petsc HOT 2
- Some files missing? HOT 2
- Question: How OPM handles well sources in multi-cell wells? HOT 4
- Installing OPM in CentOS 8.1 (RPM) HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from opm-core.