|
[Sponsors] |
June 6, 2012, 09:37 |
Question about the fvm::SuSp, and fvm::Sp
|
#1 |
Member
Charlie
Join Date: Dec 2010
Location: USA
Posts: 85
Rep Power: 15 |
Hi Foamers!
When I'm programming with the OpenFOAM, I got confused with the difference between fvm::SuSp, and fvm::Sp, could anybody tell me the difference? As far as I know, fvm::SuSp(a,b) will calculate a*b implicitly or explicitly according to the sign of a, however, when a is positive, is fvm::Sp explicit or implicit? fvm::Sp(a,b) will calculate a*b implicitly anyway. The second question is related to a phenomenon I found when I'm solving the continuity equation, the exact form should be : ddt(alpha) + div(phia, alpha) - alpha * div(phia) ==0 This supposed to be solving ddt(alpha) + U_a*grad(alpha) =0, and the implementation I made in the code is : (alpha is a scalar, volume fraction) fvScalarMatrix alphaEqn ( fvm::ddt(alpha) + fvm::div(phia, alpha, scheme) - fvm::Sp(fvc::div(phia), alpha, schemer) ); where phia = fvc::interpolation(Ua) *mesh.Sf() Interestingly, after this, I got that alpha does not change with time anymore, which means that the equation reduced to ddt(alpha) =0, which is clearly not true if we take the definition of fvm::Sp(). However, my worry is that, in the momentum equation, they are solving : fvm::ddt(Ua) + fvm::div(phia, Ua, "div(phia,Ua)") - fvm::Sp(fvc::div(phia), Ua) == ........ I'm wondering if this fvm::Sp( ) makes a difference between the scalar source and vector source. Thanks!!!! Any idea would be appreciated!! |
|
June 6, 2012, 10:06 |
Sorry for the mistake I made in the continuity equation
|
#2 |
Member
Charlie
Join Date: Dec 2010
Location: USA
Posts: 85
Rep Power: 15 |
Hi
I made a mistake, I think the exact form for continuity equation should be: ddt(alpha) + div(alpha*phia) =0 not : ddt(alpha) + U_a*grad(alpha) =0 So, forget about the continuity equation I described, however, the question about the difference between fvm::SuSp() and fvm::Sp() is still valid, Thanks! Zhen |
|
June 30, 2018, 03:29 |
I have the same doubt as well
|
#4 |
New Member
Amal
Join Date: May 2018
Posts: 3
Rep Power: 7 |
The page mentioned doesn't give much idea more than what Cheng has already mentioned.
Can someone please tell, in which case I should use Sp and where I should use SuSp? Thank you |
|
June 30, 2018, 13:03 |
|
#5 |
Senior Member
Sergei
Join Date: Dec 2009
Posts: 261
Rep Power: 21 |
Hi! Here you can find some more info on SuSp:
fvMatrix, fvOptions and SuSp: automatic implicit/explicit source-term treatment |
|
|
|