[pradeepramesh]

Dear All,

I am currently trying to solve linear and non-linear Poisson equation with implicit source terms in OpenFoam. I would like to request your feedback.

Case I : linear Poisson equation
Solve for : del^2(T) - T = 0

Source : S = -T
Implicit source: S = Su + Sp*T
Su = 0, Sp = -1

Open Foam: fvm::laplacian(T) == -fvm::Sp(Sp,T) + Su

I defined Sp and Su as "VolScalarField". I was successfully able to solve the case and compared it with the analytically results.

************************************************** ***********

Case II: Non-linear Poisson equation
Solve for : del^2(T) - T^2 = 0

Source : S = -T^2
Implicit source: S = Su + Sp*T
With taylor's series expansion:

S = To^2 + (- 2To)*T

Su = To^2, Sp = (- 2To)

I am facing errors to implement these sources terms.
How to mimic these equation in OpenFoam syntax?
How to define Su and Sp in terms of To.?

Please find the attached Snapshot for your reference.
Would be great to hear your feedback on Case II. Thank you very much and have a wonderful day.

Cheers,
Pradeep

[Cyp]

Hi,

It would have been easier if you post the error message. The implementation is similar to your first case. I guess you should program Foam:ow(T0,2) for T0^2. Maybe your compilation error comes from here.

Another point you have to be careful is the dimension of T ! Indeed, despite the case where T is dimensionless, your equation is not correct (at least you should put an dimensionedScalar factor in front of your source term in order to be homogeneous).

Best,
Cyp

[pradeepramesh]

Dear Cpy,

Thank you very much for your feedback and i highly appreciate that.

Could you please elaborate on the dimensionedScalar factor?


Case I: ----> fvm::laplacian(T)== -fvm::Sp(Sp,T) + Su

dimensionedScalar factor for source terms, you mean for Sp and Su???


Case II: ----> fvm::laplacian(T)== T*T*A

Where: A is a dummy object to make the operation compatible.


Thank you very much

Cheers,
Pradeep

[Cyp]

For both cases, I will declare A as :

dimensionedScalar A ("A", dimless/T.dimensions()/dimLength/dimLength, 1.0);

and then

fvm::laplacian(T) == -fvm::Sp(2*T0*A,T) + Foam::pow(T0,2)*A

Best,
Cyp

[pradeepramesh]

Dear Cpy,

I have updated the new syntax and with the new dimension syntax,I am facing the incompatibility error for operation. Please find the attached snapshot for your reference.

Thank you very much.

Cheers,
Pradeep

[Cyp]

Indeed it is obvious (because of the laplacian..). I modified my previous post.

[pradeepramesh]

It worked! Thank you very much Cpy

[ancolli]

Could you please explain how to modify code to solve a general nonlinear Poisson equation?
Step by Step if is it possible, because I do not fully understand