[sisetrun]

Hey Foamers,

I have run a simpleFoam simulation of an internal pipe flow.
Since the setup consists of repeated obstacles I want to cut out one repetitive element and run scalarTransportFoam on the smaller mesh. (See attached picture).

I tried to map the fields with mapFields to the a mesh of the repetitive element, but I always get the following error:

Source time: 0.1
Target time: 0.1
Create meshes

Source mesh size: 20430247      Target mesh size: 1753323


Mapping fields for time 0.1

    interpolating p
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  Foam::meshToMesh0::calculateInverseDistanceWeights() const at ??:?
#4  Foam::meshToMesh0::inverseDistanceWeights() const at ??:?
#5  ? at ??:?
#6  ? at ??:?
#7  ? at ??:?
#8  ? at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib64/libc.so.6"
#11  ? at ??:?
Segmentation fault

Mesh stats
    points:           2193964
    faces:            5623815
    internal faces:   5195678
    cells:            1753323
    faces per cell:   6.17085
    boundary patches: 7
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     1419493
    prisms:        57879
    wedges:        0
    pyramids:      0
    tet wedges:    4342
    tetrahedra:    73
    polyhedra:     271536
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            4   42311
            5   49339
            6   64796
            7   10229
            8   3530
            9   53004
           10   1244
           11   735
           12   34174
           13   177
           14   63
           15   11570
           16   6
           18   358

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology
                   inlet    50583    57911  ok (non-closed singly connected)
                  outlet    50313    57594  ok (non-closed singly connected)
                    yneg    27128    31221  ok (non-closed singly connected)
                    ypos    27116    31169  ok (non-closed singly connected)
                    xneg    27119    31159  ok (non-closed singly connected)
                    xpos    27134    31127  ok (non-closed singly connected)
                 lattice   218744   295009  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-0.005000008 -0.005000033 0.01499371) (0.00500001 0.005006101 0.0200002)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-6.950294e-17 9.097023e-16 -2.955753e-15) OK.
    Max cell openness = 4.240475e-16 OK.
    Max aspect ratio = 6.275022 OK.
    Minimum face area = 1.045066e-11. Maximum face area = 2.01948e-08.  Face area magnitudes OK.
    Min volume = 2.395032e-16. Max volume = 2.015212e-12.  Total volume = 3.813133e-07.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.1619 average: 10.9496
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.823711 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

I thought of memory problems since all works fine for a smaller element than I want to use for the postProcessing (samller than the 1.7Mio shown above). But running in parallel did not work out as well...

Can anybody tell me if this is possible in a different way and how?

Thanks a lot for your help!

Best regards,

Sebastian