|
[Sponsors] |
pimpleDyMFoam: moving cylinder in duct - crashes with high Courant Number |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 26, 2015, 09:20 |
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 |
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 |
|
July 18, 2015, 12:54 |
|
#2 |
New Member
Join Date: May 2015
Posts: 29
Rep Power: 11 |
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? |
|
|
|
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 |