Quote:
fvm::Sp( fvc::laplacian((X/Y*mag(Z)), U)? Notice I have fvc so that it solves for the laplacian explicit and then it multiplies U and that goes to the diagonal, but once again I've never tried this so that might be wrong. |
Quote:
Thanks for a quick reply. But no, this doesn't work. The error is: Code:
error: no matching function for call to ‘Sp(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >)’ Code:
note: candidate expects 2 arguments, 1 provided |
Quote:
I'm just curious if fvm::laplacian compiles? There's also SuSp https://www.cfd-online.com/Forums/op...treatment.html |
Quote:
The core issue I think is handling this term here: "X/Y * mag(Z)". |
Quote:
|
Quote:
Well I was thinking the same. But I just tried and no. Compilation fails with same errors again. Needs one more argument (:confused:). |
Quote:
|
Quote:
Compilation worked with simple fvm::laplacian in my first post. But case failed citing wrong dimensions. This was bound to happen because I was prompted to have a term in fvSchemes for a physically unrealistic quantity (e.g. laplacian of X/mag(Z)*mag(Z) instead of expected laplacian of X/mag(Y)*mag(Z).. as in the original expression). This led me think I'm not representing it right. With fvm::Sp it never managed to compile. Same error everytime. I just also tried changing the sequence of Sp and laplacian. Doesn't work. Code:
fvm::laplacian(fvc::Sp(1/Y, fvc::Sp(X,mag(Z))), U) |
Quote:
|
Quote:
Weird error: Code:
--> FOAM FATAL IO ERROR: What does that even mean!? |
Quote:
If you could post your schemes file that might help someone see something that was overlooked or maybe you could just upload the parts of your project that are necessary to see the problem, not the whole thing (unless you don't mind). Best regards! |
Thanks for your help so far nonetheless. I also tried to change the order of terms and it doesn't change the situation much.
I'll try to share more if I can soon. |
Quote:
Basically you solve for that part separately and then solve the Laplacian by taking the divergence of the gradient literally. |
Quote:
So, I was wondering if there's any way to print the dimensions of a certain quantity (with or without the value) as the solver runs? I bumped up an old thread here if you want to see what I mean. |
Quote:
Have you solved your problem? I would like to ask you if fvm::laplacian can have three terms or it just can have 2 terms. Have you tried with something like this? Code:
- fvm::laplacian Actually, I'm not sure what fvc::interpolate does but there's a product and division at the second term nut*rho/Sc .... Regards, |
All times are GMT -4. The time now is 04:31. |