|
[Sponsors] |
[waves2Foam] surfaceElevation: "cell does not contain point" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 13, 2017, 07:55 |
surfaceElevation: "cell does not contain point"
|
#1 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 11 |
Hi all -
I have a frustrating problem using the surfaceElevation utility. I use this often without problems with exactly the same locations, same domain size, but I can't seem to figure out what's happening here. Is there a way to view the wave gauge locations using surfaceElevation? I know you can do it using waveGaugesNProbes but I couldn't figure out how to set up my surfaceElevation wave gauges using that. I'm getting the following error: Code:
Reading g Reading waveProperties --> FOAM FATAL ERROR: Found cell 1719049 using face 9977105. But cell does not contain point (3.04 3.04 -14.481) From function sampledSet::getCell(const label, const point&) in file sampledSet/sampledSet/sampledSet.C at line 76. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::sampledSet::getCell(int, Foam::Vector<double> const&) const at ??:? #3 Foam::midPointSet::genSamples() at ??:? #4 Foam::midPointSet::midPointSet(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch const&, Foam::dictionary const&) at ??:? #5 Foam::sampledSet::addwordConstructorToTable<Foam::midPointSet>::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch const&, Foam::dictionary const&) at ??:? #6 Foam::sampledSet::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch const&, Foam::dictionary const&) at ??:? #7 void Foam::PtrList<Foam::sampledSet>::read<Foam::sampledSet::iNew>(Foam::Istream&, Foam::sampledSet::iNew const&) at ??:? #8 Foam::sampledSurfaceElevation::read(Foam::dictionary const&) at ??:? #9 Foam::sampledSurfaceElevation::sampledSurfaceElevation(Foam::word const&, Foam::objectRegistry const&, Foam::dictionary const&, bool) at ??:? #10 main at ??:? #11 __libc_start_main in "/lib64/libc.so.6" #12 ? at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116 Aborted (core dumped) I know this has to do with there not being a cell where it's attempting to sample but there should definitely be a cell there?? I want to view where it's trying to sample because I can't figure out why it thinks there's nothing there. I am using 160 wave gauges around a cylinder so this is just a snippet of my surfaceElevationDict: Code:
setFormat raw; interpolationScheme cellPoint; // Fields to sample. fields ( alpha.water ); sets ( loc2_gauge_0 { type midPoint; axis z; start (4.30 0 -15); end (4.30 0 15); nPoints 100; } loc2_gauge_1 { type midPoint; axis z; start (3.97 1.65 -15); end (3.97 1.65 15); nPoints 100; } loc2_gauge_2 { type midPoint; axis z; start (3.04 3.04 -15); end (3.04 3.04 15); nPoints 100; } loc2_gauge_3 { type midPoint; axis z; start (1.65 3.97 -15); end (1.65 3.97 15); nPoints 100; } ); Thanks in advance for any help on this. Ariel |
|
January 15, 2017, 12:58 |
|
#2 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,901
Rep Power: 37 |
Good evening.
@Ariel: I am not familiar with this error. I can only recommend you to track down the problem by removing one wave gauge at the time. Kind regards, Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request. Last edited by wyldckat; August 25, 2018 at 07:05. Reason: removed answers to other posts that were on the main thread |
|
January 19, 2017, 06:14 |
|
#3 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 11 |
Hi Niels,
Thanks for the reply. Is there some way to set up the waveGaugesNProbes using sets? I have my probes in pretty specific locations and I can't figure out how to do that in probeDefinitions (as in, can I define the sets within this file somehow instead of using some automated distribution?) Or is there a way to view the surfaceElevation probes using surfaceElevationDict (which I know I can do if I run waveGaugesNProbes)? I tried using userDefinedDistribution but then I got the error: Code:
"ill defined primitiveEntry starting at keyword 'sets' on line 25... Ariel **EDIT** Ok so I have fixed my problem but I still have two queries: 1) I don't understand why it was fixed, the changes I made were something like 5.06966 to 5.07, and this fixed it. I really don't understand how the tiny change would have had such an effect? Also note, I have used the same surfaceElevationDict with a range of simulations with the same computational domain size, although a slightly different mesh resolution (this case is slightly coarser than others) 2) I would still like to find a way to view my wave gauges (create a vtk file with surfaceElevationDict) Last edited by arieljeds; January 19, 2017 at 06:57. Reason: Fixed problem, although don't understand how |
|
January 19, 2017, 06:16 |
|
#4 | |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 11 |
Quote:
I found that if you shrink it too much, the wave doesn't get fully absorbed and you will see some reflection. I played around with the mu value (alphaCoefficient in relaxation zone definition) and lowered it to 1 and was able to reduce the relaxation zone to just under 1 wavelength without getting reflection. Hope that helps! |
||
January 19, 2017, 15:53 |
|
#5 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,901
Rep Power: 37 |
Hi Ariels,
I have no idea of the cause of the error. It is most likely related to the sampling function in your version of OpenFoam. An error similar to the reason that pointInMesh for snappyHexMesh is not allowed on grid lines. You can automatically create vtk-files, if you use the utility waveGaugesNProbes based on input in constant/probeDefinitions. Kind regards Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request. |
|
January 23, 2017, 10:44 |
|
#6 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 11 |
Hi Niels,
Thanks for your response. Yes I usually would use the waveGaugesNProbes to view the wave gauges but I could not figure out how to set up the locations of the probes in the same way I'm doing that in my surfaceElevationDict (I don't want them automatically set up). Is there a way to create the vtk file on its own or using surfaceElevation? Or else a better way to set up the probeDefinitions file to specify locations? Best, Ariel |
|
January 25, 2017, 13:35 |
|
#7 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,901
Rep Power: 37 |
Hi Ariel,
There is a method already in waveGaugesNProbes, where you can provide user-defined location. I think it is called userDefinedDistribution instead of lineDistribution. I do not have access to the source code right now, so cannot check, but it should be easy for you. A VTK-file will be generated, if you use this approach. Kind regards, Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Fluent UDF wrong number of cells in parallel - correct in serial | dralexpe | Fluent UDF and Scheme Programming | 7 | May 17, 2018 08:26 |
[snappyHexMesh] snappyHexMesh and cyclic boundaries | Ruli | OpenFOAM Meshing & Mesh Conversion | 2 | December 9, 2013 06:51 |
[blockMesh] error EOF in blockMesh | Ahmed Khattab | OpenFOAM Meshing & Mesh Conversion | 7 | May 17, 2012 00:37 |
Warning 097- | AB | Siemens | 6 | November 15, 2004 04:41 |
CFX4.3 -build analysis form | Chie Min | CFX | 5 | July 12, 2001 23:19 |