Solve dS_dt + div(vect(U)*f(S))=q
Hi
I am trying to solve this equation in openfoam. dS_dt + div(vect(U)*f(S))=q where f(S) is a known equation dependent of S. S cannot be explicitly extracted from f(S). up to now if I implement the equation as stated above I get an error concerning calculatedFvPatchField. Has anybody had similar problems? Thanks Andres |
A bit vague :-)
You should give some more detail on how the equation was implemented, and be specific with the error message. Best, |
Impementaiton
Hi Alberto,
I hope that the problem declaration this time is more illustrating. Implementation: The initial values of Code:
volScalarField S Code:
dimensionedScalar const1 The implementation reads Code:
volScalarField X = (S-const1)/(const2); Code:
Error: Thanks Andres |
try to replace
fvm::div(phi,f) with fvc::div(phi,f) You cant solve implicitly for both f and S, they arent really the same fields are they? |
All times are GMT -4. The time now is 15:40. |