|May 23, 2013, 04:48||
High local velocities in interDyMFoam
Join Date: Oct 2012
Posts: 33Rep Power: 6
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?
|June 28, 2013, 22:38||
Join Date: Apr 2013
Posts: 1Rep Power: 0
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!
|Thread||Thread Starter||Forum||Replies||Last Post|
|High Courant Number @ icoFoam||Artex85||OpenFOAM Running, Solving & CFD||11||February 16, 2017 14:40|
|Multiple floating objects||CKH||OpenFOAM Running, Solving & CFD||12||March 21, 2016 14:05|
|dynamic Mesh is faster than MRF????||sharonyue||OpenFOAM Running, Solving & CFD||14||August 26, 2013 07:47|
|Problem with FloatingObject||Leech||OpenFOAM Running, Solving & CFD||10||March 29, 2012 15:24|
|calculation diverge after continue to run||zhajingjing||OpenFOAM||0||April 28, 2010 04:35|