|June 14, 2012, 08:07||
Flow around cylinder, forcing term
Join Date: Feb 2012
Posts: 6Rep Power: 5
I am a new FOAM user and I am trying to model an oscillating (sinusoidal) flow around a free circular cylinder in 2D. I have some trouble, trying to add a forcing term to the UEqn. (I use a pimpleFoam solver).
The UEqn looks like this:
(Um1*omega*Foam::cos(omega*runTime.value()) + slope)*forceDir
+ fvm::div(phi, U)
if (ocorr != nOuterCorr-1)
solve(UEqn == -fvc::grad(p)); // forcing is included in UEqn
It all runs fine with no errors, however the velcity accelerates. This means that if Um1 is eg. 0.2 m/s, then in the beginning Ux oscillates between (-0.2 ; 0.2). But after some time it increases/decreases to oscillate between eg. (0 ; 0.4), (0.2 ;0.6), (0.4 ; 0.8) etc.
I'm thinking, that a pressure gradient is developing, but I don't know why..
I have zeroGradient conditions for both velocity and pressure at inlet+outlet. The top and bottom boundaries are both slip.
Does anyone have an idea of, what I have done wrong? (I use OpenFoam 1.6-ext).
Thanks in advance
Last edited by krsp; June 15, 2012 at 04:33.
|Thread||Thread Starter||Forum||Replies||Last Post|
|benchmark: flow over a circular cylinder||goodegg||Main CFD Forum||12||January 22, 2013 12:47|
|flow around a cylinder||pXYZ||Main CFD Forum||14||July 25, 2011 10:05|
|Flow induced vibration of a mobile cylinder||Hooman||Main CFD Forum||0||December 31, 2010 09:48|
|pressure gradient term in low speed flow||Atit Koonsrisuk||Main CFD Forum||2||January 10, 2002 11:52|
|Inviscid Drag at subsonic, subcritical Mach #||Axel Rohde||Main CFD Forum||1||November 19, 2001 13:19|