|
[Sponsors] |
[waves2Foam] Surface elevation sigsegv on "sampledSurfaceElevation::sampleIntegrateAndWr ite" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 8, 2016, 07:02 |
Surface elevation sigsegv on "sampledSurfaceElevation::sampleIntegrateAndWr ite"
|
#1 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 10 |
Hi, I'm trying to use the surface elevation utility and probes to calculate velocities but I'm getting this error after a few time steps:
Code:
smoothSolver: Solving for alpha.water, Initial residual = 0.000163258, Final residual = 1.37002e-13, No Iterations 2 Phase-1 volume fraction = 0.499709 Min(alpha1) = 0 Max(alpha1) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.499709 Min(alpha1) = -1.57355e-14 Max(alpha1) = 1 Relaxing time: 0.52 s DILUPBiCG: Solving for Ux, Initial residual = 0.0367039, Final residual = 4.33518e-12, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.182004, Final residual = 5.89471e-11, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.0629783, Final residual = 1.61046e-11, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.0738539, Final residual = 5.28775e-08, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 9.70499e-06, Final residual = 5.86964e-08, No Iterations 3 time step continuity errors : sum local = 2.27865e-08, global = 6.40584e-09, cumulative = 6.33201e-08 GAMG: Solving for p_rgh, Initial residual = 0.000241955, Final residual = 7.58907e-08, No Iterations 4 GAMG: Solving for p_rgh, Initial residual = 3.12518e-07, Final residual = 3.94517e-08, No Iterations 1 time step continuity errors : sum local = 1.5316e-08, global = 5.29587e-09, cumulative = 6.86159e-08 GAMG: Solving for p_rgh, Initial residual = 2.65494e-06, Final residual = 5.06483e-08, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 5.31601e-08, Final residual = 9.61977e-09, No Iterations 1 time step continuity errors : sum local = 3.73531e-09, global = 3.0406e-09, cumulative = 7.16565e-08 ExecutionTime = 73.44 s ClockTime = 74 s #0 Foam::error::printStack(Foam::Ostream&) in "/home/ariel/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigSegv::sigHandler(int) in "/home/ariel/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::sampledSurfaceElevation::sampleIntegrateAndWrite(Foam::sampledSurfaceElevation::fieldGroup<double>&) in "/home/ariel/OpenFOAM/ariel-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2FoamSampling.so" #4 Foam::OutputFilterFunctionObject<Foam::sampledSurfaceElevation>::execute(bool) in "/home/ariel/OpenFOAM/ariel-2.3.0/platforms/linux64GccDPOpt/lib/libwaves2FoamSampling.so" #5 Foam::functionObjectList::execute(bool) in "/home/ariel/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #6 Foam::Time::run() const in "/home/ariel/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #7 in "/home/ariel/OpenFOAM/ariel-2.3.0/platforms/linux64GccDPOpt/bin/waveFoam" #8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #9 in "/home/ariel/OpenFOAM/ariel-2.3.0/platforms/linux64GccDPOpt/bin/waveFoam" Segmentation fault Code:
deleteParentOutputDirectory false; regularSpectrum { callName surfaceElevation; removeDuplicate true; inputDir surfaceElevationAnyName; deltaT 0.01; tMin 10.0; //0; Default: minimum time //tMax 20; // Default: maximum time // The action list is a set of processes on the data actionList ( interpolateSurfaceElevation powerSpectraLS write2Ascii removeData ); nFreq 4; period 3; allDataSets true; } Code:
surfaceElevationAnyName { type waveGauge; pointDistribution lineDistribution; N 11; //xValues nonuniform List<scalar> 11(-3.0 -2.0 -1 -0.5 0.5 0.75 1.5 5.0 10.0 15.0 20.0); //yValues uniform 0; //zValues uniform 0.05; linestart (-1 0 0.05); lineend (3 0 0.05); add (0 0 0.45); //( 0 0.45 0); axis z; stretch 1.5; } Does anyone know what this error means or where it's coming from? |
|
January 19, 2016, 12:16 |
|
#2 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Good afternoon all,
I will try to answer some of the questions: @Ariel #1135: It is quite likely that your wave gauges are outside of your computational domain. 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 08:34. Reason: removed answers to other posts that were on the main thread |
|
January 29, 2016, 11:04 |
|
#3 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 10 |
Niels.
thank you for your response. I am now trying to fix the location of the wave gauges but am very confused by what is happening with this in the waveFlume tutorial: First of all, in this line in the controlDict file: Code:
#includeIfPresent "../waveGaugesNProbes/surfaceElevationAnyName_controlDict"; Code:
constant/probeDefinitions Code:
surfaceElevationAnyName { type waveGauge; pointDistribution lineDistribution; N 179; linestart ( 0.10 -0.3 0.05); lineend ( 17.9 -0.3 -0.05); add ( 0 0.45 0); axis y; stretch 1.; } 1) What are the add and stretch values? 2) If I want to interpolate for surface elevation at several different areas in my wave tank, say to look for surface elevation decay, how would I set up more than one lineDistribution or can I do this in the way that I would do it for a system/probeDict file? Thanks a lot for any advice, Ariel |
|
January 29, 2016, 15:49 |
|
#4 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Hi Ariel,
Note that there are generated VTK-files that show the locations of your wave gauges, when you execute waveGaugesNProbes. This makes it trivial to explore the effects of various parameters in paraView. 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. |
|
February 1, 2016, 08:33 |
|
#5 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 10 |
Hi Niels,
Thank you for your response. I am able to view the probes but am still not understanding how to manipulate some of the parameters. I have tried to turn on and off and change all features but just haven't been able to set up what I want for some reason. I am trying to calculate the surface elevation at different locations down a wave tank, I think just like the waveFlume tutorial. My wave tank is 600 meters long, the wavelength is 93 m, the depth is 10 m and gravity acts in the z direction. When I view the vtk file in parafoam, it simply doesn't look correct...as if the probes aren't stretching from the bottom to the sea surface or their location doesn't seem right. Here is my constant/probeDefinitions file: Code:
surfaceElevationAnyName { type waveGauge; pointDistribution lineDistribution; N 599; linestart ( 0.10 -3 -10); lineend (599.9 3 4); add ( 0 0 0.5); axis z; stretch 1.; } Thank you EDIT: Ok I seem to have correctly placed the probes but it was almost by luck.. - Can you please explain the "add () " parameter? - I messed around with the y-value and it didn't seem to make a difference... I don't understand why not? Last edited by arieljeds; February 1, 2016 at 08:46. Reason: Solved the problem but still not fully understanding |
|
February 1, 2016, 12:44 |
|
#6 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Hallo Ariel,
It is good that you got it working. First of all, allow me to say that you have full access to the source code, so it should be fairly straight forward to figure out, what the code it doing. That said, for the 'pointDistribution' of type 'lineDistribution', then the 'linestart', 'lineend', 'N' and 'stretch' define N points along a line (e.g. a horizontal line). The stretch parameter allows you to define a non-equidistant stretch of the points. The keyword 'axis' comes from the original sample utility and should be given as your vertical axis. The parameter 'add' is a vector (vertical), which is added to each of the points on the lineDistribution. The original and translated point (translated by 'add') define a vertical line that is your wave gauge. 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. |
|
February 1, 2016, 12:49 |
|
#7 |
Senior Member
ArielJ
Join Date: Aug 2015
Posts: 127
Rep Power: 10 |
Niels -
thank you very much for your clear explanation. I've been looking through the source code and did manage to get it working but was still struggling to really picture what was happening! Very silly, I know.. I do seem to have it working although I had a look at the free surface elevation and it's always a negative value. I will need to check where my probes are but I thought I had set them up to extend to 4 meters above the sea surface, which I expected would give negative and positive values. I am not looking at my code at the moment and will try and fix this tomorrow based on what you've said. Thanks again Ariel |
|
June 12, 2016, 10:48 |
|
#8 | |
New Member
Muhammad Nouman Khalid
Join Date: Mar 2016
Posts: 6
Rep Power: 10 |
Quote:
Hope you are doing fine. I am having the same problem as yours with probe definition. My gauges are well within the computational domain as i see them in paraFoam. but when I execute postProcessWaves2Foam it shows the same error as yours. how was you able to fix it? I would really appreciate your help. Regards Khalid |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Gmsh] Problem with Gmsh | nishant_hull | OpenFOAM Meshing & Mesh Conversion | 23 | August 5, 2015 02:09 |
SSIIM- F 36 2: computation of the vertical water surface elevation | Mummputz | Main CFD Forum | 2 | July 1, 2015 04:52 |
chart of surface elevation against time | mohammad_1986 | CFX | 1 | January 17, 2013 20:43 |
STL units and Surface elevation plot | MAB | FLOW-3D | 11 | May 9, 2010 15:55 |
ploting surface elevation | Tony Campbell | Siemens | 5 | January 29, 2005 16:09 |