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

pimpleDyMFoam: moving cylinder in duct - crashes with high Courant Number

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 26, 2015, 09:20
Default pimpleDyMFoam: moving cylinder in duct - crashes with high Courant Number
  #1
New Member
 
shokat
Join Date: Apr 2015
Location: Iran,Amol
Posts: 9
Rep Power: 11
khatereh is on a distinguished road
Hi foamers,
i started working on the flow past of the cylinder in duct,but i have a problem to use pimpleDymfoam in my case, the cylinder move in the y direction with a fixed rate as 0.1m/s.
i used pointmotionuy in 0 document and dynamicMotionSolverFvMesh in constant doc,

motionSolverLibs ( "libfvMotionSolvers.so" );

solver velocityComponentLaplacian;

velocityComponentLaplacianCoeffs
{
component y;
diffusivity directional (1000 1 0);
}
after running.it start to solve the case,but in time=0.84 it stop
the courant number is very high in this time!!!
error:
Courant Number mean: 8.746e+09 max: 4.19433e+17
Time = 0.8454

DICPCG: Solving for cellMotionUy, Initial residual = 1.59842e-05, Final residual = 9.87955e-09, No Iterations 80
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 0.0196549, No Iterations 32
time step continuity errors : sum local = 131631, global = 906.913, cumulative = -56381.3
PIMPLE: iteration 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::Field<Foam::Tensor<double> >:perator/=(Foam::UList<double> const&) at ??:?
#4 Foam::fv::gaussGrad<Foam::Vector<double> >::gradf(Foam::GeometricField<Foam::Vector<double> , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::word const&) at ??:?
#5 Foam::fv::gaussGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const at ??:?
#6 Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const at ??:?
#7 Foam::tmp<Foam::GeometricField<Foam:uterProduct< Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) at ??:?
#8 Foam::tmp<Foam::GeometricField<Foam:uterProduct< Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#9 Foam::incompressible::laminar::divDevReff(Foam::Ge ometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
#10
at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
at ??:?
Floating point exception (core dumped)
what should i do?
thank u
khatereh is offline   Reply With Quote

Old   July 18, 2015, 12:54
Default
  #2
New Member
 
Join Date: May 2015
Posts: 29
Rep Power: 11
avd28 is on a distinguished road
Start systematically.

1. Run checkMesh. Are all your cells fine?
2. Take a look at your p and U solver in fvschemes dict, try out different options.
3.Take a look at your cell sizes and your local velocities. Unless you are dealing with very small geometries, Co issue is usually tied to improper boundary conditions
4. Check your boundary conditions. Are they stable?
avd28 is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Circular cylinder vortex shedding/strouhal number too high li0012in CFX 7 April 28, 2016 20:56
foam-extend_3.1 decompose and pyfoam warning shipman OpenFOAM 3 July 24, 2014 08:14
Visualizing mesh regions with high courant number? Endel OpenFOAM 2 August 20, 2013 03:57
AMI interDyMFoam for mixer danny123 OpenFOAM Running, Solving & CFD 4 June 19, 2013 04:49
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 18:07


All times are GMT -4. The time now is 23:12.