[Petr Kazarin]

Hello guys!
I am trying to use dynamicRefineFVMesh based on velocity magnitude field.
Here is my dict:

dynamicFvMesh   dynamicRefineFvMesh;

// How often to refine
refineInterval  1;

// Field to be refinement on
field         U;

// Refine field inbetween lower..upper
lowerRefineLevel 0.1;
upperRefineLevel 13;

// If value < unrefineLevel unrefine
unrefineLevel   10;

// Have slower than 2:1 refinement
nBufferLayers   1;

// Refine cells only up to maxRefinement levels
maxRefinement   1;

// Stop refinement if maxCells reached
maxCells        2000000;

As an output i get this ERROR:

--> FOAM FATAL ERROR: 

    lookup of U from objectRegistry region0 successful
    but it is not a volScalarField, it is a volVectorField

    From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]
    in file /home/kepchuk/OpenFOAM/OpenFOAM-v1706/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 226.

FOAM aborting

However, on the OpenfoamWIKI it is said that

If a vector field is used, the dictionary will utilize the vector magnitude of that field. Only one field can be specified for mesh refinement.

How to use U as a filed in that dict?


Regards, Petr

[matejfor]

This one is easy.
just read the error message and you can figure this out yourself ..
.
.
.
.
.
.
.
.
.
.
No?
OK.
The refinement works only on the scalar fields and your velocity is vector field.
So if you want to refine with the velocity it must be velocity magnitude or velocity component ... scalar values. You may use Function object to create magnitude of U.


happy foaming
matej

[Petr Kazarin]

Sure, I tried all this stuff

[Petr Kazarin]

Matej, of course I tried whatever I could imagine before writing this post: magnitudes, functionObjects etc. However, it does see vector U. Also, Openfoam WIKI stated that

If a vector field is used, the dictionary will utilize the vector magnitude of that field. Only one field can be specified for mesh refinement.

.
It does not see any variables calculated from "outside".

I came to the conclusion that only writing the custom version of dynamicRefineFvMesh will resolve the problem.I was just hoping that there could be some tricks which could help with the issue. May be there are some..

Best regards, Petr

[matejfor]

Well well, maybe then worth reporting ( I see you are using 1706) the issue here: develop.openfoam.com

Sorry does not have time to look deeper into it at the moment.
I was just trying to be lucky with a simple hint. Pity.

Btw. out of curiosity, why are you refining with the velocity magnitude? What is the application? I can imagine you want to use gradient (which would not solve the problem either).

[chriss85]

If you can compile the solver yourself you can also just calculate a field with the appropriate indicator and use that.