# Oscillating cylinder using icoFoam blowing up

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 September 21, 2010, 07:40 Oscillating cylinder using icoFoam blowing up #1 Senior Member   Pavan Join Date: May 2009 Location: Melbourne Posts: 101 Rep Power: 9 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? utkunun likes this.

 September 22, 2010, 05:02 #2 Senior Member   Pavan Join Date: May 2009 Location: Melbourne Posts: 101 Rep Power: 9 Hmm does my acc term need to be scaled by density of the fluid?...

 September 22, 2010, 09:22 #3 Senior Member     ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 322 Rep Power: 10 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

 September 22, 2010, 09:25 #4 Senior Member   Pavan Join Date: May 2009 Location: Melbourne Posts: 101 Rep Power: 9 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...

 September 22, 2010, 09:32 #5 Senior Member     ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 322 Rep Power: 10 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

 September 22, 2010, 11:25 #6 Senior Member   Pavan Join Date: May 2009 Location: Melbourne Posts: 101 Rep Power: 9 What's wrong with them?

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post CH FLUENT 8 April 2, 2014 10:56 startingcfd Main CFD Forum 7 March 23, 2012 06:49 ammi FLUENT 0 December 18, 2006 05:06 zonexo FLUENT 3 September 14, 2006 10:13 fredD FLUENT 0 July 16, 2004 09:32

All times are GMT -4. The time now is 11:35.

 Contact Us - CFD Online - Top