View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0004245OpenFOAMBugpublic2025-06-17 11:512025-06-17 11:51
Reportercgoessni Assigned To 
PrioritynormalSeverityminorReproducibilityalways
Status newResolutionopen 
Product Versiondev 
Summary0004245: streamlines with boundaryRandom sampling, executeAtStart false and moving mesh doesn't work
DescriptionThe streamlines function object yields an error when the following conditions are met:

1. boundaryRandom seeding
2. executeAtStart false
3. moving mesh case
Steps To ReproduceAdd the following to XiFluid/kivaTest/system/functions:

streamlines
{
    type streamlines;

    libs ("libfieldFunctionObjects.so");

    // Output every
    writeControl writeTime;

    // Write format
    setFormat vtk;

    // Track forward (+U) or backward (-U) or both
    direction forward;

    // Names of fields to sample. Should contain above velocity field!
    fields (p k U);

    // Steps particles can travel before being removed
    lifeTime 10000;

    // Number of steps per cell (estimate). Set to 1 to disable subcycling.
    nSubCycle 5;

    // Cloud name to use
    cloudName particleTracks;

    // Seeding method.
    seedSampleSet
    {
        type boundaryRandom;
        patches (cylinderHead);
        nPoints 10;
    }
    executeAtStart no;
}


The solver crashes at the first execution of the streamlines function object at t = -175 CAD with the following error:

Courant Number mean: 0.000232674 max: 0.0123949
Time = -175CAD

streamlines streamlines:
    automatic track length specified through number of sub cycles : 5

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.0187947, Final residual = 8.3518e-08, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.0286253, Final residual = 1.31213e-07, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.0323516, Final residual = 8.72553e-08, No Iterations 1
DILUPBiCGStab: Solving for ft, Initial residual = 0.000139999, Final residual = 7.70679e-08, No Iterations 1
DILUPBiCGStab: Solving for ha, Initial residual = 0.12134, Final residual = 4.24599e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.642688, Final residual = 0.0635741, No Iterations 8
DICPCG: Solving for p, Initial residual = 0.126738, Final residual = 0.00890416, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.32315e-06, global = 4.25307e-07, cumulative = 3.6625e-06
DICPCG: Solving for p, Initial residual = 0.0620049, Final residual = 0.00557289, No Iterations 5
DICPCG: Solving for p, Initial residual = 0.0336373, Final residual = 7.64805e-07, No Iterations 28
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.12574e-10, global = -8.62431e-12, cumulative = 3.6625e-06
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.00243534, Final residual = 1.20503e-06, No Iterations 1
DILUPBiCGStab: Solving for k, Initial residual = 0.0316659, Final residual = 1.33305e-07, No Iterations 1
ExecutionTime = 2.70091 s ClockTime = 3 s

Piston:
    position from TDC: 0.0841156 m
    speed = 0.401661 m/s
    displacement: 1.11572e-05 m
    clearance: -1.11572e-05 m

streamlines streamlines write:


--> FOAM FATAL ERROR:
Old cell centres have not been stored

    From function virtual const pointField& Foam::polyMesh::oldCellCentres() const
    in file meshes/polyMesh/polyMesh.C at line 1405.

FOAM exiting
TagsNo tags attached.

Activities

There are no notes attached to this issue.

Issue History

Date Modified Username Field Change
2025-06-17 11:51 cgoessni New Issue