View Issue Details

IDProjectCategoryView StatusLast Update
0001194OpenFOAMBugpublic2014-03-17 15:59
ReporterSamMallinson Assigned Towill  
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version10.04
Summary0001194: Run which completes normally using writeFormat ascii crashes when using writeFormat binary
DescriptionSee the file attached:
- the case in filter3d_ascii.tar.gz runs normally

However, when I change from writeFormat ascii; to writeFormat binary; in the file system/controlDict, the run crashes.
Steps To Reproducemake a directory for the case, eg ascii

copy the file, and cd, to the directory; then tar xvf filter3d_ascii.tar.gz

you will need to edit the Makefile to reflect your local installation of OpenFOAM, ie the values of FOAMBASH and FOAMCLEAN will need to be updated

cd .. and then cp -r ascii binary

in the ascii folder
-issue the command make clean
-issue the command make & (takes ~ 30 minutes to complete on my local workstation)

in the binary folder
- edit the file system/controlDict and change the value of writeFormat from ascii to binary
-issue the command make clean
-issue the command make (takes about 1-2 minutes to crash on my local workstation)


Additional Informationcase based on lagrangian/reactingParcelFoam tutorial filter

extended to 3D by increasing z-extent and separating frontAndBack patch into separate front and back patches

in 2.2.x, this issue did not arise

TagsNo tags attached.

Activities

SamMallinson

2014-02-26 23:19

reporter  

filter3d_ascii.tar.gz (638,370 bytes)

will

2014-02-28 17:42

manager   ~0002915

Your momentum source fvOptions are specified with points directly on vertices at the mesh boundary. In the binary case, the program is failing to find the cell which contains these points. reactingParcelFoam prints the following warning, indicating this:

Selecting finite volume options model type scalarSemiImplicitSource
    Source: energySource1
    - applying source at time 0.2 for duration 2
    - selecting cells using points
--> FOAM Warning :
    From function option::setCellSet()
    in file fvOptions/fvOption.C at line 131
    Unable to find owner cell for point (2.75 0.5 0)
    - selected 0 cell(s) with volume 0

The ascii setup will generate very small perturbations in saving and re-writing the mesh. These are enough for the specified point to move enough into a cell for it to be found.

I suggest you move your points so they are unambiguously in a cell. Changing all the values of (2.75 0.5 0) in fvOptions to (2.76 0.51 0.01) makes the case run for me.

SamMallinson

2014-03-04 04:40

reporter   ~0002925

Hi Will,

Thanks for that. I can confirm that your suggestion works.

It is strange, however, that I did not encounter this problem under 2.2.x. It is not a big problem for me, I am just checking that all of our main 2.2.x cases run OK under 2.3.x.

Cheers,

Sam

Issue History

Date Modified Username Field Change
2014-02-26 23:19 SamMallinson New Issue
2014-02-26 23:19 SamMallinson File Added: filter3d_ascii.tar.gz
2014-02-28 17:42 will Note Added: 0002915
2014-03-04 04:40 SamMallinson Note Added: 0002925
2014-03-17 15:59 will Status new => resolved
2014-03-17 15:59 will Resolution open => fixed
2014-03-17 15:59 will Assigned To => will