CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Segmentation fault running waveDyMFoam solver (mod. interDyMFoam solver - waves2Foam)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 3, 2012, 08:33
Default Segmentation fault running waveDyMFoam solver (mod. interDyMFoam solver - waves2Foam)
  #1
New Member
 
Ed Ransley
Join Date: Jul 2012
Posts: 29
Rep Power: 5
Ed R is on a distinguished road
Hello,

This is my first ever forum post, so I apologise in advance for any ill forum etiquette.

I'm working on freak wave modelling and I'm using the waves2Foam library and OpenFoam 2.0.1. I have successful installed waves2Foam and compiled the solvers waveFoam and waveDyMFoam using the wiki instructions.

I've run many successful simulations using the waveFoam solver and I'm very happy with the waves2Foam library.

I'm now trying to use the waveDyMFoam solver to have a refined mesh at the free surface hoping to improve execution time for convergence. However, the solver crashes the first time it tries to refine the mesh.

I have a simple 2D wave tank with a solitary wave positioned in the centre at t=0s (Using the SolitaryFirst wave definition in my wavefoam file), it then moves from left to right and runs up on a fixed wall.

I have used the dynamicMeshDict file from the interDyMFoam tutorial for damBreakWithObstacle.

Here's what my log.waveDyMFoam file says after having read in the initial properties and selecting various transport models:

PIMPLE: Operating solver in PISO mode

time step continuity errors : sum local = 0.000158174, global = -2.98256e-09, cumulative = -2.98256e-09
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 9.09913e-08, No Iterations 18
GAMG: Solving for pcorr, Initial residual = 4.54879e-08, Final residual = 4.54879e-08, No Iterations 0
time step continuity errors : sum local = 1.43925e-11, global = 2.27339e-14, cumulative = -2.98254e-09
Courant Number mean: 0.00105706 max: 0.00636453

Starting time loop

Interface Courant Number mean: 0.000151212 max: 0.00623973
Courant Number mean: 0.00103634 max: 0.00623973
deltaT = 0.00116279
Time = 0.00116279

Selected 280 cells for refinement out of 1000.
Refined from 1000 to 2960 cells.
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::dynamicRefineFvMesh::refine(Foam::List<int> const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#4 Foam::dynamicRefineFvMesh::update() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#5
in "/home/cerg101/OpenFOAM/cerg101-2.0.1/platforms/linux64GccDPOpt/bin/waveDyMFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7
in "/home/cerg101/OpenFOAM/cerg101-2.0.1/platforms/linux64GccDPOpt/bin/waveDyMFoam"

Does anyone have any suggestions on how to resolve this problem?

Thanks,

Ed
Ed R is offline   Reply With Quote

Old   June 28, 2013, 15:32
Default Similary Problems
  #2
New Member
 
Paul B
Join Date: Sep 2012
Posts: 9
Rep Power: 4
pbarc is on a distinguished road
Hey Ed,

Did you ever figure out what was wrong with your simulation, because I am having similar problems with interDyMFoam.

Thanks,
Paul
pbarc is offline   Reply With Quote

Old   July 2, 2013, 05:55
Default
  #3
New Member
 
Ed Ransley
Join Date: Jul 2012
Posts: 29
Rep Power: 5
Ed R is on a distinguished road
Hi Paul,

Apologies for the delay, I was working with dynamic mesh a while ago and stopped as I didn't like how the initial conditions had to be applied to a coarse (unrefined) mesh. I would very much like to get the initial condions (setFields) to be reapplied after one refinement.

Anyway, YES I did resolve this problem and I have just now re-learnt what I did.

I believe the trick is to change cacheAgglomeration (in system/fvSolution) to false. (I found this on another forum post which I can't remember).

If that doesn't work. I think I might have changed the correctFluxes part of the constant/dynamicMeshDict. A working case of mine as this:

correctFluxes
(
(phi Urel)
(phiAbs U)
(phiAbs_0 U_0)
(nHatf none)
(rho*phi none)
(ghf none)
);

I hope this helps, let me now if this actually fixes your problem.

Regards,

Ed
Ed R is offline   Reply With Quote

Old   July 2, 2013, 09:39
Default RE
  #4
New Member
 
Paul B
Join Date: Sep 2012
Posts: 9
Rep Power: 4
pbarc is on a distinguished road
Thanks for the response Ed.

I made both the corrections that you suggested, but am still getting the same error for my case.

I am basically changing the damBreakWithObstacle case to a 2D case. so I shrunk the domain in one direction and applied empty bcs. However if I run normal interFoam on the case it runs just fine.

Do you have any ideas as to why interFoam would run fine, but interDyMFoam would crash after the first step?

I have uploaded my input files if you would like to look over them.

Thanks again,
Paul
Attached Files
File Type: zip 0.org.zip (2.0 KB, 8 views)
File Type: zip constant.zip (48.1 KB, 9 views)
File Type: zip system.zip (3.4 KB, 6 views)
File Type: txt log_InterDyMFoam.txt (2.0 KB, 8 views)
pbarc is offline   Reply With Quote

Old   July 2, 2013, 09:52
Default
  #5
New Member
 
Ed Ransley
Join Date: Jul 2012
Posts: 29
Rep Power: 5
Ed R is on a distinguished road
empty boundary conditions won't work with the dynamic mesh is needs the mesh to be 3D. You'll need to change these to something else. I hope that fixes it.

This is because the refinement works in a 3-dimension quadtree fashion where each cell selected for refinement is split in half in all three dimension turning one cell into 8. The case definitely must be 3D. I started with a 2D numerical wave tank but defined the front and back as symmetryPlane (seams to work if the case is symmetrical).

This is probably the reason for the error and so my other suggestions may be the answer to a different problem.

Here's the forum post regarding the cacheAgglomeration fix to a different error you may come across:
GAMG problem in interDyMFoam?

Ed
Ed R is offline   Reply With Quote

Old   July 2, 2013, 11:36
Default
  #6
New Member
 
Paul B
Join Date: Sep 2012
Posts: 9
Rep Power: 4
pbarc is on a distinguished road
Thanks for the help and info on how interDyMFoam works.

I was able to successfully run the case using the symmetry BC after following the GAMG corrections from the other forum.

Paul
pbarc is offline   Reply With Quote

Reply

Tags
dynamic mesh, solver error, wavedymfoam, waves2foam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Segmentation fault when running dieselFoam or dieselEngineFoam in parallel francesco OpenFOAM Bugs 3 June 10, 2011 09:30
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 17:55
interDyMFoam crashes with Segmentation Fault rassilon OpenFOAM 5 January 5, 2011 20:22
Statically Compiling OpenFOAM Issues herzfeldd OpenFOAM Installation 21 January 6, 2009 10:38


All times are GMT -4. The time now is 19:13.