View Issue Details

IDProjectCategoryView StatusLast Update
0003894OpenFOAMBugpublic2022-10-21 12:59
Reporterblttkgl Assigned Tohenry  
PrioritynormalSeveritycrashReproducibilityalways
Status resolvedResolutionfixed 
PlatformGNU/LinuxOSOtherOS Version(please specify)
Product Version10 
Fixed in Versiondev 
Summary0003894: Mesh refiner cannot access rho_U field from objectRegistry
DescriptionWhen refiner type topoChanger is used with a compressible solver the topoChanger cannot find the vector rhoU_0 from the object registry and throws the following error. I also found this bug report from last year which may be relevant: https://bugs.openfoam.org/view.php?id=3636

--> FOAM FATAL ERROR:

    request for volVectorField rhoU_0 from objectRegistry region0 failed
    available objects of type volVectorField are

3
(
U_0
U
rhoU
)

Steps To ReproduceRun the aachenBomb tutorial with a refiner topoChanger, using the C7H16 field for refining criteria.
TagsNo tags attached.

Activities

peksa

2022-10-09 18:17

reporter   ~0012783

I was able to reproduce the error by introducing the attached dynamicMeshDict under constant directory in the aachenBomb tutorial. As the error message suggests, fvMeshTopoChangers::refiner::refineUfs() tries to access the rhoU_0 which is not available.
dynamicMeshDict (1,681 bytes)   
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  dev
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

topoChanger
{
    type            refiner;

    libs            ("libfvMeshTopoChangers.so");

    // Refine every refineInterval timesteps
    refineInterval 1;

    // Maximum refinement level (starts from 0)
    maxRefinement 2;

    // Maximum cell limit (approximate)
    maxCells 100000000;

    // volScalarField to base refinement on
    field C7H16;

    nBufferLayers 1;

    dumpLevel true;

    lowerRefineLevel 0.0;
    upperRefineLevel 1.0;

    nBufferLayers 1;
    // Newly introduced patch points optionally get projected onto a surface
    //projectSurfaces ("fixedWalls4.stl");
    //projectPatches (fixedWalls);
    // Maximum project distance
    //projectDistance 1;

    // Fluxes to adapt. For newly created faces or split faces the flux
    // gets estimated from an interpolated volVectorField ('velocity')
    // First is name of the flux to adapt, second is velocity that will
    // be interpolated and inner-producted with the face area vector.
    correctFluxes
    (
        (phi none)
    );

}


// ************************************************************************* //
dynamicMeshDict (1,681 bytes)   

henry

2022-10-10 21:45

manager   ~0012786

Resolved by commit 30eb5e28e65d359dbe8c350b9073dcd492172904

Issue History

Date Modified Username Field Change
2022-10-05 10:56 blttkgl New Issue
2022-10-09 18:17 peksa Note Added: 0012783
2022-10-09 18:17 peksa File Added: dynamicMeshDict
2022-10-10 21:45 henry Assigned To => henry
2022-10-10 21:45 henry Status new => resolved
2022-10-10 21:45 henry Resolution open => fixed
2022-10-10 21:45 henry Fixed in Version => dev
2022-10-10 21:45 henry Note Added: 0012786