|February 14, 2013, 05:44||
Don't understand solver's behaviour
Join Date: Nov 2009
Posts: 36Rep Power: 8
today I'm facing a problem, which I can't explain. I have done a dummy calculation for a variable f. f is set everywhere to 1 (boundary/internal). Nevertheless OpenFOAM calculates values differing from 1. Can anyone explain that behaviour to me? Thank you very much in advance.
I have a normal transport equation for variable f:
f_( IOobject ( "f", mesh.time().timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ),
Boundary and internal values all set to 1.
multivariateSurfaceInterpolationScheme<scalar>::fielTable fields; fields.add(f); tmp<fv::convectionScheme<scalar> > convection ( fv::convectionScheme<scalar>::New ( mesh_, fields, phi, mesh_.divScheme("div(phi,f)") ) );
So again keep in mind, for testing purpose I set all values to 1. Now I'm trying to solve the following equations:
Info << "fmin/max = " << min(f_).value() << ", " << max(f_).value() << endl;; tmp<fvScalarMatrix> fEqn ( fvm::ddt(rho, f_) + convection->fvmDiv(phi, f_) - fvm::laplacian(turbulence_->alphaEff(), f_) ); fEqn().relax(); solve(fEqn == rho* fSource);
Info << "fmin/max = " << min(f_).value() << ", " << max(f_).value() << endl;;
f min/max = 1, 1
f min/max = 0.90082369, 1.1320949
laplacian(alphaEff,f) Gauss linear corrected;
div(phi,f) Gauss upwind;
Last edited by galap; February 14, 2013 at 08:33.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Unexpected behaviour of timeVaryingMappedFixedValue||Kopy||OpenFOAM Running, Solving & CFD||3||March 15, 2013 12:27|
|Particle behaviour in cfd||RedAdmiral||Main CFD Forum||0||November 18, 2012 10:51|
|Unstable behaviour after long period of stablility||plunge11||FLUENT||1||April 6, 2011 09:15|
|strange behaviour of GGI in parallel on axis symmetrical case||A.Devesa||OpenFOAM Running, Solving & CFD||0||April 6, 2010 03:58|
|Modelling Industrial cyclone behaviour||GŁnther Hasse||Main CFD Forum||3||October 12, 1999 19:34|