CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

High local velocities in interDyMFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   May 23, 2013, 04:48
Default High local velocities in interDyMFoam
Join Date: Oct 2012
Posts: 33
Rep Power: 5
jrrygg is on a distinguished road

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?

Kind regards,


jrrygg is offline   Reply With Quote

Old   June 28, 2013, 22:38
New Member
Join Date: Apr 2013
Posts: 1
Rep Power: 0
dawnrain is on a distinguished road
hi Jone
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;
volScalarField laplaceU
+ mag(fvc::laplacian(U.component(1)))
+ mag(fvc::laplacian(U.component(2))))
add following lines in interDyMFoam.C after pressure possion equation loop
// --- calculate laplacianU for refinement
laplaceU = (mag(fvc::laplacian(U.component(0)))
+ mag(fvc::laplacian(U.component(1)))
+ mag(fvc::laplacian(U.component(2))));
// --- 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
field alpha1;
lowerRefineLevel 0.001;
upperRefineLevel 0.999;
field laplaceU;
lowerRefineLevel 4; //test
upperRefineLevel 1e10;
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!
dawnrain is offline   Reply With Quote



Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
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
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 9 January 3, 2012 09:06
Multiple floating objects CKH OpenFOAM 10 September 21, 2011 23:13
calculation diverge after continue to run zhajingjing OpenFOAM 0 April 28, 2010 04:35

All times are GMT -4. The time now is 09:40.