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

Von Karman street simulation with icoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 7, 2008, 15:11
Default I'm trying to run a 2-D vortex
  #1
New Member
 
Anant Grewal
Join Date: Mar 2009
Posts: 6
Rep Power: 8
agrewal is on a distinguished road
I'm trying to run a 2-D vortex shedding simulation with a circular cylinder in icoFoam. The conditions are u=1 m/s, d = 0.05 m and the comp domain is 1 m in the x direction and 0.5 m in the y-direction.

The BCs are inlet, outlet and symmetric at the top and bottom.

I am using a integration time step of 1e-4 with the following schemes:

ddtSchemes
{
default CrankNicholson 0.5;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
}

laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

I also used the following solver settings:

solvers
{
p PCG
{
tolerance 1e-07;
relTol 0;
preconditioner DIC;
};
U smoothSolver
{
tolerance 1e-07;
relTol 0;
smoother GaussSeidel;
nSweeps 1;
};
}

PISO
{
momentumPredictor yes;
nCorrectors 2;
nNonOrthogonalCorrectors 2;
pRefCell 0;
pRefValue 0;
}

Before shedding starts, the Courant no blows up and the following message is seen:

#0 Foam::error::printStack(Foam:stream&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7f08420]
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&, unsigned char, int) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"ddtSche mes
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/icoFoam"
#7 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/icoFoam"
#8 main in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/icoFoam"
#9 __libc_start_main in "/lib/i686/libc.so.6"
#10 Foam::regIOobject::readIfModified() in "/home/caelinux/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/icoFoam"

Until 0.269 sec, the max Courant # is stable at 0.088, and then starts ramping up.

I use an unstructured mesh with ~9500 tetrahedral prisms. As the solution develops, the outlet region has high velocity including flow reversal that isn't transported from the cylinder. CheckMesh does not report any errors. What am I doing wrong?
Thanks
Anant
agrewal is offline   Reply With Quote

Old   February 7, 2008, 17:22
Default BTW, The Reynolds no is approx
  #2
New Member
 
Anant Grewal
Join Date: Mar 2009
Posts: 6
Rep Power: 8
agrewal is on a distinguished road
BTW, The Reynolds no is approx 350. I am reducing the time step by 10 to see if this overcomes the problem.
Anant
agrewal is offline   Reply With Quote

Old   February 9, 2008, 04:14
Default Email your case/ 0, constant a
  #3
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: South Bend, IN, USA
Posts: 688
Blog Entries: 9
Rep Power: 12
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Email your case/ 0, constant and system to me, and let me have a look.

Daniel
__________________
~
Daniel WEI
-------------
NatHaz Modeling Laboratory
Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame, USA
Email || My Personal CFD Blog
lakeat is offline   Reply With Quote

Old   February 9, 2008, 18:12
Default what's about your mesh? Is th
  #4
Senior Member
 
Marhamat Zeinali
Join Date: Mar 2009
Location: Tehran, Tehran, iran
Posts: 107
Rep Power: 8
marhamat is on a distinguished road
what's about your mesh?
Is the mesh orthogonal enough?
marhamat 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
Video about Karman vortex street simulation Marek Main CFD Forum 0 August 20, 2008 07:54
low Re (30-170) Karman vortex street Sponi FLUENT 0 March 16, 2007 14:27
Karman vortex street gellert OpenFOAM Running, Solving & CFD 15 June 17, 2005 12:19
karman vortex street michael CD-adapco 5 April 30, 2003 05:04
Von Karman Street JoHn@ Main CFD Forum 0 April 21, 2003 17:12


All times are GMT -4. The time now is 01:59.