CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Oscillating cylinder using icoFoam blowing up

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 21, 2010, 07:40
Default Oscillating cylinder using icoFoam blowing up
  #1
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
Hey guys I'm simulating a (2D) transversely oscillating cylinder using a non-inertial reference frame using the following modification to icoFoam:

Code:
scalar Pi = mathematicalConstant::pi;
dimensionedScalar accCt("accCt",dimensionSet(0,1,-2,0,0,0,0),-3.0*pow(2.0*Pi/3.0,2)*Foam::cos(2.0*Pi*runTime.value()/3.0));
dimensionedVector acc("acc",dimensionSet(0,1,-2,0,0,0,0), vector::zero);

acc.component(vector::Y) = accCt;

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu, U)
    ==
       acc
        );


        solve(UEqn ==  - fvc::grad(p));
However the run keeps blowing up (Floating Point Exception with ridiculously high velocity magnitude) at around the same time even for very small delta T... can anyone see what's wrong with my implementation?
rieuk is offline   Reply With Quote

Old   September 22, 2010, 05:02
Default
  #2
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
Hmm does my acc term need to be scaled by density of the fluid?...
rieuk is offline   Reply With Quote

Old   September 22, 2010, 09:22
Default
  #3
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 8
ata is on a distinguished road
Hi
I hope you are doing very well.
I am not sure that your solver is true. But may I ask you how implement geometry and boundary conditions?
Best regards
Good luck

Ata
ata is offline   Reply With Quote

Old   September 22, 2010, 09:25
Default
  #4
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
What do you mean you're not sure, as in you think it's wrong? Well with velocity boundary conditions I intend to use dirichlet conditions (freestream + sinusoidal cylinder velocity component) for inlet and top/bottom and zero gradient for outlet. I'm thinking groovyBC is the only way to handle that but for now I'm just using fixedValue to see if the solver is doing what it's supposed to. Right now it's just doing crazy things like creating random vortices emanating from the boundaries...
rieuk is offline   Reply With Quote

Old   September 22, 2010, 09:32
Default
  #5
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 8
ata is on a distinguished road
no no no no
Excuse me. I do not mean your code is incorrect. At least in my knowledge.
I think may be your BCs are bad.
good luck
ata is offline   Reply With Quote

Old   September 22, 2010, 11:25
Default
  #6
Senior Member
 
Pavan
Join Date: May 2009
Location: Melbourne
Posts: 101
Rep Power: 8
rieuk is on a distinguished road
What's wrong with them?
rieuk is offline   Reply With Quote

Reply

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
tutorial for solving oscillating cylinder problem CH FLUENT 8 April 2, 2014 10:56
Oscillating Cylinder in Parallel Flow: Analytical Solution startingcfd Main CFD Forum 7 March 23, 2012 06:49
problem with oscillating cylinder ammi FLUENT 0 December 18, 2006 05:06
problem with oscillating cylinder zonexo FLUENT 3 September 14, 2006 10:13
UDF for oscillating cylinder fredD FLUENT 0 July 16, 2004 09:32


All times are GMT -4. The time now is 14:06.