[TCS]

Hi everyone, I'm looking to use the interDyMFoam solver (OF Foundation 4.1).

I have a droplet defined by the volume fraction (alpha), and I want to refine both around the droplet's interface and everywhere inside. That is, I want to refine where alpha>0.01 and not coarsen for any value above 0.01 (unbounded above).

So in the end, outside the droplet should have the original mesh, which changes (adaptively) across the interface, with the mesh inside the droplet (where alpha is around 1) having the same refinement as the interface.

Does anyone know if this is possible?

[simrego]

Hi!


What if you set:
lowerRefineLevel 0.001;
upperRefineLevel 5; // max should be 1 but to be sure.


It will refine if
0.001 < alpha.water < 5, so if alpha.water > 0.001

[TCS]

Originally Posted by simrego

Hi!


What if you set:
lowerRefineLevel 0.001;
upperRefineLevel 5; // max should be 1 but to be sure.


It will refine if
0.001 < alpha.water < 5, so if alpha.water > 0.001

Hi simrego, thanks for the suggestion. I tried that and it gave me the following error which I don't understand how to deal with.

PIMPLE: iteration 1
Selected 364 cells for refinement out of 2268.
Refined from 2268 to 4816 cells.
Selected 0 split points out of a possible 364.
Execution time for mesh.update() = 0.02 s
DICPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
smoothSolver:  Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 0.130194  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.130194  Min(alpha.water) = 0  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0265856, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  ? at tensorField.C:?
#4  ? at ??:?
#5  ? at ??:?
#6  ? at ??:?
#7  __libc_start_main in "/lib64/libc.so.6"
#8  ? at ??:?
Floating point exception (core dumped)

[simrego]

There might be some other problems I think. Maybe your initial mesh is too coarse. But i'm just guessing. If it is the original dynamicRefineFvMesh (if you made no modifications in it), it should works like this without problems.

[TCS]

Originally Posted by simrego

There might be some other problems I think. Maybe your initial mesh is too coarse. But i'm just guessing. If it is the original dynamicRefineFvMesh (if you made no modifications in it), it should works like this without problems.

There was something else. My alpha=1 region was initialised such that it was in a corner of two boundaries (walls - same or separate patches, the same thing happens). If I initialise in contact with a single boundary, or free, it works fine.

Maybe this is an issue with how dynamicRefineFvMesh works?

[simrego]

I'm not sure if i understand your problem, but if you check the damBreakWithObstacle case, the alpha=1 field is also initialized in a corner and works fine.