High local velocities in interDyMFoam
I have a simulation where a high speed jet of water (38.38 ms-1) interacts with a rotating turbine bucket in interDyMFoam. I have a problem with the courant number being very low (10-8), much because of an extremely high velocity (10⁵ ms-1) in a few cells around the bucket.
I wonder if anyone else has seen this behavior in interDyMFoam, and maybe knows how I can avoid velocities of this magnitude?
Maybe you can define another parameter in creatFields.H to replace Alpha1 for refinement. Like this
Add following lines in creatFields.H
Info<< "Creating field Laplacian U\n" << endl;
add following lines in interDyMFoam.C after pressure possion equation loop
// --- calculate laplacianU for refinement
laplaceU = (mag(fvc::laplacian(U.component(0)))
// --- calculate average laplacianU
dimensionedScalar laplaceUAver =
sum(laplaceU * mesh.V()) / sum(mesh.V());
// --- laplacianU
laplaceU = laplaceU / laplaceUAver;
Logicallly, the field laplaceU will be high on the interface and area with sharp velocity gradient.
Don't forget modify costant/dynamicMeshDict, like
lowerRefineLevel 4; //test
But i am not sure whether the change will be accepted by interDyMFoam, because adaptive mesh refinement code is difficult for me. Hope you succeed!
|All times are GMT -4. The time now is 13:42.|