|
[Sponsors] |
April 1, 2015, 07:08 |
code about continuity errors
|
#1 |
New Member
nemo
Join Date: Jan 2015
Posts: 26
Rep Power: 11 |
Hello,everyone.I just start to read the code of simpleFoam,and I have some problems with the continuity errors.Here is the code:
33 volScalarField contErr(fvc::ddt(alphac) + fvc::div(alphacf*phic)); 35 scalar sumLocalContErr = runTime.deltaTValue()* 36 mag(contErr)().weightedAverage(mesh.V()).value(); 37 38 scalar globalContErr = runTime.deltaTValue()* 39 contErr.weightedAverage(mesh.V()).value(); 40 cumulativeContErr += globalContErr; I read this thread ,but I still have some questions: (1) phi=U&mesh.Sf() div(phi) = div(U) = div(U&mesh.Sf()) is this still right? According to the code,I think div(phic)=div(U),which one is right? (2) what does the alphac and alphacf mean? (3) does mag(scalar) mean the absolute value of the scalar? Thanks. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
should Courant number always be kept below 1? | wc34071209 | OpenFOAM Running, Solving & CFD | 16 | March 9, 2014 19:31 |
Micro Scale Pore, icoFoam | gooya_kabir | OpenFOAM Running, Solving & CFD | 2 | November 2, 2013 13:58 |
AMI interDyMFoam for mixer nu problem | danny123 | OpenFOAM Programming & Development | 8 | September 6, 2013 02:34 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 15:33 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 18:07 |